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. Plain Language Summary: 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. Key Points: 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 [ABSTRACT FROM AUTHOR]

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−2 cpm). 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. Plain Language Summary: 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. Key Points: 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 [ABSTRACT FROM AUTHOR]

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 T and S show gradual mixing of the upper ocean with diffusivities as high as κ ∼ 10−1 m2 s−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−3 m2 s−1 above ∼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. Plain Language Summary: 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. Key Points: 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−1 under the tropical cyclone eye and thermocline mixing in its wake [ABSTRACT FROM AUTHOR]

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. Plain Language Summary: 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. Key Points: 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 [ABSTRACT FROM AUTHOR]

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. Plain Language Summary: 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. Key Points: 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 [ABSTRACT FROM AUTHOR]

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 ωi and vertical wavenumber m provide a novel estimate of the vertical group velocity Cgz from changes in ωi with 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. Plain Language Summary: 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. Key Points: 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 [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. Key Point: Island wake generation is influenced by both low‐frequency and tidal flows [ABSTRACT FROM AUTHOR]

Abstract: 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 / c g E) 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. [ABSTRACT FROM AUTHOR]

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 (> [ABSTRACT FROM AUTHOR]