Search

A combination of in situ and remotely sensed observations are used to constrain the imprint of submesoscale turbulence in the sea surface height (SSH) field. The distribution of SSH variance across frequencies and wavenumbers is estimated by comparing an empirical model spectrum to two sets of observations. First, submesoscale SSH variance is constrained using a pair of GPS buoys spaced at about 10 km. From these data, one can estimate frequency spectra not only of SSH variance but also of the variance in the SSH difference between the buoys. The ratio between these two spectral estimates is sensitive to how much SSH variance is present in the submesoscale range and thus constrains the spectral roll-off of SSH variance in wavenumber space. Second, a combination of moored current meters and nadir altimetry is used to obtain an independent constraint. This constraint is enabled by geostrophy and the nonseparability of the wavenumber–frequency spectrum of SSH variance revealed by the GPS data. The frequency spectra of kinetic energy and SSH variance follow different power laws, and the difference constrains the spectral content in wavenumber space, allowing for a constraint without the need to actually resolve the submesoscales in space. In all four locations studied, spanning the midlatitude and subtropical ocean, these constraints indicate that the wavenumber spectral roll-off of submesoscale SSH variance is between about k−4 and k−5, where k is the wavenumber. These estimates are consistent with previous observations, model results, and theoretical predictions. They provide for a strong prior for the interpretation of upcoming high-resolution satellite data.

To close the overturning circulation, dense bottom water must upwell via turbulent mixing. Recent studies have identified thin bottom boundary layers (BLs) as locations of intense upwelling, yet it remains unclear how they interact with and shape the large-scale circulation of the abyssal ocean. The current understanding of this BL–interior coupling is shaped by 1D theory, suggesting that variations in locally produced BL transport generate exchange with the interior and thus a global circulation. Until now, however, this picture has been based on a 1D theory that fails to capture the local evolution in even highly idealized 2D geometries. The present work applies BL theory to revised 1D dynamics, which more naturally generalizes to two and three dimensions. The BL is assumed to be in quasi-equilibrium between the upwelling of dense water and the convergence of downward buoyancy fluxes. The BL transport, for which explicit formulas are presented, exerts an influence on the interior by modifying the bottom boundary condition. In 1D, this BL transport is independent of the interior evolution, but in 2D the BL and interior are fully coupled. Once interior variables and the bottom slope are allowed to vary in the horizontal, the resulting convergences and divergences in the BL transport exchange mass with the interior. This framework allows for the analysis of previously inaccessible problems such as the BL–interior coupling in the presence of an exponential interior stratification, laying the foundation for developing a full theory for the abyssal circulation.

We investigate the instabilities and transition mechanisms of Boussinesq stratified boundary layers on sloping boundaries when subjected to oscillatory body forcing parallel to the slope. We examine idealized forms of boundary layers on hydraulically smooth abyssal slopes in tranquil mid- to low-latitude regions, where low-wavenumber internal tides gently heave isopycnals up and down adiabatic slopes in the absence of mean flows, high-wavenumber internal tides, shelf breaks, resonant tide–bathymetry interactions (critical slopes) and other phenomena associated with turbulence ‘hot spots’. In non-rotating low-Reynolds-number flow, increased stratification on the downslope phase has a relaminarizing effect, while on the upslope phase we find transition-to-turbulence pathways arise from shear production triggered by gravitational instabilities. When rotation is significant (low slope Burger numbers) we find that boundary layer turbulence is sustained throughout the oscillation period, resembling stratified Stokes–Ekman layer turbulence. Simulation results suggest that oscillating boundary layers on smooth slopes at low Reynolds number ($\textit {Re}\leqslant 840$), unity Prandtl number and slope Burger numbers greater than unity do not cause significant irreversible turbulent buoyancy flux (mixing), and that flat-bottom dissipation rate models derived from the tide amplitude are accurate within an order of magnitude.

ABSTRACT: The abyssal overturning circulation is thought to be primarily driven by small-scale turbulent mixing. Diagnosed watermass transformations are dominated by rough topography “hotspots”, where the bottom-enhancement of mixing causes the diffusive buoyancy flux to diverge, driving widespread downwelling in the interior—only to be overwhelmed by an even stronger up-welling in a thin Bottom Boundary Layer (BBL). These watermass transformations are significantly underestimated by one-dimensional (1D) sloping boundary layer solutions, suggesting the importance of three-dimensional physics. Here, we use a hierarchy of models to generalize this 1D boundary layer approach to three-dimensional eddying flows over realistically rough topography. When applied to the Mid-Atlantic Ridge in the Brazil Basin, the idealized simulation results are roughly consistent with available observations. Integral buoyancy budgets isolate the physical processes that contribute to realistically strong BBL upwelling. The downwards diffusion of buoyancy is primarily balanced by upwelling along the sloping canyon sidewalls and the surrounding abyssal hills. These flows are strengthened by the restratifying effects of submesoscale baroclinic eddies and by the blocking of along-ridge thermal wind within the canyon. Major topographic sills block along-thalweg flows from restratifying the canyon trough, resulting in the continual erosion of the trough’s stratification. We propose simple modifications to the 1D boundary layer model which approximate each of these three-dimensional effects. These results provide local dynamical insights into mixing-driven abyssal overturning, but a complete theory will also require the non-local coupling to the basin-scale circulation.

27 pags., 14 figs., The cross-correlation of a diffuse or random wavefield at two points has been demonstrated to recover an empirical estimate of the Green's function under a wide variety of source conditions. Over the past two decades, the practical development of this principle, termed ambient noise interferometry, has revolutionized the fields of seismology and acoustics. Yet, because of the spatial sparsity of conventional water column and seafloor instrumentation, such array-based processing approaches have not been widely utilized in oceanography. Ocean-bottom distributed acoustic sensing (OBDAS) repurposes pre-existing optical fibers laid in seafloor cables as dense arrays of broadband strain sensors, which observe both seismic waves and ocean waves. The thousands of sensors in an OBDAS array make ambient noise interferometry of ocean waves straightforward for the first time. Here, we demonstrate the application of ambient noise interferometry to surface gravity waves observed on an OBDAS array near the Strait of Gibraltar. We focus particularly on a 3-km segment of the array on the continental shelf, containing 300 channels at 10-m spacing. By cross-correlating the raw strain records, we compute empirical ocean surface gravity wave Green's functions for each pair of stations. We first apply beamforming to measure the time-averaged dispersion relation along the cable. Then, we exploit the non-reciprocity of waves propagating in a flow to recover the depth-averaged current velocity as a function of time using a waveform stretching method. The result is a spatially continuous matrix of current velocity measurements with resolution, We thank Red Eléctrica de España for access to their cable and Aragon Photonics for providing the interrogator unit. E. F. Williams was supported by an NSF Graduate Research Fellow-ship. Z. Zhan acknowledges support from the Moore Foundation and NSF under CAREER Award 1848166. M. R. Fernández-Ruiz acknowledges support from the Spanish MICINN under contract IJC2018-035684-I. M. González-Herráez acknowledges support from projects RTI2018-097957-B-C31 and RTI2018-097957-B-C33 funded by the Spanish Ministry of Science and Innovation, project SINFOTON2-CM: P2018/NMT-4326 funded by Comunidad de Madrid and Fondo Europeo de Desarrollo Regional (FEDER) Program, and project PLEC2021-007875 funded by the Spanish Ministry of Science and Innovation MCIN/AEI/10.13039/501100011033 and by the European Union NextGenera-tionEU/PRTR program

The near-bottom mixing that allows abyssal waters to upwell tilts isopycnals and spins up flow over the flanks of mid-ocean ridges. Meso- and large-scale currents along sloping topography are subjected to a delicate balance of Ekman arrest and spin down. These two seemingly disparate oceanographic phenomena share a common theory, which is based on a one-dimensional model of rotating, stratified flow over a sloping, insulated boundary. This commonly used model, however, lacks rapid adjustment of interior flows, limiting its ability to capture the full physics of spin up and spin down of along-slope flow. Motivated by two-dimensional dynamics, the present work extends the one-dimensional model by constraining the vertically integrated cross-slope transport and allowing for a barotropic cross-slope pressure gradient. This produces a closed secondary circulation by forcing Ekman transport in the bottom boundary layer to return in the interior. The extended model can thus capture Ekman spin up and spin down physics: the interior return flow is turned by the Coriolis acceleration, leading to rapid rather than slow diffusive adjustment of the along-slope flow. This transport-constrained one-dimensional model accurately describes twodimensional mixing-generated spin up over an idealized ridge and provides a unified framework for understanding the relative importance of Ekman arrest and spin down of flow along a slope.

Mesoscale eddies contain the bulk of the ocean’s kinetic energy (KE), but fundamental questions remain on the cross-scale KE transfers linking eddy generation and dissipation. The role of submesoscale flows represents the key point of discussion, with contrasting views of submesoscales as either a source or a sink of mesoscale KE. Here, the first observational assessment of the annual cycle of the KE transfer between mesoscale and submesoscale motions is performed in the upper layers of a typical open-ocean region. Although these diagnostics have marginal statistical significance and should be regarded cautiously, they are physically plausible and can provide a valuable benchmark for model evaluation. The cross-scale KE transfer exhibits two distinct stages, whereby submesoscales energize mesoscales in winter and drain mesoscales in spring. Despite this seasonal reversal, an inverse KE cascade operates throughout the year across much of the mesoscale range. Our results are not incompatible with recent modeling investigations that place the headwaters of the inverse KE cascade at the submesoscale, and that rationalize the seasonality of mesoscale KE as an inverse cascade-mediated response to the generation of submesoscales in winter. However, our findings may challenge those investigations by suggesting that, in spring, a downscale KE transfer could dampen the inverse KE cascade. An exploratory appraisal of the dynamics governing mesoscale–submesoscale KE exchanges suggests that the upscale KE transfer in winter is underpinned by mixed layer baroclinic instabilities, and that the downscale KE transfer in spring is associated with frontogenesis. Current submesoscale-permitting ocean models may substantially understate this downscale KE transfer, due to the models’ muted representation of frontogenesis.

While the distribution of kinetic energy across spatial scales in the submesoscale range (1–100 km) has been estimated from observations, the associated time scales are largely unconstrained. These time scales can provide important insight into the dynamics of submesoscale turbulence because they help quantify to what degree the flow is subinertial and thus constrained by Earth’s rotation. Here a mooring array is used to estimate these time scales in the northeast Atlantic. Frequency-resolved structure functions indicate that energetic wintertime submesoscale turbulence at spatial scales around 10 km evolves on time scales of about 1 day. While these time scales are comparable to the inertial period, the observed flow also displays characteristics of subinertial flow that is geostrophically balanced to leading order. An approximate Helmholtz decomposition shows the order 10-km flow to be dominated by its rotational component, and the root-mean-square Rossby number at these scales is estimated to be 0.3. This rotational dominance and Rossby numbers below one persist down to 2.6 km, the smallest spatial scale accessible by the mooring array, despite substantially superinertial Eulerian evolution. This indicates that the Lagrangian evolution of submesoscale turbulence is slower than the Eulerian time scale estimated from the moorings. The observations therefore suggest that, on average, submesoscale turbulence largely follows subinertial dynamics in the 1–100-km range, even if Doppler shifting produces superinertial Eulerian evolution. Ageostrophic motions become increasingly important for the evolution of submesoscale turbulence as the scale is reduced—the root-mean-square Rossby number reaches 0.5 at a spatial scale of 2.6 km.

Hearing the heat Most of the excess heat that causes global warming is absorbed by the oceans. Quantifying that heat increase is challenging because it requires many different temperature measurements over both the vertical and horizontal extent of the oceans. Wu et al. report success in this effort through the use of a different method: They inferred temperature changes from sound waves generated by repeating earthquakes (see the Perspective by Wunsch). The travel time of these earthquakes from source to receiver reflects changes in the average water temperature that they encounter. This technique should substantially enhance our ability to monitor ocean warming. Science , this issue p. 1510 ; see also p. 1433

Weakly stratified layers over sloping topography can support a submesoscale baroclinic instability mode, a bottom boundary layer counterpart to surface mixed layer instabilities. The instability results from the release of available potential energy, which can be generated because of the observed bottom intensification of turbulent mixing in the deep ocean, or the Ekman adjustment of a current on a slope. Linear stability analysis suggests that the growth rates of bottom boundary layer baroclinic instabilities can be comparable to those of the surface mixed layer mode and are relatively insensitive to topographic slope angle, implying the instability is robust and potentially active in many areas of the global oceans. The solutions of two separate one-dimensional theories of the bottom boundary layer are both demonstrated to be linearly unstable to baroclinic instability, and results from an example nonlinear simulation are shown. Implications of these findings for understanding bottom boundary layer dynamics and processes are discussed.

The emerging view of the abyssal circulation is that it is associated with bottom-enhanced mixing, which results in downwelling in the stratified ocean interior and upwelling in a bottom boundary layer along the insulating and sloping seafloor. In the limit of slowly varying vertical stratification and topography, however, boundary layer theory predicts that these upslope and downslope flows largely compensate, such that net water mass transformations along the slope are vanishingly small. Using a planetary geostrophic circulation model that resolves both the boundary layer dynamics and the large-scale overturning in an idealized basin with bottom-enhanced mixing along a midocean ridge, we show that vertical variations in stratification become sufficiently large at equilibrium to reduce the degree of compensation along the midocean ridge flanks. The resulting large net transformations are similar to estimates for the abyssal ocean and span the vertical extent of the ridge. These results suggest that boundary flows generated by mixing play a crucial role in setting the global ocean stratification and overturning circulation, requiring a revision of abyssal ocean theories.

Submesoscale dynamics are typically intensified at boundaries and assumed to weaken below the mixed layer in the open ocean. Here, we assess both the seasonality and the vertical distribution of submesoscale motions in an open-ocean region of the northeast Atlantic. Second-order structure functions, or variance in properties separated by distance, are calculated from submesoscale-resolving ocean glider and mooring observations, as well as a 1/48° numerical ocean model. This dataset combines a temporal coverage that extends through a full seasonal cycle, a horizontal resolution that captures spatial scales as small as 1 km, and vertical sampling that provides near-continuous coverage over the upper 1000 m. While kinetic and potential energies undergo a seasonal cycle, being largest during the winter, structure function slopes, influenced by dynamical characteristics, do not exhibit a strong seasonality. Furthermore, structure function slopes show weak vertical variations; there is not a strong change in properties across the base of the mixed layer. Additionally, we compare the observations to output from a high-resolution numerical model. The model does not represent variability associated with superinertial motions and does not capture an observed reduction in submesoscale kinetic energy that occurs throughout the water column in spring. Overall, these results suggest that the transfer of mixed layer submesoscale variability down to depths below the traditionally defined mixed layer is important throughout the weakly stratified subpolar mode waters.

To close the abyssal overturning circulation, dense bottom water has to become lighter by mixing with lighter water above. This diapycnal mixing is strongly enhanced over rough topography in abyssal mixing layers, which span the bottom few hundred meters of the water column. In particular, mixing rates are enhanced over mid-ocean ridge systems, which extend for thousands of kilometers in the global ocean and are thought to be key contributors to the required abyssal water mass transformation. To examine how stratification and thus diabatic transformation is maintained in such abyssal mixing layers, this study explores the circulation driven by bottom-intensified mixing over mid-ocean ridge flanks and within ridge-flank canyons. Idealized numerical experiments show that stratification over the ridge flanks is maintained by submesoscale baroclinic eddies and that stratification within ridge-flank canyons is maintained by mixing-driven mean flows. These restratification processes affect how strong a diabatic buoyancy flux into the abyss can be maintained, and they are essential for maintaining the dipole in water mass transformation that has emerged as the hallmark of a diabatic circulation driven by bottom-intensified mixing.

In anticipation of the Surface Water and Ocean Topography (SWOT) wide-swath altimetry mission, this study reviews expectations for sea surface height (SSH) variance spectra at wavelengths of 10–100 km. Kinetic energy spectra from in situ observations and numerical simulations indicate that SSH variance spectra associated with balanced flow drop off steeply with wavenumber, with at least the negative fourth power of the wavenumber. Such a steep drop-off implies that even drastic reductions in altimetry noise yield only a modest improvement in the resolution of balanced flow. This general expectation is made concrete by extrapolating SSH variance spectra from existing altimetry to submesoscales, the results of which suggest that in the extratropics (poleward of 20° latitude) SWOT will improve the resolution from currently about 100 km to a median of 51 or 74 km, depending on whether or not submesoscale balanced flows are energetic. Internal waves, in contrast to balanced flow, give rise to SSH variance spectra that drop off relatively gently with wavenumber, so SSH variance should become strongly dominated by internal waves in the submesoscale range. In situ observations of the internal-wave field suggest that the internal-wave signal accessible by SWOT will be largely dominated by internal tides. The internal-wave continuum is estimated to have a spectral level close to but somewhat lower than SWOT’s expected noise level.

For small-scale turbulence to achieve water mass transformation and thus affect the large-scale overturning circulation, it must occur in stratified water. Observations show that abyssal turbulence is strongly enhanced in the bottom few hundred meters in regions with rough topography, and it is thought that these abyssal mixing layers are crucial for closing and shaping the overturning circulation. If it were left unopposed, however, bottom-intensified turbulence would mix away the observed mixing-layer stratification over the course of a few years. It is proposed here that the homogenizing tendency of mixing may be balanced by baroclinic restratification. It is shown that bottom-intensified mixing, if it occurs on a large-scale topographic slope such as a midocean ridge flank, not only erodes stratification but also tilts isopycnals in the bottom few hundred meters. This tilting of isopycnals generates a reservoir of potential energy that can be tapped into by submesoscale baroclinic eddies. The eddies slide dense water under light water and thus restratify the mixing layer, similar to what happens in the surface mixed layer. This restratification is shown to be effective enough to balance the homogenizing tendency of mixing and to maintain the observed mixing-layer stratification. This suggests that submesoscale baroclinic eddies may play a crucial role in providing the stratification mixing can act on, thus allowing sustained water mass transformation. Through their restratification of abyssal mixing layers, submesoscale eddies may therefore directly affect the strength and structure of the abyssal overturning circulation.

This paper revisits how the restratifying buoyancy flux generated by baroclinic mixed layer instabilities depends on environmental conditions. The frontal spindown is shown to produce buoyancy fluxes that increase significantly beyond the previously proposed and widely used scaling (f is the Coriolis parameter, Λ is the geostrophic shear, and H is the mixed layer depth), irrespective of whether the initial front is broad or narrow. This increase occurs after the initial phase of the nonlinear evolution, when the baroclinic eddies grow in size and develop velocities significantly in excess of the scaling assumption V ~ ΛH. Implications for parameterizing the restratification caused by baroclinic mixed layer instabilities in coarse-resolution models are discussed.

Baroclinic mixed-layer instabilities have recently been recognized as an important source of submesoscale energy in deep winter mixed layers. While the focus has so far been on the balanced dynamics of these instabilities, they occur in and depend on an environment shaped by atmospherically forced small-scale turbulence. In this study, idealized numerical simulations are presented that allow the development of both baroclinic instability and convective small-scale turbulence, with simple control over the relative strength. If the convection is only weakly forced, baroclinic instability restratifies the layer and shuts off convection, as expected. With increased forcing, however, it is found that baroclinic instabilities are remarkably resilient to the presence of convection. Even if the instability is too weak to restratify the layer and shut off convection, the instability still grows in the convecting environment and generates baroclinic eddies and fronts. This suggests that despite the vigorous atmospherically forced small-scale turbulence in winter mixed layers, baroclinic instabilities can persistently grow, generate balanced submesoscale turbulence, and modify the bulk properties of the upper ocean.

Spectral analysis is applied to infer the dynamics of mesoscale winds from aircraft observations in the upper troposphere and lower stratosphere. Two datasets are analyzed: one collected aboard commercial aircraft and one collected using a dedicated research aircraft. A recently developed wave–vortex decomposition is used to test the observations’ consistency with linear inertia–gravity wave dynamics. The decomposition method is shown to be robust in the vicinity of the tropopause if flight tracks vary sufficiently in altitude. For the lower stratosphere, the decompositions of both datasets confirm a recent result that mesoscale winds are consistent with the polarization and dispersion relations of inertia–gravity waves. For the upper troposphere, however, the two datasets disagree: only the research aircraft data indicate consistency with linear wave dynamics at mesoscales. The source of the inconsistency is a difference in mesoscale variance of the measured along-track wind component. To further test the observed flow’s consistency with linear wave dynamics, the ratio between tropospheric and stratospheric mesoscale energy levels is compared to a simple model of upward-propagating waves that are partially reflected at the tropopause. For both datasets, the observed energy ratio is roughly consistent with the simple wave model, but wave frequencies diagnosed from the data draw into question the applicability of the monochromatic theory at wavelengths smaller than 10 km.

Upper-ocean turbulence at scales smaller than the mesoscale is believed to exchange surface and thermocline waters, which plays an important role in both physical and biogeochemical budgets. But what energizes this submesoscale turbulence remains a topic of debate. Two mechanisms have been proposed: mesoscale-driven surface frontogenesis and baroclinic mixed-layer instabilities. The goal here is to understand the differences between the dynamics of these two mechanisms, using a simple quasi-geostrophic model. The essence of mesoscale-driven surface frontogenesis is captured by the well-known surface quasi-geostrophic model, which describes the sharpening of surface buoyancy gradients and the subsequent breakup in secondary roll-up instabilities. We formulate a similarly archetypical Eady-like model of submesoscale turbulence induced by mixed-layer instabilities. The model captures the scale and structure of this baroclinic instability in the mixed layer. A wide range of scales are energized through a turbulent inverse cascade of kinetic energy that is fuelled by the submesoscale mixed-layer instability. Major differences to mesoscale-driven surface frontogenesis are that mixed-layer instabilities energize the entire depth of the mixed layer and produce larger vertical velocities. The distribution of energy across scales and in the vertical produced by our simple model of mixed-layer instabilities compares favourably to observations of energetic wintertime submesoscale flows, suggesting that it captures the leading-order balanced dynamics of these flows. The dynamics described here in an oceanographic context have potential applications to other geophysical fluids with layers of different stratifications.

Lateral stirring is a basic oceanographic phenomenon affecting the distribution of physical, chemical, and biological fields. Eddy stirring at scales on the order of 100 km (the mesoscale) is fairly well understood and explicitly represented in modern eddy-resolving numerical models of global ocean circulation. The same cannot be said for smaller-scale stirring processes. Here, the authors describe a major oceanographic field experiment aimed at observing and understanding the processes responsible for stirring at scales of 0.1–10 km. Stirring processes of varying intensity were studied in the Sargasso Sea eddy field approximately 250 km southeast of Cape Hatteras. Lateral variability of water-mass properties, the distribution of microscale turbulence, and the evolution of several patches of inert dye were studied with an array of shipboard, autonomous, and airborne instruments. Observations were made at two sites, characterized by weak and moderate background mesoscale straining, to contrast different regimes of lateral stirring. Analyses to date suggest that, in both cases, the lateral dispersion of natural and deliberately released tracers was O(1) m2 s–1 as found elsewhere, which is faster than might be expected from traditional shear dispersion by persistent mesoscale flow and linear internal waves. These findings point to the possible importance of kilometer-scale stirring by submesoscale eddies and nonlinear internal-wave processes or the need to modify the traditional shear-dispersion paradigm to include higher-order effects. A unique aspect of the Scalable Lateral Mixing and Coherent Turbulence (LatMix) field experiment is the combination of direct measurements of dye dispersion with the concurrent multiscale hydrographic and turbulence observations, enabling evaluation of the underlying mechanisms responsible for the observed dispersion at a new level.

Midlatitude fluctuations of the atmospheric winds on scales of thousands of kilometers, the most energetic of such fluctuations, are strongly constrained by the Earth’s rotation and the atmosphere’s stratification. As a result of these constraints, the flow is quasi-2D and energy is trapped at large scales—nonlinear turbulent interactions transfer energy to larger scales, but not to smaller scales. Aircraft observations of wind and temperature near the tropopause indicate that fluctuations at horizontal scales smaller than about 500 km are more energetic than expected from these quasi-2D dynamics. We present an analysis of the observations that indicates that these smaller-scale motions are due to approximately linear inertia–gravity waves, contrary to recent claims that these scales are strongly turbulent. Specifically, the aircraft velocity and temperature measurements are separated into two components: one due to the quasi-2D dynamics and one due to linear inertia–gravity waves. Quasi-2D dynamics dominate at scales larger than 500 km; inertia–gravity waves dominate at scales smaller than 500 km.

Submesoscale (1–200 km) wavenumber spectra of kinetic and potential energy and tracer variance are obtained from in situ observations in the Gulf Stream region and in the eastern subtropical North Pacific. In the Gulf Stream region, steep kinetic energy spectra at scales between 200 and 20 km are consistent with predictions of interior quasigeostrophic–turbulence theory, both in the mixed layer and in the thermocline. At scales below 20 km, the spectra flatten out, consistent with a growing contribution of internal-wave energy at small scales. In the subtropical North Pacific, the energy spectra are flatter and inconsistent with predictions of interior quasigeostrophic–turbulence theory. The observed spectra and their dependence on depth are also inconsistent with predictions of surface quasigeostrophic–turbulence theory for the observed ocean stratification. It appears that unbalanced motions, most likely internal tides at large scales and the internal-wave continuum at small scales, dominate the energy spectrum throughout the submesoscale range. Spectra of temperature variance along density surfaces, which are not affected by internal tides, are also inconsistent with predictions of geostrophic-turbulence theories. Reasons for this inconsistency could be the injection of energy in the submesoscale range by small-scale baroclinic instabilities or modifications of the spectra by coupling between surface and interior dynamics or by ageostrophic frontal effects.