An Investigation of the Grid Sensitivity in Large-Eddy Simulations of the Stable Boundary Layer

We revisit the longstanding problem of grid sensitivity, i.e., the lack of grid convergence in large-eddy simulations (LES) of the stable boundary layer. We use a comprehensive set of LES of the well-known Global Energy and Water Cycle Experiment Atmospheric Boundary Layer Study 1 (GABLS1) case with varying grid spacings between 12.5 m and 1 m to investigate several physical processes and numerical features that are possible causes of grid sensitivity. Our results demonstrate that there are two resolution regimes in which grid sensitivity manifests differently. We find that changing the numerical advection schemes and the subgrid-scale models alters the simulation results, but the options tested do not fully address the grid-sensitivity issue. Moreover, sensitivity runs suggest that the surface boundary condition and the interaction of the surface with the near-surface flow, as well as the mixing with the free atmosphere, are unlikely to be the causes of the observed grid sensitivity. One interesting finding is that the grid sensitivity in the fine grid-spacing regime (grid spacings <= 2 m) is closely related to the reduction in the energy content of large-scale turbulence, leading to less turbulence kinetic energy and hence lower boundary-layer heights. The present work demonstrates that there is still an urgent need to address this grid-sensitivity issue in order to perform reliable LES of the stable boundary layer.