Boundary Layer Parametrization in Numerical Weather Prediction Models

Numerical weather prediction (NWP) and climate models have shown to have a challenge to correctly simulate stable boundary layers and diurnal cycles. This aim of this study is to evaluate, describe and give suggestions for improvements of the descriptions of stable boundary layers in operational NWP models. Two papers are included. Paper I focuses on the description of the surface and the interactions between the surface and the boundary layer in COAMPSR, a regional NWP model. The soil parametrization showed to be of great importance to the structure of the boundary layer. Moreover, it showed also that a low frequency of radiation calculations caused a bias in received solar energy at the surface. In paper II, the focus is on the formulation of the turbulent transport in stable boundary layers. There, an implementation of a diffusion parametrization based on the amount of turbulent kinetic energy (TKE) is tested in a single column model (SCM) version of the global NWP model Integrated Forecast System (IFS). The TKE parametrization turned out to behave similarly as the currently operational diffusion parametrization in convective regimes and neutral regimes, but showed to be less diffusive in weakly stable and stable conditions. The formulations of diffusion also turned out to be very dependent on the length scale formulation. If the turbulence and the gradients of wind temperature and wind are weak, the magnitude of turbulence can enter an oscillating mode. This oscillation can be avoided with the use of a lower limit of the length scale.Det har visat sig att det är en stor utmaning för numeriska väderprognosmodeller (NWP-modeller) att simulera stabilt skiktade atmosfäriska gränsskikt och gränsskiktets dygnscykel på ett korrekt sätt. Syftet med denna studien är att utvärdera, beskriva och ge förslag på förbättringar av beskrivningen av gränsskiktet i NWP-modeller. Studien innehåller två artiklar. Den första fokuserar på beskrivningen av markytan och interaktionen mellan marken och gränsskiktet i den regionala NWP-modellen COAMPS R . Det visade sig att beskrivningen av markytan har en signifikant inverkan på gränsskiktets struktur. Det framkom också att strålningsberäkningarna endast görs en gång i timmen vilket bland annat orsakar en bias i inkommande solinstrålning vid markytan. Den andra artikeln fokuserar på beskrivningen av den turbulenta transporten i stabila skiktade gränsskikt. En implemenering av en diffusionsparametrisering som bygger på turbulent kinetisk energy (TKE) testas i en endimensionell version av NWP-modellen Integrated Forecast System (IFS), utvecklat vid European Center for Medium Range Weather Forecasts (ECMWF). Den TKE-baserade diffussionsparametriseringen är likvärdigt med den nuvaran de operationella parametriseringen i neutrala och konvektiva gränsskikt, menär mindre diffusivt i stabila gränsskikt. Diffusionens intensitet är beroende påden turbulenta längdskalan. Vidare kan turbulensen i TKE-formuleringen hamna i ett oscillerande läge om turbulensen är svag samtidigt som temperatur- och vindgradienten är kraftig. Denna oscillation kan förhindras om längdskalans minsta tillåtna värde begränsas