Modelling of Earth thermosphere-ionosphere coupled dynamics
D. V. Kulyamin1,2, V. P. Dymnikov1,2, P. A. Ostanin3
1Marchuk Institute of Numerical Mathematics RAS, Moscow
2Fedorov Institute of Applied Geophysics, Moscow
2Moscow Institute of Physics and Technology
Email: kulyamind@mail.ru
DOI 10.24412/cl-35065-2021-1-01-33
One of the central directions in Earth system models development is inclusion of upper atmosphere and
ionospheric models. Such models create possibilities not only to study the atmospheric upper and lower layers
interaction but also to apply them to solving the practical problems, such as radio waves propagation through
ionosphere. The difficulty in constructing such models lies in variations of atmospheric characteristics and
specificity of the low-temperature plasma dynamics in the Earth's magnetic field, in addition it is necessary to
account a large amount of small gas impurities responsible for both radiation heat fluxes and ionization pro-
cesses, etc.
In our talk we discuss the problem of thermosphere-ionosphere coupled dynamics and present a new
coupled model of Earth's thermosphere and ionosphere global dynamics (altitudes 90-500 km). The model is
based on a three-dimensional thermospheric general circulation model and a dynamical model of the iono-
spheric F region, which takes into account plasma-chemical processes, ambipolar diffusion, and advective ion
transport due to neutral wind and electromagnetic drift. Realistic upper atmospheric processes main charac-
teristics in model is shown and quantitative estimates of thermosphere � ionosphere interaction are obtained
based on presented coupled model.
Eddy transport in a stably stratified atmospheric boundary layer over rough surface: a study with of nonlocal
model turbulence
L. I. Kurbatskaya
Institute of Computational Mathematics and Mathematical Geophysics SB RAS
Email: L.Kurbatskaya@ommgp.sscc.ru
DOI 10.24412/cl-35065-2021-1-01-34
The features of the atmospheric boundary layer (ABL) depend it is on atmospheric stratification (effects of
buoyancy) and the dominating mechanism of turbulence generation. The boundary layer becomes stably strat-
ified when the underlying surface is colder than the air. In these conditions turbulence can be generated by
shear and destroyed under the influence of negative buoyancy and viscosity. Because of the different shear
and buoyancy effects, turbulence in the stable boundary layer (SBL) is inhibited in comparison to the neutral
and convective boundary layers (CBL). In this study, the turbulent momentum and heat fluxes are calculated
with explicit algebraic models obtained with use of symbolic algebra from the transport equations for the
momentum and heat fluxes in the approximation of weakly equilibrium turbulence [1, 2]. The nonlocality of
the mechanism of turbulent momentum and heat transfer in the atmospheric boundary layer over a rough
surface manifests itself in the form of limited regions of countergradient of momentum and heat transfer,
which are diagnosed from analysis of balance terms in the transport equation for the variance of temperature
fluctuations and from calculation of the coefficients of turbulent momentum and heat transfer invoking model
of �gradient diffusion�. It is show is that the countergradient heat transfer in the local regions is caused by tur-
bulence diffusion or by the term of divergence of triple correlation in the balance equation for the tempera-
ture variance.
The study was carried out under state contract with ICMMG SB RAS (0251-2021-0003), and partially was financially
supported by RFBR / Russian Science Foundation (Project No. 20-01-00560 A).
References
1. Kurbatskii A.F. Modeling Nonlocal Turbulent Momentum and Heat Transfer. Nauka. Novosibirsk. 1988 (in Russian).
2. Kurbatskii A. F. and Kurbatskaya L. I. Investigation of a stable boundary layer using an explicit algebraic model of
turbulence. // Thermophysics and Aeromechanics. 2019. 26, N. 3. P. 335�350.
The use of satellite data in the problems of evaluating the characteristics of emission sources
and atmospheric parameters
A. A. Lezhenin, V. F. Raputa
Institute of Computational Mathematics and Mathematical Geophysics SB RAS
Email: lezhenin@ommfao.sscc.ru
DOI 10.24412/cl-35065-2021-1-01-35
Images from space make it possible to visually record the trajectories of smoke plumes from industrial
chimneys. Based on this information, it is possible to quickly track the spread of impurities in the atmosphere.
The paper analyzes the trajectories of smoke plumes using the equations of hydrothermodynamics and
transport of impurities in the lower atmosphere. With the use of satellite observation data in relation to emis-
sions from high-rise chimneys of thermal power plants located in the Baikal natural territory, estimates of ad-
ditional heights of rise of smoke plumes were obtained according to the previously developed method [1].
Comparison with calculations based on known relationships, including the heat flux of buoyancy and the dy-
namic impulse of emissions, is carried out [2].
The work was supported by the Ministry of Science and Higher Education of the Russian Federation, the grant
No. 075-15-2020-787 for implementation of Major scientific project "Fundamentals, methods and technologies for digital
monitoring and forecasting of the environmental situation on the Baikal natural territory".
References
1. Raputa V.F., Lezhenin A.A. Estimation of the Altitude of Smoke Plumes from Satellite Images // Atmospheric and
Oceanic Optics. 2020. V. 33, N. 5. P. 539�544.
2. Berlyand M. E. Prediction and regulation of air pollution. Leningrad: Gidrometeoizdat 1985.
One approach to restoring conditions on inner boundary in the hydrothermodynamics problem
N. R. Lezina, V. I. Agoshkov
Marchuk Institute of Numerical Mathematics of the RAS
Email: lezina@phystech.edu
DOI 10.24412/cl-35065-2021-1-01-36
In this work a domain decomposition method based on theory of optimal control and adjoint equations
[1] is considered. The initial domain is divided into subdomains by introducing an inner boundary. To formulate
interface conditions �additional unknows� are introduced to the system. The domain decomposition approach
allows to solve the problem, restore boundary conditions on inner boundary and obtain solution of the sys-
tem. The domain decomposition method is numerically studied in the problem of hydrothermodynamics mod-
elling. The numerical experiments with using the domain decomposition algorithm are presented and dis-
cussed. It could be noted that described domain decomposition method could be applied together with varia-
tional data assimilation [2].
This work was supported by Russian Foundation for Basic Research (project No 19-01-00595).