Course , General Circulation of the Earth's Atmosphere Prof. Peter Stone Section 8: Lorenz Energy Cycle
|
|
- Logan Nichols
- 6 years ago
- Views:
Transcription
1 Course.8, General Circulation of the Earth's Atmosphere Prof. Peter Stone Section 8: Lorenz Enery Cycle Enery Forms: As we saw in our discussion of the heat budet, the enery content of the atmosphere per unit mass can be expressed as the sum of four terms: E = I + Φ + LH + K I = C V T = internal enery Φ = z = potential enery K = (u + v + w ) (u + v ) = kinetic enery LH L V q = latent heat. We assume C V and L V are constants. If we assume hydrostatic equilibrium, then for the total internal enery in the column of the atmosphere we have: ρidz = ρc V Tdz = p s (z=) C V T dp; for the potential enery we have ps (z = ) ps (z=) ρφdz = ρzdz = zdp = d(zp) pdz ps RT = pdz = ρrtdz = dp ; (we choose z = to be a level below the lowest surface, z s min, and for convenience, in the presence of toporaphy, we take p = p s for z < z s. We will see later why. Note that the zero point for potential enery is arbitrary.) p s (z=) + R p s (z=) Thus ρ(i + Φ)dz = ( C V ) dp = C P T dp. Therefore in a column of the atmosphere the sum of I.E. and P.E. equals the enthalpy, and all three are equal within a constant factor. This is also sometimes called the total potential enery = TPE. This can be re-written in terms of the potential temperature, θ = T p p p Tdp = θ p R where = : 7 C P dp = θ ( + ) p dp+ ; θdp + = d ( θp + ) p + dθ;
2 interatin by parts we find for the total potential enery: p s (z=) C Tdp = C P P (+ ) p θp + p=ps,z= p s (z=) p= p + dθ What are the limits on θ? We now adopt the convention that below the surface, z < z s, θ < θ S = surface θ, and p = p s = constant, and θ at z =. (We implicitly assume that dθ >, i.e. the system is always stably stratified.) dz θ p + p s CP = and + ( ) p p + dθ = TPE However, we note that in a stably stratified atmosphere like ours much of the TPE can never be released or converted to KE. In particular, consider a typical arranement of the isentropic surfaces as shown in the first diaram below. Fiure by MIT OCW. θ 6 > θ 5... > θ. As we see there are horizontal radients and because of this PE can be released by the motions in spite of the stable stratification. Consider for example a displacement like that shown below.
3 Fiure by MIT OCW. For this displacement hih θ air is bein introduced into lower θ surroundins, so θ' >, but w ' > also, i.e., w 'θ' >, warm air is risin, and PE is bein converted into KE This is of course the mechanism of baroclinic instability. Now suppose we re-arrane the isentropes adiabatically (so no new enery is added) and while conservin mass, in such a way as to eliminate any horizontal radients. Since in H.E. p is just the amount of mass above any iven level, conservin mass means the interated or averae p above any θ surface is conserved. The averae p is p! ( θ) = p x,y,θ)da and usin this definition, one can now construct the state with no A ( horizontal radients, as shown below. 3
4 ! + >, ( p ) Fiure by MIT OCW. In this new state, since θ >, there are no displacements which can release PE, z because now for any displacement, w 'θ' <. Thus the PE of this state is unavailable. This led Lorenz (955) to define the available potential enery as the difference between the TPE in the actual atmosphere and that in the adjusted state just described. If we call this P, then for a iven system or reion, C P = P dxdy ( + ) p A Since p! is already interated over the area, we can rewrite the interal as + + A dθ ( p! p" ). This is the so-called exact formula. Note that, since dθ ( p + p! + ) > ( p" ) +, i.e. P is accordin to our definition always positive. Lorenz derived an approximation for P that is commonly used. Let p = p! + p', and similarly for other variables: 4
5 p + = p! + + p' ( + = p! + +) p' ( +) p' p! p! p! " ( +) p # ' p + = p! p! C P A ( +) P + ( ) p p! + p # ' p! p A dθ = C P dxdy p! p' dθ As a test of the accuracy of this approximation, we can take from Oort (97) the followin numbers for θ = 3 K in the annual mean: p ( equator) = 98mb, p ( 7N) = 435mb, p! ~ 78mb, p ' ~ 73mb, and the ratio of the st nelected term in the series to the first retained term is ( ) p'.4 73 = =.5 3 p! 3 78 Carryin out an analysis in θ coordinates is not very convenient, because we would have to interpolate from pressure coordinates, which would introduce error. Lorenz introduced another approximation. This is based on the approximation that both isentropes and pressure surfaces are quasi-horizontal with very small slopes: θ 4 z ο y = ~ 4 km ~ 3 and y θ 5 ο θ z km z p y ρƒu 4 = ~ ~ ~ 4 y p. p ρ z Consider for example an isentrope θ on which p has an averae p! ; now consider a point A on θ = constant. At this point, there will be a pressure deviation p'. 5
6 Fiure by MIT OCW. Now consider the θ variations on p! at the same latitude, φ. Since the θ variations are primarily vertical, or equivalently, vary primarily with pressure, we can write: p' p! θ'=. The minus sin is necessary since <, (Note that as drawn, p', θ' < ).! p θ θ Substitutin this into our interand for P, and chanin coordinates from θ to p! : dθ = d!p!p θ ( ), before interatin over x,y we have that the interand C P p! = p! θ' dp;! Note that θ' is now the variation of θ alon a constant p surface, p θ p since θ = T. p θ' T' θ θ = ; Also + = (Γ d Γ ) (use H.E.) where Γ d = θ T P Γ d p p! = θ Γ = d p! p, and θ = T ; θ p! ( Γ Γ) θ p d Substitutin this, the interand simplifies to C P p!p =!p p Γ d Γ d Γ!p Γ d Γ!p θ θ T T' d!p p!p T' T d!p C P. 6
7 T' = T (Γ Γ ) dp;! Now we substitute into the interation and absorb the minus sin by d reversin the limits of interation: P = A dxdy P S T' dp! T (Γ Γ ). d Finally Lorenz assumed that the coefficient of T' in the interand varied slowly compared to T', and replaced T and Γ by their area weihted means T! and Γ!. (i.e., 3 rd order quantities are aain bein nelected.) This is not as ood an approximation as the earlier ones. (From Oort s data I estimate ~% error, primarily from T.) Once this assumption is made we can write (notin that now p! is a dummy variable) p S P = dp T! (Γ d Γ ) A dxdyt'. Note that dxdyt' = dxdy T T! ) A A ( = dxdy ( T T T! + T! ) A dxdy T T! ; P = p S dp dxdy T T! T! (Γ Γ! ) ( ) = A ( ) The first term can be expanded in the space domain: T = T + T *, + T T = T * we can decompose P! into mean and eddy components; the mean is p S P M = T! (Γ dp d Γ! ) dxdy T T!. p S (You may verify that P M = dxdy dp T' A T! ( Γ ).) d Γ! d We have included the T! term in the mean available potential enery because it has no zonal variations. The eddy component is the remainin part, P E = p S dp T! (Γ d Γ! ) T * dydx. Now let us compare P to TPE: 7
8 p S dp dxdyt' T! (Γ d Γ! ) T' ~ p S C!! P TT(Γ d Γ) C P dp dxdyt T ~ T;! = Γd ; T' ~ ΔT; C P P (ΔT) ~ ; ΔT ~ 3 ο, T ~ 5 ο ; Γ! ~ Γ TPE d T Γ! Γ d P 3 ~ ; Hence only ~ to % of the TPE is available for conversion to TPE 5 ~ 7 kinetic enery. Balance Equations for Kinetic Enery: We have already derived this for the eddy kinetic enery, interated over the whole atmosphere, (remember, K E = u + v ), i.e.!!!! K t E dydp = ω ω Z Z dydp + p p ( v H F v H F )dydp + u uv + uω dydp + v v + vω dydp. y p y p Let us expand these terms in the space domain, u = u + u*, etc. Consider the first term: expandin, and also substitutin for H.E., Z α = = ; we have ( ω ωα α )dydp = ω * α * dydp. p ρ The term C ( P E,K E ) ω * α * dydp represents a conversion of eddy potential enery to eddy kinetic enery, i.e. in pressure coordinates in zonal planes, if warm air is risin (α* >,ω* < ) and cold air is sinkin α* <,ω* > ), then C P,K ) >, etc. The factor is necessary to make C have ( ( E E units of enery per unit time. Also, recall that this is total potential enery, i.e. internal enery plus conventional potential enery, not just the latter. 8
9 Similarly, the second term when expanded becomes simply u * F *+v * F * dydp = D K x y ( E ),! where D ( K ) represents the dissipation of eddy kinetic enery due to friction. Since v E! and F are neatively correlated, the dissipation, D(K E ) is positive. Now consider the third term, and expand it: u uv + uω dydp y p = u ( u v + u * v * ) + ( u ω + u * ω * ) dydp. y p The first and third terms can be combined and rewritten: expandin the derivatives, we have: v ω u + u u v + + u u u ω y y p p = v u + ω u = v! u, y p because of continuity. Now we have a pure diverence, and can interate by Gauss! theorem, and applyin the B.C. v nˆ =, we et no contribution from these terms. Only the eddy terms survive. The same thin happens with the fourth term in the equation, and these terms are combined to define another conversion, this time a conversion from eddy kinetic enery to mean kinetic enery: C ( K,K ) = u u * v * + u u * ω * + v v * + v v * ω * dydp E M y p y p Note that in the Q-G approximation, the vertical flux terms would disappear and the v term would be small, of order Ro. Thus when there is converence of eddy momentum flux positively correlated with the momentum itself, either zonal or meridional, then C K,K ) >. (Peixoto and Oort ( E M (99) ive the conversions in the more complicated spherical coordinates.) Thus we can now write for our eddy kinetic enery equation: t K E dydp = C ( P E,K E ) C ( K E,K M ) D ( K E ) Also, we now need merely to substitute this into our earlier result for the sum of K E and u + K M to et an equation for K M = v : 9
10 K M dydp = ω Z dydp + v! F! dydp H K E dydp; t p t Rewritin the ω term as before, and substitutin for C ( P E, K E ) and D ( K ) we et E!!!! K v H F v t ωα dydp + ω * α * dydp + dydp M dydp = H * F * dydp + C K E,K M!! = F ω α dydp + v H dydp + C ( K E,K M ). ( ) Define C P M,K M ( ) = ω α dydp; i.e. now warm air risin and cold air sinkin in meridional planes converts zonal mean total potential enery to zonal mean kinetic enery; and define D K ) =!! ( v H F dydp M = dissipation of the zonal mean flow, and we have t K M dydp = C ( K E,K M ) + C ( P M,K M ) D ( K M ). Balance equation for mean available potential enery: We need an equation for P M = P S dxdy dp T' ; T' = T T;! θ' = θ θ T! (Γ Γ! )!, d T! = T! ( ), etc. We start by writin the heat balance equation in p coordinates with θ as the dependent variable: dt dp = H;! dp p R ω = ; θ = T CP ; = dt ρ dt dt p C P T ω p θ p ω = ω RT p ρc + p θ P p p p pc P p T Substitute and divide throuh by C = C P P θ θ θ θ θ p + ω H! + u + v = = θ + v " θ. t x y p p C P t p
11 Substitute θ = θ' + θ! ( p); also H! = H! ' + H " ( ) " + "! (θ'+ θ! ) + v H (θ'+ θ! ) + (ωθ') + (ωθ! ) = p ( H ' H ). t p p p C P Because of continuity, v H θ! + (ωθ! ) = ω θ!. Substitutin this into the equation and p p! nˆ =, we have take the x,y averae, ( ~ ) ; usin the B.C., v θ + + " θ ω! p H = ; ω! = because of continuity. t p p p C P Note that this is precisely the equation on which lobal mean radiative-convective models are based (cf. Hantel, 976). Also, this equation in effect tells us that in equilibrium, the mean diabatic heatin is balanced by the mean eddy vertical heat flux, ωθ'!. Now subtract this equation from the θ'+ θ! equation:! ( ωθ' ) +! K #! θ' + v! H θ'+ " θ # p H $ ' ( ωθ' ωθ' ) + ω = t p p p C P Take the zonal mean of this equation, and multiply by G θ' where p G = p T! (Γ d Γ! ) ; G t θ' + G θ' y vθ' + G θ " G p θ' = ωθ' ωθ' + G θ' ω p C P p I II III IV V VI (! ) p H # ' θ' Now consider what will happens to each of these terms when interated over y,p. (The interation over x is implicit.) In I, in order to obtain the desired term, P M t, we have to take G inside the time derivative, i.e., we have to nelect the time variations in T! and Γ! in this term. Since these are comparable in manitude to the already nelected y variations, it is no worse an assumption. If we are doin a lobal analysis, it is a much better approximation because the seasonal chanes in the two hemispheres counteract each other. (See Oort & Peixoto, 983, Fi. 47). Substitutin for G and θ' the I term becomes: t p p!t Γ d! Γ ( ) p p T' dydp = P M. t Also note that this approximation does not affect the calculation of the mean balance from this equation.
12 Terms II and III can be combined and simplified by interatin the perfect derivatives, as follows: dydpg θ' vθ' + ωθ' y p ' (G θ' ) y p y p = dydp vθ' G θ' + G θ' θ'ω vθ' (G θ' ) ωθ = dydp ( v θ' + v * θ'* ) The first and third terms can be rewritten as ( G θ' y ) + ω θ' + ω * θ'* ( G θ' p ) v G θ' + ω G θ' = v! G θ', y p because of continuity, and this interates to zero. Notin that θ* = θ'*, we now have II + III = dydp v * θ * (G θ' ) + ω * θ * (G θ' ) and substitutin for θ', θ * y p and G, this becomes T' ω * T * p T' dydp v * T * y T! (Γ Γ! d ) + p p T! (Γ Γ = C P M,P E d ) ( ). This represents the conversion from mean to eddy available potential enery. Thus if there are neative correlations between the eddy transports and the radients of T', (actually in the ω * T * case it is not the radient of T', but the coefficient of T' does not vary stronly in p), i.e. if the transports are down radient, then C P,P ) > ( M E and P M is converted to P E (see diaram below). In effect, lonitudinal T radients are bein created from the meridional and vertical radients. Note that in the quasi eostrophic approximation, the vertical eddy flux term would not appear.
13 Fiure by MIT OCW. ( ) ωθ'! Next consider the term IV: dydpg θ' p =, because ωθ! ' and G are independent of y, θ! ' = by definition, and θ ' interates to zero over y. Now consider the term V:! =. We find p p Γ d p V = dydpg θ' ω θ! ; substitute for G, θ' and θ! (Γ Γ d ) θ! V = dydp ω α' = +C ( P M,K M ) Note that α' and α are interchaneable in the interand because the interal of ω is. The factor of comes in because we are interatin in pressure coordinates. Note that α = = RT ; α' = RT' ; if warm air rises and cold air sinks, ω and α' ρ p p in meridional planes are neatively correlated, and P M is converted to K M. Also note that we can replace α' by α because ω α! interates to zero. Finally, the VI term, upon substitutin for G and θ', becomes: VI = G K p H! ' θ' dydp C P p 3
14 H! ' = T' dydp G P M CP T" ( Γ Γ " d ) ( ) mean available potential enery is enerated if the zonal mean T and H! fluctuation are correlated, i.e. if we have heatin where T is hih and coolin where T is low. Note that # T' can aain be replaced by T because H! ' T" =. Now we put it all toether to et the P M equation: P M = G ( P M ) C ( P M, P E ) C ( P M, K M ). t Finally, a P E equation can be derived in similar fashion by derivin an equation for T * from the same thermodynamic equation. The already defined conversions for C ( P E,K E ) and C( P M,P E ) will appear aain, and the only new term is a eneration term ( M exactly analoous to G P ), i.e.! H * T * ( ) G P = dydp E " CP T) (Γ Γ ) d Thus our equation is P E = G ( P E ) + C ( P M, P E ) C( P E,K E ). t Observational Analyses: We shall follow Oort & Peixoto (983). Oort & Rasmussen (99) cite Oort and Peixoto (983) as their source, but there are minor discrepancies. We will only look at the Northern Hemisphere analyses, because of the better data base. Note however that if we do, then in eneral there are boundary fluxes across the equator which must be added to the equations above. Oort & Peixoto (974) calculated these and found that they were enerally not important except for the flux of P M in the seasonal ( M extremes, B P ), associated with the stron seasonal winter Hadley Cells and the term v G θ'. Thus we add B P ) to the P M equation, P M ( = B ( P M ) +... M The analyses are based on the years in the MIT eneral circulation library, and they calculated the four eneries, C ( P M,K M ), C( P M,P E ), and C( K,K ). Note that in E M evaluatin the first two conversions they assumed Ro <<, and only retained the leadin terms. G ( P M ) and G ( P E ) cannot be calculated because of lack of data on diabatic t 4
15 ( E E heatin, C P,K ) cannot be calculated because of the lack of data on ω * α *, and D( K ) and D K ) because of lack of data on friction. However if any one of the three M ( E G ( P ), C P,K ) and D K ) were known, all the other unknowns can be calculated as E ( E E ( E residuals, assumin a stationary state (see Fi., top). Oort & Peixoto (983) assumed values for G ( P E ) and calculated the others for the annual mean and for the December- January-February and June-July-Auust seasons. They did not explain the basis for their assumed G ( P E ) values in any of their papers. We can reard them as educated uesses. They assumed positive values, which could be explained as bein due to a positive correlation between H! * and T! * (i.e. in zonal planes) due to baroclinic eddies releasin latent heat (which is contained in H! in the enery equations), i.e. warm air movin poleward rises, leadin to condensational heatin. 5
16 Fiure 6
17 Fiure 7
18 Fi. (bottom) shows the annual mean Northern Hemisphere enery cycle, and Fi. shows the Northern Hemisphere seasonal cycles. (Note apparent small errors in the residual calculations in the annual mean.) Lookin first at the annual mean, we note:. The main input is due to G P ), caused by heatin in warmer low latitudes and ( M coolin in cooler hih latitudes. The radiative fluxes cause this both directly and indirectly by drivin low latitude moist convection and condensation.. From P M enery flows mainly to P E, because of baroclinic eddies i.e., their downradient heat transport. 3. P E enery flows mainly to K E, i.e. K E is enerated by upward eddy heat fluxes, i.e. ω * T * >. 4. This flow is mainly balanced by dissipation, i.e., mainly by friction near the surface. 5. There is however some flow of enery from K E to K M, i.e., there is an upradient transport of momentum by the eddies, on balance. In fact this is the main source of K M. 6. There is some flow of enery from K M to P M, i.e., on balance ω T <, i.e. on balance this conversion is dominated by the indirect Ferrel Cell. 7. D K ( ) are relatively small. ( E ) and G P E Fi. shows the seasonal chanes. Note:. In winter the dominant enery source is B ( P M ), i.e. P M is mainly bein enerated in the Southern Hemisphere, but a lot of it is bein transported across the equator to the Northern Hemisphere.. G ( P M ) is reater in June-July-Auust, presumably because tropical moist convection and the associated LH release is reater then. 3. The eddy converences are all much weaker in summer and thus the flow of enery is weaker in summer, and all the enery forms are less. 4. The Ferrel Cell does not dominate C ( P, K ) in December-January- M M February, because it has to contend with the very stron winter Hadley cell. Because of the difficulty in calculatin the enery cycle (calculatin G repeatedly for example) there have been very few attempts to calculate it from models. The only one I know of was with an early version of the GISS GCM (Somerville et al, 974). It used 4 lat (lon ( ) x5 ) resolution with 9 vertical levels, and a sinle January was simulated. No boundary fluxes were calculated; the Northern Hemisphere simulated enery cycle is shown in Fi. a in Somerville et al (974). Comparin with the Peixoto and Oort winter results, we note. The enery cycles are qualitatively similar.. Most of the model s eneries are too hih (except for K E ) 8
19 3. Most of the conversions, eneration and dissipation terms are too weak, except for C P M,P E ( ) ( ) is.6, close to Oort and Peixoto s assumed value. 4. Gratifyinly, G P E Other features for the GISS simulation indicate that the differences are not due to the resolution or short time base. The lare K M may be in part because the tropospheric jet is not closed. Note that the results for the Enery Cycle tell us that the EP theorem is not at all satisfied in the trosposhere, because there is a net exchane of enery between the eddies and the mean flow. C P,P ) is about four times C K,K ( M E ( E M ) 9
2.3. PBL Equations for Mean Flow and Their Applications
.3. PBL Equations for Mean Flow and Their Applications Read Holton Section 5.3!.3.1. The PBL Momentum Equations We have derived the Reynolds averaed equations in the previous section, and they describe
More informationCourse , General Circulation of the Earth's Atmosphere Prof. Peter Stone Section 4: Water Vapor Budget
Course 12.812, General Circulation of the Earth's Atmosphere Prof. Peter Stone Section 4: Water Vapor Budget Water Vapor Distribution First let us look at the distribution of specific humidity, q. The
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 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(! g. + f ) Dt. $% (" g. + f ) + f
The QG Vorticity equation The complete derivation of the QG vorticity equation can be found in Chapter 6.2.. Please read this section and keep in mind that extra approximations have been made besides those
More informationDevelopment Theory. Chapter 10
Chapter 1 Development Theory Development Theory In this section I will discuss: - further aspects of the structure and dynamics of synoptic-scale disturbances, and - derive Sutcliffe's development theory,
More informationENERGETICS. [This article treats a technical aspect of climate and weather studies; some of it is intended for readers at an advanced level.
ENERGETICS. [This article treats a technical aspect of climate and weather studies; some of it is intended for readers at an advanced level.] Atmospheric energetics, the study of the distribution and transformation
More informationCHAPTER 4. THE HADLEY CIRCULATION 59 smaller than that in midlatitudes. This is illustrated in Fig. 4.2 which shows the departures from zonal symmetry
Chapter 4 THE HADLEY CIRCULATION The early work on the mean meridional circulation of the tropics was motivated by observations of the trade winds. Halley (1686) and Hadley (1735) concluded that the trade
More informationF = ma. ATS 150 Global Climate Change Winds and Weather. Scott Denning CSU CMMAP 1. Please read Chapter 6 from Archer Textbook
Winds and Weather Please read Chapter 6 from Archer Textbook Circulation of the atmosphere and oceans are driven by energy imbalances Energy Imbalances What Makes the Wind Blow? Three real forces (gravity,
More informationToday s Lecture (Lecture 5): General circulation of the atmosphere
Climate Dynamics (Summer Semester 2017) J. Mülmenstädt Today s Lecture (Lecture 5): General circulation of the atmosphere Reference Hartmann, Global Physical Climatology (1994), Ch. 2, 3, 6 Peixoto and
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 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 information5 Shallow water Q-G theory.
5 Shallow water Q-G theory. So far we have discussed the fact that lare scale motions in the extra-tropical atmosphere are close to eostrophic balance i.e. the Rossby number is small. We have examined
More informationSimple Climate Models
Simple Climate Models Lecture 3 One-dimensional (vertical) radiative-convective models Vertical atmospheric processes The vertical is the second important dimension, because there are... Stron radients
More informationGoverning Equations and Scaling in the Tropics
Governing Equations and Scaling in the Tropics M 1 ( ) e R ε er Tropical v Midlatitude Meteorology Why is the general circulation and synoptic weather systems in the tropics different to the those in the
More informationIsentropic flows and monsoonal circulations
Isentropic flows and monsoonal circulations Olivier Pauluis (NYU) Monsoons- Past, Present and future May 20th, 2015 Caltech, Pasadena Outline Introduction Global monsoon in isentropic coordinates Dry ventilation
More informationAorangi (Mt. Cook) we typically view the state of the troposphere on "isobaric surfaces" (surfaces of constant pressure) rather
Ch3 (ctd.) EAS270_Ch3_BehaviourAtmos_B.odp JDW, EAS UAlberta, last mod. 19 Sept. 2016 Atmospheric "behaviour" Aorani (Mt. Cook) "ride" we typically view the state of the troposphere on "isobaric surfaces"
More informationCourse , General Circulation of the Earth's Atmosphere Prof. Peter Stone Section 9: Numerical Models
Course 1.81, General Circulation of the Earth's Atmosphere Prof. Peter Stone Section 9: Numerical Models Phillips Numerical Experiment: We are going to discuss Phillips (1956) original numerical experiment
More informationAnnular Mode Variability of the Atmospheric Meridional Energy Transport and Circulation
27 J O U R N A L O F T H E A T M O S P H E R I C S C I E N C E S VOLUME 72 Annular Mode Variability of the Atmospheric Meridional Enery Transport and Circulation RAY YAMADA Courant Institute of Mathematical
More informationGeopotential tendency and vertical motion
Geopotential tendency and vertical motion Recall PV inversion Knowin the PV, we can estimate everythin else! (Temperature, wind, eopotential ) In QG, since the flow is eostrophic, we can obtain the wind
More informationTropical Meridional Circulations: The Hadley Cell
Tropical Meridional Circulations: The Hadley Cell Introduction Throughout much of the previous sections, we have alluded to but not fully described the mean meridional overturning circulation of the tropics
More informationConservation of absolute vorticity. MET 171A: Barotropic Midlatitude Waves. Draw a picture of planetary vorticity
Conservation of absolute vorticity : Barotropic Midlatitude Waves Recall the important terms in the vorticity equation in the middle troposphere, near the level of non-diverence Lecture Outline 1. Conservation
More informationONLINE: MATHEMATICS EXTENSION 2 Topic 6 MECHANICS 6.3 HARMONIC MOTION
ONINE: MATHEMATICS EXTENSION Topic 6 MECHANICS 6.3 HARMONIC MOTION Vibrations or oscillations are motions that repeated more or less reularly in time. The topic is very broad and diverse and covers phenomena
More informationLecture 11: Meridonal structure of the atmosphere
Lecture 11: Meridonal structure of the atmosphere September 28, 2003 1 Meridional structure of the atmosphere In previous lectures we have focussed on the vertical structure of the atmosphere. Today, we
More informationMotion in Two Dimensions Sections Covered in the Text: Chapters 6 & 7, except 7.5 & 7.6
Motion in Two Dimensions Sections Covered in the Tet: Chapters 6 & 7, ecept 7.5 & 7.6 It is time to etend the definitions we developed in Note 03 to describe motion in 2D space. In doin so we shall find
More informationTransient and Eddy. Transient/Eddy Flux. Flux Components. Lecture 3: Weather/Disturbance. Transient: deviations from time mean Time Mean
Lecture 3: Weather/Disturbance Transients and Eddies Climate Roles Mid-Latitude Cyclones Tropical Hurricanes Mid-Ocean Eddies Transient and Eddy Transient: deviations from time mean Time Mean Eddy: deviations
More information1. The vertical structure of the atmosphere. Temperature profile.
Lecture 4. The structure of the atmosphere. Air in motion. Objectives: 1. The vertical structure of the atmosphere. Temperature profile. 2. Temperature in the lower atmosphere: dry adiabatic lapse rate.
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 informationATMO 551a Fall The Carnot Cycle
What is a arnot ycle and Why do we care The arnot ycle arnot was a French engineer who was trying to understand how to extract usable mechanical work from a heat engine, that is an engine where a gas or
More informationBerichte zur Erdsystemforschung. Daniel Hernández Deckers. Reports on Earth System Science
Warming up the atmosphere s heat engine: atmospheric energetics with higher greenhouse gas concentrations Daniel Hernández Deckers Berichte zur Erdsystemforschung 85 2011 Reports on Earth System Science
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 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 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 informationDynamics of the Zonal-Mean, Time-Mean Tropical Circulation
Dynamics of the Zonal-Mean, Time-Mean Tropical Circulation First consider a hypothetical planet like Earth, but with no continents and no seasons and for which the only friction acting on the atmosphere
More informationHydrodynamic conservation laws and turbulent friction in atmospheric circulation models
Hydrodynamic conservation laws and turbulent friction in atmospheric circulation models Erich Becker Leibniz-Institute of Atmospheric Physics, Kühlungsborn, Germany Including contributions from Ulrike
More informationThe broadest definition of a low-level jet (LLJ) is simply any lower-tropospheric maximum in the vertical profile of the horizontal winds.
2.4.3. Low-level (especially nocturnal) Jet The broadest definition of a low-level jet (LLJ) is simply any lower-tropospheric maximum in the vertical profile of the horizontal winds. A LLJ can occur under
More informationQuasi-geostrophic motion
Capter 8 Quasi-eostropic motion Scale analysis for synoptic-scale motions Simplification of te basic equations can be obtained for synoptic scale motions. Consider te Boussinesq system ρ is assumed to
More informationConvergence of DFT eigenvalues with cell volume and vacuum level
Converence of DFT eienvalues with cell volume and vacuum level Sohrab Ismail-Beii October 4, 2013 Computin work functions or absolute DFT eienvalues (e.. ionization potentials) requires some care. Obviously,
More informationThe effect of varying forcing on the transport of heat by transient eddies.
The effect of varying forcing on the transport of heat by transient eddies. LINDA MUDONI Department of Atmospheric Sciences, Iowa State University May 02 20 1. Introduction The major transport of heat
More information11 Free vibrations: one degree of freedom
11 Free vibrations: one deree of freedom 11.1 A uniform riid disk of radius r and mass m rolls without slippin inside a circular track of radius R, as shown in the fiure. The centroidal moment of inertia
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 information( δx ) = δ dx = δu, etc. Course , General Circulation of the Earth's Atmosphere Prof. Peter Stone Section 2: Analysis Techniques
Course 12.812, General Circulation of the Earth's Atmosphere Prof. Peter Stone Section 2: Analysis Techniques Pressure Coordinates: Most meteorological measurements are made in pressure coordinates, e.g.,
More informationInterhemispheric climate connections: What can the atmosphere do?
Interhemispheric climate connections: What can the atmosphere do? Raymond T. Pierrehumbert The University of Chicago 1 Uncertain feedbacks plague estimates of climate sensitivity 2 Water Vapor Models agree
More informationNOTES AND CORRESPONDENCE. On the Seasonality of the Hadley Cell
1522 JOURNAL OF THE ATMOSPHERIC SCIENCES VOLUME 60 NOTES AND CORRESPONDENCE On the Seasonality of the Hadley Cell IOANA M. DIMA AND JOHN M. WALLACE Department of Atmospheric Sciences, University of Washington,
More informationFluid Mechanics (ME 201) 1
Fluid Mechanics (ME 201) 1 9 Forces on Submered Bodies 9.1 Net Force on a Submered Plane Surfaces In this section, we will compute the forces actin on a plane surface. We consider a surface oriented at
More informationChapter 2. Turbulence and the Planetary Boundary Layer
Chapter 2. Turbulence and the Planetary Boundary Layer In the chapter we will first have a qualitative overview of the PBL then learn the concept of Reynolds averain and derive the Reynolds averaed equations.
More informationThe Circulation of the Atmosphere:
The Circulation of the Atmosphere: Laboratory Experiments (see next slide) Fluid held in an annular container is at rest and is subjected to a temperature gradient. The less dense fluid near the warm wall
More informationDynamics of the Atmosphere. General circulation of the atmosphere
12.810 Dynamics of the Atmosphere General circulation of the atmosphere 1 Spinup of the general circulation in an idealized model Fig. 1 Schneider, General circulation of the atmosphere, 2006 2 Sigma 0.2
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 informationESCI 344 Tropical Meteorology Lesson 12 Tropical Cyclones: Motion and Analysis
ESCI 344 Tropical Meteoroloy Lesson 12 Tropical Cyclones: Motion and Analysis References: A Global View of Tropical Cyclones, Elsberry (ed.) The Hurricane, Pielke Forecasters Guide to Tropical Meteoroloy,
More informationX Y. Equator. Question 1
Question 1 The schematic below shows the Hadley circulation in the Earth s tropical atmosphere around the spring equinox. Air rises from the equator and moves poleward in both hemispheres before descending
More information2.2 Differentiation and Integration of Vector-Valued Functions
.. DIFFERENTIATION AND INTEGRATION OF VECTOR-VALUED FUNCTIONS133. Differentiation and Interation of Vector-Valued Functions Simply put, we differentiate and interate vector functions by differentiatin
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 informationExam 2A Solution. 1. A baseball is thrown vertically upward and feels no air resistance. As it is rising
Exam 2A Solution 1. A baseball is thrown vertically upward and feels no air resistance. As it is risin Solution: Possible answers: A) both its momentum and its mechanical enery are conserved - incorrect.
More informationGeneral Circulation. Nili Harnik DEES, Lamont-Doherty Earth Observatory
General Circulation Nili Harnik DEES, Lamont-Doherty Earth Observatory nili@ldeo.columbia.edu Latitudinal Radiation Imbalance The annual mean, averaged around latitude circles, of the balance between the
More information1 Introduction to Governing Equations 2 1a Methodology... 2
Contents 1 Introduction to Governing Equations 2 1a Methodology............................ 2 2 Equation of State 2 2a Mean and Turbulent Parts...................... 3 2b Reynolds Averaging.........................
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 informationAtmospheric Dynamics: lecture 2
Atmospheric Dynamics: lecture 2 Topics Some aspects of advection and the Coriolis-effect (1.7) Composition of the atmosphere (figure 1.6) Equation of state (1.8&1.9) Water vapour in the atmosphere (1.10)
More informationThe Planetary Circulation System
12 The Planetary Circulation System Learning Goals After studying this chapter, students should be able to: 1. describe and account for the global patterns of pressure, wind patterns and ocean currents
More informationLecture 5: Atmospheric General Circulation and Climate
Lecture 5: Atmospheric General Circulation and Climate Geostrophic balance Zonal-mean circulation Transients and eddies Meridional energy transport Moist static energy Angular momentum balance Atmosphere
More informationHomework # 2. SOLUTION - We start writing Newton s second law for x and y components: F x = 0, (1) F y = mg (2) x (t) = 0 v x (t) = v 0x (3)
Physics 411 Homework # Due:..18 Mechanics I 1. A projectile is fired from the oriin of a coordinate system, in the x-y plane (x is the horizontal displacement; y, the vertical with initial velocity v =
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 information196 7 atmospheric oscillations:
196 7 atmospheric oscillations: 7.4 INTERNAL GRAVITY (BUOYANCY) WAVES We now consider the nature of gravity wave propagation in the atmosphere. Atmospheric gravity waves can only exist when the atmosphere
More informationANALYZE In all three cases (a) (c), the reading on the scale is. w = mg = (11.0 kg) (9.8 m/s 2 ) = 108 N.
Chapter 5 1. We are only concerned with horizontal forces in this problem (ravity plays no direct role). We take East as the +x direction and North as +y. This calculation is efficiently implemented on
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 informationPhysicsAndMathsTutor.com
. A particle of mass m is projected vertically upwards, at time t =, with speed. The particle is mv subject to air resistance of manitude, where v is the speed of the particle at time t and is a positive
More informationWater Vapor and the Dynamics of Climate Changes
Water Vapor and the Dynamics of Climate Changes Tapio Schneider California Institute of Technology (based on Rev. Geophys. article with Xavier Levine and Paul O Gorman) Water vapor dynamics in warming
More informationWire antenna model of the vertical grounding electrode
Boundary Elements and Other Mesh Reduction Methods XXXV 13 Wire antenna model of the vertical roundin electrode D. Poljak & S. Sesnic University of Split, FESB, Split, Croatia Abstract A straiht wire antenna
More informationExperiment 3 The Simple Pendulum
PHY191 Fall003 Experiment 3: The Simple Pendulum 10/7/004 Pae 1 Suested Readin for this lab Experiment 3 The Simple Pendulum Read Taylor chapter 5. (You can skip section 5.6.IV if you aren't comfortable
More informationClimate of an Earth- like Aquaplanet: the high- obliquity case and the <dally- locked case
Climate of an Earth- like Aquaplanet: the high- obliquity case and the
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 informationC. Non-linear Difference and Differential Equations: Linearization and Phase Diagram Technique
C. Non-linear Difference and Differential Equations: Linearization and Phase Diaram Technique So far we have discussed methods of solvin linear difference and differential equations. Let us now discuss
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 informationwarmest (coldest) temperatures at summer heat dispersed upward by vertical motion Prof. Jin-Yi Yu ESS200A heated by solar radiation at the base
Pole Eq Lecture 3: ATMOSPHERE (Outline) JS JP Hadley Cell Ferrel Cell Polar Cell (driven by eddies) L H L H Basic Structures and Dynamics General Circulation in the Troposphere General Circulation in the
More informationGeophysical limits to global wind power
SUPPLEMENTARY INFORMATION DOI: 1.138/NCLIMATE1683 Geophysical limits to global wind power Kate Marvel Program for Climate Model Diagnosis and Intercomparison Lawrence Livermore National Laboratory P.O.
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 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 information[ ][ ] mg = R GM g = R GM = g R 2. Definition of G: From (1) we have, 2 NEWTON S LAW OF GRAVITATION:
NEWTON S LAW OF GAVITATION: Statement: Every particle in the universe attracts every other particle with a force which is directly proportional to the product of the masses and inversely proportional to
More informationProblems of the 9 th International Physics Olympiads (Budapest, Hungary, 1976)
Problems of the 9 th International Physics Olympiads (Budapest, Hunary, 1976) Theoretical problems Problem 1 A hollow sphere of radius R = 0.5 m rotates about a vertical axis throuh its centre with an
More informationPhysics 111 P 2 A = P 1. A + mg = P 1. A + ρ( AΔh)g. Wednesday, 8-9 pm in NSC 118/119 Sunday, 6:30-8 pm in CCLIR 468.
ics Announcements day, ember 11, 011 C5: Fluids Pascal s Principle Archimede s Principle Fluid Flows Continuity Equation Bernoulli s Equation Toricelli s Theorem Wednesday, 8-9 pm in NSC 118/119 Sunday,
More informationUncertainties in (energy) budgets
Uncertainties in (enery) budets Thanks to Michael Mayer, John M. Edwards, Patrick Hyder, Andrea Storto, Marianne Pietschni, Sebastian Stichelberer, Eric Boisseson, Patrick Laloyaux, Elke Rustemeier, Markus
More informationCode_Aster. Element CABLE_GAINE
Titre : Élément CABLE_GAINE Date : 28/07/2015 Pae : 1/12 Element CABLE_GAINE Summary: The element CABLE_GAINE presented in this document is to model cables of prestressin bein able to rub or slip into
More informationTransient/Eddy Flux. Transient and Eddy. Flux Components. Lecture 7: Disturbance (Outline) Why transients/eddies matter to zonal and time means?
Lecture 7: Disturbance (Outline) Transients and Eddies Climate Roles Mid-Latitude Cyclones Tropical Hurricanes Mid-Ocean Eddies (From Weather & Climate) Flux Components (1) (2) (3) Three components contribute
More informationTest Review # 7. Combined Gas Law PV T PV T. Ideal Gas Law PV = nrt. Chemistry H-3: Form TR7.6A
Chemistry H-3: Form TR7.6A TEST 9 REVIEW Name Date Period Test Review # 7 ENERGY Calculatin Joules. When you heat a solid, it s temperature enerally oes up. There is a relationship between heat and temperature,
More informationAtmospheric Circulation
Atmospheric Circulation (WAPE: General Circulation of the Atmosphere and Variability) François Lott, flott@lmd.ens.fr http://web.lmd.jussieu.fr/~flott 1) Mean climatologies and equations of motion a)thermal,
More informationno eddies eddies Figure 3. Simulated surface winds. Surface winds no eddies u, v m/s φ0 =12 φ0 =0
References Held, Isaac M., and Hou, A. Y., 1980: Nonlinear axially symmetric circulations in a nearly inviscid atmosphere. J. Atmos. Sci. 37, 515-533. Held, Isaac M., and Suarez, M. J., 1994: A proposal
More information(a) 1m s -2 (b) 2 m s -2 (c) zero (d) -1 m s -2
11 th Physics - Unit 2 Kinematics Solutions for the Textbook Problems One Marks 1. Which one of the followin Cartesian coordinate system is not followed in physics? 5. If a particle has neative velocity
More informationThe dynamics of high and low pressure systems
The dynamics of high and low pressure systems Newton s second law for a parcel of air in an inertial coordinate system (a coordinate system in which the coordinate axes do not change direction and are
More informationMonday 7 October 2013, Class #15
Monday 7 October 2013, Class #15 Concepts for Today (Basics for Thermodynamics) Weather versus climate Lapse Rate (Adiabatic Lapse Rate) Ideal Gas Law Adiabatic Processes Potential Temperature Hydrostatic
More informationEstimating the zonal wavenumber dependence of the meridional energy transport
MO91 egree Project in Atmospheric Sciences, Oceanography and Climate, 3hp Estimating the zonal wavenumber dependence of the meridional energy transport Mattias Burtu Supervisor: Rune Grand Graversen epartment
More informationThe General Circulation of the Atmosphere
Annu. Rev. Earth Planet. Sci. 2006. 34:655 88 The Annual Review of Earth and Planetary Science is online at earth.annualreviews.org doi: 10.1146/ annurev.earth.34.031405.125144 Copyright c 2006 by Annual
More informationPHYSICS 218 SOLUTION TO HW 8. Created: November 20, :15 pm Last updated: November 21, 2004
Created: November 20, 2004 7:5 pm Last updated: November 2, 2004. Schroeder.6 (a) The three forces actin on the slab of thickness dz are ravity and the pressure from above and below. To achieve equilibrium
More informationREVIEW: Going from ONE to TWO Dimensions with Kinematics. Review of one dimension, constant acceleration kinematics. v x (t) = v x0 + a x t
Lecture 5: Projectile motion, uniform circular motion 1 REVIEW: Goin from ONE to TWO Dimensions with Kinematics In Lecture 2, we studied the motion of a particle in just one dimension. The concepts of
More informationA Statistical Generalization of the Transformed Eulerian-Mean Circulation for an Arbitrary Vertical Coordinate System
766 J O U R N A L O F T H E A T M O S P H E R I C S C I E N C E S VOLUME 68 A Statistical Generalization of the Transformed Eulerian-Mean Circulation for an Arbitrary Vertical Coordinate System OLIVIER
More informationThe Global Atmospheric Circulation in Moist Isentropic Coordinates
1JUNE 2010 P A U L U I S E T A L. 3077 The Global Atmospheric Circulation in Moist Isentropic Coordinates OLIVIER PAULUIS Courant Institute of Mathematical Sciences New York University New York New York
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 informationProblem Set 5 Solutions
Problem Set 5 Solutions Dorian Abbot APM47 0/8/04. (a) We consider π θ π. The pendulum pointin down corresponds to θ=0 and the pendulum pointin up corresponds to θ π. Define ν θ. The system can be rewritten:
More informationGrand Challenges in Global Circulation Dynamics
Grand Challenges in Global Circulation Dynamics Tapio Schneider ETH Zurich, Caltech (Source: CLAUS, http://badc.nerc.ac.uk/data/claus/) Grand Challenges in Global Circulation Dynamics Tapio Schneider ETH
More informationMODELING A METAL HYDRIDE HYDROGEN STORAGE SYSTEM. Apurba Sakti EGEE 520, Mathematical Modeling of EGEE systems Spring 2007
MODELING A METAL HYDRIDE HYDROGEN STORAGE SYSTEM Apurba Sakti EGEE 520, Mathematical Modelin of EGEE systems Sprin 2007 Table of Contents Abstract Introduction Governin Equations Enery balance Momentum
More information