Abstract. Ensemble simulation of the general atmospheric circulation of the middle and upper atmosphere up to the lower thermosphere is performed using the 3-D nonlinear mechanistic numerical model MUAM. Residual mean meridional circulation (RMC) in terms of the Transformed Eulerian Mean is calculated for the boreal winter and changes in its vertical and meridional velocity components during different phases of simulated composite stratospheric warming (SW) events are studied. The simulation results show general decrease in RMC velocity components up to 30 % during and after SW in the mesosphere and lower thermosphere of the Northern Hemisphere. There are also increases in the downward and northward velocities at altitudes 50–70 km at the northern high latitudes. Associated changes in adiabatic heating/cooling rates can contribute to heating the stratosphere and cooling the mesosphere during the composite SW. The changes in the transport of conservative species (like ozone) during SWs are estimated. Weakening of ozone fluxes at the middle latitudes of the Northern Hemisphere may reach 30 % during SWs and 30–40 % after the events at the altitudes of stratospheric maximum of ozone concentration. Such statistically confident simulations of RMC reactions on SWs at altitudes up to the lower thermosphere are performed for the first time. The study of the residual meridional circulation is useful for effective analysis of wave impacts on the mean flow and for diagnostics of the transport of atmospheric gas species in the atmosphere.