Impact of Noah-LSM Parameterizations on WRF Mesoscale Simulations: Case Study of Prevailing Summer Atmospheric Conditions over a Typical Semi-Arid Region in Eastern Spain

The current study evaluates the ability of the Weather Research and Forecasting Model (WRF) to forecast surface energy fluxes over a region in Eastern Spain. Focusing on the sensitivity of the model to Land Surface Model (LSM) parameterizations, we compare the simulations provided by the original Noah LSM and the Noah LSM with multiple physics options (Noah-MP). Furthermore, we assess the WRF sensitivity to different Noah-MP physics schemes, namely the calculation of canopy stomatal resistance (OPT_CRS), the soil moisture factor for stomatal resistance (OPT_BTR), and the surface layer drag coefficient (OPT_SFC). It has been found that these physics options strongly affect the energy partitioning at the land surface in short-time scale simulations. Aside from in situ observations, we use the Meteosat Second Generation (MSG) Spinning Enhanced Visible and Infrared Imager (SEVIRI) sensor to assess the Land Surface Temperature (LST) field simulated by WRF. Regarding multiple options in Noah-MP, WRF has been configured using three distinct soil moisture factors to control stomatal resistance (β factor) available in Noah-MP (Noah, CLM, and SSiB-types), two canopy stomatal resistance (Ball–Berry and Jarvis), and two options for surface layer drag coefficients (Monin–Obukhov and Chen97 scheme). Considering the β factor schemes, CLM and SSiB-type β factors simulate very low values of the latent heat flux while increasing the sensible heat flux. This result has been obtained independently of the canopy stomatal resistance scheme used. Additionally, the surface skin temperature simulated by Noah-MP is colder than that obtained by the original Noah LSM. This result is also highlighted when the simulated surface skin temperature is compared to the MSG-SEVIRI LST product. The largest differences between the satellite data and the mesoscale simulations are produced using the Noah-MP configurations run with the Monin–Obukhov parameterization for surface layer drag coefficients. In contrast, the Chen97 scheme shows larger surface skin temperatures than Monin–Obukhov, but at the expense of a decrease in the simulated sensible heat fluxes. In this regard, the ground heat flux and the net radiation play a key role in the simulation results.

“Gray Zone” Simulations using a Three-Dimensional Planetary Boundary Layer Parameterization in the Weather Research and Forecasting Model

Generating accurate weather forecasts of planetary boundary layer (PBL) properties is challenging in many geographical regions, oftentimes due to complex topography or horizontal variability in, for example, land characteristics. While recent advances in high-performance computing platforms have led to an increase in the spatial resolution of numerical weather prediction (NWP) models, the horizontal grid cell spacing (Δ x) of many regional-scale NWP models currently fall within or are beginning to approach the gray zone (i.e., Δ x ≈ 100 – 1000 m). At these grid cell spacings, three-dimensional (3D) effects are important, as the most energetic turbulent eddies are neither fully parameterized (as in traditional mesoscale simulations) nor fully resolved [as in traditional large eddy simulations (LES)]. In light of this modeling challenge, we have implemented a 3D PBL parameterization for high-resolution mesoscale simulations using the Weather Research and Forecasting model. The PBL scheme, which is based on the algebraic model developed by Mellor and Yamada, accounts for the 3D effects of turbulence by calculating explicitly the momentum, heat, and moisture flux divergences in addition to the turbulent kinetic energy. In this study, we present results from idealized simulations in the gray zone that illustrate the benefit of using a fully consistent turbulence closure framework under convective conditions. While the 3D PBL scheme reproduces the evolution of convective features more appropriately than the traditional 1D PBL scheme, we highlight the need to improve the turbulent length scale formulation.