Damping of inertial motions by parametric subharmonic instability in baroclinic currents

2014 ◽  
Vol 743 ◽  
pp. 280-294 ◽  
Author(s):  
Leif N. Thomas ◽  
John R. Taylor

AbstractA new damping mechanism for vertically-sheared inertial motions is described involving an inertia–gravity wave that oscillates at half the inertial frequency, $f$, and that grows at the expense of inertial shear. This parametric subharmonic instability forms in baroclinic, geostrophic currents where thermal wind shear, by reducing the potential vorticity of the fluid, allows inertia–gravity waves with frequencies less than $f$. A stability analysis and numerical simulations are used to study the instability criterion, energetics, and finite-amplitude behaviour of the instability. For a flow with uniform shear and stratification, parametric subharmonic instability develops when the Richardson number of the geostrophic current nears $Ri_{PSI}=4/3+\gamma \cos \phi $, where $\gamma $ is the ratio of the inertial to thermal wind shear magnitude and $\phi $ is the angle between the inertial and thermal wind shears at the initial time. Inertial shear enters the instability criterion because it can also modify the potential vorticity and hence the minimum frequency of inertia–gravity waves. When this criterion is met, inertia–gravity waves with a frequency $f/2$ and with flow parallel to isopycnals amplify, extracting kinetic energy from the inertial shear through shear production. The solutions of the numerical simulations are consistent with these predictions and additionally show that finite-amplitude parametric subharmonic instability both damps inertial shear and is itself damped by secondary shear instabilities. In this way, parametric subharmonic instability opens a pathway to turbulence where kinetic energy in inertial shear is transferred to small scales and dissipated.

2014 ◽  
Vol 755 ◽  
pp. 397-428 ◽  
Author(s):  
Peng Wang ◽  
James C. McWilliams ◽  
Claire Ménesguen

AbstractThe linear instability of several rotating, stably stratified, interior vertical shear flows $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}\overline{U}(z)$ is calculated in Boussinesq equations. Two types of baroclinic, ageostrophic instability, AI1 and AI2, are found in odd-symmetric $\overline{U}(z)$ for intermediate Rossby number ($\mathit{Ro}$). AI1 has zero frequency; it appears in a continuous transformation of the unstable mode properties between classic baroclinic instability (BCI) and centrifugal instability (CI). It begins to occur at intermediate $\mathit{Ro}$ values and horizontal wavenumbers ($k,l$) that are far from $l= 0$ or $k = 0$, where the growth rate of BCI or CI is the strongest. AI1 grows by drawing kinetic energy from the mean flow, and the perturbation converts kinetic energy to potential energy. The instability AI2 has inertia critical layers (ICL); hence it is associated with inertia-gravity waves. For an unstable AI2 mode, the coupling is either between an interior balanced shear wave and an inertia-gravity wave (BG), or between two inertia-gravity waves (GG). The main energy source for an unstable BG mode is the mean kinetic energy, while the main energy source for an unstable GG mode is the mean available potential energy. AI1 and BG type AI2 occur in the neighbourhood of $A-S= 0$ (a sign change in the difference between absolute vertical vorticity and horizontal strain rate in isentropic coordinates; see McWilliams et al., Phys. Fluids, vol. 10, 1998, pp. 3178–3184), while GG type AI2 arises beyond this condition. Both AI1 and AI2 are unbalanced instabilities; they serve as an initiation of a possible local route for the loss of balance in 3D interior flows, leading to an efficient energy transfer to small scales.


2001 ◽  
Vol 443 ◽  
pp. 351-376 ◽  
Author(s):  
G. M. REZNIK ◽  
R. GRIMSHAW

We consider the non-stationary dynamics of an intense localized vortex on a β-plane using a shallow-water model. An asymptotic theory for a vortex with piecewise-continuous potential vorticity is developed assuming the Rossby number to be small and the free surface elevation to be small but finite. Analogously to the well-known quasi-geostrophic model, the vortex translation is produced by a secondary dipole circulation (β-gyres) developed in the vortex vicinity and consisting of two parts. The first part (geostrophic β-gyres) coincides with the β-gyres in the geostrophic model, and the second (ageostrophic β-gyres) is due to ageostrophic terms in the governing equations. The time evolution of the ageostrophic β-gyres consists of fast and slow stages. During the fast stage the radiation of inertia–gravity waves results in the rapid development of the β-gyres from zero to a dipole field independent of the fast time variable. Correspondingly, the vortex accelerates practically instantaneously (compared to the typical swirling time) to some finite value of the translation speed. At the next slow stage the inertia–gravity wave radiation is insignificant and the β-gyres evolve with the typical swirling time. The total zonal translation speed induced by the geostrophic and ageostrophic β-gyres tends with increasing time to the speed of a steadily translating monopole exceeding (not exceeding) the drift velocity of Rossby waves for anticyclones (cyclones). This cyclone/anticyclone asymmetry generalizes the well-known finding about the greater longevity of anticyclones compared to cyclones to the case of non-stationary evolving monopoles. The influence of inertia–gravity waves upon the vortex evolution is analysed. The main role of these waves is to provide a ‘fast’ adjustment to the ‘slow’ vortex evolution. The energy of inertia–gravity waves is negligible compare to the energy of the geostrophic β-gyres. Yet another feature of the ageostrophic vortex evolution is that the area of the potential vorticity patch changes in the course of time, the cyclonic patch contracting and the anticyclonic one expanding.


2009 ◽  
Vol 66 (4) ◽  
pp. 883-901 ◽  
Author(s):  
Michael L. Waite ◽  
Chris Snyder

Abstract The atmospheric mesoscale kinetic energy spectrum is investigated through numerical simulations of an idealized baroclinic wave life cycle, from linear instability to mature nonlinear evolution and with high horizontal and vertical resolution (Δx ≈ 10 km and Δz ≈ 60 m). The spontaneous excitation of inertia–gravity waves yields a shallowing of the mesoscale spectrum with respect to the large scales, in qualitative agreement with observations. However, this shallowing is restricted to the lower stratosphere and does not occur in the upper troposphere. At both levels, the mesoscale divergent kinetic energy spectrum—a proxy for the inertia–gravity wave energy spectrum—resembles a −5/3 power law in the mature stage. Divergent kinetic energy dominates the lower stratospheric mesoscale spectrum, accounting for its shallowing. Rotational kinetic energy, by contrast, dominates the upper tropospheric spectrum and no shallowing of the full spectrum is observed. By analyzing the tendency equation for the kinetic energy spectrum, it is shown that the lower stratospheric spectrum is not governed solely by a downscale energy cascade; rather, it is influenced by the vertical pressure flux divergence associated with vertically propagating inertia–gravity waves.


1998 ◽  
Vol 354 ◽  
pp. 301-343 ◽  
Author(s):  
OLIVER BÜHLER ◽  
MICHAEL E. McINTYRE

Idealized model examples of non-dissipative wave–mean interactions, using small-amplitude and slow-modulation approximations, are studied in order to re-examine the usual assumption that the only important interactions are dissipative. The results clarify and extend the body of wave–mean interaction theory on which our present understanding of, for instance, the global-scale atmospheric circulation depends (e.g. Holton et al. 1995). The waves considered are either gravity or inertia–gravity waves. The mean flows need not be zonally symmetric, but are approximately ‘balanced’ in a sense that non-trivially generalizes the standard concepts of geostrophic or higher-order balance at low Froude and/or Rossby number. Among the examples studied are cases in which irreversible mean-flow changes, capable of persisting after the gravity waves have propagated out of the domain of interest, take place without any need for wave dissipation. The irreversible mean-flow changes can be substantial in certain circumstances, such as Rossby-wave resonance, in which potential-vorticity contours are advected cumulatively. The examples studied in detail use shallow-water systems, but also provide a basis for generalizations to more realistic, stratified flow models. Independent checks on the analytical shallow-water results are obtained by using a different method based on particle-following averages in the sense of ‘generalized Lagrangian-mean theory’, and by verifying the theoretical predictions with nonlinear numerical simulations. The Lagrangian-mean method is seen to generalize easily to the three-dimensional stratified Boussinesq model, and to allow a partial generalization of the results to finite amplitude. This includes a finite-amplitude mean potential-vorticity theorem with a larger range of validity than had been hitherto recognized.


2021 ◽  
Author(s):  
Nedjeljka Žagar ◽  
Žiga Zaplotnik ◽  
Valentino Neduhal

<p>The energy spectrum of atmospheric horizontal motions has been extensively studied in observations and numerical simulations. Its canonical shape includes a transition from the -3 power law at synoptic scale to -5/3 power law at mesoscale. The transition is taking place at scales around 500 km that can be seen as the scale where energy associated with quasi-linear inertia-gravity waves exceeds the balanced (or Rossby wave) energy. In contrast to the horizontal spectrum, the spectrum of kinetic energy of vertical motions is poorly known since the vertical motion is not an observed quantity of the global observing system and vertical kinetic energy spectra from non-hydrostatic models are difficult to validate.</p> <p>Traditionally, vertical velocities associated with the Rossby and gravity waves have been treated separately using the quasi-geostrophic omega equations and polarization relations for the stratified Boussinesq fluid in the (x,z) plane, respectively. In the tropics, the Rossby and gravity  wave regimes are difficult to separate and their frequency gap, present in the extra-tropics, is filled with the Kelvin and mixed Rossby-gravity waves. A separate treatment of the Rossby and gravity wave regimes makes it challenging to quantify energies of their vertical motions and vertical momentum fluxes. A unified treatment and wave interactions is performed by high-resolution non-hydrostatic models but their understanding requires the toolkit of theory. </p> <p>This contribution presents a unified framework for the derivation of vertical velocities of the Rossby and inertia-gravity waves and associated kinetic energy spectra. Expressions for the Rossby and gravity wave vertical velocities are derived using the normal-mode framework in the hydrostatic atmosphere that can be considered applicable up to the scale around 10 km. The derivation involves the analytical evaluation of divergence of the horizontal wind associated with the Rossby and inertia-gravity eigensolutions of the linearized primitive equations. The new framework is applied to the global analysis data of the ECMWF system. Results confirm that the tropical vertical kinetic energy spectra associated with inertia-gravity waves are on average indeed white. Deviations from the white spectrum are discussed for latitude and altitude bands.</p>


2004 ◽  
Vol 34 (12) ◽  
pp. 2761-2773 ◽  
Author(s):  
Álvaro Viúdez ◽  
David G. Dritschel

Abstract A new method for diagnosing the balanced three-dimensional velocity from a given density field in mesoscale oceanic flows is described. The method is referred to as dynamic potential vorticity initialization (PVI) and is based on the idea of letting the inertia–gravity waves produced by the initially imbalanced mass density and velocity fields develop and evolve in time while the balanced components of these fields adjust during the diagnostic period to a prescribed initial potential vorticity (PV) field. Technically this is achieved first by calculating the prescribed PV field from given density and geostrophic velocity fields; then the PV anomaly is multiplied by a simple time-dependent ramp function, initially zero but tending to unity over the diagnostic period. In this way, the PV anomaly builds up to the prescribed anomaly. During this time, the full three-dimensional primitive equations—except for the PV equation—are integrated for several inertial periods. At the end of the diagnostic period the density and velocity fields are found to adjust to the prescribed PV field and the approximate balanced vortical motion is obtained. This adjustment involves the generation and propagation of fast, small-amplitude inertia–gravity waves, which appear to have negligible impact on the final near-balanced motion. Several practical applications of this method are illustrated. The highly nonlinear, complex breakup of baroclinically unstable currents into eddies, fronts, and filamentary structures is examined. The capability of the method to generate the balanced three-dimensional motion is measured by analyzing the ageostrophic horizontal and vertical velocity—the latter is the velocity component most sensitive to initialization, and one for which a quasigeostrophic diagnostic solution is available for comparison purposes. The authors find that the diagnosed fields are closer to the actual fields than are either the geostrophic or the quasigeostrophic approximations. Dynamic PV initialization thus appears to be a promising way of improving the diagnosis of balanced mesoscale motions.


2004 ◽  
Vol 61 (22) ◽  
pp. 2692-2706 ◽  
Author(s):  
Todd P. Lane ◽  
James D. Doyle ◽  
Riwal Plougonven ◽  
Melvyn A. Shapiro ◽  
Robert D. Sharman

Abstract The characteristics and dynamics of inertia–gravity waves generated in the vicinity of an intense jet stream/ upper-level frontal system on 18 February 2001 are investigated using observations from the NOAA Gulfstream-IV research aircraft and numerical simulations. Aircraft dropsonde observations and numerical simulations elucidate the detailed mesoscale structure of this system, including its associated inertia–gravity waves and clear-air turbulence. Results from a multiply nested numerical model show inertia–gravity wave development above the developing jet/front system. These inertia–gravity waves propagate through the highly sheared flow above the jet stream, perturb the background wind shear and stability, and create bands of reduced and increased Richardson numbers. These bands of reduced Richardson numbers are regions of likely Kelvin–Helmholtz instability and a possible source of the clear-air turbulence that was observed.


2014 ◽  
Vol 757 ◽  
pp. 817-853 ◽  
Author(s):  
Callum J. Shakespeare ◽  
J. R. Taylor

AbstractDensity fronts are common features of ocean and atmosphere boundary layers. Field observations and numerical simulations have shown that the sharpening of frontal gradients, or frontogenesis, can spontaneously generate inertia–gravity waves (IGWs). Although significant progress has been made in describing frontogenesis using approximations such as quasi-geostrophy (Stone, J. Atmos. Sci., vol. 23, 1966, pp. 455–565, Williams & Plotkin J. Atmos. Sci., vol. 25, 1968, pp. 201–206) semi-geostrophy (Hoskins, Annu. Rev. Fluid Mech., vol. 14, 1982, pp. 131–151), these models omit waves. Here, we further develop the analytical model of Shakespeare & Taylor (J. Fluid Mech., vol. 736, 2013, pp. 366–413) to describe the spontaneous emission of IGWs from an initially geostrophically balanced front subjected to a time-varying horizontal strain. The model uses the idealised configuration of an infinitely long, straight front and uniform potential vorticity (PV) fluid, with a uniform imposed convergent strain across the front, similar to Hoskins & Bretherton (J. Atmos. Sci., vol. 29, 1972, pp. 11–37). Inertia–gravity waves are generated via two distinct mechanisms: acceleration of the large-scale flow and frontal collapse. Wave emission via frontal collapse is predicted to be exponentially small for small values of strain but significant for larger strains. Time-varying strain can also generate finite-amplitude waves by accelerating the cross-front flow and disrupting geostrophic balance. In both cases waves are trapped by the oncoming strain flow and can only propagate away from the frontal zone when the strain field weakens sufficiently, leading to wave emission that is strongly localised in both time and space.


Sign in / Sign up

Export Citation Format

Share Document