Search

It is generally understood that the origin of ocean diapycnal diffusivity is primarily associated with the stratified turbulence produced by breaking internal (gravity) waves (IW). However, it requires significant effort to verify diffusivity values in ocean general circulation models in any particular geographical region of the ocean due to the scarcity of microstructure measurements. Recent analyses of downscaled IW fields from an internal‐wave‐admitting global ocean simulation into higher‐resolution regional configurations northwest of Hawaii have demonstrated a much‐improved fit of the simulated IW spectra to the in‐situ profiler measurements such as the Garrett‐Munk (GM) spectrum. Here, we employ this dynamically downscaled ocean simulation to directly analyze the nature of the IW‐breaking and the wave‐turbulence cascade in this region. We employ a modified version of the Kappa Profile Parameterization (KPP) to infer what the horizontally averaged vertical profile of diapycnal diffusivity should be, and compare this to the background profile that would be employed in the ocean component of a low‐resolution coupled climate model such as the Community Earth System Model (CESM) of the US National Center for Atmospheric Research (NCAR). In pursuing this goal, we also demonstrate that the wavefield in the high‐resolution regional domain is dominated by a well‐resolved spectrum of low‐mode IWs that are predictable by solving an appropriate eigenvalue problem for stratified flow. We finally suggest a new tentative approach to improve the KPP parameterization. A much‐improved spectrum of the simulated internal wave (IW) field has recently been obtained by downscaling a global ocean model into a higher‐resolution regional configuration. The global simulation is based on the Massachusetts Institute of Technology general circulation model (MITgcm) forced by both astronomical tidal potential and surface atmospheric processes. By employing a mathematical framework to predict the structure of IWs, we first demonstrate that the interior wavefield of the high‐resolution regional domain is well dominated by a series of low‐order IW modes. Then, we address the issue as to whether the component of the K‐Profile Parameterization (KPP) associated with IW shear might be able to explain the physical origins of the background depth dependence of diapycnal diffusivity that would normally be employed in the ocean component of a modern coupled climate model. Finally, we suggest a tentative approach to further improve KPP. A representation of internal wave modes that includes dissipation is derived that explains spectra in a tidally forced numerical simulationEliminating the background component of KPP leads to a representation of the physics of internal wave breakingRealistic diapycnal diffusivity profiles can be obtained by minor adjustments to the shear component of KPP A representation of internal wave modes that includes dissipation is derived that explains spectra in a tidally forced numerical simulation Eliminating the background component of KPP leads to a representation of the physics of internal wave breaking Realistic diapycnal diffusivity profiles can be obtained by minor adjustments to the shear component of KPP

This paper examines spectra of horizontal kinetic energy (HKE) in the surface and sub‐surface ocean, with an emphasis on internal gravity wave (IGW) motions, in global high‐resolution ocean simulations. Horizontal wavenumber‐frequency spectra of surface HKE are computed over seven oceanic regions from two global simulations of the HYbrid Coordinate Ocean Model (HYCOM) and three global simulations of the Massachusetts Institute of Technology general circulation model (MITgcm). In regions with high IGW activity, high surface HKE variance in the horizontal wavenumber‐frequency spectra is aligned along IGW linear dispersion curves. For both HYCOM and MITgcm, and in almost all regions, finer horizontal resolution yields more energetic supertidal IGW continuum spectra. The ratio of high‐horizontal‐wavenumber variance in semi‐diurnal and supertidal motions relative to lower‐frequency motions, a quantity of great interest for swath altimetry, depends on the model employed and the horizontal resolution within the model, implying that quantitative predictions of the partition between low‐ and high‐frequency motions taken from particular simulations should be treated with care. The frequency‐vertical wavenumber spectra, frequency spectra, and vertical wavenumber spectra from the models are compared to spectra computed from McLane profilers at nine locations. In general, MITgcm spectra match the McLane profiler spectra more closely at high frequencies (|ω| > 4.5 cpd). In both models, vertical wavenumber spectra roll off more steeply than observations at high vertical wavenumbers (m> 10−2cpm). The vertical wavenumber spectra in such models is an important target for improvement, due to turbulence production and dissipation that takes place at high vertical wavenumbers. Recently, a small but growing number of global ocean models have begun to employ simultaneous tidal and atmospheric forcing. At the same time, increasing supercomputer power has allowed for simulations of oceanic motions with increasing accuracy, increasing feature (spatial) resolution, and more frequent time slices. Global ocean models with fine grid spacing, and simultaneous tidal and atmospheric forcing, host a vigorous spectrum of high‐frequency waves that control mixing over most of the ocean water column, and are important for many operational oceanography challenges. As an example of the latter, high‐resolution global internal wave models have been used to study the relative partition of high‐frequency versus low‐frequency motions at the small horizontal scales that will be measured by the new Surface Water Ocean Topography mission. The partition described above depends on the model employed and the grid spacing employed within that model, meaning that conclusions about the partition are dependent on the model used to estimate it. Comparisons between the models and vertically profiling instruments indicate that resolving fine scale motions in the vertical direction, where ocean mixing takes place, is not yet handled well by the models. Modeling of fine‐vertical scale motions is therefore an important future research direction. Vertical wavenumber spectra of internal gravity wave kinetic energy in two high‐resolution global models are compared to observed spectraModels under‐estimate motions at high vertical wavenumbers (small vertical scales), flagging this as a target for model improvementThe ratio of high‐ versus low‐frequency surface kinetic energy at small horizontal scales is dependent on the model and grid spacings employed Vertical wavenumber spectra of internal gravity wave kinetic energy in two high‐resolution global models are compared to observed spectra Models under‐estimate motions at high vertical wavenumbers (small vertical scales), flagging this as a target for model improvement The ratio of high‐ versus low‐frequency surface kinetic energy at small horizontal scales is dependent on the model and grid spacings employed

Wake eddies are frequently created by flow separation where ocean currents encounter abrupt topography in the form of islands or headlands. Most previous work has concentrated on wake eddy generation by either purely oscillatory (usually tidal) currents, or quasi‐steady mean flows. Here we report measurements near the point of flow separation at the northern end of the Palau island chain, where energetic tides and vertically sheared low‐frequency flows are both present. Energetic turbulence measured near the very steeply sloping ocean floor varied cubically with the total flow speed (primarily tidal). The estimated turbulent viscosity suggests a regime of flow separation and eddying wake generation for flows that directly feel this drag. Small‐scale (∼1 km), vertically sheared wake eddies of different vorticity signs were observed with a shipboard survey on both sides of the separation point, and significantly evolved over several tidal periods. The net production and export of vorticity into the wake, expected to sensitively depend on the interplay of tidal and low‐frequency currents, is explored here with a simple conceptual model. Application of the model to a 10‐month mooring record suggests that inclusion of high frequency oscillatory currents may boost the net flux of vorticity into the ocean interior by a depth dependent factor of 2 to 25. Models that do not represent the effect of these high frequency currents may not accurately infer the net momentum or energy losses felt where strong flows encounter steep island or headland topography. Island wake generation is influenced by both low‐frequency and tidal flows

The generation, propagation, and dissipation of nonlinear internal waves (NLIW) in sheared background currents is examined using 7 days of shipboard microstructure surveys and two moorings on the continental shelf offshore of Washington state. Surveys near the hypothesized generation region show semi‐diurnal (D2) energy flux is onshore and that the ratio of energy flux to group speed times energy ( F/cgE) increases sharply at the shelf break, suggesting that the incident D2 internal tide is partially reflected and partially transmitted. NLIW appear at an inshore mooring at the leading edge of the onshore phase of the baroclinic tide, consistent with nonlinear transformation of the shoaling internal tide as their generation mechanism. Of the D2 energy flux observed at the eastern extent of the generation region (133 ± 18 Wm−1), approximately 30% goes into the NLIW observed inshore (36 ± 11 Wm−1). Inshore of the moorings, 7 waves are tracked into shallow (30‐40 m) water, where a vertically sheared, southward current becomes strong. As train‐like waves propagate onshore, wave amplitudes of 25‐30 m and energies of 5 MJ decrease to 12 m and 10 kJ, respectively. The observed direction of propagation rotates from 30°N of E to ∼30°S of E in the strongly sheared region. Linear ray tracing using the Taylor‐Goldstein equation to incorporate parallel shear effects accounts for only a small portion of the observed rotation, suggesting that three‐dimensionality of the wave crests and the background currents is important here. Generation of NLIW is observed at the shelf break where the incident internal tide is partially reflected and partially transmittedApproximately 30% of the transmitted internal tide energy is transformed into NLIWObserved NLIW are either train‐like or bore‐like. Trains rotate CW up to 90 degrees and lose energy as they move onshore; bores do not

Tropical cyclones (TCs) are powered by heat fluxes across the air‐sea interface, which are in turn influenced by subsurface physical processes that can modulate storm intensity. Here, we use data from 6 profiling floats to recreate 3D fields of temperature (T), salinity (S), and velocity around the fast‐moving Super Typhoon Mangkhut (western North Pacific, September 2018). Observational estimates of vorticity (ζ) and divergence (Γ) agree with output from a 3D coupled model, while their relation to vertical velocities is explained by a linear theoretical statement of inertial pumping. Under this framework, inertial pumping is described as a linear coupling between ζand Γ, whose oscillations in quadrature cause periodic displacements in the ocean thermocline and generate near‐inertial waves (NIWs). Vertical profiles of Tand Sshow gradual mixing of the upper ocean with diffusivities as high as κ∼ 10−1m2s−1, which caused an asymmetric cold wake of sea surface temperature (SST). We estimate that ∼10% of the energy used by mixing was used to mix rainfall, therefore inhibiting SST cooling. Lastly, watermass transformation analyses suggest that κ> 3 × 10−3m2s−1above ∼110 m depth and up to 600 km behind the TC. These analyses provide an observational summary of the ocean response to fast‐moving TCs, demonstrate some advantages of ζand Γ for the study of internal wave fields, and provide conceptual clarity on the mechanisms that lead to NIW generation by winds. Near‐inertial internal waves (NIWs) are generated by winds and lead to oscillations in the internal structure of ocean currents and stratification. Turbulence induced by the vertical current shear in these waves helps sustain the upper ocean stratification and circulation. In this study, we use data from six autonomous floats deployed ahead of Super Typhoon Mangkhut to reconstruct the ocean's 3D response and compare it to output from a coupled air‐sea model. Patterns in NIW are explained using simple linear equations based on vorticity and divergence rather than current velocities, providing an alternative view of how TC winds help generate waves in the stratified ocean interior. Measurements of temperature and salinity detail how turbulence mixed rainfall and thermocline waters into the upper ocean. Our analyses indicate that turbulent mixing rates are greatest within 100 km of the typhoon eye but remain elevated hundreds of kilometers in the TC wake. Theory and observations presented here provide a comprehensive view of the ocean response to fast‐moving tropical cyclones. Float data, linear theory, and a 3D model reveal vorticity and divergence control on inertial pumping beneath a fast‐moving tropical cycloneRightward‐enhanced sea surface cooling of 1.2°C was dominated by mixing and modulated by rainfall, which suppressed cold water entrainmentEstimates of turbulent diffusivity explain sea surface cooling rates 0.1°C hr−1under the tropical cyclone eye and thermocline mixing in its wake Float data, linear theory, and a 3D model reveal vorticity and divergence control on inertial pumping beneath a fast‐moving tropical cyclone Rightward‐enhanced sea surface cooling of 1.2°C was dominated by mixing and modulated by rainfall, which suppressed cold water entrainment Estimates of turbulent diffusivity explain sea surface cooling rates 0.1°C hr−1under the tropical cyclone eye and thermocline mixing in its wake

Two global ocean models ranging in horizontal resolution from 1/12° to 1/48° are used to study the space and time scales of sea surface height (SSH) signals associated with internal gravity waves (IGWs). Frequency‐horizontal wavenumber SSH spectral densities are computed over seven regions of the world ocean from two simulations of the HYbrid Coordinate Ocean Model (HYCOM) and three simulations of the Massachusetts Institute of Technology general circulation model (MITgcm). High wavenumber, high‐frequency SSH variance follows the predicted IGW linear dispersion curves. The realism of high‐frequency motions (> 0.87cpd) in the models is tested through comparison of the frequency spectral density of dynamic height variance computed from the highest‐resolution runs of each model (1/25° HYCOM and 1/48° MITgcm) with dynamic height variance frequency spectral density computed from nine in situ profiling instruments. These high‐frequency motions are of particular interest because of their contributions to the small‐scale SSH variability that will be observed on a global scale in the upcoming Surface Water and Ocean Topography (SWOT) satellite altimetry mission. The variance at supertidal frequencies can be comparable to the tidal and low‐frequency variance for high wavenumbers (length scales smaller than ∼50 km), especially in the higher‐resolution simulations. In the highest‐resolution simulations, the high‐frequency variance can be greater than the low‐frequency variance at these scales. Two high‐resolution ocean models compare well against data in frequency spectral density of dynamic heightSea surface height frequency‐horizontal wavenumber spectral densities show high variance along internal gravity wave dispersion curvesTwo high‐resolution ocean models give different estimates of variance in high‐frequency, high wavenumber phenomena

We examine the temporal means and variability of the semidiurnal internal tide energy fluxes in 1/25° global simulations of the Hybrid Coordinate Ocean Model (HYCOM) and in a global archive of 79 historical moorings. Low‐frequency flows, a major cause of internal tide variability, have comparable kinetic energies at the mooring sites in model and observations. The computed root‐mean‐square (RMS) variability of the energy flux is large in both model and observations and correlates positively with the time‐averaged flux magnitude. Outside of strong generation regions, the normalized RMS variability (the RMS variability divided by the mean) is nearly independent of the flux magnitudes in the model, and of order 23% or more in both the model and observations. The spatially averaged flux magnitudes in observations and the simulation agree to within a factor of about 1.4 and 2.4 for vertical mode‐1 and mode‐2, respectively. The difference in energy flux computed from the full‐depth model output versus model output subsampled at mooring instrument depths is small. The global historical archive is supplemented with six high‐vertical resolution moorings from the Internal Waves Across the Pacific (IWAP) experiment. The model fluxes agree more closely with the high‐resolution IWAP fluxes than with the historical mooring fluxes. The high variability in internal tide energy fluxes implies that internal tide fluxes computed from short observational records should be regarded as realizations of a highly variable field, not as “means” that are indicative of conditions at the measurement sites over all time. The RMS variability in internal tide energy flux is large in both model and observationsMode‐1 internal tide fluxes compare much better with observationsHYCOM compares very well with internal tides fluxes from the IWAP experiment

High horizontal‐resolution ( 1/12.5°and 1/25°) 41‐layer global simulations of the HYbrid Coordinate Ocean Model (HYCOM), forced by both atmospheric fields and the astronomical tidal potential, are used to construct global maps of sea surface height (SSH) variability. The HYCOM output is separated into steric and nonsteric and into subtidal, diurnal, semidiurnal, and supertidal frequency bands. The model SSH output is compared to two data sets that offer some geographical coverage and that also cover a wide range of frequencies—a set of 351 tide gauges that measure full SSH and a set of 14 in situ vertical profilers from which steric SSH can be calculated. Three of the global maps are of interest in planning for the upcoming Surface Water and Ocean Topography (SWOT) two‐dimensional swath altimeter mission: (1) maps of the total and (2) nonstationary internal tidal signal (the latter calculated after removing the stationary internal tidal signal via harmonic analysis), with an average variance of 1.05and 0.43cm2, respectively, for the semidiurnal band, and (3) a map of the steric supertidal contributions, which are dominated by the internal gravity wave continuum, with an average variance of 0.15cm2. Stationary internal tides (which are predictable), nonstationary internal tides (which will be harder to predict), and nontidal internal gravity waves (which will be very difficult to predict) may all be important sources of high‐frequency “noise” that could mask lower frequency phenomena in SSH measurements made by the SWOT mission. Internal gravity waves and nonstationary internal tides have nonnegligible sea surface height signaturesHigh‐resolution global ocean models are beginning to resolve the internal gravity wave continuumInternal gravity waves and nonstationary internal tides may cause contamination in upcoming satellite altimeter missions

The geographical variability, frequency content, and vertical structure of near‐surface oceanic kinetic energy (KE) are important for air‐sea interaction, marine ecosystems, operational oceanography, pollutant tracking, and interpreting remotely sensed velocity measurements. Here, KE in high‐resolution global simulations (HYbrid Coordinate Ocean Model; HYCOM, and Massachusetts Institute of Technology general circulation model; MITgcm), at the sea surface (0 m) and at 15 m, are compared with KE from undrogued and drogued surface drifters, respectively. Global maps and zonal averages are computed for low‐frequency (<0.5 cpd), near‐inertial, diurnal, and semidiurnal bands. Both models exhibit low‐frequency equatorial KE that is low relative to drifter values. HYCOM near‐inertial KE is higher than in MITgcm, and closer to drifter values, probably due to more frequently updated atmospheric forcing. HYCOM semidiurnal KE is lower than in MITgcm, and closer to drifter values, likely due to inclusion of a parameterized topographic internal wave drag. A concurrent tidal harmonic analysis in the diurnal band demonstrates that much of the diurnal flow is nontidal. We compute simple proxies of near‐surface vertical structure—the ratio 0 m KE/(0 m KE + 15 m KE) in model outputs, and the ratio undrogued KE/(undrogued KE + drogued KE) in drifter observations. Over most latitudes and frequency bands, model ratios track the drifter ratios to within error bars. Values of this ratio demonstrate significant vertical structure in all frequency bands except the semidiurnal band. Latitudinal dependence in the ratio is greatest in diurnal and low‐frequency bands. It is important to map and understand ocean surface currents because they affect climate and marine ecosystems. Recent advances in global ocean models include the addition of astronomical tidal forcing alongside atmospheric forcing and the usage of more powerful computers that can resolve finer features. Here, we evaluate ocean surface currents in high‐resolution simulations of two different ocean models through comparison with observations from surface drifting buoys. We examine near‐inertial motions, forced by fast‐changing winds; semidiurnal tides, forced by the astronomical tidal potential; diurnal motions, arising from tidal and other sources; and low‐frequency currents and eddies, forced by atmospheric fields. Global patterns in the models and drifters are broadly consistent. The two models differ in their degree of proximity to drifter measurements in the near‐inertial band, most likely due to different update intervals for atmospheric forcing and in the semidiurnal band, most likely due to different damping schemes. A simple proxy for vertical structure of the currents, measured by differences in drifter flows at the surface versus 15 m depth, is tracked reasonably well by the models. Discrepancies between models and observations motivate future improvements in the models. We examine frequency content of ocean kinetic energy (KE) at the sea surface (0 m) and 15 m depth with global drifter data and two modelsNear‐surface near‐inertial and tidal KE in numerical models are sensitive to atmospheric forcing frequency and dampingThe ratio 0 m KE/(0 m KE + 15 m KE) in models lies within error bars of the observational ratios over some latitudes and frequency bands We examine frequency content of ocean kinetic energy (KE) at the sea surface (0 m) and 15 m depth with global drifter data and two models Near‐surface near‐inertial and tidal KE in numerical models are sensitive to atmospheric forcing frequency and damping The ratio 0 m KE/(0 m KE + 15 m KE) in models lies within error bars of the observational ratios over some latitudes and frequency bands

Tidal and low‐frequency flows interact with abrupt topographic features, giving rise to a variety of small‐scale phenomena including lee waves, vortical motions and turbulence. With the goal of identifying the processes that transfer momentum and energy between the larger‐scale and smaller‐scale flows, detailed tidally‐resolving shipboard surveys were conducted near a sharp and steep headland north of Velasco Reef, Palau, in the Western Pacific. We present velocity, microstructure and rapid‐cycling conductivity and temperature measurements, contextualized with a nearby 10‐month moored record, to examine wake and vortex structures. Observed flow near Velasco Reef was layered and multi‐directional in depth such that the “downstream” side of the headland depended on depth. Downstream of the headland, during a period where mean flow dominated over the tide in the top 100 m of the water column, turbulent kinetic energy dissipation rates were enhanced by up to three orders of magnitude relative to those upstream. During a period dominated by tides, the flow reversed periodically and wake eddies were associated with turbulent dissipation rates up to O10−5Wkg−1$\mathcal{O}\left(1{0}^{-5}\right)\ \mathrm{W}\hspace*{.5em}\mathrm{k}{\mathrm{g}}^{-\mathrm{1}}$. Turbulent dissipation rates in the wake increased with inverse Richardson number for both tidal and mean flows past the headland, and were larger for stronger vorticity flows. Our analysis suggests that the eddies themselves lead to shear instabilities and subsequent turbulence, either through tilting or their limited vertical extent. We associate the enhanced downstream dissipation with a combination of breaking internal tides and either these eddy processes or a multi‐scaled turbulent wake. Tall, steep underwater features such as ridges and headlands interact with tides and large currents to alter the flow downstream. Here, we present observations near Velasco Reef, Palau, an example of a sharp headland in the western Pacific. Turbulence in the wake of the headland is higher than upstream. During a time when the tides were strong the flow changed direction such that eddies were shed from both sides of the headland. These patches of rotating flow were associated with stronger turbulence than elsewhere around the headland. Shearing of the flow in the vertical direction by the eddies and the internal tides are possible mechanisms by which the turbulence is generated. Direct observations of turbulent dissipation rates show enhancement downstream of a headland during mean flow and tidal flow regimesRegions of strong vorticity in the wake are associated with the highest turbulent dissipation ratesShear instabilities are a likely mechanism for the generation of turbulence in the wake Direct observations of turbulent dissipation rates show enhancement downstream of a headland during mean flow and tidal flow regimes Regions of strong vorticity in the wake are associated with the highest turbulent dissipation rates Shear instabilities are a likely mechanism for the generation of turbulence in the wake

Two hydrographic surveys in the Gulf of Alaska and the North Pacific subtropical gyre are presented. Both surveys are roughly perpendicular to lateral temperature gradients, and were collected in the summer when there was a shallow mixed layer and a seasonal thermocline. Isopycnal displacements and horizontal velocities are dominated by internal waves. Spice anomalies along isopycnals are examined to diagnose lateral stirring mechanisms. The spectra of spice anomaly gradients along near‐surface isopycnals roughly follow power laws of ∼kx0.6±0.2(variance spectra power laws of ∼kx−1.4±0.2), and in most cases, the spectra become redder at depth. The near‐surface spectra are possibly consistent with the predictions of quasi‐geostrophic turbulence theory (when surface buoyancy effects are accounted for), but the spectra at depth are inconsistent with any quasi‐geostrophic theory. Probability distributions of spice gradients exhibit a large peak at low gradients and long tails for large gradients, symptomatic of fronts. Vertical coherence of the spice signal falls off with a decorrelation depth scale that has a maximum of about 80 m at 100 km wavelengths and depends on horizontal wavelength with a power law of approximately kx−1/2. Lateral decorrelation length scales are 20–40 km, close to the baroclinic Rossby radius. Lateral stirring occurs over large scales, with average lateral displacements of about 200 km in the upper 75 m, decreasing to 100 km at greater depths. The depth variation of the statistics indicates that time history of tracer stirring on each isopycnal is important, or that there are unconsidered depth‐dependent stirring mechanisms. Spice variance quantified along two sectionsQG/SQG theories not applicableStrong depth dependence

Internal tidal energy fluxes were obtained from June 2011 to August 2011 using underwater gliders in the South China Sea. Spraygliders profiled every ∼2 h to 500 m, which is deep enough given the shallow thermocline to compute mode‐1 fluxes from vertical mode fits to tidal displacements and currents. Westward, mode‐1 diurnal and semidiurnal fluxes exceeded 40 and 30 kW m−1. To our knowledge, these flux observations are the first from both velocity and density measurements by gliders. Fluxes compare well with a mooring near a generation site in southern Luzon Strait and a regional model. Furthermore, the zonal‐depth structure of the internal tide is obtained by binning measurements, which cover four spring‐neap cycles and over 100 km along 20°39′N. Westward phase propagation is found for currents and displacements, while roughly constant phase is found along beams. Both these features of the phase suggest a narrow‐banded internal tide. Semidiurnal energy density is largest along a raypath which coincides with generation sites on both the eastern and western ridges in Luzon Strait. Diurnal energy density is surface‐intensified consistent with relatively shallower diurnal raypaths emanating from the eastern ridge. Energy fluxes from both velocity and density measurements by glidersWestward phase propagation for currents and displacementsRoughly constant phase along beams

A near‐inertial wave stalling and breaking in a critical layer was observed for a week by a pair of autonomous velocity and density profilers on concentric paths in an ocean mesoscale eddy. Profiler observations provide estimates of the eddy's vertical vorticity and shear, quantities needed to test theories of downward near‐inertial wave (NIW) energy flux and loss from inertial wave–eddy interactions. The unique observations of the wave's intrinsic frequency ωiand vertical wavenumber mprovide a novel estimate of the vertical group velocity Cgzfrom changes in ωiwith respect to m. Shear, strain, energy flux convergence, and parameterized turbulence are all elevated near 140 m depth, near the bottom of the strongest eddy velocities. Our observations are consistent with a downgoing NIW's group velocity decreasing owing to wave–eddy interactions, providing important clues on global energetics of NIW mixing. Internal gravity waves with frequency near the Earth's Coriolis frequency, called a near‐inertial wave, was observed stalling and breaking as it propagated downward from the sea surface. The wave stalled because its downward propagation speed slowed in the presence of a background eddy that weakened with depth. We believe this is by far the clearest observation of such a process, which is an important aspect of interactions between internal gravity waves and eddies. A key finding of the paper is the first‐ever direct determination of the wave's propagation speed from our observations without assumptions or reliance on theory. A pair of EM‐APEX floats moving with a subsurface turbulent zone observe mesoscale vorticity (by Kelvin's circulation theorem) and vertical shear; observations resolve near‐inertial wave phase and group velocities, density, and shear to diagnose dynamics of NIW stalling and dissipationFirst known determination of NIW vertical group velocity computed from ∂ωi/∂m in ocean observationsAnalytic theory of mesoscale vorticity and vertical shear predicts NIW kinetic properties that compares well with observations A pair of EM‐APEX floats moving with a subsurface turbulent zone observe mesoscale vorticity (by Kelvin's circulation theorem) and vertical shear; observations resolve near‐inertial wave phase and group velocities, density, and shear to diagnose dynamics of NIW stalling and dissipation First known determination of NIW vertical group velocity computed from ∂ωi/∂m in ocean observations Analytic theory of mesoscale vorticity and vertical shear predicts NIW kinetic properties that compares well with observations

Temperature variance and kinetic energy (KE) from two global simulations of the HYbrid Coordinate Ocean Model (HYCOM; 1/12° and 1/25°) and three global simulations of the Massachusetts Institute of Technology general circulation model (MITgcm; 1/12°, 1/24°, and 1/48°), all of which are forced by atmospheric fields and the astronomical tidal potential, are compared with temperature variance and KE from a database of about 2,000 moored historical observations (MHOs). The variances are computed across frequencies ranging from supertidal, dominated by the internal gravity wave continuum, to subtidal, dominated by currents and mesoscale eddies. The most important qualitative difference between HYCOM and MITgcm, and between simulations of different resolutions, is in the supertidal band, where the 1/48° MITgcm lies closest to observations. Across all frequency bands examined, the HYCOM simulations display higher spatial correlation with the MHO than do the MITgcm simulations. The supertidal, semidiurnal, and diurnal velocities in the HYCOM simulations also compare more closely with observations than do the MITgcm simulations in a number of specific continental margin/marginal sea regions. To complement the model‐MHO comparisons, this paper also compares the surface ocean geostrophic eddy KE in HYCOM, MITgcm, and a gridded satellite altimeter product. Consistent with the model‐MHO comparisons, the HYCOM simulations have a higher spatial correlation with the altimeter product than the MITgcm simulations do. On the other hand, the surface ocean geostrophic eddy KE is too large, relative to the altimeter product, in the HYCOM simulations. Kinetic energy and temperature variance in global HYCOM and MITgcm simulations are compared with moored observationsModel resolution has large impact on the supertidal internal gravity wave continuum bandSpatial correlations with observations are higher in HYCOM than in MITgcm