Nonlocal transport and implied viscosity and diffusivity throughout the boundary layer in LES of the Southern Ocean with surface waves

Observations from the Southern Ocean Flux Station provide a wide range of wind, buoyancy, and wave (Stokes) forcing for large-eddy simulation (LES) of deep Southern Ocean boundary layers. Almost everywhere there is a nonzero angle Ω between the shear and the stress vectors. Also, with unstable forcing there is usually a depth where there is stable stratification, but zero buoyancy flux and often a number of depths above where there is positive flux, but neutral stratification. These features allow nonlocal transports of buoyancy and of momentum to be diagnosed, using either the Eulerian or Lagrangian shear. The resulting profiles of nonlocal diffusivity and viscosity are quite similar when scaled according to Monin–Obukhov similarity theory in the surface layer, provided the Eulerian shear is used. Therefore, a composite shape function is constructed that may be generally applicable. In contrast, the deeper boundary layer appears to be too decoupled from the Stokes component of the Lagrangian shear. The nonlocal transports can be dominant. The diagnosed across-shear momentum flux is entirely nonlocal and is highly negatively correlated with the across-shear component of the wind stress, just as nonlocal and surface buoyancy fluxes are related. Furthermore, in the convective limit the scaling coefficients become essentially identical, with some consistency with atmospheric experience. The nonlocal contribution to the along-shear momentum flux is proportional to (1 − cosΩ) and is always countergradient, but is unrelated to the aligned wind stress component.