Three-dimensional baroclinic instability of a Hadley cell for small Richardson number

Abstract

A three-dimensional linear stability analysis of a baroclinic flow for Richardson number Ri of order unity is presented. The model considered is a thin, horizontal, rotating fluid layer which is subjected to horizontal and vertical temperature gradients. The basic state is a Hadley cell which is a solution of the Navier–Stokes and energy equations and contains both Ekman and thermal boundary layers adjacent to the rigid boundaries; it is given in closed form. The stability analysis is also based on the Navier–Stokes and energy equations; and perturbations possessing zonal, meridional and vertical structures were considered. Numerical methods were developed for the solution of the stability problem, which results in an ordinary differential eigenvalue problem. The objectives of this work were to extend the previous theoretical work on three-dimensional baroclinic instability for small Ri to a more realistic model involving the Prandtl number σ and the Ekman number E, and to finite growth rates and a wider range of the zonal wavenumber. The study covers ranges of 0.135 [les ] Ri [les ] 1.1, 0.2 [les ] σ [les ] 5.0, and 2 × 10−4 [les ] E [les ] 2 σ 10−3. For the cases computed for E = 10−3 and σ ≠ 1, we found that conventional baroclinic instability dominates for Ri > 0.825 and symmetric baroclinic instability dominates for Ri < 0.675. However, for E [ges ] 5 × 10−4 and σ = 1 in the range 0.3 [les ] Ri [les ] 0.8, conventional baroclinic instability always dominates. Further, we found in general that the symmetric modes of maximum growth are not purely symmetric but have weak zonal structure. This means that the wavefronts are inclined at a small angle to the zonal direction. The results also show that as E decreases the zonal structure of the symmetric modes of maximum growth rate also decreases. We found that when zonal structure is permitted the critical Richardson number for marginal stability is increased, but by only a small amount above the value for pure symmetric instability. Because these modes do not substantially alter the results for pure symmetric baroclinic instability and because their zonal structure is weak, it is unlikely that they represent a new type of instability.