Many issues pioneered by Jackson Herring deal with how nonlinear interactions shape atmospheric dynamics. In this context, we analyze new direct numerical simulations of rotating stratified flows with a large-scale forcing, which is either random or quasi-geostrophic (QG). Runs were performed at a moderate Reynolds number Re and up to 1646 turn-over times in one case. We found intermittent fluctuations of the vertical velocity w and temperature θ in a narrow domain of parameters as for decaying flows. Preliminary results indicate that parabolic relations between normalized third- and fourth-order moments of the buoyancy flux ∝wθ and of the energy dissipation emerge in this domain, including for passive and active scalars, with or without rotation. These are reminiscent of (but not identical to) previous findings for other variables and systems such as oceanic and atmospheric flows, climate re-analysis data, fusion plasmas, the Solar Wind, or galaxies. For QG forcing, sharp scaling transitions take place once the Ozmidov length scale ℓOz is resolved—ℓOz being the scale after which a turbulent Kolmogorov energy spectrum likely recovers at high Re.