Abstract. When conducting meso-micro scale coupled atmospheric simulations, it is crucial to ensure an adequate treatment of gray zone or terra incognita resolutions in which a large portion of the kinetic energy is naturally produced by the momentum balance equations in the model, while the remaining part still needs to be parameterized. In this work, we conduct three multi-day, real case, full-physics atmospheric simulations that are fully coupled from the meso to the micro scale and in which the only difference is the treatment of boundary layer physics at the gray zone domain. One simulation uses a well-established parameterization, another uses its scale-aware version previously modified to accommodate gray zone resolutions, and a final one uses no parameterization at all and assumes that the gray zone domain can be run in large-eddy simulation (LES) mode. The simulated fields are cross-compared, and further compared to measurements collected during the Prince Edward Island Wind Energy Experiment. Use of LES in the gray zone domain influences the flow fields in a manner that is robust to temporal averaging. The best predictions of vertical wind shear were found for the simulations in which the gray zone is parameterized, and the inclusion of a micro scale nest run in LES mode within the gray zone domains increased the model errors by producing overly homogeneous flow fields. The parameterized simulations also produced better agreement in terms of kinetic energy spectra at the two innermost simulation domains. In the gray zone domain, the energy decays as f⁻³ throughout most of the spectral range considered. In the micro scale domain, the same is only seen in the low-frequency end of the gray zone spectral range. In the high-frequency end, the energy decay follows a f⁻¹ slope. Outside the gray zone spectral range, the micro scale simulated spectra follow the expected f^−5/3 slope and produce good agreement with measurements.