AbstractDirect numerical simulations are used to investigate potential enstrophy in stratified turbulence with small Froude numbers, large Reynolds numbers, and buoyancy Reynolds numbers ($R{e}_{b} $) both smaller and larger than unity. We investigate the conditions under which the potential enstrophy, which is a quartic quantity in the flow variables, can be approximated by its quadratic terms, as is often done in geophysical fluid dynamics. We show that at large scales, the quadratic fraction of the potential enstrophy is determined by $R{e}_{b} $. The quadratic part dominates for small $R{e}_{b} $, i.e. in the viscously coupled regime of stratified turbulence, but not when $R{e}_{b} \gtrsim 1$. The breakdown of the quadratic approximation is consistent with the development of Kelvin–Helmholtz instabilities, which are frequently observed to grow on the layerwise structure of stratified turbulence when $R{e}_{b} $ is not too small.