skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Breaking wave field statistics with a multi-layer model
The statistics of breaking wave fields are characterised within a novel multi-layer framework, which generalises the single-layer Saint-Venant system into a multi-layer and non-hydrostatic formulation of the Navier–Stokes equations. We simulate an ensemble of phase-resolved surface wave fields in physical space, where strong nonlinearities, including directional wave breaking and the subsequent highly rotational flow motion, are modelled, without surface overturning. We extract the kinematics of wave breaking by identifying breaking fronts and their speed, for freely evolving wave fields initialised with typical wind wave spectra. The $$\varLambda (c)$$ distribution, defined as the length of breaking fronts (per unit area) moving with speed $$c$$ to $$c+{\rm d}c$$ following Phillips ( J. Fluid Mech. , vol. 156, 1985, pp. 505–531), is reported for a broad range of conditions. We recover the $$\varLambda (c) \propto c^{-6}$$ scaling without wind forcing for sufficiently steep wave fields. A scaling of $$\varLambda (c)$$ based solely on the root-mean-square slope and peak wave phase speed is shown to describe the modelled breaking distributions well. The modelled breaking distributions are in good agreement with field measurements and the proposed scaling can be applied successfully to the observational data sets. The present work paves the way for simulations of the turbulent upper ocean directly coupled to a realistic breaking wave dynamics, including Langmuir turbulence, and other sub-mesoscale processes.  more » « less
Award ID(s):
2122042
PAR ID:
10441138
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
968
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We investigate wind wave growth by direct numerical simulations solving for the two-phase Navier–Stokes equations. We consider the ratio of the wave speed $$c$$ to the wind friction velocity $$u_*$$ from $$c/u_*= 2$$ to 8, i.e. in the slow to intermediate wave regime; and initial wave steepness $ak$ from 0.1 to 0.3; the two being varied independently. The turbulent wind and the travelling, nearly monochromatic waves are fully coupled without any subgrid-scale models. The wall friction Reynolds number is 720. The novel fully coupled approach captures the simultaneous evolution of the wave amplitude and shape, together with the underwater boundary layer (drift current), up to wave breaking. The wave energy growth computed from the time-dependent surface elevation is in quantitative agreement with that computed from the surface pressure distribution, which confirms the leading role of the pressure forcing for finite amplitude gravity waves. The phase shift and the amplitude of the principal mode of surface pressure distribution are systematically reported, to provide direct evidence for possible wind wave growth theories. Intermittent and localised airflow separation is observed for steep waves with small wave age, but its effect on setting the phase-averaged pressure distribution is not drastically different from that of non-separated sheltering. We find that the wave form drag force is not a strong function of wave age but closely related to wave steepness. In addition, the history of wind wave coupling can affect the wave form drag, due to the wave crest shape and other complex coupling effects. The normalised wave growth rate we obtain agrees with previous studies. We make an effort to clarify various commonly adopted underlying assumptions, and to reconcile the scattering of the data between different previous theoretical, numerical and experimental results, as we revisit this longstanding problem with new numerical evidence. 
    more » « less
  2. Atmospheric fronts embedded in extratropical cyclones are high‐impact weather phenomena, contributing significantly to mid‐latitude winter precipitation. The three vital characteristics of the atmospheric fronts, high wind speeds, abrupt change in wind direction, and rapid translation, force the induced surface waves to be misaligned with winds exclusively behind the cold fronts. The effects of the misaligned waves under atmospheric cold fronts on air‐sea fluxes remain undocumented. Using the multi‐year in situ near‐surface observations and direct covariance flux measurements from the Pioneer Array off the coast of New England, we find that the majority of the passing cold fronts generate misaligned waves behind the cold front. Once generated, the waves remain misaligned, on average, for about 8 hr. The parameterized effect of misaligned waves in a fully coupled model significantly increases the roughness length (185%), drag coefficient (19%), and air‐sea momentum flux (11%). The increased surface drag reduces the wind speeds in the surface layer. The upward turbulent heat flux is weakly decreased by the misaligned waves because of the decrease in temperature and humidity scaling parameters being greater than the increase in friction velocity. The misaligned wave effect is not accurately represented in a commonly used wave‐based bulk flux algorithm. Yet, considering this effect in the current formulation improves the overall accuracy of parameterized momentum flux estimates. The results imply that better representing a directional wind‐wave coupling in the bulk formula of the numerical models may help improve the air‐sea interaction simulations under the passing atmospheric fronts in the mid‐latitudes. 
    more » « less
  3. Wind, wave, and acoustic observations are used to test a scaling for ambient sound levels in the ocean that is based on wind speed and the degree of surface wave development (at a given wind speed). The focus of this study is acoustic frequencies in the range 1-20 kHz, for which sound is generated by the bubbles injected during surface wave breaking. Traditionally, ambient sound spectra in this frequency range are scaled by wind speed alone. In this study, we investigate a secondary dependence on surface wave development. For any given wind-speed, ambient sound levels are separated into conditions in which waves are 1) actively developing or 2) fully developed. Wave development is quantified using the non-dimensional wave height, a metric commonly used to analyze fetch or duration limitations in wave growth. This simple metric is applicable in both coastal and open ocean environments. Use of the wave development metric to scale sound spectra is first motivated with observations from a brief case study near the island of Jan Mayen (Norwegian Sea), then robustly tested with long time-series observations of winds and waves at Ocean Station Papa (North Pacific Ocean). When waves are actively developing, ambient sound levels are elevated 2-3 dB across the 1-20 kHz frequency range. This result is discussed in the context of sound generation during wave breaking and sound attenuation by persistent bubble layers. 
    more » « less
  4. Abstract The turbulent ocean surface boundary layer (OSBL) shoals during daytime solar surface heating, developing a diurnal warm layer (DWL). The DWL significantly influences OSBL dynamics by trapping momentum and heat in a shallow near‐surface layer. Therefore, DWL depth is critical for understanding OSBL transport and ocean‐atmosphere coupling. A great challenge for determining DWL depth is considering wave‐driven Langmuir turbulence (LT), which increases vertical transport. This study investigates observations with moderate wind speeds (4–7 m/s at 10 m height) and swell waves for which breaking wave effects are less pronounced. By employing turbulence‐resolving large eddy simulation experiments that cover observed wind, wave, and heating conditions based on the wave‐averaged Craik‐Lebovich equation, we develop a DWL depth scaling unifying previous approaches. This scaling closely agrees with observed DWL depths from a year‐long mooring deployment in the subtropical North Atlantic, demonstrating the critical role of LT in determining DWL depth and OSBL dynamics. 
    more » « less
  5. Abstract Energy is transferred from the atmosphere to the ocean primarily through ocean surface waves, and the majority is dissipated locally in the near‐surface ocean. Observations of turbulent kinetic energy (TKE) in the upper ocean have shown dissipation rates exceeding law‐of‐the‐wall theory by an order of magnitude. The excess near‐surface ocean TKE dissipation rate is thought to be driven primarily by wave breaking, which limits wave growth and transfers energy from the surface wave field to the wave‐affected layer of the ocean. Here, the statistical properties of breaking wave dynamics in a coastal area are extracted from visible imagery and used to estimate TKE dissipation rates due to breaking waves. The statistical properties of whitecap dynamics are quantified with Λ(c), a distribution of total whitecap crest length per unit area as a function of crest speed, and used to compute energy dissipation by breaking waves, Sds. Sdsapproximately balances elevated subsurface dissipation in young seas but accounts for only a fraction of subsurface dissipation in older seas. The wind energy input is estimated from wave spectra from polarimetric imagery and laser altimetry. Sdsbalances the wind energy input except under high winds. Λ(c)‐derived estimates of TKE dissipation rates by breaking waves compare well with the atmospheric deficit in TKE dissipation, a measure of energy input to the wave field (Cifuentes‐Lorenzen et al., 2024). These results tie the observed atmospheric dissipation deficit and enhancement in subsurface TKE dissipation to wave driven energy transport, constraining the TKE dissipation budget near the air‐sea interface. 
    more » « less