Evaluating Turbulence Parameterizations at Gray Zone Resolutions for the Ocean Surface Boundary Layer

Turbulent mixing in ocean boundary layers is often fully parameterized as a subgrid-scale process in realistic ocean simulations. However, recent submesoscale modeling studies have advanced to a horizontal grid spacing of O(10 m) that is comparable to, or even smaller than, the typical depth of the turbulent surface boundary layer. Meanwhile, efforts toward realistic large-eddy simulations (LES) nested within regional models require subdomains with similar grid spacings, where turbulent eddies are partially resolved in the mixed layer. The range of intermediate gird resolution, often known as the "gray zone", presents challenges for model configuration and analysis, including uncertainties regarding the behavior of common turbulence closures outside of their ideal use cases. In this study, we evaluate three common configurations for subgrid turbulence—k-ε, Smagorinsky, and an implicit no-closure method—in the gray zone resolutions for the ocean surface mixed layer. Results indicate the k-ε closure shows less sensitivity to grid spacing, producing accurate mean mixed-layer profiles even with partially resolved turbulence. However, it overly damps turbulent motions, significantly reducing small-scale variability that could otherwise be captured. The Smagorinsky closure and the implicit method, in contrast, exhibit higher sensitivity to grid spacing, initially performing poorly but converging toward baseline solutions at finer grids. Our findings provide guidance for submesoscale and turbulent-scale modeling, recommending Smagorinsky or implicit methods for nesting scenarios at turbulence-permitting resolutions. The k-ε closure is more suitable for high-resolution models primarily focused on accurate mean-state representations rather than explicitly resolving detailed three-dimensional turbulence.