We examine the flows induced by horizontally modulated, vertically confined (or guided), internal wavepackets in a stratified, Boussinesq fluid. The wavepacket induces both an Eulerian flow and a Stokes drift, which together determine the Lagrangian transport of passive tracers. We derive equations describing the wave-induced flows in arbitrary stable stratification and consider four special cases: a two-layer fluid, symmetric and asymmetric piecewise constant (‘top-hat’) stratification and, more representative of the ocean, exponential stratification. In a two-layer fluid, the Stokes drift is positive everywhere with the peak value at the interface, whereas the Eulerian flow is negative and uniform with depth for long groups. Combined, the net depth-integrated Lagrangian transport is zero. If one layer is shallower than the other, the wave-averaged interface displaces into that layer making the Eulerian flow in that layer more negative and the Eulerian flow in the opposite layer more positive so that the depth-integrated Eulerian transports are offset by the same amount in each layer. By contrast, in continuous stratification the depth-integrated transport due to the Stokes drift and Eulerian flow are each zero, but the Eulerian flow is singular if the horizontal phase speed of the induced flow equals the group velocity of the wavepacket, giving rise to a single resonance in uniform stratification (McIntyre, J. Fluid Mech., vol. 60, 1973, pp. 801–811). In top-hat stratification, this single resonance disappears, being replaced by multiple resonances occurring when the horizontal group velocity of the wavepacket matches the horizontal phase speed of higher-order modes. Furthermore, if the stratification is not vertically symmetric, then the Eulerian induced flow varies as the inverse squared horizontal wavenumber for shallow waves, the same as for the asymmetric two-layer case. This ‘infrared catastrophe’ also occurs in the case of exponential stratification suggesting significant backward near-surface transport by the Eulerian induced flow for modulated oceanic internal modes. Numerical simulations are performed confirming these theoretical predictions.

Numerical and observational evidence indicates that, in regions where mixed-layer instability is active, the surface geostrophic velocity is largely induced by surface buoyancy anomalies. Yet, in these regions, the observed surface kinetic energy spectrum is steeper than predicted by uniformly stratified surface quasigeostrophic theory. By generalizing surface quasigeostrophic theory to account for variable stratification, we show that surface buoyancy anomalies can generate a variety of dynamical regimes depending on the stratification’s vertical structure. Buoyancy anomalies generate longer range velocity fields over decreasing stratification and shorter range velocity fields over increasing stratification. As a result, the surface kinetic energy spectrum is steeper over decreasing stratification than over increasing stratification. An exception occurs if the near surface stratification is much larger than the deep ocean stratification. In this case, we find an extremely local turbulent regime with surface buoyancy homogenization and a steep surface kinetic energy spectrum, similar to equivalent barotropic turbulence. By applying the variable stratification theory to the wintertime North Atlantic, and assuming that mixed-layer instability acts as a narrowband small-scale surface buoyancy forcing, we obtain a predicted surface kinetic energy spectrum between k−4/3 and k−7/3, which is consistent with the observed wintertime k−2 spectrum. We conclude by suggesting a method of measuring the buoyancy frequency’s vertical structure using satellite observations.

We examine Lagrangian transport by a nonlinearly evolving vertical mode-1 internal tide in non-uniform stratification. In a companion paper (Sutherland & Dhaliwal, J. Fluid Mech., 2022, in press) it was shown that a parent internal tide can excite successive superharmonics that superimpose to form a solitary wave train. Despite this transformation, here we show that the collective forcing by the parent wave and superharmonics is effectively steady in time. Thus we derive relatively simple formulae for the Stokes drift and induced Eulerian flow associated with the waves under the assumption that the parent waves and superharmonics are long compared with the fluid depth. In all cases, the Stokes drift exhibits a mixed mode-1 and mode-2 vertical structure with the flow being in the waveward direction at the surface. If the background rotation is non-negligible, the vertical structure of the induced Eulerian flow is equal and opposite to that of the Stokes drift. This flow periodically increases and decreases at the inertial frequency with maximum magnitude twice that of the Stokes drift. When superimposed with the Stokes drift, the Lagrangian flow at the surface periodically changes from positive to negative over one inertial period. If the background rotation is zero, the induced Eulerian flow evolves non-negligibly in time and space for horizontally modulated waves: the depth below the surface of the positive Lagrangian flow becomes shallower ahead of the peak of the amplitude envelope and becomes deeper in the lee of the peak. These predictions are well-captured by fully nonlinear numerical simulations.

Three-dimensional geophysical fluids support both internal and boundary-trapped waves. To obtain the normal modes in such fluids, we must solve a differential eigenvalue problem for the vertical structure (for simplicity, we only consider horizontally periodic domains). If the boundaries are dynamically inert (e.g., rigid boundaries in the Boussinesq internal wave problem and flat boundaries in the quasigeostrophic Rossby wave problem), the resulting eigenvalue problem typically has a Sturm–Liouville form and the properties of such problems are well-known. However, when restoring forces are also present at the boundaries, then the equations of motion contain a time-derivative in the boundary conditions, and this leads to an eigenvalue problem where the eigenvalue correspondingly appears in the boundary conditions. In certain cases, the eigenvalue problem can be formulated as an eigenvalue problem in the Hilbert space L2⊕C and this theory is well-developed. Less explored is the case when the eigenvalue problem takes place in a Pontryagin space, as in the Rossby wave problem over sloping topography. This article develops the theory of such problems and explores the properties of wave problems with dynamically active boundaries. The theory allows us to solve the initial value problem for quasigeostrophic Rossby waves in a region with sloping bottom (we also apply the theory to two Boussinesq problems with a free-surface). For a step-function perturbation at a dynamically active boundary, we find that the resulting time-evolution consists of waves present in proportion to their projection onto the dynamically active boundary.

The discrete baroclinic modes of quasigeostrophic theory are incomplete and the incompleteness manifests as a loss of information in the projection process. The incompleteness of the baroclinic modes is related to the presence of two previously unnoticed stationary step-wave solutions of the Rossby wave problem with flat boundaries. These step-waves are the limit of surface quasigeostrophic waves as boundary buoyancy gradients vanish. A complete normal mode basis for quasigeostrophic theory is obtained by considering the traditional Rossby wave problem with prescribed buoyancy gradients at the lower and upper boundaries. The presence of these boundary buoyancy gradients activates the previously inert boundary degrees of freedom. These Rossby waves have several novel properties such as the presence of multiple modes with no internal zeros, a finite number of modes with negative norms, and their vertical structures form a basis capable of representing any quasigeostrophic state with a differentiable series expansion. These properties are a consequence of the Pontryagin space setting of the Rossby wave problem in the presence of boundary buoyancy gradients (as opposed to the usual Hilbert space setting). We also examine the quasigeostrophic vertical velocity modes and derive a complete basis for such modes as well. A natural application of these modes is the development of a weakly non-linear wave-interaction theory of geostrophic turbulence that takes topography into account.