Spectral distribution of wave energy dissipation by salt marsh vegetation

Spectral energy dissipation of random waves due to salt marsh vegetation (Spartina alterniflora) was analyzed using field data collected during a tropical storm. Wave data (significant wave heights up to 0.4 m in 0.8 m depth) were measured over a two-day period along a 28 m transect using 3 pressure transducers. The storm produced largely bimodal spectra on the wetland, consisting of low-frequency swell (7–10 s) and high-frequency (2–4.5 s) wind-sea. The energy dissipation varied across the frequency scales with the largest magnitude observed near the spectral peaks, above which the dissipation gradually decreased. The wind-sea energy dissipated largely in the leading section of the instrument array in the wetland, but the low-frequency swell propagated to the subsequent section with limited energy loss. Across a spectrum, dissipation did not linearly follow incident energy, and the degree of non-linearity varied with the dominant wave frequency. A rigid-type vegetation model was used to estimate the frequency-dependent bulk drag coefficient. For a given spectrum, this drag coefficient increased gradually up to the peak frequency and remained generally at a stable value at the higher frequencies. This spectral variation was parameterized by employing a frequency-dependent velocity attenuation parameter inside the canopy. This parameter had much less variability among incident wave conditions, compared to the variability of the bulk drag coefficient, allowing its standardization into a single, frequency-dependent curve for velocity attenuation inside a canopy. It is demonstrated that the spectral drag coefficient predicts the frequency-dependent energy dissipation with more accuracy than the integral coefficient.     

Short-term statistics of waves measured off Ratnagiri,eastern Arabian Sea

The short-term wave characteristics are required for design and operation of industrial facilities within the coastal areas. Water surface displacement measured using waverider buoy moored at 13 m water depth in the eastern Arabian Sea off the west coast of India have been analyzed to study the short-term statistics of waves covering full one year period. The study indicates that the values of the observed maximum wave height as a function of duration are not consistent with the theoretical expected value. There is significant variation (1.29–2.19) in the ratio between highest 1% wave and significant wave height compared to the theoretical value of 1.67. The data recorded at 13 m water depth indicates that the significant wave height is ∼8% lower than that predicted by the conventional Rayleigh distribution. The theoretical bivariate log-normal distribution represents the joint distributions of wave heights and periods for the study area.     

Waves in the nearshore waters of northern Arabian Sea during the summer monsoon

Waves at 15 m water depth in the northern Arabian Sea are measured during the summer monsoon for a period of 45 days and the characteristics are described. The significant wave height varied from 1.1 to 4.5 m with an average value of 2.5 m. 75% of the wave height at the measurement location is due to the swells arriving from the south-west and the remaining is due to the seas from south-west to north-west. Wave age of the measured data indicates that the waves in the nearshore waters of northern Arabian Sea during the summer monsoon are swells with young sea.     

Wave measurements by pressure transducers using artificial neural networks

Underwater ultrasonic acoustic transducers are frequently used in ocean wave measurements as they measure surface level using acoustic waves. However, their effectiveness can be severely affected in rough sea conditions, when bubbles generated by breaking waves interfere with their acoustic signals. When the seas are rough, one therefore often has to rely on a pressure transducer, which is generally used as a back-up for the acoustic wave gauge. A pressure transfer function is then used to obtain the surface wave information. Alternatively, the present study employed an artificial neural network to convert the pressure signal into significant wave height, significant wave period, maximum wave height, and spectral peakedness parameter using data obtained from various water depths. The results showed that, for water depths greater than 20 m, the wave parameters obtained from the artificial neural network were significantly closer to those obtained by the acoustic measurements than those obtained by using a linear pressure transfer function. Moreover, for a given water depth, the wave heights estimated by the network model from pressure data were not as good as those estimated by linear wave theory for large wave heights (above a 4 m significant wave height in this study). This can be improved if the training data set has more records with large wave heights.     

Solitary wave generation dynamics at Luzon Strait

A high resolution modeling study is undertaken, with a 2.5-dimensional nonhydrostatic model, of the generation of internal waves induced by tidal motion over the ridges in Luzon Strait. The model is forced by the barotropic tidal components K1, M2, and O1. These tidal components, along with the initial density field, were extracted from data and models. As the barotropic tide moves over the Luzon Strait sills, there is a conversion of barotropic tidal energy into baroclinic tidal energy. Depressions are generated that propagate towards the Asian Seas International Acoustics Experiment (ASIAEX) test site on the Chinese continental shelf. Nonlinear effects steepen the depressions, frequency and amplitude dispersion set in, and disintegration into large amplitude solitary waves occurs. The effects of varying the initial density field, tidal component magnitudes, as well as adding a steady background current to represent the occasional excursions of the Kuroshio Current into the strait, are considered.Depressions are generated at each of the two sills in Luzon Strait which radiate away, steepening and evolving into internal solitary wave trains. Baroclinic fluxes of available potential energy, kinetic energy and linear are calculated for various parameter combinations. The solitary wave trains produced in the simulations generally consist of large amplitude wave trains alternating with small amplitude wave trains. During strong tidal flow, Kelvin–Helmholtz type instabilities can develop over the taller double-humped sill. The solitary waves propagating towards the ASIAEX test site have been observed to reach amplitudes of 120–250 m, depending on the tidal strength. ASIAEX observations indicate amplitudes up to 150 m and the Windy Island Experiment (WISE) measurements contain magnitudes over 200 m. The model results yield solitary wave amplitudes of 70–300 m and half widths of 0.60–3.25 km, depending on parameter values. These are in the range of observations. Measurements by Klymak et al. (2006), in the South China Sea, exhibit amplitudes of 170 m, half widths of 3 km and phase speeds of 2.9 m s^?1. Model predictions indicate that the solitary waves making up the wave packet each experience different background currents with strong near surface shear.The energy in the leading soliton of the large amplitude wave trains ranges between 1.8 and 9.0 GJ m^?1. The smaller value, produced using barotropic tidal currents based on the Oregon State University data base, is the same as the energy estimated to be in a solitary wave observed by Klymak et al. (2006). Estimates of the conversion of barotropic tidal energy into radiating internal wave energy yield conversion rates ranging between 3.6% and 8.3%.     

Wave attenuation in mangroves: A quantitative approach to field observations

Coastal mangroves, dwelling at the interface between land and sea, provide an important contribution to reducing risk from coastal hazards by attenuating incident waves and by trapping and stabilizing sediments. This paper focusses on relations between vegetation densities, wave attenuation rates, sediment characteristics and sedimentation rates in mangroves. These processes were studied along two cross-shore transects through mangroves fringing estuaries in the southern Andaman region of Thailand. Volumetric vegetation densities in these mangroves were ranging up to 32‰, depending on the water depth. Generalized total wave attenuation rates increased from 0.002 m^− 1 in the sparsely vegetated forest fringes with Avicennia and Sonneratia species, up to 0.012 m^− 1 in the dense Rhizophora vegetation in the back of the forests. The total wave attenuation rates integrate effects of shoaling and energy losses due to various bio-physical interactions within the mangrove ecosystem. Wave attenuation in the mangroves is presumably dominated by energy losses due to vegetation drag, since wave attenuation due to bottom friction and viscous dissipation on the bare mudflats is significantly lower than those inside the mangrove vegetation.Additionally, wave attenuation in the mangroves was found to facilitate enhanced net sediment deposition and a gradual fining of the bed material. These findings corroborate the coastal defence function of mangroves by quantifying their contribution to wave attenuation and sediment trapping. The explicit linking of these properties to vegetation composition and structure facilitates modelling studies investigating the mechanisms determining the coastal defence capacities of mangroves.     

Wave climate scatter performance of a cycloidal wave energy converter

A lift based cycloidal wave energy converter (WEC) was investigated using potential flow numerical simulations in combination with viscous loss estimates based on published hydrofoil data. This type of wave energy converter consists of a shaft with one or more hydrofoils attached eccentrically at a radius. The main shaft is aligned parallel to the wave crests and submerged at a fixed depth. The operation of the WEC as a wave-to-shaft energy converter interacting with straight crested waves was estimated for an actual ocean wave climate. The climate chosen was the climate recorded by a buoy off the north-east shore of Oahu/Hawaii, which was a typical moderate wave climate featuring an average annual wave power P_W = 17 kWh/m of wave crest. The impact of the design variables radius, chord, span and maximum generator power on the average annual shaft energy yield, capacity factor and power production time fraction were explored. In the selected wave climate, a radius R = 5 m, chord C = 5 m and span of S = 60 m along with a maximum generator power of P_G = 1.25 MW were found to be optimal in terms of annual shaft energy yield. At the design point, the CycWEC achieved a wave-to-shaft power efficiency of 70%. In the annual average, 40% of the incoming wave energy was converted to shaft energy, and a capacity factor of 42% was achieved. These numbers exceeded the typical performance of competing renewables like wind power, and demonstrated that the WEC was able to convert wave energy to shaft energy efficiently for a range of wave periods and wave heights as encountered in a typical wave climate.     

Geometry and migration characteristics of bedforms under waves and currents: Part 2: Ripples superimposed on sandwaves

This paper presents results of a series of detailed measurements of geometric and migrating characteristics of ripples superimposed upon sandwaves under the action of combined waves and currents. Velocity measurements within the fluid, surface wave characteristics and 3D mapping of the bottom were recorded with an Acoustic Doppler Velocimeter (ADV), an acoustic water level sensor and a 32 composite element array of sub-aquatic acoustic sensors, respectively. Bottom records were statistically analyzed to obtain height, length and migration rates of ripples. Experiments examined ripple heights and wavelengths for the mobility factors (as defined in Eq. (4)) and the Reynolds wave number within the ranges 10 < ψ < 88 and 16 × 10³ < Re_w < 5 × 10⁵, respectively. Measured values were compared with laboratory and field data together with semi-empirical and analytical formulae from the literature. Good correlation was obtained when plotting measured ripple length and length in dimensionless form as a function of the Reynolds wave number Re_w. Under a given hydraulic condition, it was observed that ripples with different geometric characteristics may coexist at different locations over the sandwave. Ripple steepness is presented as a function of the Shields parameter although characterized with rather large scatter. Finally, average ripple migration speed is presented as a function of the Shields parameter and the mobility number.     

Numerical simulations of wave–current flow in an ocean basin

Wave–current flow is a phenomenon that is present in many practical engineering situations. Over the past several decades, this type of flow has been increasingly investigated under controlled laboratory conditions. This paper presents a numerical study of wave–current flow in the ocean basin of the LabOceano (COPPE/UFRJ). A homogeneous multiphase model based on the RANS equations and the k–ɛ turbulence model implemented in ANSYS-CFX code were used. A cross section of the ocean basin was represented. A regular wave with a height of 0.08 m and a period of 1.80 s (i.e., a wave steepness of H/L = 0.016), propagating on favourable currents, was simulated. The behaviour of the free surface elevation over time and the streamlines along the basin for wave and wave–current flows were presented. The numerical results were compared to the non-viscous theory given by the Rayleigh equation applied to the problem of wave–current interaction. Good agreement was found between the wave length estimated by the numerical results and the analytical solutions, with a deviation of less than 2%.     

Kinematics of breaking tsunami wavefronts: A data set from large scale laboratory experiments

This paper provides an overview of a new large scale laboratory data set on the kinematics of breaking tsunami wavefronts. The aim of the experiments was to provide an open access data set for model testing, calibration and verification, with particular emphasis on fluid kinematics in the wave breaking and run-up (swash) zones. The experiments were performed over a composite slope in the tsunami wave basin at the O. H. Hinsdale Wave Research Laboratory at Oregon State University. Data for ten different wave conditions were collected, including non-breaking and breaking waves, and both shore breaks and fully developed long bores.Surface elevation and fluid kinematics were measured with a closely spaced array of surface piercing wave gauges, non-contact ultrasonic wave gauges and four 3-D side-looking Acoustic Doppler Velocimeters. The array was traversed from the nearshore (depth = 0.2 m) to the middle and upper run-up zone, providing kinematic data at 30 cross-shore locations. Video was also recorded from 4 cameras covering the propagation, breaking and run-up zones. Surface elevation, flow velocities and the wave maker displacement were also recorded to provide offshore boundary conditions.The experiments include conditions with wave heights up to 0.55 m, notional wave periods up to 20 s and run-up lengths of up to 15.2 m on a 1/30 slope. In terms of the slope in the shoaling and breaker zones, the data correspond to Iribarren numbers in the range of 0.26–5.6. Raw, calibrated and processed data are stored with open access within the OSU Tsunami Wave Basin Experiment Notebook, which provides full access to all the wave maker control signals, data, instrument coordinates, and processing and plotting software. This paper serves as an introduction to the data set, demonstrates data quality and provides an initial analysis of some key parameters that govern the impact of tsunami events, including run-up versus offshore wave conditions and nearshore bore height, the maximum inundation depths at the original shoreline position, and the time to maximum inundation depth and flow reversal. Examples of temporal and convective accelerations and turbulent flow components are also presented to illustrate further details of the kinematics.     

Numerical study of turbulence and wave damping induced by vegetation canopies

Vegetation canopies control mean and turbulent flow structure as well as surface wave processes in coastal regions. A non-hydrostatic RANS model based on NHWAVE (Ma et al., 2012) is developed to study turbulent mixing, surface wave attenuation and nearshore circulation induced by vegetation. A nonlinear k − ϵ model accounting for vegetation-induced turbulence production is implemented to study turbulent flow within the vegetation field. The model is calibrated and validated using experimental data from vegetated open channel flow, as well as nonbreaking and breaking random wave propagation in vegetation fields. It is found that the drag-related coefficients in the k − ϵ model C_fk and C_fϵ can greatly affect turbulent flow structure, but seldom change the wave attenuation rate. The bulk drag coefficient C_D is the major parameter controlling surface wave damping by vegetation canopies. Using the empirical formula of Mendez and Losada (2004), the present model provides accurate predictions of vegetation-induced wave energy dissipation. Wave propagation through a finite patch of vegetation in the surf zone is investigated as well. It is found that the presence of a finite patch of vegetation may generate strong pressure-driven nearshore currents, with an onshore mean flow in the unvegetated zone and an offshore return flow in the vegetated zone.     

How does multiscale modelling and inclusion of realistic palaeobathymetry affect numerical simulation of the Storegga Slide tsunami?

The ∼8.15 ka Storegga submarine slide was a large (∼3000 km³), tsunamigenic slide off the coast of Norway. The resulting tsunami had run-up heights of around 10–20 m on the Norwegian coast, over 12 m in Shetland, 3–6 m on the Scottish mainland coast and reached as far as Greenland. Accurate numerical simulations of Storegga require high spatial resolution near the coasts, particularly near tsunami run-up observations, and also in the slide region. However, as the computational domain must span the whole of the Norwegian-Greenland sea, employing uniformly high spatial resolution is computationally prohibitive. To overcome this problem, we present a multiscale numerical model of the Storegga slide-generated tsunami where spatial resolution varies from 500 m to 50 km across the entire Norwegian-Greenland sea domain to optimally resolve the slide region, important coastlines and bathymetric changes. We compare results from our multiscale model to previous results using constant-resolution models and show that accounting for changes in bathymetry since 8.15 ka, neglected in previous numerical studies of the Storegga slide-tsunami, improves the agreement between the model and inferred run-up heights in specific locations, especially in the Shetlands, where maximum run-up height increased from 8 m (modern bathymetry) to 13 m (palaeobathymetry). By tracking the Storegga tsunami as far south as the southern North sea, we also found that wave heights were high enough to inundate Doggerland, an island in the southern North Sea prior to sea level rise over the last 8 ka.     

Physical model studies on wave transmission of a submerged inclined plate breakwater

This paper examines the results of physical model studies conducted in a monochromatic wave flume, to evaluate the wave transmission characteristics of a submerged plate breakwater consisting of a fixed plate of 0.50 m length and 0.003 m thickness. The model was oriented at varying inclinations and submergence. The influence of wave steepness, relative depth, relative submergence and angle of inclination on wave transmission was analysed. It was found that the horizontal plate is effective for short waves with steepness parameter higher than 5×10^?3 in relative depth grater than 0.21. The plate oriented at an angle of inclination of 60° is found to be effective for the entire ranges of wave parameters considered for the study and it reduces the wave height by about 40%.     

Distributions of freak wave heights measured in the North Sea

We present a statistical analysis of some of the largest waves occurring during 793 h of surface elevation measurements collected during 14 severe storms in the North Sea. This data contains 104 freak waves. It is found that the probability of occurrence of freak waves is only weekly dependent on the significant wave height, significant wave steepness and spectral bandwidth. The probability does show a slightly stronger dependency on the skew and kurtosis of the surface elevation data, but on removing the contribution to these measures from the presence of the freakwaves themselves, this dependency largely disappears.Distributions of extreme waves are modelled by fitting Generalised Pareto distributions, and extreme value distributions and return periods are given for freak waves in terms of the empirical fitted parameters. It is shown by comparison with these fits that both the Rayleigh distribution and the fit of Nerzic and Prevosto severely under-predict the probability of occurrence of extreme waves. For the most extreme freak wave in our data, the Rayleigh distribution over-predicts the return period by about 300 times when compared to the fitted model.     

Observations of surfzone alongshore pressure gradients onshore of an ebb-tidal delta

The relative importance of radiation stress gradients and alongshore pressure gradients to surfzone dynamics is investigated using observations of water levels, waves, and flows measured onshore of a large ebb-tidal delta. Incident wave heights measured along the ~ 11-m depth contour varied about 10% over a 1.2-km alongshore transect, resulting in alongshore wave setup differences on the order of 10 cm over the 600-m extent of the surfzone instrument array in 1.5-m depth. Despite the moderate alongshore variability in wave heights, the southerly alongshore pressure gradient, associated with the alongshore variability of wave-driven set-up, was typically twice as large as the northerly radiation stress gradient forcing, consistent with the observed southerly currents during the week-long experiment. The magnitude of the alongshore forcing and resulting alongshore velocity is reproduced by the two-dimensional depth-averaged numerical model of Shi et al. (JGR-Oceans, 2011). These observations, together with the numerical results, indicate that moderate alongshore wave height gradients (O(10^− 4)) outside the surfzone owing to alongshore variations in the offshore bathymetry can result in alongshore pressure gradients that are larger than radiation stress gradients.     

Experimental study of interactions between multi-directional focused wave and vertical circular cylinder,part II: Wave force

We describe experiments with multi-directional focused waves interacted with a vertical circular cylinder in a 3D wave basin. The focus of this study is on the run-up of multi-directional focused waves, wave forces, and wave pressures on the cylinder. Part I, the study on wave run-up, has already been presented by Li et al. (2012). In this paper, the analysis of the wave force on the vertical cylinder is presented.In this experiment, a cylinder with 0.25 m in diameter was adopted and different wave parameters, such as focused wave amplitude, peak frequency, frequency bandwidth and directional spreading index, are considered. The model scale k_pa (k_p is the wave number corresponding to peak frequency, a is the radium of the cylinder) varies from 0.32 to 0.65. The maximum forces of multi-directional focused wave on cylinder were measured and investigated. The results showed that the wave parameters have a significant influence on the wave force, and that the spatial profile of the surface of multi-directional focused wave can also affect its force on the cylinder, which is different from two-dimensional wave. In addition, the ‘secondary loading cycle’ phenomenon was also observed and discussed. In our experiments, the ‘secondary loading cycles’ occur when kA > 0.36 for all cases. While in some referred small scale experiments, the secondary load cycles are observed even for kA = 0.2, when the waves are longer enough. To larger model scale, the pronounced secondary load cycle occurs with larger wave steepness waves.     

Wave damping over artificial Posidonia oceanica meadow: A large-scale experimental study

An experimental study, conducted in the large wave flume of CIEM in Barcelona, is presented to evaluate the effects of Posidonia oceanica meadows on the wave height damping and on the wave induced velocities. The experiments were performed for irregular waves from intermediate to shallow waters with the dispersion parameter h/λ ranging from 0.09 to 0.29. Various configurations of the artificial P. oceanica meadow were tested for two stem density patterns (360 and 180 stems/m²) and for plant's height ranging from 1/3 to 1/2 of the water depth.The results for wave height attenuation are in good agreement with the analytical expressions found in literature, based on the assumption that the energy loss over the vegetated field is due to the drag forces. Based on this hypothesis, an empirical relationship for the drag coefficient related to the Reynolds number, Re, is proposed. The Reynolds number, calculated using the artificial P. oceanica leaf width as the length scale and the maximum orbital velocity over the meadow edge as the characteristic velocity scale, ranges from 1000 to 3500 and the drag coefficient C_d ranges from 0.75 to 2.0.The calculated wave heights, using the analytical expression from literature and the proposed relationship for the estimation of C_d, are in satisfactory agreement with those measured. Wave orbital velocities are shown to be significantly attenuated inside the meadow and just above the flume bed as indicated by the calculation of an attenuation parameter. Near the meadow edge, energy transfer is found in spectral wave velocities from the longer to the shorter wave period components. From the analysis it is shown that the submerged vegetation attenuates mostly longer waves.     

Modeling wave-mud interaction on the central chenier-plain coast,western Louisiana Shelf,USA

The strong coupling between hydrodynamics and seafloors on shallow muddy shelves, and resulting bed reworking, have been extensively documented. On these shelves, spectral wave transformation is driven by a complex combination of forcing mechanisms that include nonlinear wave interactions and wave energy dissipation induced by fluid-mud at a range of frequencies. Wave-mud interaction is investigated herein by using a previously validated nonlinear spectral wave model and observations of waves and near-bed conditions on a mildly-sloping seafloor off the muddy central chenier-plain coast, western Louisiana Shelf, United States. Measurements were made along a cross-shelf transect spanning 1 km between 4 and 3 m water depths. The high-resolution observations of waves and near-bed conditions suggest presence of a fluid mud layer with thickness sometimes exceeding 10 cm under strong long wave action (1 meter wave height with 7 s peak period at 4 meter depth). Spectral wave transformation is modeled using the stochastic formulation of the nonlinear Mild Slope Equation, modified to account for wave-breaking and mud-induced dissipation. The model is used in an inverse manner in order to estimate the viscosity of the fluid mud layer, which is a key parameter controlling mud-induced wave dissipation but complicated to measure in the field during major wave events. Estimated kinematic viscosities vary between 10⁻⁴-10⁻³ m²/s. Combining these results of the wave model simulations with in-depth analysis of near-bed conditions and boundary layer modeling allows for a detailed investigation of the interaction of nonlinear wave propagation and mud characteristics. The results indicate that mud-induced dissipation is most efficient when the wave-induced resuspensions of concentrations  > 10 g/L settle due to relatively small bottom stresses to form a fluid mud layer that is not as thin and viscous as a consolidated seafloor in absence of wave action but also not as thick and soft as a near-bed high concentration layer that forms during strong wave action.