Open access peer-reviewed chapter

The Global Numerical Model of the Earth’s Upper Atmosphere

Written By

Namgaladze Aleksandr, Knyazeva Maria, Karpov Mikhail, Zolotov Oleg, Martynenko Oleg, Yurik Roman, Foerster Matthias and Prokhorov Boris

Submitted: 23 March 2017 Reviewed: 20 September 2017 Published: 20 December 2017

DOI: 10.5772/intechopen.71139

From the Edited Volume

Numerical Simulations in Engineering and Science

Edited by Srinivas P. Rao

Chapter metrics overview

1,012 Chapter Downloads

View Full Metrics


The global numerical first principle 3D model of the upper atmosphere (UAM) for the heights 60–100,000 km is presented. The physical continuity, motion, heat balance and electric potential equations for the neutral, ion and electron gases and their numerical solution method are described. The numerical grids, spatial and time integration steps are given together with the boundary and initial conditions and inputs. Testing and obtained geophysical results are given for many observed situations at various levels of solar, geomagnetic and seismic activity.


  • global numerical model
  • UAM
  • equations
  • solution method
  • grid
  • integration steps
  • input
  • output
  • upper atmosphere
  • thermosphere
  • ionosphere
  • plasmasphere
  • magnetosphere
  • space weather
  • solar and geomagnetic activity

1. Introduction

The upper atmosphere is a part of the gas envelope of the Earth. It is located from the height h ≈ 60 km to several RE of the geocentric distance, where RE = 6371 km is the Earth’s radius. It is characterized by a sharp transition from the predominance of the neutral particles to the charged particles.

The upper atmosphere is divided into several height regions depending on its gas composition and dominating physical process: the thermosphere (from ~80–90 km to ~400–800 km) and the exosphere (above ~400–800 km) in relation to the neutral particles; the ionosphere, the plasmasphere and the magnetosphere in relation to the charged particles.

The ionosphere is located at the height range from ~50 km (the upper part of the middle atmosphere—mesosphere) to ~1000 km. The plasmasphere is located above the ionosphere up to the plasmapause—the geomagnetic force line with L~(2.5 – 7), where L is the geocentric distance to the geomagnetic line top expressed in units of RE. The magnetosphere is the region below the magnetopause, where the solar wind pressure balances the geomagnetic field one. The magnetopause height is ~10RE geocentric distance on the day side and ~100RE on the night one.

The upper atmosphere state is part of Space Weather. It experiences regular annual, seasonal and diurnal variations as well as disturbances caused by the solar, geomagnetic and lithosphere activities, both globally and locally. Along with the meteorological weather, the Space Weather greatly affects human activity. Such disturbances as geomagnetic storms and substorms, auroras lead to disruptions of the radio communication in the HF range up to its blackout, faulty operation of the navigation satellite systems and electronic onboard equipment of the aircrafts. They generate intense geomagnetic-induced currents in the long conducting lines (transmission facilities, telecommunication, cables, railways and oil and gas pipelines), which create failures of the automatic and relay protection systems and, as a consequence, they cause emergency shutdowns of power supply systems, etc. Therefore, monitoring and forecasting of the upper atmosphere state is an extremely important task.

Monitoring of the near-Earth space is conducted by measurements of different physical plasma parameters (temperatures, concentrations and velocities of neutral and charged components, electric and magnetic fields, etc.) at different heights and areas of the globe, both ground-based and on the board of airplanes, rockets and satellites. Despite the growing level of the technical perfectness, it is impossible to provide stable and global monitoring, and “white areas” still remain over the Earth. These areas are those where measurements cannot be conducted or experience various difficulties (especially over the oceans and near the poles). In such cases, the method of the numerical simulation becomes valuable, and the calculations of the desired parameters fill the “white gaps”. The numerical models use the basic fundamental laws of nature to describe quantitatively the near-Earth environment and/or interpret the measurements.

The Upper Atmosphere Model (UAM), described in this chapter, is a global, 3D, self-consistent, numerical model covering h = 60 − 100,000 km. It solves the system of equations describing the basic laws of mass, energy, momentum and electric charge conservation and calculates variations of the neutral (O2, N2, O, H) and charged (electrons and ions O2+, NO+, O+, H+) particles’ concentrations; temperatures of neutral gas, ions and electrons; the velocities of the particles as well as the electric field potential for any geo- and heliophysical conditions. The UAM takes into account the non-coincidence of the Earth’s geographic and geomagnetic axes.

The UAM was developed previously in Kaliningrad under the supervision of Prof. A.A. Namgaladze as the Global Self-consistent Model of the Thermosphere, Ionosphere and Protonosphere (GSM TIP) [1, 2] and further improved in Murmansk. The modern version called as UAM [3, 4] differs from the GSM TIP by implementation of algorithms for the integration with variable latitude steps and by including empirical models of the ionosphere and thermosphere to use their data for initial and boundary conditions and for the UAM testing.

Most of the modern numerical models (NCAR TIE-GCM [5], CTIM [6], CTIPe [7], GITM [8], SWMF [9, 10]) either cover the near-Earth space in very limited ranges in heights and latitudes, or perform a combination of several models without physical coupling between the parts. The UAM covers the near-Earth space as a coupled system and is still superior to all existing models by spatial coverage and resolution. This makes the UAM suitable to investigate a variety of physical processes, both globally and locally.


2. Basic equations

The fundamental conservation laws are used in the UAM: the quasi hydrodynamic equations of the continuity (1), momentum (2), heat balance (3) and electric current continuity (4).


where nα is the concentration of the α-th gas, vα is its velocity vector, Qα and Lα are its production and loss rates, respectively, ρα is the mass density of the α-th gas, Ω is the angular velocity of the Earth’s rotation, r is the radius-vector directed from the center of the Earth, Fα is the force vector acting on the unit gas volume, pα = nαkTα is the partial pressure of the α-th gas, Tα is its temperature, k is Boltzmann’s constant, g is the vector of the gravity acceleration, mα is the mass of the α-th gas particle, μαi is the reduced mass of colliding particles, vαi is the collision frequency, vαvi is the relative velocity of the α-th and i-th gases, e is the elementary charge, E and B are the electric and magnetic fields, respectively, c is the specific heat per unit volume of the α-th gas, λα is the thermal conductivity, P is the rate of heat of the α-th gas, P is its rate of cooling, P is the rate of heat exchange between the α-th gas and other particles and j is the electric current density.

Due to the changes in gas composition and the predominance of different physical processes at different heights, the UAM is divided into several computational blocks: the neutral atmosphere and lower ionosphere block; the F2-layer and plasmasphere block and the magnetosphere block. A special block calculates the electric field potential by solving the continuity equation for j. Each block covers a particular height range, uses its own particular equations, coordinate system and calculates its own set of parameters and exchanges with other blocks at each time step of numerical integration.

The block of the neutral atmosphere and lower ionosphere uses the spherical geomagnetic coordinate system and calculates the 3D variations of nn for O, O2 and N2 gases, ne, Tn and v, where ​​v → ​​ is the neutral wind velocity vector, within the height range from the model’s lower boundary to 520 km, as well as Ti and Te up to 175 km.

Concentrations n(N2) are obtained from the barometric law, n(O) and n(O2) from the continuity Eq. (1), which takes view as:


where v is the neutral wind velocity, vdn is the diffusion velocity vector, which has only a vertical component equal to the sum of molecular and turbulent diffusion velocities, Qn and Ln are the production and loss rates, respectively, taking into account the photodissociation of the molecular oxygen by the solar radiation and recombination of O and O2 in the triple collisions and radiative recombination.

The meridional and zonal components of v are obtained from the projection of Eq. (2) on the horizontal axes of the geographical coordinate system:


where ρ is the mean neutral density, p=nnnkTn is the mean partial pressure, Rnihor=niµniνninnvvihor is the horizontal projection of the neutral-ion friction force, η is the coefficient of viscosity, n and i lowercases stand for the neutral particles and ions, respectively, νni is the frequency of collisions between the neutrals and ions, vvi is the relative velocity of the neutrals and ions.

Neutral density ρ is calculated from the projection of Eq. (2) on the radii vector r directed from the center of the Earth to the particular point in space, taking into account only the gravity force and vertical component of ∇p, that is, fulfilling the hydrostatic equilibrium:


where g is the vertical projection of the sum of the gravitational and centrifugal accelerations.

By summing up for n the continuity Eq. (5) and assuming that the total number of particles remains constant, that is, nQn=nLn,we obtain:


With the calculated ρ and horizontal components of v, we calculate the vertical wind velocity from Eq. (8).

Finally, Tn is calculated from the heat balance Eq. (3) written as:


where cV and λn are the specific heat per unit volume and thermal conductivity, respectively, PnQUV, PnQJ and PnQC are the rates of the heating by the ultraviolet (UV) and extra-ultraviolet (EUV) solar radiation, Joule heating due to collisions with ions and heating by particles precipitating from the magnetosphere, respectively; PnL is the rate of the cooling due to radiation.

For concentrations of O2+, N2+ and NO+ ions in the D, E and F1 layers, the photochemical, local heating and heat exchange processes dominate over the transport processes, since the lifetime of the molecular ions is many times smaller than the transport characteristic time due to high collision frequencies of the charged particles with neutrals and each other. Thus, Eq. (1) for the molecular ions in the lower ionosphere block is written as:


where n(XY+) = ne is the total concentration of the molecular ions, Q(XY+) and L(XY+) are the production and loss rates due to the ionization by direct UV solar radiation, scatter radiation and ionization by the precipitating electrons, as well as the result of the ion-molecular reactions and dissociative recombination.

The momentum equation for the molecular ions is written as:


In the momentum equation for electrons, we neglect all terms containing me due to its smallness:


This gives the electric field of the ambipolar diffusion:


The temperatures Ti and Te are calculated in this block taking into account that up to h = 175 km the local heating dominates, and the processes of the heating transport are neglected in the heat balance Eq. (3):


where PiQJ is the rate of Joule’s heating, PiTe=PiTn are the rates of the ions’ heat exchange with electrons and neutrals, respectively, PeTi and PeTn are the rates of the electrons’ heat exchange with ions and neutrals, respectively, PeQP and PeQC are the rates of the electrons’ heating by the photoelectrons and electrons precipitating from the magnetosphere, respectively.

In the block of the F2-layer and plasmasphere n(O+), n(H+), vi, Ti and Te are calculated for h = 175 − 100,000 km. At these heights collision frequencies νin of charged particles with neutrals are much smaller than ions’ gyrofrequencies Ωi, that is, collisions do not interfere the cyclotron rotation and drift of charged particles. Ions and electrons are fully magnetized, that is, they are tied to the geomagnetic field lines and can move across the lines only under the action of an extraneous force. Because the geomagnetic field has such great influence on the behavior of ions and electrons, we use the magnetic dipole coordinate system in this block. The system of modeling Eqs. (1)(3), thus, is integrated along the field lines taking into account the electromagnetic plasma drift perpendicular to the magnetic field.

In the continuity equation, along with the production rates Qi and loss rates Li of ions O+ and H+ by photoionization, corpuscular ionization, chemical reactions between O+, O2 and N2, as well as the charge exchange between O+ and H and between O and H+, the transport of ionized components (the divergence of charged particles’ flow) is taken into account:


Here, superscripts par and per refer to directions parallel and perpendicular to the geomagnetic field lines, respectively. The operator


gives the Lagrangian time derivatives along the trajectory of the charged particle’s electromagnetic drift perpendicular to the geomagnetic field line with the velocity:


The electrons’ velocity along the field lines is calculated as:


The atomic ions’ motion equation for the i-th atomic ion is given as:


where subscripts i, j, e and n refer to O+, H+, electrons and neutrals, respectively and gpar is the projection of the sum of the gravity and centrifugal accelerations along the geomagnetic field line.

The temperatures Ti and Te are calculated from the heat balance equations written as:


where PiQJ is the rate of Joule’s heating for the ion gas; PiTj, PiTe and PiTn are the rates of the heat exchange between ions, between ions and electrons and between ions and neutrals, respectively; PeTi, PeTj and PeTn are the rates of the electrons’ heat exchange with the i-th and j-th ions as well as the neutral gas, respectively; PeQP and PeQC are the rates of the electrons’ heating by the photoelectrons and electrons precipitating from the magnetosphere, respectively.

The magnetospheric block of the UAM calculates the plasma layer ni from the continuity Eq. (21), vi from the motion Eq. (22), ion pressure pi from Eq. (23) and the field-aligned current density j (FACs) from Eq. (24) above 175 km:


where V is a half-volume of the geomagnetic field tube, γ = 5/3 is the adiabatic coefficient for the plasma layer ions and ez is the unit vector along B. V is calculated as:


where z is a distance along B, zmax is the distance to the top of the geomagnetic field line.

We assume that the pressure of the plasma layer is isotropic and constant along B. The pressure of the plasma layer electrons is neglible in comparison with pi, that is, the electrons are considered to be cold.

In the block of the electric field potential ϕ, it is calculated taking into account electric fields of magnetospheric, thermospheric (dynamo), and lower atmospheric (due to the charges transfer into the ionosphere from bellow) origins. The equation of the electric current continuity is written as:


where jm and js are the densities of the magnetosphere and lower ionosphere electric currents, and ji is the ionosphere electric current density given by Ohm’s law for the plasma:


Here, σP and σH are the Pedersen and Hall conductivities, respectively, and σpar is the conductivity along B.


where Ωi and Ωe are the ion and electron cyclotron frequencies, respectively.


where E is the electric field vector defined as the sum of the electrostatic field with the potential ϕ and the induced dynamo field:


Thus, the equation for the electric current continuity is given as:


where σ̂ is the tensor of the ionospheric electric conductivity for the coordinate system with axes along B, E and E×B:


The integration of Eq. (27a) is conducted along the electric conducting layer, from the lower boundary of the UAM up to 175 km, where we neglect the dependency of the electric field intensity on h. Above 175 km, we assume that the plasma is magnetized and the geomagnetic field lines are electrically equipotential.


3. Equations solutions: initial and boundary conditions

The integration of the model equations is carried out using a finite-difference numerical method. The near-Earth environment is considered as a discrete three-dimensional grid, and each computational block uses its own coordinate system. The derivative operators are replaced by the difference ratios, and the solution is obtained in the nodes of the numerical grid. Steps in the numerical grid are chosen depending on the task and the characteristic scale of the investigated process. Usually, the integration step in longitude is chosen as a constant value between 5° and 15°. The integration steps in latitude are variable and depend on the latitude. Near the geomagnetic equator the step is chosen in the range of 2.5–5.0° and 1° near the poles, because a tighter numerical grid is required for a precise calculation in the aurora zones and across the polar caps. Steps of the numerical integration along the height are also variable. In the blocks where a spherical coordinate system is used, the step is 3 km at the lower boundary of the UAM and increases exponentially with altitude.

The solutions of Eqs. (1)(4), which are containing derivations of coordinates, require boundary conditions: the distribution of the desired parameters at the boundaries of the numerical blocks. At the lower boundary Tn and nn are obtained from the empirical model NRLMSISE-00 (further simply MSIS [11]), v is calculated from Eqs. (6) and (8) neglecting the viscosity and ion-neutral friction, in the geostrophic approximation, when Eq. (6) contains only the Coriolis acceleration and the pressure gradient of the neutral gas.

For the upper boundary conditions in the block of the neutral atmosphere (at 520 km) the diffusion equilibrium for nn is used, v and Tn are considered to be independent of the height:


In the F2-layer and plasmasphere block the boundary conditions are defined as follows. At the altitude of 520 km, the concentration of the neutral hydrogen is set according to the empirical model of the neutral atmosphere [12]. For atomic ions ni at the bases of the geomagnetic field lines at 175 km, the production rates are equal to the loss rates:


Eqs. (14) and (15) are used for the Ti and Te at 175 km.

The model equations are integrated along the geomagnetic field lines in the areas with closed lines of force (±75° in geomagnetic latitude). In the regions with open geomagnetic field lines (poleward of 75° geomagnetic latitude) ni, as well as the heat fluxes of the ions and electrons are set equal to zero at the upper boundary of the UAM. Thus, the model reproduces the condition of the polar wind: a supersonic outflow of plasma from the F2-layer and upper ionosphere of the polar caps along the open lines of the geomagnetic field.

Solutions of those equations, which are containing time derivatives, require initial conditions: parameter distributions at the start of the numerical calculation. For quiet geophysical conditions, stationary solutions are used as initial conditions after the multiple runs of the model. The initial conditions for perturbed periods are the solutions obtained for the previous quiet days. It is also possible to use data from empirical models (NRLMSISE-00 [11], HWM-93 [13], IRI-2001 [14] or IRI-2007 [15]) as initial conditions.


4. Model inputs

The inner state of the simulated space and external forcings acting on it are characterized by the model inputs set up by the user: (1) date and time to set an initial placement of the numerical grid nodes relative to the Sun; (2) spectra of the solar UV and EUV radiation; (3) fluxes of the high energetic electrons precipitating from the magnetosphere; (4) the FACs connecting the ionosphere with the magnetosphere and/or (5) the distribution of ϕ at the boundaries of the polar cap; (6) the local js flowing through the lower boundary from below; (7) indices of the geomagnetic activity and (8) components of the interplanetary magnetic field (IMF) and solar wind.

The solar UV and EUV spectra define the coefficients of O2 dissociation and O2+, N2+, NO+ and O+ production rates due to the photoionization of the corresponding neutral components. The UV and EUV spectra dependence on the solar activity is set up according to Ref. [16]. The intensity of the night scatter radiation intensity is 5 kR for the emission with wave length λ = 121.6 nm and 5 R for the rest emission lines (102.6, 58.4 and 30.4 nm).

The precipitating electrons’ fluxes are set up at the upper boundary of the thermosphere, at 520 km, and their intensity I is written as:


where Φ and Λ are the geomagnetic latitude and longitude, respectively (Λ = 0 corresponds to the midday magnetic meridian); E is the energy of the precipitating electrons; Im(E) is the maximum intensity of the precipitating electron flux; ∆Φ and ∆Λ are the half-widths of the maximum precipitations in latitude and longitude; Φmd and Φmn are the magnetic latitudes of the maximum precipitations at the midday and midnight meridians, respectively. Specific values for the precipitating parameters in Eqs. (32) and (33) are taken from the empirical models [17, 18].

The magnetospheric sources jm in Eq. (27a) specify the distribution of FACs. The FACs of the Region 1 (R1) flow from the magnetosphere into the ionosphere on the dawn side and out of the ionosphere on the dusk side, at latitudes higher ±75°. The FACs of the Region 2 (R2) flow in the opposite direction, at areas equatorward of the R1 currents. The distribution of current densities depends on the geophysical conditions and is setup in the UAM in several different ways depending on the task, either as distribution of the FACs in the R1, 2 and the cusp region according to the model [19], or as the distribution of electric potential at the boundary of the polar cap [20] with the FACs in the cusp region and the R2.

The so-called seismogenic electric currents are the vertical electric currents switched on to simulate the ionosphere effects of various mesoscale phenomena in the lower ionosphere, such as earthquakes, thunderstorms, etc. Used as a model input, the vertical js are added locally to Eq. (27a) at the lower boundary of the UAM.

The geomagnetic activity is used in the UAM by setting up the planetary geomagnetic indices Kp and Ap, indicating global geomagnetic field disturbance, and aurora indices AE, AL and AU, where AU and AL indicate, respectively, the largest increase and lowest decrease of the Northern component of the geomagnetic field in comparison to the background (quiet) value, AE is the sum of AL and AU and characterizes the largest scale of the magnetic field during the substorm in the high-latitude regions.


5. The UAM versions

The UAM provides the possibility of integrating various empirical models and data of the upper atmosphere. The comparison of the self-consistent UAM version and the UAM versions with different combinations of the empirical models allows testing both the UAM and the empirical models.

In the UAM-MSIS version, nn and Tn are calculated directly from the MSIS [11]. The thermospheric circulation is calculated by the numerical solution of Eqs. (6) and (8) where ∇p from the MSIS is used. Finally, the MSIS is used to set up Tn and nn at the lower boundary and initial conditions as well.

In the UAM-HWM version the distribution of the horizontal thermospheric wind is calculated using the empirical model HWM-93 [13]. The vertical component of the wind velocity is calculated by the numerical solution of the continuity equation for ρ (Eq. (8)). The HWM-93 is used for the set of the initial conditions and for comparison of the theoretical model winds with observations.

The magnetospheric block of the UAM simulates the transport processes in the plasma sheet by solving the system of the equations for the plasma sheet ions (see item 2). In Ref. [21], the initial values are taken as pi = 0, 4 nPa and ni = 0, 4 cm−3, correspondingly. The program produces more or less realistic pi distribution and R2 FACs. The problem is that the obtained solution is not stable, and it falls apart after approximately 1 h.

There are several ways to set up FACs spatial-temporal distributions in the UAM, such as empirical data from the magnetic field measurements from the Dynamics Explorer 2 [22] and the Magsat satellites [23], the FACs empirical models by Papitashvili [24] in [25, 26], by Lukianova [27] in Ref. [28] and MFACE [29] in Ref. [30]. All these versions with various FACs take into account the dependence of FACs on the interplanetary magnetic field (IMF). Such methods of setting the FACs distribution allow using any other empirical data of FACs.

In the UAM version [31], the positions of the auroral oval boundaries, the values of electron flux intensities and the latitudes and longitudes of the intensity maxima were set from precipitation patterns observed by DMSP. The spectra of the precipitating electrons are assumed to be Maxwellian in this case.

The UAM-P version [32, 33], created in Potsdam, differs from the UAM by the electric field block simulation. This block uses magnetic dipole coordinates instead of spherical geographical ones within h = 80 − 526 km. It is assumed that E does not change along B inside the ionospheric current layer. After the integration, it is also represented as a 2D-distribution. It allows to exclude the conductivity parallel to the magnetic field and to keep the vertical electric field inside the current-carrying layer. As a result, the lower latitudinal and equatorial electric field distributions as well as the current system of these areas look more correctly.

The Canadian Ionosphere and Atmosphere Model (Canadian IAM or C-IAM) is comprised of the extended Canadian Middle Atmosphere (CMAM) and the UAM, currently coupled in a one-way manner [34]. This version was used to investigate the physical mechanisms responsible for forming the four-peak longitudinal structure of the 135.6 nm ionospheric emission observed by the IMAGE satellite over the tropics at 20:00 local time from March 20 to April 20, 2002. The study showed that main mechanism is driven by the diurnal eastward propagating tide with zonal number 3.


6. Results

During its development and improvement, the UAM was used to perform a number of numerical experiments aimed at testing and comparison of calculation results with measurements and other models. The UAM successfully showed its ability to reproduce the general behavior of ionospheric and thermospheric parameters such as at low and high geomagnetic and solar activity conditions. A good agreement of the numerical simulations’ results was achieved in comparison with observations by incoherent scatter radars located at various latitudes and longitudes (ISRs) [31], digital ionosonde CADI in the Voeykovo Main Geophysical Observatory [35], several chains of the ionospheric tomographical receivers [36, 37, 38], satellites, including CHAMP and GPS [3, 4, 26, 30, 39, 40, 41, 42], as well as with empirical models such as various types of IRI, MSIS, HWM [42, 43, 44, 45, 46, 47], etc. and other theoretical models [42].

Numerical modeling of the upper atmosphere behavior during substorms. These relatively short disturbances of the geomagnetic field have duration from about 0.5 to 3 h. They are generated at the magnetopause and in the magnetosphere tail via the magnetic field lines reconnection processes and connected with the polar upper atmosphere in the auroral zones via FACs including the current wedges. The substorm auroral currents are reflected by the auroral magnetic activity indexes AL, AU and AE. The UAM takes into account these indexes. This allows the modeling of the upper atmospheric behavior during substorms via the UAM simulations. The results were presented in [4, 48, 49, 50, 51, 52, 53] including the cusp and auroral oval behavior, energetic magnetospheric electron precipitations, electric fields, current wedge and internal atmospheric gravity waves generation. The main role of the thermospheric heating due to the soft electron precipitation was shown for the thermosphere substorm effects.

Numerical modeling of the upper atmosphere behavior during geomagnetic storms. These geomagnetic disturbances are global and have duration of about several days. Usually they include several substorms on the background of the geomagnetic field depression created by the ring current development. The geomagnetic activity indexes Dst and Kp characterize the geomagnetic storms. These indexes are included in the UAM inputs as well. This allows to model the upper atmospheric behavior during magnetic storms via the UAM simulations. The calculation results were presented in [25, 26, 30, 36, 37, 38, 54, 55, 56, 57, 58, 59] including the main ionospheric trough dynamics due to the ionosphere-magnetosphere convection and non-coincidence of the geographical and geomagnetic axes of the Earth. The physical mechanisms of the negative and positive F2-layer ionospheric storms (electron concentration decreases and increases, correspondingly) formation were described. The main role of the thermospheric composition (atomic to molecular neutral gas concentration ratio) and winds disturbances in the magnetic storm ionospheric effects was demonstrated in these UAM calculations.

In addition to the UAM testing and comparison with the observations for the different levels of solar and geomagnetic activities, the model has been widely used to study the ionosphere response on the local sources in the lower atmosphere, such as disturbances associated with the processes of the earthquakes’ preparation processes [26, 60, 61, 62, 63]. Numerical UAM calculations showed that the electric fields of 5–10 mV/m effects on the F2-layer plasma by the electromagnetic drift in the crossed geomagnetic field and the electric field of the seismic origin. The vertical electric current, flowing through the lower boundary of the ionosphere with the density of ~20 nA/m2, is required for the generation of the seismogenic electric field of ~10 mV/m [60, 61, 62, 63, 64]. The important role of the aerosols over the tectonic faults was underlined in this process due to the very low recombination rate for the charged aerosols.

Thus, the UAM was tested and used in many helio- and geophysical situations. Nevertheless, the amount of the UAM simulations remains to be insufficient despite the IT progress. This is related with the specifics of the geophysics as science at all. The near-Earth space environment varies due to the solar, seismic and human activities. This does not allow performing the repeated fixed experiments as in usual physical laboratories to obtain correctly the standard statistical error estimates. Moreover, the observations themselves are very limited. None of them has 3D spatial and time resolutions satisfying to the requirements of the modern technical means of the Space weather practical usage. This was well known long ago [65] and such models as UAM are aimed solving this important problem.


7. Conclusions

Further development of the UAM means a huge amount of further numerical experiments to its mathematical and physical quality. These experiments have to take into account all modern achievements of the numerical mathematics and computer science. The numerical grids, steps, various iterations, etc. should be tested to find their optimal combinations for the best stability and accuracy. The user’s manuals should be constructed, including the UAM website. The UAM prognostic features have to be improved by modeling many case studies for various helio- and geophysical situations especially for geomagnetic and seismogenic disturbances. Comparisons with ground-based and satellite observations, empirical and other theoretical models have to be made continuously. The frame approach should be widely used by including separate observational, empirical and theoretical blocks into the UAM, such as the real geomagnetic field, polar wind, plasma sheet, electric fields, lower atmosphere, aerosols, tides, etc. An international cooperation is absolutely necessary for these future scientific UAM perspectives.



We thank many people who have worked with the UAM, developed and improved it, such as A. N. Namgaladze, M. A. Volkov, E. N. Doronina, Yu. V. Romanovskaya, I. V. Korableva, Yu. A. Shapovalova, M. G. Botova, M. I. Rybakov, I. V. Artamonov, E. V. Vasilieva, V. A. Medvedeva, V. A. Shlykov, L. A. Chernyuk and many others.


  1. 1. Namgaladze AA, Korenkov YN, Klimenko VV, Karpov IV, Bessarab FS, Surotkin VA, et al. Global model of the thermosphere-ionosphere-protonosphere system. Pure and Applied Geophysics. 1988;127(2/3):219-254. DOI: 10.1007/BF00879812
  2. 2. Namgaladze AA, Korenkov YN, Klimenko VV, Karpov IV, Surotkin VA, Naumova NM.Numerical modeling of the thermosphere-ionosphere-protonosphere system. Journal of Atmospheric and Solar—Terrestrial Physics. 1991;53(11/12):1113-1124. DOI: 10.1016/0021-9169(91)90060-K
  3. 3. Namgaladze AA, Martynenko OV, Namgaladze AN. Global model of the upper atmosphere with variable latitudinal integration step. International Journal of Geomagnetism and Aeronomy. 1998;1(1):53-58
  4. 4. Namgaladze AA, Martynenko OV, Volkov MA, Namgaladze AN, Yurik RY. High-latitude version of the global numerical model of the Earth’s upper atmosphere. Proceedings of the MSTU. 1998;1(2):23-84
  5. 5. Qian L, Burns AG, Emery BA, Foster B, Lu G, Maute A, et al. The NCAR TIE-GCM. In: Huba J, Schunk R, Khazanov G, editors. Modeling the Ionosphere-Thermosphere System. 1st ed. Chichester, UK: John Wiley & Sons, Ltd; 2014. DOI: 10.1002/9781118704417.ch7
  6. 6. Fuller-Rowell TJ, Rees D, Quegan S, Moffett RJ, Codrescu MV, Millward GH. A Coupled thermosphere-ionosphere model (CTIM). In: Schunk RW, editor. STEP Handbook of Ionospheric Models. 1st ed. Utah, USA: Utah State University and SCOSTEP (Scientific Committee on Solar Terrestrial Physics); 1996. p. 217-238
  7. 7. Codrescu MV, Negrea C, Fedrizzi M, Fuller-Rowell TJ, Dobin A, Jakowsky N, et al. A real-time run of the Coupled Thermosphere Ionosphere Plasmasphere Electrodynamics (CTIPe) model. Space Weather. 2012;10(2). DOI: 10.1029/2011SW000736
  8. 8. Ridley AJ, Deng Y, Toth G. The global ionosphere–thermosphere model. Journal of Atmospheric and Solar—Terrestrial Physics. 2006;68(8):839-864. DOI: 10.1016/j.jastp.2006.01.008
  9. 9. Tóth G, Holst B, Sokolov IV, De Zeeuw DL, Gombosi TI, Fang F, et al. Adaptive numerical algorithms in space weather modeling. Journal of Computational Physics. 2012;231(3):870-903. DOI: 10.1016/
  10. 10. Tóth G, Sokolov IV, Gombosi TI, Chesney DR, Clauer CR, De Zeeuw DL, et al. Space Weather Modeling Framework: A new tool for the space science community. Journal of Geophysical Research. 2005;110(A12). DOI: 10.102 9/2005JA011126
  11. 11. Picone JM, Hedin AE, Drob DP, Aikin AC. NRLMSISE-00 empirical model of the atmosphere: Statistical comparisons and scientific issues. Journal of Geophysical Research. 2002;107(A12). DOI: 10.1029/2002JA009430
  12. 12. Jacchia LG. Thermospheric temperature, density and composition: New models. Journal of Geophysical Research. 1977;375:1-106
  13. 13. Hedin AE, Fleming EL, Manson AH, Schmidlin FJ, Avery SK, Clark RR, et al. Empirical wind model for the upper, middle and lower atmosphere. Journal of Atmospheric and Terrestrial Physics. 1996;58:1421-1447. DOI: 10.1016/0021-9169(95)00122-0
  14. 14. Bilitza D. International reference ionosphere 2000. Radio Science. 2001;36:261-275
  15. 15. Bilitza D, Reinisch BW. International reference ionosphere 2007: Improvements and new parameters. Advances in Space Research. 2008;42:599-609
  16. 16. Ivanov-Kholodny GS, Nusinov AA. Shortwave solar radiation and its impact on the atmosphere and ionosphere (in Russian). Investigations of Outer Space. 1987;26:80-154
  17. 17. Fuller-Rowell TJ, Evans DS. Height-integrated Pedersen and hall conductivity patterns inferred from the Tiros-NOAA satellite data. Journal of Geophysical Research. 1987;92(7):7606-7618
  18. 18. Hardy DA, Gussenhoven MS, Holeman E. A statistical model of auroral electron precipitation. Journal of Geophysical Research. 1985;90:4229-4248
  19. 19. Iijima T, Potemra TA. The amplitude distribution of field-aligned currents at northern high latitudes observed by Triad. Journal of Geophysical Research. 1976;81:2165-2174
  20. 20. Weimer DR, Maynard NC, Burke WJ, Liebrecht C. Polar cap potentials and the auroral electrojet indicies. Planetary and Space Science. 1990;38(9):1207-1222. DOI: 10.1016/0032-0633(90)90028-O
  21. 21. Artamonov IV, Vasilyeva EV, Medvedeva VA, Namgaladze AA. A new version of the magnetospheric block for the global numerical upper atmosphere model. In: Physics of Auroral Phenomena, XXVI Annual Apatity Seminar, Abstracts. Apatity, Russia: PGI of RAS; 2003. p. 21
  22. 22. Artamonov IV, Vasilyeva EV, Namgaladze AA, Martynenko OV. On the electric potential pattern corresponding region 2 field-aligned currents derived from the DE2 satellite measurements. In: Physics of Auroral Phenomena. Procceedings of the XXIX Annual Seminar; Apatity, Russia. 2006. p. 71-74
  23. 23. Artamonov IV, Vasilyeva EV, Martynenko OV, Namgaladze AA. Mathematical simulation of the electric potential distributions using region 2 field-aligned currents obtained from different satellite data. In: Procceedings of the 6th International Conference “Problems of Geocosmos”, Saint-Petersburg, Russia. 2006. p. 20-23
  24. 24. Papitashvili VO, Christiansen F, Neubert T. A new model of field-aligned currents derived from high-precision satellite magnetic field data. Geophysical Research Letters. 2002;29(14). DOI: 10.1029/2001GL014207
  25. 25. Förster M, Prokhorov BE, Namgaladze AA, Holschneider M. Numerical modeling of solar wind influences on the dynamics of the high-latitude upper atmosphere. Advances in Radio Science. 2012;10:299-312. DOI: 10.5194/ars-10-299-2012
  26. 26. Namgaladze AA, Förster M, Prokhorov BE, Zolotov OV. Electromagnetic drivers in the upper atmosphere: Observations and modeling. In: Bychkov V, Golubkov G, Nikitin A, editors. The Atmosphere and Ionosphere, Physics of Earth and Space Environments. 1st ed. Houten, Netherlands: Springer; 2013. p. 165-219. DOI: 10.1007/978-94-007-2914-8_4
  27. 27. Lukianova R, Christiansen F. Modeling of the global distribution of ionospheric electric fields based on realistic maps of field-aligned currents. Journal of Geophysical Research. 2006;111:A03213. DOI: 10.1029/2005JA011465
  28. 28. Knyazeva MA, Namgaladze AA, Beloushko KE. Field-aligned currents influence on the ionospheric electric fields: Modification of the upper atmosphere model. Russian Journal of Physical Chemistry B. 2015;9(5):758-763. DOI: 10.1134/S1990793115050206
  29. 29. He M, Vogt J, Lühr H, Sorbalo E, Blagau A, Le G, et al. A high-resolution model of field-aligned currents through empirical orthogonal functions analysis (MFACE). Geophysical Research Letters. 2012;39:L18105. DOI: 10.1029/2012GL053168
  30. 30. Prokhorov BE, Forster M, Namgaladze AA, Holschneider M. Using MFACE as input in the UAM to specify the MIT dynamics. Journal of Geophysical Research, Space Physics. 2014;119(8):6704-6714. DOI: 10.1002/2014JA019981
  31. 31. Namgaladze AA, Zubova YV, Namgaladze AN, Martynenko OV, Doronina EN, Goncharenko LP, et al. Modelling of the ionosphere/thermosphere behaviour during the April 2002 magnetic storms: A comparison of the UAM results with the ISR and NRLMSISE-00 data. Advances in Space Research. 2006;37(2):380-391. DOI: 10.1016/j.asr.2005.04.013
  32. 32. Prokhorov BE, Förster M, Stolle C, Lesur V, Namgalagze AA, Holschneider M. The ionospheric current system and its contribution to the Earth’s magnetic field. Geophysical Research Abstracts. 2016;18:ST4.1, EGU2016-6902
  33. 33. Prokhorov BE, Förster M, Lesur V, Lesur V, Namgalagze AA, Holschneider M, Stolle C. Using the UAM-P model to specify the additional magnetic field generated by the system of the ionospheric currents. Geophysical Research Abstracts. 2017;19:EMRP2.1/ESSI1.15/GI1.12, EGU2017-10822
  34. 34. Martynenko OV, Fomichev VI, Semeniuk K, Beagley SR, Ward WE, McConnell JC, et al. Physical mechanisms responsible for forming the 4-peak longitudinal structure of the 135.6 nm ionospheric emission: First results from the Canadian IAM. Journal of Atmospheric and Solar - Terrestrial Physics. 2014;120:51-61. DOI: 10.1016/j.jastp.2014.08.014
  35. 35. Rybakov MV, Namgaladze AA, Karpov MI. Comparison of ionospheric parameters calculated with UAM and measured at Voeykovo observatory. Geomagnetism and Aeronomy. 2016;56(5):604-609. DOI: 10.1134/S0016793216040186
  36. 36. Namgaladze AN, Evstafiev OV, Khudukon BZ, Namgaladze AA. Model interpretation of the ionospheric F-region electron density structures observed by ground-based satellite tomography at sub-auroral and auroral latitudes in Russia in January-May 1999. Annales Geophysicae. 2003;21(4):1005-1016
  37. 37. Shapovalova YA, Namgaladze AA, Namgaladze AN, Khudukon BZ. Stratification of the main ionospheric trough as an effect of the noncoincidence of the Earth’s geomagnetic and geographic axes. In: Proceedings of the XXVI Annual Apatity Seminar “Physics of Auroral Phenomena”. PGI; 2003;26:87-90
  38. 38. Korableva IV, Namgaladze AA, Namgaladze AN. High-latitude ionosphere during magnetic storms of October 26, 2003–November 1, 2003: Tomographic reconstructions and numerical modeling. Geomagnetism and Aeronomy. 2008;48(5):642-651
  39. 39. Namgaladze AA, Doronina EN, Förster M. The role of electric fields and magnetospheric electron precipitations for the formation of the equtorial total mass density minimum. In: “Physics of Auroral Phenomena”, Proc. 29 Annual Seminar. 2006. pp. 238-240
  40. 40. Doronina EN, Namgaladze AA, Förster M.A model interpretation of the CHAMP neutral mass density measurements. In: Proceedings of the 6th International Conference “Problems of Geocosmos” St. Petersburg, Russia. 2006. p. 58-61
  41. 41. Förster M, Namgaladze AA, Doronina EN, Prokhorov BE. High-latitude thermospheric winds: Satellite data and model calculations. Russian Journal of Physical Chemistry B. 2011;5(3):439-446. DOI: 10.1134/S1990793111030043
  42. 42. Shim JS, Rastaetter L, Kuznetsova M, Bilitza D, Codrescu M, Coster AJ, et al. CEDAR-GEM Challenge for Systematic Assessment of Ionosphere/Thermosphere Models in Predicting TEC during the 2006 December Storm Event. Space Weather. 05 October 2017;15:19. DOI: 10.1002/2017SW001649
  43. 43. Martynenko OV, Namgaladze AA, Namgaladze AN, Shlykov VA. Numerical modeling of the longitudinal variations in the near-earth plasma, In: Physics of Auroral Phenomena, Proceedings of the XXIII Annual Apatity Seminar, Apatity PGI-00-01-108. 2001. pp. 49-52
  44. 44. Knyazeva MA, Namgaladze AA. A model study of the seasonal and solar activity variations of the enhanced electron density regions in the night-time ionospheric F2-layer and plasmasphere. In: Physics of Auroral Phenomena, Proc. XXIX Annual Seminar, Apatity, Kola Science Centre, Russian Academy of Science. 2006. pp. 225-228
  45. 45. Knyazeva MA, Namgaladze AA. An influence of the thermospheric wind variations on the enhanced electron density regions in the night-time ionospheric F2-layer and in the plasmasphere. In: Proc. 6th International Conference “Problems of Geocosmos”, Saint-Petersburg State University. 2006. pp. 91-94
  46. 46. Knyazeva MA, Namgaladze AA, Zubova Yu V. Numerical modeling of the Weddel Sea anomaly. In: Proceedings of the 8th International Conference Problem of the Geocosmos, St. Petersburg Petrodvorets. 2010. pp. 121-125
  47. 47. Botova MG, Romanovskaya YV, Namgaladze AA. Latitudinal variations and altitude profiles of Ionospheric parameters: Comparison of theoretical and empirical model results. Russian Journal of Physical Chemistry B. 2015;9(5):764-769. DOI: 10.1134/S1990793115050164
  48. 48. Namgaladze AA, Martynenko OV, Namgaladze AN, Volkov MA, Korenkov YN, Klimenko VV, et al. Numerical simulation of an ionospheric disturbance over EISCAT using a global ionospheric model. Journal of Atmospheric and Terrestrial Physics. 1996;58(1):297-306. DOI: 10.1016/0021-9169(95)00037-2
  49. 49. Volkov MA, Namgaladze AA. Models of field-aligned currents needful to simulate the substorm variations of the electric field and other parameters observed by EISCAT. Annales Geophysicae. 1996;14(12):1356-1361. DOI: 10.1007/s00585-996-1356-0
  50. 50. Leontyev SV, Namgaladze AA, Namgaladze AN, Bogdanov NN. Thermospheric meridional winds in the vicinity of the auroral zone: Observations and modeling. Journal of Atmospheric and Solar - Terrestrial Physics. 1998;60(2):215-226. DOI: 10.1016/S1364-6826(97)00055-2
  51. 51. Namgaladze AA, Namgaladze AN, Volkov MA. Numerical modelling of the thermospheric and ionospheric effects of magnetospheric processes in the cusp region. Annales Geophysicae. 1996;14(12):1343-1355. DOI: 10.1007/s00585-996-1343-5
  52. 52. Namgaladze AA, Namgaladze AN. Numerical modelling of the thermospheric and ionospheric effects of the soft electron precipitation at the cusp. In: Proceedings of the 1st Meeting Workshop: Magnetic Reconnection at the Magnetopause and Aurora Dynamics; Verlag der Osterreichischen Akademie der Wissenschaften. 1996. p. 177-183
  53. 53. Namgaladze AA, Namgaladze AN, Volkov MA. Seasonal effects in the ionosphere-thermosphere response to the precipitation and field-aligned current variations in the cusp region. Annales Geophysicae. 1998;16(10):1283-1298. DOI: 10.1007/s00585-998-1283-3
  54. 54. Namgaladze AA, Yurik RYu. Plasmasphere state effect on the positive phase of the ionospheric storm. In: Proceedings of the XXVI Annual Apatity Seminar; Apatity, Russia. 2003. pp. 79-82
  55. 55. Namgaladze AA, Namgaladze AN, Yurik RY. Global modeling of the quiet and disturbed upper atmosphere. Physics and Chemistry of the Earth. 2000;25(5-6):533-536. DOI: 10.1016/S1464-1917(00)00071-4
  56. 56. Namgaladze AA, Foerster M, Yurik RY. Analysis of the positive ionospheric response to a moderate geomagnetic storm using a global numerical model. Annales Geophysicae. 2000;18(4):461-477. DOI: 10.1007/s00585-000-0461-8
  57. 57. Namgaladze AA, Namgaladze AN, Fadeeva Yu V, Goncharenko LP, Salah JE. Lower thermosphere and ionosphere behaviour during a strong magnetic storm of March 31, 2001: Modelling and comparison with the Millstone Hill incoherent scatter radar measurements. Procceedings of the MSTU. 2003;6(1):87-92
  58. 58. Namgaladze AA, Namgaladze AN, Martynenko OV, Doronina EN, Knyazeva MA, Zubova Yu V. Numerical modeling of the thermosphere, ionosphere and plasmasphere behaviour during the April 2002 magnetic storms. In: Proceedings of the XXVI Annual Apatity Seminar Preprint PGI. 2003. pp. 74-78
  59. 59. Förster M, Namgaladze AA, Yurik RY. Thermospheric composition changes deduced from geomagnetic storm modeling. Geophysical Research Letters. 1999;26(16):2625-2628. DOI: 10.1029/1999GL900514
  60. 60. Namgaladze AA, Klimenko MV, Klimenko VV, Zakharenkova IE. Physical mechanism and mathematical modeling of earthquake Ionospheric precursors registered in Total electron content. Geomagnetism and Aeronomy. 2009;49(2):252-262. DOI: 10.1134/S0016793209020169
  61. 61. Zolotov OV, Namgaladze AA, Zakharenkova IE, Martynenko OV, Shagimuratov II.Physical interpretation and mathematical simulation of ionospheric precursors of earthquakes at midlatitudes. Geomagnetism and Aeronomy. 2012;52(3):390-397. DOI: 10.1134/S0016793212030152
  62. 62. Karpov MI, Namgaladze AA, Zolotov OV. Modeling of Total electron content disturbances caused by electric currents between the earth and the ionosphere. Russian Journal of Physical Chemistry B. 2013;7(5):595-598. DOI: 10.1134/S1990793113050187
  63. 63. Namgaladze AA, Karpov MI. Conductivity and external electric currents in the global electric circuit. Russian Journal of Physical Chemistry В. 2015;9(4):754-757. DOI: 10.1134/S1990793115050231
  64. 64. Namgaladze АА. Earthquakes and global electrical circuit. Russian Journal of Physical Chemistry B. 2013;7(5):589-593. DOI: 10.1134/S1990793113050229
  65. 65. Nisbet JS. On the Construction and Use of the Pennsylvania State MKI Model. Ionosphere Research Laboratory; 1970 98 p

Written By

Namgaladze Aleksandr, Knyazeva Maria, Karpov Mikhail, Zolotov Oleg, Martynenko Oleg, Yurik Roman, Foerster Matthias and Prokhorov Boris

Submitted: 23 March 2017 Reviewed: 20 September 2017 Published: 20 December 2017