A relationship to describe the cumulative impact of storm clusters on beach erosion

Estimation of erosion volumes for adequate dry beach buffer zones is commonly estimated on the basis of a single extreme event, such as the 1 in 100 year storm. However, the cumulative impact of several smaller, closely spaced storms can lead to equal, if not more, dry beach loss, but this is often not quantified. Here we use a calibrated model for dune erosion, XBeach, to hindcast the cumulative erosion impact of a series of historical storms that impacted the Gold Coast, Queensland region in 1967. Over a 6-month period, four named cyclones (Dinah, Barbara, Elaine, and Glenda) and three East Coast Lows caused a cumulative erosion volume greater than the predicted 1 in 100 year event. Results presented here show that XBeach was capable of reproducing the measured dry beach erosion volume to within 21% and shoreline retreat to within 10%. The storms were then run in 17 different sequences to determine if sequencing influenced final modeled erosion volumes. It is shown that storm sequencing did not significantly affect the total eroded volumes. However, individual storm volumes were influenced by the antecedent state of the beach (i.e. prior cumulative erosion). Power-law relationships between cumulative energy density (∑ E) and eroded volume (∆V) as well as cumulative wave power ((∑ P)) and eroded volume (∆V) both explained more than 94% of the modeled dry beach erosion for the 1967 storm sequences. When the relationship was compared with observed and modeled erosion volumes for similar beaches but different storm forcing, the inclusion of pre-storm beach swash slope (β_swash) in the parameterization was found to increase the applicability of the power-law relationship over a broader range of conditions.     

Coastal erosion through integrated management: A case of Southern Thailand

This study presents how Thailand applied an integrated approach to tackle erosion problems by using a case study in Nakorn Si Thammarat province. Communities along 36 km of coastline suffered from continual erosion. Community members believed the erosion was a natural phenomenon that was intensified by human actions. Historical erosion rate estimated by overlaying aerial photographs was about 5 m per year, while LITPROF simulations suggested that approximately 5 m of beach dune would be eroded by storm waves. Stakeholders were identified based on power and legitimacy criteria. Their past attempts to mitigate the erosion were analyzed. Conflicts arose from how they selected erosion protection methods. Lessons learnt from previous management failures taught that addressing needs of the stakeholders and consulting them throughout the design process were of importance. Finally, a combination of detached nearshore breakwaters and beach nourishment was the selected protection measure and was welcome by the communities.     

Coastal defence through wave farms

The possibility of using wave farms for coastal defence warrants investigation because wave energy is poised to become a major renewable in many countries over the next decades. The fundamental question in this regard is whether a wave farm can be used to reduce beach erosion under storm conditions. If the answer to this question is positive, then a wave farm can have coastal defence as a subsidiary function, in addition to its primary role of producing carbon-free energy. The objective of this work is to address this question by comparing the response of a beach in the face of a storm in two scenarios: with and without the wave farm. For this comparison a set of ad hoc impact indicators is developed: the bed level impact (BLI), beach face eroded area (FEA), non-dimensional erosion reduction (NER), and mean cumulative eroded area (CEA); and their values are determined by means of two coupled models: a high-resolution wave propagation model (SWAN) and a coastal processes model (XBeach). The study is conducted through a case study: Perranporth Beach (UK). Backed by a well-developed dune system, Perranporth has a bar between − 5 m and − 10 m. The results show that the wave farm reduces the eroded volume by as much as 50% and thus contributes effectively to coastal protection. This synergy between marine renewable energy and coastal defence may well contribute to improving the viability of wave farms through savings in conventional coastal protection.     

The Relative Trough Froude Number for initiation of wave breaking: Theory,experiments and numerical model confirmation

It is important to accurately locate the wave breaking region for the calculation of nearshore hydrodynamics. Energy from breaking waves drives hydrodynamic phenomena such as wave set-up, set-down, wave run-up, longshore currents, rip currents, and nearshore circulation. Numerous studies have been undertaken to describe when and where wave breakings occurs. Recent development of computer resources permits the use of phase-resolving numerical models for the study of wave propagation, transformation, and nearshore hydrodynamics. This requires new types of wave breaking criteria for the numerical model. The Relative Trough Froude Number (RTFN) is a new wave breaking criterion. This model is based on the moving hydraulic jump concept, therefore it satisfies properly posed boundary-value conditions. It has been experimentally proved that a critical RTFN at the initiation of wave breaking is consistent with and without the presence of an opposing current, but previous efforts did not investigate the theory for the critical value. This paper provides a theoretical analysis and a numerical analysis to demonstrate why the RTFN theory works as a wave breaking initiation (trigger) index. The theoretical analysis provides a universal constant for the initiation of wave breaking for all water depths assuming the Miche formula properly describes the wave breaking condition. A subroutine for wave breaking in a numerical model, FUNWAVE was modified to include the RTFN trigger. The numerical model was calibrated with data from wave tank experiments, and it was found that the critical condition is very close to the theoretical number, CTFN = 1.45. A second paper (in preparation) provides details of the theory and experiments for a second criterion for termination of wave breaking. The time scale for the establishment of the breaking region i.e., between the initiation position and termination position, depends upon the additional momentum present under turbulent condition within the breaking wave. This subject is not considered herein.     

Kinematics of breaking tsunami wavefronts: A data set from large scale laboratory experiments

This paper provides an overview of a new large scale laboratory data set on the kinematics of breaking tsunami wavefronts. The aim of the experiments was to provide an open access data set for model testing, calibration and verification, with particular emphasis on fluid kinematics in the wave breaking and run-up (swash) zones. The experiments were performed over a composite slope in the tsunami wave basin at the O. H. Hinsdale Wave Research Laboratory at Oregon State University. Data for ten different wave conditions were collected, including non-breaking and breaking waves, and both shore breaks and fully developed long bores.Surface elevation and fluid kinematics were measured with a closely spaced array of surface piercing wave gauges, non-contact ultrasonic wave gauges and four 3-D side-looking Acoustic Doppler Velocimeters. The array was traversed from the nearshore (depth = 0.2 m) to the middle and upper run-up zone, providing kinematic data at 30 cross-shore locations. Video was also recorded from 4 cameras covering the propagation, breaking and run-up zones. Surface elevation, flow velocities and the wave maker displacement were also recorded to provide offshore boundary conditions.The experiments include conditions with wave heights up to 0.55 m, notional wave periods up to 20 s and run-up lengths of up to 15.2 m on a 1/30 slope. In terms of the slope in the shoaling and breaker zones, the data correspond to Iribarren numbers in the range of 0.26–5.6. Raw, calibrated and processed data are stored with open access within the OSU Tsunami Wave Basin Experiment Notebook, which provides full access to all the wave maker control signals, data, instrument coordinates, and processing and plotting software. This paper serves as an introduction to the data set, demonstrates data quality and provides an initial analysis of some key parameters that govern the impact of tsunami events, including run-up versus offshore wave conditions and nearshore bore height, the maximum inundation depths at the original shoreline position, and the time to maximum inundation depth and flow reversal. Examples of temporal and convective accelerations and turbulent flow components are also presented to illustrate further details of the kinematics.     

Numerical study of turbulence and wave damping induced by vegetation canopies

Vegetation canopies control mean and turbulent flow structure as well as surface wave processes in coastal regions. A non-hydrostatic RANS model based on NHWAVE (Ma et al., 2012) is developed to study turbulent mixing, surface wave attenuation and nearshore circulation induced by vegetation. A nonlinear k − ϵ model accounting for vegetation-induced turbulence production is implemented to study turbulent flow within the vegetation field. The model is calibrated and validated using experimental data from vegetated open channel flow, as well as nonbreaking and breaking random wave propagation in vegetation fields. It is found that the drag-related coefficients in the k − ϵ model C_fk and C_fϵ can greatly affect turbulent flow structure, but seldom change the wave attenuation rate. The bulk drag coefficient C_D is the major parameter controlling surface wave damping by vegetation canopies. Using the empirical formula of Mendez and Losada (2004), the present model provides accurate predictions of vegetation-induced wave energy dissipation. Wave propagation through a finite patch of vegetation in the surf zone is investigated as well. It is found that the presence of a finite patch of vegetation may generate strong pressure-driven nearshore currents, with an onshore mean flow in the unvegetated zone and an offshore return flow in the vegetated zone.     

Offshore wind farm impacts on surface waves and circulation in Eastern Lake Ontario

A coupled wave and hydrodynamic model was applied to the Kingston Basin of eastern Lake Ontario, a region with bathymetric variability due to channels and shoals, to assess the potential impacts on surface waves and wind-driven circulation of an offshore wind farm. The model was used to simulate a series of storm events with time-varying wind forcing and validated against wave, current and water level observations. The wind farm was simulated by adding semi-permeable structures in the surface wave model to represent the turbine monopiles, and by adding an energy loss term to the fluid momentum equations in the hydrodynamic model to represent the added drag of the monopiles on the flow. The results suggest that the wind farm would have a small influence on waves and circulation throughout the wind farm area, with spatial variability due to focussing of wave energy and re-direction of the flow. Overall, the results indicate that the wave height in coastal areas will be minimally affected with changes in significant wave height predicted to be < 3%. Larger changes to the strength of circulation occur inside the wind farm region with localized changes in current magnitude of up to 8 cm s^− 1. The results of this study may help to understand the impacts of future offshore wind farms and other offshore structures in the Great Lakes.     

Quantifying the potential impact of land cover changes due to sea-level rise on storm surge on lower Texas coast bays

In this study we investigated the impacts of potential changes of land cover due to sea-level rise (SLR) on storm surge (i.e., the rise of water above normal sea level, namely mean-sea level and the astronomical tide, caused by hurricane winds and pressure) response inside bays on the lower Texas coast. We applied a hydrodynamic and wave model (ADCIRC + SWAN) forced by hurricane wind and pressure fields to quantify the importance of SLR-induced land cover changes, considering its impacts by changing bottom friction and the transfer of wind momentum to the water column, on the peak surge inside coastal bays. The SLR increments considered, 0.5 m to 2.0 m, significantly impacted the surge response inside the bays. The contribution of land cover changes due to SLR to the surge response, on average, ranged from a mean surge increase of 2% (SLR of 0.5 m) to 15% (SLR of 2.0 m), in addition to the SLR increments. The increase in surge response strongly depended on storm condition, with larger increases for more intense storms, and geographical location. Although land cover changes had little impact on the surge increase for SLR increments lower than 1.0 m, intense storms resulted in surge increase of up to 10% even for SLR below 1.0 m, but in most cases, the geometry changes were the major factor impacting the surge response due to SLR. We also found a strong relationship between changes in bottom friction and the surge response intensification; demonstrating the importance of considering land cover changes in coastal regions that are highly susceptible to SLR when planning for climate change.     

Appraisal of storm beach buffer width for cyclonic waves

The loss of beach sand from berm and dune due to high waves and surge is a universal phenomenon associated with sporadic storm activities. To protect the development in a coastal hazard zone, hard structures or coastal setback have been established in many countries around the world. In this paper, the requirement of a storm beach buffer, being a lesser extent landward comparing with the coastal setback to ensure the safety of infrastructures, is numerically assessed using the SBEACH model for three categories of wave conditions in terms of storm return period, median sand grain size, berm width, and design water level. Two of the key outputs from the numerical calculations, berm retreat and bar formation offshore, are then analysed, as well as beach profile change. After having performed a series of numerical studies on selected large wave tank (LWT) test results with monochromatic waves using SBEACH, we may conclude that: (1) Berm erosion increases and submerged bar develops further offshore as the storm return period increases for beach with a specific sand grain size, or as the sand grain reduces on a beach under the action of identical wave condition; (2) Higher storm waves yield a large bar to form quicker and subsequently cause wave breaking on the bar crest, which can reduce the wave energy and limit the extent of the eroding berm; (3) A larger buffer width is required for a beach comprising small sand grain, in order to effectively absorb storm wave energy; and (4) Empirical relationships can be tentatively proposed to estimate the storm beach buffer width, from the input of wave conditions and sediment grain size. These results would benefit a beach nourishment project for shore protection or design of a recreational beach.     

Video and field observations of wave attenuation in a muddy surf zone

Besides the different scales within which coastal processes manifest their energy, the majority of the world's coastal regions exhibit forms of sediment heterogeneity that are physically significant. One example of a heterogeneous environment is Cassino beach, located at the southernmost part of Brazil, a fine-grain-sized sandy beach where fluid mud sporadically is transported to the nearshore and eventually onto the beach. At this site in 2005, as part of a field experiment, a video system was installed. Three years after the installation, a large mud transgression event took place in February 2008 and had 5 km of extension. In this context, the goal of the present work is to characterize the mud deposition pattern across the surf zone, describing the consequences of mud on nearshore dynamics using remote sensing techniques, beach profiles and suspended matter concentration. The surveyed beach profiles registered the deposition of fluid mud at the inner surf zone with concentrations up to 12 mg/l. The material was deposited close to the shoreline and had a cross-shore width of 100 m during the first deposition day occupying the entire water column. From surf zone time series of pixel intensity, it was possible to detect the attenuation of the surface wave spectra due to the presence of fluid mud. The combination of video techniques and field data allowed one to follow the formation of a double-layer system, where fluid mud overlays the sandy bottom. The video-based system at Cassino demonstrated that remote detection of fluid mud and quantification of its effect on the nearshore dynamics is feasible. The combination of beach profiles, measurements of suspended matter concentration and intensity timestacks allowed the analysis of the short-term evolution of the mud depositional processes.     

The effect of wave–current interactions on the storm surge and inundation in Charleston Harbor during Hurricane Hugo 1989

The effects of wave–current interactions on the storm surge and inundation induced by Hurricane Hugo in and around the Charleston Harbor and its adjacent coastal regions are examined by using a three-dimensional (3-D) wave–current coupled modeling system. The 3-D storm surge and inundation modeling component of the coupled system is based on the Princeton ocean model (POM), whereas the wave modeling component is based on the third-generation wave model, simulating waves nearshore (SWAN). The results indicate that the effects of wave-induced surface, bottom, and radiation stresses can separately or in combination produce significant changes in storm surge and inundation. The effects of waves vary spatially. In some areas, the contribution of waves to peak storm surge during Hurricane Hugo reached as high as 0.76 m which led to substantial changes in the inundation and drying areas simulated by the storm surge model.     

STWAVE simulation of Hurricane Ike: Model results and comparison to data

Hurricane Ike (2008), with its associated storm surge, caused extensive damage across parts of the northwestern Gulf Coast when it made landfall in the late hours of September 12, 2008 along the upper Texas coast at the upper end of Category 2 intensity. An extensive instrumentation effort allowed the collection of both nearshore and inland wave and water level data as Hurricane Ike passed the Louisiana coast and made landfall in Texas. This article presents the results of a validation effort for the STWAVE model and the bottom friction coefficients applied in the model with comparisons to the Hurricane Ike measured wave data. Examination of STWAVE model results as contour plots and time series of wave height and period; wave spectra at selected time steps and scatter plots of simulated versus modeled wave results allow evaluation of the model performance. STWAVE model results indicate good agreement with the measured nearshore wave data for an open water Manning ‘n’ bottom friction coefficient equal to 0.03 s/m^0.33. STWAVE model results indicate good agreement with the measured inshore wave data with Manning ‘n’ bottom friction coefficients equal to values derived from land classification data and applied in the ADCIRC model.     

Two-dimensional wave-current interaction with a locally enriched quadtree grid system

Wave-induced currents may drive nearshore mixing and transport processes, including coastal pollutant dispersion, littoral drift, and long-term morphological changes through beach erosion and accretion. In this study, a numerical model is newly developed to simulate wave climate and localized currents in complicated coastal environments. The model developed is based on a quadtree grid system. The two-dimensional hydrodynamic governing equations are solved by using an explicit Adams-Bashforth finite difference scheme. Effects of wave breaking, shoaling, refraction, diffraction, wave-current interaction, set-up and set-down, turbulent mixing, bed friction, and shoreline movement are incorporated in the model. Results are presented for set-up at a beach in a flume due to normally incident waves, and longshore currents generated by oblique waves on a plane beach.     

Sub-weekly to interannual variability of a high-energy shoreline

Sixty-one Global Positioning System (GPS), sub-aerial beach surveys were completed at 7 km long Ocean Beach, San Francisco, CA (USA), between April 2004 and March 2009. The five-year time series contains over 1 million beach elevation measurements and documents detailed changes in beach morphology over a variety of spatial, temporal, and physical forcing scales. Results show that seasonal processes dominate at Ocean Beach, with the seasonal increase and decrease in wave height being the primary driver of shoreline change. Storm events, while capable of causing large short-term changes in the shoreline, did not singularly account for a large percentage of the overall observed change. Empirical orthogonal function (EOF) analysis shows that the first two modes account for approximately three-quarters of the variance in the data set and are represented by the seasonal onshore/offshore movement of sediment (60%) and the multi-year trend of shoreline rotation (14%). The longer-term trend of shoreline rotation appears to be related to larger-scale bathymetric change. An EOF-based decomposition technique is developed that is capable of estimating the shoreline position to within one standard deviation of the range of shoreline positions observed at most locations along the beach. The foundation of the model is the observed relationship between the temporal amplitudes of the first EOF mode and seasonally-averaged offshore wave height as well as the linear trend of shoreline rotation. This technique, while not truly predictive because of the requirement of real-time wave data, is useful because it can predict shoreline position to within reasonable confidence given the absence of field data once the model is developed at a particular site.     

Numerical Modeling of the Undertow Structure and Sandbar Migration in the Surfzone

A process-based 3D numerical model for surfzone hydrodynamics and beach evolution was established. Comparisons between the experimental data and model results proved that the model could effectively describe the hydrodynamics, sediment transport feature and sandbar migration process in the surfzone with satisfactory precision. A series of numerical simulations on the wave breaking and shoaling up to a barred beach were carried out based on the model system. Analyzed from the model results, the wave-induced current system in the surfzone consists of two major processes, which are the phase-averaged undertow caused by wave breaking and the net drift caused by both of the nonlinear wave motion and surface roller effect. When storm waves come to the barred beach, the strong offshore undertow along the beach suppresses the onshore net drift, making the initial sandbar migrate to the seaside. Under the condition of calm wave environment, both the undertow and net drift flow to the shoreline at the offshore side of the sandbar, and then push the initial sandbar to the shoreline. The consideration of surface roller has significant impact on the modeling results of the sandbar migration. As the roller transfer rate increases, the sandbar moves onshore especially under the storm wave condition.     

The impact of sea level rise and climate change on inshore wave climate: A case study for East Anglia (UK)

In coastal areas, offshore wave propagation towards the shore is influenced by water depth variations, due to sea bed bathymetry, tides and surges. Considering implications of climate change both on atmospheric forcing and sea level rise, a simple methodology involving numerical modelling is implemented to compute inshore waves from 1960 to 2099. Simulations take into account five scenarios of linear sea level rise and one climatic scenario for storm surges and offshore waves. The methodology is applied to the East Anglia coast (UK). Extreme event analysis is performed to estimate climate change implication on inshore waves and the occurrence of extreme events. It is shown, for this coastal region, that wave statistics are sensitive to the trend in sea level rise, and that the climate change scenario leads to a significant increase of extreme wave heights in the northern part of the domain. For nearshore points, the increase of the mean sea level alters not only extreme wave heights but also the frequency of occurrence of extreme wave conditions.     

Estimation of infragravity waves at intermediate water depth

The accuracy of nearshore infragravity wave height model predictions has been investigated using a combination of the spectral short wave evolution model SWAN and a linear 1D SurfBeat model (IDSB). Data recorded by a wave rider located approximately 3.5 km from the coast at 18 m water depth have been used to construct the short wave frequency-directional spectra that are subsequently translated to approximately 8 m water depth with the third generation short wave model SWAN. Next the SWAN-computed frequency-directional spectra are used as input for IDSB to compute the infragravity response in the 0.01 Hz–0.05 Hz frequency range, generated by the transformation of the grouped short waves through the surf zone including bound long waves, leaky waves and edge waves at this depth. Comparison of the computed and measured infragravity waves in 8 m water depth shows an average skill of approximately 80%. Using data from a directional buoy located approximately 70 km offshore as input for the SWAN model results in an average infragravity prediction skill of 47%. This difference in skill is in a large part related to the under prediction of the short wave directional spreading by SWAN. Accounting for the spreading mismatch increases the skill to 70%. Directional analyses of the infragravity waves shows that outgoing infragravity wave heights at 8 m depth are generally over predicted during storm conditions suggesting that dissipation mechanisms in addition to bottom friction such as non-linear energy transfer and long wave breaking may be important. Provided that the infragravity wave reflection at the beach is close to unity and tidal water level modulations are modest, a relatively small computational effort allows for the generation of long-term infragravity data sets at intermediate water depths. These data can subsequently be analyzed to establish infragravity wave height design criteria for engineering facilities exposed to the open ocean, such as nearshore tanker offloading terminals at coastal locations.     

Experimental simulation of shoreface nourishments under storm events: A morphological,hydrodynamic, and sediment grain size analysis

This study focuses on barred beach shoreface nourishments physically simulated in a wave flume. The attack of a schematic storm on three different nourishments is analysed. The apex and waning storm phases lead respectively to offshore and onshore sediment transports. Nourishments in the trough and on the outer bar feed the bar and increase wave dissipation offshore. The bar acts as a wave filter and reduces shore erosion (lee effect). In contrast, nourishment on the beach face leads mostly to shore feeding and reconstruction (feeder effect). With successive nourishments, the beach face clearly becomes steeper and onshore sediment transport is reduced during moderate wave climates. The surface grain size analysis reveals marked variations. Coarser sediments are sorted on the bar and the upper beach face. These locations correspond to large wave dissipation zones during the storm apex.     

One-year along-beach variation in the maximum depth of erosion resulting from irregular shoreline morphology

Beach erosion and accretion occur across multiple time scales. Over long time scales (decades to millennia) the shoreface ravinement surface, which is recognized as a coarse lag deposit, forms at the shoreface toe as a result of wave- and current-induced erosion during shoreline transgression. Over short time scales (hours to days) the depth of sediment disturbance, which is recognized as coarse lamina and measured at the foreshore by devices and monitoring tracer beds, forms as a result of wave- and current-induced reworking during a tidal cycle. The maximum depth of erosion (MDOE), quantified here over 1 year, is modulated by processes that operate over a time scale that is between the drivers of short-term (e.g. tides and waves) and long-term (e.g. sea-level rise) beach erosion. The MDOE integrates the erosion that occurs over a discrete time interval and records the maximum depth of erosion that is likely principally induced by storms, which is difficult to quantify by other methods that rely on discrete observations (e.g. changes in elevation or movement of the mean high-water line). A novel technique for quantifying the MDOE, based on comparing the bedding and stratigraphy between cores collected at the same locations over a discrete time interval, is presented here and applied at Onslow Beach, NC, USA. This 12 km-long barrier island has irregular shoreline morphology, characterized by two embayments separated by a central headland. This shape is largely the result of variations in the depth of underlying rock strata and produces a steeper beachface at the headland than at the embayments. At each of the six sites examined along the barrier, the MDOE is found to increase from the backshore to the middle intertidal zone and is higher at the sites closer to the headland. These variations in the MDOE are likely due to the increase in average wave energy impacting the beachface from an offshore direction and steeper beaches (intermediate beach state) at the headland. Where the MDOE is within the beach facies, it is not associated with a coarsening, which is due to the heterolithic nature of the Onslow-beach strata. Where the MDOE is the contact between back-barrier and beach facies it is always associated with a coarsening and a gravel-rich lag deposit because in this case, the MDOE is an amalgamation of multiple erosional events, which is similar to the shoreface ravinement surface. Along-beach variation in the MDOE does not correspond with discrete observations of beach change over the same period and is likely a better indicator of erosion potential than long-term discrete observations, such as changes in surface elevation or the position of the mean high-water line.     

Modelling storm hydrodynamics on gravel beaches with XBeach-G

In this paper we present a process-based numerical model for the prediction of storm hydrodynamics and hydrology on gravel beaches. The model comprises an extension of an existing open-source storm-impact model for sandy coasts (XBeach), through the application of (1) a non-hydrostatic pressure correction term that allows wave-by-wave modelling of the surface elevation and depth-averaged flow, and (2) a groundwater model that allows infiltration and exfiltration through the permeable gravel bed to be simulated, and is referred to as XBeach-G. Although the model contains validated sediment transport relations for sandy environments, transport relations for gravel in the model are currently under development and unvalidated. Consequently, all simulations in this paper are carried out without morphodynamic feedback. Modelled hydrodynamics are validated using data collected during a large-scale physical model experiment and detailed in-situ field data collected at Loe Bar, Cornwall, UK, as well as remote-sensed data collected at four gravel beach locations along the UK coast during the 2012–2013 storm season. Validation results show that the model has good skill in predicting wave transformation (overall SCI 0.14–0.21), run-up levels (SCI < 0.12; median error < 10%) and initial wave overtopping (85–90% prediction rate at barrier crest), indicating that the model can be applied to estimate potential storm impact on gravel beaches. The inclusion of the non-hydrostatic pressure correction term and groundwater model is shown to significantly improve the prediction and evolution of overtopping events.