Articles | Volume 44, issue 2
https://doi.org/10.5194/angeo-44-921-2026
© Author(s) 2026. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
https://doi.org/10.5194/angeo-44-921-2026
© Author(s) 2026. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Non-Maxwellian electron distributions in the D region during artificial heating – Part 1: Model development and electron temperature
Margaretha Myrvang
CORRESPONDING AUTHOR
UiT The Arctic University of Norway, Department of Physics and Technology, Postboks 6050 Langnes, 9037 Tromsø, Norway
Björn Gustavsson
UiT The Arctic University of Norway, Department of Physics and Technology, Postboks 6050 Langnes, 9037 Tromsø, Norway
Related authors
Margaretha Myrvang and Björn Johan Gustavsson
EGUsphere, https://doi.org/10.5194/egusphere-2026-1118, https://doi.org/10.5194/egusphere-2026-1118, 2026
This preprint is open for discussion and under review for Annales Geophysicae (ANGEO).
Short summary
Short summary
This is the second of two papers. The first paper demonstrate that the electron distribution becomes non-Maxwellian during radio wave heating. Electron cooling rates are macroscopic properties of the electron gas, where electrons transfer some of their kinetic energy through collisions with neutrals and ions, and can be determined by integrating over the electron distribution. Any changes in the electron distribution affect the cooling rates.
Etienne Gavazzi, Andres Spicher, Björn Gustavsson, Juha Vierinen, James Clemmons, Robert Pfaff, and Douglas Rowland
Ann. Geophys., 44, 903–920, https://doi.org/10.5194/angeo-44-903-2026, https://doi.org/10.5194/angeo-44-903-2026, 2026
Short summary
Short summary
Auroras are caused by energetic electrons entering the upper atmosphere. For the smallest and most dynamic auroras, scientists think electrons are accelerated by waves (called Alfvén waves) thousands of kilometers above Earth. In this paper, we analyse data from a rocket that flew through auroras, applying existing and new techniques to estimate where the acceleration took place. Our results match theory, and we show how they can be used to study conditions in the near-Earth space environment.
Devin Huyghebaert, Juha Vierinen, Björn Gustavsson, Ralph Latteck, Toralf Renkwitz, Marius Zecha, Claudia C. Stephan, J. Federico Conte, Daniel Kastinen, Johan Kero, and Jorge L. Chau
Atmos. Meas. Tech., 19, 4277–4292, https://doi.org/10.5194/amt-19-4277-2026, https://doi.org/10.5194/amt-19-4277-2026, 2026
Short summary
Short summary
The phenomena of meteors occurs at altitudes of 60–120 km and can be used to measure the neutral atmosphere. We use a large high power radar system in Norway (Middle Atmosphere Alomar Radar System (MAARSY)) to measure the meteors and determine changes to the atmospheric density between the years of 2016–2023 at altitudes of 85–115 km. The same time period between years are compared, minimizing changes to the measurements due to factors other than the atmosphere.
Juha Vierinen, Dabrowka Knach, Jorge Luis Chau, Gerd Baumgarten, Devin Huyghebaert, Matthias Clahsen, Nico Pfeffer, Toralf Renkwitz, Robin Wing, Kenneth Obenberger, Björn Gustavsson, and Daniel Kastinen
EGUsphere, https://doi.org/10.5194/egusphere-2026-2857, https://doi.org/10.5194/egusphere-2026-2857, 2026
Short summary
Short summary
We observed the re-entry of a Falcon 9 rocket upper stage over Europe using a network of cameras and radars normally used to study meteors. The measurements showed how hot ionized gas formed around the debris during hypervelocity entry, producing strong radio echoes high in the atmosphere. As the number of satellites and rocket launches increases, these observations can help improve monitoring of space waste re-entering the atmosphere and its possible environmental effects.
Theresa Rexer, Björn Gustavsson, Juha Vierinen, Andres Spicher, Devin Ray Huyghebaert, Andreas Kvammen, Robert Gillies, and Asti Bhatt
Geosci. Instrum. Method. Data Syst., 15, 127–139, https://doi.org/10.5194/gi-15-127-2026, https://doi.org/10.5194/gi-15-127-2026, 2026
Short summary
Short summary
We present a second-level calibration method for electron density measurements from multi-beam incoherent scatter radars. It is based on the well-known Flat field correction method used in imaging and photography. The method improves data quality and useability as it accounts for unaccounted and unpredictable variations in the radar system. This is valuable for studies where inter-beam calibration is important such as studies of polar cap patches, plasma irregularities and turbulence.
Margaretha Myrvang and Björn Johan Gustavsson
EGUsphere, https://doi.org/10.5194/egusphere-2026-1118, https://doi.org/10.5194/egusphere-2026-1118, 2026
This preprint is open for discussion and under review for Annales Geophysicae (ANGEO).
Short summary
Short summary
This is the second of two papers. The first paper demonstrate that the electron distribution becomes non-Maxwellian during radio wave heating. Electron cooling rates are macroscopic properties of the electron gas, where electrons transfer some of their kinetic energy through collisions with neutrals and ions, and can be determined by integrating over the electron distribution. Any changes in the electron distribution affect the cooling rates.
Oliver Stalder, Björn Gustavsson, and Ilkka Virtanen
Ann. Geophys., 44, 123–135, https://doi.org/10.5194/angeo-44-123-2026, https://doi.org/10.5194/angeo-44-123-2026, 2026
Short summary
Short summary
The rapid changes in ion composition during auroral are dynamically modeled by integrating the coupled continuity equations for 15 ionospheric species. The effect of the ionospheric variation on the inversion of incoherent scatter radar (ISR) electron density profiles to differential energy spectra of precipitating electrons is studied. A systematic overestimation at high electron energies can be removed using a dynamic model. Comparisons are made with static and steady-state ionospheric models.
Kian Sartipzadeh, Andreas Kvammen, Björn Gustavsson, Njål Gulbrandsen, Magnar G. Johnsen, Devin Huyghebaert, and Juha Vierinen
Ann. Geophys., 44, 85–107, https://doi.org/10.5194/angeo-44-85-2026, https://doi.org/10.5194/angeo-44-85-2026, 2026
Short summary
Short summary
Knowledge of the charged environment in the upper atmosphere is essential for understanding space weather effects on satellites and radio communication. This environment is difficult to estimate at high latitudes, where aurora cause strong variability. We developed an artificial intelligence model to estimate this environment continuously. Our results show that the model provides reliable estimates even during auroral activity, improving monitoring of the polar upper atmosphere.
Etienne Gavazzi, Andres Spicher, Björn Gustavsson, James Clemmons, Robert Pfaff, and Douglas Rowland
Ann. Geophys., 44, 1–15, https://doi.org/10.5194/angeo-44-1-2026, https://doi.org/10.5194/angeo-44-1-2026, 2026
Short summary
Short summary
Auroral precipitation refers to energetic particles that come down into the upper part of our atmosphere, the ionosphere. There, they collide with atoms and molecules and transfer some of their energy, causing aurora. The most rapid time-variation of this energy deposition and its consequences on the ionosphere are not fully understood. We show here that one can use a new model to study auroral precipitation on sub-second timescales and advance our understanding about small-scale dynamic aurora.
Devin Huyghebaert, Björn Gustavsson, Juha Vierinen, Andreas Kvammen, Matthew Zettergren, John Swoboda, Ilkka Virtanen, Spencer M. Hatch, and Karl M. Laundal
Ann. Geophys., 43, 99–113, https://doi.org/10.5194/angeo-43-99-2025, https://doi.org/10.5194/angeo-43-99-2025, 2025
Short summary
Short summary
The EISCAT_3D radar is a new ionospheric radar under construction in the Fennoscandia region. The radar will make measurements of plasma characteristics at altitudes above approximately 60 km. The capability of the system to make these measurements at spatial scales of less than 100 m using multiple digitised signals from each of the radar antenna panels is highlighted. There are many ionospheric small-scale processes that will be further resolved using the techniques discussed here.
Yoshimasa Tanaka, Yasunobu Ogawa, Akira Kadokura, Takehiko Aso, Björn Gustavsson, Urban Brändström, Tima Sergienko, Genta Ueno, and Satoko Saita
Ann. Geophys., 42, 179–190, https://doi.org/10.5194/angeo-42-179-2024, https://doi.org/10.5194/angeo-42-179-2024, 2024
Short summary
Short summary
We present via simulation how useful monochromatic images taken by a multi-point imager network are for auroral research in the EISCAT_3D project. We apply the generalized-aurora computed tomography (G-ACT) to modeled multiple auroral images and ionospheric electron density data. It is demonstrated that G-ACT provides better reconstruction results than the normal ACT and can interpolate ionospheric electron density at a much higher spatial resolution than observed by the EISCAT_3D radar.
Thomas B. Leyser, Tima Sergienko, Urban Brändström, Björn Gustavsson, and Michael T. Rietveld
Ann. Geophys., 41, 589–600, https://doi.org/10.5194/angeo-41-589-2023, https://doi.org/10.5194/angeo-41-589-2023, 2023
Short summary
Short summary
Powerful radio waves transmitted into the ionosphere from the ground were used to study electron energization in the pumped ionospheric plasma turbulence, by detecting optical emissions from atomic oxygen. Our results obtained with the EISCAT (European Incoherent Scatter Scientific Association) facilities in northern Norway and optical detection with the ALIS (Auroral Large Imaging System) in northern Sweden suggest that long-wavelength upper hybrid waves are important in accelerating electrons.
Johann Stamm, Juha Vierinen, Björn Gustavsson, and Andres Spicher
Ann. Geophys., 41, 55–67, https://doi.org/10.5194/angeo-41-55-2023, https://doi.org/10.5194/angeo-41-55-2023, 2023
Short summary
Short summary
The study of some ionospheric events benefit from the knowledge of how the physics varies over a volume and over time. Examples are studies of aurora or energy deposition. With EISCAT3D, measurements of ion velocity vectors in a volume will be possible for the first time. We present a technique that uses a set of such measurements to estimate electric field and neutral wind. The technique relies on adding restrictions to the estimates. We successfully consider restrictions based on physics.
Daniel K. Whiter, Noora Partamies, Björn Gustavsson, and Kirsti Kauristie
Ann. Geophys., 41, 1–12, https://doi.org/10.5194/angeo-41-1-2023, https://doi.org/10.5194/angeo-41-1-2023, 2023
Short summary
Short summary
We measured the height of green and blue aurorae using thousands of camera images recorded over a 7-year period. Both colours are typically brightest at about 114 km altitude. When they peak at higher altitudes the blue aurora is usually higher than the green aurora. This information will help other studies which need an estimate of the auroral height. We used a computer model to explain our observations and to investigate how the green aurora is produced.
Mizuki Fukizawa, Takeshi Sakanoi, Yoshimasa Tanaka, Yasunobu Ogawa, Keisuke Hosokawa, Björn Gustavsson, Kirsti Kauristie, Alexander Kozlovsky, Tero Raita, Urban Brändström, and Tima Sergienko
Ann. Geophys., 40, 475–484, https://doi.org/10.5194/angeo-40-475-2022, https://doi.org/10.5194/angeo-40-475-2022, 2022
Short summary
Short summary
The pulsating auroral generation mechanism has been investigated by observing precipitating electrons using rockets or satellites. However, it is difficult for such observations to distinguish temporal changes from spatial ones. In this study, we reconstructed the horizontal 2-D distribution of precipitating electrons using only auroral images. The 3-D aurora structure was also reconstructed. We found that there were both spatial and temporal changes in the precipitating electron energy.
Johann Stamm, Juha Vierinen, and Björn Gustavsson
Ann. Geophys., 39, 961–974, https://doi.org/10.5194/angeo-39-961-2021, https://doi.org/10.5194/angeo-39-961-2021, 2021
Short summary
Short summary
Measurements of the electric field and neutral wind in the ionosphere are important for understanding energy flows or electric currents. With incoherent scatter radars (ISRs), we can measure the velocity of the ions, which depends on both the electrical field and the neutral wind. In this paper, we investigate methods to use ISR data to find reasonable values for both parameters. We find that electric field can be well measured down to 125 km height and neutral wind below this height.
Cited articles
Brunger, M. J. and Buckman, S. J.: Electron–molecule scattering cross-sections. I. Experimental techniques and data for diatomic molecules, Phys. Rep., 357, 215–458, 2002. a
Carleton, N. and Megill, L. R.: Electron energy distribution in slightly ionized air under the influence of electric and magnetic fields, Phys. Rev., 126, 2089, https://doi.org/10.1103/PhysRev.126.2089, 1962. a
Chapman, S. and Cowling, T. G.: The mathematical theory of non-uniform gases: an account of the kinetic theory of viscosity, thermal conduction and diffusion in gases, Cambridge University Press, New York, 1939. a
Dehmel, R., Fineman, M., and Miller, D.: Angular scattering of low-energy electrons by atomic and molecular oxygen, argon, and helium, Phys. Rev. A, 13, 115, https://doi.org/10.1103/PhysRevA.13.115, 1976. a
Gustavsson, B.: Time-Dependent Electron Transport I: Modelling of Supra-Thermal Electron Bursts Modulated at 5–10 Hz With Implications for Flickering Aurora, J. Geophys. Res.-Space, 127, e2019JA027608, https://doi.org/10.1029/2019JA027608, 2022. a
Gustavsson, B., Sergienko, T., Haggstrom, I., Honary, F., and Aso, T.: Simulation of high energy tail of electron distribution function, Advances in polar upper atmosphere research, 18, 1–9, 2004. a
Hedin, A. E.: Extension of the MSIS thermosphere model into the middle and lower atmosphere, J. Geophys. Res.-Space, 96, 1159–1172, 1991. a
Holstein, T.: Energy distribution of electrons in high frequency gas discharges, Phys. Rev., 70, 367, https://doi.org/10.1103/PhysRev.70.367, 1946. a
Kero, A., Bösinger, T., Pollari, P., Turunen, E., and Rietveld, M.: First EISCAT measurement of electron-gas temperature in the artificially heated D-region ionosphere, Ann. Geophys., 18, 1210–1215, https://doi.org/10.1007/s00585-000-1210-8, 2000. a
Milikh, G. and Dimant, Y. S.: Model of anomalous electron heating in the E region: 2. Detailed numerical modeling, J. Geophys. Res.-Space, 108, https://doi.org/10.1029/2002JA009527, 2003. a
Myrvang, M. and Gustavsson, B.: mmy002/Non-Maxwellian-distributions-and-cooling-rates-during-HF-ionospheric-heating: Code and data for non-Maxwellian during HF ionospheric heating (Version v1.0.0), Zenodo [software], https://doi.org/10.5281/zenodo.22069113, 2026a. a
Myrvang, M. and Gustavsson, B. J.: Non-Maxwellian electron distributions in the D region during artificial heating – Part 2: Electron cooling rates, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2026-1118, 2026b. a
Onda, K.: Rotational excitation of molecular nitrogen by electron impact, J. Phys. Soc. Jpn., 54, 4544–4554, 1985. a
Pavlov, A.: New electron energy transfer rates for vibrational excitation of N2, Ann. Geophys., 16, 176–182, https://doi.org/10.1007/s00585-998-0176-9, 1998a. a
Pavlov, A.: New electron energy transfer and cooling rates by excitation of O2, Ann. Geophys., 16, 1007–1013, https://doi.org/10.1007/s00585-998-1007-8, 1998b. a
Pavlov, A. and Berrington, K.: Cooling rate of thermal electrons by electron impact excitation of fine structure levels of atomic oxygen, Ann. Geophys., 17, 919–924, https://doi.org/10.1007/s00585-999-0919-2, 1999. a
Phelps, A. and Pack, J.: Electron collision frequencies in nitrogen and in the lower ionosphere, Phys. Rev. Lett., 3, 340, https://doi.org/10.1103/PhysRevLett.3.340, 1959. a, b
Picone, J., Hedin, A., Drob, D. P., and Aikin, A.: NRLMSISE-00 empirical model of the atmosphere: Statistical comparisons and scientific issues, J. Geophys. Res.-Space, 107, SIA–15, 2002. a
Schunk, R. W., Nagy, A. F., and Nagy, A.: Ionospheres: Physics, plasma physics, and chemistry, Cambridge University Press, ISBN 978-1-108-46210-5, 2018. a
Shoucri, M. M., Morales, G., and Maggs, J.: Ohmic heating of the polar F region by HF pulses, J. Geophys. Res.-Space, 89, 2907–2917, 1984. a
Turunen, E., Matveinen, H., Tolvanen, J., and Ranta, H.: STEP Handbook of Ionospheric Models, D-region ion chemistry model, in: STEP Handbook of Ionospheric Models, edited by: Schunk, R. W., SCOSTEP Secretariat, Boulder, Colorado, 1–25, OCLC 36598271, https://doi.org/10.1016/S1364-6826(97)00103-X, 1996. a
Viehland, L. A. and Johnsen, R.: Velocity distribution functions for ions drifting in helium and cross section for reaction of with N2 (v=0), J. Chem. Phys., 149, https://doi.org/10.1063/1.5033426, 2018. a
Short summary
This paper investigate how radio wave heating affect the electron distribution in the D region. Radio wave heating can lead to absorption of radio wave energy by electrons, thereby increasing the temperature. Electrons are cooled by inelastic collisions with the neutrals atmosphere, affecting the electron distribution since electrons lose energy by exciting different states in neutrals. Thus, electrons are redistributed to lower energies, which changes the shape of the electron distribution.
This paper investigate how radio wave heating affect the electron distribution in the D region....