Baroclinic Instability 1
Total Page:16
File Type:pdf, Size:1020Kb
BAROCLINIC INSTABILITY 1 BAROCLINIC INSTABILITY R Grotjahn, University of California, Davis, CA, USA Mathematical Formulation Copyright 2002 Elsevier Science Ltd. All Rights Reserved. The model uses quasi-geostrophic (QG) approxima- 0076-P0020 doi:10.1006/rwas.2002.0076 tions and nondimensional scaling appropriate for midlatitude frontal cyclones. Potential vorticity has Grotjahn, R contributions from the interior and from temperature University of California, gradients at rigid bottom ðz ¼ 0Þ and top ðz ¼ ZTÞ Department of Land, Air, and Water Resources, boundaries. In the QG system, PV can be written as: Davis, CA 95616-8627, USA 2 q ¼r|ffl{zffl}c þf0 þ by RV 8 9 Introduction > > <> => ½1a qc q2c qc qc 0076-P0005 Baroclinic instability refers to a process by which þ g þ e þ e Àe qz qz2 > qz qz > perturbations draw energy from the mean flow |fflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflffl} :>|fflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl}z¼0 z¼ZT ;> potential energy. The conversions of energy are TV BPV proportional to perturbation heat fluxes in both horizontal and vertical directions. The atmosphere requires heat fluxes to maintain the observed pattern where of net radiation (positive in the tropics, negative poleward of 381 N or S on an annual average). A zonal f 2L2 qlny mean meridional circulation, such as a tropical Hadley e ¼ 0 ; k ¼ D s ; and gkD qz cell, can generate these heat fluxes. However,in middle ½1b e qr qe latitudes, various factors cause eddies to accomplish g ¼ þ the bulk of the heat transport. Baroclinic instability r qz qz provides a mechanism to explain how these eddies form and evolve while incorporating the necessary where c is the horizontal velocity streamfunction, r is heat fluxes. Theoretical models of baroclinic instabil- density, g is the acceleration due to gravity, and k is the ity can simulate various observed properties of mid- static stability from the horizontal mean potential latitude eddies, including the dominant length scales, temperature. The coordinates are x eastward, y propagation speed, vertical structure, and energetics. northward, and z upward. Nondimensional length 0076-P0010 Baroclinic instability can be viewed as a shear scales are L in the horizontal and D in the vertical. f is instability. From thermal wind balance, the vertical 0 the constant part while b is the meridional derivative shear of the zonal wind is proportional to a meridional (approximated as a constant) of the Coriolis param- temperature gradient. The meridional temperature eter. gradient is proportional to the available potential An inherent horizontal length scale is the Rossby 0076-P0025 energy (APE) that the baroclinic instability mechanism radius of deformation L NHf À1 where N is the taps. Another view of baroclinic instability emphasiz- ð R ¼ 0 Þ 2 À1 es interacting potential vorticity (PV) anomalies. Brunt Va¨isa¨la¨ frequency ðN ¼ gkD Þ and H ¼ À1 Baroclinic instability is usually studied by linearizing RTg is the scale height (an inherent vertical length 2 À2 the dynamics equations and using eigenvalue or initial scale). Thus, e ¼ðLHÞ ðLRHÞ relates the assumed value techniques. These alternative views and analysis scales L and D to LR and H. procedures generally provide complementary means Quasi-geostrophic PV includes three distinct parts: 0076-P0030 to understand better baroclinic instability. absolute vorticity, which includes relative vorticity (RV); ‘‘thermal’’ vorticity (TV); and boundary PV (BPV). Positive PVis associated with an interior trough (in geopotential) and/or a warm surface temperature An Illustrative Model anomaly. When the vorticity and potential temperature con- 0076-P0035 0076-P0015 An illustrative model provides mathematical relations servation equations are combined, one obtains a time- and archetype solutions for the concepts that follow. dependent equation for QGPV conservation: rwas.2002.0076 15/4/02 13:25 Ed:: M. SHANKAR No. of pages: 10 Pgn:: thilakam 2 BAROCLINIC INSTABILITY ZZZ () q q 1 q qc 2 þ U r2c þ re qAe q rS qc qt qx r qz qz e dx dy dz ½2a qt qt 2 qz qQ qc ZZZ þ ¼ 0 qU qc qc qy qx ¼ er dx dy dz S qz qx qz |fflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl} ½5a with boundary conditions at the bottom and top ZZZ ðAz!AeÞ qc À r w dx dy dz S qz q q qc qU qc |fflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl} þ U À ¼ 0at qt qx qz qz qx ½2b ðKe!AeÞ z ¼ 0; ZT ZZZ () qK q r qc 2 qc 2 e S þ dx dy dz ‘‘Basic state’’ variables are specified: U (independent of qt qt 2 qy qx ZZZ x) is zonal wind, and Q is the interior part of the qc qc ¼À r U dx dy dz QGPV; meridional and vertical velocities are zero. S qx qy One can solve eqn [2] as an initial value problem by |fflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl}y specifying an initial streamfunction or potential ZZZ ðKz!KeÞ vorticity. qc þ r w dx dy dz 0076-P0040 An eigenvalue problem can also be formulated from S |fflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl}qz eqn [2]. A common approach assumes time and space ðA !K Þ dependence in the form: e e ½5b cðx; y; z; tÞ¼Re fgfðy; zÞ exp ½ikðx À ctÞ ½3 The volume integrals are over a closed domain. In the QG system, qc=qz is proportional to potential tem- for the ‘‘perturbation’’ streamfunction being sought. perature y, making the first term on the right-hand side This solution has zonal wavenumber k and complex of eqn [5a] proportional to a meridional heat flux, phase speed c. The growth rate is given by kImfcg.IfU while the second term is proportional to a vertical heat has no meridional variation, then one can assume a flux. The specified vertical shear, qU=qz, is propor- wavelike y dependence too: expðilyÞ. When wave- tional to the available potential energy, Az of the basic number l ¼ k, the solution is a ‘‘square wave.’’ state and is the energy upon which the baroclinic Perturbation velocities are defined as u ¼Àqc=qy instability mechanism feeds. The first term on the and v ¼Àqc=qx. right-hand side of eqn [5b] is a barotropic energy 0076-P0045 Additional simplifying approximations are often conversion. The barotropic conversion is proportional made. A particularly simple form, commonly labelled to the divergence of eddy momentum flux and also the ‘‘Eady model,’’ was described by E. T. Eady in draws energy from the mean flow. The second term on 1949. The Eady model assumes wavelike meridional the right-hand side of eqns [5a] and [5b] is the same, structure, qQ=qy ¼ 0, U ¼ z, incompressibility but with opposite sign indicating a conversion be- ðr ¼ constantÞ, and e ¼ 1. Then eqn [2a] is reduced tween Ae and Ke. simply to solving q0 ¼ 0 in the interior where the prime Example Solutions denotes the ‘‘perturbation’’ sought. The Eady eigen- value problem can be solved analytically to yield a pair This QG eigenmodel of baroclinic instability is appli- 0076-P0055 of normal modes, one growing, and one decaying, for cable to the midlatitudes. In these regions zonal flow scaled wavenumber ao 2:4. The scaled wave- increases with height reaching a maximum near the number: tropopause. Figure 1A is a representative nondimen- sional profile of U where the tropopause is at nondi- ÈÉ 2 2 À1 1=2 mensional z ¼ 1:0. The growth rate and phase speed a ¼ðk þ l Þe ½4 spectra, along with the (growing normal mode) eigenfunction structures for different k are also shown is proportional to absolute wavenumber and static in Figure 1. The growth rate has maximum value at a stability. specific value of a. The vertical structure tends to have 0076-P0050 Equations for perturbation kinetic energy, Ke and relative maxima at the surface and near the tropo- available potential energy, Ae are: pause, but it becomes progressively more bottom- rwas.2002.0076 15/4/02 13:25 Ed:: M. SHANKAR No. of pages: 10 Pgn:: thilakam BAROCLINIC INSTABILITY 3 k Im (c0) Re (c0) 2.0 0.40 0.7 1.6 0.30 0.5 1.2 0.20 z 0.8 0.3 U 0.10 0.4 Q0y 0.1 0.0 0 0 0.0 2.0 4.0 6.0 0 1.0 2.0 3.0 0 1.0 2.0 3.0 (A) (B) (C) = 0.8 = 2.0 = 3.0 2.0 2.0 2.0 A P P A A 1.6 1.6 1.6 P A P 1.2 1.2 P A 1.2 A P A z 0.8 P 0.8 0.8 P A A A P 0.4 0.4 0.4 P P P A A A 0.0 0.0 0.0 −2.0 −1.0 0 1.0 2.0 −2.0 −1.0 0 1.0 2.0 −2.0 −1.0 0 1.0 2.0 (D) (E) (F) 0076-F0001 Figure 1 Quasigeostrophic eigenanalysis. (A) Specified zonal wind U, and meridional gradient of interior potential vorticity Q0y versus scaled height. z ¼ 1 is 10 km. (B) Growth rate and (C) phase speed versus absolute wavenumber a. (D)–(F) Amplitude A, and phase P, for the growing normal mode for a ¼ 0:8, a ¼ 2:0, and a ¼ 3:0, respectively. All three modes tilt westward (upstream) with increasing height. Dimensional wavelengths depend upon scaling assumptions, but reasonable choices imply that a ¼ 0:8, a ¼ 2:0, and a ¼ 3:0 correspond to 11.0, 4.4, and 2:9Â103 km wavelengths, respectively. (Zonal and meridional scales are set equal.) The same scaling implies phase speed of 9 m s À 1 and doubling time of B1.2 days for a ¼ 2:0.