This work serves as an observation‐based exploration into the role of wave‐driven turbulence at the air‐sea interface by measuring Turbulent Kinetic Energy (TKE) dissipation rates above and below the sea surface. Subsurface ocean measurements confirm a TKE dissipation rate enhancement relative to the predicted law‐of‐the‐wall (
We present an investigation of the azimuthal bimodality of the wind-wave spectrum for waves shorter than the dominant scale comparing numerical model solutions of developing waves from idealized experiments using WAVEWATCH III (WW3). The wave solutions were forced with the “exact” Webb–Resio–Tracy (WRT) nonlinear energy fluxes and the direct interaction approximation (DIA) with three different combinations of wind input and breaking dissipation parameterizations. The WRT gives larger azimuthal bimodal amplitudes compared to the DIA regardless of wind input/dissipation. The widely used wind input/dissipation parameterizations (i.e., ST4 and ST6) generally give narrow directional distributions with relatively small bimodal amplitudes and lobe separations compared to field measurements. These biases are significantly improved by the breaking dissipation of Romero (R2019). Moreover, the ratio of the resolved cross- to downwind mean square slope is significantly lower for ST4 and ST6 compared to R2019. The overlap integral relevant for the prediction of microseisms is several orders of magnitude smaller for ST4 and ST6 compared to R2019, which nearly agrees with a semiempirical model.
Spectral gravity wave models generally cannot accurately predict the directional distribution which impacts their ability to predict the resolved down- and crosswind mean square slopes and the generation of microseisms. Our analysis shows that a directionally narrow spectral energy dissipation, accounting for long-wave–short-wave modulation, can significantly improve the directional distribution of the wind-wave spectrum by coupling to the nonlinear energy fluxes due to wave–wave interactions, which has important implications for improved predictions of the mean square slopes and the generation of microseisms.
- PAR ID:
- 10368877
- Publisher / Repository:
- American Meteorological Society
- Date Published:
- Journal Name:
- Journal of Physical Oceanography
- Volume:
- 52
- Issue:
- 7
- ISSN:
- 0022-3670
- Page Range / eLocation ID:
- p. 1549-1562
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract ε obs >ε p ), which appears to be fully supported by wave breaking highlighting the role of the transport terms in balancing the subsurface TKE budget. Simultaneous measurements of TKE dissipation rates on the atmospheric side capture a deficit relative to the law‐of‐the‐wall (ε obs <ε p ). This deficit is explained in terms of wave‐induced perturbations, with observed convergence to the law‐of‐the‐wall at 14 m above mean sea level. The deficit on the atmospheric side provides an estimate of the energy flux divergence in the wave boundary layer. An exponential function is used to integrate in the vertical and provide novel estimates of the amount of energy going into the wave field. These estimates correlate well with classic spectral input parameterizations and can be used to derive an effective wave‐scale, capturing wind‐wave coupling purely from atmospheric observations intimately tied to wave‐induced perturbations of the air‐flow. These atmospheric and oceanic observations corroborate the commonly assumed input‐dissipation balance for waves at wind speeds in the 8‐14 ms−1range in the presence of developed to young seas. At wind speeds above 14 ms−1under young seas ()observations suggest a deviation from the TKE input‐dissipation balance in the wave field. -
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
-
more » « less
Abstract Superpressure balloon data of unprecedented coverage from Loon LLC is used to investigate the seasonal and latitudinal variability of lower stratospheric gravity waves over the entire intrinsic frequency spectrum. We show that seasonal variability in both gravity wave amplitudes and spectral slopes exist for a wide range of intrinsic frequencies and provide estimates of spectral slopes in five latitudinal regions for all four seasons, in five different frequency windows. The spectral slopes can be used to infer gravity wave amplitudes of intrinsic frequencies as high as 70 cycles/day from gravity waves resolved in model and reanalysis data. We also show that a robust relationship between the phase of the quasi‐biennial oscillation and gravity wave amplitudes exists for intrinsic frequencies as high as the buoyancy frequency. These are the first estimates of seasonal and latitudinal variability of gravity wave spectral slopes and high‐frequency amplitudes and constitute a significant step toward obtaining observationally constrained gravity wave parameterizations in climate models.
-
more » « less
Abstract Drifting buoy observations of ocean surface waves in hurricanes are combined with modeled surface wind speeds. The observations include targeted aerial deployments into Hurricane Ian (2022) and opportunistic measurements from the Sofar Ocean Spotter global network in Hurricane Fiona (2022). Analysis focuses on the slope of the waves, as quantified by the spectral mean square slope. At low‐to‐moderate wind speeds (<15 m s−1), slopes increase linearly with wind speed. At higher winds (>15 m s−1), slopes continue to increase, but at a reduced rate. At extreme winds (>30 m s−1), slopes asymptote. The mean square slopes are directly related to the wave spectral shapes, which over the resolved frequency range (0.03–0.5 Hz) are characterized by an equilibrium tail (
) at moderate winds and a saturation tail ( ) at higher winds. The asymptotic behavior of wave slope as a function of wind speed could contribute to the reduction of surface drag at high wind speeds. -
more » « less
Abstract This is the first of a 2‐part series concerning remote observation and wave‐by‐wave analysis of the onset of breaking in the surf zone. In the surf zone, breaking waves drive nearshore circulation, suspend sediment, and promote air‐sea gas exchange. Nearshore wave model predictions often diverge from
in situ measurements near the break point location because common parameterizations do not account for the rapid changes that occur near the onset of breaking. This work presents extensive methodology to combine data from a line‐scanning LIDAR and thermal infrared cameras to detect breaking, classify breaker type, and measure geometric wave parameters on a wave‐by‐wave basis, which can be used to improve breaker parameterizations. Over 2,600 non‐breaking and 1,600 breaking waves are analyzed from data collected at the USACE Field Research Facility in Duck, NC, including 413 spilling and 111 plunging waves for which the onset of breaking was observed. Wave height is estimated using a spatio‐temporal method for wave tracking that preserves the sea surface elevation maximum and overcomes field of view limitations. Methods for estimating instantaneous wave speed are refined by fitting a skewed Gaussian function to each wave profile before tracking the peaks. Wave slope is estimated from a linear fit to the upper 80% of the wave face, which provides a robust metric and strong correlation with geometric wave slope defined relative to mean sea level. Finally, breaking wave face foam coverage is analyzed to assess common model assumptions about roller length for wave energy dissipation parameterizations.
