Regionale Klimamodellierung Mit Methoden der regionalen Klimamodellierung werden globale Klimadaten für einzelne Regionen räumlich verfeinert.

Entwicklung

Das Weather Research and Forecasting WRF Model (Skamarock et al., 2005, 2008, 2021) wurde gemeinsam vom National Center for Atmospheric Research (NCAR), der National Oceanic and Atmospheric Administration (NOAA), der Air Force Weather Agency (AFWA), dem Naval Research Laboratory, der University of Oklahoma, und der Federal Aviation Administration (FAA) in den USA entwickelt. WRF ist ein mesoskaliges numerisches Wettervorhersagemodell, das so entwickelt wurde, dass es sowohl den Ansprüchen operationeller Wettervorhersage, regionaler Klimasimulationen als auch atmosphärischer Forschung gerecht wird. Die WRF Softwareumgebung unterstützt 2 dynamische Kerne: das Advanced Research WRF (ARW), das am NCAR entwickelt und gepflegt wird, und das WRF Nonhydrostatic Mesoscale Model (WRF-NMM) das von den National Centers for Environmental Prediction (NCEP) und NOAA entwickelt wird. WRF wird für Anwendungen mit Gitterauflösungen von 10 m bis 50 km angewendet, dabei lassen sich Gebiete mit hoher Auflösung in die Gebiete gröberer Auflösung einbetten und entweder gleichzeitig oder nacheinander rechnen. Für die Anwendung kann vom Anwender aus einer Vielzahl von Parametrisierungen und Projektionen eine geeignete Konfiguration zusammengestellt werden. In der Klima- und Wetterforschung kommt vor allem das WRF-ARW zum Einsatz, das im Weiteren basierend auf Skamarock et al. (2008 und 2021) beschrieben wird. Details finden sich auf der WRF-Homepage, die regelmäßig aktualisiert wird.
 WRF-Homepage

Atmosphärische Dynamik und Modellkonfiguration

Bei WRF handelt es sich um ein kompressibles nichthydrostatisches Modell. Für alle skalaren Größen gelten die Erhaltungsgleichungen. Es besteht die Möglichkeit, das WRF hydrostatisch anzuwenden. Das Modell beinhaltet die prognostischen Variablen Windgeschwindigkeit (u,v,w) sowie Störungen in potentieller Temperatur, Geopotential und Oberflächendruck.

Optional können turbulente kinetische Energie und eine Vielzahl skalarer Größen (z.B. Mischungsverhältnisse) ausgegeben werden. WRF wird auf einem Arakawa C-Gitter mit geländefolgenden hybriden Sigma-Druckkoordinaten betrieben (Park et al., 2013).

Die zeitliche Integration erfolgt mit einem Runge-Kutta-Schema 2. bzw. 3. Ordnung. WRF kann in verschiedenen Projektionen betrieben werden, z.B. Lambert conformal, Mercator, und Gradnetz (auch mit rotiertem Pol). Verschiedene horizontale Gitter können ineinander geschachtelt werden: mit Wechselwirkung in einer Richtung, in zwei Richtungen und mit sich bewegenden Gebieten (z.B. bei Hurrikan-Simulationen) (Abb. 1).

WRF bietet die Wahl aus einer Vielzahl von Optionen für die Modellphysik. So gibt es z.B. verschiedene Optionen für die Mikrophysik, Landoberflächenmodelle, Grenzschicht-Schemata und Konvektionsparametrisierungen unterschiedlicher Komplexität.

Das ganze WRF-System (Abb. 2) beinhaltet neben dem eigentlichen WRF-Modell auch ein Vorbereitungspaket (WRF-Preprocessing System (WPS)) als auch ein Paket zur Aufbereitung der Ergebnisse (WRF-Postprocessing Package (WPP)). Die Aufgabe des WPS ist es, meteorologische Antriebsdaten und terrestrische Daten aufzunehmen, Initialisierungsfelder zu erzeugen und alle Daten für den geplanten Modelllauf zu formatieren. Das WRF-Modell besteht aus digitalen Filtern, die Hintergrundrauschen aus den Eingabedaten entfernen, dem dynamischen Kern (WRF-ARW oder WRF-NMM) zur 3-dimensionalen Lösung der atmosphärischen Strömungsdynamik, den physikalischen Parameterisierungen (s.u.) subskaliger Prozesse und dem Austauschmodul für physikalische Größen zwischen Parameterisierungen und dynamischem Kern. Optional ist die zusätzliche Verwendung der variationellen Datenassimilation (WRF-VAR) oder eines Modells für Atmosphärenchemie (WRF-Chem).

Abb. 2: Flussdiagramm des WRF-Systems (nach Skamarock et al., 2008)

 Vergrößerung Abb.2

Physikalische Parametrisierungen

Subskalige Prozesse werden mit Hilfe von Parameterisierungen (z.B. Mikrophysik) oder Säulenmodellen (z.B. Landoberflächenmodelle) berechnet. Dabei sind die Parameterisierungen nicht unabhängig voneinander und tauschen Größen wie z.B. Temperatur, Feuchte und Niederschlag aus.
Die Mikrophysik beinhaltet die explizite Auflösung der Prozesse, die die Bildung von Wasserdampf, Wolken und Niederschlag beschreiben. Es kann hier aus einer Vielzahl von Modellen unterschiedlicher Komplexität gewählt werden, z.B. dem Kessler-Schema (Kessler, 1969), dem Thompsen et al. -Schema (Thompsen et al. 2004) und dem Morrison et al. -2-Momenten-Schema (Morrison et al., 2008). Auch Aerosole werden z.T. inzwischen berücksichtigt (z.B. Thompson und Eidhammer 2014).

Die Konvektionsparametrisierungen werden zur Erfassung subskaliger Effekte von konvektiven und/oder flachen Wolken benötigt (d.h. wenn die horizontale Auflösung > 4km ist). Diese Schemata sollen die vertikalen Flüsse durch nicht vom Modell aufgelöste Vertikalbewegungen darstellen. Für die Konvektionsparametrisierung stehen u.a. ein modifiziertes Kain-Fritsch-Schema (Kain, 2004), ein neues Tiedke-Schema (Zhang und Wang 2017) und das Grell- Freitas-Schema (Grell und Freitas 2014) zur Verfügung.

Die Prandtl-Schicht ist die unterste Schicht der Atmosphärischen Grenzschicht. Sie beträgt in etwa 10% der Atmosphärischen Grenzschicht. Die turbulenten Austauschkoeffizienten für Wärme, Feuchte und Impuls in der Prandtl-Schicht werden in WRF mit dem sog. Oberflächenschichtschema berechnet. Hier stehen Parametrisierungen zur Wahl, welche auf der Monin-Obukkhov-Theorie (Monin und Obukhov, 1954) und Stabilitätsfunktionen basieren. Häufig wird die Parametrisierung von Jiminez et al. (2012) verwendet. Allerdings ist die Prandtl-Schicht ist das Bindeglied zwischen der Landoberfläche und der planetarischen Grenzschicht (s.u.), d. h. die Wahl muss zum ausgewählten Grenzschichtschema und Landoberflächenmodell passen.

Die Wasser- und Energieflüsse an der Landoberfläche, d.h. im Boden, in der Vegetation und zwischen Boden, Vegetation und atmosphärischer Grenzschicht werden durch Landoberflächenmodelle berechnet, die mit WRF gekoppelt sind. Es besteht die Auswahl zwischen sehr unterschiedlichen Modellen, von denen das NOAH (Chen and Dudhia, 2001) mit der Möglichkeit einer Heterogenität der Vegetation in Gitterzellen (Li et al., 2013) und seit WRF Version 3.4 das NOAH-MP (Niu et al., 2011) die größte Komplexität aufweisen. Inzwischen beinhaltet Noah-MP im WRF auch zwei Getreidewachstumsmodelle (Liu et al. 2016, Warrach-Sagi et al. 2022).

Die Parametrisierung der planetarischen Grenzschicht (PBL) ist für die vertikalen subskaligen turbulenten Flüsse in der gesamten atmosphärischen Säule oberhalb der Prandtl-Schicht verantwortlich. Wird ein PBL- Schema in WRF aktiviert, so wird die vertikale Diffusion deaktiviert, da das PBL diese beinhaltet. Die Oberflächenflüsse werden vom Oberflächenschichtschema und dem Landoberflächenmodell zur Verfügung gestellt. Das PBL Schema bestimmt daraus die Flussprofile in der gut durchmischten Grenzschicht und der stabilen Schicht und berechnet so die atmosphärischen Tendenzen für Temperatur, Feuchte (inklusive Wolken) und horizontalen Impuls in der gesamten atmosphärischen Säule. Die PBL Schemata sind eindimensional und nehmen an, dass es eine klare Trennung zwischen subskaligen und vom Modell aufgelösten Wirbeln gibt. Zur Auswahl stehen verschiedene PBL Schemata, z.B. das Yonsei University (YSU) (Hong et al. 2006), das Mellor-Yamada-Janjic (MYJ) (Janjic 1990, 2002, Mellor and Yamada 1982) und das auf turbulenter kinetischer Energie basierende Mellor-Yamada-Nakanishi-Niino (MYNN) (Nakanishi und Niino 2006, 2009, Olson 2019) PBL Schema.

Mit Hilfe von Strahlungstransportmodellen wird die kurzwellige und langwellige Strahlung berechnet. Es kann zwischen ca. 5 Modellen gewählt werden, die unterschiedliche Spektralbänder unterscheiden und unterschiedlich komplex sind, z.B. das aus dem NCAR Community Atmosphere Model (CAM 3.0) stammende Schema für Klimasimulationen. Es kann die optischen Eigenschaften verschiedener Spurengase bzw. Aerosole berücksichtigen und wechselwirkt mit im Modell aufgelösten Wolken. Das CAM-Strahlungsschema ist in Collins et al. (2004) dokumentiert. Ein weiteres Strahlungstransportmodell für die langwellige Strahlung in WRF ist das Rapid Radiative Transfer Model for GCMs (RRTMG) (Iacono et al. 2008) das die Effekte von Wasserdampf, Ozon, CO2, Spurengasen und Wolken berücksichtigt.

Ausgabegrößen

Standardmäßig werden in WRF sämtliche berechnete atmosphärischen und Landoberflächenvariablen in NetCDF ausgegeben, eine Liste findet sich z.B. in der Dokumentation von WRF (Skamarock et al., 2008). Ein Umwandlungsprogramm erlaubt die Konvertierung in CF-konformes NetCDF Format, so dass die Ergebnisse mit den CDOs (Climate Data Operators) bearbeitet werden können.