Morphodynamic evolution of a coastal lagoon entrance during swash overwash

New field observations of beach berm growth resulting in the sub-aerial closure of an intermittently closed and open lake or lagoon are presented. These perched estuarine systems frequently exhibit ephemeral entrances that respond to process forcing almost instantaneously on a geomorphic timescale, with closure by a supra-tidal berm of depositional origin. The observations were made following a mechanical opening of the entrance and show very rapid vertical growth of an initially low beach berm as a result of swash overtopping and sediment overwash. Very little simultaneous seaward progradation was observed. Sediment overwash volumes were obtained from survey data over a period of 12 days either side of spring tide and the observations also provide an accurate measure of swash sediment transport at the berm crest. A process-based parametric model is developed through robust parameterisations of wave run-up, wave run-up distributions and sediment transport and is applied to predict the total overwash transport into and infilling the estuary entrance. The model is tested against the field data and compares well with the observations. While the field data are site-specific, the proposed modelling framework represents a first step in modelling the complex growth and recovery of natural beach berms in broad-scale morphodynamic models.     

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.     

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.     

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.     

Impact of sea-level rise and coral mortality on the wave dynamics and wave forces on barrier reefs

A one-dimensional wave model was used to investigate the reef top wave dynamics across a large suite of idealized reef-lagoon profiles, representing barrier coral reef systems under different sea-level rise (SLR) scenarios. The modeling shows that the impacts of SLR vary spatially and are strongly influenced by the bathymetry of the reef and coral type. A complex response occurs for the wave orbital velocity and forces on corals, such that the changes in the wave dynamics vary reef by reef. Different wave loading regimes on massive and branching corals also leads to contrasting impacts from SLR. For many reef bathymetries, wave orbital velocities increase with SLR and cyclonic wave forces are reduced for certain coral species. These changes may be beneficial to coral health and colony resilience and imply that predicting SLR impacts on coral reefs requires careful consideration of the reef bathymetry and the mix of coral species.     

Solid-state C NMR analysis of size and density fractions of marine sediments: Insight into organic carbon sources and preservation mechanisms

Burial of organic carbon (OC) in ocean sediments acts as the ultimate long-term sink for both terrestrial and marine carbon, however, the mechanisms controlling the preservation of this carbon are poorly understood. To better understand these mechanisms, we applied solid-state ¹³C nuclear magnetic resonance (NMR) spectroscopy, along with elemental, stable carbon isotopic (δ¹³C) and lignin phenol analyses, to size and density fractions of sediments influenced by either mixed terrestrial and marine OC inputs (Washington Coast slope) or dominantly marine inputs (Mexican Margin). Elemental, isotopic and lignin analyses all reveal that within the Washington Coast sediment, the OC mixes linearly between nitrogen-poor and ¹³C-depleted, lignin-rich OC in the large and light fractions and nitrogen-rich and ¹³C-enriched, lignin-poor OC in the small and dense fractions, suggesting that this sediment contains a two-component mixture of terrestrial vascular plant- and marine-derived OC. The integral areas of each of seven NMR spectral regions in the different samples trend linearly when plotted versus δ¹³C signature, with most R² values of 0.78 or greater, demonstrating that the NMR spectra of the two sources of carbon also mix linearly between the two endmembers. The terrestrial endmember in this sediment appears to be dominated by lignin and black carbon whereas the source of the marine endmember is less clear from the NMR spectra. In contrast, all of the analyses indicate that OC in the Mexican Margin sediment fractions is homogenous and derives almost exclusively from marine sources. It appears that selective preservation of (bio)chemically recalcitrant lignin and black carbon is the primary mechanism of preservation of terrestrial OC, whereas mineral-protection is the dominant mechanism preserving marine OC in the Washington coast sediment. There is little evidence showing that either preservation mechanism functions in the Mexican Margin sediments.     

Cross-shore hydrodynamics within an unsaturated surf zone

This paper concerns the hydrodynamics induced by random waves incident on a steep beach. New experimental results are presented on surface elevation and kinematic probability density functions, cross-shore variation in wave heights, the fraction of broken waves and velocity moments. The surf zone is found to be unsaturated at incident wave frequencies, with a significant proportion of the incident wave energy remaining at the shoreline in the form of bores. Wave heights in both the outer and inner surf zones are best described by a full Rayleigh distribution [Thornton, E.B., Guza, R.T., 1983. Transformation of wave height distribution. J. Geophys. Res. 88, 5925–5938], rather than a truncated Rayleigh distribution as used by Battjes and Janssen (1978) [Battjes, J.A, Janssen, J.P., 1978. Energy loss and setup due to breaking of random waves. Proc. 16th Int. Conf. Coastal Eng. ASCE, New York, pp. 569–588]. A new parametric wave transformation model is outlined which provides explicit expressions for the fraction of broken waves and the energy dissipation rate within the surf zone. On steep beaches, the model appears to offer improved predictive capabilities over the original Battjes and Janssen model. Cross-shore variations in the velocity variance and velocity moments are best described using Linear Gaussian wave theory, with less than 20% of the velocity variance in the inner surf zone due to low frequency energy.     

Breakpoint generated surf beat induced by bichromatic wave groups

This paper presents new experimental data on 2-D surf beat generation by a time-varying breakpoint induced by bichromatic wave groups. The experimental investigation covers a broad range of wave amplitudes, short wave frequencies, group frequencies and modulation rates. The data include measurements of incident and outgoing wave amplitudes, breakpoint position, shoreline run-up and the cross-shore structure of both the short and long wave motion. Surf beat generation is shown to be in good agreement with theory [Symonds, G., Huntley, D.A., Bowen, A.J., 1982. Two dimensional surf beat: long wave generation by a time-varying breakpoint. J. Geophys. Res. 87, 492–498]. In particular, surf beat generation is dependent on the normalised surf zone width, which is a measure of the phase relationship between the seaward and shoreward breakpoint forced long waves, and linearly dependent on the short wave amplitude. The cross-shore structure of the long wave motion is also consistent with theory; at maximum and minimum surf beat generation, the mean breakpoint coincides with the nodal and anti-nodal points, respectively, for a free long wave standing at the shoreline. A numerical solution, using measured data as input, additionally shows that the phase relationship between the incident bound long wave and the outgoing breakpoint forced wave is consistent with the time-varying breakpoint mechanism.     

Direct bed shear stress measurements in bore-driven swash

Direct measurements of bed shear in the swash zone are presented. The data were obtained using a shear plate in medium and large-scale laboratory bore-driven swash and cover a wide range of bed roughness. Data were obtained across the full width of the swash zone and are contrasted with data from the inner surf zone. Estimates of the flow velocities through the full swash cycle were obtained through numerical modelling and calibrated against measured velocity data. The measured stresses and calculated flow velocities were subsequently used to back-calculate instantaneous local skin friction coefficients using the quadratic drag law. The data show rapid temporal variation of the bed shear stress through the leading edge of the uprush, which is typically two–four times greater than the backwash shear stresses at corresponding flow velocity. The measurements indicate strong temporal variation in the skin friction coefficient, particularly in the backwash. The general behaviour of the skin friction coefficient with Reynolds number is consistent with classical theory for certain stages of the swash cycle. A spatial variation in skin friction coefficient is also identified, which is greatest across the surf-swash boundary and likely related to variations in local turbulent intensities. Skin friction coefficients during the uprush are approximately twice those in the backwash at corresponding Reynolds number and cross-shore location. It is suggested that this is a result of the no-slip condition at the tip leading to a continually developing leading edge and boundary layer, into which high velocity fluid and momentum are constantly injected from the flow behind and above the tip region. Finally, the measured stress data are used to determine the asymmetry and cross-shore variation in potential sediment transport predicted by three forms of sediment transport formulae.     

The influence of seaward boundary conditions on swash zone hydrodynamics

The influence of the seaward boundary condition on the internal swash hydrodynamics is investigated. New numerical solutions of the characteristics form of the nonlinear shallow-water equations are presented and applied to describe the swash hydrodynamics forced by breaking wave run-up on a plane beach. The solutions depend on the specification of characteristic variables on the seaward boundary of the swash zone, equivalent to prescribing the flow depth or the flow velocity. It is shown that the analytical solution of Shen and Meyer [Shen, M.C., Meyer, R.E., 1963. Climb of a bore on a beach. Part 3. Runup. J. Fluid Mech. 16, 113–125] is a special case of the many possible solutions that can describe the swash flow, but one that does not appear appropriate for practical application for real waves. The physical significance of the boundary conditions is shown by writing the volume and momentum fluxes in terms of the characteristic variables. Results are presented that illustrate the dependence of internal flow depth and velocity on the boundary condition. This implies that the internal swash hydrodynamics depend on the shape and wavelength of the incident bore, which differs from the hydrodynamic similarity inherent in the analytical solution. A solution appropriate for long bores is compared to laboratory data to illustrate the difference from the analytical solution. The results are important in terms of determining overwash flows, flow forces and sediment dynamics in the run-up zone.