The Potential Vorticity Equation The Potential Vorticity Equation - - PowerPoint PPT Presentation

the potential vorticity equation the potential vorticity
SMART_READER_LITE
LIVE PREVIEW

The Potential Vorticity Equation The Potential Vorticity Equation - - PowerPoint PPT Presentation

The Potential Vorticity Equation The Potential Vorticity Equation The geopotential tendency equation is 1 f 2 2 + 2 + f 0 t = f 0 V g p p f 0 f 2 +


slide-1
SLIDE 1

The Potential Vorticity Equation

slide-2
SLIDE 2

The Potential Vorticity Equation

The geopotential tendency equation is

  • ∇2 + ∂

∂p

  • f2

σ ∂ ∂p

  • Φt = − f0Vg·∇

1 f0 ∇2Φ + f

  • + ∂

∂p

  • f2

σ Vg · ∇

  • −∂Φ

∂p

slide-3
SLIDE 3

The Potential Vorticity Equation

The geopotential tendency equation is

  • ∇2 + ∂

∂p

  • f2

σ ∂ ∂p

  • Φt = − f0Vg·∇

1 f0 ∇2Φ + f

  • + ∂

∂p

  • f2

σ Vg · ∇

  • −∂Φ

∂p

  • The second term on the right (Term (C)) may be expanded:

−Vg · ∇ ∂ ∂p

  • f2

σ ∂Φ ∂p

  • − f2

σ ∂Vg ∂p ·∇∂Φ ∂p

slide-4
SLIDE 4

The Potential Vorticity Equation

The geopotential tendency equation is

  • ∇2 + ∂

∂p

  • f2

σ ∂ ∂p

  • Φt = − f0Vg·∇

1 f0 ∇2Φ + f

  • + ∂

∂p

  • f2

σ Vg · ∇

  • −∂Φ

∂p

  • The second term on the right (Term (C)) may be expanded:

−Vg · ∇ ∂ ∂p

  • f2

σ ∂Φ ∂p

  • − f2

σ ∂Vg ∂p ·∇∂Φ ∂p But the thermal wind relationship is f0 ∂Vg ∂p = k × ∇∂Φ ∂p This is just the p-derivative of f0Vg = k × ∇Φ.

slide-5
SLIDE 5

Thus, ∂Vg/∂p is perpendicular to ∇(∂Φ/∂p) and the second term above vanishes.

2

slide-6
SLIDE 6

Thus, ∂Vg/∂p is perpendicular to ∇(∂Φ/∂p) and the second term above vanishes. The remaining term can be combined with term (B) in the tendency equation to give RHS = −f0Vg·∇ 1 f0 ∇2Φ + f + ∂ ∂p f0 σ ∂Φ ∂p

  • = −f0Vg·∇q

2

slide-7
SLIDE 7

Thus, ∂Vg/∂p is perpendicular to ∇(∂Φ/∂p) and the second term above vanishes. The remaining term can be combined with term (B) in the tendency equation to give RHS = −f0Vg·∇ 1 f0 ∇2Φ + f + ∂ ∂p f0 σ ∂Φ ∂p

  • = −f0Vg·∇q

The quantity in square brackets is called the quasi-geostrophic potential vorticity q ≡ 1 f0 ∇2Φ + f + ∂ ∂p f0 σ ∂Φ ∂p

  • 2
slide-8
SLIDE 8

Thus, ∂Vg/∂p is perpendicular to ∇(∂Φ/∂p) and the second term above vanishes. The remaining term can be combined with term (B) in the tendency equation to give RHS = −f0Vg·∇ 1 f0 ∇2Φ + f + ∂ ∂p f0 σ ∂Φ ∂p

  • = −f0Vg·∇q

The quantity in square brackets is called the quasi-geostrophic potential vorticity q ≡ 1 f0 ∇2Φ + f + ∂ ∂p f0 σ ∂Φ ∂p

  • The left side of the tendency equation may be written

LHS = f0 ∂ ∂t 1 f0 ∇2Φ + ∂ ∂p f0 σ ∂Φ ∂p

  • = f0

∂q ∂t since f does not vary with time.

2

slide-9
SLIDE 9

The tendency equation may now be written in a conserv- ative form called the quasi-geostrophic potential vorticity equation or QGPV equation: ∂ ∂t + Vg·∇

  • q = dgq

dt = 0 The potential vorticity q is comprised of three terms:

3

slide-10
SLIDE 10

The tendency equation may now be written in a conserv- ative form called the quasi-geostrophic potential vorticity equation or QGPV equation: ∂ ∂t + Vg·∇

  • q = dgq

dt = 0 The potential vorticity q is comprised of three terms:

  • The relative vorticity, ζg

3

slide-11
SLIDE 11

The tendency equation may now be written in a conserv- ative form called the quasi-geostrophic potential vorticity equation or QGPV equation: ∂ ∂t + Vg·∇

  • q = dgq

dt = 0 The potential vorticity q is comprised of three terms:

  • The relative vorticity, ζg
  • The planetary vorticity f

3

slide-12
SLIDE 12

The tendency equation may now be written in a conserv- ative form called the quasi-geostrophic potential vorticity equation or QGPV equation: ∂ ∂t + Vg·∇

  • q = dgq

dt = 0 The potential vorticity q is comprised of three terms:

  • The relative vorticity, ζg
  • The planetary vorticity f
  • The stretching vorticity, (∂/∂p)[(f0/σ)∂Φ/∂p].

3

slide-13
SLIDE 13

The tendency equation may now be written in a conserv- ative form called the quasi-geostrophic potential vorticity equation or QGPV equation: ∂ ∂t + Vg·∇

  • q = dgq

dt = 0 The potential vorticity q is comprised of three terms:

  • The relative vorticity, ζg
  • The planetary vorticity f
  • The stretching vorticity, (∂/∂p)[(f0/σ)∂Φ/∂p].

The equation states that q is conserved following the geo- strophic flow.

3

slide-14
SLIDE 14

The tendency equation may now be written in a conserv- ative form called the quasi-geostrophic potential vorticity equation or QGPV equation: ∂ ∂t + Vg·∇

  • q = dgq

dt = 0 The potential vorticity q is comprised of three terms:

  • The relative vorticity, ζg
  • The planetary vorticity f
  • The stretching vorticity, (∂/∂p)[(f0/σ)∂Φ/∂p].

The equation states that q is conserved following the geo- strophic flow. Note that q is completely determined once the three-dim- ensional distribution of geopotential Φ is given.

3

slide-15
SLIDE 15

The QGPV equation may be used to predict the evolution

  • f atmospheric flows in midlatitudes.

4

slide-16
SLIDE 16

The QGPV equation may be used to predict the evolution

  • f atmospheric flows in midlatitudes.

Or, in plain language, to make Weather Forecasts.

4

slide-17
SLIDE 17

The QGPV equation may be used to predict the evolution

  • f atmospheric flows in midlatitudes.

Or, in plain language, to make Weather Forecasts. Indeed an even simpler equation, ignoring the stretching term, was used for the first computer forecasts on the ENIAC computer (Electronic Numerical Integrator and Computer) in 1950.

4

slide-18
SLIDE 18

The QGPV equation may be used to predict the evolution

  • f atmospheric flows in midlatitudes.

Or, in plain language, to make Weather Forecasts. Indeed an even simpler equation, ignoring the stretching term, was used for the first computer forecasts on the ENIAC computer (Electronic Numerical Integrator and Computer) in 1950. This was the equation

d dt(ζ + f) = 0

for the conservation of absolute vorticity.

4

slide-19
SLIDE 19

Exercise

5

slide-20
SLIDE 20

Exercise

An idealized geopotential field is given at time t = 0 by Φ = Φ0 − f0¯ uy + A sin(kx − mp) where Φ0, ¯ u and A are functions of p.

5

slide-21
SLIDE 21

Exercise

An idealized geopotential field is given at time t = 0 by Φ = Φ0 − f0¯ uy + A sin(kx − mp) where Φ0, ¯ u and A are functions of p. This represents a mean westerly flow ¯ u > 0 distorted by a wave disturbance of amplitude A.

5

slide-22
SLIDE 22

Exercise

An idealized geopotential field is given at time t = 0 by Φ = Φ0 − f0¯ uy + A sin(kx − mp) where Φ0, ¯ u and A are functions of p. This represents a mean westerly flow ¯ u > 0 distorted by a wave disturbance of amplitude A. (a) Compute the geostrophic wind components as functions

  • f x and p.

5

slide-23
SLIDE 23

Exercise

An idealized geopotential field is given at time t = 0 by Φ = Φ0 − f0¯ uy + A sin(kx − mp) where Φ0, ¯ u and A are functions of p. This represents a mean westerly flow ¯ u > 0 distorted by a wave disturbance of amplitude A. (a) Compute the geostrophic wind components as functions

  • f x and p.

(b) Ignoring the β-effect, compute the geostrophic vorticity and divergence.

5

slide-24
SLIDE 24

Exercise

An idealized geopotential field is given at time t = 0 by Φ = Φ0 − f0¯ uy + A sin(kx − mp) where Φ0, ¯ u and A are functions of p. This represents a mean westerly flow ¯ u > 0 distorted by a wave disturbance of amplitude A. (a) Compute the geostrophic wind components as functions

  • f x and p.

(b) Ignoring the β-effect, compute the geostrophic vorticity and divergence. (c) Compute the variations in the temperature field due to the wave.

5

slide-25
SLIDE 25

(d) Compute the vorticity advection and temperature advec- tion.

6

slide-26
SLIDE 26

(d) Compute the vorticity advection and temperature advec- tion. (e) Using the geopotential tendency equation, describe how the pressure at a point upstream from a trough and down- stream from a ridge is expected to change.

6

slide-27
SLIDE 27

(d) Compute the vorticity advection and temperature advec- tion. (e) Using the geopotential tendency equation, describe how the pressure at a point upstream from a trough and down- stream from a ridge is expected to change. (f) Using the omega equation, describe the pattern of vertical velocity associated with the wave disturbance.

6

slide-28
SLIDE 28

x

The ENIAC Integrations

(ENIAC: Electronic Numerical Integrator and Computer)

7

slide-29
SLIDE 29

Electronic Computer Project, 1946 (under direction of John von Neumann)

Von Neumann’s idea: Weather forecasting was a scientific problem par excellence for solution using a large computer.

8

slide-30
SLIDE 30

Electronic Computer Project, 1946 (under direction of John von Neumann)

Von Neumann’s idea: Weather forecasting was a scientific problem par excellence for solution using a large computer. The objective of the project was to study the problem of predicting the weather by simulating the dynamics of the atmosphere using a digital electronic computer.

8

slide-31
SLIDE 31

Electronic Computer Project, 1946 (under direction of John von Neumann)

Von Neumann’s idea: Weather forecasting was a scientific problem par excellence for solution using a large computer. The objective of the project was to study the problem of predicting the weather by simulating the dynamics of the atmosphere using a digital electronic computer. A Proposal for funding listed three “possibilities”:

8

slide-32
SLIDE 32

Electronic Computer Project, 1946 (under direction of John von Neumann)

Von Neumann’s idea: Weather forecasting was a scientific problem par excellence for solution using a large computer. The objective of the project was to study the problem of predicting the weather by simulating the dynamics of the atmosphere using a digital electronic computer. A Proposal for funding listed three “possibilities”:

  • 1. Entirely new methods of weather prediction by calcula-

tion will have been made possible;

8

slide-33
SLIDE 33

Electronic Computer Project, 1946 (under direction of John von Neumann)

Von Neumann’s idea: Weather forecasting was a scientific problem par excellence for solution using a large computer. The objective of the project was to study the problem of predicting the weather by simulating the dynamics of the atmosphere using a digital electronic computer. A Proposal for funding listed three “possibilities”:

  • 1. Entirely new methods of weather prediction by calcula-

tion will have been made possible;

  • 2. A new rational basis will have been secured for the plan-

ning of physical measurements and field observations;

8

slide-34
SLIDE 34

Electronic Computer Project, 1946 (under direction of John von Neumann)

Von Neumann’s idea: Weather forecasting was a scientific problem par excellence for solution using a large computer. The objective of the project was to study the problem of predicting the weather by simulating the dynamics of the atmosphere using a digital electronic computer. A Proposal for funding listed three “possibilities”:

  • 1. Entirely new methods of weather prediction by calcula-

tion will have been made possible;

  • 2. A new rational basis will have been secured for the plan-

ning of physical measurements and field observations;

  • 3. The first step towards influencing the weather by rational

human intervention will have been made.

8

slide-35
SLIDE 35

“Conference on Meteorology”

A “Conference on Meteorology” was arranged in the Insti- tute for Advanced Studies (IAS), Princeton on 29–30 Au- gust, 1946. Participants included:

9

slide-36
SLIDE 36

“Conference on Meteorology”

A “Conference on Meteorology” was arranged in the Insti- tute for Advanced Studies (IAS), Princeton on 29–30 Au- gust, 1946. Participants included:

  • Carl Gustav Rossby

9

slide-37
SLIDE 37

“Conference on Meteorology”

A “Conference on Meteorology” was arranged in the Insti- tute for Advanced Studies (IAS), Princeton on 29–30 Au- gust, 1946. Participants included:

  • Carl Gustav Rossby
  • Jule Charney

9

slide-38
SLIDE 38

“Conference on Meteorology”

A “Conference on Meteorology” was arranged in the Insti- tute for Advanced Studies (IAS), Princeton on 29–30 Au- gust, 1946. Participants included:

  • Carl Gustav Rossby
  • Jule Charney
  • George Platzman

9

slide-39
SLIDE 39

“Conference on Meteorology”

A “Conference on Meteorology” was arranged in the Insti- tute for Advanced Studies (IAS), Princeton on 29–30 Au- gust, 1946. Participants included:

  • Carl Gustav Rossby
  • Jule Charney
  • George Platzman
  • Norman Phillips

9

slide-40
SLIDE 40

“Conference on Meteorology”

A “Conference on Meteorology” was arranged in the Insti- tute for Advanced Studies (IAS), Princeton on 29–30 Au- gust, 1946. Participants included:

  • Carl Gustav Rossby
  • Jule Charney
  • George Platzman
  • Norman Phillips
  • Ragnar Fjørtoft

9

slide-41
SLIDE 41

“Conference on Meteorology”

A “Conference on Meteorology” was arranged in the Insti- tute for Advanced Studies (IAS), Princeton on 29–30 Au- gust, 1946. Participants included:

  • Carl Gustav Rossby
  • Jule Charney
  • George Platzman
  • Norman Phillips
  • Ragnar Fjørtoft
  • Arnt Eliassen

9

slide-42
SLIDE 42

“Conference on Meteorology”

A “Conference on Meteorology” was arranged in the Insti- tute for Advanced Studies (IAS), Princeton on 29–30 Au- gust, 1946. Participants included:

  • Carl Gustav Rossby
  • Jule Charney
  • George Platzman
  • Norman Phillips
  • Ragnar Fjørtoft
  • Arnt Eliassen
  • Joe Smagoinsky

9

slide-43
SLIDE 43

“Conference on Meteorology”

A “Conference on Meteorology” was arranged in the Insti- tute for Advanced Studies (IAS), Princeton on 29–30 Au- gust, 1946. Participants included:

  • Carl Gustav Rossby
  • Jule Charney
  • George Platzman
  • Norman Phillips
  • Ragnar Fjørtoft
  • Arnt Eliassen
  • Joe Smagoinsky
  • Phil Thompson

9

slide-44
SLIDE 44

The ENIAC

10

slide-45
SLIDE 45

The ENIAC

The ENIAC (Electronic Nu- merical Integrator and Com- puter) was the first multi- purpose programmable elec- tronic digital computer. It had:

  • 18,000 vacuum tubes
  • 70,000 resistors
  • 10,000 capacitors
  • 6,000 switches

11

slide-46
SLIDE 46

The ENIAC

The ENIAC (Electronic Nu- merical Integrator and Com- puter) was the first multi- purpose programmable elec- tronic digital computer. It had:

  • 18,000 vacuum tubes
  • 70,000 resistors
  • 10,000 capacitors
  • 6,000 switches

Power Consumption: 140 kWatts

11

slide-47
SLIDE 47

The ENIAC: Technical Details.

ENIAC was a decimal machine. No high-level language. Assembly language. Fixed-point arithmetic: −1 < x < +1. 10 registers, that is, Ten words of high-speed memory. Function Tables: 624 6-digit words of “ROM”, set on ten-pole rotary switches. “Peripheral Memory”: Punch-cards. Speed: FP multiply: 2ms (say, 500 Flops). Access to Function Tables: 1ms. Access to Punch-card equipment: You can imagine!

12

slide-48
SLIDE 48

Evolution of the Project:

  • Plan A: Integrate the Primitive Equations

Problems similar to Richardson’s would arise

13

slide-49
SLIDE 49

Evolution of the Project:

  • Plan A: Integrate the Primitive Equations

Problems similar to Richardson’s would arise

  • Plan B: Integrate baroclinic Q-G System

Too computationally demanding

13

slide-50
SLIDE 50

Evolution of the Project:

  • Plan A: Integrate the Primitive Equations

Problems similar to Richardson’s would arise

  • Plan B: Integrate baroclinic Q-G System

Too computationally demanding

  • Plan C: Solve barotropic vorticity equation

Very satisfactory initial results

13

slide-51
SLIDE 51

Charney, Fjørtoft, von Neumann

14

slide-52
SLIDE 52

Charney, et al., Tellus, 1950.

15

slide-53
SLIDE 53

Charney, et al., Tellus, 1950.

Absolute Vorticity

  • =

Relative Vorticity

  • +

Planetary Vorticity

  • η = ζ + f .

15

slide-54
SLIDE 54

Charney, et al., Tellus, 1950.

Absolute Vorticity

  • =

Relative Vorticity

  • +

Planetary Vorticity

  • η = ζ + f .

The atmosphere is treated as a single layer, and the flow is assumed to be nondivergent. Absolute vorticity is con- served:

d(ζ + f) dt = 0.

15

slide-55
SLIDE 55

Charney, et al., Tellus, 1950.

Absolute Vorticity

  • =

Relative Vorticity

  • +

Planetary Vorticity

  • η = ζ + f .

The atmosphere is treated as a single layer, and the flow is assumed to be nondivergent. Absolute vorticity is con- served:

d(ζ + f) dt = 0.

This equation looks deceptively simple. But it is nonlinear: ∂ζ ∂t + V · ∇(ζ + f) = 0 .

15

slide-56
SLIDE 56

Charney, et al., Tellus, 1950.

Absolute Vorticity

  • =

Relative Vorticity

  • +

Planetary Vorticity

  • η = ζ + f .

The atmosphere is treated as a single layer, and the flow is assumed to be nondivergent. Absolute vorticity is con- served:

d(ζ + f) dt = 0.

This equation looks deceptively simple. But it is nonlinear: ∂ζ ∂t + V · ∇(ζ + f) = 0 . In more detail: ∂ ∂t[∇2ψ − Fψ] +

  • ∂ψ

∂x ∂∇2ψ ∂y − ∂ψ ∂y ∂∇2ψ ∂x

  • + β∂ψ

∂x = 0

15

slide-57
SLIDE 57

Solution method for QGBVE

∂ζ ∂t = −J(ψ, ζ + f)

16

slide-58
SLIDE 58

Solution method for QGBVE

∂ζ ∂t = −J(ψ, ζ + f)

  • 1. Compute Jacobian

16

slide-59
SLIDE 59

Solution method for QGBVE

∂ζ ∂t = −J(ψ, ζ + f)

  • 1. Compute Jacobian
  • 2. Step forward (Leapfrog scheme)

16

slide-60
SLIDE 60

Solution method for QGBVE

∂ζ ∂t = −J(ψ, ζ + f)

  • 1. Compute Jacobian
  • 2. Step forward (Leapfrog scheme)
  • 3. Solve Poisson equation for ψ (Fourier expansion)

16

slide-61
SLIDE 61

Solution method for QGBVE

∂ζ ∂t = −J(ψ, ζ + f)

  • 1. Compute Jacobian
  • 2. Step forward (Leapfrog scheme)
  • 3. Solve Poisson equation for ψ (Fourier expansion)
  • 4. Go to (1).

16

slide-62
SLIDE 62

Solution method for QGBVE

∂ζ ∂t = −J(ψ, ζ + f)

  • 1. Compute Jacobian
  • 2. Step forward (Leapfrog scheme)
  • 3. Solve Poisson equation for ψ (Fourier expansion)
  • 4. Go to (1).
  • Timestep : ∆t = 1 hour (2 and 3 hours also tried)
  • Gridstep : ∆x = 750 km (approximately)
  • Gridsize : 18 x 15 = 270 points
  • Elapsed time for 24 hour forecast: About 24 hours.

16

slide-63
SLIDE 63

Solution method for QGBVE

∂ζ ∂t = −J(ψ, ζ + f)

  • 1. Compute Jacobian
  • 2. Step forward (Leapfrog scheme)
  • 3. Solve Poisson equation for ψ (Fourier expansion)
  • 4. Go to (1).
  • Timestep : ∆t = 1 hour (2 and 3 hours also tried)
  • Gridstep : ∆x = 750 km (approximately)
  • Gridsize : 18 x 15 = 270 points
  • Elapsed time for 24 hour forecast: About 24 hours.

Forecast involved punching about 25,000 cards. Most of the elapsed time was spent handling these.

16

slide-64
SLIDE 64

ENIAC Algorithm

17

slide-65
SLIDE 65

ENIAC: First Computer Forecast

18

slide-66
SLIDE 66

Richardson’s Reaction

19

slide-67
SLIDE 67

Richardson’s Reaction

  • “Allow me to congratulate you and your collaborators on

the remarkable progress which has been made in Prince- ton.

19

slide-68
SLIDE 68

Richardson’s Reaction

  • “Allow me to congratulate you and your collaborators on

the remarkable progress which has been made in Prince- ton.

  • “I have today made a tiny psychological experiment on

the diagrams in your Tellus paper (he had asked his wife Dorothy to compare the charts and to decide whether the initial analysis or the 24 hour forecast more closely resembled the verify- ing analysis).

19

slide-69
SLIDE 69

Richardson’s Reaction

  • “Allow me to congratulate you and your collaborators on

the remarkable progress which has been made in Prince- ton.

  • “I have today made a tiny psychological experiment on

the diagrams in your Tellus paper (he had asked his wife Dorothy to compare the charts and to decide whether the initial analysis or the 24 hour forecast more closely resembled the verify- ing analysis).

  • My wife’s opinions were that the ‘forecast’ has it on av-

erage, but only slightly.

19

slide-70
SLIDE 70

Richardson’s Reaction

  • “Allow me to congratulate you and your collaborators on

the remarkable progress which has been made in Prince- ton.

  • “I have today made a tiny psychological experiment on

the diagrams in your Tellus paper (he had asked his wife Dorothy to compare the charts and to decide whether the initial analysis or the 24 hour forecast more closely resembled the verify- ing analysis).

  • My wife’s opinions were that the ‘forecast’ has it on av-

erage, but only slightly.

  • “This is . . . an enormous scientific advance on the single,

and quite wrong, result in which . . . [Richardson (1922)] ended.”

19