Original paper

Evaluation of the COSMO-SC turbulence scheme in a shear-driven stable boundary layer

Buzzi, Matteo; Rotach, Mathias W.; Holtslag, Matthias; Holtslag, Albert A.M.

Meteorologische Zeitschrift Vol. 20 No. 3 (2011), p. 335 - 350

published: Jun 1, 2011

DOI: 10.1127/0941-2948/2011/0050

BibTeX file

ArtNo. ESP025012003008, Price: 29.00 €

Download preview PDF Buy as PDF

Abstract

The performance of the COSMO single column turbulence scheme (a TKE scheme with a 1.5 order turbulence closure at the hierarchy level 2.5 following Mellor and Yamada) is investigated in the framework of the first GABLS intercomparison case. This is an idealized shear-driven stable boundary layer case with no advection. Overall the COSMO model performs reasonably well compared to the other participating models and the reference Large Eddy Simulations. However, the modification of some model parameters, together with the prescribed high vertical resolution, reveals a problem of numerical stability in the turbulence scheme: for the investigated shear-driven stable boundary layer the vertical diffusivities show unrealistic oscillations. This model deficiency, which has previously been described in literature, is explored in quite substantial detail and possible solutions are evaluated. It is found that under the given conditions the numerical description of the vertical wind gradients is crucial for the stability of the turbulence scheme. It is shown that for the determination of vertical gradients information from grid points beyond the immediately neighboring model levels must be incorporated - as it is common practice in the treatment of horizontal gradients - in order to obtain a numerically stable turbulence scheme. As a first approach vertical wind gradients are filtered using a 5-point filter prior to the evaluation of the stability functions. This approach yields to the overall best performance among all those tested and found in literature. The simulations additionally show that the use of a too high minimum diffusion coefficient (which is introduced in the model in order to avoid too low mixing) leads to losing important structures of the planetary boundary layer, such as the low level jet or a near-surface temperature inversion.

Kurzfassung

Die Vorhersagegüte des Turbulenzschemas im COSMO-Modell (ein TKE-basiertes Schema mit einer Turbulenzschliessung der Ordnung 1.5 auf dem Hierarchie-Niveau 2.5 gemäß Mellor und Yamada) wird anhand der ersten GABLS -Fallstudie zur stabilen Grenzschicht untersucht. Im Vergleich mit den anderen teilnehmenden Modellen und den Large Eddy Simulationen als Referenz erreicht das Modell recht gute Ergebnisse. Allerdings führt die Modifikation einiger Modellparameter, kombiniert mit der sehr hohen vertikalen Auflösung, die das GABLS Experiment vorschreibt, zur Ausbildung einer numerischen Instabilität. In der durch Windscherung bestimmten Grenzschicht zeigen die vertikalen Diffusionskoeffizienten unrealistische Oszillationen. Dieses problematische Modellverhalten, das schon früher in der Literatur beschrieben wurde, wird im Detail untersucht und verschiedene Lösungsmöglichkeiten werden diskutiert und ausgetestet. Es wird gezeigt, dass - analog zur numerischen Behandlung horizontaler Gradienten - bei der Bestimmung der vertikalen Gradienten Informationen aus benachbarten Modellschichten einbezogen werden müssen, um numerische Stabilität des Turbulenzschemas für die untersuchten Grenzschichtbedingungen zu gewährleisten. Von den untersuchten Lösungsmöglichkeiten zeigt die Filterung der vertikalen Gradienten der Windgeschwindigkeit mit einem 5-Punkte-Filter die besten Ergebnisse. Die Simulationen zeigen außerdem, dass die Limitierung des Diffusionskoeffizienten, wie sie in vielen operationellen Modellen zur Verhinderung von Entkoppelung unter sehr stabilen Bedingungen verwendet wird, zu einer schlechten Simulation wichtiger Charakteristika der stabilen Grenzschicht, wie des “low level jet” oder einer ausgeprägten Temperaturinversion nahe der Oberfläche, führt.