Regime-dependent turbulence length scale formulation for NWP models based on turbulence kinetic energy, shear and stratification

This paper presents a new Turbulence Length Scale (TLS) formulation based on turbulence kinetic energy, shear, and stratification. Starting from the well-known parcel-based method, several updates are introduced to obtain (i) a smooth and stability-dependent near-surface transition towards the similarity theory-based TLS, (ii) a regime-dependent TLS decrease near the Planetary Boundary Layer (PBL) top and (iii) a non-vanishing (constant) upper-air TLS to support the jet-stream-related turbulence. The new formulation’s TLS profile better matches its LES-diagnosed counterpart in idealized stratocumulus and very stable PBL cases, thus improving the PBL structure and decay rate of turbulence in aloft layers. Simultaneously, the performance in convective cases is nearly neutral to the reference formulation. Aiming its implementation in the 3D operational numerical weather prediction model, the new TLS formulation is evaluated in inversion and mesoscale convective system cases. Thereby, a spatial distribution of cloudiness and precipitation better corresponds to products derived from meteorological satellites and radars. Furthermore, the comparison with atmospheric soundings and tower measurements points to improved wind speed and temperature profiles within the PBL and turbulent fluxes in the surface layer. The preliminary objective verification suggests (i) the essential role of so-called k scaling, (ii) a nearly neutral upper-air performance with the reference formulation, (iii) a decrease in RMSE of cloudiness and random error of 2-m temperature (1.5–2.0%) and (iv) a potential to improve the 10-m wind speed prediction for categories above the 90th percentile (3.0–6.5%). Finally, we discuss possible further upgrades of the turbulence scheme used in our model.