Dynamik und Klimatologie von diabatischen Rossby-Wellen

Der Wintersturm Lothar zog am 26. Dezember 1999 über Europa und richtete in Frankreich, in Deutschland, in der Schweiz und in Österreich ungewöhnlich hohe Schäden an. Lothar entstand aus einer diabatischen Rossby Welle (DRW) und erreichte erst wenige Stunden vor dem europäischen Kontinent Orkanstärke. DRWs weisen ein interessantes atmosphärisches Strömungsmuster auf. Sie bestehen aus einer positiven PV-Anomalie in der unteren Troposphäre, die sich in einer Region mit starkem meridionalen Temperaturgradient befindet. Die positive PV-Anomalie löst eine zyklonale Strömung aus, dadurch wird östlich der PV-Anomalie warme Luft aus dem Süden herantransportiert. Während des Aufstieg der warmen Luft finden diabatische Prozesse statt, die zur Bildung einer neuen positiven PV-Anomalie in der unteren Troposphäre (PVA) führen. DRWs entstehen unabhängig von PV-Anomalien an der Tropopause. Falls sie jedoch mit ihnen in Wechselwirkung treten, kann - wie im Falle von Lothar - eine explosive Zyklogenese daraus resultieren. Im ersten Teil wird die Dynamik einer DRW am Beispiel des Wintersturms Lothar untersucht. Es wird insbesondere auf das Potential einer DRW zur explosiven Zyklogenese eingegangen. Im zweiten Teil wird das Aufretreten von DRWs in ECMWF-Vorhersagen untersucht. Es werden Unterschiede zwischen DRWs und anderen PV-Anomalien in der unteren Troposphäre hervorgehoben. Die Dynamik von DRWs wird mit Hilfe eines ECMWF-"Ensemble Prediction System" (EPS) des Wintersturms Lothar untersucht. Die 50 Modellläufe des EPS starten am 24. Dezember 1999 um 12 UTC und reichen bis zum 26. Dezember 1999 um 12 UTC. Nur 16 der 50 Modellläufe sagen einen ähnlich starken Sturm wie Lothar vorher. 10 Modellläufen sagen am 26. Dezember keine Zyklone mehr vorher. Die Ausprägung der baroklinen Zone, in der sich die DRW befindet, ist ausschlaggebend für die Intensität der DRW. Weitere wichtige Parameter sind der Feuchtegehalt der unteren Troposphäre und der latente Wärmefluss über dem Ozean. Diejenigen DRWs, die sich zu am 25. Dezember um 12 UTC näher als 400 km am Tropopausenjet befinden, entwickeln sich zu einer starken Zyklone. Alle anderen lösen sich auf oder bleiben schwache Zyklonen. Es ist schwierig, diabatische Prozesse in Wettervorhersagemodellen abzubilden, dementsprechend treten Schwierigkeiten bei der Vorhersage von PVAs auf. In den operationellen ECMWF-Vorhersagen von Juni 2004 bis Mai 2005 werden mit Hilfe eines Tracking- Algorithmus PVAs im Nordpazifik und Nordatlantik bestimmt und in fünf Kategorien eingeteilt. Die fünf Kategorien unterscheiden sich in ihrer Häufigkeit, ihrer Zugbahn und ihrer Gestalt. Im Nordpazifik entstehen doppelt so viele PVAs wie im Nordatlantik. Durchschnittlich werden im Winter weniger PVAs gefunden als im Sommer. Die Baroklinität und die Geschwindigkeit des Tropopausenjets ist in der Nähe von DRWs besonders hoch. Verglichen mit anderen PVAs weisen DRWs eine ähnliche Verteilung des reduzierten Bodendrucks auf. DRWs können in etwa gleich gut vorhergesagt werden wie andere PVAs.

Diskontinuierliche Galerkin-Verfahren für die operationelle Wettervorhersage

In dieser Arbeit wird ein neuer Dynamikkern entwickelt und in das bestehendernnumerische Wettervorhersagesystem COSMO integriert. Für die räumlichernDiskretisierung werden diskontinuierliche Galerkin-Verfahren (DG-Verfahren)rnverwendet, für die zeitliche Runge-Kutta-Verfahren. Hierdurch ist ein Verfahrenrnhoher Ordnung einfach zu realisieren und es sind lokale Erhaltungseigenschaftenrnder prognostischen Variablen gegeben. Der hier entwickelte Dynamikkern verwendetrngeländefolgende Koordinaten in Erhaltungsform für die Orographiemodellierung undrnkoppelt das DG-Verfahren mit einem Kessler-Schema für warmen Niederschlag. Dabeirnwird die Fallgeschwindigkeit des Regens, nicht wie üblich implizit imrnKessler-Schema diskretisiert, sondern explizit im Dynamikkern. Hierdurch sindrndie Zeitschritte der Parametrisierung für die Phasenumwandlung des Wassers undrnfür die Dynamik vollständig entkoppelt, wodurch auch sehr große Zeitschritte fürrndie Parametrisierung verwendet werden können. Die Kopplung ist sowohl fürrnOperatoraufteilung, als auch für Prozessaufteilung realisiert.rnrnAnhand idealisierter Testfälle werden die Konvergenz und die globalenrnErhaltungseigenschaften des neu entwickelten Dynamikkerns validiert. Die Massernwird bis auf Maschinengenauigkeit global erhalten. Mittels Bergüberströmungenrnwird die Orographiemodellierung validiert. Die verwendete Kombination ausrnDG-Verfahren und geländefolgenden Koordinaten ermöglicht die Behandlung vonrnsteileren Bergen, als dies mit dem auf Finite-Differenzenverfahren-basierendenrnDynamikkern von COSMO möglich ist. Es wird gezeigt, wann die vollernTensorproduktbasis und wann die Minimalbasis vorteilhaft ist. Die Größe desrnEinflusses auf das Simulationsergebnis der Verfahrensordnung, desrnParametrisierungszeitschritts und der Aufteilungsstrategie wirdrnuntersucht. Zuletzt wird gezeigt dass bei gleichem Zeitschritt die DG-Verfahrenrnaufgrund der besseren Skalierbarkeit in der Laufzeit konkurrenzfähig zurnFinite-Differenzenverfahren sind.

Stabilität und Transport in der planetaren Grenzschicht - Untersuchungen mit Radiosonden- und Spurengasmessungen während PARADE

Die vorliegende Arbeit untersucht die Struktur und Zusammensetzung der untersten Atmosphäre im Rahmen der PARADE-Messkampagne (PArticles and RAdicals: Diel observations of the impact of urban and biogenic Emissions) am Kleinen Feldberg in Deutschland im Spätsommer 2011. Dazu werden Messungen von meteorologischen Grundgrößen (Temperatur, Feuchte, Druck, Windgeschwindigkeit und -richtung) zusammen mit Radiosonden und flugzeuggetragenen Messungen von Spurengasen (Kohlenstoffmonoxid, -dioxid, Ozon und Partikelanzahlkonzentrationen) ausgewertet. Ziel ist es, mit diesen Daten, die thermodynamischen und dynamischen Eigenschaften und deren Einfluss auf die chemische Luftmassenzusammensetzung in der planetaren Grenzschicht zu bestimmen. Dazu werden die Radiosonden und Flugzeugmessungen mit Lagrangeschen Methoden kombiniert und es wird zwischen rein kinematischen Modellen (LAGRANTO und FLEXTRA) sowie sogenannten Partikeldispersionsmodellen (FLEXPART) unterschieden. Zum ersten Mal wurde im Rahmen dieser Arbeit dabei auch eine Version von FLEXPART-COSMO verwendet, die von den meteorologischen Analysefeldern des Deutschen Wetterdienstes angetrieben werden. Aus verschiedenen bekannten Methoden der Grenzschichthöhenbestimmung mit Radiosondenmessungen wird die Bulk-Richardson-Zahl-Methode als Referenzmethode verwendet, da sie eine etablierte Methode sowohl für Messungen und als auch Modellanalysen darstellt. Mit einer Toleranz von 125 m, kann zu 95 % mit mindestens drei anderen Methoden eine Übereinstimmung zu der ermittelten Grenzschichthöhe festgestellt werden, was die Qualität der Grenzschichthöhe bestätigt. Die Grenzschichthöhe variiert während der Messkampagne zwischen 0 und 2000 m über Grund, wobei eine hohe Grenzschicht nach dem Durchzug von Kaltfronten beobachtet wird, hingegen eine niedrige Grenzschicht unter Hochdruckeinfluss und damit verbundener Subsidenz bei windarmen Bedingungen im Warmsektor. Ein Vergleich zwischen den Grenzschichthöhen aus Radiosonden und aus Modellen (COSMO-DE, COSMO-EU, COSMO-7) zeigt nur geringe Unterschiede um -6 bis +12% während der Kampagne am Kleinen Feldberg. Es kann allerdings gezeigt werden, dass in größeren Simulationsgebieten systematische Unterschiede zwischen den Modellen (COSMO-7 und COSMO-EU) auftreten. Im Rahmen dieser Arbeit wird deutlich, dass die Bodenfeuchte, die in diesen beiden Modellen unterschiedlich initialisiert wird, zu verschiedenen Grenzschichthöhen führt. Die Folge sind systematische Unterschiede in der Luftmassenherkunft und insbesondere der Emissionssensitivität. Des Weiteren kann lokale Mischung zwischen der Grenzschicht und der freien Troposphäre bestimmt werden. Dies zeigt sich in der zeitlichen Änderung der Korrelationen zwischen CO2 und O3 aus den Flugzeugmessungen, und wird im Vergleich mit Rückwärtstrajektorien und Radiosondenprofilen bestärkt. Das Einmischen der Luftmassen in die Grenzschicht beeinflusst dabei die chemische Zusammensetzung in der Vertikalen und wahrscheinlich auch am Boden. Diese experimentelle Studie bestätigt die Relevanz der Einmischungsprozesse aus der freien Troposphäre und die Verwendbarkeit der Korrelationsmethode, um Austausch- und Einmischungsprozesse an dieser Grenzfläche zu bestimmen.