1. Introduction
Estimation of the total water level (TWL) at the shoreline is an important asset for coastal engineers and those involved in coastal zone management and engineering design. For instance, the TWL describes one of the key components in forecasting tools for the assessment of coastal flood risk or storm impact intensity [
1,
2,
3]. The empirical formulae commonly used to design coastal structures, such as sea walls or rubble mound breakwaters, also rely on the determination of the maximum water level [
4]. The calculation of the TWL has thus been the subject of many studies [
5,
6] aiming in particular to improve the estimation of the wave-induced runup [
7,
8,
9,
10,
11,
12], which is one of the primary contributions to the TWL with tide and atmospheric surge.
Wave runup is composed of a mean time component, the wave setup, and a time-varying component, the swash [
13]. The setup depends on an increase in mean sea level at the wave period scale that balances the onshore component of the momentum flux of the waves in the breaking and surf zones [
14]. The swash is composed of a short-wave (SW) or incident wave component, corresponding to high frequency oscillations of the water level in the frequency band between 0.04 and 0.25 Hz, and an infragravity (IG) component corresponding to the contribution of long waves with frequency ranging between 0.002 and 0.04 Hz. Therefore, the accurate determination of the runup contribution to the TWL requires tackling a series of challenges associated with the processes of transformation of short waves from intermediate to shallow waters together with the interaction of bound and free long waves. In addition, the respective contribution of SW and IG waves will depend on the type of beach. In the case of a dissipative beach, the dynamics of the swash zone will be dominated by IG waves, whereas for intermediate to reflective beaches both types of waves will contribute to the TWL at the shoreline [
15].
One type of approach to estimate the runup consists of applying empirical formulations derived either from laboratory data [
16,
17,
18,
19,
20] or field observations [
13,
21,
22,
23,
24]. These formulae have the advantage of providing an estimation of the runup essentially based on the knowledge of offshore bulk wave parameters, such as the significant wave height
the peak wave period
, and the beach geometry as the foreshore beach slope
. This type of approach can easily be implemented into coastal risk forecasting tools based on fast and low cost computations. However, their application to beaches with complex 3D features is usually limited [
7,
10,
23]. Indeed, comparing several empirical formulations Atkinson et al. [
23] showed that the most accurate models give a relative error of
of up to 25%. Thus, it is often necessary to develop site-specific runup formulations [
25,
26], which require a significant measurement effort to cover a wide range of oceanographic conditions at a given site [
25]. Furthermore, it is often hazardous to collect data in natural environments, especially during extremely energetic events [
10,
27], which is probably the reason for the sparse existence of runup data from extreme events. Another limitation of empirical formulae is that they do not provide any information on the physical processes that control the wave-induced water level at the shoreline.
The limits of empirical formulae can be overcome through application of process-based deterministic numerical wave models. For instance, a phase-and-depth resolving model based on the Reynolds-averaged Navier–Stokes equations was recently used to study the sources of runup variability on planar beaches [
28]. However, the application of this type of wave model is mainly intended for in-depth studies of physical processes that control wave transformations and their interactions with coastal structures [
29]. Indeed, the high resolution and long computation time limit their application to real beach configurations. Phase-resolving and depth-integrated models offer a promising alternative. This type of approach allows to account for the main processes of wave transformation in intermediate and shallow waters, including dispersion and nonlinear effects, while requiring an acceptable computation time. For instance, the SWASH model [
30], a widely used nonhydrostatic nonlinear shallow water model, was applied in 1D mode on an urbanized field site [
11] and in 2D mode on a natural open sandy beach [
10] to compute storm-induced runup. The COULWAVE model [
31], a weakly dispersive and fully nonlinear Boussinesq-type model, was used to investigate wave processes in a fringing coral reef environment at two atoll sites in the western tropical Pacific [
32]. The
BOSZ model [
33], a weakly dispersive and weakly nonlinear Boussinesq-type model was used to compute wave setup induced by energetic breaking waves at a fringing reef site in Hawai’i [
34]. The model was incorporated into a full model suite for coastal inundation [
35] and later used for probabilistic mapping of storm-induced coastal inundation under climate change scenarios [
36]. Both studies involved large computational domains with millions of cells. Most of the previously cited studies demonstrate the ability of phase-resolving depth-averaged models to compute the cross-shore sea and swell waves, the IG waves, and wave-induced setup. This type of approach also succeeds in correctly estimating the
exceedance runup value (
), which is usually used as an indicator of storm impact intensity. However, few studies have focused on the detailed computation of time-varying swash dynamics.
Accurate measurements of water level oscillation at the shoreline under real conditions are usually difficult to perform. It requires one to instrument a thin layer of water, usually during energetic wave conditions in a changing environment. Laboratory data can offer the advantage of providing synchronized high temporal and spatial resolution of wave transformations and wave-driven water level oscillations under controlled conditions. Free surface elevations can be measured using resistance or capacitance wave gauges distributed along a cross-shore transect. The runup oscillation on the beach face is usually measured using a long capacitance wire gauge mounted normal to the beach slope at a fixed height above the bottom. For example, the data collected during the GLOBEX project [
37] was used to validate the application of SWASH to compute the runup variability under dissipative conditions corresponding to irregular waves breaking over a gentle slope [
38] for three different incident wave conditions. The relative runup errors ranged between 1% and 11%. Laboratory data of free surface elevation for eight gauges and runup oscillations were used to provide a comprehensive and detailed methodology for sensitivity analysis, calibration, and validation of the SWASH model for its application to the computation of runup oscillations over fringing reefs [
39]. The fully nonlinear and dispersive Boussinesq-type model FUNWAVE-TVD was tested with laboratory data to assess its ability to predict the cross-shore evolution of significant wave heights in the SW and IG frequency bands, and the runup spectrum for irregular waves propagating over a laboratory scale fringing reef [
40]. Detailed data from a laboratory experiment for waves breaking over submerged reef [
41] were also used to validate the computation of runup over a steep-sided coastal structure with a Boussinesq-type model [
42]. As an alternative to commonly used measuring devices such as pressure sensors or runup wires, the use of LiDAR scanners in coastal research is becoming increasingly common in both field [
43,
44,
45] and laboratory [
46] experiments. The use of a LiDAR scanner provides a continuous description of the area of interest, as opposed to point-by-point measurements with previously cited measuring devices. A single instrument is required to cover a relatively large area. Moreover, LiDAR scanners allow for remote measurements, thus providing data in a nonintrusive way. This can be critical for studies of small-scale processes where surface piercing instruments can lead to obstruction or disturbance of the flow.
In the present work, the Boussinesq-type wave model
BOSZ [
33,
34] is compared to an extensive runup laboratory data set based on LiDAR data obtained during the DynaRev large-scale experiment [
46]. The study focuses on the ability of a depth-integrated phase-resolving model to accurately compute both the setup and the contribution of SW and IG waves of the time-varying water elevation at the shoreline in intermediate and reflective beach configurations. Furthermore, this work presents a detailed sensitivity analysis of the computed results to the model settings, including the influence of phase distribution of the incident short waves and the definition of the threshold value for runup determination.
The outline of the paper is as follows. In
Section 2, the model used in this study is described together with the laboratory data set. The comparisons between the model computations and measurements of spectral wave characteristics, free surface elevation time series, and runup components are presented in
Section 3. A discussion of best practices for proper setup of a phase-resolving wave model with the objective to compute runup oscillations over intermediate and reflective beaches is presented in
Section 4. Finally, concluding remarks are drawn in
Section 5.
5. Conclusions
As phase-resolving depth-integrated models are gaining importance for runup studies, precise validations under controlled conditions and recommendations for best practice are required. In this work, a high-quality laboratory data set is used to investigate the capability and sensitivity of the Boussinesq-type model BOSZ for the computation of nearshore wave transformations, including swash processes over intermediate and reflective beaches. The data set includes accurate LiDAR measurements of the free surface elevation in the surf and swash zone as well as shoreline elevation oscillations.
Wave transformations are accurately captured with low cross-shore errors for both the significant wave height and the wave setup . Time series from the numerical model output of shoreline elevation oscillations as well as swash spectra show a satisfying agreement with the laboratory data. The statistical runup quantity is successfully computed with relative errors of less than . The IG swash is well predicted with errors smaller than . The SW swash and shoreline setup are reasonably well predicted with errors of less than .
The discussion evaluates the sensitivity of the results to the model settings for general numerical computations of wave runup by depth-integrated phase-resolving models. Multiple computations with different breaking indices show that in the range of
to
, the numerical results show little variability overall. Some parameters, such as the grid size or the threshold depth defining the edge of the runup tongue, are found to have a significant impact on the results, and thus the performance of the model. A nondimensional parameter is proposed to find the optimal grid size to improve numerical accuracy. A depth threshold of 10 cm, consistent with other studies [
11,
22], is found to be the most appropriate value for systematic comparisons of numerical and laboratory data to prevent small changes in the beach profile from having a disproportional impact. For model/data comparison, a free surface time series is used as the boundary condition for the numerical model, thus providing information of both amplitude and phase angles. Computations with different sets of random phases demonstrate that accurate replication of the laboratory data can only be achieved when the exact phases are not known. Moreover, a significant variability is observed among the runs with different random phases due to the sensitivity of the IG swash to the initial phase distribution. For beaches with high influence of IG energy, i.e., intermediate to dissipative beaches, the variability of the runup can be significant. If the goal is to reproduce particular runup values, which were previously measured in the laboratory or field, the lack of information of the incident wave phases can contribute to substantial uncertainty. It should also be noted that special attention is necessary when data from laboratory experiments such as the one in this study are used. The generation of irregular waves in a laboratory environment essentially relies on the same technique as that which is used in phase-resolving models, i.e., a wave spectrum is decomposed into a series of individual waves. Laboratory data are inevitably subject to the same problem related to the waves’ phases as phase-resolving wave models.
Overall, the phase-resolving depth-integrated BOSZ model shows satisfying capabilities in modeling irregular wave runup on intermediate and reflective beaches. It proves that this type of numerical model can be a powerful tool for coastal risks assessment and hazard mitigation projects. The sensitivity analysis performed provides guidelines on how to utilize the model and, more generally, any phase-resolving depth-integrated model to find the best accuracy at the lowest computational cost and ensure quality results for runup modeling studies.