Einführung in die Meteorologie I

Clemens Simmer

Einführung

in die Meteorologie I

- Teil IV: Meteorologische

Zustandsvariable -

Gliederung der Vorlesung

0 Allgemeines I Einführung

II Zusammensetzung und Aufbau der Atmosphäre III Strahlung

IV Die atmosphärischen Zustandsvariablen V Thermodynamik der Atmosphäre

--- VI Dynamik der Atmosphäre

VII Atmosphärische Grenzschicht

VIII Synoptische Meteorologie

IV Die atmosphärischen Zustandsvariablen

IV.1 Luftdruck

IV.2 Windgeschwindigkeit und turbulente Transporte IV.3 Temperatur

IV.4 Feuchte

IV.2 Windgeschwindigkeit

1. Allgemeines

2. TurbulenteTransporte, Reynolds-Mittelung und logarithmisches Windprofil

3. Windmessung

• In der Atmosphäre finden durch die Luftbewegung ständig Transporte von Eigenschaften (z.B. Masse, Wasserdampf, Wärme) statt.

• Diese Transporte finden i. a. auf allen Skalen statt.

– Die Bewegung einzelner Moleküle transportieren,

– so wie auch kleine Wirbel beginnend mit Millimeter großen Auslenkungen in einer Zigarettenrauchfahne

– bis zu hundert Meter großen Wirbeln, die in Rauchfahnen von Schloten sichtbar werden,

– ebenso wie ganze Kubikkilometer Luft, die durch großskalige Druckgradienten bewegt werden.

• Legt man eine räumliche Skala fest, auf der man die Zustandsgrößen (als Mittel) betrachtet – wie z.B. in numerischen Modellen mit

Gitterweiten von Kilometern - so unterscheidet man:

– skalige Transporte, wenn sie z.B. durch die mittleren Größen auf dieser Skala (Mittelung z.B. über 1 km³, über 10 Minuten, etc.) bestimmt werden, und

– subskalige Transporte, die auf kleineren Skalen stattfinden. Diese subskaligen Transporte von Energie (Wärme, Wasserdampf, Masse, Impuls…) sind insbesondere an der Grenzfläche zwischen Untergrund (Erdboden, Wasserflächen) und der Atmosphäre bedeutend.

IV.2.2 Turbulente Transporte

Transportflussdichten

• In der Atmosphäre beziehen wir Transporte gern auf Einheitsflächen oder –querschnitte; d.h. wir möchten wissen, wieviel von einer

beliebigen Eigenschaft (Masse, Energie, Impuls) pro Sekunde durch eine Einheitsfläche hindurch geht.

• Dies kennen wir bereits von der Strahlung: Den Energiefluss durch Strahlung durch eine horizontale Einheitsfläche haben wir als

Strahlungsflussdichte F bezeichnet:

• Analog gibt es eine Massenflussdichte F

:

• ...und eine Impulsflussdichte F

(die wir schon als Druck kennen):

F = Strahlungsenergie

Quadratmeter ´ Sekunde mit éë ùû= F W

m

= J m

s

F

= Masse

Quadratmeter ´ Sekunde mit [ ] F

^{= kg} m

s

F

= Masse ´ Geschwindigkeit

Quadratmeter ´ Sekunde mit [ ] F

^{= kgms}

m

s = kg

ms

= Pa

Flussdichten durch Massentransporte

• Anders als beim Strahlungstransport sind die Transporte von Impuls,

Wärmeenergie (=kinetische Energie der Moleküle) oder Wasserdampf immer auch mit Massentransporten verbunden.

• Die Einheitenanalyse zeigt, dass wir massengebundene Transporte durch das Produkt von Massenflussdichte ρv, [ρv]=(kg/m

)(m/s)=kg/(m

s) und die ent-

sprechenden massenspezifischen Größe χ (d.h. Größe pro kg) schreiben können:

Wasserdampf: éë ùû = F

kg

m

s = kg kg

kg m

m

s Þ F

= r

v = r

r r v = q r v , éë ùû = q

spezifische Feuchte Impuls: [ ] F

^{= kg m s} _m

s = kg m s

kg

kg m

m s = m

s kg m

m

s Þ F

= v r v , [ ] v

mass. spez. Impuls Wärmeenergie: [ ] F

⁼ _m ^J

s = J

kg kg m

m

s Þ F

= e r v mit [ ] e ⁼

massenspezif. Energie e = c

T innere Energie mit

c

spezifische Wärme bei konst.

Volumen ( = 717J / (KgK ))

Allgemein: F

= cr v mit c = e,v , q,... massenspezifische Größen

Transporte: skalig und subskalig

���

Transporte/Flussdichten erfolgen skalig und/oder subskalig (Turbulenz).

Annahme: Die Luft links sei wärmer als rechts. Dann

können Wirbel an der Grenz- fläche Wärme von links nach rechts transportieren, auch wenn die mittlere Wind-

geschwindigkeit Null ist.

Wie beschreiben wir die skaligen und subskaligen Flussdichteanteile?

→ Reynolds-Mittelung

F

= cr v

Turbulenz und Reynolds Mittelung (1)

Mittelungsintervall Δt, Δx

Betrachte eine beliebige Eigenschaft ε als

Funktion der Zeit t oder des Ortes x und mittle diese über Intervalle in der Zeit Δt oder bezüglich des Ortes Δx.

 

x t ,

 '

Rechenregeln:  + d =  + d , d ' =  d ' = 0 , d = ?

 º lim

N®¥

1

N 

i=1

å

⁼ _Dt ^{1  (t ) dt}

Dt

ò Zeitmittelwert oder = 1

Dx  (x) dx

Dx

ò Raummittelwert

Turbulenz und Reynolds Mittelung (2)

• Produkte von meteorologischen Größen treten vielfach auf, z.B. im Advektionsterm bei der Eulerzerlegung.

• Für die Flussdichte einer Eigenschaft χ müssen wir ebenfalls über ein

Produkt - χρv - mitteln. Dazu schreiben wir alle drei Größen als Summe von Mittelwert und Abweichung, also , und damit

• Wir berücksichtigen weiter, dass die Luftdichte ρ im Vergleich zu anderen Größen in der Atmosphäre lokal kaum variiert; diese Approximation heißt Boussinesq-Approximation (Joseph Boussinesq, 1842 – 1929):

. Sie trägt der Tatsache Rechnung, dass die Luft wenig kompressibel ist (verhindert „störende“ Schallwell. bei num. Modellen).

 =  +  '

F

= r v c = ( r + r ')(v + v ')( c + c ')

11

Turbulenz und Reynolds Mittelung (3)

Was bedeutet der letzte Term?

t v‘(t)χ‘(t)

v(t)

t

t χ(t)

Δt

negativer turbulenter χ-Transport

Gemittelte Flussdichte r v c @ r v c + r v' c '

v ¢ c ¢ = 1

Dt v ¢ (t ) c ¢ (t ) dt

Dt

ò ⁼ _Dt ¹ ( ^{v (t ) - v} ) ( ^{c (t) - c} ) ^dt

Dt

ò

 0

¢

¢ c v

v

c '

v

c '

In Bodennähe verschwindet die mittlere Vertikalgeschwindigkeit (Massenerhaltung). Vertikaltransporte können dann nur noch turbulent erfolgen:

Besonderheit:

vertikaler Transport in Bodennähe

Beispiele:

• vertikaler turbulenter Horizontalimpulstransport τ (= Reibung)

• vertikaler turbulenter Wasserdampftransport E (= Verdunstung)

• vertikaler turbulenter Wärmetransport H

' ' da 0

w w w

r c r @ c =

t = r wv @ r w ' v' , E º r wq @ r w ' q' , H = r c

wT @ c

r w 'T ' mit v spezifischer Horizontalimpuls, q = r

r spezifische Feuchte und '

'

e Flussdicht

Gemittelte r v c @ r v c + r v c

• Nahe der Erdoberfläche (< 100m) nimmt der Horizontalwind von 0 am Boden stark mit der Höhe zu (bis zum geostrophischen Wind) und ist durch diese Windscherung oft turbulent (sog. Mischungsschicht oder Prandtl Schicht (Ludwig Prandtl, 1875-1953)).

• In dieser Schicht wird der Impuls der Luft durch den Impulsfluss τ - verursacht durch den vertikalen Austausch von Luftpaketen - zum Boden transportiert.

• Unter stationären Verhältnissen muss mit der Höhe konstant sein – denn Impulsflusszu- bzw. –abnahmen würden dort jeweils zu Geschwindigkeitsab- bzw. –zunahmen führen; die Kovarianz von w’ und v’ ist also höhenkonstant.

• Betrachten wir den vertikalen Austausch von Luftpaketen als Ursache des Impulsflusses, so muss – wenn der Wind zum Boden hin abnimmt - immer ein positives u’ mit einem negativen w’ zusammenhängen und umgekehrt, also τ  0.

• Wir führen eine Geschwindigkeit u* (Reibungsgeschwindigkeit) ein, welche die miteinander korrelierten (gleich- oder gegensinnig variierende) Geschwindig-

keitsabweichungen u’ und w’ quantifiziert (nach Messungen gilt etwa u*=1,25σ

mit σ

die Standardabweichung von v bzw. w), womit wir schreiben können:

Logarithmisches Windprofil (1)

t = r w ' v' = - r u _* ² = const

t = r w' v'

Logarithmisches Windprofil (2)

• Das vertikale Windprofil in der Prandtl Schicht erschließen wir uns wiederum durch eine Dimensionsanalyse.

• Die Änderung der Geschwindigkeit mit der Höhe ∂u/∂z sollte abhängen

– der Stärke des Impulsflusses und damit von u

– und von der Höhe über Grund z .

• Es gibt nur eine Möglichkeit, wie man nach Einheiten u

und z mit ∂u/∂z in Verbindung bringen kann:

(Theodore von Karman, 1881-1963)

• Integration liefert schließlich das logarithmische Windprofil:

¶u

¶z = u

k z mit k eine dimensionslose Konstante, die experimentell zu ~0,4 bestimmt wurde. Dies ist die von-Karman Konstante .

¶u = du = u

k

dz

z = u

k ^{dln z , u( z ) =} u

k ( ^{ln z - ln z}

) ^{mit u( z}

^{) = 0}

u z

Logarithmisches Windprofil (3)

• Das logarithmische Windprofil

gilt nur, wenn die Wärmeflüsse (Turbulenz durch Aufheizen oder Abkühlen des Bodens) klein sind – sog. neutrale Verhältnisse (Abb.).

• u

und z

lassen sich durch zwei Messungen des Horizontalwindes bestimmen.

• z

hängt mit der Bodenrauhigkeit zusammen und beträgt zwischen wenigen Millimetern (Grass) bis zu Metern (Wald).

Messungen über Grass (aus Blackadar, 1996) 15

u(z) = u

k ^ln

z

z

mit u

= t

r Reibungsgeschwindigkeit

und z

Rauhigkeitslänge

Struktur atmosphä- rische

Grenz-

schicht

Mittelung der prognostischen Grundgleichungen - Schließung

Die linken Seiten der prognostischen Grundgleichungen haben die Form:

• Die gemittelten (über Gitterelement, über Zeitintervall) prognostischen Grundgleichungen erhalten alle einen turbulenten Zusatzterm, die

Divergenz turbulenter Flüsse der die Anzahl der Unbekannten erhöht.

• Das Finden der nun notwendigen zusätzlichen Bestimmungs-

gleichungen wird Schließung (des Gleichungssystems) genannt.

d c

dt = ¶ c

¶t + v×  

( ) Ñ ^c

Zusammenfassung Flussdichten F

der massenspezifischen Eigenschaft c in Abhängigkeit von der Massenflussdichte r v

F

= cr v mit c = 1,v,e ,q ,...

Gemittelte Flussdichte

Reynolds Mittelung + Boussinesc Approximation r v c @ r v c + r v ' c '

vertikale Flussdichten am Boden:

Horizontalimpuls t = r wv @ r w' v' Wasserdampf E º r wq @ r w' q'

fühlbare Wärme H = r c

wT @ c

r w'T '

In der unteren Grenzschicht:

t = rw 'v ' = -ru

= const

¶u

¶z = u

k z aus Dimensionsgründen mit k ~0,4 von-Karman Konstante.

u(z) = u

k ^ln

z

z

mit z

Rauhigkeitslänge

Übungen zu IV.2.2

Folgende vertikale Windgeschwindigkeiten w, Temperaturen T und spezifische Feuchten q seien an einer zeitlich hochauflösenden Messstation (in 2 m

Höhe) über ein Zeitintervall von 10 Sekunden gemessen worden.

a) Berechne mittels R

oder Python die vertikalen turbulenten Flüsse der

(fühlbaren) Wärme H, des Wasserdampfes E und des Horizontalimpulses M (der Druck sei 1000 hPa):

b) Wie groß sind u

und z

?

c) Wie hoch ist die Windgeschwindigkeit in 10 m Höhe?

Zeit (s) 1 2 3 4 5 6 7 8 9 10

w (m/s) 0,5 0,4 0,1 -0,3 -0,6 -0,2 0,4 0,4 0,1 -0,3 v (m/s) 2,1 2,4 1,5 1,4 1,9 2,2 1,9 0,9 0,5 2,0 q (g/kg) 10,5 11,2 12,0 11,4 10,3 12,0 13,0 13,2 12,6 11,9

T (°C) 20,5 21,0 20,9 20,8 20,9 21,0 21,0 20,8 20,9 21,0