2 Transport of heat, momentum and potential vorticity
|
|
- Alberta Day
- 6 years ago
- Views:
Transcription
1 Transport of heat, momentum and potential vorticity. Conventional mean momentum equation We ll write the inviscid equation of onal motion (we ll here be using log-pressure coordinates with = H ln p=p, with H constant) as and the continuity equation du dt fv = x () r u = () where () = p=gh. We introduce the onal mean of any variable a as a(y; ; t) = L Z L a(x; y; ; t) dx ; (3) (where we assume periodic behavior in x such that a(l; y; ; t) = a(; y; ; t) note that this means that (a=x) = ) and then the eddy component is It follows that, by de nition, a (x; y; ; t) = a(x; y; ; t) a(y; ; t) : a = : (4) Note that the de nition of the eddy variable is contingent on the de nition of mean. Then the mean of (), invoking (), becomes u t + v u y + w u f v = r u u ; (5) thus introducing the concept of eddy momentum ux u u = u v ; u w (it has no onal component in our onal mean formulation). The way in which the eddies in uence the mean state is evidently equivalent to a body force r u u acting on the mean onal motion. This ux term clearly represents the advection, by the eddy component of the meridional velocity, of the eddy component of onal momentum. The onal momentum per unit mass is u; so the advective ux of momentum across a coordinate surface (of constant y and ) is uu, whence its onal mean is uu = uu + u u ; hence it is natural to regard u u as the eddy ux of momentum. It is important to realise, however, that this de nition of eddy momentum ux is a consequence of our de nition (3) of what is meant by mean, speci cally, that the average is taken at at constant pressure (constant in log-p coordinates).
2 One could easily de ne other means: since transport is fundamentally a Lagrangian process, one might try taking some kind of Lagragian mean, but we will not do that here. One can take a di erent integration path in (3) the way we have done it is to integrate at constant latitude and height (i.e.; constant pressure) by using a di erent coordinate system, in which case the whole problem can look rather di erent, as we are about to see.. Wave, mean ow interaction in isentropic coordinates In (x; y; ) coordinates, the inviscid onal equation of motion is u t + uu x + v u y + _ u fv = M x (6) and that of continuity t + x (u) + y (v) + _ = (7) where = g p= is the isentropic coordinate density, and M = c p T + g = + (8) is the Montgomery potential, with (p) = c p (p=p ) the Exner function and = g. (Note that is real height here, not log-pressure height.) It follows from (8), with some juggling and the assumption of hydrostatic balance, that Now our de nition of onal mean is M = : (9) a(y; ; t) = L Z L a(x; y; ; t) dx where the integral is now, of course, taken at constant, and the eddy term is a (x; y; ; t) = a(x; y; ; t) a(y; ; t) : If you are interested in the Lagrangian perspective, there are some additional notes on this on the class web page. The conceptually powerful generalied Lagrangian mean theory of wave, mean ow, interaction can be found in Andrews & McIntyre, J. Fluid Mech., 978 (though it is not for the faint hearted!). The derivation of the -coordinate equations is discussed in Appendix 3A in Andrews, Holton and Leovy [987] (at the end of Chapter 3).
3 It makes sense to de ne mean meridional velocities in a density-weighted sense: v = (v) ; _ = _ () since then the mean continuity equation has no eddy ux terms: t + y (v ) + _ = : The momentum equation (6), on multiplication by, use of (7), and onal averaging, gives t (u)+ y (v u)+ _ u f v = y (v) u _ u M x : () The rst two terms on the rhs are easily identifable as advection by the eddy mass uxes, but what about the third term? Since M=x =, it is an eddy term M x = M x and, in fact, the term can be written, using the de nition of and (8) M x = p M g x = g pm x g p M x : But M= = (p), so p M x = p p (p) = pd x dp x = Z p d x dp dp = ; hence In total, therefore, () becomes t (u)+ y (v u)+ _ u M x = g pm = x f v = y h(v) u i M p g x : _ M u p g x ; () 3
4 there is an eddy momentum ux whose y and components are given by (v) u ; _ M u p : (3) g x The additional, non-advective, term in (3) represents the form stress acting on the isentropic surface. If one has, e.g., surface topography (x; y) there is an inviscid stress on the earth s surface given by psin = p=x, where 3 is the angle of the local surface with respect to the horiontal; accordingly there is a form drag of p=x acting on the atmosphere above. Now, p M x = p M x = p d p dp x + g x = p d p + gp dp x x since we just established that the rst term is ero. Hence M p g x = p x : = gp x ; Since =x here is just the slope of the isentropic surface, we can therefore immediately identify this extra term as the form drag acting on an isentropic surface due to the presence of the eddies. So, in a di erent coordinate system, eddy transport of momentum can look quite di erent, physically as well as mathematically..3 QG momentum, heat, and PV transport.3. The quasigeostrophic equations on a -plane We consider motion with characteristic horiontal length scale L, height scale h, velocity scale U, time scale L=U on a -plane for which the Coriolis parameter is f = f + y. We make the assumptions that: (i) the Rossby number Ro = U=f L is small, (ii) L=f Ro (which, in a spherical context, means that L << a, where a is the Earth s radius), (iii) the isentropic slopes j x j=j j and j y j=j j are Ro (h=l), (otherwise vertical motions would not be small), and 3 Under the shallow atmosphere" approximation which all these equations assume, the slope is small, so =x = tan ' sin : 4
5 (iv) the static stability is, to leading order, a function of only [this stems from (iii)], so we write = () + ~, where ~ is O(Ro) smaller than. We then de ne a background geopotential, in hydrostatic balance with the background : = g c p T = H : Under these assumptions, we nd that the leading order equations give geostrophic balance, in which we write the geostrophic velocities [i:e:, the leading order terms in a Rossby number expansion of (u; v; w)] as u = y ; v = x ; w = ; (4) where is the geostrophic streamfunction, which we de ne to be Hydrostatic balance then gives us = f H f = f H [ = [ ()] =f. (5) ()] = f H ~. (6) At next order, we obtain our quasigeostrophic equations. The equations of motion are D g u yv f v a = G (x), D g v + yu + f u a = G (y) (7), where G represents frictional or other forces per unit mass, D g is the time derivative following the geostrophic ow D g t + u x + v y and (u a ; v a ; w a ) is the ageostrophic velocity (i:e:, the di erence between the actual velocity and the geostrophic one). Similarly, the thermodynamic equation becomes D g + w a ; = c p Q = () J, (8) where J is the diabatic heating rate per unit volume (and Q = J =c p is the heating rate expressed in units of K per unit time). From these, we can readily derive the equation for quasigeostrophic potential vorticity, q : D g q = X, (9) where q = f + y + v x u y + f! ~ ; = f + y + () 5
6 and where with X = G (y) x x + y G y (x) + f + N = H ; J ; f N, (), () being the square of the buoyancy frequency. Eq. (9) tells us that, for conservative ow (G =, J =, whence X = ), q is conserved following the geostrophic ow. When the ow is not conservative, X represents the local sources or sinks of q arising from viscous and/or diabatic e ects..3. PV uxes and the Eliassen-Palm theorem Consider perturbations to a steady, onally-uniform basic state U = U(y; ) ; = (y; ) ; = (y; ) where The basic state PV is: y = U ; H y = f U. Q(y; ) = f + y + yy + f N = f + y +. (3) Write = + (x; y; ; t) where is a small perturbation. Now, v = x and q = = xx + yy + f N therefore v q = x xx + x yy + f x N Simple manipulation yields the following identities: ; (4) : (5) 6
7 . x xx = h ( x) i = ; x. x yy = x y y = x y = x y y y ; y xy h y i x and x f N = f N x = f N x = f N x : f N x f N h( ) i x Therefore, from (5), where F (y) F = = F () v q = r F. (6)! x y u f = v x =N f v : (7) = F is known as the ELIASSEN-PALM ux. Note that the northward component of F is (minus) the northward eddy ux of onal momentum, while the vertical component is proportional to the northward eddy ux of heat, v T : Now, assume that the perturbations are small, such that the QGPV equation (9) can be linearied to give ( t + U x ) q + v Q y = X ; multiply by q and average: If we de ne q + v q Q y = v X : (8) t A = q =Q y and D = v X =Q y ; 7
8 then, using (6), A t + r F = D. (9) Eq. (9) is the ELIASSEN-PALM RELATION. It is a conservation law for onally-averaged wave activity whose density is A. Note that D! for conservative ow. The signi cance of this relation is that it gives us a measure of the ux of wave activity through wave propagation. For example, if the waves are conservative (D = ) then A must increase with time wherever F is convergent and decrease wherever it is divergent. Thus F is a meaningful measure of the propagation of wave activity from one place to another. This becomes most obvious for almost-plane waves (WKB theory) when, as we shall see, where c g is group velocity, in which case F = c g A A t + r c g A = D. However, note that F remains valid as a measure of the ux of wave activity even when WKB theory is not valid and we cannot even de ne group velocity. [N.B. there are some subtleties to these arguments if Q y changes sign anywhere, since A is then not positive de nite (and so increasing A does not necessarily mean increasing wave amplitude). Indeed, if this occurs, the basic state may be unstable the Charney-Stern necessary condition for instability can in fact be readily obtained from (9).].3.3 The Eliassen-Palm theorem For waves which are steady (A t = ), of small amplitude and conservative (D = ), the ux F is nondivergent. This is, through (6), the same thing as saying that the northward ux of quasigeostrophic potential vorticity vanishes under these conditions [a result we could of course have obtained from (8) without involving F]:.3.4 Potential vorticity transport and the nonacceleration theorem We now consider the problem of how eddies (in this case, quasigeostrophic eddies) impact the onal mean circulation. We return to eddies which may not be small and to the QGPV budget, which becomes, on onal averaging q t + (v q ) y = X. (3) Note that, unlike (5), there is no mean advection term in (3). This is because there is no advection by w in this quasigeostrophic case and v = x =. Similarly, the eddy ux contains no vertical component, because of QG scalings. 8
9 The dynamical in uence of the eddies on mean potential vorticity, therefore, is described entirely 4 by the northward ux v q. Given the equivalence between pv q and rf, we could make the same statement about r F which informs us immediately that F is telling us something about wave transport as well as propagation. Now, we know from the Eliassen-Palm theorem that if the waves are everywhere (I) of steady amplitude, (II) conservative, and (III) of small amplitude, (so that terms of O( 3 ) in the wave activity budget are negligible) then F is nondivergent and v q =. Under these conditions, therefore, the budget equation for onally-averaged QGPV is independent of the waves (if we assume that X is so independent). Now, q = f + ( ) therefore t = qt = X : If q t is independent of the waves, then so is t ; since is an elliptic operator, the solution of the above equation for t invokes boundary conditions. If, however, we invoke the further condition that (IV) the boundary conditions on t are independent of the waves then t is everywhere independent of the waves. Since u = y and = (f T =g), it then follows that, under conditions (I)-(IV), the tendencies of the mean geostrophic ow, mean temperature and mean QGPV are independent of the waves. This is known as the nonacceleration theorem and conditions (I)-(IV) are sometimes known as nonacceleration conditions..3.5 Zonal mean momentum and heat budgets: Conventional Eulerian-mean approach In a similar fashion, we can take the onal mean of the quasigeostrophic momentum and heat equations. The onal mean of the rst of (7) gives us u t f v a = G (x) (u v ) y, (3) 4 One should add the caveat that, in principle, the eddies could also in uence X (e.g., eddies could modify mean precipitation). 9
10 where we have used v = x = (the mean ageostrophic wind is not ero, however). Similarly, the heat equation (8) gives t + w a = () J (v ) y. (3) The mean ow equations are closed by the continuity equation v a;y + (w a) = (33) and the thermal wind equation f u = H ~ y. (34) (The two latter equations are linear and so the onal averaging is trivial). The evolution of the onal mean state in the presence of eddies is therefore speci ed by (3) - (34). In this case, the e ects of wave transport are manifested in two terms the convergence of the eddy uxes of momentum u v and heat v (strictly, the eddy ux of heat is c p v T, but we shall refer to v in this way for simplicity). Both these terms force the mean ow equations and it is important to note that the whole system is coupled, i:e:, the heat uxes can impact on the mean winds just as much as can the momentum uxes. Thermal wind balance requires this to be true. Consider, for example, an upward propagating wave with v 6= but u v =. The mean state could not respond with a changing mean temperature only; thermal wind balance demands a corresponding change in u. From (3) this could only be achieved through an ageostrophic meridional circulation, which would impact on both the momentum and heat budgets. Thus, the waves will not only drive u t and t, but also v a and w a (except in the unlikely case where the eddy forcing terms conspire not to disturb thermal wind balance). To put the same statements into mathematics, what (3) - (34) give us is a set of 4 equations in the 4 unknowns u t, t, v a and w a in terms of the two eddy driving terms. In general, when both the eddy ux terms are nonero, there is no simple way of saying which forcing achieves what response. Moreover, the central role of the potential vorticity ux obvious in the PV budget is not at all obvious here. Indeed, we have seen from the mean potential vorticity budget that under nonacceleration conditions u t and t, must be ero under these conditions. What must (and does) happen under such circumstances is that the eddies induce ageostrophic mean motions whose e ects in (3) and (3) exactly balance the eddy ux terms. All in all, this approach to the mean momentum and heat budgets is not giving us much insight into what is going on.
11 .4 Transformed Eulerian-mean theory The most important lesson to learn from the isentropic-coordinate analysis we did earlier is that the de nition of mean, and consequently of concepts like eddy uxes, is non-unique. We can take this message to heart by allowing ourselves some exibility in de ning what we mean by the mean meridional circulation. We begin by noting that, from (33), we may de ne an ageostrophic mean streamfunction a such that (v a ; w a ) = ( a ), a y. (35) Now, we look for a more revealing way of de ning mean circulation. To make things as simple as possible, we insist that our modi ed circulation be nondivergent also, so we write ( (v ; w ) = ), ; (36) y where the new streamfunction is If we substitute this into (3), we obtain = a + c : t + w = () J y (v ) c y. Noting that = (), it follows that if we make the choice c = v, so that the streamfunction of the so-called residual circulation is we obtain the thermodynamic equation = a v, (37) t + w = () J. (38) We have thus succeeded in deriving a mean heat equation in which there are no explicit eddy terms; heat is transported solely through the mean vertical residual" motion. It might be thought, of course, that the eddy terms are still there, implicit in w ; but this was also true of w a which, as noted earlier, is in general in uenced by the eddies. What we have done is to rede ne this in uence so as to put the mean heat budget into its simplest possible form.
12 is We now need to complete our transformed system of equations. The continuity equation v y + (w ) = (39) [we arranged this by (36)]. The thermal wind equation stays as before, i:e: : f u = H ~ y : (4) The momentum equation is less trivial. We need to replace v a using (36); the result is u t f v = G x + r F, (4) where F is the Eliassen-Palm ux. This transformation which is nothing more than a di erent way of writing the same equations makes the role of the eddies look rather di erent. We now have, in (38) - (4), a set of equations for v ; w ; u=t and =t in which the only term representing eddy forcing is a term r F, which appears as an e ective body force (per unit mass) in the mean momentum equation. [We could, for example, rede ne G x to absorb this term]. It is clear, therefore, that under nonacceleration conditions when rf = and the boundary conditions are independent of wave-dependent terms, v, w, u=t and =t are independent of the waves. When nonacceleration conditions are not satis ed, the transformed equations o er a more transparent approach to the problem simply because the single term represented by the e ective body force rf entirely summaries the eddy forcing of the mean state (there being no thermal eddy forcing to confuse the issue). In fact, this formulation gives us another interpretation of F: as an eddy ux of (negative, i:e:, easterly) momentum which, because of the properties we have just described, is a more reliable measure of the wave transport of momentum than u v alone. Moreover, we now see that the Eliassen-Palm ux gives us a uni ed picture of wave propagation (through its role in wave activity conservation) and transport (through its interpretation as a momentum ux). This perspective is conceptually very powerful, as we shall see. Moreover, through the relationship (6) we immediately recover the central role of PV uxes in the eddy forcing of the mean state..4. Relationship to isentropic mean perspective Finally, note again that, while the perspective which regards F as a momentum ux may seem to be the result of mere mathematical juggling, it is rather more than that. It might seem odd to have what looks like an eddy heat ux appearing in the vertical momentum ux in this formulation. But it can be shown (see relevant Problem Set) that, for small-amplitude eddies under adiabatic conditions:
13 . this heat ux term corresponds to the form stress term that appeared in the isentropic coordinate formulation, and. that the residual circulation corresponds to the density-weighted mean circulation is isentropic coordinates. [It is also, under these assumptions, directly related to the Lagrangian mean motion.] Thus, aside from the mathematical advantages of the TEM set of equations over the conventional set, the formaultion is also describing the physics rather well. [You might well ask: why not use the isentropic coordinate version? The answer is that, for wave problems, isentropic coordinates are a bit unwieldy.].4. Boundary conditions on the TEM equations Finally, a short but important note about boundary conditions on the TEM ow. At a at lower boundary, the boundary condition w = translates straightforwardly, under conventional averaging, to w =. For the TEM equations, however, (36) and (37) tells us that! w = v y which is nonero whenever there is a heat ux along the boundary. So, in this respect, the TEM framework is the more complicated (though it is actually telling us something about the dynamics near the boundary see Problem Set). As we shall see, on a lower boundary with topography, the opposite may be true. In general w is nonero on such a boundary. 3
Transformed Eulerian Mean
Chapter 15 Transformed Eulerian Mean In the last few lectures we introduced some fundamental ideas on 1) the properties of turbulent flows in rotating stratified environments, like the ocean and the atmosphere,
More information8 3D transport formulation
8 3D transport formulation 8.1 The QG case We consider the average (x; y; z; t) of a 3D, QG system 1. The important distinction from the zonal mean case is that the mean now varies with x, or (more generally)
More informationEliassen-Palm Theory
Eliassen-Palm Theory David Painemal MPO611 April 2007 I. Introduction The separation of the flow into its zonal average and the deviations therefrom has been a dominant paradigm for analyses of the general
More information7 The General Circulation
7 The General Circulation 7.1 The axisymmetric state At the beginning of the class, we discussed the nonlinear, inviscid, axisymmetric theory of the meridional structure of the atmosphere. The important
More informationEliassen-Palm Cross Sections Edmon et al. (1980)
Eliassen-Palm Cross Sections Edmon et al. (1980) Cecily Keppel November 14 2014 Eliassen-Palm Flux For β-plane Coordinates (y, p) in northward, vertical directions Zonal means F = v u f (y) v θ θ p F will
More informationQuasi-geostrophic ocean models
Quasi-geostrophic ocean models March 19, 2002 1 Introduction The starting point for theoretical and numerical study of the three dimensional large-scale circulation of the atmosphere and ocean is a vorticity
More information1 The circulation of a zonally symmetric atmosphere. 1.1 Angular momentum conservation and its implications
1 The circulation of a zonally symmetric atmosphere We will begin our attempts to understand the big picture of the structure of the atmosphere by asking about what theory predicts if we ignore eddies
More informationChapter 2. Quasi-Geostrophic Theory: Formulation (review) ε =U f o L <<1, β = 2Ω cosθ o R. 2.1 Introduction
Chapter 2. Quasi-Geostrophic Theory: Formulation (review) 2.1 Introduction For most of the course we will be concerned with instabilities that an be analyzed by the quasi-geostrophic equations. These are
More information10 Shallow Water Models
10 Shallow Water Models So far, we have studied the effects due to rotation and stratification in isolation. We then looked at the effects of rotation in a barotropic model, but what about if we add stratification
More informationChapter 3. Stability theory for zonal flows :formulation
Chapter 3. Stability theory for zonal flows :formulation 3.1 Introduction Although flows in the atmosphere and ocean are never strictly zonal major currents are nearly so and the simplifications springing
More informationStratospheric Dynamics and Coupling with Troposphere and Mesosphere
WDS'13 Proceedings of Contributed Papers, Part III, 6 66, 13. ISBN 978-8-7378-5-8 MATFYZPRESS Stratospheric Dynamics and Coupling with Troposphere and Mesosphere P. Šácha Charles University in Prague,
More information1 Climatological balances of heat, mass, and angular momentum (and the role of eddies)
1 Climatological balances of heat, mass, and angular momentum (and the role of eddies) We saw that the middle atmospheric temperature structure (which, through thermal wind balance, determines the mean
More informationQuasi-Geostrophic ω-equation. 1. The atmosphere is approximately hydrostatic. 2. The atmosphere is approximately geostrophic.
Quasi-Geostrophic ω-equation For large-scale flow in the atmosphere, we have learned about two very important characteristics:. The atmosphere is approximately hydrostatic.. The atmosphere is approximately
More informationGoal: Use understanding of physically-relevant scales to reduce the complexity of the governing equations
Scale analysis relevant to the tropics [large-scale synoptic systems]* Goal: Use understanding of physically-relevant scales to reduce the complexity of the governing equations *Reminder: Midlatitude scale
More informationThe General Circulation of the Atmosphere: A Numerical Experiment
The General Circulation of the Atmosphere: A Numerical Experiment Norman A. Phillips (1956) Presentation by Lukas Strebel and Fabian Thüring Goal of the Model Numerically predict the mean state of the
More informationChapter 9. Geostrophy, Quasi-Geostrophy and the Potential Vorticity Equation
Chapter 9 Geostrophy, Quasi-Geostrophy and the Potential Vorticity Equation 9.1 Geostrophy and scaling. We examined in the last chapter some consequences of the dynamical balances for low frequency, nearly
More information6 Two-layer shallow water theory.
6 Two-layer shallow water theory. Wewillnowgoontolookatashallowwatersystemthathastwolayersofdifferent density. This is the next level of complexity and a simple starting point for understanding the behaviour
More informationIntroduction to Isentropic Coordinates:! a new view of mean meridional & eddy circulations" Cristiana Stan
Introduction to Isentropic Coordinates:! a new view of mean meridional & eddy circulations" Cristiana Stan School and Conference on the General Circulation of the Atmosphere and Oceans: a Modern Perspective!
More informationESCI 343 Atmospheric Dynamics II Lesson 11 - Rossby Waves
ESCI 343 Atmospheric Dynamics II Lesson 11 - Rossby Waves Reference: An Introduction to Dynamic Meteorology (4 rd edition), J.R. Holton Atmosphere-Ocean Dynamics, A.E. Gill Fundamentals of Atmospheric
More informationFour ways of inferring the MMC. 1. direct measurement of [v] 2. vorticity balance. 3. total energy balance
Four ways of inferring the MMC 1. direct measurement of [v] 2. vorticity balance 3. total energy balance 4. eliminating time derivatives in governing equations Four ways of inferring the MMC 1. direct
More information1/18/2011. Conservation of Momentum Conservation of Mass Conservation of Energy Scaling Analysis ESS227 Prof. Jin-Yi Yu
Lecture 2: Basic Conservation Laws Conservation Law of Momentum Newton s 2 nd Law of Momentum = absolute velocity viewed in an inertial system = rate of change of Ua following the motion in an inertial
More information2 A: The Shallow Water Equations
2 A: The Shallow Water Equations 2.1 Surface motions on shallow water Consider two-dimensional (x-z) motions on a nonrotating, shallow body of water, of uniform density, as shown in Fig. 1 below. The ow
More information1. tangential stresses at the ocean s surface due to the prevailing wind systems - the wind-driven circulation and
Chapter 9 The wind-driven circulation [Hartmann, Ch. 7 (7.4-7.5)] There are two causes of the circulation of the ocean: 1. tangential stresses at the ocean s surface due to the prevailing wind systems
More information2.5 Shallow water equations, quasigeostrophic filtering, and filtering of inertia-gravity waves
Chapter. The continuous equations φ=gh Φ=gH φ s =gh s Fig..5: Schematic of the shallow water model, a hydrostatic, incompressible fluid with a rigid bottom h s (x,y), a free surface h(x,y,t), and horizontal
More informationConservation of Mass Conservation of Energy Scaling Analysis. ESS227 Prof. Jin-Yi Yu
Lecture 2: Basic Conservation Laws Conservation of Momentum Conservation of Mass Conservation of Energy Scaling Analysis Conservation Law of Momentum Newton s 2 nd Law of Momentum = absolute velocity viewed
More information7 Balanced Motion. 7.1 Return of the...scale analysis for hydrostatic balance! CSU ATS601 Fall 2015
7 Balanced Motion We previously discussed the concept of balance earlier, in the context of hydrostatic balance. Recall that the balanced condition means no accelerations (balance of forces). That is,
More informationLecture #2 Planetary Wave Models. Charles McLandress (Banff Summer School 7-13 May 2005)
Lecture #2 Planetary Wave Models Charles McLandress (Banff Summer School 7-13 May 2005) 1 Outline of Lecture 1. Observational motivation 2. Forced planetary waves in the stratosphere 3. Traveling planetary
More information9 Rossby Waves. 9.1 Non-divergent barotropic vorticity equation. CSU ATS601 Fall (Holton Chapter 7, Vallis Chapter 5)
9 Rossby Waves (Holton Chapter 7, Vallis Chapter 5) 9.1 Non-divergent barotropic vorticity equation We are now at a point that we can discuss our first fundamental application of the equations of motion:
More informationChapter 9. Barotropic Instability. 9.1 Linearized governing equations
Chapter 9 Barotropic Instability The ossby wave is the building block of low ossby number geophysical fluid dynamics. In this chapter we learn how ossby waves can interact with each other to produce a
More informationIn two-dimensional barotropic flow, there is an exact relationship between mass
19. Baroclinic Instability In two-dimensional barotropic flow, there is an exact relationship between mass streamfunction ψ and the conserved quantity, vorticity (η) given by η = 2 ψ.the evolution of the
More informationDynamics of the Atmosphere. Large-scale flow with rotation and stratification
12.810 Dynamics of the Atmosphere Large-scale flow with rotation and stratification Visualization of meandering jet stream Upper level winds from June 10th to July 8th 1988 from MERRA Red shows faster
More informationATMOSPHERIC AND OCEANIC FLUID DYNAMICS
ATMOSPHERIC AND OCEANIC FLUID DYNAMICS Fundamentals and Large-scale Circulation G E O F F R E Y K. V A L L I S Princeton University, New Jersey CAMBRIDGE UNIVERSITY PRESS An asterisk indicates more advanced
More informationThe Eady problem of baroclinic instability described in section 19a was shown to
0. The Charney-Stern Theorem The Eady problem of baroclinic instability described in section 19a was shown to be remarkably similar to the Rayleigh instability of barotropic flow described in Chapter 18.
More informationDiagnosis of a Quasi-Geostrophic 2-Layer Model Aaron Adams, David Zermeño, Eunsil Jung, Hosmay Lopez, Ronald Gordon, Ting-Chi Wu
Diagnosis of a Quasi-Geostrophic 2-Layer Model Aaron Adams, David Zermeño, Eunsil Jung, Hosmay Lopez, Ronald Gordon, Ting-Chi Wu Introduction For this project we use a simple two layer model, which is
More informationAtmosphere, Ocean and Climate Dynamics Answers to Chapter 8
Atmosphere, Ocean and Climate Dynamics Answers to Chapter 8 1. Consider a zonally symmetric circulation (i.e., one with no longitudinal variations) in the atmosphere. In the inviscid upper troposphere,
More information7 The Quasi-Biennial Oscillation (QBO)
7 The Quasi-Biennial Oscillation (QBO) (Reviewed by Baldwin et al., Rev. Geophys., 001) Previously we noted the remarkable quasi-periodic reversal of zonal winds in the tropical stratosphere, the quasi-biennial
More informationROSSBY WAVE PROPAGATION
ROSSBY WAVE PROPAGATION (PHH lecture 4) The presence of a gradient of PV (or q.-g. p.v.) allows slow wave motions generally called Rossby waves These waves arise through the Rossby restoration mechanism,
More informationPAPER 333 FLUID DYNAMICS OF CLIMATE
MATHEMATICAL TRIPOS Part III Wednesday, 1 June, 2016 1:30 pm to 4:30 pm Draft 21 June, 2016 PAPER 333 FLUID DYNAMICS OF CLIMATE Attempt no more than THREE questions. There are FOUR questions in total.
More informationBALANCED FLOW: EXAMPLES (PHH lecture 3) Potential Vorticity in the real atmosphere. Potential temperature θ. Rossby Ertel potential vorticity
BALANCED FLOW: EXAMPLES (PHH lecture 3) Potential Vorticity in the real atmosphere Need to introduce a new measure of the buoyancy Potential temperature θ In a compressible fluid, the relevant measure
More informationChapter 5. Shallow Water Equations. 5.1 Derivation of shallow water equations
Chapter 5 Shallow Water Equations So far we have concentrated on the dynamics of small-scale disturbances in the atmosphere and ocean with relatively simple background flows. In these analyses we have
More informationNew variables in spherical geometry. David G. Dritschel. Mathematical Institute University of St Andrews.
New variables in spherical geometry David G Dritschel Mathematical Institute University of St Andrews http://www-vortexmcsst-andacuk Collaborators: Ali Mohebalhojeh (Tehran St Andrews) Jemma Shipton &
More information8 Mechanisms for tropical rainfall responses to equatorial
8 Mechanisms for tropical rainfall responses to equatorial heating More reading: 1. Hamouda, M. and Kucharski, F. (2019) Ekman pumping Mechanism driving Precipitation anomalies in Response to Equatorial
More informationUsing simplified vorticity equation,* by assumption 1 above: *Metr 430 handout on Circulation and Vorticity. Equations (4) and (5) on that handout
Rossby Wave Equation A. Assumptions 1. Non-divergence 2. Initially, zonal flow, or nearly zonal flow in which u>>v>>w. 3. Initial westerly wind is geostrophic and does not vary along the x-axis and equations
More informationOn the Instability of the African Easterly Jet and the Generation of AfricanWaves: Reversals of the Potential Vorticity Gradient
On the Instability of the African Easterly Jet and the Generation of AfricanWaves: Reversals of the Potential Vorticity Gradient Jen-Shan Hsieh * and Kerry H. Cook Department of Earth and Atmospheric Sciences,
More informationThe counter-propagating Rossby-wave perspective on baroclinic instability. I: Mathematical basis
Q. J. R. Meteorol. Soc. (2004), 130, pp. 211 231 doi: 10.1256/qj.02.184 The counter-propagating Rossby-wave perspective on baroclinic instability. I: Mathematical basis By E. HEIFETZ 1, C. H. BISHOP 2,
More information第四章 : 中纬度的经向环流系统 (III) 授课教师 : 张洋. - Ferrel cell, baroclinic eddies and the westerly jet
第四章 : 中纬度的经向环流系统 (III) - Ferrel cell, baroclinic eddies and the westerly jet 授课教师 : 张洋 2016. 10. 24 Outline Review! Observations! The Ferrel Cell!! Review: baroclinic instability and baroclinic eddy life
More informationTurbulence in the Atmosphere and Oceans
Turbulence in the Atmosphere and Oceans Instructors: Raffaele Ferrari and Glenn Flierl Course description The course will present the phenomena, theory, and modeling of turbulence in the Earth s oceans
More informationNote that Rossby waves are tranverse waves, that is the particles move perpendicular to the direction of propagation. f up, down (clockwise)
Ocean 423 Rossby waves 1 Rossby waves: Restoring force is the north-south gradient of background potential vorticity (f/h). That gradient can be due to either the variation in f with latitude, or to a
More informationOn the effect of forward shear and reversed shear baroclinic flows for polar low developments. Thor Erik Nordeng Norwegian Meteorological Institute
On the effect of forward shear and reversed shear baroclinic flows for polar low developments Thor Erik Nordeng Norwegian Meteorological Institute Outline Baroclinic growth a) Normal mode solution b) Initial
More information2. Baroclinic Instability and Midlatitude Dynamics
2. Baroclinic Instability and Midlatitude Dynamics Midlatitude Jet Stream Climatology (Atlantic and Pacific) Copyright 26 Emily Shuckburgh, University of Cambridge. Not to be quoted or reproduced without
More informationPV Thinking. What is PV thinking?
PV Thinking CH 06 What is PV thinking? Two main approaches to solving fluid flow problems:. We can integrate the momentum, continuity and thermodynamic equations (the primitive equations) directly.. In
More information( ) where the phase! is given by! = kx + mz!"t. We also know
GFD I, Final Exam Solutions 3/7/1 Parker MacCready 1.(a) The expression for the pressure perturbation is found from the vertical momentum equation: Z-MOM w t! 1! p' z b which may be rearranged to give:
More informationLecture 23. Vertical coordinates General vertical coordinate
Lecture 23 Vertical coordinates We have exclusively used height as the vertical coordinate but there are alternative vertical coordinates in use in ocean models, most notably the terrainfollowing coordinate
More informationA Hamiltonian Numerical Scheme for Large Scale Geophysical Fluid Systems
A Hamiltonian Numerical Scheme for Large Scale Geophysical Fluid Systems Bob Peeters Joint work with Onno Bokhove & Jason Frank TW, University of Twente, Enschede CWI, Amsterdam PhD-TW colloquium, 9th
More informationA mechanistic model study of quasi-stationary wave reflection. D.A. Ortland T.J. Dunkerton NorthWest Research Associates Bellevue WA
A mechanistic model study of quasi-stationary wave reflection D.A. Ortland T.J. Dunkerton ortland@nwra.com NorthWest Research Associates Bellevue WA Quasi-stationary flow Describe in terms of zonal mean
More informationThe Hydrostatic Approximation. - Euler Equations in Spherical Coordinates. - The Approximation and the Equations
OUTLINE: The Hydrostatic Approximation - Euler Equations in Spherical Coordinates - The Approximation and the Equations - Critique of Hydrostatic Approximation Inertial Instability - The Phenomenon - The
More informationGFD 2012 Lecture 1: Dynamics of Coherent Structures and their Impact on Transport and Predictability
GFD 2012 Lecture 1: Dynamics of Coherent Structures and their Impact on Transport and Predictability Jeffrey B. Weiss; notes by Duncan Hewitt and Pedram Hassanzadeh June 18, 2012 1 Introduction 1.1 What
More information( ) (9.1.1) Chapter 9. Geostrophy, Quasi-Geostrophy and the Potential Vorticity Equation. 9.1 Geostrophy and scaling.
Chapter 9 Geostrophy, Quasi-Geostrophy and the Potential Vorticity Equation 9.1 Geostrophy and scaling. We examined in the last chapter some consequences of the dynamical balances for low frequency, nearly
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationModeling Large-Scale Atmospheric and Oceanic Flows 2
Modeling Large-Scale Atmospheric and Oceanic Flows V. Zeitlin known s of the Laboratoire de Météorologie Dynamique, Univ. P. and M. Curie, Paris Mathematics of the Oceans, Fields Institute, Toronto, 3
More informationWave Mean-Flow Interactions in Moist Baroclinic Life Cycles
JULY 2017 Y A M A D A A N D P A U L U I S 2143 Wave Mean-Flow Interactions in Moist Baroclinic Life Cycles RAY YAMADA Department of Earth, Atmospheric and Planetary Sciences, Massachusetts Institute of
More informationChapter 1. Governing Equations of GFD. 1.1 Mass continuity
Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for
More informationIntroduction to Isentropic Coordinates: a new view of mean meridional & eddy circulations. Cristiana Stan
Introduction to Isentropic Coordinates: a new view of mean meridional & eddy circulations Cristiana Stan School and Conference on the General Circulation of the Atmosphere and Oceans: a Modern Perspective
More informationGravity Waves. Lecture 5: Waves in Atmosphere. Waves in the Atmosphere and Oceans. Internal Gravity (Buoyancy) Waves 2/9/2017
Lecture 5: Waves in Atmosphere Perturbation Method Properties of Wave Shallow Water Model Gravity Waves Rossby Waves Waves in the Atmosphere and Oceans Restoring Force Conservation of potential temperature
More informationModelling the atmosphere. Hennie Kelder University of Technology Eindhoven
Modelling the atmosphere Hennie Kelder University of Technology Eindhoven Content Basics of the atmosphere Atmospheric dynamics Large scale circulation Planetary waves Brewer-Dobson circulation Some Key
More informationLecture 10a: The Hadley Cell
Lecture 10a: The Hadley Cell Geoff Vallis; notes by Jim Thomas and Geoff J. Stanley June 27 In this short lecture we take a look at the general circulation of the atmosphere, and in particular the Hadley
More informationCirculation and Vorticity
Circulation and Vorticity Example: Rotation in the atmosphere water vapor satellite animation Circulation a macroscopic measure of rotation for a finite area of a fluid Vorticity a microscopic measure
More informationDefinition of a Generalized Diabatic Circulation Based on a Variational Approach
ISSN 000-4338, Izvestiya, Atmospheric and Oceanic Physics, 2007, Vol. 43, No. 4, pp. 436 44. Pleiades Publishing, Ltd., 2007. Original Russian Text A.S. Medvedev, 2007, published in Izvestiya AN. Fizika
More informationQuasi-geostrophic system
Quasi-eostrophic system (or, why we love elliptic equations for QGPV) Charney s QG the motion of lare-scale atmospheric disturbances is overned by Laws of conservation of potential temperature and potential
More informationThe Hadley Circulation and the Weak Temperature Gradient Approximation
1744 JOURNAL OF THE ATMOSPHERIC SCIENCES The Hadley Circulation and the Weak Temperature Gradient Approximation L. M. POLVANI AND A. H. SOBEL Department of Applied Physics and Applied Mathematics, and
More informationDynamics Rotating Tank
Institute for Atmospheric and Climate Science - IACETH Atmospheric Physics Lab Work Dynamics Rotating Tank Large scale flows on different latitudes of the rotating Earth Abstract The large scale atmospheric
More informationt tendency advection convergence twisting baroclinicity
RELATIVE VORTICITY EQUATION Newton s law in a rotating frame in z-coordinate (frictionless): U + U U = 2Ω U Φ α p U + U U 2 + ( U) U = 2Ω U Φ α p Applying to both sides, and noting ω U and using identities
More information3.1.1 The Primitive Equations in Log-Pressure Coordinates on the Sphere
Chapter 3 Basic Dynamics 3.1 Introduction The atmosphere is a shallow envelope of compressible gas surrounding an approximately spherical, rotating planet. The equations of motion in a rotating frame for
More informationρ x + fv f 'w + F x ρ y fu + F y Fundamental Equation in z coordinate p = ρrt or pα = RT Du uv tanφ Dv Dt + u2 tanφ + vw a a = 1 p Dw Dt u2 + v 2
Fundamental Equation in z coordinate p = ρrt or pα = RT Du uv tanφ + uw Dt a a = 1 p ρ x + fv f 'w + F x Dv Dt + u2 tanφ + vw a a = 1 p ρ y fu + F y Dw Dt u2 + v 2 = 1 p a ρ z g + f 'u + F z Dρ Dt + ρ
More informationwhere p oo is a reference level constant pressure (often 10 5 Pa). Since θ is conserved for adiabatic motions, a prognostic temperature equation is:
1 Appendix C Useful Equations Purposes: Provide foundation equations and sketch some derivations. These equations are used as starting places for discussions in various parts of the book. C.1. Thermodynamic
More informationNonlinear Balance on an Equatorial Beta Plane
Nonlinear Balance on an Equatorial Beta Plane David J. Raymond Physics Department and Geophysical Research Center New Mexico Tech Socorro, NM 87801 April 26, 2009 Summary Extension of the nonlinear balance
More informationLecture 12: Angular Momentum and the Hadley Circulation
Lecture 12: Angular Momentum and the Hadley Circulation September 30, 2003 We learnt last time that there is a planetary radiative drive net warming in the tropics, cooling over the pole which induces
More informationThe Effect of a Hadley Circulation on the Propagation and Reflection of Planetary Waves in a Simple One-Layer Model
1536 JOURNAL OF THE ATMOSPHERIC SCIENCES The Effect of a Hadley Circulation on the Propagation and Reflection of Planetary Waves in a Simple One-Layer Model J. GAVIN ESLER Program in Atmospheres, Oceans
More informationChapter 1. Introduction
Chapter 1. Introduction In this class, we will examine atmospheric phenomena that occurs at the mesoscale, including some boundary layer processes, convective storms, and hurricanes. We will emphasize
More information1/27/2010. With this method, all filed variables are separated into. from the basic state: Assumptions 1: : the basic state variables must
Lecture 5: Waves in Atmosphere Perturbation Method With this method, all filed variables are separated into two parts: (a) a basic state part and (b) a deviation from the basic state: Perturbation Method
More informationVacillations in a Shallow-Water Model of the Stratosphere
1174 JOURNAL OF THE ATMOSPHERIC SCIENCES VOLUME 61 Vacillations in a Shallow-Water Model of the Stratosphere PING-PING RONG AND DARRYN W. WAUGH Department of Earth and Planetary Sciences, Johns Hopkins
More informationChapter 5. Sound Waves and Vortices. 5.1 Sound waves
Chapter 5 Sound Waves and Vortices In this chapter we explore a set of characteristic solutions to the uid equations with the goal of familiarizing the reader with typical behaviors in uid dynamics. Sound
More informationOcean dynamics: the wind-driven circulation
Ocean dynamics: the wind-driven circulation Weston Anderson March 13, 2017 Contents 1 Introduction 1 2 The wind driven circulation (Ekman Transport) 3 3 Sverdrup flow 5 4 Western boundary currents (western
More informationControl Volume. Dynamics and Kinematics. Basic Conservation Laws. Lecture 1: Introduction and Review 1/24/2017
Lecture 1: Introduction and Review Dynamics and Kinematics Kinematics: The term kinematics means motion. Kinematics is the study of motion without regard for the cause. Dynamics: On the other hand, dynamics
More informationLecture 1: Introduction and Review
Lecture 1: Introduction and Review Review of fundamental mathematical tools Fundamental and apparent forces Dynamics and Kinematics Kinematics: The term kinematics means motion. Kinematics is the study
More informationDynamics and Kinematics
Geophysics Fluid Dynamics () Syllabus Course Time Lectures: Tu, Th 09:30-10:50 Discussion: 3315 Croul Hall Text Book J. R. Holton, "An introduction to Dynamic Meteorology", Academic Press (Ch. 1, 2, 3,
More informationLecture 3. Turbulent fluxes and TKE budgets (Garratt, Ch 2)
Lecture 3. Turbulent fluxes and TKE budgets (Garratt, Ch 2) The ABL, though turbulent, is not homogeneous, and a critical role of turbulence is transport and mixing of air properties, especially in the
More informationGeophysics Fluid Dynamics (ESS228)
Geophysics Fluid Dynamics (ESS228) Course Time Lectures: Tu, Th 09:30-10:50 Discussion: 3315 Croul Hall Text Book J. R. Holton, "An introduction to Dynamic Meteorology", Academic Press (Ch. 1, 2, 3, 4,
More informationg (z) = 1 (1 + z/a) = 1
1.4.2 Gravitational Force g is the gravitational force. It always points towards the center of mass, and it is proportional to the inverse square of the distance above the center of mass: g (z) = GM (a
More information3. Midlatitude Storm Tracks and the North Atlantic Oscillation
3. Midlatitude Storm Tracks and the North Atlantic Oscillation Copyright 2006 Emily Shuckburgh, University of Cambridge. Not to be quoted or reproduced without permission. EFS 3/1 Review of key results
More informationChapter 4. Gravity Waves in Shear. 4.1 Non-rotating shear flow
Chapter 4 Gravity Waves in Shear 4.1 Non-rotating shear flow We now study the special case of gravity waves in a non-rotating, sheared environment. Rotation introduces additional complexities in the already
More informationGeostrophic and Quasi-Geostrophic Balances
Geostrophic and Quasi-Geostrophic Balances Qiyu Xiao June 19, 2018 1 Introduction Understanding how the atmosphere and ocean behave is important to our everyday lives. Techniques such as weather forecasting
More information2 Equations of Motion
2 Equations of Motion system. In this section, we will derive the six full equations of motion in a non-rotating, Cartesian coordinate 2.1 Six equations of motion (non-rotating, Cartesian coordinates)
More informationQuasi-equilibrium Theory of Small Perturbations to Radiative- Convective Equilibrium States
Quasi-equilibrium Theory of Small Perturbations to Radiative- Convective Equilibrium States See CalTech 2005 paper on course web site Free troposphere assumed to have moist adiabatic lapse rate (s* does
More informationA Truncated Model for Finite Amplitude Baroclinic Waves in a Channel
A Truncated Model for Finite Amplitude Baroclinic Waves in a Channel Zhiming Kuang 1 Introduction To date, studies of finite amplitude baroclinic waves have been mostly numerical. The numerical models,
More informationThe Euler Equation of Gas-Dynamics
The Euler Equation of Gas-Dynamics A. Mignone October 24, 217 In this lecture we study some properties of the Euler equations of gasdynamics, + (u) = ( ) u + u u + p = a p + u p + γp u = where, p and u
More informationAPPENDIX B. The primitive equations
APPENDIX B The primitive equations The physical and mathematical basis of all methods of dynamical atmospheric prediction lies in the principles of conservation of momentum, mass, and energy. Applied to
More informationAtmospheric dynamics and meteorology
Atmospheric dynamics and meteorology B. Legras, http://www.lmd.ens.fr/legras III Frontogenesis (pre requisite: quasi-geostrophic equation, baroclinic instability in the Eady and Phillips models ) Recommended
More informationThe direct stratosphere troposphere interaction
Submitted to the Journal of Atmospheric Sciences special issue on Jets, June 2006 The direct stratosphere troposphere interaction MAARTEN H. P. AMBAUM 1 Department of Meteorology, University of Reading,
More information2. Meridional atmospheric structure; heat and water transport. Recall that the most primitive equilibrium climate model can be written
2. Meridional atmospheric structure; heat and water transport The equator-to-pole temperature difference DT was stronger during the last glacial maximum, with polar temperatures down by at least twice
More information