On the existence of extreme waves and the Stokes conjecture with vorticity

This is a study of singular solutions of the problem of traveling gravity water waves on flows with vorticity. We show that, for a certain class of vorticity functions, a sequence of regular waves converges to an extreme wave with stagnation points at its crests. We also show that, for any vorticity function, the profile of an extreme wave must have either a corner of 120° or a horizontal tangent at any stagnation point about which it is supposed symmetric. Moreover, the profile necessarily has a corner of 120° if the vorticity is nonnegative near the free surface.

Steady periodic water waves with constant vorticity: regularity and local bifurcation

This paper studies periodic traveling gravity waves at the free surface of water in a flow of constant vorticity over a flat bed. Using conformal mappings the free-boundary problem is transformed into a quasilinear pseudodifferential equation for a periodic function of one variable. The new formulation leads to a regularity result and, by use of bifurcation theory, to the existence of waves of small amplitude even in the presence of stagnation points in the flow.

Baroclinic waves with parameterized effects of moisture interpreted using Rossby wave components

A theoretical framework is developed for the evolution of baroclinic waves with latent heat release parameterized in terms of vertical velocity. Both wave–conditional instability of the second kind (CISK) and large-scale rain approaches are included. The new quasigeostrophic framework covers evolution from general initial conditions on zonal flows with vertical shear, planetary vorticity gradient, a lower boundary, and a tropopause. The formulation is given completely in terms of potential vorticity, enabling the partition of perturbations into Rossby wave components, just as for the dry problem. Both modal and nonmodal development can be understood to a good approximation in terms of propagation and interaction between these components alone. The key change with moisture is that growing normal modes are described in terms of four counterpropagating Rossby wave (CRW) components rather than two. Moist CRWs exist above and below the maximum in latent heating, in addition to the upper- and lower-level CRWs of dry theory. Four classifications of baroclinic development are defined by quantifying the strength of interaction between the four components and identifying the dominant pairs, which range from essentially dry instability to instability in the limit of strong heating far from boundaries, with type-C cyclogenesis and diabatic Rossby waves being intermediate types. General initial conditions must also include passively advected residual PV, as in the dry problem.

Equatorial waves in opposite QBO phases

A methodology for identifying equatorial waves is used to analyze the multilevel 40-yr ECMWF Re-Analysis (ERA-40) data for two different years (1992 and 1993) to investigate the behavior of the equatorial waves under opposite phases of the quasi-biennial oscillation (QBO). A comprehensive view of 3D structures and of zonal and vertical propagation of equatorial Kelvin, westward-moving mixed Rossby–gravity (WMRG), and n = 1 Rossby (R1) waves in different QBO phases is presented. Consistent with expectation based on theory, upward-propagating Kelvin waves occur more frequently during the easterly QBO phase than during the westerly QBO phase. However, the westward-moving WMRG and R1 waves show the opposite behavior. The presence of vertically propagating equatorial waves in the stratosphere also depends on the upper tropospheric winds and tropospheric forcing. Typical propagation parameters such as the zonal wavenumber, zonal phase speed, period, vertical wavelength, and vertical group velocity are found. In general, waves in the lower stratosphere have a smaller zonal wavenumber, shorter period, faster phase speed, and shorter vertical wavelength than those in the upper troposphere. All of the waves in the lower stratosphere show an upward group velocity and downward phase speed. When the phase of the QBO is not favorable for waves to propagate, their phase speed in the lower stratosphere is larger and their period is shorter than in the favorable phase, suggesting Doppler shifting by the ambient flow and a filtering of the slow waves. Tropospheric WMRG and R1 waves in the Western Hemisphere also show upward phase speed and downward group velocity, with an indication of their forcing from middle latitudes. Although the waves observed in the lower stratosphere are dominated by “free” waves, there is evidence of some connection with previous tropical convection in the favorable year for the Kelvin waves in the warm water hemisphere and WMRG and R1 waves in the Western Hemisphere, which is suggestive of the importance of convective forcing for the existence of propagating coupled Kelvin waves and midlatitude forcing for the existence of coupled WMRG and R1 waves.

Generation of inertia-gravity waves in the rotating thermal annulus by a localised boundary layer instability

Waves with periods shorter than the inertial period exist in the atmosphere (as inertia-gravity waves) and in the oceans (as Poincaré and internal gravity waves). Such waves owe their origin to various mechanisms, but of particular interest are those arising either from local secondary instabilities or spontaneous emission due to loss of balance. These phenomena have been studied in the laboratory, both in the mechanically-forced and the thermally-forced rotating annulus. Their generation mechanisms, especially in the latter system, have not yet been fully understood, however. Here we examine short period waves in a numerical model of the rotating thermal annulus, and show how the results are consistent with those from earlier laboratory experiments. We then show how these waves are consistent with being inertia-gravity waves generated by a localised instability within the thermal boundary layer, the location of which is determined by regions of strong shear and downwelling at certain points within a large-scale baroclinic wave flow. The resulting instability launches small-scale inertia-gravity waves into the geostrophic interior of the flow. Their behaviour is captured in fully nonlinear numerical simulations in a finite-difference, 3D Boussinesq Navier-Stokes model. Such a mechanism has many similarities with those responsible for launching small- and meso-scale inertia-gravity waves in the atmosphere from fronts and local convection.

Equivalence of weak formulations of the steady water waves equations

We prove the equivalence of three weak formulations of the steady water waves equations, namely: the velocity formulation, the stream function formulation and the Dubreil-Jacotin formulation, under weak Hölder regularity assumptions on their solutions.

The interaction of flexural-gravity waves with periodic geometries

A periodic structure of finite extent is embedded within an otherwise uniform two-dimensional system consisting of finite-depth fluid covered by a thin elastic plate. An incident harmonic flexural-gravity wave is scattered by the structure. By using an approximation to the corresponding linearised boundary value problem that is based on a slowly varying structure in conjunction with a transfer matrix formulation, a method is developed that generates the whole solution from that for just one cycle of the structure, providing both computational savings and insight into the scattering process. Numerical results show that variations in the plate produce strong resonances about the ‘Bragg frequencies’ for relatively few periods. We find that certain geometrical variations in the plate generate these resonances above the Bragg value, whereas other geometries produce the resonance below the Bragg value. The familiar resonances due to periodic bed undulations tend to be damped by the plate.

The Stokes conjecture for waves with vorticity

We study stagnation points of two-dimensional steady gravity free-surface water waves with vorticity. We obtain for example that, in the case where the free surface is an injective curve, the asymptotics at any stagnation point is given either by the “Stokes corner flow” where the free surface has a corner of 120°, or the free surface ends in a horizontal cusp, or the free surface is horizontally flat at the stagnation point. The cusp case is a new feature in the case with vorticity, and it is not possible in the absence of vorticity. In a second main result we exclude horizontally flat singularities in the case that the vorticity is 0 on the free surface. Here the vorticity may have infinitely many sign changes accumulating at the free surface, which makes this case particularly difficult and explains why it has been almost untouched by research so far. Our results are based on calculations in the original variables and do not rely on structural assumptions needed in previous results such as isolated singularities, symmetry and monotonicity.

The importance of friction in mountain wave drag amplification by Scorer parameter resonance

A mechanism for amplification of mountain waves, and their associated drag, by parametric resonance is investigated using linear theory and numerical simulations. This mechanism, which is active when the Scorer parameter oscillates with height, was recently classified by previous authors as intrinsically nonlinear. Here it is shown that, if friction is included in the simplest possible form as a Rayleigh damping, and the solution to the Taylor-Goldstein equation is expanded in a power series of the amplitude of the Scorer parameter oscillation, linear theory can replicate the resonant amplification produced by numerical simulations with some accuracy. The drag is significantly altered by resonance in the vicinity of n/l_0 = 2, where l_0 is the unperturbed value of the Scorer parameter and n is the wave number of its oscillation. Depending on the phase of this oscillation, the drag may be substantially amplified or attenuated relative to its non-resonant value, displaying either single maxima or minima, or double extrema near n/l_0 = 2. Both non-hydrostatic effects and friction tend to reduce the magnitude of the drag extrema. However, in exactly inviscid conditions, the single drag maximum and minimum are suppressed. As in the atmosphere friction is often small but non-zero outside the boundary layer, modelling of the drag amplification mechanism addressed here should be quite sensitive to the type of turbulence closure employed in numerical models, or to computational dissipation in nominally inviscid simulations.

A linear model for the structure of turbulence beneath surface water waves

The structure of turbulence in the ocean surface layer is investigated using a simplified semi-analytical model based on rapid-distortion theory. In this model, which is linear with respect to the turbulence, the flow comprises a mean Eulerian shear current, the Stokes drift of an irrotational surface wave, which accounts for the irreversible effect of the waves on the turbulence, and the turbulence itself, whose time evolution is calculated. By analysing the equations of motion used in the model, which are linearised versions of the Craik–Leibovich equations containing a ‘vortex force’, it is found that a flow including mean shear and a Stokes drift is formally equivalent to a flow including mean shear and rotation. In particular, Craik and Leibovich’s condition for the linear instability of the first kind of flow is equivalent to Bradshaw’s condition for the linear instability of the second. However, the present study goes beyond linear stability analyses by considering flow disturbances of finite amplitude, which allows calculating turbulence statistics and addressing cases where the linear stability is neutral. Results from the model show that the turbulence displays a structure with a continuous variation of the anisotropy and elongation, ranging from streaky structures, for distortion by shear only, to streamwise vortices resembling Langmuir circulations, for distortion by Stokes drift only. The TKE grows faster for distortion by a shear and a Stokes drift gradient with the same sign (a situation relevant to wind waves), but the turbulence is more isotropic in that case (which is linearly unstable to Langmuir circulations).

On the momentum fluxes associated with mountain waves in directionally sheared flows

The direct impact of mountain waves on the atmospheric circulation is due to the deposition of wave momentum at critical levels, or levels where the waves break. The first process is treated analytically in this study within the framework of linear theory. The variation of the momentum flux with height is investigated for relatively large shears, extending the authors’ previous calculations of the surface gravity wave drag to the whole atmosphere. A Wentzel–Kramers–Brillouin (WKB) approximation is used to treat inviscid, steady, nonrotating, hydrostatic flow with directional shear over a circular mesoscale mountain, for generic wind profiles. This approximation must be extended to third order to obtain momentum flux expressions that are accurate to second order. Since the momentum flux only varies because of wave filtering by critical levels, the application of contour integration techniques enables it to be expressed in terms of simple 1D integrals. On the other hand, the momentum flux divergence (which corresponds to the force on the atmosphere that must be represented in gravity wave drag parameterizations) is given in closed analytical form. The momentum flux expressions are tested for idealized wind profiles, where they become a function of the Richardson number (Ri). These expressions tend, for high Ri, to results by previous authors, where wind profile effects on the surface drag were neglected and critical levels acted as perfect absorbers. The linear results are compared with linear and nonlinear numerical simulations, showing a considerable improvement upon corresponding results derived for higher Ri.

Mountain waves in two-layer sheared flows: Critical-level effects, wave reflection, and drag enhancement

Internal gravity waves generated in two-layer stratified shear flows over mountains are investigated here using linear theory and numerical simulations. The impact on the gravity wave drag of wind profiles with constant unidirectional or directional shear up to a certain height and zero shear above, with and without critical levels, is evaluated. This kind of wind profile, which is more realistic than the constant shear extending indefinitely assumed in many analytical studies, leads to important modifications in the drag behavior due to wave reflection at the shear discontinuity and wave filtering by critical levels. In inviscid, nonrotating, and hydrostatic conditions, linear theory predicts that the drag behaves asymmetrically for backward and forward shear flows. These differences primarily depend on the fraction of wavenumbers that pass through their critical level before they are reflected by the shear discontinuity. If this fraction is large, the drag variation is not too different from that predicted for an unbounded shear layer, while if it is small the differences are marked, with the drag being enhanced by a considerable factor at low Richardson numbers (Ri). The drag may be further enhanced by nonlinear processes, but its qualitative variation for relatively low Ri is essentially unchanged. However, nonlinear processes seem to interact constructively with shear, so that the drag for a noninfinite but relatively high Ri is considerably larger than the drag without any shear at all.