Turbulent winds and temperature fronts in large-eddy simulations of the stable atmospheric boundary layer

The nighttime high-latitude stably stratified atmospheric boundary layer (SBL) is computationally simulated using high-Reynolds number large-eddy simulation on meshes varying from 200³ to 1024³ over 9 physical hours for surface cooling rates C_r = [0.25, 1] K h^-1. Continuous weakly stratified turbulence is maintained for this range of cooling, and the SBL splits into two regions depending on the location of the low-level jet (LLJ) and . Above the LLJ, turbulence is very weak and the gradient Richardson number is nearly constant: . Below the LLJ, small scales are dynamically important as the shear and buoyancy frequencies vary with mesh resolution. The heights of the SBL and Ri noticeably decrease as the mesh is varied from 200³ to 1024³. Vertical profiles of the Ozmidov scale L_o show its rapid decrease with increasing C_r with L_o < 2 over a large fraction of the SBL for high cooling. Flow visualization identifies ubiquitous warm-cool temperature fronts populating the SBL. The fronts span a large vertical extent, tilt forward more so as the surface cooling increases, and propagate coherently. In a height-time reference frame, an instantaneous vertical profile of temperature appears intermittent, exhibiting a staircase pattern with increasing distance from the surface. Observations from CASES-99 also display these features. Conditional sampling based on linear stochastic estimation is used to identify coherent structures. Vortical structures are found upstream and downstream of a temperature front, similar to those in neutrally stratified boundary layers, and their dynamics are central to the front formation.