Abstract
Numerical simulations of a convective boundary layer (CBL) are performed to investigate model behavior in the terra incognita, also known as the gray zone. The terra incognita of the CBL refers to a range of model grid spacing that is comparable to the size of the most energetic convective eddies, which are on the order of the boundary layer depth. Using the Rayleigh–Bénard thermal instability as reference, a set of idealized simulations is used to show that gray zone modeling is not only a numerical challenge, but also poses dynamical difficulties. When the grid spacing falls within the CBL gray zone, grid-dependent convection can occur. The size of the initial instability structures is set by the grid spacing rather than the natural state of the flow. This changes higher-order flow statistics and poses fundamental difficulties for gray zone modeling applications.