Grid resolution dependency of land surface heterogeneity effects on boundary‐layer structure

• Published in: Quarterly journal of the Royal Meteorological Society 2022-01, Vol.148 (742), p.141-158
• Main Authors: Poll, S., Shrestha, P., Simmer, C.
• Format: Article
• Language: English
• Subjects: boundary layer; Boundary layers; boundary‐layer scheme; Eddies; grey zone; Height; Heterogeneity; Land cover; land surface heterogeneity; Mixing length; Modelling; numerical methods and NWP; Oscillations; Prediction models; Resolution; secondary circulations; Soil moisture; Surface properties; Thermal circulation; Turbulence; turbulence parametrization; Turbulent mixing; Weather forecasting; Wind; Wind speed
• ISSN: 0035-9009; 1477-870X
• DOI: 10.1002/qj.4196

Land surface heterogeneity exerts a strong control on atmospheric boundary‐layer (ABL) evolution by spatially varying the distribution and partitioning of surface energy fluxes and triggering secondary circulations. The representation of this physical process in numerical weather prediction (NWP) models is especially affected in the terra incognita as the model grid resolution approaches the length‐scale of the largest eddies in the boundary layer. We explore these effects for a mesoscale strip‐like land surface inhomogeneity in land cover, soil moisture or a superposition of both embedded in an elsewhere homogeneous landscape. The study is conducted with the numerical weather prediction model ICON (ICOsahedral Nonhydrostatic), using the default operational level 2.5 Mellor–Yamada turbulence closure (MY) and a large‐eddy simulation (LES) configuration as a benchmark. While simulations with the default ABL scheme approach the LES reference when refining the spatial grid towards finer resolution, the model generates artificial circulations leading to ABL height oscillations when the horizontal grid resolution (∆x) approaches the ABL height zi. The effect of these model‐induced circulations on the state of the boundary layer is even present with weak thermal heterogeneity (∆H) under low background wind speed (vx) but diminishes with increasing background wind speed. The tuning of the asymptotic turbulent mixing length‐scale l∞ in the operational ABL scheme helps in reducing the amplitude of the oscillations, thereby reducing the artificially induced circulations due to thermal heterogeneity which might act as unintentional trigger for clouds and precipitation. Based on the tuned synthetic model data from sensitivity runs, we propose a new parametrization for a 2‐D l∞ as a function of ∆H, zi/∆x and vx, which is otherwise held as a constant in the ABL scheme. This study analyses the effects of land surface inhomogeneities in the terra incognita range. While simulations with the Mellor–Yamada turbulence closure approach the LES reference simulations when refining the spatial grid towards the LES resolution, the model generates artificial circulations when the horizontal grid resolution approaches the length‐scale of the largest eddies in the boundary layer. These circulations are effectively attenuated by a proposed boundary‐layer parametrization tuning following the land surface heterogeneity.