Ensemble simulation of the atmospheric general circulation at altitudes up to the lower thermosphere is performed using the 3-D nonlinear mechanistic numerical model MUAM. The residual mean meridional circulation (RMC), which is the superposition of the mean Eulerian and wave-induced eddy components, is calculated for the boreal winter. Changes in the vertical and meridional RMC velocity components are analysed at different stages of a simulated composite sudden stratospheric warming (SSW) event averaged over 19 model runs. The simulation results show a general decrease in RMC velocity components up to 30 % during and after SSW in the mesosphere and lower thermosphere of the Northern Hemisphere. There are also increases in the downward and northward velocities at altitudes of 20–50 km at the northern polar latitudes during SSW. Associated vertical transport and adiabatic heating can contribute to warming the stratosphere and downward shifting of the stratopause during the composite SSW. The residual mean and eddy mass fluxes are calculated for different SSW stages. It is shown that before the SSW, planetary wave activity creates wave-induced eddy circulation cells in the northern upper stratosphere, which are directed upwards at middle latitudes, northward at high latitudes and downwards near the North Pole. These cells increase heat transport and adiabatic heating in the polar region. During SSW, the region of upward eddy vertical velocity is shifted to high latitudes, but the velocity is still downward near the North Pole. After SSW, upward eddy-induced fluxes span the entire polar region, producing upward transport and adiabatic cooling of the stratosphere and providing the return of the stratopause to higher altitudes. The obtained statistically significant results on the evolution of RMC and eddy circulation at different SSW stages at altitudes up to the lower thermosphere can be useful for a better understanding the mechanisms of planetary wave impacts on the mean flow and for the diagnostics of the transport of conservative tracers in the atmosphere.

Using the general circulation model of the middle and upper atmosphere (MUAM), a number of numerical scenarios were implemented to study the impact of individual planetary waves (PWs) on the global atmospheric circulation, including zonal wind, temperature, and residual meridional circulation (RMC). The calculations were performed for the winter conditions of the Northern Hemisphere (January–February). We show the contribution to the formation of the dynamic and temperature regimes of the MUAM made by equatorial Kelvin waves propagating to the east, as well as atmospheric normal modes (NMs) with periods from 4 to 16 d. In particular, it is demonstrated that the impact of a 5 d PW and an ultra-fast Kelvin wave (UFKW) can change the speed of circulation flows by up to 6 % in the areas of their amplitude maxima. At the same time, this effect can be significantly enhanced in certain periods of time. The presented research results are important for a deeper understanding of the mechanisms of large-scale atmospheric interactions. Despite the obviousness and simplicity of the problem, such work has not been carried out yet.

Using a nonlinear model of the general circulation of the middle and upper atmosphere (MUAM), spatio-temporal structures of planetary waves (PWs) during boreal winter were studied. Modeling of global atmospheric circulation was performed for January–February. Despite the tropospheric PW sources shaped in the model, the phenomenon of 16-day PW excitation arise out of internal atmospheric sources in the southern lower thermosphere was discovered. To explain the observed phenomenon, a number of numerical experiments were carried out according to different scenarios with a selective turning (on/off) tropospheric sources of PW individual modes (having periods of 4–16 days) in the model. Also, the evolution of perturbed potential enstrophy for a 16-day PW, as well as the contribution of nonlinear interactions between individual PW to it, was studied. This made it possible for the first time to demonstrate explicitly the process of generation a secondary 16-day PW as a result of the nonlinear interconnection of 4-day and 5-day PWs. Graphical Abstract

The behavior of planetary waves and their influence on the global circulation of the Northern Hemisphere during different El Niño types is studied. Three sets of five boreal winters were chosen for each El Niño type: Modoki I and II and canonical El Niño. Based on data of the Japanese 55-year Reanalysis and the Modern-Era Retrospective Analysis for Research and Applications, the spatio-temporal structure of planetary waves and the residual mean circulation were analyzed. The results show that the canonical El Niño type is characterized by the weakest wave activity in March. It is also demonstrated that warming of the polar stratosphere, accompanied by maximizing wave activity and weakening of the zonal wind, may lead to earlier stratospheric polar vortex collapse and the early spring transition under Modoki I conditions. This study is the next step in understanding of the so-called long-range teleconnections, consisting of the propagation of a signal from the tropical El Niño Southern Oscillation source into the polar stratosphere.

Numerical simulation of the general atmospheric circulation was performed to estimate changes in amplitudes of the westward-travelling planetary waves (PWs) at altitudes from the Earth’s surface up to 300 km under different solar activity (SA) levels. The three-dimensional nonlinear mechanistic model of circulation of the middle and upper atmosphere “MUAM” was used. The atmospheric general circulation and PW amplitudes were calculated based on ensembles containing 16 model runs for conditions corresponding to low and high SA. PWs having periods of 4–10 days were considered. Comparison with the data of digital ionosondes showed that the MUAM model is capable of reproducing the considered PW modes at thermospheric heights. It is shown that under high SA conditions, PW amplitudes are significantly larger in the thermosphere and smaller in the middle atmosphere. The observed PW structures are influenced not only by changes in atmospheric refractive index and Eliassen–Palm flux but also by varying PW reflection in the lower thermosphere, which can change proportions of the wave energy transferred from the lower atmosphere to the upper layers and reflected downwards.

Numerical modeling of primary cosmic ray protons’ transport through the Earth’s atmosphere was performed for the energy spectra of solar energetic particle events (SEPs). Several events in the last three solar cycles were considered. A comparative analysis of the characteristics of coronal mass ejections and primary proton fluxes was carried out. The main results were quantitative estimates of the calculated atmospheric ionization count rate for a wide range of altitudes (from sea level up to 98 km). The difference in the influence of solar protons on the Earth’s atmosphere is considered for seven SEPs divided into three groups with similar solar sources (X-flare magnitude and coordinates) but with different characteristics of accelerated particle fluxes. The data obtained in this work are very important for future studies of radio wave propagation, atmospheric chemistry and climate change.