The QG Vorticity Equation The QG Vorticity Equation The - - PowerPoint PPT Presentation

the qg vorticity equation the qg vorticity equation
SMART_READER_LITE
LIVE PREVIEW

The QG Vorticity Equation The QG Vorticity Equation The - - PowerPoint PPT Presentation

The QG Vorticity Equation The QG Vorticity Equation The quasi-geostrophic vorticity is g = k V g = 1 2 f 0 This enables g to be computed immediately once the geopo- tential is known. The QG Vorticity Equation The


slide-1
SLIDE 1

The QG Vorticity Equation

slide-2
SLIDE 2

The QG Vorticity Equation

The quasi-geostrophic vorticity is ζg = k · ∇ × Vg = 1 f0 ∇2Φ This enables ζg to be computed immediately once the geopo- tential is known.

slide-3
SLIDE 3

The QG Vorticity Equation

The quasi-geostrophic vorticity is ζg = k · ∇ × Vg = 1 f0 ∇2Φ This enables ζg to be computed immediately once the geopo- tential is known. It also means that the geopotential can be deduced from the vorticity by inverting the Laplacian operator.

slide-4
SLIDE 4

The QG Vorticity Equation

The quasi-geostrophic vorticity is ζg = k · ∇ × Vg = 1 f0 ∇2Φ This enables ζg to be computed immediately once the geopo- tential is known. It also means that the geopotential can be deduced from the vorticity by inverting the Laplacian operator. This invertibility principle holds in a much more general context and is a central tenet of the theory of balanced flows.

slide-5
SLIDE 5

The QG Vorticity Equation

The quasi-geostrophic vorticity is ζg = k · ∇ × Vg = 1 f0 ∇2Φ This enables ζg to be computed immediately once the geopo- tential is known. It also means that the geopotential can be deduced from the vorticity by inverting the Laplacian operator. This invertibility principle holds in a much more general context and is a central tenet of the theory of balanced flows. Since the Laplacian of a function tends to have a minimum where the function has a maximum, and vice-versa, Positive Vorticity is associated with Low Pressure that is, low values of the geopotential, and Negative Vorticity is associated with High Pressure

slide-6
SLIDE 6

We write the quasi-geostrophic momentum equation in com- ponent form dgug dt − f0va − βy vg = 0 dgvg dt + f0ua + βy ug = 0

2

slide-7
SLIDE 7

We write the quasi-geostrophic momentum equation in com- ponent form dgug dt − f0va − βy vg = 0 dgvg dt + f0ua + βy ug = 0 Subtracting the y-derivative of the first equation from the x-derivative of the first, we get the vorticity equation dgζg dt = −f0∇·Va − βvg

2

slide-8
SLIDE 8

We write the quasi-geostrophic momentum equation in com- ponent form dgug dt − f0va − βy vg = 0 dgvg dt + f0ua + βy ug = 0 Subtracting the y-derivative of the first equation from the x-derivative of the first, we get the vorticity equation dgζg dt = −f0∇·Va − βvg Note that a term δgζg arising from the total time derivative vanishes.

2

slide-9
SLIDE 9

We write the quasi-geostrophic momentum equation in com- ponent form dgug dt − f0va − βy vg = 0 dgvg dt + f0ua + βy ug = 0 Subtracting the y-derivative of the first equation from the x-derivative of the first, we get the vorticity equation dgζg dt = −f0∇·Va − βvg Note that a term δgζg arising from the total time derivative vanishes. Exercise: Verify the derivation of the vorticity equation. Expand d/dt and proceed as indicated above.

2

slide-10
SLIDE 10

The β-term may be written −βvg = −Vg·∇f

3

slide-11
SLIDE 11

The β-term may be written −βvg = −Vg·∇f The divergence of the ageostrophic wind may be replaced by the vertical gradient of vertical velocity using the continuity equation.

3

slide-12
SLIDE 12

The β-term may be written −βvg = −Vg·∇f The divergence of the ageostrophic wind may be replaced by the vertical gradient of vertical velocity using the continuity equation. Then ∂ζg ∂t = −Vg·∇(ζg + f) + f0 ∂ω ∂p

3

slide-13
SLIDE 13

The β-term may be written −βvg = −Vg·∇f The divergence of the ageostrophic wind may be replaced by the vertical gradient of vertical velocity using the continuity equation. Then ∂ζg ∂t = −Vg·∇(ζg + f) + f0 ∂ω ∂p This equation means that the local rate of change of geo- strophic vorticity is determined by the sum of two terms:

3

slide-14
SLIDE 14

The β-term may be written −βvg = −Vg·∇f The divergence of the ageostrophic wind may be replaced by the vertical gradient of vertical velocity using the continuity equation. Then ∂ζg ∂t = −Vg·∇(ζg + f) + f0 ∂ω ∂p This equation means that the local rate of change of geo- strophic vorticity is determined by the sum of two terms:

  • The advection of the absolute vorticity by the geostrophic

wind

3

slide-15
SLIDE 15

The β-term may be written −βvg = −Vg·∇f The divergence of the ageostrophic wind may be replaced by the vertical gradient of vertical velocity using the continuity equation. Then ∂ζg ∂t = −Vg·∇(ζg + f) + f0 ∂ω ∂p This equation means that the local rate of change of geo- strophic vorticity is determined by the sum of two terms:

  • The advection of the absolute vorticity by the geostrophic

wind

  • The stretching or shrinking of fluid columns (divergence

effect).

3

slide-16
SLIDE 16

The advection itself is the sum of two terms −Vg·∇(ζg + f) = −Vg·∇ζg − βvg the advection of relative vorticity and the advection of plan- etary vorticity respectively.

4

slide-17
SLIDE 17

The advection itself is the sum of two terms −Vg·∇(ζg + f) = −Vg·∇ζg − βvg the advection of relative vorticity and the advection of plan- etary vorticity respectively. For wave-like disturbances in the mid-latitude westerlies, these two terms tend to be of opposite sign so that they counteract each other.

4

slide-18
SLIDE 18

The advection itself is the sum of two terms −Vg·∇(ζg + f) = −Vg·∇ζg − βvg the advection of relative vorticity and the advection of plan- etary vorticity respectively. For wave-like disturbances in the mid-latitude westerlies, these two terms tend to be of opposite sign so that they counteract each other. We can estimate their relative sizes:

  • Vg·∇ζg

βvg

  • ∼ V a

L2f = Ro L/a ∼ 1 so the two terms are of comparable size.

4

slide-19
SLIDE 19

5

slide-20
SLIDE 20

Let us consider a wave-like disturbance, as shown above.

6

slide-21
SLIDE 21

Let us consider a wave-like disturbance, as shown above. In Region I — upstream of the 500 hPa trough — the geo- strophic wind is flowing from the relative vorticity minimum (at the ridge) towards the relative vorticity maximum (at the trough) so that −Vg·∇ζg < 0.

6

slide-22
SLIDE 22

Let us consider a wave-like disturbance, as shown above. In Region I — upstream of the 500 hPa trough — the geo- strophic wind is flowing from the relative vorticity minimum (at the ridge) towards the relative vorticity maximum (at the trough) so that −Vg·∇ζg < 0. However, since vg < 0 the flow is from higher planetary vorticity values to lower values, so −βvg > 0.

6

slide-23
SLIDE 23

Let us consider a wave-like disturbance, as shown above. In Region I — upstream of the 500 hPa trough — the geo- strophic wind is flowing from the relative vorticity minimum (at the ridge) towards the relative vorticity maximum (at the trough) so that −Vg·∇ζg < 0. However, since vg < 0 the flow is from higher planetary vorticity values to lower values, so −βvg > 0. Hence, in Region 1

6

slide-24
SLIDE 24

Let us consider a wave-like disturbance, as shown above. In Region I — upstream of the 500 hPa trough — the geo- strophic wind is flowing from the relative vorticity minimum (at the ridge) towards the relative vorticity maximum (at the trough) so that −Vg·∇ζg < 0. However, since vg < 0 the flow is from higher planetary vorticity values to lower values, so −βvg > 0. Hence, in Region 1

  • The advection of relative vorticity decreases ζg

6

slide-25
SLIDE 25

Let us consider a wave-like disturbance, as shown above. In Region I — upstream of the 500 hPa trough — the geo- strophic wind is flowing from the relative vorticity minimum (at the ridge) towards the relative vorticity maximum (at the trough) so that −Vg·∇ζg < 0. However, since vg < 0 the flow is from higher planetary vorticity values to lower values, so −βvg > 0. Hence, in Region 1

  • The advection of relative vorticity decreases ζg
  • The advection of planetary vorticity increases ζg

6

slide-26
SLIDE 26

7

slide-27
SLIDE 27

Similar arguments, with signs reversed, apply in Region II.

8

slide-28
SLIDE 28

Similar arguments, with signs reversed, apply in Region II. Therefore, the relative vorticity advection tends to move the vorticity pattern, and hence the troughs and ridges, downstream or eastward.

8

slide-29
SLIDE 29

Similar arguments, with signs reversed, apply in Region II. Therefore, the relative vorticity advection tends to move the vorticity pattern, and hence the troughs and ridges, downstream or eastward. On the other hand, the planetary vorticity advection tends to move the vorticity pattern upstream or westward, causing the wave to regress.

8

slide-30
SLIDE 30

Similar arguments, with signs reversed, apply in Region II. Therefore, the relative vorticity advection tends to move the vorticity pattern, and hence the troughs and ridges, downstream or eastward. On the other hand, the planetary vorticity advection tends to move the vorticity pattern upstream or westward, causing the wave to regress. Since the two terms are of comparable magnitude, either may dominate, depending on the particular case.

8

slide-31
SLIDE 31

Let us consider an idealized streamfunction on a midlati- tude β-plane comprising a zonally averaged part and a wave disturbance Φ = Φ0 − f0¯ uy + f0A sin kx cos ℓy

9

slide-32
SLIDE 32

Let us consider an idealized streamfunction on a midlati- tude β-plane comprising a zonally averaged part and a wave disturbance Φ = Φ0 − f0¯ uy + f0A sin kx cos ℓy The parameters Φ0, ¯ u and A depend only on pressure and the wavenumbers are k = 2π/Lx and ℓ = 2π/Ly.

9

slide-33
SLIDE 33

Let us consider an idealized streamfunction on a midlati- tude β-plane comprising a zonally averaged part and a wave disturbance Φ = Φ0 − f0¯ uy + f0A sin kx cos ℓy The parameters Φ0, ¯ u and A depend only on pressure and the wavenumbers are k = 2π/Lx and ℓ = 2π/Ly. The geostrophic winds are ug = − 1 f0 ∂Φ ∂y = ¯ u + u′

g = ¯

u + ℓA sin kx sin ℓy vg = + 1 f0 ∂Φ ∂x = + v′

g =

+ kA cos kx cos ℓy

9

slide-34
SLIDE 34

Let us consider an idealized streamfunction on a midlati- tude β-plane comprising a zonally averaged part and a wave disturbance Φ = Φ0 − f0¯ uy + f0A sin kx cos ℓy The parameters Φ0, ¯ u and A depend only on pressure and the wavenumbers are k = 2π/Lx and ℓ = 2π/Ly. The geostrophic winds are ug = − 1 f0 ∂Φ ∂y = ¯ u + u′

g = ¯

u + ℓA sin kx sin ℓy vg = + 1 f0 ∂Φ ∂x = + v′

g =

+ kA cos kx cos ℓy The geostrophic vorticity is ζg = 1 f0 ∇2Φ = −(k2 + ℓ2)A sin kx cos ℓy

9

slide-35
SLIDE 35

It is easily shown that the advection of relative vorticity by the wave component vanishes, u′

g

∂ζg ∂x + v′

g

∂ζg ∂y = 0

10

slide-36
SLIDE 36

It is easily shown that the advection of relative vorticity by the wave component vanishes, u′

g

∂ζg ∂x + v′

g

∂ζg ∂y = 0 Thus, the advection of relative vorticity reduces to −Vg·∇ζg = −¯ u ∂ζg ∂x = k¯ u(k2 + ℓ2)A cos kx cos ℓy

10

slide-37
SLIDE 37

It is easily shown that the advection of relative vorticity by the wave component vanishes, u′

g

∂ζg ∂x + v′

g

∂ζg ∂y = 0 Thus, the advection of relative vorticity reduces to −Vg·∇ζg = −¯ u ∂ζg ∂x = k¯ u(k2 + ℓ2)A cos kx cos ℓy The advection of planetary vorticity is −βvg = −βkA cos kx cos ℓy

10

slide-38
SLIDE 38

It is easily shown that the advection of relative vorticity by the wave component vanishes, u′

g

∂ζg ∂x + v′

g

∂ζg ∂y = 0 Thus, the advection of relative vorticity reduces to −Vg·∇ζg = −¯ u ∂ζg ∂x = k¯ u(k2 + ℓ2)A cos kx cos ℓy The advection of planetary vorticity is −βvg = −βkA cos kx cos ℓy The total vorticity advection is −Vg·∇(ζg + f) = [¯ u(k2 + ℓ2) − β]kA cos kx cos ℓy

10

slide-39
SLIDE 39

Repeat: The total vorticity advection is −Vg·∇(ζg + f) = [¯ u(k2 + ℓ2) − β]kA cos kx cos ℓy

11

slide-40
SLIDE 40

Repeat: The total vorticity advection is −Vg·∇(ζg + f) = [¯ u(k2 + ℓ2) − β]kA cos kx cos ℓy For relatively short wavelengths (L ≪ 3, 000 km) the advec- tion of relative vorticity dominates. For planetary-scale waves (L ∼ 10, 000 km) the β-term dominates and the waves regress.

11

slide-41
SLIDE 41

Repeat: The total vorticity advection is −Vg·∇(ζg + f) = [¯ u(k2 + ℓ2) − β]kA cos kx cos ℓy For relatively short wavelengths (L ≪ 3, 000 km) the advec- tion of relative vorticity dominates. For planetary-scale waves (L ∼ 10, 000 km) the β-term dominates and the waves regress. Thus, as a general rule, short-wavelength synoptic-scale dis- turbances should move eastward in a westerly flow. Long planetary waves regress or remain stationary.

11

slide-42
SLIDE 42

The Tendency Equation

12

slide-43
SLIDE 43

The Tendency Equation

Although the vertical velocity plays an essential role in the dynamics, the evolution of the geostrophic circulation can be determined without explicitly determining the distribu- tion of ω.

12

slide-44
SLIDE 44

The Tendency Equation

Although the vertical velocity plays an essential role in the dynamics, the evolution of the geostrophic circulation can be determined without explicitly determining the distribu- tion of ω. The vorticity equation is ∂ζg ∂t = −Vg·∇(ζg + f) + f0 ∂ω ∂p

12

slide-45
SLIDE 45

The Tendency Equation

Although the vertical velocity plays an essential role in the dynamics, the evolution of the geostrophic circulation can be determined without explicitly determining the distribu- tion of ω. The vorticity equation is ∂ζg ∂t = −Vg·∇(ζg + f) + f0 ∂ω ∂p Recalling that the vorticity and geopotential are related by ζg = (1/f0)∇2Φ and reversing the order of differentiation, we get 1 f0 ∇2Φt = −Vg·∇ 1 f0 ∇2Φ + f

  • + f0

∂ω ∂p

Note: Φt ≡ ∂Φ/∂t. Holton uses χ.

12

slide-46
SLIDE 46

The thermodynamic equation is ∂ ∂t + Vg · ∇ ∂Φ ∂p

  • + σω = −κ ˙

Q p

13

slide-47
SLIDE 47

The thermodynamic equation is ∂ ∂t + Vg · ∇ ∂Φ ∂p

  • + σω = −κ ˙

Q p Let us multiply by f0/σ and differentiate with respect to p: ∂ ∂p f0 σ ∂Φt ∂p

  • = − ∂

∂p f0 σ Vg · ∇ ∂Φ ∂p

  • − f0

∂ω ∂p − f0 ∂ ∂p

  • κ ˙

Q σp

  • 13
slide-48
SLIDE 48

The thermodynamic equation is ∂ ∂t + Vg · ∇ ∂Φ ∂p

  • + σω = −κ ˙

Q p Let us multiply by f0/σ and differentiate with respect to p: ∂ ∂p f0 σ ∂Φt ∂p

  • = − ∂

∂p f0 σ Vg · ∇ ∂Φ ∂p

  • − f0

∂ω ∂p − f0 ∂ ∂p

  • κ ˙

Q σp

  • We now ignore the effects of diabatic heating and set ˙

Q = 0.

13

slide-49
SLIDE 49

The thermodynamic equation is ∂ ∂t + Vg · ∇ ∂Φ ∂p

  • + σω = −κ ˙

Q p Let us multiply by f0/σ and differentiate with respect to p: ∂ ∂p f0 σ ∂Φt ∂p

  • = − ∂

∂p f0 σ Vg · ∇ ∂Φ ∂p

  • − f0

∂ω ∂p − f0 ∂ ∂p

  • κ ˙

Q σp

  • We now ignore the effects of diabatic heating and set ˙

Q = 0. It is simple to eliminate ω by addition of the thermodynamic and vorticity equations as expressed above. We then get

  • ∇2 + ∂

∂p

  • f2

σ ∂ ∂p

  • Φt
  • A

= − f0Vg·∇ 1 f0 ∇2Φ + f

  • B

+ ∂ ∂p

  • f2

σ Vg · ∇

  • −∂Φ

∂p

  • C

13

slide-50
SLIDE 50

Again, ∇2 + ∂ ∂p

  • f2

σ ∂ ∂p

  • Φt
  • A

= − f0Vg·∇ 1 f0 ∇2Φ + f

  • B

+ ∂ ∂p

  • f2

σ Vg · ∇

  • −∂Φ

∂p

  • C

14

slide-51
SLIDE 51

Again, ∇2 + ∂ ∂p

  • f2

σ ∂ ∂p

  • Φt
  • A

= − f0Vg·∇ 1 f0 ∇2Φ + f

  • B

+ ∂ ∂p

  • f2

σ Vg · ∇

  • −∂Φ

∂p

  • C

This is the geopotential tendency equation. It provides a relationship between

14

slide-52
SLIDE 52

Again, ∇2 + ∂ ∂p

  • f2

σ ∂ ∂p

  • Φt
  • A

= − f0Vg·∇ 1 f0 ∇2Φ + f

  • B

+ ∂ ∂p

  • f2

σ Vg · ∇

  • −∂Φ

∂p

  • C

This is the geopotential tendency equation. It provides a relationship between Term A: The local geopotential tendency Φt

14

slide-53
SLIDE 53

Again, ∇2 + ∂ ∂p

  • f2

σ ∂ ∂p

  • Φt
  • A

= − f0Vg·∇ 1 f0 ∇2Φ + f

  • B

+ ∂ ∂p

  • f2

σ Vg · ∇

  • −∂Φ

∂p

  • C

This is the geopotential tendency equation. It provides a relationship between Term A: The local geopotential tendency Φt Term B: The advection of vorticity

14

slide-54
SLIDE 54

Again, ∇2 + ∂ ∂p

  • f2

σ ∂ ∂p

  • Φt
  • A

= − f0Vg·∇ 1 f0 ∇2Φ + f

  • B

+ ∂ ∂p

  • f2

σ Vg · ∇

  • −∂Φ

∂p

  • C

This is the geopotential tendency equation. It provides a relationship between Term A: The local geopotential tendency Φt Term B: The advection of vorticity Term C: The vertical shear of temperature advection.

14

slide-55
SLIDE 55

Term (A) involves second derivatives with respect to spatial variables of the geopotential tendency Φt. For sinusoidal variations, this is typically proportional to −Φt.

15

slide-56
SLIDE 56

Term (A) involves second derivatives with respect to spatial variables of the geopotential tendency Φt. For sinusoidal variations, this is typically proportional to −Φt. Term (B) is proportional to the advection of absolute vor-

  • ticity. For the upper troposphere it is usually the dominant

term.

15

slide-57
SLIDE 57

Term (A) involves second derivatives with respect to spatial variables of the geopotential tendency Φt. For sinusoidal variations, this is typically proportional to −Φt. Term (B) is proportional to the advection of absolute vor-

  • ticity. For the upper troposphere it is usually the dominant

term. For short waves we have seen that the relative vorticity ad- vection dominates the planetary vorticity advection. With a ridge to the west and a trough to the east, this term is then negative.

15

slide-58
SLIDE 58

Term (A) involves second derivatives with respect to spatial variables of the geopotential tendency Φt. For sinusoidal variations, this is typically proportional to −Φt. Term (B) is proportional to the advection of absolute vor-

  • ticity. For the upper troposphere it is usually the dominant

term. For short waves we have seen that the relative vorticity ad- vection dominates the planetary vorticity advection. With a ridge to the west and a trough to the east, this term is then negative. Thus, Term (B) makes Φt positive, so that a ridge tends to develop and, associated with this, the vorticity becomes negative.

15

slide-59
SLIDE 59

Term (A) involves second derivatives with respect to spatial variables of the geopotential tendency Φt. For sinusoidal variations, this is typically proportional to −Φt. Term (B) is proportional to the advection of absolute vor-

  • ticity. For the upper troposphere it is usually the dominant

term. For short waves we have seen that the relative vorticity ad- vection dominates the planetary vorticity advection. With a ridge to the west and a trough to the east, this term is then negative. Thus, Term (B) makes Φt positive, so that a ridge tends to develop and, associated with this, the vorticity becomes negative. Term (B) acts to transport the pattern of geopotential. However, since Vg·∇ζ = 0 on the trough and ridge axes, this term does not cause the wave to amplify or decay.

15

slide-60
SLIDE 60

The means of amplification or decay of midlatitude waves is contained in Term (C). This term is proportional to minus the rate of change of temperature advection with respect to pressure.

16

slide-61
SLIDE 61

The means of amplification or decay of midlatitude waves is contained in Term (C). This term is proportional to minus the rate of change of temperature advection with respect to pressure. It is therefore related to plus the rate of change of temper- ature advection with respect to height. This is called the differential temperature advection.

16

slide-62
SLIDE 62

The means of amplification or decay of midlatitude waves is contained in Term (C). This term is proportional to minus the rate of change of temperature advection with respect to pressure. It is therefore related to plus the rate of change of temper- ature advection with respect to height. This is called the differential temperature advection. The magnitude of the temperature (or thickness) advection tends to be largest in the lower troposphere, beneath the 500 hPa trough and ridge lines in a developing baroclinic wave.

16

slide-63
SLIDE 63

17

slide-64
SLIDE 64
  • Below the 500 hPa ridge, there is warm advection as-

sociated with the advancing warm front. This increases thickness and builds the upper level ridge.

18

slide-65
SLIDE 65
  • Below the 500 hPa ridge, there is warm advection as-

sociated with the advancing warm front. This increases thickness and builds the upper level ridge.

  • Below the 500 hPa trough, there is cold advection as-

sociated with the advancing cold front. This decreases thickness and deepens the upper level trough.

18

slide-66
SLIDE 66
  • Below the 500 hPa ridge, there is warm advection as-

sociated with the advancing warm front. This increases thickness and builds the upper level ridge.

  • Below the 500 hPa trough, there is cold advection as-

sociated with the advancing cold front. This decreases thickness and deepens the upper level trough. Thus in contrast to term (B), term (C) is dominant in the lower troposphere; but its effect is felt at higher levels.

18

slide-67
SLIDE 67
  • Below the 500 hPa ridge, there is warm advection as-

sociated with the advancing warm front. This increases thickness and builds the upper level ridge.

  • Below the 500 hPa trough, there is cold advection as-

sociated with the advancing cold front. This decreases thickness and deepens the upper level trough. Thus in contrast to term (B), term (C) is dominant in the lower troposphere; but its effect is felt at higher levels. In words, we may write the geopotential tendency equation: Falling Pressure

  • Positive

Vorticity Advection

  • +
  • Differential

Temperature Advec’n

  • 18
slide-68
SLIDE 68

Tendency due to vorticity advection Tendence due to diff’l temperature advection

19

slide-69
SLIDE 69

The Omega Equation

20

slide-70
SLIDE 70

The Omega Equation

We will now eliminate the geopotential tendency by com- bining the momentum and thermodynamic equations, and

  • btain an equation for the vertical velocity ω.

20

slide-71
SLIDE 71

The Omega Equation

We will now eliminate the geopotential tendency by com- bining the momentum and thermodynamic equations, and

  • btain an equation for the vertical velocity ω.

The thermodynamic equation for adiabatic flow is ∂ ∂t + Vg · ∇ ∂Φ ∂p

  • + σω = 0

20

slide-72
SLIDE 72

The Omega Equation

We will now eliminate the geopotential tendency by com- bining the momentum and thermodynamic equations, and

  • btain an equation for the vertical velocity ω.

The thermodynamic equation for adiabatic flow is ∂ ∂t + Vg · ∇ ∂Φ ∂p

  • + σω = 0

We now write it as ∂Φt ∂p = −Vg · ∇ ∂Φ ∂p

  • − σω

20

slide-73
SLIDE 73

The Omega Equation

We will now eliminate the geopotential tendency by com- bining the momentum and thermodynamic equations, and

  • btain an equation for the vertical velocity ω.

The thermodynamic equation for adiabatic flow is ∂ ∂t + Vg · ∇ ∂Φ ∂p

  • + σω = 0

We now write it as ∂Φt ∂p = −Vg · ∇ ∂Φ ∂p

  • − σω

We take the Laplacian of this and obtain

∇2∂Φt ∂p = −∇2

  • Vg · ∇

∂Φ ∂p

  • − σ∇2ω

20

slide-74
SLIDE 74

Recall that the vorticity equation may be written ∂ζg ∂t = −Vg·∇(ζg + f) + f0 ∂ω ∂p

21

slide-75
SLIDE 75

Recall that the vorticity equation may be written ∂ζg ∂t = −Vg·∇(ζg + f) + f0 ∂ω ∂p Multiply by f0 and use f0ζg = ∇2Φ: ∇2Φt = −f0Vg·∇( 1 f0 ∇2Φ + f) + f2 ∂ω ∂p

21

slide-76
SLIDE 76

Recall that the vorticity equation may be written ∂ζg ∂t = −Vg·∇(ζg + f) + f0 ∂ω ∂p Multiply by f0 and use f0ζg = ∇2Φ: ∇2Φt = −f0Vg·∇( 1 f0 ∇2Φ + f) + f2 ∂ω ∂p Now differentiate with respect to pressure:

∇2∂Φt ∂p = −f0 ∂ ∂p

  • Vg·∇

1 f0 ∇2Φ + f

  • + f 2

∂2ω ∂p2

21

slide-77
SLIDE 77

Recall that the vorticity equation may be written ∂ζg ∂t = −Vg·∇(ζg + f) + f0 ∂ω ∂p Multiply by f0 and use f0ζg = ∇2Φ: ∇2Φt = −f0Vg·∇( 1 f0 ∇2Φ + f) + f2 ∂ω ∂p Now differentiate with respect to pressure:

∇2∂Φt ∂p = −f0 ∂ ∂p

  • Vg·∇

1 f0 ∇2Φ + f

  • + f 2

∂2ω ∂p2

We now have two equations with identical expressions for the tendency Φt. So we can subtract one from the other to obtain a diagnostic equation for the vertical velocity.

21

slide-78
SLIDE 78
  • σ∇2 + f2

∂2 ∂p2

  • ω
  • A

= − f0 ∂ ∂p

  • −Vg·∇

1 f0 ∇2Φ + f

  • B

+ ∇2

  • Vg · ∇
  • −∂Φ

∂p

  • C

22

slide-79
SLIDE 79
  • σ∇2 + f2

∂2 ∂p2

  • ω
  • A

= − f0 ∂ ∂p

  • −Vg·∇

1 f0 ∇2Φ + f

  • B

+ ∇2

  • Vg · ∇
  • −∂Φ

∂p

  • C

This is the omega equation, a diagnostic relationship for the vertical velocity ω.

22

slide-80
SLIDE 80
  • σ∇2 + f2

∂2 ∂p2

  • ω
  • A

= − f0 ∂ ∂p

  • −Vg·∇

1 f0 ∇2Φ + f

  • B

+ ∇2

  • Vg · ∇
  • −∂Φ

∂p

  • C

This is the omega equation, a diagnostic relationship for the vertical velocity ω. It provides a relationship between

22

slide-81
SLIDE 81
  • σ∇2 + f2

∂2 ∂p2

  • ω
  • A

= − f0 ∂ ∂p

  • −Vg·∇

1 f0 ∇2Φ + f

  • B

+ ∇2

  • Vg · ∇
  • −∂Φ

∂p

  • C

This is the omega equation, a diagnostic relationship for the vertical velocity ω. It provides a relationship between Term A: The vertical velocity

22

slide-82
SLIDE 82
  • σ∇2 + f2

∂2 ∂p2

  • ω
  • A

= − f0 ∂ ∂p

  • −Vg·∇

1 f0 ∇2Φ + f

  • B

+ ∇2

  • Vg · ∇
  • −∂Φ

∂p

  • C

This is the omega equation, a diagnostic relationship for the vertical velocity ω. It provides a relationship between Term A: The vertical velocity Term B: The differential advection of vorticity

22

slide-83
SLIDE 83
  • σ∇2 + f2

∂2 ∂p2

  • ω
  • A

= − f0 ∂ ∂p

  • −Vg·∇

1 f0 ∇2Φ + f

  • B

+ ∇2

  • Vg · ∇
  • −∂Φ

∂p

  • C

This is the omega equation, a diagnostic relationship for the vertical velocity ω. It provides a relationship between Term A: The vertical velocity Term B: The differential advection of vorticity Term C: The temperature advection.

22

slide-84
SLIDE 84

Term (A), the left side of the equation, involves spatial sec-

  • nd derivatives of ω.

For sinusoidal variations it is proportional to the negative

  • f ω and is thus related directly to the vertical velocity w.

23

slide-85
SLIDE 85

Term (A), the left side of the equation, involves spatial sec-

  • nd derivatives of ω.

For sinusoidal variations it is proportional to the negative

  • f ω and is thus related directly to the vertical velocity w.

Term (B) is the change with pressure of the advection of absolute vorticity, that is, the differential vorticity advec- tion.

23

slide-86
SLIDE 86

Term (A), the left side of the equation, involves spatial sec-

  • nd derivatives of ω.

For sinusoidal variations it is proportional to the negative

  • f ω and is thus related directly to the vertical velocity w.

Term (B) is the change with pressure of the advection of absolute vorticity, that is, the differential vorticity advec- tion. Term (C) is the Laplacian of minus the temperature advec- tion, and is thus proportional to the advection of tempera- ture.

23

slide-87
SLIDE 87

Term (A), the left side of the equation, involves spatial sec-

  • nd derivatives of ω.

For sinusoidal variations it is proportional to the negative

  • f ω and is thus related directly to the vertical velocity w.

Term (B) is the change with pressure of the advection of absolute vorticity, that is, the differential vorticity advec- tion. Term (C) is the Laplacian of minus the temperature advec- tion, and is thus proportional to the advection of tempera- ture. In words, we may write the omega equation as follows Rising Motion

  • Differential

Vorticity Advection

  • +

Temperature Advection

  • 23
slide-88
SLIDE 88

Idealized baroclinic wave. Solid: 500 hPa geopotential con-

  • tours. Dashed: 1000 hPa contours. Regions of strong ver-

tical motion due to differential vorticity advection are indi- cated.

24

slide-89
SLIDE 89

Term (B) is −f0 ∂ ∂p

  • −Vg·∇

1 f0 ∇2Φ + f

  • ∝ ∂

∂z

  • −Vg·∇
  • ζg + f
  • so it is proportional to differential vorticity advection.

25

slide-90
SLIDE 90

Term (B) is −f0 ∂ ∂p

  • −Vg·∇

1 f0 ∇2Φ + f

  • ∝ ∂

∂z

  • −Vg·∇
  • ζg + f
  • so it is proportional to differential vorticity advection.
  • At the surface Low , the advection of vorticity is small

25

slide-91
SLIDE 91

Term (B) is −f0 ∂ ∂p

  • −Vg·∇

1 f0 ∇2Φ + f

  • ∝ ∂

∂z

  • −Vg·∇
  • ζg + f
  • so it is proportional to differential vorticity advection.
  • At the surface Low , the advection of vorticity is small
  • Above this, there is strong positive vorticity advection at

500 hPa

25

slide-92
SLIDE 92

Term (B) is −f0 ∂ ∂p

  • −Vg·∇

1 f0 ∇2Φ + f

  • ∝ ∂

∂z

  • −Vg·∇
  • ζg + f
  • so it is proportional to differential vorticity advection.
  • At the surface Low , the advection of vorticity is small
  • Above this, there is strong positive vorticity advection at

500 hPa

  • Therefore, the differential vorticity advection is positive

25

slide-93
SLIDE 93

Term (B) is −f0 ∂ ∂p

  • −Vg·∇

1 f0 ∇2Φ + f

  • ∝ ∂

∂z

  • −Vg·∇
  • ζg + f
  • so it is proportional to differential vorticity advection.
  • At the surface Low , the advection of vorticity is small
  • Above this, there is strong positive vorticity advection at

500 hPa

  • Therefore, the differential vorticity advection is positive
  • This indices an upward vertical velocity

25

slide-94
SLIDE 94

Term (B) is −f0 ∂ ∂p

  • −Vg·∇

1 f0 ∇2Φ + f

  • ∝ ∂

∂z

  • −Vg·∇
  • ζg + f
  • so it is proportional to differential vorticity advection.
  • At the surface Low , the advection of vorticity is small
  • Above this, there is strong positive vorticity advection at

500 hPa

  • Therefore, the differential vorticity advection is positive
  • This indices an upward vertical velocity
  • Correspondingly, w < 0 above the surface High Pressure

25

slide-95
SLIDE 95

Term (B) is −f0 ∂ ∂p

  • −Vg·∇

1 f0 ∇2Φ + f

  • ∝ ∂

∂z

  • −Vg·∇
  • ζg + f
  • so it is proportional to differential vorticity advection.
  • At the surface Low , the advection of vorticity is small
  • Above this, there is strong positive vorticity advection at

500 hPa

  • Therefore, the differential vorticity advection is positive
  • This indices an upward vertical velocity
  • Correspondingly, w < 0 above the surface High Pressure
  • We assume the scale is short enough that relative vortic-

ity advection dominates planetary vorticity advection. Conclusion: Differential vorticity advection implies: Rising motion above the surface low Subsidence above the surface High.

25

slide-96
SLIDE 96

Idealized baroclinic wave. Solid: 500 hPa geopotential con-

  • tours. Dashed: 1000 hPa contours. Regions of strong ver-

tical motion due to temperature advection are indicated.

26

slide-97
SLIDE 97

Term (C) is +∇2

  • Vg · ∇
  • −∂Φ

∂p

  • ∝ −Vg · ∇T

so it is propotrional to the temperature advection.

27

slide-98
SLIDE 98

Term (C) is +∇2

  • Vg · ∇
  • −∂Φ

∂p

  • ∝ −Vg · ∇T

so it is propotrional to the temperature advection.

  • Ahead of the surface Low there is warm advection

27

slide-99
SLIDE 99

Term (C) is +∇2

  • Vg · ∇
  • −∂Φ

∂p

  • ∝ −Vg · ∇T

so it is propotrional to the temperature advection.

  • Ahead of the surface Low there is warm advection
  • Therefore, Term (C) is positive

27

slide-100
SLIDE 100

Term (C) is +∇2

  • Vg · ∇
  • −∂Φ

∂p

  • ∝ −Vg · ∇T

so it is propotrional to the temperature advection.

  • Ahead of the surface Low there is warm advection
  • Therefore, Term (C) is positive
  • So, there is Rising Motion ahead of the Low centre.

27

slide-101
SLIDE 101

Term (C) is +∇2

  • Vg · ∇
  • −∂Φ

∂p

  • ∝ −Vg · ∇T

so it is propotrional to the temperature advection.

  • Ahead of the surface Low there is warm advection
  • Therefore, Term (C) is positive
  • So, there is Rising Motion ahead of the Low centre.
  • Behind the Low, the cold front is associated with cold

advection

27

slide-102
SLIDE 102

Term (C) is +∇2

  • Vg · ∇
  • −∂Φ

∂p

  • ∝ −Vg · ∇T

so it is propotrional to the temperature advection.

  • Ahead of the surface Low there is warm advection
  • Therefore, Term (C) is positive
  • So, there is Rising Motion ahead of the Low centre.
  • Behind the Low, the cold front is associated with cold

advection

  • Hence, there is subsidence at the 500 hPa trough

27

slide-103
SLIDE 103

Vertical motion due to differential vorticity advection. Vertical motion due to temperature advection.

28

slide-104
SLIDE 104

Summary

29

slide-105
SLIDE 105

Summary

For synoptic scale motions, the flow is approximately in geo- strophic balance.

29

slide-106
SLIDE 106

Summary

For synoptic scale motions, the flow is approximately in geo- strophic balance. For purely geostrophic flow, the horizontal velocity is de- termined by the geopotential field.

29

slide-107
SLIDE 107

Summary

For synoptic scale motions, the flow is approximately in geo- strophic balance. For purely geostrophic flow, the horizontal velocity is de- termined by the geopotential field. The QG system allows us to determine both the geostrophic and ageostrophic components of the flow.

29

slide-108
SLIDE 108

Summary

For synoptic scale motions, the flow is approximately in geo- strophic balance. For purely geostrophic flow, the horizontal velocity is de- termined by the geopotential field. The QG system allows us to determine both the geostrophic and ageostrophic components of the flow. The vertical velocity is also determined by the geopotential field.

29

slide-109
SLIDE 109

Summary

For synoptic scale motions, the flow is approximately in geo- strophic balance. For purely geostrophic flow, the horizontal velocity is de- termined by the geopotential field. The QG system allows us to determine both the geostrophic and ageostrophic components of the flow. The vertical velocity is also determined by the geopotential field. This vertical velocity is just that required to ensure that the vorticity remains geostrophic and the temperature remains in hydrostatic balance.

29

slide-110
SLIDE 110

Summary

For synoptic scale motions, the flow is approximately in geo- strophic balance. For purely geostrophic flow, the horizontal velocity is de- termined by the geopotential field. The QG system allows us to determine both the geostrophic and ageostrophic components of the flow. The vertical velocity is also determined by the geopotential field. This vertical velocity is just that required to ensure that the vorticity remains geostrophic and the temperature remains in hydrostatic balance. The Tendence Equation allows us to predict the evolution

  • f the mass field and the diagnostic relationships then yield

all the other fields.

29

slide-111
SLIDE 111
  • Recap. on Φt and ω Equations

30

slide-112
SLIDE 112
  • Recap. on Φt and ω Equations

The Tendency Equation is Falling Pressure

  • Positive

Vorticity Advection

  • +
  • Differential

Temperature Advec’n

  • 30
slide-113
SLIDE 113
  • Recap. on Φt and ω Equations

The Tendency Equation is Falling Pressure

  • Positive

Vorticity Advection

  • +
  • Differential

Temperature Advec’n

  • The Omega Equation is

Rising Motion

  • Differential

Vorticity Advection

  • +

Temperature Advection

  • 30
slide-114
SLIDE 114
  • Recap. on Φt and ω Equations

The Tendency Equation is Falling Pressure

  • Positive

Vorticity Advection

  • +
  • Differential

Temperature Advec’n

  • The Omega Equation is

Rising Motion

  • Differential

Vorticity Advection

  • +

Temperature Advection

  • Note the complimentarity between these two equations.

30

slide-115
SLIDE 115

Exercise:

Study a chart of the 850 hPa or 700 hPa temperature and geopotential.

31

slide-116
SLIDE 116

Exercise:

Study a chart of the 850 hPa or 700 hPa temperature and geopotential.

  • Pick out the areas of maximum baroclinicity.

31

slide-117
SLIDE 117

Exercise:

Study a chart of the 850 hPa or 700 hPa temperature and geopotential.

  • Pick out the areas of maximum baroclinicity.
  • How are they related to surface Frontal zones?

31

slide-118
SLIDE 118

Exercise:

Study a chart of the 850 hPa or 700 hPa temperature and geopotential.

  • Pick out the areas of maximum baroclinicity.
  • How are they related to surface Frontal zones?
  • Identify where warm and cold advection are taking place.

31

slide-119
SLIDE 119

Exercise:

Study a chart of the 850 hPa or 700 hPa temperature and geopotential.

  • Pick out the areas of maximum baroclinicity.
  • How are they related to surface Frontal zones?
  • Identify where warm and cold advection are taking place.
  • Identify regions of vorticity advection.

31

slide-120
SLIDE 120

Exercise:

Study a chart of the 850 hPa or 700 hPa temperature and geopotential.

  • Pick out the areas of maximum baroclinicity.
  • How are they related to surface Frontal zones?
  • Identify where warm and cold advection are taking place.
  • Identify regions of vorticity advection.
  • Draw deductions about the vertical velocity.

31

slide-121
SLIDE 121

Exercise:

Study a chart of the 850 hPa or 700 hPa temperature and geopotential.

  • Pick out the areas of maximum baroclinicity.
  • How are they related to surface Frontal zones?
  • Identify where warm and cold advection are taking place.
  • Identify regions of vorticity advection.
  • Draw deductions about the vertical velocity.
  • How is the vertical velocity correlated with the geopoten-

tial field?

31