Maximal Heights of Nearshore Storm Waves and Resultant Onshore Flow Velocities

Storm waves, after breaking or overtopping, generate strong onshore flows that do significant mechanical work, including eroding and transporting large boulders. The waves can be amplified on approach, and the flows themselves may be further intensified by local topographic effects. These processes are currently poorly parameterised, but are of great importance for understanding the interactions between waves and coasts. We present a highly generalised equation for estimating maximal coastal wave heights and consequent onshore flow velocities. Although very approximate, this method contains no embedded assumptions, and thus provides a more realistic first-order check of storm wave capabilities than previous approaches. Initial analysis suggests that exceptional wave impacts may generate onshore flow velocities up to six times greater than expected from previous approaches. Although the probability of occurrence in any given storm is very low, the possibility of such extreme values cannot be ignored, especially when interpreting ancient deposits of large boulders. The equations presented here can be used as a first-order test for coastal boulder deposits currently interpreted as tsunami deposits, to evaluate whether a storm-wave origin should be reconsidered. This approach could also be employed at coasts in general, to evaluate long-term probabilities of damaging flows, as a component of coastal risk analysis.


INTRODUCTION
The behaviour of large waves as they approach steep coasts is one of the great unknowns in marine science. In contrast to waves that shoal gradually across smooth bathymetry, and run up gentle slopes-which are well-described and amenable to mathematical modelling-waves at rocky coasts are far less studied. They are difficult to measure, and equally difficult to model (Dodet et al., 2018;Herterich and Dias, 2019). Thus, there are few data for storm-wave runup in those environments. This matters, because coastal wave energy fluxes are predicted to increase (Mentaschi et al., 2017) and reliable risk analysis for coastal populations is urgently needed (e.g., Switzer et al., 2014;Esteban et al., 2016;Breivik et al., 2017). It is therefore important to understand inshore transmission of storm-wave energy. Extreme wave behaviour can most readily be observed at exposed rocky coasts.
Coastal boulder deposits have received attention as markers for extreme wave inundation events (e.g., Nott, 2003a;Lorang, 2011;Cox et al., 2018). However, there has been controversy about the kinds of power required to emplace supratidal megagravel, i.e., boulders with an intermediate axis longer than 4.1 m (e.g., Switzer and Burston, 2010;Engel and May, 2012;Kennedy et al., 2016;Weiss and Sheremet, 2017). Some workers argue that storm waves are insufficiently powerful to move massive boulders, which are therefore often interpreted to be tsunami deposits (Scicchitano et al., 2007;Scheffers and Kinis, 2014;Roig-Munar et al., 2019a). But abundant recent evidence shows that storm waves can and do move extremely large rocks (e.g., May et al., 2015;Kennedy et al., 2017;Cox et al., 2018;Cox, 2019). New applications of surface-exposure dating to coastal boulders are revealing chronologies that conflict with tsunami interpretations (e.g. Brill et al., 2020), and new models of responses to wave loading are revealing the true power of storms (e.g., Park et al., 2018). Thus, parameterising wave and flow amplification at steep coasts can help close the interpretive gap for coastal boulder deposits, and also provide a framework for coastal risk analysis.
Previous attempts to calculate coastal wave heights have employed a hindcasting method, using masses of emplaced boulders to calculate the minimum flow velocity and thence wave height needed to transport them (Nott, 2003b;Pignatelli et al., 2009;Barbano et al., 2010;Benner et al., 2010;Engel and May, 2012;Mottershead et al., 2014). This approach however, embeds preconceptions about onshore flow hydrodynamics, relying on the presumption of a fixed Froude number depending on wave type (Kennedy et al., 2016;Cox et al., 2020).
Here, we take a different tack. Examining factors that contribute to wave height amplification and onshore bore velocity, we estimate maximal onshore wave impacts that might occur, given offshore significant wave height H s . These estimates can be compared to existing boulder masses, but do not depend on them; and further, can be applied to structural design specifications at any coastline. Extreme waves and flow velocities are low-probability events, but appraising their impacts provides an additional lens for examining coastal boulder deposits, and also long term potential for storm damage at coasts.

ESTIMATING COASTAL WAVE-HEIGHT MAXIMA
Approaching the shore, waves interact with other waves, sea-bottom topography, and coastal geometry. Many waves are dampened, others amplified-sometimes tremendously (Herterich and Dias, 2019). The highest possible coastal wavecrest elevation will occur if the largest wave in a sea undergoes maximum amplification during shoaling and runup. Clearly, this will be a rare occurrence, with very low probability in any given storm. Over many storms, however the likelihood increases. Thus, in long-term evaluations of coastal change or coastal risk, maximal events should not be ignored.
Coastal wave heights depend on complex interplay of surface water with wind gust and fetch, as well as interactions with bathymetry and topography. The maximum offshore wave height H max can be approximated by 2H s (e.g., Krogstad, 1985;Cattrell et al., 2018), but seas evolve as they transit the shelf and shoal toward land, so that H max obtained from buoy data may not match maximum wave height at the coast. Wave energy attenuates across the continental shelf, typically causing large swell height to decrease by about half (Ardhuin et al., 2003;Janjić et al., 2018), but amplification commonly occurs near shore, as shoaling waves interact with sea bed and coast (Viotti and Dias, 2014;Akrish et al., 2016;Brennan et al., 2017). Predicting coastal wave heights is therefore complicated.
Here we consider only pre-impact unbroken waves, and horizontal inland bore velocity. We ignore storm surge, tide level, and infragravity effects, which for this analysis represent short-term sea-level increases. In particular too, we ignore jetting, by which waves impacting cliffs can generate spray heights far greater than the incident wave height (e.g., Peregrine, 2003;Lugni et al., 2006). Although important for erosion (e.g., Bredmose et al., 2015), and for quarrying cliff-top boulders (e.g., Herterich et al., 2018), most of that force is directed vertically. Part of the jet may collapse on the cliff-top platform, but with relatively little inshore transmission of mass or energy. Horizontal flows generated by breaking or overtopping waves are the main agents of onshore sediment movement and structural damage; and our analysis therefore focuses on waves likely to produce them.

Factors Contributing to Wave Height at the Coast
Processes affecting incoming waves can be given as variables in an expression relating offshore H s to the crest elevation of the largest wave impacting the coast η max : C max is the multiplier relating H max to H s (H max = C max H s ), and has a value ≥ 1 (H max being either measured or statistically inferred). C atten describes the proportional attenuation as waves cross the continental shelf, and thus has a maximum value of 1. C amp describes nearshore amplification, as shoaling waves interact with bathymetry. Because η max is the sea surface elevation above still-water level (i.e. not crest-trough height), then in the case where there is no amplification relative to H s , C amp = 0.5 (i.e., η max is given by local amplitude, or half the inshore wave height). Values > 0.5 denote increasing amplification. C amp = 1 would be the case where sea surface elevation is at double the amplitude of the incoming wave before amplification. We discuss and quantify each of these variables below. Theoretically, as outlined above, C max ∼ 2 (because H max ∼ 2H s ); and C atten ∼ 0.5. At first blush, therefore, the two factors cancel, suggesting that offshore H s may be a reasonable predictor of the upper end of the coastal wave height spectrum. But we emphasise that this will vary. If, for example, H s data come from an inshore buoy, or the continental shelf is narrow, there will be little attenuation, so for our illustration we choose C max = 2 and C atten = 1.
Choosing a value for C amp is difficult, because complexities abound near the coast (Hansom et al., 2008;Didenkulova and Pelinovsky, 2011;Slunyaev et al., 2011;Carbone et al., 2013). Reflection can double wave height by simple constructive interference (C amp = 1); but non-linear interactions, and complex interplay with irregular bathymetry and other waves, FIGURE 1 | Strongly amplified wave overtopping the harbour wall at Kalk Bay, Cape Town, South Africa, on June 3rd, 2015. The greatest H s offshore of Cape Town on that day was ∼ 4-5 m (South African Weather Service, document reference: ops-fcast-CT-083). The lighthouse is 6 m tall; the continuous green-water wall is almost 12 m high, i.e., more than twice the offshore H s . Photo copyright Rob Tarr (https://www.sapeople.com/2015/07/11/massive-waves-south-africa-kalk-bayharbour-waves/), used with permission.
Here, we adopt a conservative value of C amp = 1, based on numerous measurements of inshore waves > 2H s (i.e., "rogue waves, " e.g., Kharif and Pelinovsky, 2003;Gemmrich and Thomson, 2017). Photographs document dramatic amplifications for which there are no direct measurements, but which are clearly exceptional waves of substantial size (e.g., Figure 1). In some places and at some times-e.g., locations with unusual topography, providing resonances with certain wave spectra (e.g., Brennan et al., 2017;Herterich and Dias, 2019)-C amp may be > 1 (e.g., Didenkulova et al., 2006;Didenkulova and Anderson, 2010). At most sites, however, it should be less. And we emphasise that users should pick coefficient values appropriate for individual sites.

Overtopping Wave Height: Taking Coastal Elevation Into Account
At sea level, η max fully describes impinging wave amplitude. At steep or cliffed coasts, however, only the overtopping portion of the wave will propagate inland (Figure 2). Therefore, coastal elevation Z must be subtracted from η max to calculate overtopping wave height h O : On low-lying shores, Z ∼ 0 so that h O is essentially η max , and the full height of the wave comes ashore to generate the overland bore. Where there is a cliff or a seawall, Z is the elevation of that barrier. If the coast is steeply sloping but not vertical, the simplest approach is to use the elevation at the point of interest: e.g., in the case of boulder movement or infrastructure impact, Z is the elevation at that location. We do not apply any correction for energy dissipation as water moves inland and against gravity, but such calculations exist (e.g., Noormets et al., 2004;Ogawa et al., 2015;Poate et al., 2018), and can be incorporated into these equations.

ESTIMATING MAXIMUM BORE VELOCITY ONSHORE
Onshore flow takes many forms, including swash uprush following breaking and dryland bores formed by collapse of overtopping waves (Figure 2). The processes are complex, FIGURE 2 | Highly stylised diagram of wave amplification and flow acceleration at the coast. This is not intended as a realistic depiction of a wave storm sea, but simply illustrates components discussed in the text. Offshore wave spectra can include large individual waves (H max ). Interaction with bathymetry and with steep coastal topography may further amplify the largest waves, resulting in an extreme surface elevation at the coast η max , with the elevation above the land surface given by h O . Collapse on the coastal platform or cliff top yields a bore with velocity U bore , which may be focussed and accelerated by topography, generating high-energy flows (U max ) with the capability of moving large boulders.
multiphase, and difficult to model, but field and lab observations show that storm waves can generate flows hydrodynamically similar to tsunami bores (Kennedy et al., 2016;Wüthrich et al., 2018). Dam-break theory (Ritter, 1892) relates horizontal velocity of the bore front U bore to wave height H: For a horizontal dry bed, analytical solutions give C dam = 2 (Ritter, 1892;Chanson, 2006), but experiments and measurements show that the value actually varies considerably in natural settings (Chanson, 2006;Hansom et al., 2008;Hu et al., 2015). This makes sense, because hydrodynamically, C dam is simply the Froude number Fr 1 . Wave-generated flows can have Fr > 2 (Holland et al., 1991;Hansom et al., 2008) and even Fr ≥ 4 (Matias et al., 2016). Likewise, tsunami flows with lower Fr (< 2) are also common (e.g., Yeh, 1991;Imamura et al., 2008;Matsutomi and Okamoto, 2010). We therefore cast C dam as a variable (> 0), with value depending on local topographic and roughness conditions. We can restate the Ritter equation in terms of Equation (2): where U bore is onshore flow velocity, and h O can be thought of as the instantaneous depth of water above the shore platform (e.g., Hu et al., 2015).
where D is the depth of the flow.

Factors that Amplify Flow Velocity
Variability in real-world bore velocities is substantially greater than predicted by classic dam-break theory. Lab and field measurements in many environments-from overwash of sandy and gravel barriers to greenwater overtopping of platforms-have yielded C dam -equivalent values ranging from 1.2 to 10 (Holland et al., 1991;Cox and Ortega, 2002;Ryu et al., 2007;Donnelly, 2008;Chang et al., 2011;Matias et al., 2014Matias et al., , 2016Chuang et al., 2015;Song et al., 2015;Wüthrich et al., 2018). These coefficients are both higher and lower than the theoretical Ritter (1892) value of 2. Environmental interactions clearly regulate flow behaviour, sometimes weakening but often accelerating it. A major mechanism for flow-front acceleration, commonly overlooked, is topographic focusing after wave-breaking and during bore propagation (Denny et al., 2003). This funneling and acceleration can be expressed by a parameter C topo (≥ 1), yielding an expression for maximum onshore velocity: Which can be restated as: C topo is difficult to estimate, and data for surf-zone flow velocities under high-energy conditions are rare. But the few measurements available indicate that strong flow magnifications do occur during runup on uneven rocky coasts. For example, Yeh et al. (1994) reported that convergence of refracted wave fronts produced a three-fold increase in tsunami bore runup velocity. Denny et al. (2003) showed that collision of refracted bore fronts, channeled and focused by topography, can double flow speeds, i.e., C topo = 2. Similarly, Carrasco et al. (2014), studying a breakwater with large roughness elements, reported overtopping flows with U max of 11 m/s, related to 5 m waves. Plugged into Equation (6) and assuming C dam = 2, this yields C topo ∼ 4.5. These are only two studies, and without more data it's hard to know how representative they are; but whereas the Carrasco et al. (2014) study was short term and at one location, the Denny et al. (2003) study encompassed a 1.5-year time series of field data from 221 sites, backed up with wave tank experiments, both approaches yielding velocity-amplification factors ∼ 2. Recognising that Denny et al. (2003) may not have captured the maximum possible amplification, and pointing out that more studies of this phenomenon are needed, we choose C topo = 2 as a reasonable preliminary value.

A UNIFIED EQUATION RELATING OFFSHORE WAVE HEIGHT AND ONSHORE FLOW VELOCITY
Combining Equations (2) and (6), we can relate U max to H s to give a theoretical maximum velocity for wave-induced onshore flow: With coefficient values as chosen above (C topo = 2, C dam = 2, C max = 2, C atten = 1, C amp = 1) and for locations at sea level (Z = 0), Equation (7) becomes: Equation (8) means that for a given H s -and if all amplification factors are maximised-resultant extreme onshore velocity could be almost six times that of a shallow-water wave with speed gH s . Such values deliberately push the envelope of what is possible and we emphasise that velocities this high relative to the generating wave heights would be exceedingly rare.

COMPARISON WITH PREVIOUS APPROACHES
Previous studies calculating onshore flow velocities have generally used them to hindcast nominal heights for the boregenerating waves. The most widely applied approach in coastal boulder studies is that of Nott (1997Nott ( , 2003b, which has driven thinking about likely heights and transport capacities of coastal wave flows (see Kennedy et al., 2019, their Table 1), so we address it in detail here. Attempting to discriminate between tsunami and storm waves as mechanisms of boulder transport, Nott (1997Nott ( , 2003b introduced hydrodynamic equations linking clast characteristics to threshold horizontal flow velocity U bore . Starting with boulder dimensions and density, he derived by force balance the flow velocity necessary for transport. Once this velocity is calculated, the dam-break expression (Equation 3) can be reorganised and used to solve for wave height H-but only if the value of C dam is known. Nott (2003b) assumed that the Ritter coefficient C dam = 2 would apply to all tsunami flows, and that flows generated by storm waves would have the velocity of a shallow water wave (i.e., C dam = 1). He therefore defined a 'wave parameter' δ (equivalent to the square of the coefficient in the Ritter equation, or C 2 dam in the terminology of Equation 3). Thus, Nott's (2003b) formulation for wave height is: Nott's presumption of specific δ values (1 for storm waves and 4 for tsunami) opened a door to using boulder measurements to calculate nominal wave heights required for transport, with modifications for different pre-transport settings (i.e., whether the boulder was initially submerged, sitting on dry land, or still part of bedrock). Subsequent studies incorporated Nott's (1997;2003b) approach in a range of equations attempting to relate flow speeds to wave heights (see Table 1). Such equations are now often employed to hindcast wave heights for both modern and ancient coastal boulder deposits. In particular, they have frequently been used to argue that coastal boulders could not have been emplaced by storm waves, because the equations call for waves too large relative to measured or inferred H s (e.g., Kennedy et al., 2007;Mastronuzzi et al., 2007;Scicchitano et al., 2007;Barbano et al., 2010;Boulton and Whitworth, 2017;Roig-Munar et al., 2019b). Flaws in the Nott equations have been detailed by many (see Cox et al., 2020 and references therein). A key problem is that Nott's (2003) wave parameter δ is actually equivalent to Fr 2 (Kennedy et al., 2016;Cox et al., 2020): Nott's approach requires Fr = 1 for storm waves and Fr = 2 for tsunami, unvarying. But given the complexity of wave-coast interactions, and the non-linearity of wave breaking and overtopping (Cooker and Peregrine, 1992;Carbone et al., 2013), there is no basis for assuming that storm-wave flows have lower Fr than tsunami. Nott's δ is not a constant, but a variable, which will differ from one wave-generated flow to another. Therefore, Equation's (9) central premise fails. In particular Nott's (1997;2003b) assumption that storm-wave bore velocity never exceeds shallow-water wave speed-hard-wired into the equations in Table 1-neglects local effects (and is contradicted by many observations, discussed previously). Thus, coastal storm-wave heights and flow velocities calculated by equations based on the Nott assumptions are not meaningful.
In contrast, the relationships laid out in Equation (7) are explicitly variable, and make no assumptions about hydrodynamic state. No coefficient values are stipulated. The equation is designed to be flexible and adaptable, providing an estimate of the most extreme conditions possible for a given storm and coastal setting.

USES OF THESE EQUATIONS
This analysis indicates that storm waves can generate flows with δ-equivalent values larger than those asserted in the  Benner et al. (2010) Submerged and subaerial The threshold flow speed U bore is calculated from boulder dimensions and then linked to a minimum wave height necessary for boulder transport Min. H through the relation Min. H = U 2 bore /(δg). Equation terms are δ, wave parameter (see text for explanations); a, b, c, long, intermediate, and short boulder dimensions; V, volume; ρs, boulder density; ρw, water density;ü, instantaneous flow acceleration; C d , drag coefficient; C L , lift coefficient; Cm, inertia coefficient; µ, static friction coefficient; Q, empirical coefficient; g, gravity; θ , angle of the bed slope. W, weight of heaviest boulder in study; Z, boulder elevation; U bore , velocity threshold for boulder motion; Min. H, nominal minimum storm wave height required (Min. H = U 2 bore /g), from the cited study; Max. Hs, largest Hs at the boulder location (used as nominal maximum wave height in the cited study); Umax, bore velocity calculated from Equation (7). Storm-wave transport was discounted in these studies because their calculated Min. H was greater than the known Max. Hs. However, applying Equation (7), with boulder elevations and Max. Hs from the cited studies (and using non-maximal coefficients: see text for details), shows that bores with sufficient velocity to move coastal megagravel (Umax ≥ U bore ) can be produced by storm waves consistent with local wave climate. Nott Approach. For example, if H s = 10 m (i.e., conditions seen on sub-decadal timescales on high-energy coasts: e.g., Komar and Allan, 2007;Chen and Curcic, 2016;Masselink et al., 2016) an extreme onshore velocity U max could be as great as 56 m/s according to Equation (8). This is more than 5 times the Nott U bore corresponding value of 10 m/s, and would imply a δequivalent value of 10 (i.e., ten times the Nott-asserted stormwave value for δ).
It's extremely unlikely that all factors of Equation (7) would be maximized for any single wave event, and we do not imply that flows with outsize velocities occur during any given storm. But it is probable that some of these factors act on some large waves on a fairly regular basis, producing local runup velocities much greater than simple dam-break theory would predict. Didenkulova et al. (2006) reported that coastal waves with H ≥ 2H s occur much more frequently than they are observed. We suggest that large coastal wave events and consequent high-velocity onshore flows occur more commonly than predicted from classical theories.
For illustration, we examined coastal boulders that have been interpreted as tsunami deposits via the argument that storm waves had insufficient power to move them. Table 2 lists examples of such boulders from the literature, along with computed threshold velocities and upper limits on available wave energy. We used the published data (boulder dimensions, elevations, and wave climate information) as inputs for Equation (7), but we dialled down its coefficient values: we supposed waves that are 20% attenuated (C atten = 0.8), where H max is 80% greater than H s (C max = 1.8), and amplification occurs in the nearshore environment (C amp = 0.8). We also assumed that the resulting flow behaved like a classic dam break flow (C dam = 2) and could be accelerated 70% by topographic focusing (C topo = 1.7). The point of Table 2 is to show that when amplification factors are considered-even if they are sub-maximal-bores potentially generated by storm waves are more than sufficient to move very large objects. Does this mean that storm waves have moved those boulders? Maybe, although that cannot be proved without additional evidence. The main result, however, is that the possibility of displacement by storm waves cannot be discarded on the basis of hydrodynamic equations alone. Any large storm wave can be rendered exceptional if the amplification factors align constructively.
As another illustration, let us consider the case of the largest boulder known to have been moved by storm waves (620 t, at 3 m above high water: Cox et al., 2018). Bore velocity required to slide this rock is around 11 m/s from equations in Nandasena et al. (2011) (velocity equations of Nott, 2003b return a higher value, ∼ 17 m/s, but do not adequately consider boulder geometry). The maximum H s in the winter that this boulder moved was 14.7 m  but the Nott approach would require wave heights more than twice this. Boulders of this size (as well as smaller ones) are routinely interpreted as tsunami transported. However, per Equation (7), this boulder could have been moved by a wave with a height equivalent to H s = 14.7 m (C atten = 1, C max = 1, C amp = 0.5) and which produced a dam-break-style flow (C dam = 2) that had no topographic focusing (C topo = 1). Moreover, it could also have been moved by a much smaller wave, say ∼ 10 m on approach to the coast (C atten = 0.5, C max = 1.4), that was amplified inshore by 50% (C amp = 0.75), and then collapsed to form a dam-break flow (C dam = 2) that was mildly topographically focused (C topo = 1.5). Many other numerical combinations are possible, and we cannot know which ones may have occurred. But the point is that there are abundant ways to amplify storm waves sufficiently to produce bores strong enough to do extreme hydrodynamic work onshore.
We emphasize that the equations developed here are highly generalised and return indicator values only. They use a simple linear augmentation of height and velocity, whereas non-linear interactions certainly happen: non-linearity could yield extreme dampening, or additional amplification. The coefficients can only be approximate, because it is not possible to better constrain them. Despite recent advances in modelling (e.g., Weiss and Diplas, 2015;Zainali and Weiss, 2015;Zhang et al., 2017), the dynamics of high-energy coasts still are so poorly known that rough guesses are the best we can do. But coupled with reliable hydrodynamic estimates of the velocity required to move a boulder, this approach permits evaluation of wave conditions needed for transport. Given H s , for example, one could use Equation (7) to estimate the combined total of the amplification factors required to match a threshold boulder-entertainment velocity. This would permit initial evaluation of the feasibility or likelihood of storm-wave transportation.
Exceptional events cannot be ignored when considering longterm storm histories, or future projections. Equations (2) and (6) provide a simple way to evaluate the kinds of amplifications that could generate remarkably large waves and accelerated bores. The longer the time period under consideration, the more important it is to consider the possibility that maximal waves might have occurred. In the case of coastal boulder deposits, the timescale is commonly centuries to millennia (e.g., Barbano et al., 2010;Costa et al., 2011;Shah-hosseini et al., 2013;Rixhon et al., 2018), and it is generally decades to centuries in forward-looking risk analyses for coasts. It takes only one extraordinary wave to move a giant boulder, break a seawall, or tear out a building.

IN CONCLUSION: TOWARD A NEW AWARENESS OF EXTREME COASTAL WAVE BEHAVIOUR
Amplification of wave height and onshore flow velocity should be taken into account when trying to determine upper limits of storm energy. Maximal waves with extreme onshore flow velocities will be rare events. But in considering whether boulders could be storm-transported, the possibility of extreme wave events-however statistically unlikely-cannot be discarded, and must be actively considered. This is especially true both for interpreting ancient deposits, and for coastal forward planning.
But onshore damage does not necessarily require the largest possible waves and the maximum possible velocities. Lesser waves-still larger than H s or even H max of a spectrum-may occur with greater frequency than the maximal monsters of Equation (8). The bottom line is that when considering stormdamage probabilities, and wave-coast interactions in general, waves with excess height and excess onshore velocity relative to offshore wave heights must be expected. This analysis is rough, but the roughness is actually part of its value. It shines a light on an important truth: the equations are imprecise by necessity, because basic information is lacking. We hope that by unpacking the processes and exposing how poorly quantified they are, we can spur further research into the upper limits of onshore wave heights and onshore flows. As necessary data are gathered and understanding grows, better constraints and more robust values will be generated. In addition, every location is different, and so people applying this approach to specific sites should carefully consider the best local estimates for these variables. And we emphasise that we are not predicting wave heights and bore velocities for specific storms, but adopting a probabilistic method for evaluating possible extreme wave impacts over long time periods.

DATA AVAILABILITY STATEMENT
All datasets generated for this study are included in the article.

AUTHOR CONTRIBUTIONS
NB wrote the initial draft. NB and RC collaboratively wrote and revised the current paper.