The observation-based source terms available in the third-generation wave model WAVEWATCH III (i.e., the ST6 package for parameterizations of wind input, wave breaking, and swell dissipation terms) are recalibrated and verified against a series of academic and realistic simulations, including the fetch/duration-limited test, a Lake Michigan hindcast, and a 1-yr global hindcast. The updated ST6 not only performs well in predicting commonly used bulk wave parameters (e.g., significant wave height and wave period) but also yields a clearly improved estimation of high-frequency energy level (in terms of saturation spectrum and mean square slope). In the duration-limited test, we investigate the modeled wave spectrum in a detailed way by introducing spectral metrics for the tail and the peak of the omnidirectional wave spectrum and for the directionality of the two-dimensional frequency–direction spectrum. The omnidirectional frequency spectrum E ( f ) from the recalibrated ST6 shows a clear transition behavior from a power law of approximately f ⁻⁴ to a power law of about f ⁻⁵ , comparable to previous field studies. Different solvers for nonlinear wave interactions are applied with ST6, including the Discrete Interaction Approximation (DIA), the more expensive Generalized Multiple DIA (GMD), and the very expensive exact solutions [using the Webb–Resio–Tracy method (WRT)]. The GMD-simulated E ( f ) is in excellent agreement with that from WRT. Nonetheless, we find the peak of E ( f ) modeled by the GMD and WRT appears too narrow. It is also shown that in the 1-yr global hindcast, the DIA-based model overestimates the low-frequency wave energy (wave period T > 16 s) by 90%. Such model errors are reduced significantly by the GMD to ~20%.

The authors present airborne observations of fetch-limited waves during strong offshore winds in the Gulf of Tehuantepec. The measurements, collected over a wide range of fetches, include one- and two-dimensional surface wavenumber spectra and turbulent fluxes in winds up to 25 m s−1. The evolution of the wave spectra is in good agreement with the fetch relationships from previous observations. The tails of the observed one-dimensional k1 spectra, in the dominant wave direction, exhibit a power law over a wide range of wavenumbers. The authors present the first quantification of the transition between the equilibrium and saturation ranges for the omnidirectional spectrum in the wavenumber domain. The energy density within the equilibrium range shows a dependence on the wave age. At wavenumbers higher than the spectral peak, the width of the spectrum in the direction orthogonal to the dominant waves is nearly independent of the wave age. However, in the azimuthal direction, the spreading of the spectrum decreases with increasing effective wave age. The bimodal directional distribution, characterized by the lobe amplitude and separation, shows a consistent scaling with both parameters collapsing when scaled by the square root of the wave age. The one-dimensional fetch-limited k1 spectrum is well parameterized with dependence on the effective fetch and friction velocity. At higher wavenumbers within the saturation range, although the one-dimensional saturation in the dominant wave direction is independent of the wind forcing, the saturation in the crosswind direction is weakly dependent on the effective wave age and on average 30% larger than that in the downwind direction. The results are discussed in the context of previous observations and current numerical wind-wave prediction models.

A semiempirical determination of the spectral dependence of the energy dissipation due to surface wave breaking is presented and then used to propose a model for the spectral dependence of the breaking strength parameter b, defined in the O. M. Phillips’s statistical formulation of wave breaking dynamics. The determination of the spectral dissipation is based on closing the radiative transport equation for fetch-limited waves, measured in the Gulf of Tehuantepec Experiment, by using the measured evolution of the directional spectra with fetch, computations of the four-wave resonant interactions, and three models of the wind input source function. The spectral dependence of the breaking strength is determined from the Kleiss and Melville measurements of the breaking statistics and the semiempirical spectral energy dissipation, resulting in b = b(k, cp/u*), where k is the wavenumber and the parametric dependence is on the wave age, cp/u*. Guided by these semiempirical results, a model for b(k, cp/u*) is proposed that uses laboratory data from a variety of sources, which can be represented by b = a(S − S0)n, where S is a measure of the wave slope at breaking, a is a constant, S0 is a threshold slope for breaking, and 2.5 < n < 3 is a power law consistent with inertial wave dissipation scaling and laboratory measurements. The relationship between b(S) in the laboratory and b(k) in the field is based on the relationship between the saturation and mean square slope of the wave field. The results are discussed in the context of wind wave modeling and improved measurements of breaking in the field.

Wave–current interaction can result in significant inhomogeneities of the ocean surface wave field, including modulation of the spectrum, wave breaking rates, and wave statistics. This study presents novel airborne observations from two experiments: 1) the High-Resolution Air–Sea Interaction (HiRes) experiment, with measurements across an upwelling jet off the coast of Northern California, and 2) an experiment in the Gulf of Mexico with measurements of waves interacting with the Loop Current and associated eddies. The significant wave height and slope varies by up to 30% because of these interactions at both sites, whereas whitecap coverage varies by more than an order of magnitude. Whitecap coverage is well correlated with spectral moments, negatively correlated with the directional spreading, and positively correlated with the saturation. Surface wave statistics measured in the Gulf of Mexico, including wave crest heights and lengths of crests per unit surface area, show good agreement with second-order nonlinear approximations, except over a focal area. Similarly, distributions of wave heights are generally bounded by the generalized Boccotti distribution, except at focal regions where the wave height distribution reaches the Rayleigh distribution with a maximum wave height of 2.55 times the significant wave height, which is much larger than the standard classification for extreme waves. However, theoretical distributions of spatial statistics that account for second-order nonlinearities approximately bound the observed statistics of extreme wave elevations. The results are discussed in the context of improved models of breaking and related air–sea fluxes.

A large-eddy simulation (LES) model, which adopts wave-averaged equations with vortex force, is used to investigate Langmuir turbulence and ocean boundary layer (OBL) dynamics in high-wind hurricane conditions. The temporally evolving spatially asymmetric wind and wave Stokes drift velocity imposed in the LES are generated by a spectral wave prediction model adapted to Hurricane Frances traveling at a speed of 5.5 m s−1. The potency of Langmuir turbulence depends on the turbulent Langmuir number, the wind–Stokes drift alignment, and the depth scale of the Stokes profile Ds relative to the OBL depth h. At the time of maximum winds, large-scale vigorous coherent cells develop on the right-hand side of the storm under the inertially rotating winds; the Stokes drift velocity is well tuned to the surface winds. Much weaker cells develop on the left-hand side of the storm, partly because of reduced Stokes production. With misaligned winds and waves the vertical momentum fluxes can be counter to the gradient of Stokes drift, and the cell orientation tracks the direction of the mean Lagrangian shear. The entrainment flux is increased by 20% and the sea surface temperature is 0.25 K cooler on the right-hand side of the storm in the presence of Langmuir turbulence. Wave effects impact entrainment when the ratio Ds/|h| > 0.75. Because of wind–wave asymmetry Langmuir cells add quantitatively to the left–right asymmetry already understood for hurricanes due to resonance. And the transient evolution of the OBL cannot be understood simply in terms of equilibrium snapshots.

This study describes a model of Phillips' Λ(c) distribution, which is the expected length of breaking fronts (per unit surface area) moving with velocity c to c+dc, providing a framework for coupled atmosphere‐wave‐ocean models to explicitly account for wave breaking related air‐sea fluxes. The model of Λ depends on the spectral saturation, based on Gaussian statistics of the lengths of crest exceeding wave slope criteria, and long wave‐short wave modulation. A wave breaking dissipation function based on Λ was implemented in the model WaveWatchIII. The wave solutions are consistent with the observations, including several metrics of the spectrum and Λ(c) distributions. The whitecap coverage derived from Λ reproduces recent parameterizations saturating at high winds. The wave breaking variability due to wave‐current interaction is significant at submesoscales (order 1 km or smaller). The wave breaking model can be further developed to model gas transfer coefficients and aerosol production.