This study summarizes the revision performed on the surface layer formulation of the Weather Research and Forecasting (WRF) model. A first set of modifications are introduced to provide more suitable similarity functions to simulate the surface layer evolution under strong stable/unstable conditions. A second set of changes are incorporated to reduce or suppress the limits that are imposed on certain variables in order to avoid undesired effects (e.g., a lower limit in u * ). The changes introduced lead to a more consistent surface layer formulation that covers the full range of atmospheric stabilities. The turbulent fluxes are more (less) efficient during the day (night) in the revised scheme and produce a sharper afternoon transition that shows the largest impacts in the planetary boundary layer meteorological variables. The most important impacts in the near-surface diagnostic variables are analyzed and compared with observations from a mesoscale network.

The Weather Research and Forecasting (WRF) model presents a high surface wind speed bias over plains and valleys that constitutes a limitation for the increasing use of the model for several applications. This study attempts to correct for this bias by parameterizing the effects that the unresolved topographic features exert over the momentum flux. The proposed parameterization is based on the concept of a momentum sink term and makes use of the standard deviation of the subgrid-scale orography as well as the Laplacian of the topographic field. Both the drag generated by the unresolved terrain and the possibility of an increase in the speed of the flow over the mountains and hills, where it is herein shown that WRF presents a low wind speed bias, are considered in the scheme. The surface wind simulation over a complex-terrain region that is located in the northeast of the Iberian Peninsula is improved with the inclusion of the new parameterization. In particular, the underestimation of the wind speed spatial variability resulting from the mentioned biases is corrected. The importance of selecting appropriate grid points to compare with observations is also examined. The wind speed from the nearest grid point is not always the most appropriate one for this comparison, nearby ones being more representative. The new scheme not only improves the climatological winds but also the intradiurnal variations at the mountains, over which the default WRF shows limitations in reproducing the observed wind behavior. Some advantages of the proposed formulation for wind-resource evaluation are also discussed.

Mesoscale numerical weather prediction models using fine-grid [O(1) km] meshes for weather forecasting, environmental assessment, and other applications capture aspects of larger-than-grid-mesh size, convectively induced secondary circulations (CISCs) such as cells and rolls that occur in the convective planetary boundary layer (PBL). However, 1-km grid spacing is too large for the simulation of the interaction of CISCs with smaller-scale turbulence. The existence of CISCs also violates the neglect of horizontal gradients of turbulent quantities in current PBL schemes. Both aspects—poorly resolved CISCs and a violation of the assumptions behind PBL schemes—are examples of what occurs in Wyngaard’s “terra incognita,” where horizontal grid spacing is comparable to the scale of the simulated motions. Thus, model CISCs (M-CISCs) cannot be simulated reliably. This paper describes how the superadiabatic layer in the lower convective PBL together with increased horizontal resolution allow the critical Rayleigh number to be exceeded and thus allow generation of M-CISCs like those in nature; and how the M-CISCs eventually neutralize the virtual temperature stratification, lowering the Rayleigh number and stopping their growth. Two options for removing M-CISCs while retaining their fluxes are 1) introducing nonlocal closure schemes for more effective removal of heat from the surface and 2) restricting the effective Rayleigh number to remain subcritical. It is demonstrated that CISCs are correctly handled by large-eddy simulation (LES) and thus may provide a way to improve representation of them or their effects. For some applications, it may suffice to allow M-CISCs to develop, but account for their shortcomings during interpretation.

T he ChAllenge. Climate and weather forecasting applications share a common ancestry and build on the same physical principles. Nevertheless, climate research and numerical weather prediction (NWP) are commonly seen as different disciplines. The emerging concept of "seamless prediction" forges weather forecasting and climate change studies into a single framework. At the same