Transcript Zusammenfassung

VI.1.1 Divergenz und Massenerhaltung
Die Divergenz eines Windfeldes quantifiziert das Zusammen(Konvergenz, negative Divergenz) oder Auseinanderströmen
(Divergenz) der Luft.
u v w Bei Beschränkung auf die horizontalen
div v    v 
 
Windkomponenten wird der
x y z Zusammenhang zwischen Form des
Strömungsfeldes und Divergenz
u v
div vH    vH 

unmittelbar deutlich.
x y
y
t=0
t=t1
<0
x
>0
<0
1
Eulersche und Lagrangesche
Kontinuitätsgleichung
Euler‘sche Zerlegung für ρ:
Euler‘sche Kont‘gleichung:
Umrechnung:
 
d   

 v  
dt t



   v 
t
d
     v   v  
dt
 
 
  
  v      v  v  
aus Produktreg el anwenden
auf   v 
Lagrange‘sche Kont‘gleichung
 
d
    v
dt
2
Sonderfall: Inkompressibles Medium
• Ein Medium ist inkompressibel, wenn man es weder zusammenpressen
noch auseinander ziehen kann (z.B. näherungsweise Wasser). Dabei kann
es durchaus seine Form verändern oder im Inneren inhomogen sein
(veränderliche Dichte, z.B. eine Wasser-Öl-Mischung).
• Auch Luft kann für bestimmte Betrachtungen in guter Näherung als
inkompressibel angenommen werden. Dann gibt es z.B.
 keine Ausdehnung beim Aufsteigen,
 keine Schallwellen (Vereinfachung der Numerik bei Modellen).
• Man macht daher die Annahme der Inkompressibilität oft bei der
Beschreibung der Strömungsprozesse bei relativ geringen und langsamen
Vertikalauslenkungen, z.B. Strömungen in der Grenzschicht.
• Es gilt dann offensichtlich:
beachte aber:
 
d
 0   v  0
dt

 0 !!!
t
dicht
dünn
3
Rotation eines Vektorfeldes
- Vektor-Produkt des Nabla-Operators mit einem Vektor -

 
 
  v  rot v  



x

y

z

i
 u 
   
   v   x
 w  u
  
y
Ist die Vertikalgeschwindigkeit w=0 und
sind u und v vertikal konstant, dann gilt
offensichtlich:

j

y
v
 w v 


 
k  y z     xi
 u w   



    eta
z

z x
   zeta


w
v u
 



 x y 


.
 v u 
  
  v  k   k   
 x y 
Da die Luftströmung großskalig i.w. horizontal ist, hat ς
(Vorticity) eine große Bedeutung in der Meteorologie.
Offensichtlich ist die
Rotation aus der
Zeichenebene zum Beobachter gerichtet. Sie
wird als zyklonal
(Zyklone!) bezeichnet.
Die Rotation ist ein
x axialer Vektor.
 v u 
  k   v    
 x 4 y 
Vorticitygleichung
   f
mit  absolute Vorticity

relative Vorticity
f  2 sin  Coriolisparameter

 
Pol

• f ist die Rotation um die lokale Vertikale, die durch
die Erddrehung erzeugt wird (NH positiv, SH negativ).
• Ist der Drehsinn der Relativbewegung so wie der
Drehsinn der Erde, nennt man diesen zyklonal;
zyklonal heißt also auf der NH gegen Uhrzeigersinn
auf der SH im Uhrzeigersinn.
• Für die absolute Vorticity lässt sich unter barotropen
Verhältnissen (keine Vertikaländerungen des
Horizontalwindes und keine horizontalen
Dichteänderungen) ableiten (barotrope
Vorticitygleichung):
𝒅𝜼
𝒅𝒕
= −𝜼𝜵𝒉 ∙
  𝒗𝒉 .
 z  k   • Konvergenz erhöht die Vorticity und Divergenz
Äquator
reduziert sie.
• Bei zusätzlicher Divergenzfreiheit
führt
eine Nordwärtsbewegung eines Tief auf der NH
zu seiner Abschwächung
und
Südwärtsbewegung zu
einer
Verstärkung.
5
Natürliches Koordinatensystem
• Zur Untersuchung von Strömungen ist es oft nützlich anstatt des
starren und ortsfesten kartesischen Koordinatensystems ein
Koordinatensystem zu verwenden, das an die Strömung selbst
gebunden ist.
• Betrachtet man einen sehr kleinen Ausschnitt aus einer
beliebigen dreidimensionalen Strömung, so lässt sich dieser als
ein Teil eines Kreisbogens auffassen.
• Ein geeignetes Koordinatensystem wird dann festgelegt durch

drei Einheitsvektoren in Richtung
n0
- des Windrichtungsvektors (𝒔𝟎)
- des Vektors senkrecht dazu nach links in der
Strömungsebene (dieser ist dann parallel zur Richtung zum
hypothetischen Kreismittelpunkt) (𝒏𝟎)
- der Normalen auf der Ebene des Kreises (𝒌).
 


v v  
s0    , s0  n0  k
v v
  
s0 , n0 , k bilden ein Rechtssyst em

s0

n0

s0
6
Krümmungs- und Scherungsvorticity (b)
y
y
+
+
a
b
x


V

n
Scherungsvorticity
x

V
Rs
Krümmungsvorticity
7
VI.1.3 Stromlinien und Trajektorien
• Stromlinien sind Momentaufnahmen eines
Geschwindigkeitsfeldes. An jedem Punkt
bewegt sich zu diesem Zeitpunkt die Luft
parallel zu den Stromlinien.
• Trajektorien repräsentieren den Weg eines
Teilchens über eine Zeitspanne
8
 2

u  U  const , v  A cos
( x  ct ) 
 

Beispiel (1):
𝜆
Stromlinie für t=0
0.5
S2
y'
Trajektorie
S1
Trajektorie mit Start bei
0.0
S3
y 0  0 , x0  0 , t 0  0
-0.5
0.0
0.5
1.0
1.5
x'
2.0
Die Trajektorie hat hier eine größere Amplitude als die Stromlinie, da c
und U in die gleiche Richtung gehen, und entsprechend auch eine
längere Wellenlänge.
In der Abbildung wurden x und y mit λ normiert (→x‘, y‘) und U=A
und c=0,3U gesetzt.
9
Reibungskraft (4)
Berechnung der Nettokraft (=Nettoimpulsflussdichte x Fläche) in x-Richtung:
K R ,x   xz ( z0  z / 2)x y   xz ( z0  z / 2)x y

fR ,x 
 xz

x y z über  xz ( z0  z / 2)   xz ( z0 )  xz z / 2
z
z
V
K R ,x
m

 xz V
1  xz

Reibungsbeschleunigung nach x
z m  z
Laminare und turbulente Strömungen (Einsetzen von τ)
1  xz 1   u 
fR,x 



 z
 z  z 
laminar (   )
turbulent (   K ( z ))

1   u  1  
u 



 z  z   z 
z 
 2u
 2
z
mit ,  



1  
u 

K
(
z
)
 z 
z 

 
u 
K
(
z
)
z 
z 
dynamische, bzw. molekulare Viskosität (  1,5 · 10 5 m2 / s)
K
turbulenter Diffusionskoeffizient (
1 m2 / s)
10
Bewegungsgleichung für die Atmosphäre im
Inertialsystem


dv a
1 
1 
  p  g N    
dt


In der Bewegungsgleichung für das Inertialsystem treten Coriolis- und
Zentrifugalbeschleunigung nicht auf!
Ein brauchbares Inertialsystem ist ein in der Sonne verankertes
Koordinatensystem, das seine Achsen starr am Fixsternhimmels
ausrichtet.
11
Coriolisbeschleunigung
- formal (4) -



    
d v dv a d  


 r  2

v 


r


dt
dt 
dt


IV
V


I
II
III
I. Scheinbare Beschleunigung relativ zur Erdoberfläche
II. Beschleunigung im Inertialsystem (= Summe der angreifenden Kräfte)
III. Beschleunigung durch Änderung der Erdrotation (Herbsttag 0,05 s kürzer als
Sommertag, i.a. aber vernachlässigbar)
IV. Coriolisbeschleunigung
V. Zentrifugalbeschleunigung
12
Navier-Stokes-Gleichung (2)



  
  
dv v
1 

 (v  )v   p  gk  2  v  fR
dt
t

komponentenweise
du u
u
u
u
1 p

u
v w

 fv  2w cos   f R , x
dt t
x
y
z
 x
dv v
v
v
v
1 p
 u v w

 fu
 f R, y
dt t
x
y
z
 y
dw w
w
w
w
1 p

u
v
w

g
 2u cos   f R , z
dt
t
x
y
z
 z
gekoppelte nichtlineare Diff‘gleichungen 2. Ordnung
13
Skalenanalyse der horizontalen
Bewegungsgleichung
du
dt
1 p

 fv  2 cos  w  fR,x
 x
dv
dt
1 p

 fu
 y
U/T 1/ p/L
10-4
10-3
1 p
fv 
 x
1 p
fu  
 y
 fR,y
fU
fW
10-3
10-6
-
m/s2
...Coriolisbeschleunigung und
Druckgradientbeschleunigung
heben sich gegenseitig auf!
14
fv 
1 p
 x
Geostrophischer
Wind
 v 
0 u 
, fu  
1 p
 y
1
 
   
, f  u   f  0    v   fk  v h    h p

 0 
 1  0 
 
   
oder mittels Division durch -f und Multiplikation von links mit k 
0  v 
1
   
0


u

v

k  h p
h
   
f
 1  0 
   
T
p 3 p
geostrophischer Wind:

1  
vg 
k  h p
ρf
1 p
1 p
ug  
, vg 
ρf y
ρf x
FP,H
p 2 p
vg
p 1 p
FC,H
p
H
15
Fallunterscheidung und Bezeichnungen
Je nach wirkenden Kräften ergeben sich unterschiedliche Bewegungssysteme, die im
folgenden diskutiert werden.
Druckgradient
CoriolisBeschl.
Reibung
Zentrifugalbeschleunigung
geostrophischer
Wind
synoptische Systeme
Gradientwind
zyklostrophischer
Wind
Staubteufel
Trägheitskreis
Grenzschichtstrahlstrom
antitriptischer
Wind

s:
0 
Äquator
1 p
 fR,s
 s
 v h2
1 p
n:
 
 fv h
R
 n
16
Gradientwind – gekrümmte
Annahmen
Zusätzliche Annahme
• Stationarität
Stromlinien
- keine Reibung
• keine Bahnbeschleunigung

s:
1 p
0 
 fR,s
 s
 v h2
1 p
n:
 
 fv h
R
 n
s:
n:
0 
1 p
 s
vh2
1 p
 
 fvh
R
 n
R0
T

n


1  1 p vh2 
 vh   
 
f   n R 
  fv

g


R0
 p  0
s n
1  1 p vh2 
vh  
 

f   n R 
vh  vG  vg
 Stromlinien || Isobaren
H
p
 n  0
n

 1 p vh2 
1
s vh  
 

f   n R 
vh  vG  vg
Im Hoch wirkt Coriolis entgegen der
Im Tief kompensieren Coriolis und
Zentrifugalbeschleunigung gemeinsam den Zentrifugalbeschleunigung, daher höhere
17
Geschwindigkeit bei gleichem Druckgradient!
Druckgradient.
Diskussion - Besonderheit bei Hochs
Diskussion
2
fR
 fR  R p
vG  
   
2
 n
 2 
• Anormale Fälle werden auf der synoptischen Skala nicht beobachtet, da
Druckgradient die primäre Bewegungsursache ist.
• Anormale Fälle können nur auf sehr kleiner Skala durch Trägheitseffekte
auftreten (Staubteufel, Badewanne)
2
• Besonderheit des Hochs (R<0^∂p/∂n<0)
fR
R p
 fR 
vG 
   
(Wurzelargument muss positiv sein):
 2 
2
 n
f R
R p
 fR 
  
 
 n
 2 
 2
p
f2

 R
n
4
2
2

R p
 

 n

Druckgradient muss zum Zentrum abnehmen.
Hochs sind flach. Tiefs haben diese Beschränkung nicht.
18

s:
Zyklostrophischer
Wind
Zusätzliche Annahmen:
1 p
0 
 fR,s
 s
- keine Reibung
- keine Coriolisbeschleunigung (z. B. Äquatornähe,
kleiner Krümmungsradius, z.B. Staubteufel)
 v h2
1 p
n:
 
 fv h
R
 n
s:
0
n:
v h2

R

1 p
 s

 Stromlinien || Isobaren
1 p
 n
 R und
p
entgegengesetzte
n
Vorzeichen
Isobare und Stromlinie
vH
n
n
T
T
FP
Z
FP
vH
antizyklonal R < 0
p/ n > 0 , d.h. Tief
zyklonal R > 0
p/ n < 0 , d.h. Tief
Z
Welche Luftdruckdifferenz
herrscht in einem typischen
Staubteufel (Außenrand zu
Zentrum) mit 1 m
Durchmesser und einer
Windgeschwindigkeit am
Rand von 20 m/s?
19

s:
Zusätzliche Annahmen
- keine Reibung
- kein Druckgradient
Stromlinie
n
FC
Z
vH
0 
1 p
 fR,s
 s
Trägheitskreis n(1)
v
1 p
:
 
 fv
2
h
s:
0 
n:
v h2

R
 n
R
h
0
 fv h  v h  fR , also R  0 antizyklonal
vh
 f
R
also breitenabhängige Winkelgeschwindigkeit,

gleich der doppelten "Erdrotationsgeschwindigkeit"
0°
20°
43,3
60°
90°
f in 10-4s-1
0
0,5
1
1,26
1,46
Umlaufzeit, T=2π|R|/vh =2π/f, in Stunden
∞
35
17,5
13,8
12
|R| bei vh=10 m/s
∞
200
100
79
69
Als solche in der Atmosphäre kaum direkt beobachtet. Im Ozean
dagegen sind diese Trägheitsschwingungen durchaus häufig.
20
Trägheitskreis (2)Grenzschichtstrahlstrom
Der Trägheitskreis taucht
aber in der Form des sogenannten Grenzschichtstrahlstroms auf:
- Ausgangspunkt ist der subgeostrophische Wind in der Grenzschicht bedingt durch
Reibung an der Erdoberfläche. Stabilisiert sich die Luft durch Ausbleiben der Heizung
vom Boden in der Nacht, so reduziert sich die Reibung.
- Nehmen wir an, dass die Reibung plötzlich entfällt. Bei gegebenem Druck-gradient
wird dieser dann nicht durch die Coriolisbeschleunigung ausgeglichen – der Wind
beschleunigt zum Druckgefälle hin, wodurch die Coriolisbeschleunigung zunimmt.
- Der Wind beschleunigt, und zwar solange die Windrichtung eine Komponente zum
tiefen Druck hat, da die Resultierende von Coriolis- und Druckgradient-beschleunigung
eine Komponente in Richtung der Windrichtung hat.
- Ist der Wind parallel zu den Isobaren, so ist er stärker als der geostrophische Wind, er
ist supergeostrophisch.
- Die Coriolisbeschleunigung ist nun aber stärker als der Druckgradient, er dreht den
Windvektor zum hohen Druck. Die Resultierende „bremst“ dann den Wind. Dies geht
so lange bis die Coriolisbeschleunigung kleiner als der Druckgradient ist und wieder
eine Linksbeschleunigung wirkt….
Um dies quantitativ zu beschreiben müssen wir wieder zur
im x,y,z-System zurück.
Bewegungsgleichung
21
Richtung der Reibung unter Einfluss der
Coriolisbeschleunigung
Annahmen:
1. stationäre horizontale Strömung
2. gradlinige Isobaren (keine Zentrifuglabeschleunigung)

1

 h p  fk  vh  f R  0

fp
 fk vg
f R  fk   vh  vg  also f R  vh  vg  vag
Die Reibungsbeschleunigung steht senkrecht auf dem
ageostrophischen Wind - also nicht parallel zum Windvektor.
T

vh

fR

fC
H
22
Bausteine der modernen
Wettervorhersage
1. Online-Datensammlung
2. Datenassimilation -> aktueller Zustand der Atmosphäre
 Verschmelzen von Beobachtungen und „alter“ Vorhersage
 Methoden
-
Nudging
3-dimensionale variationelle Datenassimilation
4-dimensionale variationelle Datenassimilation
Ensemble-basierte Datenassimilation
…
3. Vorhersagelauf mit Modell
-
deterministische Vorhersage
Ensemble-Vorhersage
4. Interpretation der Modellausgabe
– Model Output Statistics (MOS)
23
Höhenkarten
• sind Topographien von isobaren Flächen, angegeben in geopotentiellen
Metern (gpm) h=(g/g0)z
– absolute Topographien, z.B. 850 hPa, 700 hPa, 500 hPa, 300 hPa, … enthalten
• h850, h700, … als Isolinien (sog. Isohypsen) in gpd(eka)m
• Isothermen
• relevante Messwerteintragungen (Radiosonden,
Flugzeuge, Satellit) als reduziertes Stationsmodell
– relative Topographien, z.B. h300 – h700
• geben Informationen über die mittlere virtuelle
Temperatur in den Schichten (niedrige Höhendifferenz
= kalt, große Höhendifferenz = warm, siehe später)
24
Der thermische Wind
- Zusammenfassung -
1
vg :
k  H p
f
vg
g

k   H Tv
z Tv f
vthermisch
g
 vg 
k   H Tv z
Tv f
𝑣𝑔,𝑢𝑛𝑡𝑒𝑛
Der thermische Wind
(= Änderung des
geostrophischen Windes mit
der Höhe durch einen
horizontalen
Temperaturgradienten) „weht“
um ein Kaltluftgebiet, wie der
geostrophische Wind um das
Tief.
𝑣𝑔,𝑜𝑏𝑒𝑛
T
𝑣𝑡ℎ𝑒𝑟𝑚𝑖𝑠𝑐ℎ
H
W
T
K
H
25
Barotrope und barokline Felder
•
barotrop:
Isoflächen von Druck und Temperatur sind parallel zueinander
  pTv  0


v g
 ln p
0
geostrophischer Wind mit der Höhe konstant
•
baroklin:
Isoflächen von Druck und Temperatur sind gegeneinander geneigt
  pTv  0


v g
 ln p
0
geostrophischer Wind ändert sich mit der Höhe
26
Gegenüberstellung von thermischen
und dynamischen Druckgebilden
H
T
T
warm
H
warm kalt warm
Thermische Tiefs
und Hochs
kalt
Divergenz
kalt
Konvergenz
Dynamische Tiefs
und Hochs werden
durch StrömungsT
strukturen (Divergenzen und Konvergenzen) in der Höhe
angetrieben.
Die resultierende Strömung am Boden verändert dann aber
wieder auch die Strömung in der Höhe.
H
27
Allgemeine
Vorticitygleichung
(2)




 u v   w v w u 
d 
1   p  p 
  2 


 f       f     


   x y   x z y z 
dt   


x

y

y

x
 


 
  
 



Tiltingterm
Solenoidterm
Divergenzterm
   w v w u  1   p  p 
d
  2 

    h  vh  


dt
 x z y z    x y y x 
Absolute Vorticity η (bzw. relative Vorticity ζ, wenn sich die Breite
nur wenig ändert) wird also erzeugt durch:
1. Horizontale Konvergenz
2. Kombination von horizontaler Änderung des Vertikalwindes
mit einer vertikalen Änderung des Horizontalwindes
3. Schneiden von Isolinien von Druck und Temperatur (Sonderfall
barokliner Verhältnisse).
28
Barotrope
Rossby-Wellen
(3)
𝒅𝜼/𝒅𝒕 = 𝒅𝝇/𝒅𝒕 + 𝒅𝒇/𝒅𝒕
= 𝒅𝝇/𝒅𝒕 + 𝒗 𝒅𝒇/𝒅𝒚 = 𝟎
Durch Breitenänderung initiierte
Drehbewegung der Strömung
N
λ
𝜍<0
S
𝜍=0
Initialstörung
𝜍>0
𝜍>0
𝜍=0
η=f
df/dt<0
df/dt>0
da
also
also
ς=0
𝜍=0
dς/dt>0
dς/dt<0
df/dt<0
also
dς/dt>0
29
Barotrope Rossby-Wellen – Ausbreitung
(2)
• Rossby-Wellen wandern also mit einer Geschwindigkeit, die von der
Strömungsgeschwindigkeit u0 und der Wellenlänge λ abhängt.
c  u0 
•
•
•
•

, mit k Wellenzah l (k  2π ,  Wellenlän ge)
λ
k²
d.h. die Wellen pflanzen sich mit Grundstromgeschwindigkeit u0 aus, aber
vermindert um β/k².
Je kürzer die Wellen, desto schneller wandern sie in Richtung des
Grundstroms (also nach Osten).
Bei 45° und λ > 7000 km Wellenlänge wandern Die Wellen bei einer
Grundstromgeschwindigkeit ū = 10 m/s nach Westen. Oft sind die langen
Wellen quasi-stationär.
Genauer: Alle Rossby-Wellen laufen bezogen auf ein mitdriftendes Partikel im
Grundstrom (also Grundstrom abziehen) nach Westen, und zwar je länger die
Welle, desto schneller (k~1/λ).
Wichtig: Rossby-Wellen erfordern neben der Erdrotation auch die
Kugelgestalt der Erde (β-Effekt)!
30
Barokline Rossby-Wellen - Schema (2)
 u v 
d
 


dt

x
y 

Aus dem
Divergenz/Konvergenzmuster ergibt sich
Aufsteigen auf der
Trogvorderseite und
Absteigen auf der
Trogrückseite.
(aus Roedel, 1994)
• Da die Geschwindigkeiten in der Höhe höher sind als darunter in
Bodennähe, überkompensieren die „Vergenzen“ in der Höhe die
„Vergenzen“ in Bodennähe.
• Daraus folgen Druckfall (Tief) auf der Trogvorderseite und
Druckanstieg (Hoch) auf der Trogrückseite.
31