skip to main content


Title: Surface and Sub‐Surface Kinetic Energy Wavenumber‐Frequency Spectra in Global Ocean Models and Observations
Abstract

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.

 
more » « less
NSF-PAR ID:
10511168
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Oceans
Volume:
129
Issue:
6
ISSN:
2169-9275
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The internal wave (IW) continuum of a regional ocean model is studied in terms of the vertical spectral kinetic energy (KE) fluxes and transfers at high vertical wavenumbers. Previous work has shown that this model permits a partial representation of the IW cascade. In this work, vertical spectral KE flux is decomposed into catalyst, source, and destination vertical modes and frequency bands of nonlinear scattering, a framework that allows for the discernment of different types of nonlinear interactions involving both waves and eddies. Energy transfer within the supertidal IW continuum is found to be strongly dependent on resolution. Specifically, at a horizontal grid spacing of 1/48°, most KE in the supertidal continuum arrives there from lower-frequency modes through a single nonlinear interaction, whereas at 1/384° and with sufficient vertical resolution KE transfers within the supertidal IW continuum are comparable in size to KE transfer from lower-frequency modes. Additionally, comparisons are made with existing theoretical and observational work on energy pathways in the IW continuum. Induced diffusion (ID) is found to be associated with a weak forward frequency transfer within the supertidal IW continuum. ID is also limited to the highest vertical wavenumbers and is more sensitive to resolution relative to spectrally local interactions. At the same time, ID-like processes involving high-vertical-wavenumber near-inertial and tidal waves as well as low-vertical-wavenumber eddy fields are substantial, suggesting that the processes giving rise to a Garrett–Munk-like spectra in the present numerical simulation and perhaps the real ocean may be more varied than in idealized or wave-only frameworks.

     
    more » « less
  2. Abstract

    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.

     
    more » « less
  3. Abstract

    Recent work demonstrates that high‐resolution global models forced simultaneously by atmospheric fields and the astronomical tidal potential contain a partial internal (gravity) wave (IW) spectral continuum. Regional simulations of the MITgcm forced at the horizontal boundaries by a global run that carries a partial IW continuum spectrum are performed at the same grid spacing as the global run and at finer grid spacings in an attempt to fill out more of the IW spectral continuum. Decreasing only the horizontal grid spacing from 2 to 0.25 km greatly improves the frequency spectra and slightly improves the vertical wavenumber spectra of the horizontal velocity. Decreasing only the vertical grid spacing by a factor of 3 does not yield any significant improvements. Decreasing both horizontal and vertical grid spacings yields the greatest degree of improvement, filling the frequency spectrum out to 72 cpd. Our results suggest that improved IW spectra in regional models are possible if they are run at finer grid spacings and are forced at their lateral boundaries by remotely generated IWs. Additionally, consistency relations demonstrate that improvements in the spectra are indeed due to the existence of IWs at higher frequencies and vertical wavenumbers when remote IW forcing is included and model grid spacings decrease. By being able to simulate an IW spectral continuum to 0.25 km scales, these simulations demonstrate that one may be able to track the energy pathways of IWs from generation to dissipation and improve the understanding of processes such as IW‐driven mixing.

     
    more » « less
  4. Abstract Hundreds of full-depth temperature and salinity profiles collected by Deepglider autonomous underwater vehicles (AUVs) in the North Atlantic reveal robust signals in eddy isopycnal vertical displacement and horizontal current throughout the entire water column. In separate glider missions southeast of Bermuda, subsurface-intensified cold, fresh coherent vortices were observed with velocities exceeding 20 cm s −1 at depths greater than 1000 m. With vertical resolution on the order of 20 m or less, these full-depth glider slant profiles newly permit estimation of scaled vertical wavenumber spectra from the barotropic through the 40th baroclinic mode. Geostrophic turbulence theory predictions of spectral slopes associated with the forward enstrophy cascade and proportional to inverse wavenumber cubed generally agree with glider-derived quasi-universal spectra of potential and kinetic energy found at a variety of locations distinguished by a wide range of mean surface eddy kinetic energy. Water-column average spectral estimates merge at high vertical mode number to established descriptions of internal wave spectra. Among glider mission sites, geographic and seasonal variability implicate bottom drag as a mechanism for dissipation, but also the need for more persistent sampling of the deep ocean. Significance Statement Relative to upper-ocean measurements of temperature, salinity, and velocity, deep ocean measurements (below 2000 m) are fewer in number and more difficult to collect. Deep measurements are needed, however, to explore the nature of deep ocean circulation contributing to the global redistribution of heat and to determine how upper-ocean behavior impacts or drives deep motions. Understanding of geographic and temporal variability in vertical structures of currents and eddies enables improved description of energy pathways in the ocean driven by turbulent interactions. In this study, we use newly developed autonomous underwater vehicles, capable of diving to the seafloor and back on a near daily basis, to collect high-resolution full ocean depth measurements at various locations in the North Atlantic. These measurements reveal connections between surface and deep motions, and importantly show their time evolution. Results of analyzing these vertical structures reveal the deep ocean to regularly “feel” events in the upper ocean and permit new comparisons to deep motions in climate models. 
    more » « less
  5. The effects of horizontal resolution and wave drag damping on the semidiurnal M2 tidal energetics are studied for two realistically-forced global HYbrid Coordinate Ocean Model (HYCOM) simulations with 41 layers and horizontal resolutions of 8 km (H12) and 4 km (H25). In both simulations, the surface tidal error is minimized by tuning the strength of the linear wave drag, which is a parameterization of the surface-tide energy conversion to the unresolved baroclinic wave modes. In both simulations the M2 surface tide error with TPXO8-atlas, an altimetry constrained model, is 2.6 cm. Compared to H12, the surface tide energy conversion to the resolved vertical modes is increased by 50% in H25. This coincides with an equivalent reduction in the tuned loss of energy from the surface tide to the wave drag. For the configurations studied here, the horizontal and not the vertical resolution is the factor limiting the number of vertical modes that are resolved in most of the global ocean: modes 1–2 in H12 and modes 1–5 in H25. The wave drag also dampens the resolved internal tides. The 40% reduction in wave-drag strength does not result in a proportional increase in the mode-1 energy density in H25. In the higher-resolution simulations, topographic mode-scattering and wave–wave interactions are better resolved. This allows for an energy flux out of mode 1 to the higher modes, mitigating the need for an internal tide damping term. The HYCOM simulations are validated with analytical conversion models and altimetry-inferred sea-surface height, fluxes, and surface tide dissipation. H25 agrees best with these data sets to within 10%. To facilitate the comparison of stationary tide signals extracted from time series with different durations, we successfully apply a spatially-varying correction factor. 
    more » « less