Tidal-cycle changes in oscillation ripples on the inner part of an estuarine sand flat

Oscillation ripples form on subaqueous sand beds when wave-generated, near-bottom water motions are strong enough to move sand grains. The threshold of grain motion is the lower bound of the regime of oscillation ripples and the onset of sheet flow is the upper bound. Based on the relation between ripple spacing and orbital diameter, three types of symmetrical ripples occur within the ripple regime. In the lower part of the ripple regime (orbital ripples), spacing is proportional to orbital diameter; in the upper part (anorbital ripples) spacing is independent of orbital diameter. Between these regions occurs a transitional region (suborbital ripples).

Oscillation ripples develop on a sandy tidal flat in Willapa Bay, Washington, as a result of waves traversing the area when it is submerged. Because wave energy is usually low within the bay, the ripples are primarily orbital in type. This means that their spacing should respond in a systematic way to changes in wave conditions. During the high-water parts of some tidal cycles, ripples near the beach decrease in spacing during the latter stage of the ebb tide while ripples farther offshore do not change. Observations made over several tidal cycles show that the zone of active ripples shifts on- or offshore in response to different wave conditions.

Detailed bed profiles and current measurements taken during the high-water part of spring tides show the manner in which the oscillation ripples change with changes in orbital diameter. Changes in ripple spacing at the study site could be correlated with changes in orbital diameter in the manner suggested by the criterion for orbital ripples. However, there appeared to be a lag time between a decrease in orbital diameter and the corresponding decrease in ripple spacing. Absence of change during a tidal cycle could be attributed to orbital velocities below the threshold for grain motion that negated the effects of changes in orbital diameter.

Because changes in sand-flat ripples depend both upon changes in orbital diameter and upon the magnitude of the orbital velocity, exposed ripples were not necessarily produced during the preceding high tide. In fact, some ripples may have been just produced, while others, farther offshore, may have been produced an unknown number of tides earlier. Therefore, when interpreting past wave conditions over tidal flats from low-tide ripples, one must remember that wave periods have to be short enough to produce velocities greater than the threshold velocity for the orbital diameters calculated from the observed ripple spacings.     

Ripples in intertidal mud—a conceptual explanation

On protected mudflats and along sheltered tidal channel margins, wave- and current-generated ripples are frequently observed on surficial and subsurface mud beds, although such bedforms are generally not thought to occur in cohesive sediments. In this paper, examples of such ripple marks in the German Wadden Sea (back-barrier tidal flats of Spiekeroog island) and also along the west coast of Korea (Baeksu tidal flats) are documented and analyzed. The mud ripples are 5–8 cm in spacing and 0.3–0.8 cm in height, and are composed of slightly sandy to virtually pure mud (80–98% mud content). For the Spiekeroog study area, a comparison of in situ particle-size measurements of suspended matter and of dispersed mud collected from the ripples shows that the former consists of low-density flocs which are considerably larger than the constituent grains of the latter. To assess local wave effects, near-bed orbital velocities and orbital diameters were calculated on the basis of standard wave theory using estimated wave parameters at the time of the study (June 2004) as well as wave data recorded nearby within the back-barrier tidal basin. The relationships between grain size, morphometric ripple parameters, and the near-bed orbital diameter show the wave-generated mud ripples to be of the orbital post-vortex type. It is demonstrated that only short-period shoaling (intermediate water depth) waves with periods of 1.5–2.5 s and heights of 0.1–0.5 m are able to generate and maintain such ripples. Corresponding near-bed orbital velocities range from 8–32 cm s^–1 and near-bed orbital diameters from 6.25–10 cm. It can be anticipated that increased current shear and turbulence associated with higher and longer waves prevent ripple formation due to the resuspension of settled mud, and the breakdown of suspended flocs and aggregates into smaller particles which then tend to remain in suspension. The most plausible explanation for the formation of the mud ripples is that mud flocs and aggregates deposited from suspension around high-water slack tide under moderate weather conditions initially respond as single (non-cohesive) particles which are hydraulically equivalent to ambient very fine sands. During exposure at low tide, gradual loss of water transforms the rippled mud into increasingly more cohesive mud drapes which are more resistant to erosion. Unless destroyed during high-energy events, the mud ripples may remain intact long enough to become buried and thereby preserved. Indeed, occasional but persistent observations of ripples in sub-Recent to ancient mudrocks document their preservation potential.     

Modeling of the propagation of transient waves of moderate steepness

A theoretical approach is applied to predict the propagation and transformation of nonlinear water waves. A semi-analytical solution was derived by applying an eigenfunction expansion method. The solution is applied to analyze the effect of wave frequencies and wave steepness on the propagation of nonlinear waves. The main attention is paid to the wave profile, the wave energy spectrum, and the changes of wave profile and energy spectrum due to the interaction of wave components in a wave train. The results show that for waves of low steepness the nonlinear wave effects and effects associated with the interaction of water waves in a wave train are of secondary importance. For waves of moderate steepness and steep waves the effects associated with the interactions between waves in a wave train are becoming significant and a train of initially sinusoidal waves may drastically change its form within a short distance from its original position. The evolution of wave components has substantial effects on the wave spectrum. A train of initially very narrow-banded waves changes its one-peak spectrum to a multi-peak one in a fairly short period of time. Laboratory experiments were conducted in a wave flume to verify theoretical approaches. The free-surface elevation recorded by a system of wave gauges was compared with the results provided by the semi-analytical solution. Theoretical results are in a fairly good agreement with experimental data. A reasonable agreement between theoretical results and experimental data is observed often even for relatively steep waves.     

Generation and propagation of transient nonlinear waves in a wave flume

A semi-analytical nonlinear wavemaker model is derived to predict the generation and propagation of transient nonlinear waves in a wave flume. The solution is very efficient and is achieved by applying eigenfunction expansions and FFT. The model is applied to study the effect of the wavemaker and its motion on the generation and propagation of nonlinear waves. The results indicate that the linear wavemaker theory may be applied to predict only the generation of waves of low steepness for which the nonlinear terms in the kinematic wavemaker boundary condition and free-surface boundary conditions are of secondary importance. For waves of moderate steepness and steep waves these nonlinear terms have substantial effects on wave profile and wave spectrum just after the wavemaker. A wave spectrum corresponding to a sinusoidally moving wavemaker possesses a multi-peak form with substantial nonlinear components, which disturbs or may even exclude physical modeling in wave flumes. The analysis shows that the widely recognized weakly nonlinear wavemaker theory may only be applied to describe the generation and propagation of waves of low steepness. This is subject to further restrictions in shallow and deep waters because the kinematic wavemaker boundary condition as well as the nonlinear interaction of wave components and the evolution of wave energy spectrum is not properly described by weakly nonlinear wavemaker theory. Laboratory experiments were conducted in a wave flume to verify the nonlinear wavemaker model. The comparisons show a reasonable agreement between predicted and measured free-surface elevation and the corresponding amplitudes of Fourier series. A reasonable agreement between theoretical results and experimental data is observed even for fairly steep waves.     

Estimation of the wave-related ripple characteristics and induced bed shear stress

A large data set on ripples was collected and examined. A set of new formulas for the prediction of the ripple characteristics is proposed with an emphasis on the disappearance of the ripples. The ripple wavelength was observed to be proportional to the bottom wave excursion but also to be a function of the grain-related Shields parameter and wave period parameter introduced by Mogridge et al. (1994). The ripple steepness was found to be nearly constant for orbital ripples, and with a sharp decrease for suborbital ripples. Two empirical functions are added including the effects of the critical Shields parameters (inception of transport and inception of sheet flow), i.e. giving the boundaries for the ripple existence's domain. The proposed formulas yield better prediction capabilities compared to the previously published formulas, especially when ripples are washed out. The effect of the ripple characteristics on the roughness height and the calculation of the bed shear stress is also discussed. It appeared that the bed shear stress calculation is more sensitive to the empirical coefficient a_r introduced in the estimation of the ripple-induced roughness height or to the limits of existence of the ripples than the ripple characteristics themselves.     

Laboratory investigation on wave transmission through two rows of perforated hollow piles

Experimental investigations on perforated hollow piles in two rows were conducted in a two dimensional wave flume. The influence of water depth, incident wave steepness, clear spacing between the piles and the spacing of pile rows on transmission coefficient have been studied. The effect of staggering of piles in rows is investigated. The results are also compared with the results of experiments on piles without perforations. The investigations have revealed that perforated piles attenuate more wave energy than non-perforated piles. The transmission coefficient K_t decreases as the wave steepness increases for both non-perforated and perforated piles. For non-perforated piles as relative clear spacing between the piles (b/D) decreases, for waves of higher steepness, K_t decreases while for perforated piles as b/D decreases, K_t is decreasing for all the steepness considered. As the relative clear spacing between the pile rows (B/D) increases K_t initially decreases till B/D is around one and later it starts increasing for both non-perforated and perforated piles. Staggering of piles has little effect on K_t. It is also found that water depth has insignificant influence on transmission coefficient at higher steepness for both perforated and non-perforated piles. Wave period alone does not directly influence transmission coefficient K_t.     

Statistical characteristics of individual waves in laboratory wind waves

On the basis of data on the statistical characteristics of individual waves in laboratory wind waves reported in part I of this series, a self-consistent similarity regime is found to exist among properties of the individual waves, such as the nondimensional frequency, the wave number, the phase speed, and the steepness. Also, it is shown that forms of past empirical formulas for the development of the peak wave can be derived starting from the 3/2-power law, as an extension of the persent laboratory experimental data. In the derivation, only values of the coefficient of the 3/2-power law, and the fraction of momentum transferred from the wind retained by the wind waves, remain on an empirical basis.     

Boundary roughness and bedforms in the surf zone

Hydrodynamical models of the nearshore system frequently assume that a single friction coefficient is sufficient to represent flow conditions at a point in the surf zone. Furthermore, models attempting to relate bed configuration to surf zone flows have relied primarily upon the wave orbital velocity as an indicator of potential bedforms, and thus as the control on boundary roughness. The data presented here point out potential errors arising from either of these approaches. The results of a field experiment conducted at Wendake Beach, Ontario, show that at a single location in an active surf zone, the Darcy-Weisbach friction coefficient, f, varied by approximately 250% (in this case between 0.016 and 0.041).

It is also shown that existing bedform models, based upon primary wave motions alone, do not accurately predict conditions at this study site. For a relatively constant wave orbital velocity and velocity asymmetry, it is found that changes in bed roughness, as a result of bedform development, are reflected mainly in the vertical profile of the long-shore current velocity. A sequence of bedforms, from oscillatory ripples through flat bed, is inferred from the data, and found to be supported by diver observations and preserved primary sedimentary structures.     

Wave Breaking Phenomena of Irregular Waves Combined with Opposing Current

- Experimental study and theoretical analysis show that the critical value of relative wave height (H/ d)b given by Goda and the critical wave steepness (H/ L)b given by Michell and Miche can be adopted as the spilling breaking indices of regular waves. According to the same principle, a systematic theoretical analysis and experiment of irregular waves have been done by the authors in order to solve the breaking problem of irregular waves. It is indicated that the authors' method for determining wave breaking of regular waves can also be used for irregular waves.     

Parametric modelling of joint probability density distributions for steepness and asymmetry in deep water waves

Traditional wave steepness s = H/L does not define steep asymmetric waves uniquely. Three additional parameters characterising single zero-downcross waves in a time series are crest front steepness, vertical asymmetry factor and horizontal asymmetry factor. Parametric models for joint probability density distributions for deep water waves are presented. The joint distributions are for crest front steepness-wave height, vertical asymmetry factor-wave height, total wave steepness-wave height and wave height-wave period. The parametric models are estimated from zero-downcross analysis of wave data obtained from measurements at sea on the Norwegian continental shelf. The results of the analysis presented here can be used in the estimation of the probabilities of occurrence of steep asymmetric waves and breaking waves in deep water. Thus the results are useful for the practical naval architect and ocean engineer who are considering unusual events in the sea, the associated accidents or responses and the probability of occurrence of such events.     

粗糙度与风浪特征量关系的研究

在实验室风浪槽中观测风浪和风速,研究粗糙度与波面特征量的关系,发现风浪谱宽度增加,粗糙度增大。在窄谱时,粗糙度随谱宽的增加变化不明显,当波陡降低,粗糙度降低;在宽谱时,当谱宽度增加,即使波陡降低,粗糙度仍可增大。这一结果表明,波陡不足以完全决定粗糙度。当风浪波龄增加,粗糙度呈下降趋势,但由于谱宽度对粗糙度的影响,当波龄增加,部分波浪可有较大的粗糙度。由于这一因素,在粗糙度与波龄关系的观测结果中,数据点的散落不完全由观测误差造成。     

Numerical Investigation of Run-ups on Cylinder in Steep Regular Wave

The run-up on offshore structures induced by the steep regular wave is a highly nonlinear flow with a free surface.This article focuses on the investigation of the steep regular wave run-up on a single vertical cylinder by solving the Navier-Stokes equations. A numerical wave tank is established based on the open-source package to simulate the wave scattering induced by a vertical cylinder. The VOF method is applied to capture the large deformation and breaking of the free surface. The numerical model is validated by experimental results. The relative wave run-ups on the front face and the back face along the centerline of a cylinder are analyzed. The changes of the relative run-ups with the wave steepness, the relative diameter and the relative depth are studied. It is found that the relative run-ups on the front face and the back face of the cylinder depend mainly on the wave steepness and the relative diameter,while the dependence on the relative depth is weak. The empirical formulae are proposed to calculate the relative run-ups in terms of the wave steepness of incident regular waves and the relative diameter of a cylinder.     

Steepness and asymmetry of extreme waves and the highest waves in deep water

Traditional wave steepness s=H/L does not define steep asymmetric waves in a random sea uniquey. Three additional parameters characterising single zero-downcross waves in a time series are crest front steepness, vertical asymmetry factor and horizontal asymmetry factor. Results for steepness and asymmetry from zero-downcross analysis of wave data obtained from full scale measurements in deep water on the Norwegian continental shelf in 58 time series are presented. The analysis demonstrates clearly the asymmetry of both “extreme waves” and the highest waves. The period and height of the highest waves are also given together with their correlation to spectral parameters. The measured maximum wave heights are also compared with predicted values of maximum wave heights showing good agreement.     

The limiting wave height in wind-induced wave trains

The vertical acceleration threshold concept has been applied to evaluate the limiting wave height in the train of wind-induced waves propagated over a horizontal bottom. This concept yields very simple computation of the probability of breaking for stochastic sea in deep and finite water depths. The computations confirmed the available field and laboratory observations that the limiting wave steepness in the deep water is lower than the steepness predicted by Stokes. For shallow water depth, the limiting wave height is smaller than 0.55h. This conclusion is consistent with field as well as wave tank observations.     

The steepness and shape of wind waves

Variations are found in the shape and the steepness of wind-generated surface gravity waves between very young waves, such as seen in a laboratory tank, and larger waves of various wave ages encountered at sea as the result of wind stress over larger fetches. These differences in the characteristic shape of wind waves are presented as a function of the wave age. The wave steepness is also expressed as a function of wave age, the measurement of which is consistent with the 3/2-power law connecting wave height and characteristic period, normalized by the air friction velocity.     

Wave parameter gradients along the wave ray

Computer methods have been used to track selected deep-water waves along their orthogonals to the surf zone. Three bottom orbital parameters were computed at each step: diameter, maximum velocity and acceleration. These three parameters were then plotted along the several wave rays, thereby providing shoreward gradients.The three gradients have essentially the same geometry, so that any one can be taken as an indicator of the other two. Bottom slope changes are responsible for departure of the gradients from linearity. The bottom orbital gradients are quite different from gradients of wave height and wave length. These results from coasting waves cannot be applied directly to the case of forced, or wind-driven, waves.The non-linear nature of the bottom orbital gradients indicates that the near-shore bottom should exhibit a crude banding, more-or-less parallel with the coast, with different sediment transport and ripple-mark geometry from band to band.     

不同间距大直径圆筒防波堤波压力试验研究

通过开展物理模型试验,进行不规则波作用下不同间距大直径圆筒防波堤的波峰压力分布规律研究,分析了圆筒相对间距、相对水深、波陡等主要影响因素对圆筒波压力的影响规律,并通过波压力试验值与各家理论公式计算值的比较,进行不同间距圆筒波压力的计算方法研究。基于直立墙结构和大尺度孤立墩柱结构波压力公式给出修正系数以及横向衰减函数的拟合公式,用以计算圆筒周身不同位置处的波压力,从而为工程设计和理论研究提供参考。     

波群内单个波的波陡分布与波破碎

对波群内单个波的波陡分布和波破碎进行了实验研究。研究结果是,波群中波动的最大振幅出现在波群前部而不是出现在波群中央,这种不对称性导致波群前部单个波出现大波陡的概率大于后部单个波出现大波陡的概率;进一步的波破碎统计发现波群前部单个波破碎的频率是后部单个波破碎频率的4倍。因此认为,波群结构的不对称性能够导致单个波发生破碎的...     

The wave plus current flow over vortex ripples at an arbitrary angle

This work concerns the wave plus current flow over a sand bed covered by vortex ripples, with the current and the waves coming from different angles. Experiments were performed in a basin, where current and waves were perpendicular, in order to determine the conditions (current strength) leading to a regular ripple pattern formation. Numerical simulations were conducted changing the direction between the waves and the current from 0° to 90° and the ratio between the current strength and the wave orbital velocity from 0.2 to 1.5. Close to the bed, the current aligns parallel to the ripple crests, leading to a veering current profile with the vertical coordinate. The current-related friction coefficient was calculated. It was found that it decreases as the angle approaches 90°, while it increases for decreasing values of the current with a trend that can be described by a power law.