Ephemeral sand waves in the hurricane surf zone

Airborne bathymetric LIDAR observations along the Florida panhandle after Hurricane Dennis (2005) show the first unequivocal observations of surf-zone sand wave trains.

These are found in depths of 5m along the trough of the hurricane bar, where hindcasts show strong longshore currents only during severe storms. The waves extend over tens of kilometers of coast after Dennis but are absent from the same area in four other datasets. Observed wavelength to water depth ratios are comparable to river dunes and tidal sand waves but height to depth ratios are smaller, with the largest wave heights around 0.1 times the water depth. The sand wave generation mechanism is hypothesized to be from wind-and-wave-induced longshore currents, which were hindcast to be large during Dennis, with destruction from water wave orbital velocities.     

SWAN predictions of waves observed in shallow water onshore of complex bathymetry

SWAN model predictions, initialized with directional wave buoy observations in 550-m water depth offshore of a steep, submarine canyon, are compared with wave observations in 5.0-, 2.5-, and 1.0-m water depths. Although the model assumptions include small bottom slopes, the alongshore variations of the nearshore wave field caused by refraction over the steep canyon are predicted well over the 50 days of observations. For example, in 2.5-m water depth, the observed and predicted wave heights vary by up to a factor of 4 over about 1000 m alongshore, and wave directions vary by up to about 10°, sometimes changing from south to north of shore normal. Root-mean-square errors of the predicted wave heights, mean directions, periods, and radiation stresses (less than 0.13 m, 5°, 1 s, and 0.05 m³/s² respectively) are similar near and far from the canyon. Squared correlations between the observed and predicted wave heights usually are greater than 0.8 in all water depths. However, the correlations for mean directions and radiation stresses decrease with decreasing water depth as waves refract and become normally incident. Although mean wave properties observed in shallow water are predicted accurately, nonlinear energy transfers from near-resonant triads are not modeled well, and the observed and predicted wave energy spectra can differ significantly at frequencies greater than the spectral peak, especially for narrow-band swell.     

Effect of roughness formulation on the performance of a coupled wave,hydrodynamic, and sediment transport model

A variety of algorithms are available for parameterizing the hydrodynamic bottom roughness associated with grain size, saltation, bedforms, and wave–current interaction in coastal ocean models. These parameterizations give rise to spatially and temporally variable bottom-drag coefficients that ostensibly provide better representations of physical processes than uniform and constant coefficients. However, few studies have been performed to determine whether improved representation of these variable bottom roughness components translates into measurable improvements in model skill. We test the hypothesis that improved representation of variable bottom roughness improves performance with respect to near-bed circulation, bottom stresses, or turbulence dissipation. The inner shelf south of Martha’s Vineyard, Massachusetts, is the site of sorted grain-size features which exhibit sharp alongshore variations in grain size and ripple geometry over gentle bathymetric relief; this area provides a suitable testing ground for roughness parameterizations. We first establish the skill of a nested regional model for currents, waves, stresses, and turbulent quantities using a uniform and constant roughness; we then gauge model skill with various parameterization of roughness, which account for the influence of the wave-boundary layer, grain size, saltation, and rippled bedforms. We find that commonly used representations of ripple-induced roughness, when combined with a wave–current interaction routine, do not significantly improve skill for circulation, and significantly decrease skill with respect to stresses and turbulence dissipation. Ripple orientation with respect to dominant currents and ripple shape may be responsible for complicating a straightforward estimate of the roughness contribution from ripples. In addition, sediment-induced stratification may be responsible for lower stresses than predicted by the wave–current interaction model.     

Modeling drift currents in a shallow basin with allowance for variations in tangential stresses induced by wind waves

A coupling model for calculating wind-driven currents and waves in a shallow basin with allowance for current-wave interactions is introduced. The model is constructed on the basis of the three-dimensional σ-coordinate model of currents [3] and the SWAN (Simulating Waves Nearshore) spectral wave model [4]. The effect of waves on currents is taken into account in the coefficients of surface and bottom friction through roughness parameters. Results of combined modeling of stationary fields of currents and waves generated by spatially homogeneous wind are correlated with the corresponding results of separate modeling for a cylindrical basin of constant depth and the water area of Lake Donuzlav (the northwestern coast of the Crimea). The allowance for the effect of waves during calculation of tangential wind stresses in the model of currents is shown to be among major factors intensifying water circulation and forming spatial inhomogeneities of the vortex type. In addition, some cases of local decreases in tangential wind stresses are revealed; they appear when the lake is penetrated from the side of the open sea by relatively long waves, which significantly decrease the roughness of the water surface.     

Coupled modeling of currents and wind waves in the Kerch Strait

We present a numerical model of the dynamics of the Kerch Strait allowing one to perform the coordinated analysis of the fields of currents and wind waves. The model includes the spectral wave module and the hydrodynamic block of currents. The influence of waves on the currents is taken into account in the hydrodynamic block both via the surface and bottom tangential stresses and via the radiation stresses. In order to take into account the inverse influence of currents upon the waves, we use the fields of currents and sea level from the hydrodynamic block in the wave module. The specific features of the structure of currents and wind waves in the strait are studied for the typical wave situations. The results of the coupled and separate simulation are compared and the importance of taking into account the mechanisms of interaction between waves and currents in the analysis of the dynamic processes in the strait is demonstrated. __________ Translated from Morskoi Gidrofizicheskii Zhurnal, No. 5, pp. 3–20, September–October, 2007.     

Modeling wave�Ccurrent bottom boundary layers beneath shoaling and breaking waves

The boundary layer characteristics beneath waves transforming on a natural beach are affected by both waves and wave-induced currents, and their predictability is more difficult and challenging than for those observed over a seabed of uniform depth. In this research, a first-order boundary layer model is developed to investigate the characteristics of bottom boundary layers in a wave–current coexisting environment beneath shoaling and breaking waves. The main difference between the present modeling approach and previous methods is in the mathematical formulation for the mean horizontal pressure gradient term in the governing equations for the cross-shore wave-induced currents. This term is obtained from the wave-averaged momentum equation, and its magnitude depends on the balance between the wave excess momentum flux gradient and the hydrostatic pressure gradient due to spatial variations in the wave field of propagating waves and mean water level fluctuations. A turbulence closure scheme is used with a modified low Reynolds number k-ε model. The model was validated with two published experimental datasets for normally incident shoaling and breaking waves over a sloping seabed. For shoaling waves, model results agree well with data for the instantaneous velocity profiles, oscillatory wave amplitudes, and mean velocity profiles. For breaking waves, a good agreement is obtained between model and data for the vertical distribution of mean shear stress. In particular, the model reproduced the local onshore mean flow near the bottom beneath shoaling waves, and the vertically decreasing pattern of mean shear stress beneath breaking waves. These successful demonstrations for wave–current bottom boundary layers are attributed to a novel formulation of the mean pressure gradient incorporated in the present model. The proposed new formulation plays an important role in modeling the boundary layer characteristics beneath shoaling and breaking waves, and ensuring that the present model is applicable to nearshore sediment transport and morphology evolution.     

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…     

Quadtree grid numerical model of nearshore wave–current interaction

This paper describes an adaptive quadtree-based 2DH wave–current interaction model for evaluating nearly horizontal wave-induced currents in the surf-zone. The model accounts for wave breaking, shoaling, refraction, diffraction, wave–current interaction, set-up and set-down, mixing processes (turbulent diffusion), bottom frictional effects, and movement of the land–water interface at the shoreline. The wave period- and depth-averaged governing equations, which conserve mass, momentum, energy and wave action, are discretised explicitly by means of an Adams–Bashforth second-order finite difference technique on adaptive hierarchical staggered quadtree grids. Grid adaptation is achieved through seeding points distributed according to flow criteria (e.g. local current gradients). The model is verified for nearshore circulation at a sinusoidal beach and nearshore currents at a multi-cusped beach. Reasonable agreement is obtained with experimental data from da Silva Lima [da Silva Lima, S.S.L., 1981. Wave-induced Nearshore Currents. PhD Thesis, Department of Civil Engineering, University of Liverpool] and Borthwick et al. [Borthwick, A.G.L., Foote, Y.L.M., Ridehalgh, A., 1997. Nearshore measurements at a cusped beach in the UK Coastal Research Facility, Coastal Dynamics '97, Plymouth, 953–962]. The modelling approach presented herein should be useful in simulating nearshore processes in complicated natural coastal domains. Of particular value is the local grid enrichment capability, which permits refined modelling of important localised flow behaviour such as rip currents and surf-zone circulation systems.     

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.     

Simulation of Random Waves and Associated Laminar Bottom Shear Stresses

This work presents a new approach for simulating the random waves in viscous fluids and the associated bottom shear stresses. By generating the incident random waves in a numerical wave flume and solving the unsteady two-dimensional Navier-Stokes equations and the fully nonlinear free surface boundaiy conditions for the fluid flows in the flume, the viscous flows and laminar bottom shear stresses induced by random waves axe determined. The deterministic spectral amplitude method implemented by use of the fast Fourier transform algorithm was adopted to generate the incident random waves. The accuracy of the numerical scheme is confirmed by comparing the predicted wave spectrum with the target spectrum and by comparing the nanlerical transfer function between the shear stress and the surface elevation with the theoretical transfer function. The maximum bottom shear stress caused by random waves, computed by this wave model, is compared with that obtained by Myrhaug＇ s model （1995）. The transfer function method is also employed to determine the maximum shear stress, and is proved accurate.     

Numerical study on water waves and wave-induced longshore currents in Obaky coastal water

In this paper, the water waves and wave-induced longshore currents in Obaky coastal water which is located at the Mediterranean coast of Turkey were numerically studied. The numerical model is based on the parabolic mild-slope equation for coastal water waves and the nonlinear shallow water equation for the wave-induced currents. The wave transformation under the effects of shoaling, refraction, diffraction and breaking is considered, and the wave provides radiation stresses for driving currents in the model. The numerical results for the water wave-induced longshore currents were validated by the measured data to demonstrate the efficiency of the numerical model. Then the water waves and longshore currents induced by the waves from main directions were numerically simulated and analyzed based on the numerical results. The numerical results show that the movement of the longshore currents was different while the wave propagated to a coastal zone from different directions.     

On the vertical distribution of

The aim of this paper is to present an analytical expression for the vertical distribution of the correlation between the horizontal ( ) and vertical ( ) wave velocity components. This quantity, , which appears explicitly in the time-averaged momentum balance equations, has been shown to play an important role in the vertical distribution of wave-induced currents.The proposed formulation for is based on an identity that relates the effective (wave) shear stress to the effective (wave) normal stresses ( ² and ²) and to the vorticity of the oscillatory flow gw. This general expression has been applied to simplified situations and has been shown to degenerate into other existing formulations with comparable simplifying assumptions, viz. irrotational waves in shallow water over an arbitrary bottom topography and breaking waves over a horizontal bottom.The model has also been applied to the case of waves interacting with a depth-varying current over a horizontal bottom, in which preliminary results have been obtained for a simplified situation invoking linear (small-amplitude) wave theory.     

Seabed shear stresses under irregular waves plus current from Monte Carlo simulations of parameterized models

Shear stresses on a rough seabed under irregular waves plus current are calculated. Parameterized models valid for regular waves plus current have been used in Monte Carlo simulations, assuming the wave amplitudes to be Rayleigh-distributed. Numerical estimates of the probability distribution functions are presented. For waves only, the shear stress maxima follow a Weibull distribution, while for waves plus current, both the maximum and time-averaged shear stresses are well represented by a three-parameter Weibull distribution. The behaviour of the maximum shear stresses under a wide range of wave-current conditions has been investigated, and it appears that under certain conditions, the current has a significant influence on the maximum shear stresses. Results of comparison between predictions and measurements of the maximum bottom shear stresses from laboratory and field experiments are presented.     

Transport of Sediments Entrapped by Topographic Waves

In the Boussinesq approximation, for topographic waves entrapped by a sloping bottom, we determine mean currents induced by a wave due to nonlinearity with regard for turbulent viscosity and diffusion. We determine the thickness of the bottom boundary layer, the vertical turbulent exchange coefficients, and turbulent stresses on the upper boundary of the boundary layer depending on the parameters of the wave. In the diffusion approximation, we find the vertical distribution of the concentration of sediments suspended by the wave and the flow rates of sediments along and perpendicular to the isobaths. __________ Translated from Morskoi Gidrofizicheskii Zhurnal, No. 5, pp. 13–24, September–October, 2005.     

Momentum balance in shoaling gravity waves: Comment on ’shoaling surface gravity waves cause a force and a torque on the bottom’ by K. E. Kenyon

In a recent paper, Kenyon (2004) proposed that the wave-induced energy flux is generally not conserved, and that shoaling waves cause a mean force and torque on the bottom. That force was equated to the divergence of the wave momentum flux estimated from the assumption that the wave-induced mass flux is conserved. This assumption and conclusions are contrary to a wide body of observations and theory. Most importantly, waves propagate in water, so that the momentum balance generally involves the mean water flow. Although the expression for the non-hydrostatic bottom force given by Kenyon is not supported by observations, a consistent review of existing theory shows that a smaller mean wave-induced force must be present in cases with bottom friction or wave reflection. That force exactly balances the change in wave momentum flux due to bottom friction and the exchange of wave momentum between incident and reflected wave components. The remainder of the wave momentum flux divergence, due to shoaling or wave breaking, is compensated by the mean flow, with a balance involving hydrostatic pressure forces that arise from a change in mean surface elevation that is very well verified by observations.     

Wave-induced Benthic Velocity Variations in Shallow Waters

Waves propagating from deep water into shallow coastal areas produce oscillatory currents near the sea bottom. The magnitude of these currents depend upon the period and amplitude of the incoming waves, and the dissipation mechanism such as wave breaking and bottom friction. Field experiments in a gently shoaling bay, i.e. Cleveland Bay, Northern Australia, showed that there is a broad band of water at around 6 m depth, where the benthic surge velocities are maximum. Both further inshore and offshore, the bottom velocities were less than at 6 m depth, contrary to the normal expectation that the velocities should increase as the water becomes shallower. A new and computationally efficient wave model was developed and was able to reproduce experimental results for waves above 50 cm wave height, but not for small waves (wave height about 30 cm). One implication of this higher band of benthic surge velocities may be to produce high water turbidities in this region. Turbidity data from Cleveland Bay is consistent with this hypothesis.     

Downward transfer of momentum by wind-driven waves

A model for the downward transfer of wind momentum is derived for growing waves. It is shown that waves, which grow due to an uneven pressure distribution on the water surface or a wave-coherent surface shear stress have horizontal velocities out of phase with the surface elevation. Further, if the waves grow in the x-direction, while the motion is perhaps time-periodic at any fixed point, the Reynolds stresses associated with the organized motion are positive. This is in agreement with several field and laboratory measurements which were previously unexplained, and the new theory successfully links measured wave growth rates and measured sub-surface Reynolds stresses. Wave coherent air pressure (and/or surface shear stress) is shown to change the speed of wave propagation as well as inducing growth or decay. From air pressure variations that are in phase with the surface elevation, the influence on the waves is simply a phase speed increase. For pressure variations out of phase with surface elevation, both growth (or decay) and phase speed changes occur. The theory is initially developed for long waves, after which the velocity potential and dispersion relation for linear waves in arbitrary depth are given. The model enables a sounder model for the transfer to storm surges or currents of momentum from breaking waves in that it does not rely entirely on ad-hoc turbulent diffusion. Future models of atmosphere-ocean exchanges should also acknowledge that momentum is transferred partly by the organized wave motion, while other species, like heat and gasses, may rely totally on turbulent diffusion. The fact that growing wind waves do in fact not generally obey the dispersion relation for free waves may need to be considered in future wind wave development models.     

The vertical structure of combined wave–current flow

A simple numerical model, based on the Reynolds stress equations and k–ε turbulence closure scheme, is developed for the coastal wave and current bottom boundary layer. The current friction velocity is introduced to account for the effect of currents on waves. The implicit Crank–Nicolson finite difference method discretizes the governing equations. Vertical changing step grids with the constant ratio for two adjacent spatial steps are used together with the equal time steps in the modeling. Vertical profiles of mean current velocity and wave velocity amplitude are obtained. These modeled results are compared with the laboratory experimental data of Van Doorn [1981. Experimental investigation of near bottom velocities in water waves with and without a current. Report M1423, Delft Hydraulics Laboratory, Delft, The Netherlands; 1982. Experimenteel onderzoek naar het snelheidsveld in de turbulente bodemgrenslaag in een oscillerende stroming in een golftunnel. Report M1562, Delft Hydraulics Laboratory, Delft, The Netherlands]. It has been shown that modeled and observed (Van Doorn, T., 1981. Experimental investigation of near bottom velocities in water waves with and without a current. Report M1423, Delft Hydraulics Laboratory, Delft, The Netherlands; 1982. Experimenteel onderzoek naar het snelheidsveld in de turbulente bodemgrenslaag in een oscillerende stroming in een golftunnel. Report M1562, Delft Hydraulics Laboratory, Delft, The Netherlands) mean velocity profiles within the wave and current bottom boundary layer are in better agreement than outside. Modeled and observed (Van Doorn, T., 1981. Experimental investigation of near bottom velocities in water waves with and without a current. Report M1423, Delft Hydraulics Laboratory, Delft, The Netherlands) wave velocity amplitude profiles within the wave and current bottom boundary layer are in better agreement than outside. Modeled wave velocity amplitudes are in good agreement with the laboratory experimental data of Van Doorn [1982. Experimenteel onderzoek naar het snelheidsveld in de turbulente bodemgrenslaag in een oscillerende stroming in een golftunnel. Report M1562, Delft Hydraulics Laboratory, Delft, The Netherlands].     

Wave-tide-surge coupled simulation for typhoon Maemi

The main task of this study focuses on studying the effect of wave-current interaction on currents, storm surge and wind wave as well as effects of current induced wave refraction and current on waves by using numerical models which consider the bottom boundary layer and sea surface roughness parameter for shallow and smooth bed area around Korean Peninsula. The coupled system (unstructured-mesh SWAN wave and ADCIRC) run on the same unstructured mesh. This identical and homogeneous mesh allows the physics of wave-circulation interactions to be correctly resolved in both models. The unstructured mesh can be applied to a large domain allowing all energy from deep to shallow waters to be seamlessly followed. There is no nesting or overlapping of structured wave meshes, and no interpolation is required. In response to typhoon Maemi (2003), all model components were validated independently, and shown to provide a faithful representation of the system’s response to this storm. The waves and storm surge were allowed to develop on the continental shelf and interact with the complex nearshore environment. The resulting modeling system can be used extensively for prediction of the typhoon surge. The result show that it is important to incorporate the wave-current interaction effect into coastal area in the wave-tide-surge coupled model. At the same time, it should consider effects of depth-induced wave breaking, wind field, currents and sea surface elevation in prediction of waves. Specially, we found that: (1) wave radiation stress enhanced the current and surge elevation otherwise wave enhanced nonlinear bottom boundary layer decreased that, (2) wind wave was significantly controlled by sea surface roughness thus we cautiously took the experimental expression. The resulting modeling system can be used for hindcasting (prediction) the wave-tide-surge coupled environments at complex coastline, shallow water and fine sediment area like areas around Korean Peninsula.     

Numerical investigation of the influence of various factors on the propagation and transformation of solitary waves in the sea

For various stratifications and different types of bottom patterns we study the transformations of solitary perturbations of density appearing in the depth of the sea. In the two-dimensional case, under the assumption that the average dynamic characteristics weakly vary in time as compared with the wave characteristics, we deduce the equations for mean currents and waves taking into account vertical and horizontal viscosity and the diffusion of density. Numerical examples show that the stratification, bottom topography, nonlinearity, mean currents, and dissipation strongly affect both the process of splitting of a solitary wave into wave trains and their amplitude and length. The wave currents exhibit the oscillatory (train-like) character. It is emphasized that, in the case of propagation of solitary perturbations of density with dissipation, it is also important to take into account the combined influence of nonlinearity, currents, bottom topography, and stratification. Translated by Peter V. Malyshev and Dmitry V. Malyshev