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]

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]

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]

Unprecedented quantities of heat are entering the Pacific sector of the Arctic Ocean through Bering Strait, particularly during summer months. Though some heat is lost to the atmosphere during autumn cooling, a significant fraction of the incoming warm, salty water subducts (dives beneath) below a cooler fresher layer of near-surface water, subsequently extending hundreds of kilometers into the Beaufort Gyre. Upward turbulent mixing of these sub-surface pockets of heat is likely accelerating sea ice melt in the region. This Pacific-origin water brings both heat and unique biogeochemical properties, contributing to a changing Arctic ecosystem. However, our ability to understand or forecast the role of this incoming water mass has been hampered by lack of understanding of the physical processes controlling subduction and evolution of this this warm water. Crucially, the processes seen here occur at small horizontal scales not resolved by regional forecast models or climate simulations; new parameterizations must be developed that accurately represent the physics. Here we present novel high resolution observations showing the detailed process of subduction and initial evolution of warm Pacific-origin water in the southern Beaufort Gyre. Warming ocean water plays a significant role in accelerating Arctic sea ice melt. Here the authors present detailed observations of warm water of Pacific origin entering and diving beneath the Arctic ocean surface, and explore the dynamical processes governing its evolution. [ABSTRACT FROM AUTHOR]

Despite sufficient wind forcing, internal waves in the South China Sea do not exhibit the strong near-inertial wave (NIW) peak that is typical in most of the world oceans. Using data from 10 contemporaneous moorings deployed in summer 2011, we show that strong isopycnal vertical tidal displacements transfer most of the near-inertial (NI) kinetic energy (KE) to frequencies higher than the inertial frequency in an Eulerian reference frame. Transforming to an isopycnal-following reference frame increases the KE at NI frequencies, suggesting the presence of NIWs. However, the projection onto a semi-Lagrangian coordinate system still underestimates the expected NI peak. To fully resolve NIWs requires the use of time-dependent vertical wavenumber–frequency spectra because the intrinsic frequency of the NIWs varies substantially, owing to Doppler shifting by lateral mesoscale flows. Here, we show NIW intrinsic frequency variations of ±0.2 cpd within few days, of similar magnitude as the observed variations of relative vorticity associated with the meandering Kuroshio. [ABSTRACT FROM AUTHOR]

The spatial, temporal, and directional characteristics of shear are examined in the upper 1400 m of the North Pacific during late spring with an array of five profiling moorings deployed from 25° to 37°N (1330 km) and simultaneous shipboard transects past them. The array extended from a regime of moderate wind generation at the north to south of the critical latitude 28.8°N, where parametric subharmonic instability (PSI) can transfer energy from semidiurnal tides to near-inertial motions. Analyses are done in an isopycnal-following frame to minimize contamination by Doppler shifting. Approximately 60% of RMS shear at vertical scales >20m (and 80% for vertical scales >80 m) is contained in near-inertial motions. An inertial back-rotation technique is used to index shipboard observations to a common time and to compute integral time scales of the shear layers. Persistence times are O(7) days at most moorings but O(25) days at the critical latitude. Simultaneous shipboard transects show that these shear layers can have lateral scales ≥100 km. Layers tend to slope downward toward the equator north of the critical latitude and are more flat to its south. Phase between shear and strain is used to infer lateral propagation direction. Upgoing waves are everywhere laterally isotropic. Downgoing waves propagate predominantly equatorward north and south of the critical latitude but are isotropic near it. Broadly, results are consistent with wind generation north of the critical latitude and PSI near it-and suggest a more persistent and laterally coherent near-inertial wave field than previously thought. [ABSTRACT FROM AUTHOR]

Shipboard ADCP and towed CTD measurements are presented of a near-inertial internal gravity wave radiating away from a zonal jet associated with the Subtropical Front in the North Pacific. Three-dimensional spatial surveys indicate persistent alternating shear layers sloping downward and equatorward from the front. As a result, depth-integrated ageostrophic shear increases sharply equatorward of the front. The layers have a vertical wavelength of about 250 m and a slope consistent with a wave of frequency 1.01 f. They extend at least 100 km south of the front. Time series confirm that the shear is associated with a downward-propagating near-inertial wave with frequency within 20% of f. A slab mixed layer model forced with shipboard and NCEP reanalysis winds suggests that wind forcing was too weak to generate the wave. Likewise, trapping of the near-inertial motions at the low-vorticity edge of the front can be ruled out because of the extension of the features well south of it. Instead, the authors suggest that the wave arises from an adjustment process of the frontal flow, which has a Rossby number about 0.2-0.3. [ABSTRACT FROM AUTHOR]

The downward propagation of near-inertial internal waves following winter storms is examined in the context of a 2-yr record of velocity in the upper 800 m at Ocean Station Papa. The long time series allow accurate estimation of wave frequency, whereas the continuous data in depth allow separation into upward- and downward-propagating components. Near-inertial kinetic energy (KEin) dominates the record. At all measured depths, energy in downgoing motions exceeds that of upward-propagating motions by factors of 3-7, whereas KEin is elevated by a factor of 3-5 in winter relative to summer. The two successive winters are qualitatively similar but show important differences in timing and depth penetration. Energy is seen radiating downward in a finite number of wave groups, which are tagged and catalogued to determine the vertical group velocity cgz, which has a mean of about 1.5 × 10−4 m s−1 (13 m day−1). Case studies of three of these are presented in detail. Downward energy flux is estimated as cgz × KEin (i) by summing over the set of events, (ii) from time series near the bottom of the record, and (iii) from the wavenumber-frequency spectrum and the dispersion relationship. These estimates are compared to the work done on near-inertial motions in the mixed layer by the wind, which is directly estimated from mixed layer near-inertial currents and winds measured from a surface buoy 10 km away. All three methods yield similar values, indicating that 12%-33% of the energy input into the mixed layer transits 800 m toward the deep sea. This simple picture neglects lateral energy flux carried by the first few vertical modes, which was not measured. The substantial deep penetration implies that near-inertial motions may play a role in mixing the deep ocean, but the strong observed variability calls for a need to better understand the role of lateral mesoscale structures in modulating the vertical propagation. [ABSTRACT FROM AUTHOR]

Observations are reported of the semidiurnal ( M2) internal tide across Kaena Ridge, Hawaii. Horizontal velocity in the upper 1000-1500 m was measured during eleven ~240-km-long shipboard acoustic Doppler current profiler (ADCP) transects across the ridge, made over the course of several months. The M2 motions are isolated by means of harmonic analysis and compared to numerical simulations using the Princeton Ocean Model (POM). The depth coverage of the measurements is about 3 times greater than similar past studies, offering a substantially richer view of the internal tide beams. Sloping features are seen extending upward north and south from the ridge and then downward from the surface reflection about ±40 km from the ridge crest, closely matching theoretical M2 ray paths and the model predictions. [ABSTRACT FROM AUTHOR]

The relative strength and spatiotemporal structure of near-inertial waves (NIW) and internal tides (IT) are examined in the context of recent moored observations made 19 km south of Mendocino Escarpment, an abrupt ridge/step feature in the eastern Pacific. In addition to strong internal tide generation, steps and ridges give rise to the possibility of ''shadowing,'' wherein near-inertial energy is prevented from reaching depths beneath a characteristic intersecting the ridge top. A combination of two moored profilers and a long-range acoustic Doppler current profiler (ADCP) yielded velocity and shear measurements from 100 to 3640 m (60 m above bottom) and isopycnal depth, strain, and overturn-inferred turbulence dissipation rate from 1000 to 3640 m. Sampling intervals (20 min in the upper 1000 m and 1.5 h below that) were fast enough to minimize aliasing of higher-frequency internal-wave motions. The 67-day-long record is easily sufficient to isolate NIW and IT via bandpass filtering and to capture low-frequency variability in all quantities. No near-inertial shadowing was observed. Instead, energetic near-inertial waves were observed at all depths, radiating both upward and downward. A strong upward internal tide beam, showing a pronounced spring--neap cycle, was also seen near the expected depth. Case studies of each of these are presented in depth and isopycnal-following coordinates. Except for immediately above the bottom and in the ''beam,'' where IT kinetic energy shows marked peaks, kinetic energy in the two bands is within a factor of 2 of each other. However, because of the redder NIW vertical wavenumber spectrum, NIW shear exceeded IT shear at all depths by a factor of 2--4. Dissipation rate was strongly enhanced in the bottom 1000 m and in the depth range of the internal tide beam. However, except for very near the bottom and possibly in one NIW event, no clear phase relationship was observed between dissipation rate and wave shear or strain, suggesting that turbulence occurs through a cascade process rather than by direct breaking at most locations. [ABSTRACT FROM AUTHOR]

Westward-propagating atmospheric easterly waves contribute to much of the variability of the low-level wind fields within the northeast tropical Pacific. With the dominant period of these waves (3–5 days) close to the local inertial period (2.4 days at 12°N to 5.7 days at 5°N), there is the expectation that the associated winds may resonantly force mixed layer inertial motions in this region. The authors test this hypothesis using a simple slab model and roughly 4½ yr of wind data from four NOAA Tropical Atmosphere Ocean/Eastern Pacific Investigation of Climate Processes (TAO/EPIC) buoys along 95°W at 12°, 10°, 8°, and 5°N. The degree of resonance is determined by comparing model simulations using observed wind stress with simulations forced with reversed-rotation wind stress. Results strongly indicate that Pacific easterly waves (PEWs) resonantly force inertial motions in the region. This resonance shows both significant seasonality and latitudinal dependence that appears to be related to the meridional position and intensity of the PEWs. North of the zonal axis of the mean track of the PEWs, the low-level winds associated with the waves rotate predominantly clockwise with time and resonantly force mixed layer inertial motions. South of this axis, the winds rotate counterclockwise, resulting in dissonant (antiresonant) forcing. As this axis migrates annually from about 4°N during the boreal winter/spring to a maximum northerly position of about 8°–10°N in the late boreal summer/early fall, the region of strongest resonance follows, consistently remaining to its north. Model output suggests that resonant forcing results in roughly 10%–25% greater net annual flux of kinetic energy from the wind to mixed layer inertial motions than in neutral or nonresonant conditions. This finding has strong implications for mixed layer properties, air–sea coupling, and the generation of near-inertial internal waves. [ABSTRACT FROM AUTHOR]

Temporal and spatial patterns of near-inertial kinetic energy (KEmoor) are investigated in a database of 2480 globally distributed, moored current-meter records (deployed on 690 separate moorings) and compared with the distribution of wind-forced mixed-layer energy flux FML. By computing KEmoor using short (30 day) multitaper spectral windows, the seasonal cycle is resolved. Clear winter enhancement by a factor of 4–5 is seen in the Northern Hemisphere for latitudes 25°–45° at all depths <4500 m, in close agreement with the magnitude, phase, and latitudinal dependence of the seasonal cycle of FML. In the Southern Hemisphere, data coverage is poorer, but a weaker seasonal cycle (a factor of 2) in both KEmoor and FML is still resolvable between 35° and 50°. When Wentzel–Kramers–Brillouin (WKB) scaled using climatological buoyancy-frequency profiles, summer KEmoor is approximately constant in depth while showing a clear decrease by a factor of 4–5 from 500 to 3500 m in winter. Spatial coverage is sufficient in the Northern Hemisphere to resolve broad KEmoor maxima in the western portion of each ocean basin in winter, generally collocated with FML maxima associated with storm forcing. The ratio of depth-integrated KEmoor to FML gives a replenishment time scale, which is about 10 days in midlatitudes, consistent with 1) previous estimates of the dissipation time scale of the internal wave continuum and 2) the presence of a seasonal cycle. Its increase to ≈70–80 days at lower latitudes is a possible signature of equatorward propagation of near-inertial waves. The seasonal modulation of the magnitude of KEmoor, its similarity to that in FML, and the depth decay and western intensification in winter but not in summer are consistent with a primarily wind-forced near-inertial field for latitudes poleward of ≈25°. [ABSTRACT FROM AUTHOR]

Using a set of 80 globally distributed time series of near-inertial and semidiurnal energy E and energy flux F computed from historical moorings, the group velocity ĉg ≡ FE-1 is estimated. For a single free wave, observed group speed |ĉg| should equal that expected from linear wave theory. For comparison, the latitude dependence of perceived group speed for perfectly standing waves is also derived. The latitudinal dependence of observed semidiurnal |ĉg| closely follows that expected for free waves at all latitudes, implying that 1) low-mode internal tides obey linear theory and 2) standing internal-tidal waves are rare in the deep ocean for latitudes equatorward of about 35°. At higher latitudes, standing waves cannot be easily distinguished from free waves using this method because their expected group speeds are similar. Near-inertial waves exhibit scattered |ĉg| values consistent with the passage of events generated at various latitudes, with implied frequencies ω ≈ 1.05–1.25 × f, as typically observed in frequency spectra. [ABSTRACT FROM AUTHOR]

Extending an earlier attempt to understand long-range propagation of the global internal-wave field, the energy E and horizontal energy flux F are computed for the two gravest baroclinic modes at 80 historical moorings around the globe. With bandpass filtering, the calculation is performed for the semidiurnal band (emphasizing M2 internal tides, generated by flow over sloping topography) and for the near-inertial band (emphasizing wind-generated waves near the Coriolis frequency). The time dependence of semidiurnal E and F is first examined at six locations north of the Hawaiian Ridge; E and F typically rise and fall together and can vary by over an order of magnitude at each site. This variability typically has a strong spring–neap component, in addition to longer time scales. The observed spring tides at sites northwest of the Hawaiian Ridge are coherent with barotropic forcing at the ridge, but lagged by times consistent with travel at the theoretical mode-1 group speed from the ridge. Phase computed from 14-day windows varies by approximately ±45° on monthly time scales, implying refraction by mesoscale currents and stratification. This refraction also causes the bulk of internal-tide energy flux to be undetectable by altimetry and other long-term harmonic-analysis techniques. As found previously, the mean flux in both frequency bands is O(1 kW m-1), sufficient to radiate a substantial fraction of energy far from each source. Tidal flux is generally away from regions of strong topography. Near-inertial flux is overwhelmingly equatorward, as required for waves generated at the inertial frequency on a β plane, and is winter-enhanced, consistent with storm generation. In a companion paper, the group velocity, ĉg ≡ FE-1, is examined for both frequency bands. [ABSTRACT FROM AUTHOR]

Energy flux is a fundamental quantity for understanding internal wave generation, propagation, and dissipation. In this paper, the estimation of internal wave energy fluxes 〈u′p′〉 from ocean observations that may be sparse in either time or depth are considered. Sampling must be sufficient in depth to allow for the estimation of the internal wave–induced pressure anomaly p′ using the hydrostatic balance, and sufficient in time to allow for phase averaging. Data limitations that are considered include profile time series with coarse temporal or vertical sampling, profiles missing near-surface or near-bottom information, moorings with sparse vertical sampling, and horizontal surveys with no coherent resampling in time. Methodologies, interpretation, and errors are described. For the specific case of the semidiurnal energy flux radiating from the Hawaiian ridge, errors of ∼10% are typical for estimates from six full-depth profiles spanning 15 h. [ABSTRACT FROM AUTHOR]

The energy flux from the wind to inertial mixed layer motions is computed for all oceans from 50°S to 50°N for the years 1996–99. The wind stress, τ, is computed from 6-h, 2.5°-resolution NCEP–NCAR global reanalysis surface winds. The inertial mixed layer response, u[sub I] , and the energy flux, Π = τ · u[sub I] , are computed using a slab model. The validity of the reanalysis winds and the slab model is demonstrated by direct comparison with wind and ADCP velocity records from NDBC buoys. (At latitudes > 50°, the inertial response is too fast to be resolved by the reanalysis wind 6-h output interval.) Midlatitude storms produce the greatest fluxes, resulting in broad maxima near 40° latitude during each hemisphere's winter, concentrated in the western portion of each basin. Northern Hemisphere fluxes exceed those in the Southern Hemisphere by about 50%. The global mean energy flux from 1996 to 1999 and 50°S to 50°N is (0.98 ± 0.08) × 10[sup -3] W m[sup -2] , for a total power of 0.29 TW (1 TW = 10[sup 12] W). This total is the same order of magnitude as recent estimates of the global power input to baroclinic M[sub 2] tidal motions, suggesting that wind-generated near-inertial waves may play an important role in the global energy balance. [ABSTRACT FROM AUTHOR]