Articles | Volume 17, issue 16
https://doi.org/10.5194/gmd-17-6401-2024
© Author(s) 2024. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
https://doi.org/10.5194/gmd-17-6401-2024
© Author(s) 2024. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Development and technical paper
| 
29 Aug 2024
Development and technical paper |
| 29 Aug 2024
| 29 Aug 2024
Physics-motivated cell-octree adaptive mesh refinement in the Vlasiator 5.3 global hybrid-Vlasov code
Department of Physics, University of Helsinki, Helsinki, Finland
Markus Battarbee
Department of Physics, University of Helsinki, Helsinki, Finland
Yann Pfau-Kempf
Department of Physics, University of Helsinki, Helsinki, Finland
Minna Palmroth
Department of Physics, University of Helsinki, Helsinki, Finland
Finnish Meteorological Institute, Space and Earth Observation Centre, Helsinki, Finland
Related authors
Urs Ganse, Yann Pfau-Kempf, Hongyang Zhou, Liisa Juusola, Abiyot Workayehu, Fasil Kebede, Konstantinos Papadakis, Maxime Grandin, Markku Alho, Markus Battarbee, Maxime Dubart, Leo Kotipalo, Arnaud Lalagüe, Jonas Suni, Konstantinos Horaites, and Minna Palmroth
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2024-101, https://doi.org/10.5194/gmd-2024-101, 2024
Revised manuscript accepted for GMD
Short summary
Short summary
Vlasiator is a kinetic space-plasma model that simulates the behaviour of plasma, solar wind and magnetic fields in near-Earth space. So far, these simulations had been run without any interaction wtih the ionosphere, the uppermost layer of Earth's atmosphere. In this manuscript, we present the new methods that add an ionospheric electrodynamics model to Vlasiator, coupling it with the existing methods and presenting new simulation results of how space Plasma and Earth's ionosphere interact.
Tuomas Häkkilä, Maxime Grandin, Markus Battarbee, Monika E. Szeląg, Markku Alho, Leo Kotipalo, Niilo Kalakoski, Pekka T. Verronen, and Minna Palmroth
Ann. Geophys. Discuss., https://doi.org/10.5194/angeo-2024-7, https://doi.org/10.5194/angeo-2024-7, 2024
Revised manuscript under review for ANGEO
Short summary
Short summary
We study the atmospheric impact of auroral electron precipitation, by the novel combination of both magnetospheric and atmospheric modelling. We first simulate fluxes of auroral electrons, and then use these fluxes to model their atmospheric impact. We find an increase of up to 200 % in thermospheric odd nitrogen, and a corresponding decrease in stratospheric ozone of around 0.7 %. The produced auroral electron precipitation is realistic, and shows the potential for future studies.
Markku Alho, Giulia Cozzani, Ivan Zaitsev, Fasil Tesema Kebede, Urs Ganse, Markus Battarbee, Maarja Bussov, Maxime Dubart, Sanni Hoilijoki, Leo Kotipalo, Konstantinos Papadakis, Yann Pfau-Kempf, Jonas Suni, Vertti Tarvus, Abiyot Workayehu, Hongyang Zhou, and Minna Palmroth
Ann. Geophys., 42, 145–161, https://doi.org/10.5194/angeo-42-145-2024, https://doi.org/10.5194/angeo-42-145-2024, 2024
Short summary
Short summary
Magnetic reconnection is one of the main processes for energy conversion and plasma transport in space plasma physics, associated with plasma entry into the magnetosphere of Earth and Earth’s substorm cycle. Global modelling of these plasma processes enables us to understand the magnetospheric system in detail. However, finding sites of active reconnection from large simulation datasets can be challenging, and this paper develops tools to find magnetic topologies related to reconnection.
Yann Pfau-Kempf, Konstantinos Papadakis, Markku Alho, Markus Battarbee, Giulia Cozzani, Lauri Pänkäläinen, Urs Ganse, Fasil Kebede, Jonas Suni, Konstantinos Horaites, Maxime Grandin, and Minna Palmroth
Ann. Geophys. Discuss., https://doi.org/10.5194/angeo-2024-26, https://doi.org/10.5194/angeo-2024-26, 2024
Preprint under review for ANGEO
Short summary
Short summary
Flux ropes are peculiar structures of twisted magnetic field occurring in many regions of space, near Earth and other planets, at the Sun, and in astrophysical objects. We developed a new way of detecting flux ropes in large supercomputer simulations of near-Earth space and we use it to follow the evolution of flux ropes for long distances past the Earth in the flow direction. This will be useful in future studies as these flux ropes are involved in the transport of matter and energy in space.
Urs Ganse, Yann Pfau-Kempf, Hongyang Zhou, Liisa Juusola, Abiyot Workayehu, Fasil Kebede, Konstantinos Papadakis, Maxime Grandin, Markku Alho, Markus Battarbee, Maxime Dubart, Leo Kotipalo, Arnaud Lalagüe, Jonas Suni, Konstantinos Horaites, and Minna Palmroth
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2024-101, https://doi.org/10.5194/gmd-2024-101, 2024
Revised manuscript accepted for GMD
Short summary
Short summary
Vlasiator is a kinetic space-plasma model that simulates the behaviour of plasma, solar wind and magnetic fields in near-Earth space. So far, these simulations had been run without any interaction wtih the ionosphere, the uppermost layer of Earth's atmosphere. In this manuscript, we present the new methods that add an ionospheric electrodynamics model to Vlasiator, coupling it with the existing methods and presenting new simulation results of how space Plasma and Earth's ionosphere interact.
Tuomas Häkkilä, Maxime Grandin, Markus Battarbee, Monika E. Szeląg, Markku Alho, Leo Kotipalo, Niilo Kalakoski, Pekka T. Verronen, and Minna Palmroth
Ann. Geophys. Discuss., https://doi.org/10.5194/angeo-2024-7, https://doi.org/10.5194/angeo-2024-7, 2024
Revised manuscript under review for ANGEO
Short summary
Short summary
We study the atmospheric impact of auroral electron precipitation, by the novel combination of both magnetospheric and atmospheric modelling. We first simulate fluxes of auroral electrons, and then use these fluxes to model their atmospheric impact. We find an increase of up to 200 % in thermospheric odd nitrogen, and a corresponding decrease in stratospheric ozone of around 0.7 %. The produced auroral electron precipitation is realistic, and shows the potential for future studies.
Markku Alho, Giulia Cozzani, Ivan Zaitsev, Fasil Tesema Kebede, Urs Ganse, Markus Battarbee, Maarja Bussov, Maxime Dubart, Sanni Hoilijoki, Leo Kotipalo, Konstantinos Papadakis, Yann Pfau-Kempf, Jonas Suni, Vertti Tarvus, Abiyot Workayehu, Hongyang Zhou, and Minna Palmroth
Ann. Geophys., 42, 145–161, https://doi.org/10.5194/angeo-42-145-2024, https://doi.org/10.5194/angeo-42-145-2024, 2024
Short summary
Short summary
Magnetic reconnection is one of the main processes for energy conversion and plasma transport in space plasma physics, associated with plasma entry into the magnetosphere of Earth and Earth’s substorm cycle. Global modelling of these plasma processes enables us to understand the magnetospheric system in detail. However, finding sites of active reconnection from large simulation datasets can be challenging, and this paper develops tools to find magnetic topologies related to reconnection.
Jonas Suni, Minna Palmroth, Lucile Turc, Markus Battarbee, Giulia Cozzani, Maxime Dubart, Urs Ganse, Harriet George, Evgeny Gordeev, Konstantinos Papadakis, Yann Pfau-Kempf, Vertti Tarvus, Fasil Tesema, and Hongyang Zhou
Ann. Geophys., 41, 551–568, https://doi.org/10.5194/angeo-41-551-2023, https://doi.org/10.5194/angeo-41-551-2023, 2023
Short summary
Short summary
Magnetosheath jets are structures of enhanced plasma density and/or velocity in a region of near-Earth space known as the magnetosheath. When they propagate towards the Earth, these jets can disturb the Earth's magnetic field and cause hazards for satellites. In this study, we use a simulation called Vlasiator to model near-Earth space and investigate jets using case studies and statistical analysis. We find that jets that propagate towards the Earth are different from jets that do not.
Konstantinos Papadakis, Yann Pfau-Kempf, Urs Ganse, Markus Battarbee, Markku Alho, Maxime Grandin, Maxime Dubart, Lucile Turc, Hongyang Zhou, Konstantinos Horaites, Ivan Zaitsev, Giulia Cozzani, Maarja Bussov, Evgeny Gordeev, Fasil Tesema, Harriet George, Jonas Suni, Vertti Tarvus, and Minna Palmroth
Geosci. Model Dev., 15, 7903–7912, https://doi.org/10.5194/gmd-15-7903-2022, https://doi.org/10.5194/gmd-15-7903-2022, 2022
Short summary
Short summary
Vlasiator is a plasma simulation code that simulates the entire near-Earth space at a global scale. As 6D simulations require enormous amounts of computational resources, Vlasiator uses adaptive mesh refinement (AMR) to lighten the computational burden. However, due to Vlasiator’s grid topology, AMR simulations suffer from grid aliasing artifacts that affect the global results. In this work, we present and evaluate the performance of a mechanism for alleviating those artifacts.
Vertti Tarvus, Lucile Turc, Markus Battarbee, Jonas Suni, Xóchitl Blanco-Cano, Urs Ganse, Yann Pfau-Kempf, Markku Alho, Maxime Dubart, Maxime Grandin, Andreas Johlander, Konstantinos Papadakis, and Minna Palmroth
Ann. Geophys., 39, 911–928, https://doi.org/10.5194/angeo-39-911-2021, https://doi.org/10.5194/angeo-39-911-2021, 2021
Short summary
Short summary
We use simulations of Earth's magnetosphere and study the formation of transient wave structures in the region where the solar wind first interacts with the magnetosphere. These transients move earthward and play a part in the solar wind–magnetosphere interaction. We show that the transients are a common feature and their properties are altered as they move earthward, including an increase in temperature, decrease in solar wind speed and an alteration in their propagation properties.
Andrei Runov, Maxime Grandin, Minna Palmroth, Markus Battarbee, Urs Ganse, Heli Hietala, Sanni Hoilijoki, Emilia Kilpua, Yann Pfau-Kempf, Sergio Toledo-Redondo, Lucile Turc, and Drew Turner
Ann. Geophys., 39, 599–612, https://doi.org/10.5194/angeo-39-599-2021, https://doi.org/10.5194/angeo-39-599-2021, 2021
Short summary
Short summary
In collisionless systems like space plasma, particle velocity distributions contain fingerprints of ongoing physical processes. However, it is challenging to decode this information from observations. We used hybrid-Vlasov simulations to obtain ion velocity distribution functions at different locations and at different stages of the Earth's magnetosphere dynamics. The obtained distributions provide valuable examples that may be directly compared with observations by satellites in space.
Minna Palmroth, Savvas Raptis, Jonas Suni, Tomas Karlsson, Lucile Turc, Andreas Johlander, Urs Ganse, Yann Pfau-Kempf, Xochitl Blanco-Cano, Mojtaba Akhavan-Tafti, Markus Battarbee, Maxime Dubart, Maxime Grandin, Vertti Tarvus, and Adnane Osmane
Ann. Geophys., 39, 289–308, https://doi.org/10.5194/angeo-39-289-2021, https://doi.org/10.5194/angeo-39-289-2021, 2021
Short summary
Short summary
Magnetosheath jets are high-velocity features within the Earth's turbulent magnetosheath, separating the Earth's magnetic domain from the solar wind. The characteristics of the jets are difficult to assess statistically as a function of their lifetime because normally spacecraft observe them only at one position within the magnetosheath. This study first confirms the accuracy of the model used, Vlasiator, by comparing it to MMS spacecraft, and then carries out the first jet lifetime statistics.
Minna Palmroth, Maxime Grandin, Theodoros Sarris, Eelco Doornbos, Stelios Tourgaidis, Anita Aikio, Stephan Buchert, Mark A. Clilverd, Iannis Dandouras, Roderick Heelis, Alex Hoffmann, Nickolay Ivchenko, Guram Kervalishvili, David J. Knudsen, Anna Kotova, Han-Li Liu, David M. Malaspina, Günther March, Aurélie Marchaudon, Octav Marghitu, Tomoko Matsuo, Wojciech J. Miloch, Therese Moretto-Jørgensen, Dimitris Mpaloukidis, Nils Olsen, Konstantinos Papadakis, Robert Pfaff, Panagiotis Pirnaris, Christian Siemes, Claudia Stolle, Jonas Suni, Jose van den IJssel, Pekka T. Verronen, Pieter Visser, and Masatoshi Yamauchi
Ann. Geophys., 39, 189–237, https://doi.org/10.5194/angeo-39-189-2021, https://doi.org/10.5194/angeo-39-189-2021, 2021
Short summary
Short summary
This is a review paper that summarises the current understanding of the lower thermosphere–ionosphere (LTI) in terms of measurements and modelling. The LTI is the transition region between space and the atmosphere and as such of tremendous importance to both the domains of space and atmosphere. The paper also serves as the background for European Space Agency Earth Explorer 10 candidate mission Daedalus.
Markus Battarbee, Thiago Brito, Markku Alho, Yann Pfau-Kempf, Maxime Grandin, Urs Ganse, Konstantinos Papadakis, Andreas Johlander, Lucile Turc, Maxime Dubart, and Minna Palmroth
Ann. Geophys., 39, 85–103, https://doi.org/10.5194/angeo-39-85-2021, https://doi.org/10.5194/angeo-39-85-2021, 2021
Short summary
Short summary
We investigate local acceleration dynamics of electrons with a new numerical simulation method, which is an extension of a world-leading kinetic plasma simulation. We describe how large supercomputer simulations can be used to initialize the electron simulations and show numerical stability for the electron method. We show that features of our simulated electrons match observations from Earth's magnetic tail region.
Maxime Dubart, Urs Ganse, Adnane Osmane, Andreas Johlander, Markus Battarbee, Maxime Grandin, Yann Pfau-Kempf, Lucile Turc, and Minna Palmroth
Ann. Geophys., 38, 1283–1298, https://doi.org/10.5194/angeo-38-1283-2020, https://doi.org/10.5194/angeo-38-1283-2020, 2020
Short summary
Short summary
Plasma waves are ubiquitous in the Earth's magnetosphere. They are responsible for many energetic processes happening in Earth's atmosphere, such as auroras. In order to understand these processes, thorough investigations of these waves are needed. We use a state-of-the-art numerical model to do so. Here we investigate the impact of different spatial resolutions in the model on these waves in order to improve in the future the model without wasting computational resources.
Markus Battarbee, Xóchitl Blanco-Cano, Lucile Turc, Primož Kajdič, Andreas Johlander, Vertti Tarvus, Stephen Fuselier, Karlheinz Trattner, Markku Alho, Thiago Brito, Urs Ganse, Yann Pfau-Kempf, Mojtaba Akhavan-Tafti, Tomas Karlsson, Savvas Raptis, Maxime Dubart, Maxime Grandin, Jonas Suni, and Minna Palmroth
Ann. Geophys., 38, 1081–1099, https://doi.org/10.5194/angeo-38-1081-2020, https://doi.org/10.5194/angeo-38-1081-2020, 2020
Short summary
Short summary
We investigate the dynamics of helium in the foreshock, a part of near-Earth space found upstream of the Earth's bow shock. We show how the second most common ion in interplanetary space reacts strongly to plasma waves found in the foreshock. Spacecraft observations and supercomputer simulations both give us a new understanding of the foreshock edge and how to interpret future observations.
Lucile Turc, Vertti Tarvus, Andrew P. Dimmock, Markus Battarbee, Urs Ganse, Andreas Johlander, Maxime Grandin, Yann Pfau-Kempf, Maxime Dubart, and Minna Palmroth
Ann. Geophys., 38, 1045–1062, https://doi.org/10.5194/angeo-38-1045-2020, https://doi.org/10.5194/angeo-38-1045-2020, 2020
Short summary
Short summary
Using global computer simulations, we study properties of the magnetosheath, the region of near-Earth space where the stream of particles originating from the Sun, the solar wind, is slowed down and deflected around the Earth's magnetic field. One of our main findings is that even for idealised solar wind conditions as used in our model, the magnetosheath density shows large-scale spatial and temporal variation in the so-called quasi-parallel magnetosheath, causing varying levels of asymmetry.
Harriet George, Emilia Kilpua, Adnane Osmane, Timo Asikainen, Milla M. H. Kalliokoski, Craig J. Rodger, Stepan Dubyagin, and Minna Palmroth
Ann. Geophys., 38, 931–951, https://doi.org/10.5194/angeo-38-931-2020, https://doi.org/10.5194/angeo-38-931-2020, 2020
Short summary
Short summary
We compared trapped outer radiation belt electron fluxes to high-latitude precipitating electron fluxes during two interplanetary coronal mass ejections (ICMEs) with opposite magnetic cloud rotation. The electron response had many similarities and differences between the two events, indicating that different acceleration mechanisms acted. Van Allen Probe data were used for trapped electron flux measurements, and Polar Operational Environmental Satellites were used for precipitating flux data.
Milla M. H. Kalliokoski, Emilia K. J. Kilpua, Adnane Osmane, Drew L. Turner, Allison N. Jaynes, Lucile Turc, Harriet George, and Minna Palmroth
Ann. Geophys., 38, 683–701, https://doi.org/10.5194/angeo-38-683-2020, https://doi.org/10.5194/angeo-38-683-2020, 2020
Short summary
Short summary
We present a comprehensive statistical study of the response of the Earth's space environment in sheath regions prior to interplanetary coronal mass ejections. The inner magnetospheric wave activity is enhanced in sheath regions, and the sheaths cause significant changes to the outer radiation belt electron fluxes over short timescales. We also show that non-geoeffective sheaths can result in a significant response.
Markus Battarbee, Urs Ganse, Yann Pfau-Kempf, Lucile Turc, Thiago Brito, Maxime Grandin, Tuomas Koskela, and Minna Palmroth
Ann. Geophys., 38, 625–643, https://doi.org/10.5194/angeo-38-625-2020, https://doi.org/10.5194/angeo-38-625-2020, 2020
Short summary
Short summary
The structure and medium-scale dynamics of Earth's bow shock and how charged solar wind particles are reflected by it are studied in order to better understand space weather effects. We use advanced supercomputer simulations to model the shock and reflected ions. We find that the thickness of the shock depends on solar wind conditions but also has small-scale variations. Charged particle reflection is shown to be non-localized. Magnetic fields are important for ion reflection.
Theodoros E. Sarris, Elsayed R. Talaat, Minna Palmroth, Iannis Dandouras, Errico Armandillo, Guram Kervalishvili, Stephan Buchert, Stylianos Tourgaidis, David M. Malaspina, Allison N. Jaynes, Nikolaos Paschalidis, John Sample, Jasper Halekas, Eelco Doornbos, Vaios Lappas, Therese Moretto Jørgensen, Claudia Stolle, Mark Clilverd, Qian Wu, Ingmar Sandberg, Panagiotis Pirnaris, and Anita Aikio
Geosci. Instrum. Method. Data Syst., 9, 153–191, https://doi.org/10.5194/gi-9-153-2020, https://doi.org/10.5194/gi-9-153-2020, 2020
Short summary
Short summary
Daedalus aims to measure the largely unexplored area between Eart's atmosphere and space, the Earth's
ignorosphere. Here, intriguing and complex processes govern the deposition and transport of energy. The aim is to quantify this energy by measuring effects caused by electrodynamic processes in this region. The concept is based on a mother satellite that carries a suite of instruments, along with smaller satellites carrying a subset of instruments that are released into the atmosphere.
Emilia Kilpua, Liisa Juusola, Maxime Grandin, Antti Kero, Stepan Dubyagin, Noora Partamies, Adnane Osmane, Harriet George, Milla Kalliokoski, Tero Raita, Timo Asikainen, and Minna Palmroth
Ann. Geophys., 38, 557–574, https://doi.org/10.5194/angeo-38-557-2020, https://doi.org/10.5194/angeo-38-557-2020, 2020
Short summary
Short summary
Coronal mass ejection sheaths and ejecta are key drivers of significant space weather storms, and they cause dramatic changes in radiation belt electron fluxes. Differences in precipitation of high-energy electrons from the belts to the upper atmosphere are thus expected. We investigate here differences in sheath- and ejecta-induced precipitation using the Finnish riometer (relative ionospheric opacity meter) chain.
Maxime Grandin, Markus Battarbee, Adnane Osmane, Urs Ganse, Yann Pfau-Kempf, Lucile Turc, Thiago Brito, Tuomas Koskela, Maxime Dubart, and Minna Palmroth
Ann. Geophys., 37, 791–806, https://doi.org/10.5194/angeo-37-791-2019, https://doi.org/10.5194/angeo-37-791-2019, 2019
Short summary
Short summary
When the terrestrial magnetic field is disturbed, particles from the near-Earth space can precipitate into the upper atmosphere. This work presents, for the first time, numerical simulations of proton precipitation in the energy range associated with the production of aurora (∼1–30 keV) using a global kinetic model of the near-Earth space: Vlasiator. We find that nightside proton precipitation can be regulated by the transition region between stretched and dipolar geomagnetic field lines.
Antti Lakka, Tuija I. Pulkkinen, Andrew P. Dimmock, Emilia Kilpua, Matti Ala-Lahti, Ilja Honkonen, Minna Palmroth, and Osku Raukunen
Ann. Geophys., 37, 561–579, https://doi.org/10.5194/angeo-37-561-2019, https://doi.org/10.5194/angeo-37-561-2019, 2019
Short summary
Short summary
We study how the Earth's space environment responds to two different amplitude interplanetary coronal mass ejection (ICME) events that occurred in 2012 and 2014 by using the GUMICS-4 global MHD model. We examine local and large-scale dynamics of the Earth's space environment and compare simulation results to in situ data. It is shown that during moderate driving simulation agrees well with the measurements; however, GMHD results should be interpreted cautiously during strong driving.
Liisa Juusola, Sanni Hoilijoki, Yann Pfau-Kempf, Urs Ganse, Riku Jarvinen, Markus Battarbee, Emilia Kilpua, Lucile Turc, and Minna Palmroth
Ann. Geophys., 36, 1183–1199, https://doi.org/10.5194/angeo-36-1183-2018, https://doi.org/10.5194/angeo-36-1183-2018, 2018
Short summary
Short summary
The solar wind interacts with the Earth’s magnetic field, forming a magnetosphere. On the night side solar wind stretches the magnetosphere into a long tail. A process called magnetic reconnection opens the magnetic field lines and reconnects them, accelerating particles to high energies. We study this in the magnetotail using a numerical simulation model of the Earth’s magnetosphere. We study the motion of the points where field lines reconnect and the fast flows driven by this process.
Minna Palmroth, Heli Hietala, Ferdinand Plaschke, Martin Archer, Tomas Karlsson, Xóchitl Blanco-Cano, David Sibeck, Primož Kajdič, Urs Ganse, Yann Pfau-Kempf, Markus Battarbee, and Lucile Turc
Ann. Geophys., 36, 1171–1182, https://doi.org/10.5194/angeo-36-1171-2018, https://doi.org/10.5194/angeo-36-1171-2018, 2018
Short summary
Short summary
Magnetosheath jets are high-velocity plasma structures that are commonly observed within the Earth's magnetosheath. Previously, they have mainly been investigated with spacecraft observations, which do not allow us to infer their spatial sizes, temporal evolution, or origin. This paper shows for the first time their dimensions, evolution, and origins within a simulation whose dimensions are directly comparable to the Earth's magnetosphere. The results are compared to previous observations.
Xochitl Blanco-Cano, Markus Battarbee, Lucile Turc, Andrew P. Dimmock, Emilia K. J. Kilpua, Sanni Hoilijoki, Urs Ganse, David G. Sibeck, Paul A. Cassak, Robert C. Fear, Riku Jarvinen, Liisa Juusola, Yann Pfau-Kempf, Rami Vainio, and Minna Palmroth
Ann. Geophys., 36, 1081–1097, https://doi.org/10.5194/angeo-36-1081-2018, https://doi.org/10.5194/angeo-36-1081-2018, 2018
Short summary
Short summary
We use the Vlasiator code to study the characteristics of transient structures that exist in the Earth's foreshock, i.e. upstream of the bow shock. The structures are cavitons and spontaneous hot flow anomalies (SHFAs). These transients can interact with the bow shock. We study the changes the shock suffers via this interaction. We also investigate ion distributions associated with the cavitons and SHFAs. A very important result is that the arrival of multiple SHFAs results in shock erosion.
Liisa Juusola, Yann Pfau-Kempf, Urs Ganse, Markus Battarbee, Thiago Brito, Maxime Grandin, Lucile Turc, and Minna Palmroth
Ann. Geophys., 36, 1027–1035, https://doi.org/10.5194/angeo-36-1027-2018, https://doi.org/10.5194/angeo-36-1027-2018, 2018
Short summary
Short summary
The Earth's magnetic field is shaped by the solar wind. On the dayside the field is compressed and on the nightside it is stretched as a long tail. The tail has been observed to occasionally undergo flapping motions, but the origin of these motions is not understood. We study the flapping using a numerical simulation of the near-Earth space. We present a possible explanation for how the flapping could be initiated by a passing disturbance and then maintained as a standing wave.
Minna Palmroth, Sanni Hoilijoki, Liisa Juusola, Tuija I. Pulkkinen, Heli Hietala, Yann Pfau-Kempf, Urs Ganse, Sebastian von Alfthan, Rami Vainio, and Michael Hesse
Ann. Geophys., 35, 1269–1274, https://doi.org/10.5194/angeo-35-1269-2017, https://doi.org/10.5194/angeo-35-1269-2017, 2017
Short summary
Short summary
Much like solar flares, substorms occurring within the Earth's magnetic domain are explosive events that cause vivid auroral displays. A decades-long debate exists to explain the substorm onset. We devise a simulation encompassing the entire near-Earth space and demonstrate that detailed modelling of magnetic reconnection explains the central substorm observations. Our results help to understand the unpredictable substorm process, which will significantly improve space weather forecasts.
Antti Lakka, Tuija I. Pulkkinen, Andrew P. Dimmock, Adnane Osmane, Ilja Honkonen, Minna Palmroth, and Pekka Janhunen
Ann. Geophys., 35, 907–922, https://doi.org/10.5194/angeo-35-907-2017, https://doi.org/10.5194/angeo-35-907-2017, 2017
Short summary
Short summary
We studied the impact on global MHD simulations from different simulation initialisation methods. While the global MHD code used is GUMICS-4 we conclude that the results might be generalisable to other codes as well. It is found that different initialisation methods affect the dynamics of the Earth's space environment by creating differences in momentum transport several hours afterwards. These differences may even grow as a response to rapid solar wind condition changes.
Yann Pfau-Kempf, Heli Hietala, Steve E. Milan, Liisa Juusola, Sanni Hoilijoki, Urs Ganse, Sebastian von Alfthan, and Minna Palmroth
Ann. Geophys., 34, 943–959, https://doi.org/10.5194/angeo-34-943-2016, https://doi.org/10.5194/angeo-34-943-2016, 2016
Short summary
Short summary
We have simulated the interaction of the solar wind – the charged particles and magnetic fields emitted by the Sun into space – with the magnetic field of the Earth. The solar wind flows supersonically and creates a shock when it encounters the obstacle formed by the geomagnetic field. We have identified a new chain of events which causes phenomena in the downstream region to eventually cause perturbations at the shock and even upstream. This is confirmed by ground and satellite observations.
P. T. Verronen, M. E. Andersson, A. Kero, C.-F. Enell, J. M. Wissing, E. R. Talaat, K. Kauristie, M. Palmroth, T. E. Sarris, and E. Armandillo
Ann. Geophys., 33, 381–394, https://doi.org/10.5194/angeo-33-381-2015, https://doi.org/10.5194/angeo-33-381-2015, 2015
Short summary
Short summary
Electron concentrations observed by EISCAT radars can be reasonable well represented using AIMOS v1.2 satellite-data-based ionization model and SIC D-region ion chemistry model. SIC-EISCAT difference varies from event to event, probably because the statistical nature of AIMOS ionization is not capturing all the spatio-temporal fine structure of electron precipitation. Below 90km, AIMOS overestimates electron ionization because of proton contamination of the satellite electron detectors.
D. Pokhotelov, S. von Alfthan, Y. Kempf, R. Vainio, H. E. J. Koskinen, and M. Palmroth
Ann. Geophys., 31, 2207–2212, https://doi.org/10.5194/angeo-31-2207-2013, https://doi.org/10.5194/angeo-31-2207-2013, 2013
A. T. Aikio, T. Pitkänen, I. Honkonen, M. Palmroth, and O. Amm
Ann. Geophys., 31, 1021–1034, https://doi.org/10.5194/angeo-31-1021-2013, https://doi.org/10.5194/angeo-31-1021-2013, 2013
Related subject area
Solar-terrestrial science
New routine NLTE15µmCool-E v1.0 for calculating the non-local thermodynamic equilibrium (non-LTE) CO2 15 µm cooling in general circulation models (GCMs) of Earth's atmosphere
The Vlasiator 5.2 Ionosphere – Coupling a magnetospheric hybrid-Vlasov simulation with a height-integrated ionosphere model
Daily INSOLation (DINSOL-v1.0): an intuitive tool for classrooms and specifying solar radiation boundary conditions
SSolar-GOA v1.0: a simple, fast, and accurate Spectral SOLAR radiative transfer model for clear skies
Application of CCM SOCOL-AERv2-BE to cosmogenic beryllium isotopes: description and validation for polar regions
UBER v1.0: a universal kinetic equation solver for radiation belts
Azimuthal averaging–reconstruction filtering techniques for finite-difference general circulation models in spherical geometry
Improved forecasting of thermospheric densities using multi-model ensembles
Accounting for anthropic energy flux of traffic in winter urban road surface temperature simulations with the TEB model
Calculations of the integral invariant coordinates I and L* in the magnetosphere and mapping of the regions where I is conserved, using a particle tracer (ptr3D v2.0), LANL*, SPENVIS, and IRBEM
A simple parameterization of the short-wave aerosol optical properties for surface direct and diffuse irradiances assessment in a numerical weather model
Verification of SpacePy's radial diffusion radiation belt model
LANL*V2.0: global modeling and validation
Alexander Kutepov and Artem Feofilov
Geosci. Model Dev., 17, 5331–5347, https://doi.org/10.5194/gmd-17-5331-2024, https://doi.org/10.5194/gmd-17-5331-2024, 2024
Short summary
Short summary
Infrared CO2 cooling of the middle and upper atmosphere is increasing. We developed a new routine for very fast and accurate calculations of this cooling in general circulation models. The new algorithm accounts for non-local thermodynamic equilibrium and is about 1000 times faster than the standard matrix algorithms. It is based on advanced techniques for non-equilibrium emission calculations in stellar atmospheres, which so far have not been used in Earth’s and planetary atmospheres.
Urs Ganse, Yann Pfau-Kempf, Hongyang Zhou, Liisa Juusola, Abiyot Workayehu, Fasil Kebede, Konstantinos Papadakis, Maxime Grandin, Markku Alho, Markus Battarbee, Maxime Dubart, Leo Kotipalo, Arnaud Lalagüe, Jonas Suni, Konstantinos Horaites, and Minna Palmroth
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2024-101, https://doi.org/10.5194/gmd-2024-101, 2024
Revised manuscript accepted for GMD
Short summary
Short summary
Vlasiator is a kinetic space-plasma model that simulates the behaviour of plasma, solar wind and magnetic fields in near-Earth space. So far, these simulations had been run without any interaction wtih the ionosphere, the uppermost layer of Earth's atmosphere. In this manuscript, we present the new methods that add an ionospheric electrodynamics model to Vlasiator, coupling it with the existing methods and presenting new simulation results of how space Plasma and Earth's ionosphere interact.
Emerson D. Oliveira
Geosci. Model Dev., 16, 2371–2390, https://doi.org/10.5194/gmd-16-2371-2023, https://doi.org/10.5194/gmd-16-2371-2023, 2023
Short summary
Short summary
Based on the Milankovitch cycle theory, the Daily INSOLation (DINSOL-v1.0) program simulates the incoming solar radiation at the top of the atmosphere, such as the PMIP boundary conditions. Still, users can simulate hypothetical cases by freely setting the Earth's orbital parameters. The program is recommended for educational purposes (from a user-friendly interface) or to prepare data for simplified climate models (from command lines). The program is supported on Linux and Windows.
Victoria Eugenia Cachorro, Juan Carlos Antuña-Sanchez, and Ángel Máximo de Frutos
Geosci. Model Dev., 15, 1689–1712, https://doi.org/10.5194/gmd-15-1689-2022, https://doi.org/10.5194/gmd-15-1689-2022, 2022
Short summary
Short summary
This work describes the features of a simple, fast, accurate, and physically based spectral radiative transfer model (SSolar-GOA) in the solar wavelength range under clear skies. The model is intended for a wide community of users for many different applications, was designed to be easily replicated, and has sufficient accuracy. The validation of the model was carried out through extensive comparison with simulated spectra from the LibRadtran and with direct and global spectral measurements.
Kseniia Golubenko, Eugene Rozanov, Gennady Kovaltsov, Ari-Pekka Leppänen, Timofei Sukhodolov, and Ilya Usoskin
Geosci. Model Dev., 14, 7605–7620, https://doi.org/10.5194/gmd-14-7605-2021, https://doi.org/10.5194/gmd-14-7605-2021, 2021
Short summary
Short summary
A new full 3-D time-dependent model, based on SOCOL-AERv2, of beryllium atmospheric production, transport, and deposition has been developed and validated using directly measured data. The model is recommended to be used in studies related to, e.g., atmospheric dynamical patterns, extreme solar particle storms, long-term solar activity reconstruction from cosmogenic proxy data, and solar–terrestrial relations.
Liheng Zheng, Lunjin Chen, Anthony A. Chan, Peng Wang, Zhiyang Xia, and Xu Liu
Geosci. Model Dev., 14, 5825–5842, https://doi.org/10.5194/gmd-14-5825-2021, https://doi.org/10.5194/gmd-14-5825-2021, 2021
Short summary
Short summary
Earth’s Van Allen belts are studied by solving particular kinds of equations that could be notoriously difficult when different physical processes are acting together. In this article, we describe a numerical code that can solve these equations with unprecedented freedom from the numerous restrictions of existing models, even the ones that no other can solve. The abilities of our code could mean a breakthrough in Van Allen belt studies from the diffusive into the non-diffusive transport regime.
Tong Dang, Binzheng Zhang, Jiuhou Lei, Wenbin Wang, Alan Burns, Han-li Liu, Kevin Pham, and Kareem A. Sorathia
Geosci. Model Dev., 14, 859–873, https://doi.org/10.5194/gmd-14-859-2021, https://doi.org/10.5194/gmd-14-859-2021, 2021
Short summary
Short summary
This paper describes a numerical treatment (ring average) to relax the time step in finite-difference schemes when using spherical and cylindrical coordinates with axis singularities. The ring average is used to develop a high-resolution thermosphere–ionosphere coupled community model. The technique is a significant improvement in space weather modeling capability, and it can also be adapted to more general finite-difference solvers for hyperbolic equations in spherical and polar geometries.
Sean Elvidge, Humberto C. Godinez, and Matthew J. Angling
Geosci. Model Dev., 9, 2279–2292, https://doi.org/10.5194/gmd-9-2279-2016, https://doi.org/10.5194/gmd-9-2279-2016, 2016
Short summary
Short summary
This paper presents the first known application of multi-model ensembles to the forecasting of the thermosphere. A multi-model ensemble (MME) is a method for combining different, independent models. The main advantage of using an MME is to reduce the effect of model errors and bias, since it is expected that the model errors will, at least partly, cancel. This paper shows that use of MMEs for forecasting thermospheric densities can reduce errors by 60 %.
A. Khalifa, M. Marchetti, L. Bouilloud, E. Martin, M. Bues, and K. Chancibaut
Geosci. Model Dev., 9, 547–565, https://doi.org/10.5194/gmd-9-547-2016, https://doi.org/10.5194/gmd-9-547-2016, 2016
Short summary
Short summary
An experimental study was conducted to quantify the anthropic energy flux of traffic impact on RST in the winter season. It indicated an RST increase by 1 °C to 3 °C with respect to the absence of traffic. Additional work was undertaken so as to evaluate to which extent an accurate description of traffic might improve the TEB numerical model when dedicated to RST simulations. Two approaches to traffic integration in this model were detailed and tested.
K. Konstantinidis and T. Sarris
Geosci. Model Dev., 8, 2967–2975, https://doi.org/10.5194/gmd-8-2967-2015, https://doi.org/10.5194/gmd-8-2967-2015, 2015
Short summary
Short summary
The 2nd & 3rd adiabatic invariants (in particular their proxies I & L*) are commonly used to characterize charged particle motion in a magnetic field. However care should be taken when calculating them, as the assumption of their conservation is not valid everywhere in the Earth’s magnetosphere. In this paper we compare calculations of I and L* using LANLstar, SPENVIS, IRBEM and a 3D particle tracer, and we map the areas in the Earth’s magnetosphere where I & L* can be assumed to be conserved.
J. A. Ruiz-Arias, J. Dudhia, and C. A. Gueymard
Geosci. Model Dev., 7, 1159–1174, https://doi.org/10.5194/gmd-7-1159-2014, https://doi.org/10.5194/gmd-7-1159-2014, 2014
D. T. Welling, J. Koller, and E. Camporeale
Geosci. Model Dev., 5, 277–287, https://doi.org/10.5194/gmd-5-277-2012, https://doi.org/10.5194/gmd-5-277-2012, 2012
J. Koller and S. Zaharia
Geosci. Model Dev., 4, 669–675, https://doi.org/10.5194/gmd-4-669-2011, https://doi.org/10.5194/gmd-4-669-2011, 2011
Cited articles
Battarbee, M., Hannuksela, O. A., Pfau-Kempf, Y., von Alfthan, S., Ganse, U., Jarvinen, R., Leo, Suni, J., Alho, M., lturc, Ilja, tvbrito, and Grandin, M.: fmihpc/analysator: v0.9, Zenodo [code], https://doi.org/10.5281/zenodo.4462515, 2021. a
Berger, M. J. and Jameson, A.: Automatic adaptive grid refinement for the Euler equations, AIAA J., 23, 561–568, https://doi.org/10.2514/3.8951, 1985. a
Childs, H., Brugger, E., Whitlock, B., Meredith, J., Ahern, S., Pugmire, D., Biagas, K., Miller, M. C., Harrison, C., Weber, G. H., Krishnan, H., Fogal, T., Sanderson, A., Garth, C., Bethel, E. W., Camp, D., Rubel, O., Durant, M., Favre, J. M., and Navratil, P.: High Performance Visualization–Enabling Extreme-Scale Scientific Insight, edited by: Bethel, E. W., Childs, H., and Hansen, C., 1st Edn., Chapman and Hall/CRC, 520 pp., https://doi.org/10.1201/b12985, 2012. a
Devine, K., Boman, E., Heapby, R., Hendrickson, B., and Vaughan, C.: Zoltan Data Management Service for Parallel Dynamic Applications, Comput. Sci. Eng., 4, 90–97, https://doi.org/10.1109/5992.988653, 2002. a, b, c
Dubart, M., Ganse, U., Osmane, A., Johlander, A., Battarbee, M., Grandin, M., Pfau-Kempf, Y., Turc, L., and Palmroth, M.: Resolution dependence of magnetosheath waves in global hybrid-Vlasov simulations, Ann. Geophys., 38, 1283–1298, https://doi.org/10.5194/angeo-38-1283-2020, 2020. a, b, c
Ganse, U., Koskela, T., Battarbee, M., Pfau-Kempf, Y., Papadakis, K., Alho, M., Bussov, M., Cozzani, G., Dubart, M., George, H., Gordeev, E., Grandin, M., Horaites, K., Suni, J., Tarvus, V., Kebede, F. T., Turc, L., Zhou, H., and Palmroth, M.: Enabling technology for global 3D + 3V hybrid-Vlasov simulations of near-Earth space, Phys. Plasmas, 30, 042902, https://doi.org/10.1063/5.0134387, 2023. a, b, c, d
Gombosi, T. I., Chen, Y., Glocer, A., Huang, Z., Jia, X., Liemohn, M. W., Manchester, W. B., Pulkkinen, T., Sachdeva, N., Shidi, Q. A., Sokolov, I. V., Szente, J., Tenishev, V., Toth, G., van der Holst, B., Welling, D. T., Zhao, L., and Zou, S.: What sustained multi-disciplinary research can achieve: The space weather modeling framework, J. Space Weather Space Clim., 11, 42, https://doi.org/10.1051/swsc/2021020, 2021. a
Graham, D. B., Khotyaintsev, Y. V., André, M., Vaivads, A., Chasapis, A., Matthaeus, W. H., Retinò, A., Valentini, F., and Gershman, D. J.: Non-Maxwellianity of Electron Distributions Near Earth's Magnetopause, J. Geophys. Res.-Space, 126, e29260, https://doi.org/10.1029/2021JA029260, 2021. a
Hoilijoki, S., Palmroth, M., Walsh, B. M., Pfau-Kempf, Y., von Alfthan, S., Ganse, U., Hannuksela, O., and Vainio, R.: Mirror modes in the Earth's magnetosheath: Results from a global hybrid-Vlasov simulation, J. Geophys. Res.-Space, 121, 4191–4204, https://doi.org/10.1002/2015JA022026, 2016. a
Honkonen, I.: fmihpc/dccrg: dccrg, Github [code], https://github.com/fmihpc/dccrg, 2023. a
Honkonen, I., von Alfthan, S., Sandroos, A., Janhunen, P., and Palmroth, M.: Parallel grid library for rapid and flexible simulation development, Comput. Phys. Commun., 184, 1297–1309, https://doi.org/10.1016/j.cpc.2012.12.017, 2013. a, b
Janhunen, P., Palmroth, M., Laitinen, T., Honkonen, I., Juusola, L., Facsko, G., and Pulkkinen, T. I.: The GUMICS-4 global MHD magnetosphere-ionosphere coupling simulation, J. Atmos. Sol.-Terr. Phys., 80, 48–59, https://doi.org/10.1016/j.jastp.2012.03.006, 2012. a, b, c
Kempf, Y., Pokhotelov, D., von Alfthan, S., Vaivads, A., Palmroth, M., and Koskinen, H. E. J.: Wave dispersion in the hybrid-Vlasov model: Verification of Vlasiator, Phys. Plasmas, 20, 112–114, https://doi.org/10.1063/1.4835315, 2013. a
Kotipalo, L.: AMR Test Configuration, University of Helsinki [data set], Matemaattis-luonnontieteellinen tiedekunta, https://doi.org/10.23729/f7f9d95d-1e23-49e3-9c35-1d8e26c47bf7, 2023. a, b
Müller, J., Simon, S., Motschmann, U., Schüle, J., Glassmeier, K.-H., and Pringle, G. J.: A.I.K.E.F.: Adaptive hybrid model for space plasma simulations, Comput. Phys. Commun., 182, 946–966, https://doi.org/10.1016/j.cpc.2010.12.033, 2011. a
Nishikawa, K., Duţan, I., Köhn, C., and Mizuno, Y.: PIC methods in astrophysics: simulations of relativistic jets and kinetic physics in astrophysical systems, Living Reviews in Computational Astrophysics, 7, 1, https://doi.org/10.1007/s41115-021-00012-0, 2021. a
Palmroth, M., Ganse, U., Pfau-Kempf, Y., Battarbee, M., Turc, L., Brito, T., Grandin, M., Hoilijoki, S., Sandroos, A., and von Alfthan, S.: Vlasov methods in space physics and astrophysics, Living Reviews in Computational Astrophysics, 4, 1, https://doi.org/10.1007/s41115-018-0003-2, 2018. a, b
Papadakis, K., Pfau-Kempf, Y., Ganse, U., Battarbee, M., Alho, M., Grandin, M., Dubart, M., Turc, L., Zhou, H., Horaites, K., Zaitsev, I., Cozzani, G., Bussov, M., Gordeev, E., Tesema, F., George, H., Suni, J., Tarvus, V., and Palmroth, M.: Spatial filtering in a 6D hybrid-Vlasov scheme to alleviate adaptive mesh refinement artifacts: a case study with Vlasiator (versions 5.0, 5.1, and 5.2.1), Geosci. Model Dev., 15, 7903–7912, https://doi.org/10.5194/gmd-15-7903-2022, 2022. a
Pfau-Kempf, Y., Battarbee, M., Ganse, U., Hoilijoki, S., Turc, L., von Alfthan, S., Vainio, R., and Palmroth, M.: On the Importance of Spatial and Velocity Resolution in the Hybrid-Vlasov Modeling of Collisionless Shocks, Front. Phys., 6, 44, https://doi.org/10.3389/fphy.2018.00044, 2018. a, b, c
Pfau-Kempf, Y., von Alfthan, S., Ganse, U., Battarbee, M., Kotipalo, L., Koskela, T., Honkonen, I., Sandroos, A., Papadakis, K., Alho, M., Zhou, H., Palmu, M., Grandin, M., Suni, J., Pokhotelov, D., and Horaites, K.: fmihpc/vlasiator: Vlasiator 5.3, Zenodo [code], https://doi.org/10.5281/zenodo.10600112, 2024. a, b
Rembiasz, T., Obergaulinger, M., Cerdá-Durán, P., Ángel Aloy, M., and Müller, E.: On the Measurements of Numerical Viscosity and Resistivity in Eulerian MHD Codes, Astrophys. J., 230, 18, https://doi.org/10.3847/1538-4365/aa6254, 2017. a
Stout, Q. F., De Zeeuw, D. L., Gombosi, T. I., Groth, C. P. T., Marshall, H. G., and Powell, K. G.: Adaptive Blocks: A High Performance Data Structure, in: Proceedings of the 1997 ACM/IEEE Conference on Supercomputing, SC '97, 1–10, Association for Computing Machinery, New York, NY, USA, ISBN 0897919858, https://doi.org/10.1145/509593.509650, 1997. a
Swisdak, M.: Quantifying gyrotropy in magnetic reconnection, Geophys. Res. Lett., 43, 43–49, https://doi.org/10.1002/2015GL066980, 2016. a
von Alfthan, S., Pokhotelov, D., Kempf, Y., Hoilijoki, S., Honkonen, I., Sandroos, A., and Palmroth, M.: Vlasiator: First global hybrid-Vlasov simulations of Earth's foreshock and magnetosheath, J. Atmos. Sol.-Terr. Phys., 120, 24–35, https://doi.org/10.1016/j.jastp.2014.08.012, 2014. a, b, c
Wang, X., Chen, Y., and Tóth, G.: Global Magnetohydrodynamic Magnetosphere Simulation With an Adaptively Embedded Particle-In-Cell Model, J. Geophys. Res.-Space, 127, e2021JA030091, https://doi.org/10.1029/2021JA030091, 2022. a
Zerroukat, M. and Allen, T.: A three-dimensional monotone and conservative semi-Lagrangian scheme (SLICE-3D) for transport problems, Q. J. Roy. Meteor. Soc., 138, 1640–1651, https://doi.org/10.1002/qj.1902, 2012. a
Zhang, W., Myers, A., Gott, K., Almgren, A., and Bell, J.: AMReX: Block-structured adaptive mesh refinement for multiphysics applications, Int. J. High Perform. Comput. Appl., 35, 508–526, https://doi.org/10.1177/10943420211022811, 2021. a
Short summary
This paper examines a method called adaptive mesh refinement in optimization of the space plasma simulation model Vlasiator. The method locally adjusts resolution in regions which are most relevant to modelling, based on the properties of the plasma. The runs testing this method show that adaptive refinement manages to highlight the desired regions with manageable performance overhead. Performance in larger-scale production runs and mitigation of overhead are avenues of further research.
This paper examines a method called adaptive mesh refinement in optimization of the space plasma...
