Cross-shore hydrodynamics within an unsaturated surf zone

This paper concerns the hydrodynamics induced by random waves incident on a steep beach. New experimental results are presented on surface elevation and kinematic probability density functions, cross-shore variation in wave heights, the fraction of broken waves and velocity moments. The surf zone is found to be unsaturated at incident wave frequencies, with a significant proportion of the incident wave energy remaining at the shoreline in the form of bores. Wave heights in both the outer and inner surf zones are best described by a full Rayleigh distribution [Thornton, E.B., Guza, R.T., 1983. Transformation of wave height distribution. J. Geophys. Res. 88, 5925–5938], rather than a truncated Rayleigh distribution as used by Battjes and Janssen (1978) [Battjes, J.A, Janssen, J.P., 1978. Energy loss and setup due to breaking of random waves. Proc. 16th Int. Conf. Coastal Eng. ASCE, New York, pp. 569–588]. A new parametric wave transformation model is outlined which provides explicit expressions for the fraction of broken waves and the energy dissipation rate within the surf zone. On steep beaches, the model appears to offer improved predictive capabilities over the original Battjes and Janssen model. Cross-shore variations in the velocity variance and velocity moments are best described using Linear Gaussian wave theory, with less than 20% of the velocity variance in the inner surf zone due to low frequency energy.     

Laboratory study of wave and turbulence characteristics in narrow-band irregular breaking waves

Wave elevations and water particle velocities were measured in a laboratory surf zone created by the breaking of a narrow-band irregular wave train on a 1/35 plane slope. The incident waves form wave groups that are strongly modulated. It is found that the waves that break close to the shoreline generally have larger wave-height-to-water-depth ratios before breaking than the waves that break farther offshore. After breaking, the wave-height-to-water-depth ratio for the individual waves approaches a constant value in the inner surf zone, while the standard deviation of the wave period increases as the still water depth decreases. In the outer surf zone, the distribution of the period-averaged turbulent kinetic energy is closely correlated to the initial wave heights, and has a wider variation for narrow-band waves than for broad-band waves. In the inner surf zone, the distribution of the period-averaged turbulent kinetic energy is similar for narrow-band waves and broad-band waves. It is found that the wave elevation and turbulent kinetic energy time histories for the individual waves in a wave group are qualitatively similar to those found in a spilling regular wave. The time-averaged transport of turbulent kinetic energy by the ensemble-averaged velocity and turbulence velocity under the irregular breaking waves are also consistent with the measurements obtained in regular breaking waves. The experimental results indicate that the shape of the incident wave spectrum has a significant effect on the temporal and spatial variability of wave breaking and the distribution of turbulent kinetic energy in the outer surf zone. In the inner surf zone, however, the distribution of turbulent kinetic energy is relatively insensitive to the shape of the incident wave spectrum, and the important parameters are the significant wave height and period of the incident waves, and the beach slope.     

Discussion of “Wave setup and setdown generated by obliquely incident waves” by T.-W. Hsu et al., Coastal Engineering, 53, 865–877, 2006

In the recent paper by Tai-Wen Hsu, John R.-C. Hsu, Wen-Kai Weng, Swun-Kwang Wang, and Shan-Hwei Ou (Coastal Engineering, 53, 865–877, 2006), the authors derived theoretical formulations for calculating the wave setup and setdown induced by obliquely incident waves on a beach. The derivation of an expression for setdown contains errors which would lead to an imbalance in longshore momentum flux outside the surfzone. We correct their derivation and give results in terms of the radiation stress concept in a general case including an oblique wave incidence. We also point out that the correct form of wave setdown is important to describe the zero-net force in the momentum balance outside the surfzone.     

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.     

Implementation and modification of a three-dimensional radiation stress formulation for surf zone and rip-current applications

Regional Ocean Modeling System (ROMS v 3.0), a three-dimensional numerical ocean model, was previously enhanced for shallow water applications by including wave-induced radiation stress forcing provided through coupling to wave propagation models (SWAN, REF/DIF). This enhancement made it suitable for surf zone applications as demonstrated using examples of obliquely incident waves on a planar beach and rip current formation in longshore bar trough morphology (Haas and Warner, 2009). In this contribution, we present an update to the coupled model which implements a wave roller model and also a modified method of the radiation stress term based on Mellor (2008, 2011a,b,in press) that includes a vertical distribution which better simulates non-conservative (i.e., wave breaking) processes and appears to be more appropriate for sigma coordinates in very shallow waters where wave breaking conditions dominate. The improvements of the modified model are shown through simulations of several cases that include: (a) obliquely incident spectral waves on a planar beach; (b) obliquely incident spectral waves on a natural barred beach (DUCK'94 experiment); (c) alongshore variable offshore wave forcing on a planar beach; (d) alongshore varying bathymetry with constant offshore wave forcing; and (e) nearshore barred morphology with rip-channels. Quantitative and qualitative comparisons to previous analytical, numerical, laboratory studies and field measurements show that the modified model replicates surf zone recirculation patterns (onshore drift at the surface and undertow at the bottom) more accurately than previous formulations based on radiation stress (Haas and Warner, 2009). The results of the model and test cases are further explored for identifying the forces operating in rip current development and the potential implication for sediment transport and rip channel development. Also, model analysis showed that rip current strength is higher when waves approach at angles of 5° to 10° in comparison to normally incident waves.     

Reply to “Discussion of wave setup and setdown generated by obliquely incident waves” by F. Shi and J.T. Kirby

The authors are grateful to the derivation and correction addressed by Shi and Kirby (2008) for the calculation of the wave setup and setdown induced by obliquely incident waves. The mathematical equations for the expression of setup and setdown are rederived in which the velocity potential is expanded to the second-order on the still water lever. The correct form of setdown presented by Shi and Kirby is confirmed in the present paper.     

Surf zone dynamics simulated by a Boussinesq type model. Part II: surf beat and swash oscillations for wave groups and irregular waves

This is the second of three papers on the modelling of various types of surf zone phenomena. In the first paper the general model was described and it was applied to study cross-shore motion of regular waves in the surf zone. In this paper, part II, we consider the cross-shore motion of wave groups and irregular waves with emphasis on shoaling, breaking and runup as well as the generation of surf beats. These phenomena are investigated numerically by using a time-domain Boussinesq type model, which resolves the primary wave motion as well as the long waves. As compared with the classical Boussinesq equations, the equations adopted here allow for improved linear dispersion characteristics and wave breaking is modelled by using a roller concept for spilling breakers. The swash zone is included by incorporating a moving shoreline boundary condition and radiation of short and long period waves from the offshore boundary is allowed by the use of absorbing sponge layers. Mutual interaction between short waves and long waves is inherent in the model. This allows, for example, for a general exchange of energy between triads rather than a simple one-way forcing of bound waves and for a substantial modification of bore celerities in the swash zone due to the presence of long waves. The model study is based mainly on incident bichromatic wave groups considering a range of mean frequencies, group frequencies, modulation rates, sea bed slopes and surf similarity parameters. Additionally, two cases of incident irregular waves are studied. The model results presented include transformation of surface elevations during shoaling, breaking and runup and the resulting shoreline oscillations. The low frequency motion induced by the primary-wave groups is determined at the shoreline and outside the surf zone by low-pass filtering and subsequent division into incident bound and free components and reflected free components. The model results are compared with laboratory experiments from the literature and the agreement is generally found to be very good. Finally the paper includes special details from the breaker model: time and space trajectories of surface rollers revealing the breakpoint oscillation and the speed of bores; envelopes of low-pass filtered radiation stress and surface elevation; sensitivity of surf beat to group frequency, modulation rate and bottom slope is investigated. Part III of this work (Sørensen et al., 1998) presents nearshore circulations induced by the breaking of unidirectional and multi-directional waves.     

Numerical experiments of swash oscillations on steep and gentle beaches

A weakly non-linear Boussinesq model with a slot-type shoreline boundary is used to simulate swash oscillations on beaches. Numerical simulations of swash were compared with laboratory measurements and in general good agreement found (less than 15% root-mean-square error of surface elevation except in regular waves). A series of numerical experiments on shoreline movement were then performed for a range of beach slopes and incident wave conditions. The resulting swash characteristics are then discussed in terms of their physical nature and spectral properties. On steep slopes, both individual bores and infragravity waves are equally significant in driving the swash while infragravity waves alone drive them on mild slopes. Swash excursions on any given slope are found to be highest when individual bores from a partially saturated surf zone ride on top of low-frequency waves. This is confirmed by the relationship found between swash excursion and wave groupiness in the surf zone. Swash excursions increase with increasing incident wave energy, even in fully saturated surf zones. However, a poor correlation is found between swash excursion and the surf similarity parameter due to the involvement of infragravity wave energy in the swash.     

Surf zone dynamics simulated by a Boussinesq type model. Part I. Model description and cross-shore motion of regular waves

This is the first of three papers on the modelling of various types of surf zone phenomena. In this first paper, part I, the model is presented and its basic features are studied for the case of regular waves. The model is based on two-dimensional equations of the Boussinesq type and it features improved linear dispersion characteristics, possibility of wave breaking, and a moving boundary at the shoreline. The moving shoreline is treated numerically by replacing the solid beach by a permeable beach characterized by an extremely small porosity. Run-up of nonbreaking waves is verified against the analytical solution for nonlinear shallow water waves. The inclusion of wave breaking is based on the surface roller concept for spilling breakers using a geometrical determination of the instantaneous roller thickness at each point and modelling the effect of wave breaking by an additional convective momentum term. This is a function of the local wave celerity, which is determined interactively. The model is applied to cross-shore motions of regular waves including various types of breaking on plane sloping beaches and over submerged bars. Model results comprise time series of surface elevations and the spatial variation of phase-averaged quantities such as the wave height, the crest and trough elevations, the mean water level, and the depth-averaged undertow. Comparisons with physical experiments are presented. The phaseaveraged balance of the individual terms in the momentum and energy equation is determined by time-integration and quantities such as the cross-sectional roller area, the radiation stress, the energy flux and the energy dissipation are studied and discussed with reference to conventional phase-averaged wave models. The companion papers present cross-shore motions of breaking irregular waves, swash oscillations and surf beats (part II) and nearshore circulations induced by breaking of unidirectional and multidirectional waves (part III).     

Quasi-3D modelling of nearshore currents

Existing concepts of wave-induced nearshore current models, in the cross-shore vertical plane (2DV) and depth-integrated (2DH), are combined to a quasi-3D mathematical model. This combination is tested for reproducing correct results in 2DV and 2DH situations. The importance of the various contributions to the wave-induced secondary circulation in the vertical plane is investigated for realistic parameter ranges, which leads to the conclusion that both the non-breaking and the breaking fraction of a random wave field in the surf zone generate important secondary currents.Additional computations show the relevance of a 3D-approach of nearshore currents, even in seemingly simple situations like a plane sloping beach with obliquely incident waves.     

Wave setup/setdown and longshore current on non-planar beaches

Wave-induced setup/setdown and longshore currents are examined theoretically for non-planar, concave-up beaches. A beach profile in which the still-water-depth is proportional to the horizontal distance offshore raised to the $\text{2}\text{3}$ power is examined in detail. The total mean-water-depth, which includes the sum of the still-water-depth for this $\text{2}\text{3}$ power beach profile plus the wave-induced setup/setdown, may be approximated shoreward of the breaker line in a best least-squares sense by a profile that is proportional to the horizontal distance offshore raised to the $\text{1}\text{2}$ power. The longshore current profile for this non-planar beach is found to differ significantly from that predicted analytically for a planar beach. The longshore current velocity does not vanish at the shoreline as in the planar beach case. In addition, the peak velocities and the total longshore transport of water are found to be less than for the corresponding planar beach cases. For a given concave-up beach profile, the influence of the lateral mixing increases as the wave height decreases.     

潮汐影响下海滩前滨波浪传播耗能过程分析

文章基于长乐海滩前滨剖面的实测波浪数据, 通过统计分析以及谱分析的方法, 探讨了潮汐过程中长乐海滩波浪参数及耗能过程的变化规律。结果表明, 观测期间内波浪以混合浪为主, 各测点谱型较宽, 存在多峰振荡现象。向岸传播过程中, 波能耗散的形式为窄频域向宽频域转变, 能量分布趋于分散, 高频波能减小, 低频波能反而有所上升, 波浪破碎后生成长重力波。破波带内的能量衰减与波浪传播距离具有良好的相关性, 破碎波能在破波带内大约衰减了98.3%。潮汐水位对波浪具有明显的调制作用。入射波能随潮汐水位的增加而有所增加, 且水位越高, 入射波能分布越分散。破波带内的有效波高和潮汐水位具有显著的正相关关系。潮汐过程中固定测点的波谱变化与波浪沿剖面的波谱变化具有明显的相似性。     

Magnitude and cross-shore distribution of bed return flow measured on natural beaches

Field measurements of cross-shore currents 0.25 m from the bed were made on two natural beaches under a range of incident wave conditions. The results indicated the presence of a relatively strong, offshore-directed mean current, both within and seaward of the surf zone. Typical velocities within the surf zone were of the order of 0.2–0.3 m/s. This bed return flow, or “undertow”, represents a mass conservation response, returning water seaward that was initially transported onshore in the upper water column, primarily above the trough of the incident waves. The measurements demonstrated that the bed return flow velocity increases with the incident wave height. In addition, the crossshore distribution of the bed return flow is characterised by a mid-surf zone maximum, which exhibits a strong decrease in velocity towards the shoreline and a more gradual decay in the offshore direction. Several bed return flow models based on mass continuity were formulated to predict the cross-shore distribution of the bed return flow under an irregular wave field and were compared with the field data. Best agreement was obtained using shallow water linear wave theory, after including the mass transport associated with unbroken waves. The contribution of the unbroken waves enables net offshore-directed bottom currents to persist outside the region of breaking waves, providing a mechanism, other than rip currents, to transport sediment offshore beyond the surf zone.     

The effect of groundwater on topographic changes in a gravel beach

In recent years, several attempts to stabilize the beach by control of the percolation of water have been proposed. However, morphodynamics in the surf zone is still not clear because of the complexity of wave actions and sediment transport. Especially, there is a little research on gravel beach morphodynamics including wave breaking in the surf zone. The present study investigates experimentally how groundwater level influences topographic changes in a gravel beach and simulates numerically the wave fields and flow patterns in the surf zone, considering the porosity of the media and the presence of groundwater. In experiments, water-level control tank was designed to control the simulated groundwater elevation and the wave flume was divided into two parts to maintain a constant mean water level. The experimental results show that the berm formed in the upper portion of the shoreline moves up the beach as the groundwater level falls and the lower the groundwater level, the steeper the beach surface. The numerical model was developed to clarify these features capable of simulating the difference of groundwater and mean water level. Numerical results showed different flow patterns due to the groundwater elevation; wave run-up weakens and wave run-down strengthens by the seaward currents caused by elevated groundwater. These deformations of the flow pattern explain well how the beach profile is affected by the groundwater elevation.     

Morphodynamics of a bar-trough surf zone

A field study was made of the distinguishing morphodynamic processes operating in a surf zone which perennially exhibits accentuated bar-trough topography (the “longshore-bar-trough” and “rhytmic-bar-and-beach” states as described by Wright and Short, 1984). Characteristic features of the morphology include a shallow bar with a steep shoreward face, a deep trough, and a steep beach face. This morphology, which is favored by moderate breaker heights and small tidal ranges, strongly controls the coupled suite of hydrodynamic processes. In contrast to fully dissipative surf zones, the bar-trough surf zone is not at all saturated and oscillations at incident wave frequency remain dominant from the break point to the subaerial beach. The degree of incident wave groupiness does not change appreciably across the surf zone. Infragravity standing waves which, in dissipative surf zones, dominate the inshore energy, remain energetically secondary and occur at higher frequencies in the bar trough surf zone. Analyses of the field data combined with numerical simulations of leaky mode and edge wave nodal—antinodal positions over observed surf-zone profiles, indicate that the frequencies which prevail are favored by the resonant condition of antinodes over the bar and nodes in the trough. Standing waves which would have nodes over the bar are suppressed. Sediment resuspension in the surf zone appears to be largely attributable to the incident waves which are the main source of bed shear stress. In addition, the extra near-bottom eddy viscosity provided by the reformed, non-breaking waves traversing the trough significantly affects the vertical velocity profile of the longshore current. Whereas the bar is highly mobile in terms of onshore—offshore migration rates, the beach face and inner regions of the trough are remarkably stable over time.     

Statistical characteristics of long waves nearshore

The random long wave runup on a beach of constant slope is studied in the framework of the rigorous solutions of the nonlinear shallow water theory. These solutions are used for calculation of the statistical characteristics of the vertical displacement of the moving shoreline and its horizontal velocity. It is shown that probability characteristics of the runup heights and extreme values of the shoreline velocity coincide in the linear and nonlinear theory. If the incident wave is represented by a narrow-band Gaussian process, the runup height is described by a Rayleigh distribution. The significant runup height can also be found within the linear theory of long wave shoaling and runup. Wave nonlinearity nearshore does not affect the Gaussian probability distribution of the velocity of the moving shoreline. However the vertical displacement of the moving shoreline becomes non-Gaussian due to the wave nonlinearity. Its statistical moments are calculated analytically. It is shown that the mean water level increases (setup), the skewness is always positive and kurtosis is positive for weak amplitude waves and negative for strongly nonlinear waves. The probability of the wave breaking is also calculated and conditions of validity of the analytical theory are discussed. The spectral and statistical characteristics of the moving shoreline are studied in detail. It is shown that the probability of coastal floods grows with an increase in the nonlinearity. Randomness of the wave field nearshore leads to an increase in the wave spectrum width.     

Morphodynamic zone variability on a microtidal barred beach

This paper builds on the work of Masselink [Masselink, G., 1993. Simulating the effects of tides on beach morphodynamics. J. Coast. Res. SI 15, 180–197.] on the use of the residence times of shoaling waves, breaking waves and swash/backwash motions across a cross-shore profile to qualitatively understand temporal beach behaviour. We use a data set of in-situ measurements of wave parameters (height and period) and water depth, and time-exposure video images overlooking our single-barred intertidal measurement array at Egmond aan Zee (Netherlands) to derive boundaries between the shoaling zone, the surf zone and the swash zone. We find that the boundaries are functional dependencies of the local relative wave height on the local wave steepness. This contrasts with the use of constant relative wave heights or water levels in earlier work. We use the obtained boundaries and a standard cross-shore wave transformation model coupled to an inner surf zone bore model to show that large (> 5) relative tide ranges (RTR, defined as the ratio tide range–wave height) indicate shoaling wave processes across almost the entire intertidal profile, with surf processes dominating on the beach face. When the RTR is between 2 and 5, surf processes dominate over the intertidal bar and the lower part of the beach face, while swash has the largest residence times on the upper beach face. Such conditions, associated with surf zone bores propagating across the bar around low tide, were observed to cause the intertidal bar to migrate onshore slowly and the upper beach face to steepen. For RTR values less than about 2, surf zone processes dominate across the intertidal bar, while the dominance of swash processes now extends across most of the beach face. The surf zone processes were now observed to lead to offshore bar migration, while the swash eroded the upper beach face.     

Large-wave simulation of three-dimensional,cross-shore and oblique,spilling breaking on constant slope beach

In this work, the large-wave simulation (LWS) method is adapted for application in spilling wave breaking over a constant slope beach. According to LWS, large scales of velocities, pressure and free-surface elevation are numerically resolved, while the corresponding unresolved scale effects are taken into consideration by a subgrid scale (SGS) model for wave and eddy stresses. The model may be not fully applicable in very shallow water, close to the shoreline, where the unresolved, turbulent, free-surface oscillation is of the same order with the water depth. Time integration of the Euler equations is achieved by a two-stage fractional scheme, combined with a hybrid scheme for spatial discretization, consisting of finite difference and pseudospectral approximation methods. Model parameters are calibrated by comparison to available experimental data of free-surface elevation and velocities in the surf zone for cross-shore incoming waves. The action of the wave SGS stresses in the outer coastal and surf zones initiates breaking and generates appropriate vorticity, in the form of an eddy structure (surface roller), at the breaking wavefront. At incipient breaking, both advection and gravity contribute to the vorticity flux at the free surface, while only after the full development of the surface roller, the effect of advection becomes stronger. The SGS model is also utilized to simulate propagation, refraction and breaking of oblique incoming waves. The gradual breaking and dissipation of wave crestlines and the surface roller structure along the breaking wavefront are automatically captured without any empirical input, such as data for the roller shape or the wave propagation angle at breaking.     

On Radiation Boundary Conditions and Wave Transformation Across the Surf Zone

The purpose of this paper is to extend the validity of Li's parabolic model (1994) by incorporating a combined energy factor in the mild-slope equation and by improving the traditional radiation boundary conditions. With wave breaking and energy dissipation expressed in a direct form in the equation, the proposed model could provide an efficient numerical scheme and accurate predictions of wave transformation across the surf zone. The radiation boundary conditions are iterated in the model without use of approximations. The numerical predictions for wave height distributions across the surf zone are compared with experimental data over typical beach profiles. In addition, tests of waves scattering around a circular pile show that the proposed model could also provide reasonable improvement on the radiation boundary conditions for large incident angles of waves.     

Experimental and Numerical Study of Wave-Induced Long-Shore Currents on A Mild Slope Beach

1.Introduction Long shorecurrenthasbecomethesubjectofextensiveworldwideresearchformanyyears.The purposeofthestudyistounderstandandpredicttheprocessesofsedimenttransport,shorelineevolu tionandpollutanttransportinthenear shorezoneundertheactionofwaves.Goda(2001)regarded thesuspendedsedimentastheprincipalloadinsedimenttransportandtheeffectoflong shorecurrents wasemphasized.Theeffectofwavesonpollutanttransportwasstudiedbynumericalmodelandfield experiment(TaoandHan,2002;Rodriguezetal.,1995),andth…