A fully nonlinear non-dispersive energy balance for surfzone waves is derived based on the nonlinear shallow water equations to study the nearshore dynamics of infragravity (IG) waves. Based on simulations of waves on a relatively moderate and mild beach slope with a non-hydrostatic wave-flow model (SWASH), the new theory shows that spatial gradients in IG energy flux are nearly completely balanced by the combined effect of bottom stresses and predominantly nonlinear triad interactions. The new balance confirms many features of existing weakly nonlinear theories, and yields an improved description in the inner surfzone where waves become highly nonlinear. A gain of IG energy flux throughout the shoaling and outer surfzones is driven by triad interactions between IG waves and pairs of sea-swell (SS) waves. The IG energy flux decreased in the inner surfzone, primarily through an energy cascade to the swell-band and superharmonic frequencies where wave energy is ultimately dissipated. Dissipation by bottom friction was weak on both slopes. The IG wave breaking, characterized by triads between three IG or two IG waves and one SS wave, was significant only deep inside the surfzone of the mild slope. Even though IG waves broke on the mild slope, nonlinear interactions between IG waves and pairs of SS waves were responsible for at least half of the net IG flux loss. @en

A 2-layer non-hydrostatic model with improved dispersive behaviour is presented. Due to the assumption of a constant non-hydrostatic pressure distribution in the lower layer, the dispersive behaviour is improved without much additional computational time. A comparison with linear wave theory showed that this 2-layer model gives a better result for the dispersion relation and shoaling of waves in intermediate water. This means that the 2-layer model is applicable in shallow and intermediate water depths (up to relative depths kh equals 4), whereas the 1-layer model is only applicable in shallow water depths (kh smaller than 1). Three laboratory experiments, including a fringing reef and a barred beach, were used to validate the presented mode for different hydrodynamic conditions. Based on these results, it can be concluded that the 2-layer model can be applied to accurately simulate the bulk wave height and spectral properties. The low frequency wave height, the setup and in particular the second order statistics contain more scatter, but the model accurately captured the general trend. Furthermore, the model showed good results for complex bathymetries in shallow to intermediate water.@en

Wave-induced currents are an ubiquitous feature in coastal waters that can spread material over the surf zone and the inner shelf. These currents are typically under resolved in non-hydrostatic wave-flow models due to computational constraints. Specifically, the low vertical resolutions adequate to describe the wave dynamics – and required to feasibly compute at the scales of a field site – are too coarse to account for the relevant details of the three-dimensional (3D) flow field. To describe the relevant dynamics of both wave and currents, while retaining a model framework that can be applied at field scales, we propose a two grid approach to solve the governing equations. With this approach, the vertical accelerations and non-hydrostatic pressures are resolved on a relatively coarse vertical grid (which is sufficient to accurately resolve the wave dynamics), whereas the horizontal velocities and turbulent stresses are resolved on a much finer subgrid (of which the resolution is dictated by the vertical scale of the mean flows). This approach ensures that the discrete pressure Poisson equation – the solution of which dominates the computational effort – is evaluated on the coarse grid scale, thereby greatly improving efficiency, while providing a fine vertical resolution to resolve the vertical variation of the mean flow. This work presents the general methodology, and discusses the numerical implementation in the SWASH wave-flow model. Model predictions are compared with observations of three flume experiments to demonstrate that the subgrid approach captures both the nearshore evolution of the waves, and the wave-induced flows like the undertow profile and longshore current. The accuracy of the subgrid predictions is comparable to fully resolved 3D simulations – but at much reduced computational costs. The findings of this work thereby demonstrate that the subgrid approach has the potential to make 3D non-hydrostatic simulations feasible at the scale of a realistic coastal region.@en

Wave forecasting in ocean and coastal waters commonly relies on spectral models based on the spectral action balance equation. These models assume that different wave components are statistically independent and as a consequence cannot resolve wave interference due to statistical correlation between crossing waves, as may be found in, for instance, a focal zone. This study proposes a statistical model for the evolution of wave fields over non-uniform currents and bathymetry that retains the information on the correlation between different wave components. To this end, the quasi-coherent model (Smit & Janssen, J. Phys. Oceanogr., vol. 43, 2013, pp. 1741-1758) is extended to allow for wave-current interactions. The outcome is a generalized action balance model that predicts the evolution of the wave statistics over variable media, while preserving the effect of wave interferences. Two classical examples of wave-current interaction are considered to demonstrate the statistical contribution of wave interferences: (1) swell field propagation over a jet-like current and (2) the interaction of swell waves with a vortex ring. In both examples cross-correlation terms lead to development of prominent interference structures, which significantly change the wave statistics. Comparison with results of the SWAN model demonstrates that retention of cross-correlation terms is essential for accurate prediction of wave statistics in shear-current-induced focal zones.@en