Computation of the zonally-averaged circulation driven by heating due to radiation and turbulence

In a zonally-averaged dynamic model the steady state circulation of the middle atmosphere between 10 and 110 km altitude has been calculated for solstice. To investigate the combined effect of turbulent heat conduction and dissipation of eddy kinetic energy on the mean circulation, the dissipative heating has been parameterized in terms of the buoyancy term modified by a residual Richardson number. It is shown that turbulence will result in net heating of the mesopause region to be consistent with a zero mass flux through a pressure surface. It is also demonstrated that the combined effect of turbulent heat conduction and dissipation can modify the mean circulation remarkably if the Richardson number is made latitude-dependent.