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

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

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.     

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.     

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.     

Wave groupiness variations in the nearshore

This paper proposes a new definition of the groupiness factor, GF, based on the envelope of the incident-wave time series. It is shown that an envelope-based GF has several important advantages over the SIWEH-based groupiness factor, including objective criteria for determining the accuracy of the envelope function and well-defined numerical limits.Using this new GF, the variability of incident wave groupiness in the field is examined both temporally, in unbroken waves at a fixed location, and spatially, in a cross-shore array through the surf zone. Contrary to previous studies using the SIWEH-based GF, results suggest that incident wave groupiness may not be an independent parameter in unbroken waves; through a wide range of spectral shapes, from swell to storm waves, the groupiness did not vary significantly. As expected, the groupiness decreases rapidly as waves break through the surf zone, although significant wave height variability persists even through a saturated surf zone. The source of this inner surf zone groupiness is not identified; however, this observation implies that models of long wave generation must account for nonsteady radiation stress gradients landward of some narrow zone near the mean breakpoint.     

浅滩地形上波浪破碎对低频波能放大的影响

低频长波对港湾共振、泥沙输运、波浪爬高与越浪等过程有重要影响。采用非静压模型SWASH模拟了不规则波在浅滩上的传播及破碎过程,重点探讨了滩顶短波破碎程度对低频波能演变的影响。模拟结果显示,当滩顶短波处于临界破碎状态时,低频波能流沿浅滩持续增长,浅滩对低频长波的放大效应（滩后与滩前低频波能流之比）显著;当滩顶短波轻微破碎、破波仅发生在外破波区时,低频波能流的沿滩增长率进一步变大;当滩顶短波剧烈破碎、破波延续到内破波区时,低频波能流沿滩先增长而后转为衰减,滩后放大率较临界破碎时明显减小。研究结果表明,浅滩顶部水深改变了短波破碎程度,进而影响低频长波的演化过程,浅滩上长波总体增长率随滩顶水深的减小呈现先增大后减小的规律,在短波轻微破碎时最大。     

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.     

Modelling of periodic wave transformation in the inner surf zone

In this paper, we analyse the ability of the nonlinear shallow-water (NSW) equations to predict wave distortion and energy dissipation of periodic broken waves in the inner surf zone. This analysis is based on the weak-solution theory for conservative equations. We derive a new one-way model, which applies to the transformation of non-reflective periodic broken waves on gently sloping beaches. This model can be useful to develop breaking-wave parameterizations (in particular broken-wave celerity expression) in both time-averaged wave models and time-dependent Boussinesq-type models. We also derive a new wave set-up equation which provides a simple and explicit relation between wave set-up and energy dissipation. Finally, we compare numerical simulations of both, the NSW model and the simplified one-way model, with spilling wave breaking experiments and we find a good agreement.     

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.     

Rip currents under obliquely incident wind waves and tidal longshore currents

A field experiment on the nature of rip currents was conducted on the Dutch coast, which differs from previous rip current study sites because it is a wind-sea dominated environment with mostly obliquely incident waves and tidally-driven longshore currents. During the experiment three distinct flow patterns, obtained with GPS tracked drifter instruments, were observed: (1) a locally governed circulation cell, (2) an offshore current that is deflected shore parallel outside the surf zone and (3) a meandering longshore current. The transition from rip currents (flow patterns 1 and 2) to meandering longshore currents (flow pattern 3) occurred gradually within the tidal cycle with longshore currents prevalent at mid to high tide. Rip currents at this site appeared at depressions in the surf zone bar and typically occurred when the water level fell below NAP (equivalent to MSL), even in the presence of obliquely incident waves and tidally driven longshore currents. Hindcast simulations of the drifter experiments were performed with the numerical model XBeach and showed good agreement with field observations. The model was subsequently used to investigate the influence of tidal water level fluctuations, longshore currents and obliquely incident waves on rip currents.Rip currents were initiated when the water level dropped below a specific threshold with the magnitude of the rip current associated with the water level. The strength of the tidal current and its orientation with respect to the incident waves governed the offshore extent and orientation of the rip current. In contrast to other studies that suggest that rip currents solely occur under shore normal (or slightly oblique waves), in this study both observations and numerical model simulations indicate that rip currents can exist under large angles of wave incidence, when the rip channel is sufficiently wide and the wave height is small.     

Non-linear propagation of unidirectional wave fields over varying topography

The present study aims at investigating the non-linear triad interaction process affecting shoaling surface gravity wave fields. The triad interaction phenomenon being enhanced towards the shore, the domain of study is extended up to the surf zone. Three 1D non-linear wave models (one phase-resolving and two phase-averaged spectral models) have been implemented and compared to laboratory experiments performed in a wave flume. This set of models includes two existing models and a new one which has been developed in the frame of this work. The models include a breaking dissipation term based on the parametrical model of Battjes and Janssen [Battjes, J.A., Janssen, P.A.E.M., 1978. Energy loss and set-up due to breaking of random waves. Proc. 16th Int. Conf. Coastal Eng. (ASCE), Vol. 1, pp. 569–587.]. The investigations concern the evolution of variance spectra, spectral significant wave height and mean period over a barred bathymetric profile. In addition, the performances of the different models are analysed by computing the spectral source term for triad interactions. We found that all models are able to reproduce the main features of non-linear mechanisms affecting a wave field in the near-shore zone. The phase-resolving model gives the most accurate results for non-breaking situations. It correctly reproduces the non-linear coupling effect in decreasing water depths due to wave–wave interactions, as well as the harmonic release after a bar. However, the model is computationally time-consuming. The CPU time is considerably reduced using phase-averaged models. They give satisfactorily results on harmonic generation. However, they do not reproduce the release of harmonics as water depth increases. In breaking conditions, the variance spectra undergo significant changes under the combined effects of non-linear energy transfers and dissipation. The depth-induced wave breaking model included in the equations provides a good estimate of the energy decay in the surf zone.     

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.     

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.     

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.     

珊瑚岸礁破碎带附近波浪演化实验研究

通过波浪水槽实验对珊瑚岸礁破碎带附近波浪演变规律开展研究,实验采用了概化的岸礁模型,测试了4种礁坪水深、4种礁前斜坡坡度和一系列入射波高的组合工况。对破碎带宽度和破碎带附近波浪的入射、反射、透射以及能量耗散进行了测量分析,透射波的计算考虑了礁坪上高次谐波的影响。结果表明:礁坪水深和入射深水波高的比值(即礁坪相对水深)是影响岸礁破碎带附近波浪演化的关键参数,而礁前斜坡坡度的影响在本文测量的范围内可以忽略不计。破碎带宽度与礁坪上浅水波波长为同一数量级,并与礁坪相对水深成反比;透射系数随礁坪相对水深的增大呈线性增长,而反射系数的变化却无类似规律;岸礁能够削弱超过50%入射波能,礁坪相对水深越小,波浪破碎造成的能量耗散越大。     

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.     

Turbulence and energy dissipations of surf-zone spilling breakers

A laboratory study on the turbulence and wave energy dissipations of spilling breakers in a surf zone is presented. Instantaneous velocity fields of propagating breaking waves on a 1/20 slope were measured using Particle Image Velocimetry (PIV). Due to the large region of the evolving wave breaking generated turbulent flow, seven PIV fields of view (FOVs) were mosaicked to form a continuous flow field in the surf zone. Mean and turbulence quantities were extracted by ensemble averaging 25 repeated instantaneous measurements at each FOV. New results for distribution and evolution of turbulent kinetic energy, mean flow energy, and total energy across the surf zone were obtained from analyzing the data. The turbulence dissipation rate was estimated based on several different approaches. It was found that the vertical distribution of the turbulence dissipation rate decays exponentially from the crest level to the bottom. The resulting energy budget and energy flux were also calculated. The calculated total energy dissipation rate was compared to that based on a bore approximation. It was found that the ratio of turbulence dissipation rate to total energy dissipation rate was about 0.01 in the outer surf zone and increased to about 0.1 after the breaking waves transformed into developed turbulent bores in the inner surf zone.     

珊瑚礁破碎带附近波浪演化和波生流实验研究

为了研究珊瑚岸礁破碎带附近波浪演化和波生流特性,通过水槽实验对规则波浪作用下珊瑚岸礁上沿礁分布的水位和流速进行了详细的测量。在典型卷破波条件下,测试了礁冠存在与不存在的两种情况。实验结果表明,多重波浪反射作用引起岸礁上形成不完全驻波,而破碎带附近的浅化作用则产生高次谐波,波浪破碎所耗散的波能主要来源于主频波,礁坪上透射波成分中二次谐波与主频波的能量相当;礁冠的存在引起破碎带宽度减小、礁坪上增水变大以及礁坪上各次谐波变小;礁冠不存在时,岸礁上波生流的沿礁分布与平直海岸相似,而礁冠的存在一定程度上阻碍了礁坪上水体向外海的回流。研究成果将丰富和发展珊瑚礁水动力学理论,并为岛礁工程的建设和维护提供一定的理论参考。     

Validation of the three-wave quasi-kinetic approximation for the spectral evolution in shallow water

This paper aims at validating the three-wave quasi-kinetic approximation for the spectral evolution of weakly nonlinear gravity waves in shallow water. The problem is investigated using a one-dimensional numerical wave propagation model, formulated in the spectral representation. This model includes both a nonlinear triad interactions term and a wave breaking dissipation term. Some numerical tests were carried out in order to show the importance of using the triad nonlinear term in wave propagation spectral models, particularly to describe both behavior of the spectral integral parameters and of the spectral shape evolution in shallow water depth. Furthermore; a comparison against different set of experimental observations was carried out. Comparing the numerical results with the experimental observations made it possible to show the modeling efficiency of the three-wave quasi-kinetic approximation.     

Wave breaking velocity effects in depth-integrated models

A simple model for predicting the velocities under breaking waves in depth-integrated models is developed. A velocity modification due to wave breaking is formulated based on a specific exponential profile, which is then added to the numerically predicted, depth-integrated velocity profile. This modification is superficial in that it does not directly change the hydrodynamic calculations inside the depth-integrated model. The modifications can be employed in any of the numerous Boussinesq-type models, and is not dependant on the use of a particular breaking dissipation scheme. Horizontal velocity profiles, both mean and instantaneous, are compared with experimental data in the surf zone. The comparisons show good agreement, markedly better than the un-modified results, and on par with published numerical results from sophisticated models.