Sediment transport and beach morphodynamics induced by free long waves,bound long waves and wave groups

New laboratory data are presented on the influence of free long waves, bound long waves and wave groups on sediment transport in the surf and swash zones. As a result of the very significant difficulties in isolating and identifying the morphodynamic influences of long waves and wave groups in field conditions, a laboratory study was designed specifically to enable measurements of sediment transport that resolve these influences. The evolution of model sand beaches, each with the same initial plane slope, was measured for a range of wave conditions, firstly using monochromatic short waves. Subsequently, the monochromatic conditions were perturbed with free long waves and then substituted with bichromatic wave groups with the same mean energy flux. The beach profile changes and net cross-shore transport rates were extracted and compared for the different wave conditions, with and without long waves and wave groups. The experiments include a range of wave conditions, e.g. high-energy, moderate-energy, low-energy waves, which induce both spilling and plunging breakers and different turbulent intensities, and the beaches evolve to form classical accretive, erosive, and intermediate beach states. The data clearly demonstrate that free long waves influence surf zone morphodynamics and promote increased onshore sediment transport during accretive conditions and decreased offshore transport under erosive conditions. In contrast, wave groups, which can generate both forced and free long waves, generally reduce onshore transport during accretive conditions and increase offshore transport under erosive conditions. The influence of the free long waves and wave groups is consistent with the concept of the relative fall velocity, H/w_sT, as a dominant parameter controlling net beach erosion or accretion. Free long waves tend to reduce H/w_sT, promoting accretion, while wave groups tend to increase the effective H/w_sT, promoting erosion.     

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.     

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.     

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).     

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.     

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.     

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.     

Morphodynamics of intertidal bar morphology on a macrotidal beach under low-energy wave conditions, North Lincolnshire, England

The morphological changes of multiple intertidal bars (ridges) on a macrotidal beach were examined under low-energy wave conditions during a spring-to-spring tidal cycle. The morphological response was coupled to the tidal water level variations and related residence times for swash processes and surf (breaking waves and bores) over the cross-shore profile. Spring tides induced a large spatial variation in water lines and small residence times for distinct processes. Neap tides narrowed the intertidal area and increased the time for certain processes to work on the sediment at one location. The observed morphological changes could be coupled to the stagnation of processes at a certain bar crest position. The action of surf (breaking waves and bores) played the major role in the onshore migration of the intertidal bars and the simultaneous erosion of the seaward flank. Swash action, responsible for the generation and migration of intertidal bars in microtidal settings, was not the dominant process in causing the observed morphological changes. Intertidal ridges on macrotidal beaches cannot be considered swash bars as suggested by most previous investigations into these morphological features.     

卷破波作用下沙质岸滩剖面形态判别式初探

通过波浪水槽实验,开展不同类型波浪作用下的沙质岸滩演化规律研究工作。本次实验研究不考虑比尺,采用1:10与1:20组成的复合沙质斜坡对岸滩进行概化,选取规则波和椭圆余弦波两种典型波浪作用,对波浪的传播、变形和破碎、上爬、回落过程以及波浪作用前后沙质岸滩床面地形进行了观测,探讨波浪作用下沙质岸滩剖面演化规律。本文实验工况中,规则波作用下,岸滩剖面呈现出沙坝剖面和滩肩剖面,椭圆余弦波作用下的岸滩剖面均呈滩肩形态,发现岸滩剖面形态不仅与波浪作用类型、强度、周期等因素相关,还与波浪破碎的强度等因素有关。通过对实验过程中现象的进行观察和分析,引入了卷破波水舌冲击角的概念。对波浪卷破破碎后形成的水流挟沙运动与岸滩剖面形态的关系进行定性分析,对水舌冲击角与Irribarren参数之间的关系进行定量分析,基于Irribarren参数与岸滩剖面形态的关系初步建立了波浪作用下沙质岸滩剖面形态判别关系式。通过本文实验结果和前人实验结果对趋势线进行拟合,求得其判别系数,判别式能够较好地划分淤积型岸滩、侵蚀型岸滩及过渡型岸滩三种岸滩形态。     

粤东后江湾冲流带滩面高频振动及动力作用分析

利用涌浪影响下短时段内的冲流带滩面高频高程数据和碎波带波流资料,在奇异谱分析(SSA)的基础上,以比研究了不同形态滩面的冲淤变化趋势、趋势分布形状、冲淤变化周期和冲淤变化强度,以及同一条剖面不同桩点间各因素间的变化关系;用交叉谱方法探索了每分钟滩面高频冲淤变化与碎波带长重力波间的作用关系.分析结果表明,滩角韵律地形引起的冲流分流作用促进了滩脊向滩谷的泥沙转运,冲流带滩面存在明显的长重力波频段的周期性冲淤振动,滩面冲淤振动强度由滩面下部向上部递减,碎波带长重力波对滩面高频冲淤变化起重要作用.     

粤东后江湾冲流带滩面高频振动及动力作用分析

利用涌浪影响下短时段内的冲流带滩面高频高程数据和碎波带波流资料,在奇异谱分析(SSA)的基础上,对比研究了不同形态滩面的冲淤变化趋势、趋势分布形状、冲淤变化周期和冲淤变化强度,以及同一条剖面不同桩点间各因素间的变化关系;用交叉谱方法探索了每分钟滩面高频冲淤变化与碎波带长重力波间的作用关系。分析结果表明,滩角韵律地形引起的冲流分流作用促进了滩脊向滩谷的泥沙转运,冲流带滩面存在明显的长重力波频段的周期性冲淤振动,滩面冲淤振动强度由滩面下部向上部递减,碎波带长重力波对滩面高频冲淤变化起重要作用。     

Modeling sediment transport in the swash zone: A review

A critical review of conceptual and mathematical models developed in recent decades on sediment transport in the swash zone is presented. Numerous studies of the hydrodynamics and sediment transport in the swash zone in recent years have pointed out the importance of swash processes in terms of science advancement and practical applications. Evidently, the hydrodynamics of the swash zone are complex and not fully understood. Key hydrodynamic processes include both high-frequency bores and low-frequency infragravity motions, and are affected by wave breaking and turbulence, shear stresses and bottom friction. The prediction of sediment transport that results from these complex and interacting processes is a challenging task. Besides, sediment transport in this oscillatory environment is affected by high-order processes such as the beach groundwater flow. Most relationships between sediment transport and flow characteristics are empirical, based on laboratory experiments and/or field measurements. Analytical solutions incorporating key factors such as sediment characteristics and concentration, waves and coastal aquifer interactions are unavailable. Therefore, numerical models for wave and sediment transport are widely used by coastal engineers. This review covers mechanisms of sediment transport, important forcing factors, governing equations of wave-induced flow, groundwater interactions, empirical and numerical relations of cross-shore and longshore sediment transport in the swash zone. Major advantages and shortcomings of various numerical models and approaches are highlighted and reviewed. These will provide coastal modelers an impetus for further detailed investigations of fluid and sediment transport in the swash zone.     

Measurements and modelling of the advection of suspended sediment in the swash zone by solitary waves

Novel laboratory experiments and numerical modelling have been performed to study the advection scales of suspended sediment in the swash zone. An experiment was designed specifically to measure only the sediment picked up seaward of the swash zone and during bore collapse. The advection scales and settling of this sediment were measured during the uprush along a rigid sediment-free beach face by a sediment trap located at varying cross-shore positions. Measurements were made using a number of repeated solitary broken waves or bores. Approximately 25% of the pre-suspended sediment picked up by the bores reaches the mid-swash zone (50% of the horizontal run-up distance), indicating the importance of the sediment advection in the lower swash zone. The pre-suspended sediment is sourced from a region seaward of the shoreline (still water line) which has a width of about 20% of the run-up distance. An Eulerian–Lagrangian numerical model is used to model the advection scales of the suspended sediment. The model resolves the hydrodynamics by solving the non-linear shallow water equations in an Eulerian framework and then solves the advection–diffusion equation for turbulence and suspended sediment in a Lagrangian framework. The model provides good estimates of the measured mass and distribution of sediment advected up the beach face. The results suggest that the correct modelling of turbulence generation prior to and during bore collapse and the advection of the turbulent kinetic energy into the lower swash is important in resolving the contribution of pre-suspended sediment to the net sediment transport in the swash zone.     

Process-based model for nearshore hydrodynamics, sediment transport and morphological evolution in the surf and swash zones

Hydrodynamics and sediment transport in the nearshore zone were modeled numerically taking into account turbulent unsteady flow. The flow field was computed using the Reynolds Averaged Navier–Stokes equations with a k–ε turbulence closure model, while the free surface was tracked using the Volume-Of-Fluid technique. This hydrodynamical model was supplemented with a cross-shore sediment transport formula to calculate profile changes and sediment transport in the surf and swash zones. Based on the numerical solutions, flow characteristics and the effects of breaking waves on sediment transport were studied. The main characteristic of breaking waves, i.e. the instantaneous sediment transport rate, was investigated numerically, as was the spatial distribution of time-averaged sediment transport rates for different grain sizes. The analysis included an evaluation of different values of the wave friction factor and an empirical constant characterizing the uprush and backwash. It was found that the uprush induces a larger instantaneous transport rate than the backwash, indicating that the uprush is more important for sediment transport than the backwash. The results of the present model are in reasonable agreement with other numerical and physical models of nearshore hydrodynamics. The model was found to predict well cross-shore sediment transport and thus it provides a tool for predicting beach morphology change.     

Physical modeling of the cross-shore short-term evolution of protected and unprotected beach nourishments

This paper illustrates the results of an experimental investigation (model-to-prototype length ratio equal to 12) carried out to reproduce the cross-shore evolution of nourished sandy beaches. New two-dimensional experiments were performed to study the short-term response of the cross-shore profile for both “soft” (unprotected) and “mixed” (protected by submerged breakwaters) beach fill projects. Due to the simplified reproduction of prototype conditions in a two-dimensional geometry, only cross-shore sediment transport is considered. The results are related to the immediate post-nourishment evolution and far from beach fill boundaries where long-shore gradients of long-shore sediment transport are likely to be negligible. Three different pseudo-random wave trains were generated in order to simulate both accretive and erosive conditions. A fourth wave train, characterised by time-varying incident wave spectrum was generated for the investigation of the beach response to simplified storm time evolution. Dimensionless experimental results are given in terms of wave parameters, key features of cross-shore profile evolution and sediment transport rates. Furthermore, being highly resolved in both time and space, experimental data are suitable for mathematical model validation. It was observed that submerged breakwater switches erosive conditions to slightly accretive, at least within the tested experimental range.     

On the cross-shore profile change of gravel beaches

This paper investigates cross-shore profile changes of gravel beaches, with particular regard to discussing the tendency for onshore transport and profile steepening in the swash zone. The discussion includes observed morphological changes on a gravel beach from experimental investigations at the Large Wave Flume (GWK) in Hanover, Germany. During the tests all the profile changes occurred in the swash zone, resulting in erosion below the still water line (SWL) and formation of a berm above the SWL. We investigate the profile evolution evaluating the transport rates from a bed load sediment transport formulation coupled with velocities calculated from a set of Boussinesq equations that have been validated for its use in the surf and swash zones [Lynett, P.J., Wu, T.-R., and Liu, L.-F., P., 2002. Modelling wave runup with depth-integrated equations. Coastal Engineering, 46, 89–107; Otta, A.K., and Pedrozo-Acuña, A., 2004. Swash boundary and cross-shore variation of horizontal velocity on a slope. In: J.M. Smith (Editor), Proceedings 29th International Conference on Coastal Engineering. World Scientific, Lisbon, Portugal, pp. 1616–1628]. We discuss the influence of bottom friction on the predicted profiles, using reported friction factors from experimental studies. It is shown that the use of a different friction factor within a realistic range in each phase of the swash (uprush and backwash) improves prediction of the beach profiles, although quantitative agreement between the measured and computed profile evolutions is not satisfactory. Furthermore, if the friction factor and the transport efficiency (C) of the sediment transport formulation are kept the same in the uprush and backwash, accurate representation of profile evolution is not possible. Indeed, the features of the predicted profiles are reversed. However, when the C parameter is set larger during the uprush than during the backwash, the predicted profiles are closer to the observations. Differences between the predicted profiles from setting non-identical C-values and friction factors for the swash phase, are believed to be linked to both the infiltration effects on the flow above the beachface and the more accelerated flow in the uprush.     

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.     

Swash zone sediment fluxes: Field observations

This paper describes newly obtained, high-frequency observations of beach face morphological change over numerous tidal cycles on a macrotidal sandy beach made using a large array of ultrasonic altimeters. These measurements enable the net cross-shore sediment fluxes associated with many thousands of individual swash events to be quantified. It is revealed that regardless of the direction of net morphological change on a tidal time scale, measured net fluxes per event are essentially normally distributed, with nearly equal numbers of onshore and offshore-directed events. The majority of swash events cause net cross-shore sediment fluxes smaller than ± 50 kg m^− 1 and the mean sediment flux per swash event is only O(± 1 kg m^− 1) leading to limited overall morphological change. However, much larger events which deposit or remove hundreds of kilograms of sand per meter width of beach occur at irregular intervals throughout the course of a tide. It was found that swash–swash interactions tend to increase the transport potential of a swash event and the majority of the swash events that cause these larger values of sediment flux include one or more interactions. The majority of the larger sediment fluxes were therefore measured in the lower swash zone, close to the surf/swash boundary where swash–swash interactions are most common. Despite the existence of individual swash events that can cause fluxes of sediment that are comparable to those observed on a tidal time scale, frequent reversals in transport direction act to limit net transport such that the beach face volume remains in a state of dynamic equilibrium and does not rapidly erode or accrete.     

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.