Measuring run-up on a natural beach

Field experiments have been performed to evaluate and intercompare two techniques for measuring run-up on natural beaches, resistance wires and films. Simultaneous deployment of wire sensors shows a low error (< 5%) in electronics gain, but a strong sensitivity to the elevation of the wires above the beach face. On a low slope (β ～ 0.02) beach, with incident wind waves of moderate height (H ～ 1 m), differences of only a few cm in the wire elevation cause variance differences as large as 25%, in otherwise identical sensors. Replicate digitizations of the same run-up film show variance differences as large as 20%, with an average deviation from the mean variance of 8%.Use of the film and resistance wire sensors on the same run-up field showed small differences in the mean swash elevation (i.e., set-up), but an 83% difference in swash variance. Much further work is needed to determine the dependence of sensor differences on beach slope, porosity, camera elevation and other factors.     

Empirical parameterization of setup,swash, and runup

Using shoreline water-level time series collected during 10 dynamically diverse field experiments, an empirical parameterization for extreme runup, defined by the 2% exceedence value, has been developed for use on natural beaches over a wide range of conditions. Runup, the height of discrete water-level maxima, depends on two dynamically different processes; time-averaged wave setup and total swash excursion, each of which is parameterized separately. Setup at the shoreline was best parameterized using a dimensional form of the more common Iribarren-based setup expression that includes foreshore beach slope, offshore wave height, and deep-water wavelength. Significant swash can be decomposed into the incident and infragravity frequency bands. Incident swash is also best parameterized using a dimensional form of the Iribarren-based expression. Infragravity swash is best modeled dimensionally using offshore wave height and wavelength and shows no statistically significant linear dependence on either foreshore or surf-zone slope. On infragravity-dominated dissipative beaches, the magnitudes of both setup and swash, modeling both incident and infragravity frequency components together, are dependent only on offshore wave height and wavelength. Statistics of predicted runup averaged over all sites indicate a − 17 cm bias and an rms error of 38 cm: the mean observed runup elevation for all experiments was 144 cm. On intermediate and reflective beaches with complex foreshore topography, the use of an alongshore-averaged beach slope in practical applications of the runup parameterization may result in a relative runup error equal to 51% of the fractional variability between the measured and the averaged slope.     

Probability distributions for wave runup on beaches

Measured probability distributions of shoreline elevation, swash height (shoreline excursion length) and swash maxima and minima from a wide range of beach types are compared to theoretical probability distributions. The theoretical distributions are based on assumptions that the time series are weakly steady-state, ergodic and a linear sum of random variables. Despite the swash process being inherently non-linear, results indicate that these assumptions are not overly restrictive with respect to modeling exceedence statistics in the upper tail of the probability distribution. The RMS-errors for a range of exceedence level statistics (50, 10, 5, 2, and 1%) were restricted to < 10 cm (and often < 5 cm) for all of the swash variables that were investigated. The results presented here provide the basis for further refinement of coastal inundation modeling as well as stochastic-type morphodynamic modeling of beach response to waves. Further work is required, however, to relate the parameters of swash probability distributions to wave conditions further offshore.     

礁面大糙率存在下孤立波传播变形及爬高实验研究

本文采用圆柱体阵列来模拟珊瑚礁面的大糙率,通过波浪水槽实验研究礁面糙率对孤立波传播变形及岸滩爬高的影响.结果表明,粗糙礁面的存在显著削弱了礁坪上孤立波的首峰和礁后岸滩反射造成的次峰,同时降低了波浪在珊瑚礁面的传播速度;垂直于岸线方向沿礁相对波高随着入射波增大而减小,随着礁坪水深的增大而增大,粗糙礁面上波高沿礁的衰减更为...     

Geometric characteristics of storm-beach profiles caused by inclined waves

A theoretical analysis shows that the geometric characteristics of a storm-beach profile is governed by a modified Iribarren number which includes the effects among the factors of beach slope, breaking wave angle and wave steepness. A series of experiments have been conducted in a three-dimensional movable bed model on the conditions of two different beach slopes and two incident wave angles as well as several erosive wave steepnesses. Based on the experimental data, the relative importance of each factor involved in the parameter is discussed. The empirical relationships between the geometric characteristics of a storm-beach profile and the modified Iribarren number are proposed through regression analysis.     

On low-frequency waves in the surf and swash

Low-frequency waves in the surf and swash zones on various beach slopes are discussed using numerical simulations. Simulated surface elevations of both primary waves and low-frequency waves across the surf zone were first compared with experimental data and good agreement found. Low-frequency wave characteristics are then discussed in terms of their physical nature and their relationship to the primary wave field on a series of sea bottom slopes. Unlike primary waves, low-frequency wave energy increases towards the shoreline. Low-frequency waves in the surf and swash are a function of incident waves and the sea bottom slope and hence the saturation level of the surf zone. Wave energy on a gently sloping beach is dominated by low-frequency waves while primary waves play a significant role on a steep beach. Low-frequency wave radiation from the surf zone on a given beach depends on primary wave frequency and beach slope. However, a very poor correlation was found between surf similarity parameter and low-frequency wave radiation.     

Simulation and prediction of swash oscillations on a steep beach

This paper presents numerical simulations and analytical predictions of key aspects of swash oscillations on a steep beach. Simulations of the shoreline displacement based on bore run-up theory are found to give excellent agreement with recent experimental data for regular waves, wave groups and random waves. The theory is used to derive parameters that predict the onset of swash saturation and the spectral characteristics of the saturated shoreline motion. These parameters are again in good agreement with the measured laboratory data and are also consistent with previous experimental data. Simulation of irregular wave run-up using a series of overlapping monochromatic swash events is found to reproduce typical features of swash oscillations and can accurately describe both the low and high frequency spectral characteristics of the swash zone. In particular, the low frequency components of the run-up can be modelled directly using a sequence of incident short wave bores, with no direct long wave input to the numerical simulations. This suggests that wave groupiness must be accounted for when modelling shoreline oscillations.     

Lagrangian measurements and modelling of fluid advection in the inner surf and swash zones

New laboratory and field data are presented on fluid advection into the swash zone. The data illustrate the region of the inner surf zone from which sediment can be directly advected into the swash zone during a single uprush, which is termed the advection length. Experiments were conducted by particle tracking in a Lagrangian reference frame, and were performed for monochromatic breaking waves, solitary bores, non-breaking solitary waves and field conditions. The advection length is normalised by the run-up length to give an advection ratio, A, and different advection ratios are identified on the basis of the experimental data. The data show that fluid enters the swash zone from a region of the inner surf zone that can extend a distance seaward of the bore collapse location that is approximately equal to half of the run-up length. This region is about eight times wider than the region predicted by the classical swash solution of Shen and Meyer [Shen, M.C., Meyer, R.E., 1963. Climb of a bore on a beach. Part 3. Runup. Journal of Fluid Mechanics 16, 113–125], as illustrated by Pritchard and Hogg [Pritchard, D., Hogg, A.J., 2005. On the transport of suspended sediment by a swash event on a plane beach. Coastal Engineering 52, 1–23]. Measured advection ratios for periodic waves show no significant trend with Iribarren number, consistent with self-similarity in typical swash flows. The data are compared to recent characteristic solutions of the non-linear shallow water wave (NLSW) equations and both finite difference and finite volume solutions of the NLSW equations.     

RANS v2–f simulation of a swash event: Detailed flow structure

The flow structure of a swash event over a uniform slope is studied using a RANS-VOF numerical model coupled with a v²–f turbulence closure. The model is compared with experimental data of recent laboratory experiments. The ability of the turbulence modelling for simulating swash flow and the evolution of the computed bed shear stress during run-up and run-down are investigated. The agreement between numerical results and measured data, such as water depth, depth-averaged velocity and bed shear stress is very good during run-up. Main discrepancies are found during run-down. The paper also examines the aeration of the water layer in the swash flow, taking advantage of the PLIC method for computing the air–water interfaces. Air is continuously entrapped in the swash front and released at its rear during run-up. A detailed analysis indicates that the flow reversal is initiated near the bottom at the outer boundary of the swash zone and progresses landward. The study highlights the asymmetry between run-up and run-down. During run-up, the swash front propagation determines the turbulence properties and the bed shear stress profile on the beach, whereas the flow properties are more homogeneously distributed in the swash area during run-down.     

Run-up of solitary waves on slopes with different profiles

The problem of sea-wave run-up on a beach is discussed within the framework of exact solutions of a nonlinear theory of shallow water. Previously, the run-up of solitary waves with different forms (Gaussian and Lorentzian pulses, a soliton, special-form pulses) has already been considered in the literature within the framework of the same theory. Depending on the form of the incident wave, different formulas were obtained for the height of wave run-up on a beach. A new point of this study is the proof of the universality of the formula for the maximum height of run-up of a solitary wave on a beach for the corresponding physical choice of the determining parameters of the incident wave, so that the effect of difference in form is eliminated. As a result, an analytical formula suitable for applications, in particular, in problems related to tsunamis, has been proposed for the height of run-up of a solitary wave on a beach.     

Swash overtopping and sediment overwash on a truncated beach

New experimental laboratory data are presented on swash overtopping and sediment overwash on a truncated beach, approximating the conditions at the crest of a beach berm or inter-tidal ridge-runnel. The experiments provide a measure of the uprush sediment transport rate in the swash zone that is unaffected by the difficulties inherent in deploying instrumentation or sediment trapping techniques at laboratory scale. Overtopping flow volumes are compared with an analytical solution for swash flows as well as a simple numerical model, both of which are restricted to individual swash events. The analytical solution underestimates the overtopping volume by an order of magnitude while the model provides good overall agreement with the data and the reason for this difference is discussed. Modelled flow velocities are input to simple sediment transport formulae appropriate to the swash zone in order to predict the overwash sediment transport rates. Calculations performed with traditional expressions for the wave friction factor tend to underestimate the measured transport. Additional sediment transport calculations using standard total load equations are used to derive an optimum constant wave friction factor of f_w = 0.024. This is in good agreement with a broad range of published field and laboratory data. However, the influence of long waves and irregular wave run-up on the overtopping and overwash remains to be assessed. The good agreement between modelled and measured sediment transport rates suggests that the model provides accurate predictions of the uprush sediment transport rates in the swash zone, which has application in predicting the growth and height of beach berms.     

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.     

On the evolution and run-up of breaking solitary waves on a mild sloping beach

This paper presents new laboratory experiments carried out in a supertank (300 m × 5 m × 5.2 m) of breaking solitary waves evolution on a 1:60 plane beach. The measured data are employed to re-examine existing formulae that include breaking criterion, amplitude evolution and run-up height. The properties of shoreline motion, underwater particle velocity and scale effect on run-up height are briefly discussed. Based on our analyses, it is evidently found that there exist five zones during a wave amplitude evolution course on the present mild slope. A simple formula which is capable of predicting maximum run-up height for a breaking solitary wave on a uniform beach with a wide range of beach slope (1:15–1:60) is also proposed. The calculated results from the present model agree favorably with available laboratory data, indicating that our method is compatible with other predictive models.     

Infra-gravity wave generation by the shoaling wave groups over beaches

A physical parameter, μb, which was used to meet the forcing of primary short waves to be off-resonant before wave breaking, has been considered as an applicable parameter in the infra-gravity wave generation. Since a series of modulating wave groups for different wave conditions are performed to proceed with the resonant mechanism of infra-gravity waves prior to wave breaking, the amplitude growth of incident bound long wave is assumed to be simply controlled by the normalized bed slope, βb. The results appear a large dependence of the growth rate, α, of incident bound long wave, separated by the three-array method, on the normalized bed slope, βb. High spatial resolution of wave records enables identification of the cross-correlation between squared short-wave envelopes and infra-gravity waves. The cross-shore structure of infra-gravity waves over beaches presents the mechanics of incident bound- and outgoing free long waves with the formation of free standing long waves in the nearshore region. The wave run-up and amplification of infra-gravity waves in the swash zone appear that the additional long waves generated by the breaking process would modify the cross-shore structure of free standing long waves. Finally, this paper would further discuss the contribution of long wave breaking and bottom friction to the energy dissipation of infra-gravity waves based on different slope conditions.     

Spectral characteristics of random wave run-up

The time series of shoreline variations (run-up variations) due to random waves have been measured on uniform sloping beaches with slopes ranging from $\text{1}\text{5}$ to $\text{1}\text{30}$ and the energy spectra of the variations (run-up spectra) have been examined. The main characteristics of run-up spectra obtained from the experimental results are as follows: (1) a phenomenon of energy saturation is seen in a high frequency region; and (2) the spectral energy densities are independent of offshore incident wave energy. In the saturation region, the run-up spectra show f⁻⁴ dependence and tan ⁴θ dependence (f: frequency, tanθ: beach slope). Only in a low frequency region, the energy densities increase with increasing incident wave energy. In addition to the experimental study, it is shown by numerical simulations that if run-up variations are formed by parabolas induced by bores running up and down on the beach surface, the spectra of the variations show f⁻⁴ dependence, and the low frequency run-up energy densities increase with increasing running-up velocities of bores.     

The influence of seaward boundary conditions on swash zone hydrodynamics

The influence of the seaward boundary condition on the internal swash hydrodynamics is investigated. New numerical solutions of the characteristics form of the nonlinear shallow-water equations are presented and applied to describe the swash hydrodynamics forced by breaking wave run-up on a plane beach. The solutions depend on the specification of characteristic variables on the seaward boundary of the swash zone, equivalent to prescribing the flow depth or the flow velocity. It is shown that the analytical solution of Shen and Meyer [Shen, M.C., Meyer, R.E., 1963. Climb of a bore on a beach. Part 3. Runup. J. Fluid Mech. 16, 113–125] is a special case of the many possible solutions that can describe the swash flow, but one that does not appear appropriate for practical application for real waves. The physical significance of the boundary conditions is shown by writing the volume and momentum fluxes in terms of the characteristic variables. Results are presented that illustrate the dependence of internal flow depth and velocity on the boundary condition. This implies that the internal swash hydrodynamics depend on the shape and wavelength of the incident bore, which differs from the hydrodynamic similarity inherent in the analytical solution. A solution appropriate for long bores is compared to laboratory data to illustrate the difference from the analytical solution. The results are important in terms of determining overwash flows, flow forces and sediment dynamics in the run-up zone.     

An experimental study on geometric characteristics of beach erosion profiles

The beach profile and sediment transport are very important factors in the design of coastal structures, and the beach profile is mainly affected by a number of parameters, such as wave height and period, beach slope, and the material properties of the bed. In this study, considering wave height (H₀=6.5, 11.5, 16, 20, 23, 26 and 30 cm), wave period (T=1.46 and 2.03 s), beach slope (m=1/10 and 1/15) and mean sediment diameter (d₅₀=0.18, 0.26, 0.33 and 0.40 mm), an experimental investigation of coastal erosion profile (storm profile) was carried out in a wave flume using regular waves, and geometric characteristics of erosion profile were determined by the resultant erosion profile. Dimensional and non-dimensional equations were obtained by using linear and non-linear regression methods through the experimental data and were compared with previously developed equations in the literature. The results have shown that the experimental data fitted well to the proposed equations with respect to the previously developed equations.     

Run-up of long solitary waves of different polarities on a plane beach

We study the run-up of long solitary waves of different polarities on a beach in the case of composite bottom topography: a plane sloping beach transforms into a region of constant depth. We confirm that nonlinear wave deformation of positive polarity (wave crest) resulting in an increase in the wave steepness leads to a significant increase in the run-up height. It is shown that nonlinear effects are most strongly pronounced for the run-up of a wave with negative polarity (wave trough). In the latter case, the run-up height of such waves increases with their steepness and can exceed the amplitude of the incident wave.     

Run-up of long waves on a beach: The influence of the incident wave form

The influence of the incident wave form on the extreme (maximal) characteristics of a wave at a beach (run-up and draw-down heights, run-up and draw-down velocities, and the breaking parameter) is studied. It is suggested to use in the calculations the definition of wavelength at a level of 2/3 of the maximal height, which to a certain degree correlates with the definition of the significant wavelength accepted in oceanology. Such a definition allows us to unify the relations for extreme run-up characteristics so that the influence of the incident wave form becomes insignificant. The obtained universal relations can be used for the estimates of run-up characteristics when the exact information about the form of the incident wave is not available.     

Numerical analysis of run-up oscillations under dissipative conditions

This paper presents laboratory and numerical simulations of run-up induced by irregular waves breaking on a gentle-sloping planar beach. The experimental data are well reproduced by a numerical model based on the nonlinear shallow water equations. By extending the incoming wave conditions considered in the laboratory experiments, the model is applied to study the run-up variability under highly energetic incoming conditions. The numerical results support the idea that, for cases characterized by the same incident peak frequency, infragravity run-up increases almost linearly with the offshore significant wave height. Moreover, the most energetic conditions lead to an upper limit of the swash similarity parameter of about 1.8.