Abstract:In the numerical gray zone of the convective boundary layer (CBL), the horizontal resolution is comparable to the size of organized convective circulation. As turbulence becomes partially resolved, gridscale variations of the subgrid-scale (SGS) turbulent fluxes become significant compared to the mean. Previously, such variations have often been ignored in scale-adaptive planetary boundary layer (PBL) schemes developed for the gray zone. This study investigates these variations with respect to height and resolution based on large-eddy simulations. It is found that SGS fluxes exhibit maximum variability at the center of the gray zone, where the resolved and the SGS mean fluxes are approximately equal. A simple analytical model is used to associate such characteristic variations to the nonlinear interactions of the dominant energy-containing mode of CBL turbulence. Examination of the horizontal distribution of the SGS fluxes reveals their preferential location over the updraft edges surrounding the core. A priori analysis further suggests the ability of a scale-similarity closure to reproduce the unique spatial patterns of the SGS fluxes at gray zone resolutions. Four scale-adaptive PBL schemes are evaluated focusing on their representations of the modeled SGS flux variability. Their shared shortcomings as a result of their gradient diffusion-based formulation are exposed. This study suggests that a mixed model consisting of a scale-adaptive PBL scheme to represent the mean, and a scale-similarity component to account for gridscale variability to be advantageous for the gray zone. Significance StatementAs Gresho and Lee's famous quote on numerical schemes goes, "Don't suppress the wiggles. They're telling you something!" In the numerical gray zone, where the grid spacing is comparable to the characteristic length scale of the flow, turbulence is partially resolved. The "wiggles" (or gridscale variability) reflecting the resolved heterogeneity of turbulent fluxes is an outstanding feature of the gray zone. However, they are often overlooked as error bars to the mean. This study uncovers the significance of these error bars, and characterizes and explains their variations with height and grid spacing. A suitable model that captures such variations is investigated to help build better planetary boundary layer schemes.

On the Departure from Monin–Obukhov Surface Similarity and Transition to the Convective Mixed Layer

On the departure from Monin-Obukhov surface similarity and transition to the convective mixed layer

A Novel Similarity Approach for Describing the Bulk Shear in the Atmospheric Surface Layer

The mean wind and potential temperature flux profiles in convective boundary layers

Moving Beyond Monin–Obukhov Similarity Theory in Modelling Wind-Speed Profiles in the Lower Atmospheric Boundary Layer under Stable Stratification

Stable Boundary Layers with Subsidence: Scaling and Similarity of the Steady State

The Universality of the Normalized Vertical Velocity Variance in Contrast to the Horizontal Velocity Variance in the Convective Boundary Layer

Analytical model coupling Ekman and surface layer structure in atmospheric boundary layer flows

An Improved Surface Boundary Condition for Large-Eddy Simulations Based on Monin–Obukhov Similarity Theory: Evaluation and Consequences for Grid Convergence in Neutral and Stable Conditions

Stable Boundary Layers with Subsidence: Scaling and Similarity of the Truly Steady State

Connection between mass flux transport and eddy diffusivity in convective atmospheric boundary layers

The local similarity theory, including the "spectral" Prandtl mixing length. Dissipation of energies and structural parameters in the atmospheric CBL

A Local Similarity Function for Katabatic Flows Derived from Field Observations over Steep‐ and Shallow‐Angled Slopes

Variations of Subgrid-Scale Turbulent Fluxes in the Dry Convective Boundary Layer at Gray Zone Resolutions

Large Eddies Regulate Turbulent Flux Gradients in Coupled Stable Boundary Layers

Near‐surface wind profiles from numerical model predictions. Part I: Algorithms and comparisons with wind profile based on Monin–Obukhov similarity theory

A Cautionary Note on the Use of Monin–Obukhov Similarity Theory in Very High-Resolution Large-Eddy Simulations

The Effect of Coupling Between CLUBB Turbulence Scheme and Surface Momentum Flux on Global Wind Simulations

Distinct Turbulence Structures in Stably Stratified Boundary Layers with Weak and Strong Surface Shear

Quasi-stationary State of Convective Boundary Layer over Heterogeneous Flat Surfaces