ST – Solar-Terrestrial Sciences

EGU22-901 | Presentations | ST2.4 | Arne Richter Award for Outstanding ECS Lecture

Ion dynamics in the inner magnetosphere during Van Allen Probe Era 

Chao Yue

Ion dynamics are controlled by the energy-dependent source, transport, energization, and loss processes. Systematic changes in the ion dynamics are essential to understand the ring current variations in the inner magnetosphere. The Van Allen Probes mission, which orbits near the equatorial plane inside the geosynchronous orbit, has a wide energy coverage with high energy resolution and state-of-the-art ion composition instrumentation. It provides a great opportunity to investigate plasma dynamics. In this talk, I will present some of our recent studies on the ion dynamics of different populations and species as well as the related plasma wave activity during geomagnetic quiet and active times.

How to cite: Yue, C.: Ion dynamics in the inner magnetosphere during Van Allen Probe Era, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-901,, 2022.

EGU22-1115 | Presentations | MAL13 | Hannes Alfvén Medal Lecture

Our Heliosphere and its Interstellar Interaction: A Solar Cycle of Global Observations and Discoveries 

David McComas

The supersonic solar wind and its embedded magnetic field continuously flow outward in all directions from the sun. This magnetized plasma inflates a bubble – the heliosphere – in the very-local interstellar medium (VLISM). The Interstellar Boundary Explorer (IBEX) mission launched in late 2008 and has been continuously returning 3-D global images of Energetic Neutral Atoms (ENAs) that derive from ion populations in the heliosheath and beyond. Now spanning more than a full solar cycle, IBEX’s all-sky maps and observations uniquely inform the global outer heliosphere and its evolving interstellar interaction. Insights from IBEX, in concert with in situ observations by the two Voyager spacecraft, which were transiting two different trajectories through the outer boundaries of the heliosphere contemporaneously with IBEX, have led to a true scientific revolution in our understanding of the outer heliosphere and its interstellar interaction. This Hannes Alfvén Medal Lecture will summarize some of the many discoveries and “firsts” from the IBEX mission and their implications for the outer heliosphere and VLISM. Finally, we will also look forward to the promise of the even more advanced Interstellar Mapping and Acceleration Probe (IMAP) mission, which is under development and slated to launch in 2025.

How to cite: McComas, D.: Our Heliosphere and its Interstellar Interaction: A Solar Cycle of Global Observations and Discoveries, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-1115,, 2022.

EGU22-8388 | Presentations | MAL13 | ST Division Outstanding ECS Award Lecture

A Revised Collection of Sunspot Group Numbers: Context and Future Improvements 

Víctor Carrasco

Rudolf Wolf, first director of the Zürich Observatory around mid-19th century, recovered a large number of sunspot observations made by astronomers several solar cycles back. Based on that database, he defined the relative sunspot number from the number of sunspot groups and individual sunspots. He extended his daily and monthly series until 1749, whereas his yearly series to 1700 (Clette et al. 2014). Nowadays, the World Data Center Sunspot Index and Long-term Solar Observations is the responsible to maintain this sunspot number index. At the end of the 20th century, Hoyt and Schatten (1998) compiled more sunspot observations made by astronomers since the beginning of the 17th century. Thus, they created the group sunspot number index from the number of sunspot groups. Unlike the relative sunspot number, their series starts in 1610.

More recently, several works have detected some problems both in these two indices and the databases. For example, Vaquero et al. (2016) published a revised collection of sunspot group numbers correcting some of the mistakes found in the Hoyt and Schatten database, in addition to incorporate other unknown sunspot records. Currently, there is an ongoing global effort to improve the weakness of the database and recalibrate the indices. Some remarkable improvements to be carried out in future versions of the sunspot number databases have been made regarding the earliest sunspot observations recorded by astronomers such as Galileo and Scheiner, inter alia. Then, corrections of significant mistakes detected in the sunspot counting assigned to these observers in the existing databases are proposed as well as the incorporation of telescopic sunspot records made by the earliest observers not included in these databases.

The sunspot number series is the index including the longest direct solar observation set to study the long-term solar activity evolution and its influence on the Earth. Therefore, we need that the databases, in which these indices are based, are free of problematic observations and, moreover, to improve their observational coverage before mid-19th century. Thus, we will understand better past, present and future solar activity.


Clette, F., Svalgaard, L., Vaquero, J.M., Cliver, E.W.: 2014, Revisiting the Sunspot Number. A 400-Year Perspective on the Solar Cycle, SSRv 186, 35. DOI: 10.1007/s11214-014-0074-2.

Hoyt, D.V., Schatten, K.H.: 1998, Group sunspot numbers: a new solar activity reconstruction. Solar Phys. 179, 189. DOI: 10.1023/A:1005007527816.

Vaquero, J.M., Svalgaard, L., Carrasco, V.M.S., Clette, F., Lefèvre, L., Gallego, M.C., Arlt, R., Aparicio, A.J.P., Richard, J.-G., Howe, R.: 2016, A Revised Collection of Sunspot Group Numbers, Sol. Phys. 291, 3061. DOI: 10.1007/s11207-016-0982-2.

How to cite: Carrasco, V.: A Revised Collection of Sunspot Group Numbers: Context and Future Improvements, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-8388,, 2022.

ST1 – The Sun and Heliosphere

EGU22-307 | Presentations | ST1.1

Joint Instabilities of Sheared Flows and Magnetic Fields 

Patrick Lewis, David Hughes, and Evy Kersale

Shear flows and magnetic fields are ubiquitous in astrophysical bodies such as stars and accretion discs. Furthermore,
the interaction between flows and magnetic field plays a key role in the dynamics of plasma fusion devices. Typically,
the flows and magnetic field are both sheared, and it is therefore a problem of fundamental importance to understand
the instabilities that may occur in such a system.

In the absence of magnetic field, the linear stability of a viscous sheared flow is governed by the Orr-Sommerfeld
equation; this is one of the classic problems of hydrodynamics. At the other limit, there are somewhat analogous
instabilities of a fluid of finite electrical conductivity containing a static sheared magnetic field. These are related to
the classical tearing modes that have received considerable attention in both the astrophysical and plasma physics

In general though, the fluid flow and the magnetic field will both be important players. Previous studies have investigated
configurations which have served as models for systems such as the magnetotail and solar surges. While these
investigations have been fruitful, the prescription of the basic field and flow, while physically motivated, have been
chosen somewhat arbitrarily. It is therefore of interest to consider the instability problem within this more general

Motivated astrophysically, such as by the dynamics in the solar tachocline, here we consider a self-consistent problem
in which both instabilities can occur. In particular, we consider the stability of equilibrium states arising from the
shearing of a uniform magnetic field by a forced transverse flow. The problem is governed by three non-dimensional
parameters: the Chandrasekhar number, and the flow and magnetic Reynolds numbers. In opposite limits of parameter
space, we recover the predictions of the aforementioned classical problems. As we move through this three-dimensional
parameter space, a range of interactions are possible: We demonstrate the stabilisation of a purely hydrodynamic
instability through the magnetic field, show the existence of a joint instability outlining the physical mechanisms at
play, and demonstrate that under certain conditions, hydrodynamically-stable parallel shear flows lead to instability
growth rates that exceed those of static tearing modes. To conclude, we elucidate the consequences of considering
the linear stability of an evolving background state and show that a quasi-static approach may not be meaningful. In
these circumstances, it therefore becomes essential to perform a stability analysis of a time-varying basic state.



How to cite: Lewis, P., Hughes, D., and Kersale, E.: Joint Instabilities of Sheared Flows and Magnetic Fields, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-307,, 2022.

EGU22-1496 | Presentations | ST1.1

Modelling magneto-hydro-static equilibria in quiet Sun regions 

Thomas Wiegelmann and Maria Madjarska

The solar magnetic field is measured routinely only
in the solar photosphere. While reliable measurements
of the photospheric magnetic field vector are only available
in active regions, in the quiet Sun at present only
 the vertical component of the magnetic field can be obtained accurately.
 To derive magnetic field structures throughout the solar atmosphere, from
 the chromosphere to the corona, we extrapolate these photospheric measurements
into the upper photosphere, chromosphere and corona with a magneto-hydro-static
model. We optimize free model parameters by comparing the modelled magnetic
field lines with structures observed in solar images. The comparison is
done automatically with a number of quantitative measurements and
the optimal model parameters are found with the help of
a downhill simplex minimization. This newly developed modelling approach can
provide an accurate and deep understanding of the magnetic field structures
that extend to any height in the solar atmosphere.

How to cite: Wiegelmann, T. and Madjarska, M.: Modelling magneto-hydro-static equilibria in quiet Sun regions, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-1496,, 2022.

EGU22-1598 | Presentations | ST1.1

Oblique instabilities driven by pickup ion ring-beam distributions in the outer heliosheath 

Ameneh Mousavi, Kaijun Liu, and Sina Sadeghzadeh

The energetic neutral atom (ENA) ribbon observed by the Interstellar Boundary Explorer (IBEX) spacecraft is believed to originate from the pickup ions in the outer heliosheath. The outer heliosheath pickup ions generally have a ring-beam velocity distribution at a certain pickup angle, α, the angle at which these ions are picked up by the interstellar magnetic field. The pickup ion ring-beam distributions can drive unstable waves of different propagation angles with respect to the background interstellar magnetic field, θ. Previous studies of the outer heliosheath pickup ion dynamics were mainly focused on ring-like pickup ion distributions with α≈90° and/or the parallel- and anti-parallel-propagating unstable waves (θ=0°and 180°). The present study carries out linear kinetic instability analysis to investigate both the parallel and oblique unstable modes (0°≤θ≤180°) driven by ring-beam pickup ion distributions of different pickup angles between 0° and 90°. Our linear analysis reveals that ring-beam pickup ions can excite mirror waves as well as oblique left-helicity waves and their harmonics. The maximum growth rate of the mirror mode increases with increasing α. On the other hand, the wavenumber and growth rate of the most unstable oblique left-helicity modes are consistent with the unstable modes of 0°and 180° examined in our earlier work.

How to cite: Mousavi, A., Liu, K., and Sadeghzadeh, S.: Oblique instabilities driven by pickup ion ring-beam distributions in the outer heliosheath, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-1598,, 2022.

EGU22-1705 | Presentations | ST1.1

Characteristic ion length scales for four types of interplanetary shocks 

Byeongseon Park, Alexander Pitna, Jana Safrankova, and Zdenek Nemecek

The interaction between interplanetary (IP) shock and solar wind has been studied for the understanding of the energy dissipation mechanism within the collisionless plasma. The power spectra of the magnetic field exhibit breaks, where steepening of these spectra occurs. These breaks have been observed and also regarded as a threshold distinguishing the kinetic range from the inertial range of turbulence. Different heating processes can be related mainly to two characteristic ion length scales — ion inertial length and ion gyroradius. We attempt to establish the relation between these length scales and the spectral break. Data for four different types of IP shocks (fast forwards, fast reverse, slow forwards, slow reverse) measured for 2 hours (one hour for up and downstream plasma) by WIND at 1 AU were used. Continuous wavelet transform for the estimation of the power spectra of measured magnetic field was employed. Spectral breaks were determined by fitting 2-segment piecewise linear function around the expected break position in log-log space. Preliminary analysis of these spectral breaks and the characteristic length scales in fast shocks yields results consistent with the previous studies. Additionally, we extended this analysis towards slow shocks and obtained similar results. While the level of power enhancement of the magnetic field due to fast shocks reaches the order of ten on average, only the order of one was shown for slow shocks. The level of the compression of the characteristic spatial scales, however, is approximately similar for fast and slow shocks.

How to cite: Park, B., Pitna, A., Safrankova, J., and Nemecek, Z.: Characteristic ion length scales for four types of interplanetary shocks, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-1705,, 2022.

EGU22-1895 | Presentations | ST1.1

Towards a realistic evaluation of transport coefficients in non-equilibrium plasmas from space 

Edin Husidic, Marian Lazar, Klaus Scherer, Horst Fichtner, and Stefaan Poedts

While space plasmas are largely considered to be nearly collisionless, at relatively low heliocentric distances, that is, below 1 AU, particle-particle collisions still play an important role in the transport of matter, momentum, and energy. A way to quantify these processes macroscopically, e.g., in fluid models, is within classical transport theory, where fluxes and their sources are linearly related by transport coefficients. In the solar wind context, of particular interest are the observed velocity distributions of plasma particles with Kappa-distributed suprathermal tails, conditioned not only by binary collisisions, but also by their interaction with plasma waves and turbulence. We present first derivations of the main transport coefficients based on regularised Kappa distributions (RKDs), which, unlike standard Kappa distributions (SKDs), enable a macroscopic description of non-equilibrium plasmas without mathematical divergences or physical inconsistencies. All transport coefficients are finite, well defined for all values of κ > 0, and markedly enhanced in the presence of suprathermal electrons. The results indicate that for low values of κ, that is, below the SKD poles, the transport coefficients can be many orders of magnitudes higher than the corresponding Maxwellian limits, which can lead to significant underestimations if suprathermal electrons are ignored. Moreover, we show the importance of an adequate Kappa modeling of suprathermal populations by contrasting our results to other modified interpretations that underestimate the effects of suprathermals.

How to cite: Husidic, E., Lazar, M., Scherer, K., Fichtner, H., and Poedts, S.: Towards a realistic evaluation of transport coefficients in non-equilibrium plasmas from space, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-1895,, 2022.

EGU22-2814 | Presentations | ST1.1

Hybrid Particle-In-Cell model with Adaptive Mesh Refinement 

Nicolas Aunai, Roch Smets, Philip Deegan, Andrea Ciardi, and Alexis Jeandet

Collisionless magentized plasmas a priori need to be evolved using Vlasov-Maxwell kinetic formalism.
However the tremendous number of spatial and temporal scales involved in phenomena of interest makes it prohibitive, from a computational standpoint.
Fully kinetic particle in cell and single fluid MHD codes are commonly used at very small or very large scales.
The hybrid formalism, treating ions kinetically and electrons as a fluid, is in principle advantageous to fill the gap between these two extremities.
However, a correct treatment of critical regions such as reconnection X-lines require a good resolution of sub-ion dissipative scales, which
still constitute a major challenge if aiming at simulating meso/macro scale systems.
This work presents a new code, named PHARE, which successfully implements the adaptive mesh refinement mechanism in a hybrid particle-in-cell code.
Such a code is able to dynamically focus the resolution in critical regions while others not only have a coarser spatial resolution, but are also
evolved much less often thanks to a recursive time stepping procedure.
Adopting a patch based AMR mechanism, the code architecture is made so that the specific solver/physical model that is solved at a given refined level
is abstracted, thus giving the opportunity to handling multi-formalisms AMR patch hierarchies, where, for instance, coarsest levels are solved in MHD while
dynamicall created refined levels are solved within the Hybrid framework.

How to cite: Aunai, N., Smets, R., Deegan, P., Ciardi, A., and Jeandet, A.: Hybrid Particle-In-Cell model with Adaptive Mesh Refinement, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-2814,, 2022.

EGU22-2862 | Presentations | ST1.1 | Highlight

The ESA Virtual Space Weather Modelling Centre 

Stefaan Poedts and the VSWMC-P3 team

The ESA Virtual Space Weather Modelling Centre (VSWMC) project was defined as a long term project including different successive parts. Parts 1 and 2 were completed in the first 4-5 years and designed and developed a system that enables models and other components to be installed locally or geographically distributed and to be coupled and run remotely from the central system. A first, limited version went operational in May 2019 under the H-ESC umbrella on the ESA SSA SWE Portal. It is similar to CCMC but interactive (no runs on demand) and the models are geographically distributed and coupled over the internet.

The goal of the ESA project "Virtual Space Weather Modelling Centre - Part 3" (2019-2021) was to further develop the Virtual Space Weather Modelling Centre, building on the Part 2 prototype system and focusing on the interaction with the ESA SSA SWE system. The objectives and scope of this new project include maintaining the current operational system, the efficient integration of 11 new models and many new model couplings, including daily automated end-to-end (Sun to Earth) simulations, the further development and wider use of the coupling toolkit  and front-end GUI, making the operational system more robust and user-friendly. The VSWMC-Part 3 project finished recently.

The 11 new models that have been integrated are Wind-Predict (a global coronal model from CEA, France), the Coupled Thermosphere/Ionosphere Plasmasphere (CTIP) model, Multi-VP (another global coronal model form IRAP/CNRS, France), the BIRA Plasma sphere Model of electron density and temperatures inside and outside the plasmasphere coupled with the ionosphere (BPIM, Belgium), the SNRB  (also named SNB3GEO) model for electron fluxes at geostationary orbit (covering the GOES 15 energy channels >800keV and >2MeV) and the SNGI geomagnetic indices Kp and Dst models (University of Sheffield, UK), the SPARX Solar Energetic Particles transport model (University of Central Lancashire, UK), Spenvis DICTAT tool for s/c internal charging analysis (BISA, Belgium), the Gorgon magnetosphere model (ICL, UK), and the Drag Temperature Model (DTM) and operations-focused whole atmosphere model MCM being developed in the H2020 project SWAMI. Many new couplings have also been implemented and a dynamic coupling facility has been installed. Moreover, Daily runs are implemented of two model chains covering the whole Sun-to-Earth domain. The results of these daily runs are made available to all VSWMC users.

We will provide an overview of the state-of-the-art, including the new available model couplings and daily model chain runs, and demonstrate the system.

How to cite: Poedts, S. and the VSWMC-P3 team: The ESA Virtual Space Weather Modelling Centre, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-2862,, 2022.

The acceleration of high-energy charged particles is common both in heliophysics and astrophysics. Although the diffusive shock acceleration (DSA) has been the well-accepted standard mechanism for particle acceleration at shocks, the fundamental issue is that DSA does not predict the number of accelerated particles. In other words, it relies on an unprescribed injection mechanism that provides a seed population from which the particle acceleration proceeds. Resolving the so-called injection problem is more challenging for electrons than ions because scattering low-energy electrons requires high-frequency waves, which are usually much lower in intensity than low-frequency fluctuations.


We have proposed stochastic shock drift acceleration (SSDA) as a plausible electron injection mechanism that can take place within the transition layer of quasi-perpendicular shocks [Katou & Amano, 2019]. The energy gain mechanism of SSDA is essentially the same as the conventional shock drift acceleration (SDA), but the presence of stochastic pitch-angle scattering makes the acceleration more efficient. We will show that the theoretical predictions nicely explain in-situ observations by Magnetospheric MultiScale (MMS) spacecraft [Amano et al. 2020]. Recent fully kinetic Particle-In-Cell (PIC) simulation results will also be shown, in which we found signatures of SSDA [Matsumoto et al. 2017, Kobzar et al. 2021]. We will also present an extended theoretical model that unifies SSDA and DSA. The model predicts a wide range of the energy spectrum from sub-relativistic and relativistic energies. The particle acceleration in the sub-relativistic energy will be dominated by SSDA, which has a spectral index steeper than the standard DSA. Under certain conditions, the particle acceleration mechanism may smoothly transition from SSDA to DSA, and the spectral index approaches the canonical DSA prediction. Therefore, the model can consistently describe the whole particle acceleration process, including the injection by SSDA and the main acceleration to cosmic-ray energies by DSA. We argue that the electron injection scheme through SSDA is realized preferentially at shocks with higher Alfven Mach numbers defined in the Hoffmann-Teller frame.

How to cite: Amano, T.: Electron Injection via Stochastic Shock Drift Acceleration: Theory, Simulation, and Observation, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-3277,, 2022.

EGU22-3304 | Presentations | ST1.1 | Highlight

Recent progress on developments of the magnetohydrostatic extrapolation 

Xiaoshuai Zhu, Thomas Wiegelmann, and Bernd Inhester

The magnetohydrostatic (MHS) extrapolation is developed to model the three-dimensional magnetic fields and plasma of the solar atmosphere with the measured vector magnetogram used as boundary condition. It differs from the nonlinear force-free field (NLFFF) extrapolation in that it takes into account plasma forces include pressure gradient and gravity. In this presentation, I will report the recent progress in two aspects on developments of the MHS extrapolation. The first one is the development of a preprocessing method to deal with the non-MHS vector magetograms. The reason of doing this is that there are a small number of the vector magnetograms which are not consistent with the MHS equilibria. The second aspect is the combination of the MHS extrapolation and the NLFFF extrapolation to improve the efficiency of the computation.

How to cite: Zhu, X., Wiegelmann, T., and Inhester, B.: Recent progress on developments of the magnetohydrostatic extrapolation, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-3304,, 2022.

EGU22-3345 | Presentations | ST1.1

Fluid modeling of collisionless plasmas with a non-local kinetic closure 

Taiki Jikei and Takanobu Amano

Modeling of collisionless plasmas can be divided into two categories: fluid models and kinetic models. Generally speaking, fluid models require less computational resources than kinetic models, so they are suited for large-scale simulations. However, conventional fluid models such as MHD ignores wave-particle interaction (WPI). It has been pointed out that WPI affects microscopic and macroscopic dynamics and should not be ignored even in MHD scales. This creates a demand for a fluid model of collisionless plasma that takes into account WPI effects.

We have developed a fluid model called the cyclotron resonance closure (CRC) model [1]. The CRC model reproduces linear cyclotron resonance effects using a non-local closure method similar to the Landau closure model. The CRC model reproduces the linear cyclotron resonance and linear growth of ion temperature anisotropy instabilities qualitatively correct. We have also shown that the quasilinear relaxation of temperature anisotropy via resonant waves incorporated in the CRC model.

 Another example of a kinetic fluid model is the well-known Chew-Goldberger-Low (CGL) model. The CGL model is used to analyze low-frequency waves in collisionless plasmas. The CGL model enriched by finite Larmor radius correction and Landau closure predicts the growth rate of firehose instability with reasonable accuracy. However, the CGL model cannot reproduce cyclotron resonance effects such as cyclotron damping and electromagnetic ion cyclotron (EMIC) instability because of the low-frequency assumption.

We will discuss some basic concepts of these kinetic fluid models and their range of application, especially in nonlinear simulation. The CRC model is not limited by frequency (at least up to cyclotron frequency) and can be used for both EMIC and parallel firehose instabilities but need improvement for quantitative agreement with fully kinetic models.

The CGL model can be very accurate in linear analysis of low-frequency waves. We compared the CRC and the CGL model using a simulation of an initially parallel firehose unstable system [2]. We found that the low-frequency approximation of the CGL model fails in some parameters after the appearance of high-frequency oscillation in the nonlinear stage. Also, the absence of cyclotron damping in the CGL model results in a quasi-steady final state that is not consistent with marginal stability analysis.

We conclude that a kinetic fluid model that does not make the low-frequency approximation should be considered instead of the CGL-based approach. The CRC model is a candidate for such a model that can be used in a wide range of parameters.



1. T. Jikei and T. Amano (2021), A non-local fluid closure for modeling cyclotron resonance in collisionless magnetized plasmas, Physics of Plasmas, 28, 042105

2. T. Jikei and T. Amano (2022), Critical comparison of collisionless fluid models: Nonlinear simulations of parallel firehose instability, Physics of Plasmas, Accepted

How to cite: Jikei, T. and Amano, T.: Fluid modeling of collisionless plasmas with a non-local kinetic closure, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-3345,, 2022.

EGU22-3911 | Presentations | ST1.1

Flux tube dependent propagation of Alfvén waves in the solar corona 

Chaitanya Prasad Sishtla, Jens Pomoell, Emilia Kilpua, Simon Good, Farhad Daei, and Minna Palmroth
Alfvén wave turbulence has emerged as an important heating mechanism to accelerate the solar wind. The generation of this turbulent heating is dependent on the presence and subsequent interaction of counter-propagating alfvén waves. This requires us to understand the propagation and evolution of alfvén waves in the solar wind in order to develop an understanding of the relationship between turbulent heating and solar wind parameters. In this paper we aim to study the  response  of  the  solar  wind  upon  injecting  monochromatic  single frequency alfvén waves at the base of the corona for various magnetic flux tube geometries. We use an ideal magnetohydrodynamic (MHD) model using an adiabatic equation of state. An alfvén pump wave is injected into the quiet solar wind by perturbing the transverse magnetic field and velocity components. The alfvén waves were found to be reflected due to the development of the parametric decay instability (PDI). Further investigation revealed that the PDI was suppressed both by efficient reflections at low frequencies as well as magnetic flux tube geometries.

How to cite: Sishtla, C. P., Pomoell, J., Kilpua, E., Good, S., Daei, F., and Palmroth, M.: Flux tube dependent propagation of Alfvén waves in the solar corona, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-3911,, 2022.

EGU22-4064 | Presentations | ST1.1

Mating of two current circuits – A new type of plasma instability 

Gerhard Haerendel

In the context of the onset of substorms, observations had revealed the frequent occurrence of an unidentified trigger process. Inspections of such onsets had led to the discovery that high-beta plasma at the magnetosphere-tail boundary became suddenly unstable when so-called auroral streamers lined up closely with that boundary. The manifestation was the appearance of bead-like auroral structures preceding the auroral breakup and subject to nonlinear growth. Highly resolved video coverage of a few events showed that the fast moving beads moved opposite to the convective motions on either side. This led to the proposal [Haerendel & Frey 2021] that the trigger of the instability was the formation of a new current circuit, by a non-MHD process, in the gap between the two adjacent circuits of the high-beta plasma boundary and the auroral streamer. Observational evidence and a model of the current structure will be presented.

How to cite: Haerendel, G.: Mating of two current circuits – A new type of plasma instability, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-4064,, 2022.

EGU22-4892 | Presentations | ST1.1 | Highlight

The 3D Dynamics of Solar Flare Magnetic Reconnection 

Spiro Antiochos, Joel Dahlin, and Richard DeVore

Solar eruptive events (SEE) consisting of a massive filament eruption, intense X-class flare, and a fast CME are the most powerful manifestations of explosive energy release in our solar system and the primary drivers of highly destructive space weather at Earth and in interplanetary space. These giant events, which have global scale of solar radii, allow us to study in great detail fundamental space plasma processes such as magnetic reconnection that are important to many cosmic phenomena. Both 2.5D and 3D numerical simulations have shown that the fast energy release is due to reconnection in a large-scale current sheet that forms in the corona, but the 3D dynamics of the reconnection are far from understood.  The greatest challenge to understanding SEEs is their extreme rate of energy release, and for some events, the amazing efficiency at converting magnetic energy into high-energy particle energy. We present new ultra-high-resolution 3D simulations of flare reconnection using the adaptive-mesh-refinement code ARMS. We find that the reconnection dynamics are dominated by 3D magnetic islands, and show that the islands should have clear observational signatures, especially in the so-called flare ribbons that are commonly observed in the chromosphere. We discuss the central role of the islands for understanding the multiscale coupling at the heart of reconnection, the fast energy release rate, and the high efficiency of particle acceleration.


This work was supported by the NASA Living With a Star Program and by the NASA DRIVE Center Program.


How to cite: Antiochos, S., Dahlin, J., and DeVore, R.: The 3D Dynamics of Solar Flare Magnetic Reconnection, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-4892,, 2022.

EGU22-5498 | Presentations | ST1.1

Simulating the FIP effect in coronal loops using a multi-species kinetic-fluid model. 

Nicolas Poirier, Michael Lavarra, Alexis Rouillard, Pierre-Louis Blelly, Victor Réville, Andrea Verdini, Marco Velli, Eric Buchlin, and Mikel Indurain

We investigate abundance variations of heavy ions in coronal loops. We develop and exploit a multi-species model of the solar atmosphere (called IRAP’s Solar Atmospheric Model: ISAM) that solves for the transport of neutral and charged particles from the chromosphere to the corona. We investigate the effect of different mechanisms that could produce the First Ionization Potential (FIP) effect. We compare the effects of the thermal, friction and ponderomotive force. The propagation, reflection and dissipation of Alfvén waves is solved using two distinct models, the first one from Chandran et al. (2011) and the second one that is a more sophisticated turbulence model called Shell-ATM. ISAM solves a set of 16-moment transport equations for both neutrals and charged particles. Protons and electrons are heated by Alfvén waves, which then heat up the heavy ions via collision processes. We show comparisons of our results with other models and observations, with an emphasis on FIP biases. This work was funded by the European Research Council through the project SLOW_SOURCE - DLV-819189.

How to cite: Poirier, N., Lavarra, M., Rouillard, A., Blelly, P.-L., Réville, V., Verdini, A., Velli, M., Buchlin, E., and Indurain, M.: Simulating the FIP effect in coronal loops using a multi-species kinetic-fluid model., EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-5498,, 2022.

EGU22-5546 | Presentations | ST1.1 | Highlight

A simulation of flare-driven coronal rain 

Wenzhi Ruan, Yuhao Zhou, and Rony Keppens

Coronal rains are cool materials (~10,000 K) that appear at hot corona. They are frequently observed in non-flaring loops of active regions and recently observed in flaring loops at gradual phases. Hot coronal loops (~10 MK) are often produced in flare events due to magnetic reconnection. The hot flare loops gradually recover to typical coronal temperature due to thermal conduction and radiative loss, during which condensation can happen due to thermal instability. Here we demonstrate how the rains formed in a flare loop with a two-and-a-half dimensional magnetohydrodynamic simulation. We simulate a flare event from pre-flare phase all the way to gradual phase and successfully reproduce coronal rains. We find that thermal conduction and radiative losses alternately dominate the cooling of the flare loop. We find that runaway cooling and rain formation also induce the appearance of dark post-flare loop systems, as observed in extreme ultraviolet (EUV) channels.

How to cite: Ruan, W., Zhou, Y., and Keppens, R.: A simulation of flare-driven coronal rain, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-5546,, 2022.

EGU22-5633 | Presentations | ST1.1

Electron cyclotron maser model of solar radio zebras 

Jan Benáček and Marian Karlický
Solar radio zebras, occurring as fine structures in radiograms during Type IV bursts, are an excellent tool for plasma diagnostics during solar flares. The main model of zebras is the electron cyclotron maser instability based on double plasma resonance between electron cyclotron and plasma frequencies and an unstable wave in the presence of an unstable type of velocity distribution is necessary, e.g., a loss-cone. The radio emission occurs in the electromagnetic Z-mode along the magnetic field or at the first harmonic of the X-mode in the perpendicular direction. However, it is still unclear where and how the instability evolves and how the locally captured electrostatic waves are converted to escaping radio waves. To obtain the instability evolution, we calculated its growth rates and saturation energies as functions of cyclotron-to-plasma frequency ratio, loss-cone density, cold background  temperature, hot electron thermal velocity, and loss-cone angle by using analytical calculations and particle-in-cell simulations. We found that the growth rates and saturation energies form maxima, approximately located at the harmonic numbers of cyclotron frequency. The maxima shift to lower frequencies with increasing the plasma temperature, they broaden and decrease with increasing the harmonic number. We also estimated electromagnetic energy densities in the emission region and the conversion efficiency to the radio waves.

How to cite: Benáček, J. and Karlický, M.: Electron cyclotron maser model of solar radio zebras, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-5633,, 2022.

Ubiquitous small-scale vortical motions in the solar atmosphere are thought to play an important role in local heating of the quiet chromosphere and corona. However, an unambiguous method for their detection in observations and numerical simulations has not been found yet.
We aim at developing a robust method for the automated identification of vortices. Local and global rotations in the flow should be considered as both are necessary for the detection of coherent vortical structures. Moreover, the use of a threshold should be avoided to not exclude slow vortices in the identification process.
We present a new method that combines the rigor of mathematical criteria and the global perspective of morphological techniques. The core of the method is the estimation of the center of rotation for every point that presents some degree of local rotation in the flow. For that, we employ the Rortex criterion and the morphology of the neighboring velocity field. We then identify coherent vortical structures by clustering the estimated centers of rotation.
The application of the method to synthetic velocity fields demonstrates its reliability and accuracy. A first statistical study is performed on realistic numerical simulations of the solar atmosphere carried out with the radiative magneto-hydrodynamical code CO5BOLD. We counted on average 0.8 Mm-2 swirls in the photosphere and 1.9 Mm-2 at the bottom of the chromosphere. The average radius varies between 59 km and 72 km. Compared to previous studies, our analysis reveals more and smaller vortical motions in the simulated solar atmosphere. Moreover, we find that 84 % of the swirls in the photosphere show twists in the magnetic field lines compatible with torsional Alfvén waves. 

How to cite: Canivete Cuissa, J. R. and Steiner, O.: An innovative and automated vortex identification method based on the estimation of the center of rotation with application to solar simulations, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-5744,, 2022.

EGU22-6060 | Presentations | ST1.1

Coronal rain in randomly heated arcades 

Xiaohong Li, Rony Keppens, and Yuhao Zhou

Adopting the MPI-AMRVAC code, we present a 2.5-dimensional magnetohydrodynamic simulation, which includes thermal conduction and radiative cooling, to investigate the formation and evolution of the coronal rain phenomenon. We perform the simulation in initially linear force-free magnetic fields that host chromospheric, transition-region, and coronal plasma, with turbulent heating localized on their footpoints. Due to thermal instability, condensations start to occur at the loop top, and rebound shocks are generated by the siphon inflows. Condensations fragment into smaller blobs moving downwards, and as they hit the lower atmosphere, concurrent upflows are triggered. Larger clumps show us clear coronal rain showers as dark structures in synthetic EUV hot channels and as bright blobs with cool cores in the 304 Å channel, well resembling real observations. Following coronal rain dynamics for more than 10 hr, we carry out a statistical study of all coronal rain blobs to quantify their widths, lengths, areas, velocity distributions, and other properties. The coronal rain shows us continuous heating–condensation cycles, as well as cycles in EUV emissions. Compared to the previous studies adopting steady heating, the rain happens faster and in more erratic cycles. Although most blobs are falling downward, upward-moving blobs exist at basically every moment. We also track the movement of individual blobs to study their dynamics and the forces driving their movements. The blobs have a prominence-corona transition-region-like structure surrounding them, and their movements are dominated by the pressure evolution in the very dynamic loop system.

How to cite: Li, X., Keppens, R., and Zhou, Y.: Coronal rain in randomly heated arcades, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-6060,, 2022.

EGU22-6120 | Presentations | ST1.1 | Highlight

Modeling Reconnection and Turbulence in the Magnetosphere on Kinetic Scales 

Wieslaw M. Macek and Szymon Gorka

We consider magnetic turbulence using observations from the Magnetospheric Multiscale (MMS) mission on kinetic (ions and electron) scales, which are far shorter than the scales characteristic for description of plasma by magnetohydrodynamic (MHD) theory. We have shown that a break of the magnetic spectral exponent to about -5.5 agrees with the predictions of kinetic theory (-16/3), see Ref. [1]. It is worth noting that the unprecedented very high (millisecond) resolution of the magnetic field instrument allowed to grasp the mechanism of reconnection in the magnetotail on kinetic scales, Ref. [2]. As expected from numerical simulations, we have verified that when the field lines and plasma become decoupled a large reconnecting electric field related to the Hall current (1–10 mV/m) is responsible for fast reconnection in the ion diffusion region both at the magnetopause and in the magnetotail regions. Although inertial accelerating forces remain moderate (1–2 mV/m), the electric fields resulting from the divergence of the full electron pressure tensor provide the main contribution to the generalized Ohm’s law at the neutral sheet (of the order of 10 mV/m), cf. [3]. This illustrates that when ions decouple electron physics dominates. The results obtained on kinetic scales may be useful for better understanding the physical mechanisms governing reconnection processes in various magnetized space and laboratory plasmas.

Acknowledgments. This work has been supported by the National Science Center, Poland (NCN), through grant No. 2021/41/B/ST10/00823.


1. Macek, W. M., Krasinska, A., Silveira, M. V. D., Sibeck, D. G., Wawrzaszek, A., Burch, J. L., & Russell, C. T. 2018, Magnetospheric Multiscale observations of turbulence in the magnetosheath on kinetic scales, Astrophys. J. Lett., 864, L29,

2. Macek, W. M., Silveira, M. V. D., Sibeck, D. G., Giles, B.L., & Burch, J. L. 2019a, Magnetospheric Multiscale mission observations of reconnecting electric fields in the magnetotail on kinetic scales, Geophys. Res. Lett., 46, 10,295—10,302,

3. Macek, W. M., Silveira, M. V. D., Sibeck, D. G., Giles, B.L., & Burch, J. L. 2019, Mechanism of reconnection on kinetic scales based on Magnetospheric Multiscale mission observations, Astrophys. J. Lett., 885, L26,


How to cite: Macek, W. M. and Gorka, S.: Modeling Reconnection and Turbulence in the Magnetosphere on Kinetic Scales, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-6120,, 2022.

EGU22-6171 | Presentations | ST1.1

Two-fluid numerical model towards reproducing the observed solar wave cutoffs. 

Błażej Kuźma, Kris Murawski, and Zdzisław Musielak

With use of JOANNA code we developed a numerical model that describes a partially ionized plasma described by a set of two-fluid equations for ion + electron and neutral fluids, coupled by ion-neutral collisions. We used this model to simulate a quiet region of the solar atmosphere with granules spontaneously generated in the photosphere by the underlying solar convective motions. We found that such ongoing granulation excites a wide range of waves propagating into the upper atmospheric layers, with their cutoff frequencies strongly depending on height above the photosphere. We report for the first time numerically obtained cutoff frequencies that are consistent with the cutoff frequencies computed by Stark & Musielak (1993), Kraśkiewicz et al. (2019) and Wójcik et al. (2019). What is even more important, our results remain in agreement with the observational data of Wiśniewska et al. (2016) and Kayshap et al. (2018). As the exact analytical formula for two-fluid cutoff frequencies has not been found up to date, the numerical simulations are crucial tool to answer the ongoing question about impact of different physical processes on cutoffs and their variation in the solar atmosphere.

How to cite: Kuźma, B., Murawski, K., and Musielak, Z.: Two-fluid numerical model towards reproducing the observed solar wave cutoffs., EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-6171,, 2022.

EGU22-6875 | Presentations | ST1.1 | Highlight

Data-Driven MHD Simulations on Magnetic Flux Rope Eruptions 

Yang Guo, Mingde Ding, Pengfei Chen, Chun Xia, Rony Keppens, Ze Zhong, Jinhan Guo, and Yiwei Ni

Solar eruptions such as flares and coronal mass ejections could cause disastrous space weather. To understand and predict these eruptive activities, we have to combine multi-wavelength observations and numerical simulations. Recently, data-driven magnetohydrodynamic (MHD) simulations have provided a series of new findings in studying the accumulation of electric current and magnetic energy in active regions, in explaining magnetic flux rope eruptions and coronal mass ejections. We briefly review the progress in this field and introduce one way to realize data-driven MHD simulation, including processing magnetic field observational data, inversion of velocity field and electric field, models as initial conditions and subsequent dynamic simulations. Finally, we will look into the future of the data-driven simulations and point out several methods to improve the simulation results. 

How to cite: Guo, Y., Ding, M., Chen, P., Xia, C., Keppens, R., Zhong, Z., Guo, J., and Ni, Y.: Data-Driven MHD Simulations on Magnetic Flux Rope Eruptions, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-6875,, 2022.

EGU22-7233 | Presentations | ST1.1

Three-dimensional particle-in-cell simulation of high-speed plasma jets interacting with Mercury’s magnetosphere 

Gabriel Voitcu, Marius Echim, Eliza Teodorescu, and Costel Munteanu

The transport of high-speed plasma jets (or clouds, streams, blobs, plasmoids) across magnetic discontinuities/shocks is a key process for planetary magnetic environments. Recently, a large number of localized magnetic structures detected in-situ by MESSENGER in the Hermean magnetosheath has been reported in the literature. These structures are similar to the high-speed plasma jets identified in the Earth’s magnetosheath. Due to some limitations in the plasma measurements on-board MESSENGER, only the magnetic signature has been studied. The BepiColombo mission provides a great opportunity to further advance this type of investigation. In this paper we use 3d3v electromagnetic particle-in-cell simulations to study the transport and entry of high-speed plasma jets into Mercury’s magnetosphere. The physical setup is adapted to simulate the kinetic effects and their role on the dynamics of localized plasma structures propagating from the magnetosheath toward the Hermean magnetopause. The magnetospheric field of planet Mercury is provided by the KT17 model, while the high-speed plasma jets are defined as 3D finite-size elements with their bulk velocity pointing towards the dayside magnetopause. We investigate the space and time evolution of the plasma jets prior, during and after their impact on the Hermean magnetopause. We analyse the parallel and perpendicular dynamics with respect to the background magnetic field and emphasize key physical processes for the propagation of high-speed plasma jets across transverse magnetic fields. Our simulations shall support the future exploitation of in-situ data from BepiColombo.

How to cite: Voitcu, G., Echim, M., Teodorescu, E., and Munteanu, C.: Three-dimensional particle-in-cell simulation of high-speed plasma jets interacting with Mercury’s magnetosphere, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-7233,, 2022.

EGU22-7377 | Presentations | ST1.1

Similarity and difference of Alfvén waves' propagation and evolution in the slow and fast solar wind of inner heliosphere 

Jiansen He, Liping Yang, Die Duan, Xingyu Zhu, and Chuanpeng Hou
Parker Solar Probe detected the Alfvénic slow solar wind in the inner heliosphere. Before that, although different spacecraft had measured that slow solar wind sometimes has Alfvénic characteristics, they did not attract extensive and strong attention and discussion. The critical question is, how does the propagation and evolution of Alfvénic waves perform through the same or different processes in the fast and slow solar wind? To study this problem, we simulate the formation of high and low-speed solar wind and the propagation of evolving Alfvén waves therein from a global perspective. Compared with one-dimensional or multi-dimensional simulations with a limited range of latitude and longitude, the advantage of global simulation is that it provides a self-consistent model of fast and slow solar wind coexisting in different flow tubes. Based on this model, we study and evaluate the effects of the expansion, bending, and non-uniformity across the flux tube on the propagation of evolving Alfvén wave. The varying characteristics during the propagation consist of wave amplitude, wave-vector anisotropy, wave mode conversion, etc. As the critical interface to distinguish the sub-Alfvénic and super-Alfvénic solar wind, which is also the vital interface to distinguish the corona and interplanetary space, Alfvén surface is another important aspect of our research. We study the propagation characteristics of Alfvén waves inside and outside the non-spherical interface. In addition, we also discuss the possible relationship between the propagation and evolution of the Alfvén wave and the formation and development of switchback.

How to cite: He, J., Yang, L., Duan, D., Zhu, X., and Hou, C.: Similarity and difference of Alfvén waves' propagation and evolution in the slow and fast solar wind of inner heliosphere, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-7377,, 2022.

EGU22-7843 | Presentations | ST1.1

How are the energetics of transverse waves in a coronal loop affected by a transition region? 

Timothy Duckenfield, Gabriel Pelouze, and Tom Van Doorsselaere

Transverse waves of a coronal loop, also known as kink waves, are commonplace in the solar atmosphere, both in their standing and propagating form. Such waves may be important in the energy + mass cycles of the solar corona. Furthermore, recent numerical studies of coronal loops have shown that the plasma heating from dissipation of these waves is sufficient to overcome radiative cooling. However, the addition of a dense mass reservoir at the end(s) of the loop in the form of a chromosphere and transition region can alter the energetics of the wave and its evolution compared to a purely coronal loop.
In this talk, I will outline current progress in 3D MHD simulations of a solar loop incorporating a chromosphere, transition region and coronal component in a stratified, thermally conducting atmosphere. Transverse waves are induced from a driver in the chromosphere, showing these waves are able to penetrate the transition region. This result is important for the decay-less oscillation regime, in which very small amplitude transverse oscillations are seen to persist for many periods, despite the (presumable) action of damping mechanisms such as resonant absorption. The fact that decay-less oscillations may be driven down in the chromosphere supports the notion that decay-less oscillations are powered from below.
When the loop is sufficiently driven, the motion of the coronal plasma leads to small scales generated from Kelvin Helmholtz instability eddies, and these deformations are regions of enhanced heating. I will discuss the simulation results on how the wave energy is dissipated into heat; the relationship between the driver and the heating; and the extent to which the entire loop is heated, and compare with the purely coronal case. 

How to cite: Duckenfield, T., Pelouze, G., and Van Doorsselaere, T.: How are the energetics of transverse waves in a coronal loop affected by a transition region?, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-7843,, 2022.

EGU22-8175 | Presentations | ST1.1 | Highlight

The magnetopause discontinuity: a MMS study. 

Giulio Ballerini, Gerard Belmont, Laurence Rezeau, and Francesco Califano

The magnetopause boundary seems to escape the general classification of discontinuities since it mixes characteristics of shocks (magnetic field magnitude increase) and those typical for the rotational discontinuities (magnetic rotation), whereas it is very often described as a tangential discontinuity. As, the main issue is the amount of matter/momentum/energy from the solar wind and entering into the magnetosphere, the solution cannot be simply achieved by assuming the discontinuity as strictly tangential, everywhere and at all times. Here we propose to study the magnetopause boundary as a "quasi-tangential" discontinuity, with the normal magnetic component Bn small but not null since even small departures from the standard hypothesis of a zero Bn can lead to noticeable changes in the global properties. In that aim, we look into the MMS database for a large number of magnetopause crossings. For each case we will determine what are the most important features (non-planarity, non-stationarity, Hall effect, pressure anisotropy and agyrotropy) that allow the discontinuity to escape the general classification, i.e. to noticeably change the form of the conservation laws on which the theory of discontinuities is based for non-strictly tangential discontinuities. We put a special emphasis on the refined methods that can be used for determining the spatial gradients from four spacecraft data and on the accuracy that can be attained by these methods.

How to cite: Ballerini, G., Belmont, G., Rezeau, L., and Califano, F.: The magnetopause discontinuity: a MMS study., EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-8175,, 2022.

EGU22-8824 | Presentations | ST1.1

Importance of accurate numerical implementaion of electron inertial terms in hybrid-kinetic simulations  of collisionless plasma turbulence 

Neeraj Jain, Patricio A. Muñoz, Jörg Büchner, Meisam Tabriz, and Markus Rampp

A comprehensive understanding of the turbulent dissipation of magnetic energy in collisionless space and astrophysical plasmas, an unsolved problem yet, requires efficient kinetic simulations of collisionless plasma turbulence. Fully kinetic simulations (where all the plasma species, ions and electrons, are treated kinetically) of collisionless plasma turbulence covering a full range of kinetic scales (from ion to electron scales) are computationally demanding.  Hybrid-kinetic simulations (ions treated as kinetic species and electrons as fluid and therefore ignoring electron kinetic effects) with inertia-less electron fluid, although less demanding computationally,  can not address the physics at the electron scales. Hybrid-kinetic simulations with inertial electron fluid are applicable all the way down to electron scales (still ignoring electron kinetic effects) with computational demands in between those for hybrid-kinetic simulations with inertia-less electron fluid and fully kinetic simulations, and, therefore have begun to attract significant interest recently.

The majority of hybrid-kinetic codes, solving either the Vlasov equation for the ions by an Eulerian method (called Vlasov-hybrid codes) or the equations of motion for ion macro-particles by the Lagrangian Particle-in-Cell (PIC) method (called PIC-hybrid codes), numerically implement the electron inertial terms of  the electron fluid equations under varying approximations which are not necessarily valid at electron scales. In hybrid-kinetic codes, electric field is calculated from either the generalized Ohm's law or an elliptic partial differential equation. In the former case, the non-stationary electron inertial term (time derivative of electron bulk velocity) in the generalized Ohm's law is neglected (approximation A1). In the latter case, a part of the electron inertial term involving cross partial derivatives of electric field is  neglected in comparison to the other part involving second order partial derivatives to obtain a simpler elliptic partial differential equation for electric field (approximation A2). For two dimensional collisionless plasma turbulence, we assess the validity of the two approximations of electron inertial terms used in hybrid-kinetic codes for the calculation of electric field. We employ our recently parallelized three-dimensional PIC-hybrid code CHIEF, which numerically implements the electron inertial terms without any of these approximations, in a quasi-two dimensional setup for the simulations of collisionless plasma turbulence.  We find that the approximation A1 impacts the accuracy of the results at the electron scales and therefore may lead to physically incorrect results. Approximation A2, on the other hand, is found to be invalid from ion to electron scales. We conclude that the simulation results obtained using the hybrid-kinetic codes with an approximate numerical implementation of the electron inertial term need to be interpreted with caution.  

How to cite: Jain, N., Muñoz, P. A., Büchner, J., Tabriz, M., and Rampp, M.: Importance of accurate numerical implementaion of electron inertial terms in hybrid-kinetic simulations  of collisionless plasma turbulence, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-8824,, 2022.

EGU22-8919 | Presentations | ST1.1

Modelling the propagation of slow magnetoacoustic waves in a multi-stranded coronal loop 

Krishna Prasad Sayamanthula and Tom Van Doorsselaere

The cross-field thermal structure of a coronal loop is one of the critical parameters useful to distinguish the major heating theories. In a recent study we are able to isolate two thermal components of a coronal loop using observations of propagating slow magnetoacoustic waves in two different temperature channels. In order to properly interpret these observations and identify the actual cross-field thermal structure, we develop multiple three-dimensional magnetohydrodynamic models of coronal loop. In each of these models slow magnetoacoustic waves are driven by perturbing the plasma at one end and the corresponding multi-wavelength propagation characteristics are studied by applying forward modelling techniques. We compare the wave propagation properties between different models including a monolithic and a multi-stranded model with those from observations to draw some important inferences which will be discussed in this talk.

How to cite: Sayamanthula, K. P. and Van Doorsselaere, T.: Modelling the propagation of slow magnetoacoustic waves in a multi-stranded coronal loop, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-8919,, 2022.

EGU22-10585 | Presentations | ST1.1 | Highlight

Properties of cosmic ray test particles in global MHD simulation of the heliosphere 

Shuichi Matsukiyo, Kotaro Yoshida, Haruichi Washimi, and Tohru Hada

Heliosphere is a bubble occupied by the solar wind plasma and magnetic field in the local interstellar space. The motion of galactic cosmic rays (GCRs) invading into the heliosphere are strongly affected by the electromagnetic structures of the heliosphere. The statistical behavior of the GCRs near and inside the heliosphere have been conventionally studied by many authors using the diffusion convection model [e.g., Moraal (2013)].

  In this study we investigate the behavior of GCRs invading into the heliosphere in the level of particle trajectory. We conduct test particle simulations of GCRs by using the electromagnetic fields obtained from a global MHD simulation of the heliosphere. The MHD simulation assumes steady solar wind and interstellar wind. GCR protons are initially distributed outside the heliosphere and their motions in the steady virtual heliosphere are calculated by using the Buneman-Boris method. Depending on their initial energy, various types of particle motions, current sheet drift, polar drift, spiral motion, shock drift, Fermi-like acceleration, linear motion, resonantly scattered motion, mirror reflection by bottleneck interstellar field, are observed. We further discuss some statistics of the particles reached at the inner boundary (=50AU from the sun) of the simulation domain.

How to cite: Matsukiyo, S., Yoshida, K., Washimi, H., and Hada, T.: Properties of cosmic ray test particles in global MHD simulation of the heliosphere, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-10585,, 2022.

The solar wind has historically been considered as a far simpler medium than the solar or magnetospheric plasma. From the beginning of the space era through to XXI century, it has been supposed that all solar wind structures freely expand and can be modeled in 2D. Correspondingly, current sheets have been pictured as thin planar objects carrying the electric current and separating magnetic fields of different directions. Before 2010th, Petscheck magnetic reconnection was considered as the only possible mechanism transforming the magnetic energy into the thermal energy at current sheets in the solar wind. Acceleration of particles to suprathermal energies was thought to be impossible there because of too slow inflow speeds. Interpretations of observations totally followed the theoretical dominant paradigm mainly because of the insufficiency of observational material. Things began to change when more and more theoretical and observational studies in magnetospheric and solar physics appeared pointing to the complex character of magnetic reconnection. In particular, ideas about stochastic or turbulent reconnection at current sheets in realistic space plasmas become dominating. In turn, observations of the fine structure of current sheets in the solar wind as well as evidence for local acceleration of energetic particles found with help of modern missions, including Parker Solar Probe, allow re-considering views on solar wind current sheets and better understanding physics of the processes associated with magnetic reconnection.   

How to cite: Khabarova, O.: Magnetic reconnection and particle acceleration in the solar wind: theory, observations and opinions, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-10641,, 2022.

EGU22-10807 | Presentations | ST1.1 | Highlight

Homologous Coronal Mass Ejections Caused by Recurring Formation and Disruption of Current Sheet within a Sheared Magnetic Arcade 

Chaowei Jiang, Xinkai Bian, and Xueshang Feng

The Sun often produces coronal mass ejections with similar structure repeatedly from the same source region, and how these homologous eruptions are initiated remains an open question. Here, by using a new magnetohydrodynamic simulation, we show that homologous solar eruptions can be efficiently produced by recurring formation and disruption of coronal current sheet as driven by continuously shearing of the same polarity inversion line within a single bipolar configuration. These eruptions are initiated by the same mechanism, in which an internal current sheet forms slowly in a gradually sheared bipolar field and reconnection of the current sheet triggers and drives the eruption. Each of the eruptions does not release all the free energy but with a large amount left in the post-flare arcade below the erupting flux rope. Thus, a new current sheet can be more easily formed by further shearing of the post-flare arcade than by shearing a potential field arcade, and this is favorable for producing the next eruption. Furthermore, it is found that the new eruption is stronger since the newly formed current sheet has a larger current density and a lower height. In addition, our results also indicate the existence of a magnetic energy threshold for a given flux distribution, and eruption occurs once this threshold is approached.

How to cite: Jiang, C., Bian, X., and Feng, X.: Homologous Coronal Mass Ejections Caused by Recurring Formation and Disruption of Current Sheet within a Sheared Magnetic Arcade, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-10807,, 2022.

EGU22-10833 | Presentations | ST1.1 | Highlight

Can the ultraviolet bursts be generated in the low solar chromosphere? 

Lei Ni

Ultraviolet (UV) bursts and Ellerman bombs (EBs) are small scale magnetic reconnection events in the highly stratified low solar atmosphere. The plasma density, reconnection mechanisms and radiative cooling/transfer process are very different at different atmospheric layers. EBs are believed to form in the up photosphere or low chromosphere. It is still not clear whether UV bursts have to be generated at a higher atmospheric layer than the EBs or UV bursts and EBs can actually both appear in the low chromosphere.   We numerically studied the low beta magnetic reconnection process around the solar temperature minimum region (TMR) by including more realistic physical diffusions and radiative cooling models. We aim to find out if UV bursts can be formed in the low chromosphere and investigate the dominant mechanism that transfer the magnetic energy to heat in an UV burst.The single-fluid MHD code NIRVANA was used to perform simulations. The time dependent ionization degrees of Hydrogen and Helium are included in the code, which lead to the more realistic magnetic diffusion caused by electron-neutral collision and ambipolar diffusion. The more realistic radiative cooling model from Carlsson& Leenaarts 2012 is included in the simulations. The high resolution simulation results indicate that the plasmas in the reconnection region can be heated above 20,000 K as long as the reconnection magnetic fields reach above 500 G, which further proves that UV bursts can be generated in the dense low chromosphere. When the reconnection magnetic fields are stronger than 900 G, the width of the synthesized Si IV 1394 A line profile with multiple peaks can reach above 100 km s-1, which is consistent with the usually observed broad-line-width UV bursts. The dominate mechanism that converts magnetic energy to heat in an UV burst in the low chromosphere is Joule heating that is contributed by magnetic diffusion caused by electron-ion collision in the reconnection region. The average power density that is converted to the thermal energy in the reconnection region is about 100 erg cm-3 s-1, which is comparable to the average power density of the released heat in an UV burst.


How to cite: Ni, L.: Can the ultraviolet bursts be generated in the low solar chromosphere?, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-10833,, 2022.

EGU22-10945 | Presentations | ST1.1

Doppler shifts of spectral lines formed in the solar transition region and corona 

Yajie Chen, Hardi Peter, Damien Przybylski, Hui Tian, and Jiale Zhang

Emission lines formed in the transition region and corona show dominantly redshifts and blueshifts, respectively. We investigate the Doppler shifts in a 3D radiation magnetohydrodynamic (MHD) model of the quiet Sun and compare these to observed properties. We concentrate on Si IV 1394 Å originating in the transition region and examine the Doppler shifts of several other spectral lines at different formation temperatures. We construct a radiation MHD model extending from the upper convection zone to the lower corona using the MURaM code. In this quiet Sun model the magnetic field is self-consistently maintained by the action of a small-scale dynamo alone. We synthesize the profiles of several optically thin emission lines, formed at temperatures from the transition region into the corona. We investigate the spatial structure and coverage of red- and blueshifts and how this changes with line-formation temperature. The model successfully reproduces the observed change of average net Doppler shifts from red- to blueshifted from the transition region into the corona. In particular, the model shows a clear imbalance of area coverage of red- vs. blueshifts in the transition region of ca. 80% to 20%, even though it is even a bit larger on the real Sun. We isolate that (at least) four processes generate the systematic Doppler shifts in our model, including pressure enhancement in the transition region, transition region brightenings unrelated to coronal emission, boundaries between cold and hot plasma, and siphon-type flows. We show that there is not the single process that is responsible for the observed net Doppler shifts in the transition region and corona. Because current 3D MHD models do not yet fully capture the evolution of spicules, one of the key ingredients of the chromosphere, most probably these have still to be added to the list of processes responsible for the persistent Doppler shifts.

How to cite: Chen, Y., Peter, H., Przybylski, D., Tian, H., and Zhang, J.: Doppler shifts of spectral lines formed in the solar transition region and corona, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-10945,, 2022.

EGU22-11960 | Presentations | ST1.1

A new global nonlinear force-free coronal magnetic-field extrapolation code implemented on a Yin Yang grid. 

Argyrios Koumtzis and Thomas Wiegelmann

The solar magnetic field dominates and structures the
coronal plasma and detailed insights are important to understand almost
all physical processes. While direct routine measurements
of the coronal magnetic field are not available, we have to extrapolate
the photospheric vector field measurements into the corona. To do so,
we developed a new code that performs state-of-the-art nonlinear force-free
magnetic field extrapolations in spherical geometry.
Our new implementation is based on an optimization principle and is
able to reconstruct the magnetic field in the entire corona, including
the polar regions.
Because of the nature of the finite-difference numerical scheme used in
the past, extrapolation close to polar regions was computationally
inefficient. In the current code, the so-called Yin Yang grid is used.
Both the speed and accuracy of the code is improved compared to previous
implementations. We tested our new code with a well known
semi-analytical model (Low and Lou solution). This new Yin and Yang
implementation is timely because
the Solar Orbiter mission is expected to provide reliable
vector magnetograms also in the polar regions within the following years.
Thus, this code can be used in the future when these synoptic magnetograms
are available to model the magnetic field of the solar corona for the entire
Sun including the polar regions.

How to cite: Koumtzis, A. and Wiegelmann, T.: A new global nonlinear force-free coronal magnetic-field extrapolation code implemented on a Yin Yang grid., EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-11960,, 2022.

EGU22-12161 | Presentations | ST1.1 | Highlight

Excitations of Alfven Wave by 3D Patchy and Intermittent Magnetic Reconnection 

Liping Yang, Jiansen He, Xueshang Feng, and Ming Xiong

Alfven waves make a central role in energy transfer of the solar atmosphere and heliosphere, with the potential to heat corona, accelerate solar wind, and drive Alfvenic turbulence. It has long been suggested that magnetic reconnection can generate Alfven waves through a relaxation of a highly curved reconnected field lines. Here, with a high-resolution simulation of 3D magnetic reconnection under the solar corona environment, we study this sketch and find that Alfven waves, whose features resemble to those of the Alfven Waves observed in the solar atmosphere, are continually and energetically excited by reconnection mainly through two ways. One involves the current sheet which experiences patchy and intermittent reconnections, and the other refers to the turbulence forming in the outflow regions. The Pointing flux carried by the excited upward-propagating Alfven waves can satisfy the requirements of plasma heating in the corona. This has implications for self-consistent considerations of energy budgets in the solar atmosphere and heliosphere.

How to cite: Yang, L., He, J., Feng, X., and Xiong, M.: Excitations of Alfven Wave by 3D Patchy and Intermittent Magnetic Reconnection, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-12161,, 2022.

EGU22-12196 | Presentations | ST1.1

On the evolution equations of field gradient tensor invariants in MHD Theory 

Virgilio Quattrociocchi, Giuseppe Consolini, Massimo Materassi, Tommaso Alberti, and Ermanno Pietropaolo

The understanding of the evolution of turbulence in space plasmas requires the investigation of magnetic field and plasma structures and their time evolution. Indeed, the interactions among different topological multi scale structures have been shown to play an important role tu understand the energy transfer across scales and dissipation. A possible approach to this issue is the study of the geometrical invariants of magnetic and velocity field gradient tensors from a Lagrangian point of view. In the early 1980 a series of works (Vielliefosse, 1982 and 1984) discussed a nonlinear homogeneous evolution equation for the velocity gradient tensor in fluid dynamics.  Here, we derive the evolution equations of the geometrical invariants of the magnetic and velocity field gradient tensors in the case of magneto-hydrodynamic theory and discuss their application to the analysis of magneto-hydrodynamic turbulence in space plasmas.

How to cite: Quattrociocchi, V., Consolini, G., Materassi, M., Alberti, T., and Pietropaolo, E.: On the evolution equations of field gradient tensor invariants in MHD Theory, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-12196,, 2022.

EGU22-12533 | Presentations | ST1.1

MHD wave modes of solar magnetic flux tubes with the realistic cross-section 

Anwar Aldhafeeri, Viktor Fedun, Istvan Ballai, and Gary Verth


Current and near-future high resolution solar observations indicate that theoretical modelling of the magnetohydrodynamic (MHD) modes in magnetic waveguides with realistic structure and shape now becomes an imperative necessity. 

It was recently shown that even a magnetic structure with elliptical shape, that corresponds to the weak irregularity, may significantly influence the spatial structure of MHD mode in comparison to the mode structure obtained from the modelling which is based on the magnetic flux tube shape with cylindrical cross-section (Aldhafeeri et al., ApJ, 2021). An inaccurate model used for describing waves may lead to the misinterpretation of observational data. 

The expressions for the linear MHD perturbations of a magnetic flux tube are derived by assuming zero value of the vertical component of the velocity perturbation at the boundary of the magnetic flux tube, which is in good agreement with observations. The governing equation for the vertical velocity perturbation was solved by taking into account the observed realistic shape of the sunspot umbra. With these conditions the proposed model is applicable for the analysis of slow body modes under photospheric conditions. 

Our results show that under solar photospheric conditions the conditions of continuity of the component of radial velocity and pressure at the boundary are enough to be imposed, enabling us to use Cartesian coordinates with varsity numerical methods to model the MHD modes with their realistic cross-sectional shape.


How to cite: Aldhafeeri, A., Fedun, V., Ballai, I., and Verth, G.: MHD wave modes of solar magnetic flux tubes with the realistic cross-section, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-12533,, 2022.

EGU22-12904 | Presentations | ST1.1

Three-dimensional Particle-in-Cell (PIC) simulations of non-dipolar minimagnetospheres and comparison to the experiment. 

Andrey Divin, Ildar Shaikhislamov, Igor Paramonik, Marina Rumenskikh, Daniil Korovinskiy, and Jan Deca

Magnetospheres are formed when plasma flow interacts with an external source of the magnetic field. Objects (with a size comparable to or less than the ion inertial length or ion gyroradius) formed by relatively weak magnetic field sources are called minimagnetospheres since they are rather different from more common large planetary magnetospheres. 

Moon surface has regions called Lunar Magnetic Anomalies (LMAs) where the remanent magnetization of the Lunar crust provides sources of the magnetic field strong enough to stand off the solar wind. These fields produce minimagnetospheres of the size of the several ion inertial lengths (or gyroradii) or below, typically having weakly structured topology and mostly non-dipolar nature. In this study, we combine numerical simulations and laboratory experiments to investigate ion scattering and basic properties of a non-dipolar minimagnetosphere.  A series of laboratory experiments were carried out on the KI-1 facility (Novosibirsk, Russia) to investigate minimagnetosphere properties for the case of a quadrupolar magnetic field source. The experiment consists of a vacuum chamber, of 5 m length and 1.2 m diameter (with a residual pressure of ~10-7 Torr) filled with a moving plasma. The quadrupolar magnetic field is generated by two coils connected in anti-parallel. The experimental results are supported by the Particle-in-Cell (PIC) three-dimensional simulations using code iPIC3D which capture the full kinetic behavior of the interaction. 

We report several important results based on both the experiment and numerical simulations: 1) a majority of particles is reflected by the Hall electric field formed due to the formation of the magnetopause electron current; however, a hotter portion of the inbound distribution also experiences magnetic deflection closer to the B field source; 2) reflecting electrostatic potential is smaller in the quadrupolar case (if compared to dipolar minimagnetosphere); 3) numerical simulations reproduce well the ion reflection pattern seen in the laboratory experiment, but simulations show slightly less reflected ions which might be attributed to unsteady processes developing during each laboratory run.

How to cite: Divin, A., Shaikhislamov, I., Paramonik, I., Rumenskikh, M., Korovinskiy, D., and Deca, J.: Three-dimensional Particle-in-Cell (PIC) simulations of non-dipolar minimagnetospheres and comparison to the experiment., EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-12904,, 2022.

EGU22-915 | Presentations | NP6.3

Quantifying the electron scattering by electrostatic fluctuations in the Earth’s bow shock 

Sergey Kamaletdinov, Ivan Vasko, Anton Artemyev, and Rachel Wang

Collisionless shocks are known to be natural sources of suprathermal particles, but the mechanism resulting in acceleration of thermal electrons to suprathermal energies still remains elusive. The problem is, that the Diffusive Shock Acceleration (DSA) becomes efficient only for suprathermal electrons, which fluxes in the far upstream region are relatively low. Recent studies have shown that the so-called Stochastic Shock Drift Acceleration (SSDA) mechanism can potentially provide the necessary pre-acceleration of incoming thermal electrons to suprathermal energies. In this mechanism, electrons are temporarily kept trapped in the shock transition region due to magnetic mirror reflection by the magnetic ramp and pitch-angle scattering of electrons trying to escape upstream by wave turbulence. Spacecraft measurements showed that broadband electrostatic turbulence is always present in the Earth’s bow shock, but its efficiency in scattering suprathermal electrons has not been estimated up to date. In this study we have quantified the electron scattering by the broadband electrostatic turbulence and, specifically, by electrostatic solitary waves (ion holes) substantially contributing to this turbulence in the Earth’s bow shock. Adopting the solitary wave and turbulence parameters typical of the Earth’s bow shock, we obtain quasi-linear scattering rates and compare these scattering rates to the results of test-particle simulations. This analysis showed that scattering of suprathermal electrons by the osberved electrostatic turbulence is relatively well estimated by the quasi-linear approach. We estimated the quasi-linear scattering rates at various energies and pitch-angles and demonstrated that the electrostatic turbulence in the Earth’s bow shock can provide pre-acceleration of thermal electrons from a few tens of eV to a few hundred eV via the SSDA mechanism.

This work was supported by the Russian Scientific Foundation, Project No. 19–12-00313

How to cite: Kamaletdinov, S., Vasko, I., Artemyev, A., and Wang, R.: Quantifying the electron scattering by electrostatic fluctuations in the Earth’s bow shock, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-915,, 2022.

EGU22-1485 | Presentations | NP6.3

Plasmoid-dominated turbulent reconnection in symmetric and asymmetric systems 

Seiji Zenitani, Momoka Yamamoto, and Takahiro Miyoshi

In magnetohydrodynamics (MHD), magnetic reconnection has been discussed by three theoretical models: Sweet--Parker reconnection, Petschek reconnection, and plasmoid-dominated turbulent reconnection. Among these models, properties of plasmoid-dominated reconnection remain unclear, because it was discovered only recently. In this talk, we explore basic properties of plasmoid-dominated reconnection in a low-beta plasma such as in a solar corona, by using large-scale MHD simulations [1]. We have found that the system becomes highly complex due to repeated formation of plasmoids and shocks. We have further found that the reconnection rate goes higher than previously thought. Next we explore influence of asymmetry in background plasma densities in plasmoid-dominated reconnection. We have found that the average reconnection rate follows Cassak-Shay's hybrid relation [2]. Many signatures become asymmetric across the reconnection layer, and plasmas inside the plasmoids start to swirl in specific directions. Formation processes of these vortices and a potential extension of our numerical survey will be discussed.

[1] S. Zenitani and T. Miyoshi, Astrophys. J. Lett., 894, L7 (2020)
[2] P. A. Cassak and M. A. Shay, Phys. Plasmas, 14, 102114 (2007)

How to cite: Zenitani, S., Yamamoto, M., and Miyoshi, T.: Plasmoid-dominated turbulent reconnection in symmetric and asymmetric systems, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-1485,, 2022.

EGU22-1577 | Presentations | NP6.3

Physical Regimes of 2D MHD turbulent reconnection in different Lundquist numbers 

Haomin Sun, Yan Yang, Quanming Lu, San Lu, Minping Wan, and Rongsheng Wang

Using two-dimensional (2D) MHD simulations in different Lundquist numbers , we investigate physical regimes of turbulent reconnection and the role of turbulence in enhancing the reconnection rate. Turbulence is externally injected into the system with varying strength. External driven turbulence contributes to the conversion of magnetic energy to kinetic energy flowing out of the reconnection site and thus enhances the reconnection rate. The plasmoids formed in high Lundquist numbers contribute to the fast reconnection rate as well. Moreover, an analysis of the power of turbulence implies its possible association with the generation of plasmoids. Additionally, the presence of turbulence has great impact on the magnetic energy conversion and may be favorable for the Kelvin-Helmholtz (K-H) instability in the magnetic reconnection process.

How to cite: Sun, H., Yang, Y., Lu, Q., Lu, S., Wan, M., and Wang, R.: Physical Regimes of 2D MHD turbulent reconnection in different Lundquist numbers, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-1577,, 2022.

EGU22-1740 | Presentations | NP6.3

Local slope of magnetic field power spectrum in inertial and kinetic ranges of solar wind turbulence 

Alexander Pitna, Jana Safrankova, Zdenek Nemecek, Gilbert Pi, Luca Franci, and Byeongseon Park

Solar wind, a supersonic flow of plasma embedded in the magnetic field, exhibits turbulent behavior. The character of turbulent fluctuations has been investigated through low cadence measurements of particle distribution function and high cadence magnetic field measurements. One of the most frequently adopted approach in the analysis of the ‘measured’ time series of any particular quantity is the estimation of its power spectral density (PSD). The shape of the PSD then may infer which physical mechanisms govern the evolution of turbulent fluctuations. Generally, every ‘measured’ time series is ‘noisy’ and it differs from the ‘true’ one (measured by an ideal instrument). In turn, the shape of PSD is affected as well. In this paper, we focus on a special case where the signal and noise are independent, i.e., the noise is additive and therefore, the PSD of measured signal can be expressed as a sum of ‘true’ and ‘noise’ PSDs. Moreover, we define a so-called local slope in the framework of continuous wavelet transform as the finite difference derivative between the two consequent values of a global PSD. Employing this technique, we show that the noise of magnetic field measurements of the MFI instrument on board the Wind spacecraft is additive. Finally, we applied the technique to measurements of the Parker Solar Probe close to the Sun. Our preliminary results suggest that our technique may lead to a more accurate estimations of the kinetic range spectral indices.

How to cite: Pitna, A., Safrankova, J., Nemecek, Z., Pi, G., Franci, L., and Park, B.: Local slope of magnetic field power spectrum in inertial and kinetic ranges of solar wind turbulence, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-1740,, 2022.

The measurements of the longitudinal velocity were performed in an open-circuit suction wind tunnel installed at the laboratory of the Max-Planck Institute for Dynamics and Self-Organization in Gottingen, using hot wire anemometer at different positions in turbulent flow generated by a traditional fractal square grid (FSG) and by a spaced fractal square grid (SFSG) with similar physical properties have shown that the self-similarity is present. The statistical description of this complex turbulent system was performed using Extended Self Similarity (ESS). We propose a complementary methodology suitable for non-homogeneous turbulence based on the analysis of the energy transfer hierarchy. The signature of the non-homogeneous characteristics of a turbulent field, indicated by nonlocal dynamics, is separated from those usually assigned as being only due to the intermittency. We propose a physical interpretation of the observed scale independence of the relative scaling exponents in such non-homogeneous flows by means of the compensation effect of the energy transfer on the difference between the strong coherent turbulent events and the background less intense turbulence. This procedure is able to distinguish whether the intermittency arises from the small scales or is linked to coherent structures. The practical interest of this type of turbulent excitation concerns several fields of aeronautical and space application and energy or environmental problems of noise reduction of mixers in combustion or for the numerical models of prediction of the dispersion of pollutants in the atmosphere.

How to cite: Ben Mahjoub, O. and Ouadoud, A.: Intermittency in turbulence generated by traditional fractal square grid and spaced fractal square grid, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-1997,, 2022.

Turbulence, the self-generated turbulence by plasmas and magnetic field collective interaction, has been found to play an essential role in energizing charged particles in the large-scale reconnecting current sheet in the major solar eruption.

The typical large-scale CME/Flare events involve sudden bursts of particle acceleration from the sudden release of magnetic energy in a few minutes to a few tens of minutes. The X-rays emission and gamma rays burst produced by the combined result from the interactions of electrons, hydrogen, helium, and other heavier ions. Space and laboratory researchers are more inclined to believe that turbulence acceleration is belonged to shock acceleration. Solar and astrophysics researchers are more inclined to believe that turbulence acceleration is an independent acceleration mechanism that belongs to the flare acceleration. The evidence in both theories and observations from solar atmosphere activities shows that the acceleration is related to nonlinear resonant wave-particle interaction (e.g., Landau acceleration). So far, many-particle acceleration models consider turbulence acceleration as an effective way of generating energetic electrons, protons, and heavier ions. However, the detailed role of turbulence in this process remains unclear. More effort needs to invest in looking into particle accelerations by turbulence that occurs over a large range of the scale in space from the inertial scale of individual particles to the MHD scale.

In this work, applying the statistical treatment of plasma physics, combing with filter theory of turbulence, the actual ratio of the proton mass to the electron mass, and mass-to-charge ratios, we investigate the interaction of charged particles with the turbulent electric field and magnetic field in the large-scale CME/flare current sheet by applying the

We found the significant Langmuir turbulence acceleration (LTA) through the nonlinear resonant wave-particle interaction in the diffusion region via tracking the trajectories and analyzing the energy spectrum of energetic protons and electrons. The results show that protons and electrons could be efficiently accelerated simultaneously and that the way of LTA is similar to that of the shock acceleration}} but is much more efficient than the shock acceleration. This indicates that large-scale reconnection is a good candidate for the mechanism for the efficient acceleration of protons and electrons in the major solar eruption.

The acceleration of heavy ion considered Helium (3He/4He) and other heavy elements in 3He-rich flares burst would explore in the follow-up work series.


How to cite: Zhu, B., Li, Y., and Lin, J.: Investigations of Particle Accelerations by Turbulent Magnetic Reconnection in Large-Scale CME/Flare Current Sheet: I. Protons and Electrons, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-2116,, 2022.

EGU22-2649 | Presentations | NP6.3

Electromagnetic energy conversion by various processes in turbulent plasmas observed by MMS 

Thanapon Aiamsai, Peera Pongkitiwanichakul, Rungployphan Kieokaew, David Ruffolo, and Theerasarn Pianpanit

A key issue in space plasma physics is how electromagnetic energy is converted to plasma particle energy and heat. Electromagnetic energy conversion generally involves turbulence and/or instabilities. With the Magnetospheric Multiscale (MMS) mission data, we investigate such energy conversion in turbulent plasmas, separating the plasma currents from various drift motions and other processes and assessing their contributions. For example, we have explored the roles of curvature drift, gradient drift, particle inertia drift and perpendicular magnetization currents. We will discuss their roles and related mechanisms in turbulent plasmas. This research has been supported in part by grant RTA6280002 from Thailand Science Research and Innovation, by DPST scholarship grant ,and by grant RGNS 63-045 from Office of the Permanent Secretary, Ministry of Higher Education, Science, Research and Innovation.

How to cite: Aiamsai, T., Pongkitiwanichakul, P., Kieokaew, R., Ruffolo, D., and Pianpanit, T.: Electromagnetic energy conversion by various processes in turbulent plasmas observed by MMS, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-2649,, 2022.

EGU22-3198 | Presentations | NP6.3

Instabilities and Turbulence in Two-Dimensional Magnetohydrodynamics 

Jonathan Tessier, Francis J. Poulin, and David W. Hughes

The solar tachocline is a dynamically important thin region in the Sun, located between the convective and radiative zones and characterised by strong shear in both the radial and latitudinal directions. Furthermore, it is believed to play a key role in the solar dynamo process through the shearing of a poloidal field into a stronger toroidal component. Motivated by the dynamics of the tachocline we have conducted a detailed numerical exploration of the dynamics of sheared MHD turbulence.

Specifically, we have implemented a parallelized numerical model using the "shenfun" Python library to solve the nonlinear two-dimensional Magnetohydrodynamic (MHD) equations to study the dynamics of unstable jets and turbulence in astrophysical plasmas. In particular, we study details of how the jet becomes unstable and the resulting cascade of energy in the case of MHD turbulence. In addition to studying the evolution of the physical quantities, we also investigate the evolution of the spectral slopes and spectral fluxes. As has been found in previous studies of MHD turbulence, a very weak large-scale magnetic field can play a key dynamical role through its amplification on small scales. For extremely weak fields, the behaviour is essentially hydrodynamic. However, once the field is dynamic, the nature of the resulting MHD solution is very different. We are able to classify the various flows and quantify the nature of the solutions in the two regimes.

How to cite: Tessier, J., Poulin, F. J., and Hughes, D. W.: Instabilities and Turbulence in Two-Dimensional Magnetohydrodynamics, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-3198,, 2022.

EGU22-3354 | Presentations | NP6.3

High-Speed Jets in Earth’s Magnetosheath Downstream of the Quasi-Parallel Shock: A Two-Dimensional Global Hybrid Simulation 

Jin Guo, San Lu, Quanming Lu, Yu Lin, Xueyi Wang, Yufei Hao, Kai Huang, Rongsheng Wang, and Xinliang Gao

High-speed jets (HSJs) occur frequently in Earth’s magnetosheath downstream of the quasi-parallel bow shock. They have great impacts on the magnetosheath and the magnetosphere. Using a two-dimensional global hybrid simulation, we investigate the formation and evolution of the HSJs with an IMF cone angle of 0°. The quasi-parallel shock is near the subsolar point, and the HSJs begin to appear in the quasi-parallel magnetosheath with a parallel (perpendicular) scale size of about 1RE (0.2RE). These HSJs then converge, leading to the formation of a large-scale HSJ with a parallel (perpendicular) scale size of 6RE (1.2RE). Some long HSJs, with a large parallel but small perpendicular scale size, are formed at the quasi-parallel bow shock and extend toward the quasi-perpendicular magnetosheath along with the background magnetosheath flow. Moreover, these long HSJs can cause filamentary structures in the magnetosheath.

How to cite: Guo, J., Lu, S., Lu, Q., Lin, Y., Wang, X., Hao, Y., Huang, K., Wang, R., and Gao, X.: High-Speed Jets in Earth’s Magnetosheath Downstream of the Quasi-Parallel Shock: A Two-Dimensional Global Hybrid Simulation, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-3354,, 2022.

EGU22-3546 | Presentations | NP6.3

Overshoot dependence on the shock parameters 

Natalia Borodkova, Olga Sapunova, Victor Eselevich, Georgy Zastenker, and Yury Yermolaev

The structure of the solar wind plasma flow downstream of the ramp of the interplanetary and bow shocks was studied based on the BMSW plasma spectrometer installed onboard the SPEKTR-R spacecraft. Particular attention was paid to the overshoot region, where correlated oscillations of the ion flux and magnetic field are observed. They are formed by two populations of ions: the inflowing solar wind and the beam of coherent gyrating ions. Based on the statistical analysis it was shown that overshoots form both in supercritical and subcritical shocks. It is found that maximum values of the overshoot amplitudes are significantly influenced by the angle between the shock normal and magnetic field vectors, Mach number, plasma and magnetic compression at the shock front. It was established that the oscillation wavelength determined from the magnetic field measurements onboard the WIND spacecraft, on average, coincides with the oscillation wavelength determined from the ion flux on the SPEKTR-R, while the rates of relaxation of these oscillations can greatly differ. It was also shown that the estimates of the overshoot wavelength good correlate with the convected ion gyroradius.

How to cite: Borodkova, N., Sapunova, O., Eselevich, V., Zastenker, G., and Yermolaev, Y.: Overshoot dependence on the shock parameters, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-3546,, 2022.

Previous numerical works on electron/ion foreshocks observed upstream of a curved shock have been already performed within a self-consistent approach based on 2D PIC simulation (Savoini et Lembege, 2010, 2013, 2015), but are restricted to a supercritical regime only. Present two dimensional PIC (Particle in cell) simulations are used in order to analyze the features of a curved shock and associated foreshocks in a subcritical regime. In order to investigate the dynamic of each electron and ion backstreaming populations, we used test-particles in a pre-computed electromagnetic field (issued from 2D PIC simulations) which allows us to define precisely the characteristic of each population in terms of initial velocity and/or their upstream position to the  θBn angle (angle between the local shock normal and the interplanetary magnetic field IMF). Then, results allow to clarify the following questions: what is the impact of the subcritical regime (i) on the persistence of each electron/ion foreshock respectively ?, (ii) in the case the persistence is confirmed,  how the location (along the curved front) and the angular direction of each foreshock edge are affected ?, and (iii) how the mapping of upstream local  distribution functions are impacted ? Preliminary results will be presented and compared with those already obtained for a supercritical shock.

How to cite: Savoini, P. and Lembege, B.: Analysis of a curved shock front microstructures and associated electron/ion foreshock for a subcritical shock regime, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-3631,, 2022.

EGU22-3868 | Presentations | NP6.3

Model of KAW contribution to cross-magnetopause ion transport 

Alexander Lukin, Anton Artemyev, and Anatoly Petrukovich

Magnetosheath ion transport across the night-side magnetopause can be contributed by ion cross-field diffusion due to wave-particle scattering. In this presentation we focus on such scattering mechanism for the most intense magnetosheath wave emission, kinetic Alfven waves (KAWs). These waves carry a finite field-aligned electric field and potentially can accelerate particles along magnetic field lines. In the fast plasma flows these waves are usually observed as a wide Doppler-shifted electromagnetic spectrum characterized by strong electric fields in high wave-number range. Dense frequency spectrum leads to overlapping of particles resonances with waves and causes particle diffusion in pitch-angle and energy space. We investigate particles diffusion caused by interactions with KAW turbulence in a realistic model of the Earth flank magnetopause with nonuniform ambient magnetic field fitting the tangential discontinuity. The KAW spectrum is determined by a sum of a several thousand plane waves with different frequencies and propagation angles. We estimate diffusion coefficients as function of ion pitch-angle and energy for different distances from the magnetopause and discuss the expected cross-field transport rate for this model.

How to cite: Lukin, A., Artemyev, A., and Petrukovich, A.: Model of KAW contribution to cross-magnetopause ion transport, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-3868,, 2022.

EGU22-4040 | Presentations | NP6.3

Spectral features and energy cascade of kinetic scale plasma turbulence 

Giuseppe Arrò, Francesco Califano, and Giovanni Lapenta

Solar wind (SW) in situ observations of plasma turbulence show that the turbulent magnetic field spectrum follows a Kolmogorov-like scaling ∼k-5/3 at large MHD scales and steepens at ion scales where a different power law develops with a scaling exponent varying between -2 and -4, depending on SW conditions. Recent satellite measurements revealed the presence of a second spectral break around electron scales where the magnetic field spectrum shows an exponential falloff described by the so called exp model ∼k-8/3exp(-ρek), where ρe is the electron gyroradius. This model was tested on a large number of magnetic spectra at various distances from the Sun (from 0.3 to 1 AU) and appears to be a solid feature of turbulent magnetic field fluctuations at kinetic scales [1]. 

Using a fully kinetic energy conserving particle-in-cell (PIC) simulation of freely decaying plasma turbulence we study the spectral properties of the turbulent cascade at kinetic scales. Consistently with satellite observations, we find that the magnetic field spectrum follows the kexp(-λk) law at sub-ion scales, with an exponential range developing around kρe≈1. The same exponential falloff is observed also in the electron velocity spectrum but not in the ion velocity spectrum that drops like a power law without reaching electron scales. We investigate the development of these spectral features by analyzing the high-pass filtered electromagnetic work J·E and pressure-strain interaction -P:∇u of both the ions and the electrons. Our analysis shows that the magnetic field dynamics at kinetic scales is mainly driven by the electrons that are responsible for the formation of the exponential range. In particular, we see that at fully developed turbulence the magnetic field energy is dissipated by a two-stage mechanism lead by the electrons that first subtract energy from the magnetic field and then convert it into internal energy at electron scales through the pressure-strain interaction, that accounts for the electron heating [2].



[1] Alexandrova, O., Jagarlamudi, V. K., Hellinger, P., Maksimovic, M., Shprits, Y., & Mangeney, A. (2021). Spectrum of kinetic plasma turbulence at 0.3–0.9 astronomical units from the Sun. Physical Review E, 103(6), 063202.

[2] Arrò, G., Califano, F., & Lapenta, G. (2021). Spectral properties and energy cascade at kinetic scales in collisionless plasma turbulence. arXiv preprint arXiv:2112.12753.

How to cite: Arrò, G., Califano, F., and Lapenta, G.: Spectral features and energy cascade of kinetic scale plasma turbulence, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-4040,, 2022.

EGU22-4447 | Presentations | NP6.3

Electron heating scales in quasi-perpendicular shocks 

Andreas Johlander, Andrew Dimmock, Yuri Khotyaintsev, Daniel Graham, and Ahmad Lalti

Collisionless shock waves are important for particle heating and acceleration in space. Electron heating at shocks is a combination of adiabatic heating due to large-scale electric and magnetic fields and scattering by high-frequency oscillations. Electron heating and scattering at the shock is still poorly understood but the scales at which heating happens can hint to which physical processes are taking place. Here, we study electron heating scales with the Magnetospheric Multiscale (MMS) spacecraft at Earth’s quasi-perpendicular bow shock. We utilize the small tetrahedron formation and rapid plasma measurements of MMS to directly measure the electron temperature gradient inside the shock. From this, we reconstruct the electron temperature profile inside the shock ramps of a number of shock crossings with varying shock parameters. We find that most of the electron temperature increase takes place on a scale of tens of electron inertial lengths. Further, we investigate the electron distribution functions and attempt to disentangle the effects of the large-scale adiabatic heating and scattering by high-frequency waves.

How to cite: Johlander, A., Dimmock, A., Khotyaintsev, Y., Graham, D., and Lalti, A.: Electron heating scales in quasi-perpendicular shocks, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-4447,, 2022.

EGU22-5068 | Presentations | NP6.3

Measurement of short  (Debye length) scale electrostatic waves with MMS EDP instrument:  Problems and possible mitigation 

Ahmad Lalti, Yuri Khotyaintsev, Daniel Graham, Konrad Steinvall, and Andreas Johlander

High frequency electrostatic oscillations are one of the most fundamental players in energy conversion in collisionless plasmas. Whether at collisionless shocks, turbulence energy cascades or reconnection, small scale Debye length processes  are at the heart of irreversible energy exchange between particles and fields. MMS is one of the most advanced still active spacecraft, with high resolution field and particle instruments. The electric field instrument (EDP) on board of MMS is formed of 3 axial double probes positioned in a perpendicular configuration allowing for the measurement of the 3D electric field. In this study we probe the limitations of the EDP instrument in measuring Debye-scale electrostatic oscillations. In particular we show that at such small wavelengths the electric field is attenuated due to the finite probe-to-probe separation. Furthermore, we propose a method to correct for the electric field attenuation based on the single spacecraft interferometry technique which will allow us to properly determine the observed wave modes.

How to cite: Lalti, A., Khotyaintsev, Y., Graham, D., Steinvall, K., and Johlander, A.: Measurement of short  (Debye length) scale electrostatic waves with MMS EDP instrument:  Problems and possible mitigation, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-5068,, 2022.

EGU22-5192 | Presentations | NP6.3

Identifying active magnetic reconnection in simulations and in situ observations of plasma turbulence using magnetic flux transport 

Tak Chu Li, Yi-Hsin Liu, Yi Qi, and Christopher T. Russell

For decades, magnetic reconnection has been suggested to play an important role in the dynamics and energetics of plasma turbulence by spacecraft observations, numerical simulations and theory. Reliable approaches to study reconnection in turbulence are essential to advance this frontier topic of plasma physics. A new method based on magnetic flux transport (MFT) has been recently developed to identify reconnection activity in turbulent plasmas. Applications to gyrokinetic simulations of two- and three-dimensional (2D and 3D) plasma turbulence, and MMS observations of reconnection events in the magnetosphere have demonstrated the capability and accuracy of MFT in identifying active reconnection in turbulence. In 2D, MFT identifies multiple active reconnection X-lines; two of them have developed bi-directional electron and ion outflow jets, observational signatures for reconnection, while one of the X-line does not have bi-directional electron or ion outflow jets, beyond the category of electron-only reconnection recently discovered in the turbulent magnetosheath. In 3D, plentiful reconnection X-lines are identified through MFT, and a new picture of reconnection in turbulence results. In space, MMS observations have provided first evidence for MFT signatures of active reconnection under varying plasma conditions throughout the Earth's magnetosphere. MFT is applicable to in situ measurements by spacecraft missions, including PSP and Solar Orbiter, and laboratory experiments.

How to cite: Li, T. C., Liu, Y.-H., Qi, Y., and Russell, C. T.: Identifying active magnetic reconnection in simulations and in situ observations of plasma turbulence using magnetic flux transport, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-5192,, 2022.

EGU22-5222 | Presentations | NP6.3

Fluctuations of the solar wind ion flux near the Earth bow shock 

Olga Sapunova, Natalia Borodkova, Yuri Yermolaev, and Georgii Zastenker

In our study we analyzing fluctuations of the solar wind ion flux associated with the Earth bow shock using data obtained by the BMSW experiment, installed onboard the SPEKTR-R satellite. The high time resolution of the spectrometer (0.031 s for the plasma flux magnitude and direction and 1.5 s for velocity, temperature, and density) makes it possible to study fine structures in detail.

From 2011 to 2019 SPEKTR-R satellite crossed the Earth bow shock many times. In our work we analyzed more than 200 bow shock crossings including multiple ones. More than half of them had fluctuations near the Earth bow shock front.

It was shown that in 25% of events the frequencies of ion flux fluctuations were in the range of 3-4 Hz. In 5-7% of events the frequencies of ion flux fluctuations lay in the interval of 5-6 Hz. Just few cases had frequencies of ion flux fluctuations equal or more than 7 Hz. In other cases the frequencies of ion flux fluctuations were lower than 3 Hz or no fluctuations were observed at all.

We also observed low-frequencies fluctuations about 0.1 Hz and lower. These fluctuations were also visible by the 1.5 s plasma parameters: protons density and velocity; He++ (alpha particles) density and velocity (including helium abundance).

How to cite: Sapunova, O., Borodkova, N., Yermolaev, Y., and Zastenker, G.: Fluctuations of the solar wind ion flux near the Earth bow shock, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-5222,, 2022.

EGU22-5560 | Presentations | NP6.3

Statistical study of the ripples and reformation in the collisionless shocks using MMS observation 

Ajay Lotekar, Yuri Khotyaintsev, Daniel Graham, Andrew Dimmock, Ahmad Lalti, and Andreas Johlander

Collisionless shocks are ubiquitous throughout the universe in near-Earth and astrophysical plasma environments. The behavior of collisionless shocks in terms of their structure and energy dissipation has been the subject of extensive research over many decades, but many open questions remain. Recent studies have demonstrated that the Earth's bow shock can exhibit ripples that propagate along the shock surface. However, their occurrence, dependence on shock parameters, and their role in shock dynamics is still under investigation. One signature of rippling is the presence of phase space holes in reduced ion distribution (integrated along the tangential plane of the shock). Such ion phase space holes are also observed in association with the shock reformation. It is unclear at what part of the parametric space these ion phase space holes are expected. In this study, we have focused on characterizing ion phase space holes at the Earth’s bow shock using MMS observations. We analyze more than 500 shock crossings observed by the MMS spacecraft and establish a systematic procedure to find the shocks exhibiting phase space holes. We investigate the key shock physical processes responsible for the existence of these phase space holes (e.g. ripples and reformation) and study the association to shock parameters such as Mach number and geometry. We present the first statistical study of this nature, and these results are important to understanding the non-stationary behavior of collisionless shocks. 

How to cite: Lotekar, A., Khotyaintsev, Y., Graham, D., Dimmock, A., Lalti, A., and Johlander, A.: Statistical study of the ripples and reformation in the collisionless shocks using MMS observation, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-5560,, 2022.

EGU22-6084 | Presentations | NP6.3

Suprathermal populations and small scale fluctuations in the solar wind 

Marian Lazar, Rodrigo Lopez, Hamd Shaaban, Stefaan Poedts, Horst Fichtner, and Peter Yoon

In the last decade, studies of solar wind plasma have shown that suprathermal populations (up to a few keV) are closely linked to wave turbulence and fluctuations at small (or kinetic) scales. We aim to identify those types of wave fluctuations observed at these scales scales, for which existing theories predict a major implication in particle acceleration and formation of suprathermal tails in the velocity distributions of plasma particles. On the other hand, it is currently believed that fluctuation power (magnetic, density, velocity) measured at ion scales and lower are generated by the turbulent cascade but also wave instabilities. Therefore, we also intend to discuss a number of recent results describing the kinetic instabilities driven by the anisotropy of velocity distributions (e.g., temperature anisotropy, field-aligned drifts), and how are these instabilities influenced by the suprathermal populations. These results help to understand the energy exchanges between particles and electromagnetic fields, not only in the solar wind but also in the coronal plasma ejections, with consequences for the space weather and terrestrial magnetosphere.

How to cite: Lazar, M., Lopez, R., Shaaban, H., Poedts, S., Fichtner, H., and Yoon, P.: Suprathermal populations and small scale fluctuations in the solar wind, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-6084,, 2022.

EGU22-6378 | Presentations | NP6.3 | Highlight

Kinetic signatures of magnetic reconnection in the global hybrid-Vlasov and local particle-in-cell simulations.  

Ivan Zaitsev, Andrey Divin, Urs Ganse, Yann Pfau-Kempf, Markus Battarbee, Markku Alho, Jonas Suni, Maxime Grandin, Lucile Turc, Giulia Cozzani, Maarja Bussov, Maxime Dubart, Harriet George, Konstantinos Horaites, Konstantinos Papadakis, Talgat Manglayev, Vertti Tarvus, Honyang Zhou, and Minna Palmroth

Magnetic reconnection is the energy converter in space plasma that releases magnetic energy into the kinetic energy of particles. We study the magnetotail reconnection in the first 3D global magnetospheric hybrid-Vlasov simulation performed with Vlasiator code. We also performed a simulation of symmetric magnetic reconnection in particle-in-cell technique with the iPIC3D code to compare ion kinetic signatures of reconnection for both hybrid-Vlasov and fully-kinetic approaches. Despite the relatively coarse spatial resolution in the global 3D hybrid-Vlasov model, we are able to recognize the most distinguished reconnection features: ion demagnetization, non-gyrotropic ion acceleration and energy dissipation. Using the well-known signatures of the different subregions of symmetric magnetic reconnection we are able to identify ion diffusion regions, separatrices and reconnection jet fronts in the global simulation. Guided by the measure of the ion perpendicular slippage, we identify ion diffusion regions where ion non-gyrotropic crescent-type distributions are formed. These distinguishable features are nicely visible in the PIC simulation data as well. Separatrix regions are visible as the layers containing the potential Hall electric field at the boundaries of accelerated outflow. Reconnection jet fronts in the global simulation are highlighted at the positions where the energy dissipation peaks. Three-dimensional effects affecting the extending of the reconnection characteristics in the equatorial plane are discussed.

How to cite: Zaitsev, I., Divin, A., Ganse, U., Pfau-Kempf, Y., Battarbee, M., Alho, M., Suni, J., Grandin, M., Turc, L., Cozzani, G., Bussov, M., Dubart, M., George, H., Horaites, K., Papadakis, K., Manglayev, T., Tarvus, V., Zhou, H., and Palmroth, M.: Kinetic signatures of magnetic reconnection in the global hybrid-Vlasov and local particle-in-cell simulations. , EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-6378,, 2022.

EGU22-6722 | Presentations | NP6.3

Asymmetric magnetic reconnection between two coalescing flux ropes as squeezed by convergent flows 

Qiaowen Luo, Jiansen He, and Jun Cui

We identify two coalescing flux ropes as squeezed by convergent ion flows near the magnetopause from the MMS observations. According to the electron distributions, we find that one flux rope is closer to the magnetosphere, while the other is closer to the magnetosheath. A current sheet with magnetic field reversal is found to sit at the interface between the two colliding flux ropes, and have magnetic reconnection occurring in the ion diffusion region (IDR). Due to the density asymmetry of flux ropes, the embedded magnetic reconnection event with a significant guide field component shows a large asymmetry in energy conversion across the reconnection site. On the side where the flux rope is closer to the magnetosphere with low density, we find that electrons gained energy from electromagnetic fields resulting in a parallel heating effect. In contrast, ions are found to obtain the energy from electromagnetic fields on the other side of the reconnection current sheet, where the flux rope is near the magnetosheath with high density.

How to cite: Luo, Q., He, J., and Cui, J.: Asymmetric magnetic reconnection between two coalescing flux ropes as squeezed by convergent flows, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-6722,, 2022.

EGU22-7069 | Presentations | NP6.3

Interplay between magnetic reconnection and flapping instabilities in the magnetotail: global hybrid-Vlasov simulation of the Earth’s magnetosphere 

Giulia Cozzani, Urs Ganse, Yann Pfau-Kempf, Markku Alho, Jonas Suni, Maxime Grandin, Lucile Turc, Ivan Zaitsev, Maarja Bussov, Maxime Dubart, Harriet George, Konstantinos Horaites, Talgat Manglayev, Konstantinos Papadakis, Vertti Tarvus, Honyang Zhou, and Minna Palmroth

Magnetic reconnection is a fundamental process in plasma and a major cause of energy conversion and transport by means of magnetic field topology reconfiguration. It takes place in thin plasma sheets, where energy is often explosively converted from the magnetic field to plasma heating and particle energization. Magnetic reconnection in Earth’s magnetotail is thought to play a crucial role in geomagnetic storms and substorms, one of the most explosive phenomena in the context of Earth’s magnetosphere. Several other current sheet-related processes, such as the ballooning instability, tearing instability, and a variety of flapping instabilities, can occur in the magnetotail, and the interplay between magnetic reconnection and these current sheet instabilities is largely unexplored. In this study, we investigate the interplay between magnetic reconnection and other instabilities taking place in the magnetotail current sheet, using a hybrid-Vlasov simulation that provides a three-dimensional description of the global coupled solar wind-magnetosphere system down to the ion-kinetic scale. In particular, we identify and characterize the flapping instability that develops in the magnetotail midnight sector and we discuss its dynamics in relation to magnetotail magnetic reconnection.

How to cite: Cozzani, G., Ganse, U., Pfau-Kempf, Y., Alho, M., Suni, J., Grandin, M., Turc, L., Zaitsev, I., Bussov, M., Dubart, M., George, H., Horaites, K., Manglayev, T., Papadakis, K., Tarvus, V., Zhou, H., and Palmroth, M.: Interplay between magnetic reconnection and flapping instabilities in the magnetotail: global hybrid-Vlasov simulation of the Earth’s magnetosphere, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-7069,, 2022.

EGU22-7354 | Presentations | NP6.3

Hybrid-Vlasov simulations of ion velocity distribution functions within Kelvin-Helmholtz vortices 

Vertti Tarvus, Lucile Turc, Hongyang Zhou, Giulia Cozzani, Urs Ganse, Yann Pfau-Kempf, Markku Alho, Markus Battarbee, Maarja Bussov, Maxime Dubart, Harriet George, Maxime Grandin, Konstantinos Horaites, Talgat Manglayev, Konstantinos Papadakis, Jonas Suni, Ivan Zaitsev, and Minna Palmroth

The Kelvin-Helmholtz instability (KHI) is a ubiquitous fluid instability in space plasmas. At the flanks of Earth's magnetopause, the KHI can typically develop during periods of northward interplanetary magnetic field, and it drives the solar wind-magnetosphere mass/energy transfer in the absence of dayside magnetic reconnection. We use local 2D-3V hybrid-Vlasov simulations to study the ion velocity distribution functions (VDFs) associated with the KHI in a magnetopause-like setup. Our results indicate that when the KHI enters the non-linear stage, the ion VDFs in the region perturbed by the instability become increasingly non-Maxwellian. The degree of non-Maxwellianity increases along with the magnitude of the density jump across the KHI boundary. We assess the impact of the non-Maxwellian ion VDFs on the development of the KHI, and compare the simulated VDFs with those observed by the Magnetospheric Multiscale Mission.

How to cite: Tarvus, V., Turc, L., Zhou, H., Cozzani, G., Ganse, U., Pfau-Kempf, Y., Alho, M., Battarbee, M., Bussov, M., Dubart, M., George, H., Grandin, M., Horaites, K., Manglayev, T., Papadakis, K., Suni, J., Zaitsev, I., and Palmroth, M.: Hybrid-Vlasov simulations of ion velocity distribution functions within Kelvin-Helmholtz vortices, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-7354,, 2022.

EGU22-7371 | Presentations | NP6.3

Uniturbulence due to non-linear damping of surface Alfvén waves 

Rajab Ismayilli and Tom Van Doorsselaere

We consider a simple 1-D planar equilibrium model with piece-wise constant density. We completed analytical computations in incompressible MHD. First, we derived mathematical formulas for the wave energy density, the rate of energy dissipation, and the energy cascade damping time. Following that, we developed an analytical model to estimate the damping time for the evolution of uniturbulence in surface Alfvén waves. According to the derived equation, the damping time is inversely proportional to the perpendicular wavenumber and the amplitude of the surface Alfvén waves. Next, we determined the numerical energy dissipation rate using the Fourier transform through numerical simulations. Finally, we approximated the damping time using the fundamental mode of a perpendicular wavenumber.
Consequently, we compared our theoretical model to a series of 3D ideal MHD simulations and observed a remarkable resemblance. The numerical findings demonstrate, in particular, that the damping time is inversely related to the density contrast and amplitude of surface Alfvén waves. Besides, we studied third-order structure-function (Yaglom's law) for Uniturbulence. We compared Yaglom's law (predicted energy dissipation) statistics obtained through simulation with our analytical model. In addition, we estimated the inertial range of the turbulent flow.

How to cite: Ismayilli, R. and Van Doorsselaere, T.: Uniturbulence due to non-linear damping of surface Alfvén waves, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-7371,, 2022.

EGU22-8160 | Presentations | NP6.3

Exploring 2.5D magnetic reconnection due to Rayleigh-Taylor induced turbulence in solar prominences 

Madhurjya Changmai and Rony Keppens

The internal dynamics of solar prominences have been observed for many decades to be highly complex, many of which also indicate the possibility of turbulence. Prominences represent large-scale, dense condensations suspended against gravity at great heights within the solar atmosphere. It is therefore of no surprise that the fundamental process of the Rayleigh-Taylor (RT) instability has been suggested as the potential mechanism for driving the dynamics and turbulence remarked upon within observations. We use the open-source MPI-AMRVAC code to construct an extremely high-resolution, 2.5D fully-resistive magnetohydrodynamic model, and employ it to explore the turbulent nature of RT-induced magnetic reconnection processes within solar prominences. The intermittent events of heating and energy dissipation are caused by magnetic reconnection. Furthermore, the strength of the mean magnetic field directed into the 2D plane, and its alignment with the plane itself, creates a system with varying turbulent behaviour. Based on low plasma beta (magnetic pressure dominant) evolution near the chromosphere and a higher value (plasma pressure dominant) evolution within the corona, the stratified numerical model generates different fluctuation statistics. Hence, we find the turbulent dynamics and prominence reconnection events to differ distinctly from those elsewhere within the solar corona.

How to cite: Changmai, M. and Keppens, R.: Exploring 2.5D magnetic reconnection due to Rayleigh-Taylor induced turbulence in solar prominences, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-8160,, 2022.

EGU22-8705 | Presentations | NP6.3

Fully kinetic simulations of the near-Sun solar wind plasma: turbulence, reconnection, and particle heating 

Luca Franci, Emanuele Papini, Alfredo Micera, Lorenzo Matteini, Julia Stawarz, Giovanni Lapenta, David Burgess, Petr Hellinger, Simone Landi, Andrea Verdini, and Victor Montagud-Camps

We model the development of plasma turbulence in the near-Sun solar wind with high-resolution fully-kinetic particle-in-cell (PIC) simulations, initialised with plasma conditions measured by Parker Solar Probe during its first solar encounter (ion and electron plasma beta ≤ 1 and a large amplitude of the turbulent fluctuations). The power spectra of the plasma and electromagnetic fluctuations are characterized by multiple power-law intervals, with a transition and a considerable steepening in correspondence of the electron scales. In the same range of scales, the kurtosis of the magnetic fluctuations is observed to further increase, hinting at a higher level of intermittency. We observe a number of electron-only reconnection events, which are responsible for an increase of the electron temperature in the direction parallel to the ambient field. The total electron temperature, however, exhibits only a small increase due to the cooling of electrons in the perpendicular direction, leading to a strong temperature anisotropy. We also analyse the power spectra of the different terms of the electric field in the generalised Ohm’s law, their linear and nonlinear components, and their alignment, to get a deeper insight on the nature of the turbulent cascade. Finally, we compare our results with those from hybrid simulations with the same parameters, as well as with spacecraft observations.

How to cite: Franci, L., Papini, E., Micera, A., Matteini, L., Stawarz, J., Lapenta, G., Burgess, D., Hellinger, P., Landi, S., Verdini, A., and Montagud-Camps, V.: Fully kinetic simulations of the near-Sun solar wind plasma: turbulence, reconnection, and particle heating, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-8705,, 2022.

In supercritical shocks a substantial fraction of ions is reflected at the steep shock ramp. The beam of reflected ions carries a considerable amount of energy and momentum. As a consequence, different plasma populations can co-exist within the same foot region, which constitutes a source of micro-instabilities excited by the relative drifts between incoming ions, reflected ions, and electrons across the ambient magnetic field B. With the help of a spectral periodic 2D PIC code, we investigate the resulting micro-turbulence. Three different waves with different frequency/wave number ranges can be excited simultaneously: Bernstein waves and whistler waves near the lower-hybrid frequency as well as the electron cyclotron frequency. The present work is a 2D extension of a previous analysis (Muschietti et Lembege, Ann. Geophys. 2017) and allows to self-consistently include the mutual interaction between the different instabilities/waves which propagate in different directions with respect to Bo and are at different stages of their respective linear/nonlinear phases. In order to clarify their intricate synergies, a new filtering procedure (low or high pass filter of a given wave number range) has been developed. Taking thus advantage of the spectral nature of the code, we can include/exclude at will the impact of a given instability on the other ones. We have performed several times the simulation with exactly the same initial conditions yet with different filtering ranges. The procedure allows us to illuminate the role played by each instability in the scenario when all are included. Recent results will be presented. 

How to cite: Muschietti, L., Lembege, B., and Decyk, V.: How to define the interplay between different instabilities excited within the foot of a supercritical shock : 2D PIC simulations, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-8962,, 2022.

The Jovian magnetosphere is loaded internally with material from the volcanic moon of Io, which is ionised and brought into co-rotation forming the Io plasma torus. Plasma is removed from the torus mainly via ejection as energetic neutrals and by bulk transport into sink regions in the outer magnetosphere.

There are two physical processes that are implicated in the bulk transport process, these are diffusion and the radial-interchange (RI) instability. The latter is analogous to the Rayleigh-Taylor instability, but with centrifugal force replacing gravity. This allows magnetic flux tubes containing hot, tenuous plasma to exchange places with tubes containing cool, dense plasma, moving material from the inner to outer magnetosphere whilst returning magnetic flux to the planet. Observational data does not currently provide strong evidence to favour either process and indeed they may be non-linearly coupled. Furthermore, current state-of-the-art simulations do not permit an understanding of non-linear phases of the instability nor the effect of magnetosphere-ionosphere coupling on small length scales.

In order to examine the bulk transport process we have developed a full hybrid kinetic ion, fluid-electron plasma model in 2.5-dimensions, JERICHO. The technique of hybrid modelling allows for probing of plasma motions from the scales of planetary-radii down to the ion-inertial length, considering constituent ion species kinetically as charged particles and forming the electrons into a single magnetised fluid continuum. This allows for insights into particle motions on spatial scales below the size of the magnetic flux tubes. We are particularly interested in exploring a) bulk transport on spatial scales not currently accessible with other state-of-the-art models; b) the relative contributions from diffusive motions against those from RI instabilities; and c) non-linear effects generated by RI instabilities and the impact of these on plasma transport from the inner to outer magnetosphere. In this presentation we will examine the latest simulation results from JERICHO, initialised with a range of Jovian parameters, examining the evolution of the RI instability on differing spatial and temporal scales.

How to cite: Wiggs, J. and Arridge, C.: Examining Radial-Interchange in the Jovian Magnetosphere using JERICHO: a Kinetic-Ion, Fluid-Electron Hybrid Model, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-9788,, 2022.

EGU22-9885 | Presentations | NP6.3

MMS bow shock crossings database 

Yuri Khotyaintsev, Ahmad Lalti, Andrew P. Dimmock, Andreas Johlander, and Daniel B. Graham

Identifying collisionless shock crossings in data sent from spacecraft has so far been done manually. It is a tedious job that shock physicists have to go through if they want to conduct case studies or perform statistical studies. We use a machine learning approach to automatically identify shock crossings from the Magnetospheric Multiscale (MMS) spacecraft. We compile a database of those crossings including various spacecraft related and shock related parameters for each event. Furthermore, we show that the shocks in the database have properties that are spread out both in real space and parameter space. We also present a possible science application of the database by looking for correlations between ion acceleration efficiency at shocks and different shock parameters such as the shock geometry and the Mach number.

How to cite: Khotyaintsev, Y., Lalti, A., Dimmock, A. P., Johlander, A., and Graham, D. B.: MMS bow shock crossings database, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-9885,, 2022.

EGU22-10339 | Presentations | NP6.3

The combined effect of electron and proton firehose instabilities for the solar wind plasma conditions 

Rodrigo A. López, Alfredo Micera, Marian Lazar, Shaaban M. Shaaban, Stefaan Poedts, and Giovanni Lapenta

In the absence of collision, kinetic instabilities triggered by velocity space anisotropies of plasma particles play an essential role in limiting the deviations from isotropy. For example, in the solar wind, firehose instabilities may inhibit the growth of the temperatures in the direction parallel to the background magnetic field, counterbalancing the effect of the expansion. Electron and proton firehose instabilities can be triggered depending on the plasma parameters and the different branches within (periodic and aperiodic). Despite the significant difference between electron and proton spatial and temporal scales, both modes can work together to alter the dynamic of the plasma.
We use a fully kinetic 2D semi-implicit particle-in-cell simulation, iPic3D, to study the evolution and interplay of firehose instabilities triggered by electrons and protons when both species are anisotropic. The aperiodic electron firehose instability remains largely unaffected by the proton anisotropy and saturates rapidly at low-level fluctuations. On the other hand, the presence of anisotropic electrons has a considerable impact on the proton firehose modes, especially on the aperiodic branch, shifting the onset of the instability and boosting the saturation levels of the fluctuations. Anisotropic electrons contribute to more effective regulation of the proton anisotropy.

How to cite: López, R. A., Micera, A., Lazar, M., Shaaban, S. M., Poedts, S., and Lapenta, G.: The combined effect of electron and proton firehose instabilities for the solar wind plasma conditions, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-10339,, 2022.

EGU22-10688 | Presentations | NP6.3

Multi-Spacecraft Observations of Interplanetary Shocks 

Oksana Kruparova, Vratislav Krupar, and Adam Szabo

Interplanetary (IP) shocks provide us with a unique opportunity to extensively investigate properties of collisionless shocks using in situ measurements under a wide range of upstream conditions. Here we report a case study of several IP shock crossings observed by the Wind, Solar and Heliospheric Observatory (SOHO), Advanced Composition Explorer (ACE), and Deep Space Climate Observatory (DSCOVR) spacecraft. By applying a simple timing method to multipoint measurements, we are able to investigate their characteristic spatiotemporal features. We assume that an IP shock can be represented by a moving plane with a constant velocity, when observed at closely separated points in space and time. We compared IP shock parameters obtained with the timing method with those obtained using the magnetic coplanarity, the mixed mode methods, and Rankine-Hugoniot jump relations.


How to cite: Kruparova, O., Krupar, V., and Szabo, A.: Multi-Spacecraft Observations of Interplanetary Shocks, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-10688,, 2022.

EGU22-11028 | Presentations | NP6.3

Enernization of Alpha Particles in the Solar Wind Magnetic Reconnection 

Die Duan, Jiansen He, Xingyu Zhu, Rui Zhuo, Ziqi Wu, and Liu Yang
The acceleration and heating of solar wind particles by energy dissipation in the magnetic reconnection is an important problem in space physics. Although alpha particles are the second most abundant element in the solar wind, their dynamical behavior in the magnetic reconnection of the solar wind is not well understood. Using the high energy (1500~3000 eV) part of the SWA/PAS instrument on board the Solar Orbiter, we study the kinetic behavior of alpha particles in a magnetic-reconnetion exhaust region within a heliospheric current sheet. In this event, protons and alpha particles have similar bulk velocities. Alpha particles are accelerated and form a jet in the exhaust region. The counter-stream distribution of alpha particles is observed inside the exhuast region, which changes the direction from parallel to perpendicular to the magnetic field direction when the magnetic field is reversed. In addition, a pair of the slow shock/rotational discontinuity is observed in the exhuast region. The exhaust region is heated and bounded by the slow shocks , while the accelerated plasma jet is bounded by the rotational discontinuities.

How to cite: Duan, D., He, J., Zhu, X., Zhuo, R., Wu, Z., and Yang, L.: Enernization of Alpha Particles in the Solar Wind Magnetic Reconnection, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-11028,, 2022.

EGU22-11524 | Presentations | NP6.3

Turbulence-driven magnetic reconnection and the magnetic correlation length in collisionless plasma turbulence 

Julia E. Stawarz, Jonathan P. Eastwood, Tai Phan, Imogen L. Gingell, Prayash S. Pyakurel, Michael A. Shay, Sadie L. Robertson, Christopher T. Russell, and Olivier Le Contel

Observations of Earth’s magnetosheath from the Magnetospheric Multiscale (MMS) mission have provided an unprecedented opportunity to examine the detailed structure of the multitude of thin current sheets that are generated by plasma turbulence, revealing that a novel form of magnetic reconnection, which has come to be known as electron-only reconnection, can occur within magnetosheath turbulence. These electron-only reconnection events occur at thin electron-scale current sheets and have super-Alfvénic electron jets that can approach the electron Alfvén speed; however, they do not appear to have signatures of ion jets. It is thought that electron-only reconnection can occur when the length of the reconnecting current sheets along the outflow direction is short enough that the ions cannot fully couple to the newly reconnected magnetic field lines before they fully relax. In this work, we examine how the correlation length of the magnetic fluctuations in a turbulent plasma, which constrains the length of the current sheets that can be formed by the turbulence, impacts the nature of turbulence-driven magnetic reconnection. Using observations from MMS, we systematically examine 60 intervals of magnetosheath turbulence – identifying 256 small-scale reconnection events, both with and without ion jets. We demonstrate that the properties of the reconnection events transition to become more consistent with electron-only reconnection when the magnetic correlation length of the turbulence is below ~20 ion inertial lengths. We further discuss the implications of the results in the context of other turbulent plasmas by considering observations of turbulent fluctuations in the solar wind.

How to cite: Stawarz, J. E., Eastwood, J. P., Phan, T., Gingell, I. L., Pyakurel, P. S., Shay, M. A., Robertson, S. L., Russell, C. T., and Le Contel, O.: Turbulence-driven magnetic reconnection and the magnetic correlation length in collisionless plasma turbulence, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-11524,, 2022.

EGU22-11660 | Presentations | NP6.3

Magnetosheath jets at the magnetopause: reconnection onset conditions 

Adrian LaMoury, Heli Hietala, Jonathan Eastwood, Laura Vuorinen, and Ferdinand Plaschke

Magnetosheath jets are localised pulses of high dynamic pressure plasma observed in Earth’s magnetosheath. They are believed to form from the interaction between the solar wind and ripples in Earth’s collisionless bow shock, before propagating into the turbulent magnetosheath. Upon impacting the magnetopause, jets can influence magnetospheric dynamics. In particular, previous studies have suggested that, by virtue of their internal magnetic field orientations, jet impacts may be able to trigger local magnetic reconnection at the magnetopause. This is most notable during traditionally unfavourable solar wind conditions, such as intervals of northward interplanetary magnetic field. This idea has been supported by a small number of case studies and simulations. We present a large statistical study into the properties of jets near the magnetopause. We examine the components of the magnetic reconnection onset condition – the competing effects of magnetic shear angle and plasma beta – to determine how jets may affect magnetopause reconnection in a statistical sense. We find that, due to their increased beta, jet plasma is typically not favourable to reconnection, often more so than the non-jet magnetosheath. Most jets do contain some reconnection-favourable plasma, however, suggesting that jets may be able to both trigger and suppress magnetopause reconnection. We complement this with new case studies of jets interacting with the magnetopause.

How to cite: LaMoury, A., Hietala, H., Eastwood, J., Vuorinen, L., and Plaschke, F.: Magnetosheath jets at the magnetopause: reconnection onset conditions, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-11660,, 2022.

EGU22-11807 | Presentations | NP6.3 | Highlight

Agyrotropy patterns in 3D small-scall turbulent reconnection 

Jeffersson Andres Agudelo Rueda, Daniel Verscharen, Robert T. Wicks, Christopher J. Owen, Andrew P. Walsh, and Kai Germaschewski
Turbulence and magnetic reconnection are at the core of the long-standing problem of energy dissipation in collisionless plasmas. More than two decades of research on magnetic reconnection have led us to understand the characteristic plasma flows and particle agyrotropy patterns present in collisionless reconnection events. However, it is still not clear what the agyrotropy patterns associated with reconnection events are that form in a turbulent cascade. In this work, we use an explicit fully kinetic particle-in-cell code to study the plasma particles’ agyrotropy associated with three-dimensional small-scale magnetic reconnection events generated by anisotropic and Alfvénic decaying turbulence. We select one reconnection event involving two reconnecting flux ropes. Although we observe similarities with agyrotropy patterns known from two-dimensional steady-state reconnection events, the agyrotropy patterns in our event are more complex. This has further implications for the energy transfer channels available in three-dimensional turbulent reconnection.

How to cite: Agudelo Rueda, J. A., Verscharen, D., Wicks, R. T., Owen, C. J., Walsh, A. P., and Germaschewski, K.: Agyrotropy patterns in 3D small-scall turbulent reconnection, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-11807,, 2022.

EGU22-11945 | Presentations | NP6.3 | Highlight

Turbulence modified by velocity shear in coronal mass ejection sheaths 

Juska Soljento, Simon Good, Adnane Osmane, and Emilia Kilpua

Fast coronal mass ejections (CMEs) drive shock waves ahead of them. The turbulent sheath region between the shock and the CME itself contains magnetic field and velocity fluctuations on a broad spectrum of frequencies. In this work we aim to characterise the direction and source of solar wind fluctuations at MHD fluid scales in CME-driven sheaths near Earth. One possible source for these fluctuations is velocity shear, which are common occurrences in CME-driven sheaths. Here we first identify velocity shear as it occurs and then relate that to signatures of new fluctuations being created locally in the sheath. Turbulence parameters such as cross helicity, residual energy, Elsasser ratio, and Alfvén ratio are calculated, and they are correlated against large-scale signatures of velocity shear. Findings indicate a clear association between velocity shear and locally generated fluctuations, as well as a balance in the directionality of these new fluctuations, i.e., they tend to propagate equally towards and away from the Sun. In contrast, most solar wind is typically dominated by anti-sunward fluctuations.

How to cite: Soljento, J., Good, S., Osmane, A., and Kilpua, E.: Turbulence modified by velocity shear in coronal mass ejection sheaths, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-11945,, 2022.

The kinetic turbulence generated by accelerated particles in a reconnecting current sheet (RCS) with X- and O-nullpoints is considered. The simulations of magnetic reconnection using particle-in-cell (PIC) approach is carried out in a thin current sheet with 3D magnetic field topology affected by tearing instability that leads to a formation of two large magnetic islands . The model utilises a strong guiding field that leads to separation of the particles of opposite charges, generation of a strong polarisation electric field across the RCS and suppression of kink instability in the ’out-of-plane’ direction. The accelerated particles of the same charge entering an RCS from the opposite edges are shown accelerated to different energies forming the ‘bump-in-tail’ velocity distributions that, in turn, can generates plasma turbulence in different locations. The turbulence produced by either electron or proton beams is identified from the energy spectra of electromagnetic field fluctuations in the phase and frequency domains.

The spectral index of the power spectrum In a wavenumber space of the turbulent magnetic field near the ion inertial length approaches -2.7. The collective turbulence power spectra are consistent with the high-frequency fluctuations of perpendicular electric field, or upper hybrid waves, to occur in a vicinity of X-nullpoints, with the Langmuir waves  generated by accelerated electrons which can be converted to  Bernstein waves when electron beams become moving across the magnetic field lines. The frequency spectra of high and low-frequency waves are explored in the kinetic turbulence in parallel and perpendicular directions to the local magnetic field showing noticeable lower hybrid turbulence. The implication of finding for observations is also discussed.

How to cite: Zharkova, V. and Xia, Q.: Kinetic turbulence generated by accelerated particles in a reconnecting current sheet with magnetic islands, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-12159,, 2022.

EGU22-12840 | Presentations | NP6.3

Studying multi-beam ion temperature inside a collisionless reconnection plasmoid by means of Gaussian Mixture Model 

Igor Paramonik, Andrey Divin, Ivan Zaitsev, and Vladimir Semenov

Being ubiquitous energy converter is space plasmas, magnetic reconnection releases stored magnetic energy into kinetic energy of particles. Magnetic reconnection involves several particle acceleration mechanisms which form beams directed parallel to the magnetic field. It was recently demonstrated analytically that in the presence of complicated velocity space structures, the definition of higher moments (like thermal pressure) should be extended to cover such multibeam distributions. In practice, the number of beams at each spatial point of interest is not know a priori. With the aim to automatically reveal the information about the beams generated in the reconnection process, we applied an unsupervised machine learning algorithm (Gaussian Mixture Model, GMM) to the 2.5D Particle-in-Cell simulations of collisionless magnetic reconnection. We studied the ion distributions inside a plasmoid and found that the multibeam ion temperature  within the reconnected outflow deviates significantly from the standard ion temperature (calculated as the 2nd moment of the ion distribution function). In particular, the regions of the strong parallel heating contain in fact relatively cold counterstreaming beams and the overestimation of parallel temperature in this case could be as high as 10. In the current study, we make an attempt to figure out how long the multi-beam regime exists without significant thermalization inside a plasmoid formed by two adjacent X-lines.

How to cite: Paramonik, I., Divin, A., Zaitsev, I., and Semenov, V.: Studying multi-beam ion temperature inside a collisionless reconnection plasmoid by means of Gaussian Mixture Model, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-12840,, 2022.

Cosmogenic radionuclides concentrations are predominantly determined by the solar activity and space weather around the Earth, forming an important source of cosmic-origin background radiation in the terrestrial environment. The highest values of such radiation are observed during the solar minima because the penetrability of the Earth’s magnetosphere is greatest at that time. Beryllium 7Be binds to aerosols and is transported within a few years to the Earth’s surface. Its concentrations are higher during the spring and summer months when the stratospheric 7Be penetrates the troposphere as a result of the exchange of air masses between the troposphere and stratosphere. We compare periods of strong solar and geomagnetic storms with periods of very low solar activity in the longitudinal view during the years 1986 – 2020.

For a better understanding of the process dynamics, in our work we investigate the coupling of concentrations of the cosmogenic radionuclide 7Be (time series of activity concentration of 7Be in aerosols) to space weather parameters around the Earth (Kp planetary index, disturbance storm time Dst, proton density, proton flux), proxy parameters of the solar activity (intensity of solar radio flux, relative sunspot number), stratospheric dynamics parameters (temperature, zonal component of wind, O3), and aggregates of strong atmospheric frontal transition. The beryllium radionuclide 7Be concentration was evaluated by the corresponding activity in aerosols on a weekly basis at the National Radiation Protection Institute Monitoring Section in Prague.

We also perform the case study of cosmogenic radionuclide 7Be concentrations during the period of strong solar and geomagnetic storm in November 2021 with the ERA5 reanalysis data, and Aeolus satellite lidar wind measurements.

How to cite: Podolská, K., Kozubek, M., Hýža, M., and Šindelářová, T.: The effect of space weather, proxy parameters of solar activity, and stratospheric phenomena on the concentration of cosmogenic radionuclide 7Be (in the Czech Republic), EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-4659,, 2022.

EGU22-5539 | Presentations | GI6.4

Measurements of cosmic rays by a mini neutron monitor aboard the German research vessel Polarstern. 

Bernd Heber, Sasa Banjac, Sönke Burmeister, Martin Zoska, Hanna Giese, Konstantin Herbst, Lisa Romaneehsen, Carolin Schwerdt, Dutoit Stauss, Carsten Wallmann, Adrian Vogt, and Michael Walter

Galactic cosmic rays (GCRs) consist of energetic electrons and nuclei which are a direct sample of material from far beyond the solar system. Measurements by various particle detectors have shown that the intensity varies on different timescales, caused by the Sun’s activity and geomagnetic variation. Interplanetary disturbances cause space weather effects which warrant a more detailed study. Many studies on GCR intensity decreases is based on the analysis of ground-based neutron monitors and muon telescopes. Their measurements depend on the geomagnetic position, and the processes in the Earth's atmosphere. In order to get a better understanding of the geomagnetic filter over the solar cycle, the Christian-Albrechts-Universität zu Kiel, DESY Zeuthen, and the North-West University in Potchefstroom, South Africa agreed on a regular monitoring of the GCR intensity as a function of latitude, by installing a portable device aboard the German research vessel Polarstern in 2012. The vessel is ideally suited for this research campaign because it covers extensive geomagnetic latitudes (i.e. goes from the Arctic to the Antarctic) at least once per year. Here we present the measurements for different latitude surveys including the periods of solar maximum in 2014 and solar minimum in 2019. 

The Kiel team received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No 870405. The team would like to thank the crew of the Polarstern and the AWI for supporting our research campaign.

How to cite: Heber, B., Banjac, S., Burmeister, S., Zoska, M., Giese, H., Herbst, K., Romaneehsen, L., Schwerdt, C., Stauss, D., Wallmann, C., Vogt, A., and Walter, M.: Measurements of cosmic rays by a mini neutron monitor aboard the German research vessel Polarstern., EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-5539,, 2022.

EGU22-6238 | Presentations | GI6.4 | Highlight

Challenges and solutions for cosmic-ray neutron sensing in heterogeneous soil moisture situations related to irrigation practices 

Cosimo Brogi, Heye Reemt Bogena, Markus Köhli, Harrie-Jan Hendricks Franssen, Olga Dombrowski, Vassilios Pisinaras, Anna Chatzi, Kostantinos Babakos, Jannis Jakobi, Patrizia Ney, and Andreas Panagopoulos

Water availability is a key challenge in agriculture, especially given the expected increase of droughts related to climate change. Soil moisture (SM) sensors can be used to collect information on water availability in a reliable and accurate way. However, due to their very small measuring volume, the installation of multiple sensors is required. In addition, in-situ sensors may need to be removed during field management and connecting cables are often damaged by rodents and other wilderness animals. Hence, the demand for SM sensors that do not have such limitations will increase in the upcoming years. A promising non-invasive technique to monitor SM is cosmic-ray neutron sensing (CRNS), which is based on the negative correlation between fast neutrons originating from cosmic radiation and SM content. With its large measuring footprint of ~130-210m, CRNS can efficiently cover the field-scale. However, heterogeneous agricultural management (e.g., irrigation) can lead to abrupt SM differences, which pose a challenge for the analysis of CRNS data. Here, we investigate the effects of small-scale soil moisture patterns on the CRNS signal by using both modelling approaches and field studies. The neutron transport model URANOS was used to simulate the neutron signal of a CRNS station located in irrigated plots of different sizes (from 1 to 8 ha) with different soil moisture (from 5 and 50 Vol.%) inside and outside such a plot. A total of 400 different scenarios were simulated and the response functions of multiple detector types were further considered. In addition, two CRNS with Gadolinium shielding were installed in two irrigated apple orchards of ~1.2 ha located in the Pinios Hydrologic Observatory (Greece) in the context of the H2020 ATLAS project. Reference soil moisture was determined using 25 SoilNet stations, each with 6 SM sensors installed in pairs at 5, 20 and 50 cm depth and water potential sensors at 20 cm depth. The orchards were also equipped with two Atmos41 climate stations and eight water meters for irrigation monitoring. The CRNS were calibrated using either soil samples or the SM measured by the SoilNet network. In the URANOS simulations, the percentage of neutrons detected by the CRNS that are representative of an irrigated plot varied between 45 and 90% and was strongly influenced by both the dimension and SM of the irrigated plot. As expected, the CRNS footprint decreased considerably with increasing SM but did not appear to be influenced by the plot dimension. SM variation within the irrigated plot strongly affected the neutron energy at detection, which was not the case for SM variations outside the plot. The instrumented fields corroborated the URANOS findings and the performance of the local CRNS was dependent on a) the timing and intensity of irrigation and precipitation, b) the CRNS calibration strategy, and c) the management of the surrounding fields. These results provide novel and meaningful information on the impact of horizontal SM patterns on CRNS measurements, which will help to make CRNS more useful in irrigated agriculture.

How to cite: Brogi, C., Bogena, H. R., Köhli, M., Hendricks Franssen, H.-J., Dombrowski, O., Pisinaras, V., Chatzi, A., Babakos, K., Jakobi, J., Ney, P., and Panagopoulos, A.: Challenges and solutions for cosmic-ray neutron sensing in heterogeneous soil moisture situations related to irrigation practices, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-6238,, 2022.

EGU22-6430 | Presentations | GI6.4

Utilizing Cosmic Ray data as input for neutron-based soil moisture measurement 

Hanna Giese, Bernd Heber, Konstantin Herbst, and Martin Schrön

Neutrons on Earth interact with the soil and are substantially moderated by hydrogen atoms. Since the reflected neutron flux is a function of the soil water content, cosmic-ray neutron measurements above the ground can be used to estimate the average field soil moisture. Thus, if the local incoming neutron flux and the abundance of nearby hydrogen pools are known, the reflected neutron flux could be modeled and compared to observed detector count rates. However, the incoming neutrons are secondaries produced by interacting energetic Galactic Cosmic Rays (GCRs) in the atmosphere. The total neutron flux on the ground depends on the solar modulation-dependent GCR flux, the geomagnetic position, and the altitude within the atmosphere. So far, measurements of either the Jungfraujoch neutron monitor (NM) or a NM of similar cutoff rigidity have been used and altered to estimate the neutron flux at the position of each neutron detector. In this contribution we present a new method based on the Dorman function to directly compute the local neutron flux using remote neutron monitor data.

We received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No. 870405

How to cite: Giese, H., Heber, B., Herbst, K., and Schrön, M.: Utilizing Cosmic Ray data as input for neutron-based soil moisture measurement, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-6430,, 2022.

EGU22-8872 | Presentations | GI6.4

Relationship between the time series of cosmic ray data and aerosol optical properties: (Case study: southern Italy, 2016-2020) 

Faezeh Karimian Sarakhs, Fabio Madonna, Marco Rosoldi, and Salvatore De Pasquale


High energy Cosmic Ray (CR) particles are capable of ionizing the Earth’s atmosphere, which leads to changes in the atmospheric physical and chemical properties. One of the most important effects of interactions between the CR particles and atmospheric molecules is the formation of aerosol and its subsequent condensation nuclei processes. These interactions are known with considerable uncertainty yet and may translate into even bigger uncertainties in future climate predictions. Laser Detection and Ranging (LIDAR) is currently the best suited technology to retrieve aerosol optical and microphysical properties is also used for the atmosphere correction of high energy cosmic ray observatory data. LIDAR measurements are available from single stations or from networks at continental scale like the European Aerosol Research LIdar NETwork (EARLINET). Sun photometer data are the most suitable complement to LIDAR measurements for the study of aerosol properties due to the extensive coverage of their measurements available through the AErosol RObotic NETwork (AERONET) network. The purpose of this study is to find the correlation between the aerosol properties and the CR data. The aerosol properties retrieved from two databases for the period of 2016-2020: I) the multi-wavelength LIDAR system Potenza EArlinet Raman Lidar (PEARL) which operates at the CNR-IMAA (Tito Scalo (Italy) and contributes to the EAELINET); and II) the AERONET sun photometer data from the stations located at Southern Italy i.e. Potenza (40.60° N, 15.72° E, 820m), Naples (40.83° N, 14.30° E, 50 m) and Lecce (40.33° N, 18.11° E, 30m). whereas, the CR data made available in Italy from the Extreme Energy Events project ( Air mass back-trajectories were used to confirm the observed aerosol types and support the correlation study. Our study showed promising results in understanding the relationship between cosmic ray and aerosol properties.

Keywords: Cosmic Ray, Aerosol, Lidar, Sun Photometer, Back-trajectory

How to cite: Karimian Sarakhs, F., Madonna, F., Rosoldi, M., and De Pasquale, S.: Relationship between the time series of cosmic ray data and aerosol optical properties: (Case study: southern Italy, 2016-2020), EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-8872,, 2022.

EGU22-9264 | Presentations | GI6.4

Smart Scintillating Neutron Detectors for Soil Moisture Monitoring 

Patrick Stowell and the COSMIC-SWAMP and STFC Food Network+ Collaborations

Cosmic ray neutron sensing has been shown to be a powerful method for continuously monitoring soil moisture over large areas. This technique relies on the detection of albedo cosmic ray neutrons coming from from the soil to infer the local hydrogen content. Cosmic ray neutron sensing is well-suited for hydrological monitoring in the field sizes typically seen on smallholder farms. The ongoing development of new lower-cost neutron detector instrumentation and processing tools will help to further support the adoption of this novel technique within the agricultural industry.

In this presentation I will discuss recent efforts at Durham University (UK) to develop low-cost cosmic ray neutron detectors to support soil moisture monitoring in the agriculture sector. These systems rely on lithium and boron-based scintillator foils for thermal neutron detection. Recent pilot studies in collaboration with the COSMOS-UK network have shown that the detected neutron rate in these sensors correlates well with results obtained from traditional gaseous systems. Work is now underway to improve the robustness of these scintillator systems for use in agricultural and civil engineering applications. 

In addition, I will present a new international research network, COSMIC-SWAMP, which is looking at the integration of cosmic ray neutron sensors with managed irrigation sites in Brazil. By combining low-cost neutron probes with a smart water management platform (SWAMP), this research network is looking at using cosmic ray neutrons to perform data-driven irrigation control over large areas. The instrumentation being considered for COSMIC-SWAMP will be presented before discussing the future plans for the network.

How to cite: Stowell, P. and the COSMIC-SWAMP and STFC Food Network+ Collaborations: Smart Scintillating Neutron Detectors for Soil Moisture Monitoring, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-9264,, 2022.

EGU22-9438 | Presentations | GI6.4

Geomagnetic field shielding over the past 100 000 years 

Monika Korte, Jiawei Gao, and Sanja Panovska

The geomagnetic field prevents energetic particles, such as galactic cosmic rays, from directly interacting with the Earth's atmosphere. The geomagnetic field is not static but constantly changing, and over the last 100,000 years several geomagnetic excursions occurred. During geomagnetic field excursions, the field strength is significantly decreased and the field morphology is controlled by non-dipole components, and more cosmic ray particles can access the Earth's atmosphere. Paleomagnetic field models provide a global view of the long-term geomagnetic field evolution, however, with individual spatial and temporal resolution. Here, we reconstruct the geomagnetic shielding effect over the last 100,000 years by calculating the geomagnetic field cutoff rigidity using four global paleomagnetic field models, i.e., GGF100k, GGFSS70, LSMOD.2, and CALS10k.2. We find that the non-dipole components of the geomagnetic field are not negligible for estimating the long-term geomagnetic shielding effect, in particular during excursions. Our results indicate that cosmic ray flux, impact area, and cosmic ray radiation intensity increase strongly during the excursions. Our results provide the possibility to accurately estimate the cosmogenic isotope production rate and cosmic radiation dose rate covering the last 100,000 years.

How to cite: Korte, M., Gao, J., and Panovska, S.: Geomagnetic field shielding over the past 100 000 years, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-9438,, 2022.

EGU22-11230 | Presentations | GI6.4

Space or climate? Disentangling cosmogenic and climatic drivers of present-day tritium (3H) in global precipitation 

Stefan Terzer-Wassmuth, Luis J. Araguas-Araguas, Lorenzo Copia, and Jodie A. Miller

The generation of cosmogenic tritium (3H) through spallation of 14N in the upper atmosphere and a its decay (half-life of 12.32 y) are the two main processes resulting in the global steady-state inventory of tritium in the hydrosphere of approximately 2.95 kg. Various mechanisms of scavenging of stratospheric 3H into the troposphere, such as stratosphere-troposphere transports (STTs) during the so-called “spring leak”, or the tropospheric distribution by means of the Brewer-Dobson circulation, have been described to explain the observed spatial and seasonal distribution of present-day tritium levels in global precipitation. Following thermonuclear weapons testing prior to the Preliminary Test Ban Treaty in 1963, the natural 3H input signal was overlaid by the so-called “bomb peak”. This characteristic tritium pulse has been used for decades in nuclear and hydrological sciences, with 3H values in Vienna, the reference northern hemisphere station of the IAEA-WMO Global Network of Isotopes in Precipitation (GNIP), peaking in 1963 at approximately 400 Bq L-1. Since the year 2000, this 3H pulse has dissipated in the northern hemisphere, and 3H levels at the Vienna monitoring site have reached their natural background value of ca. 1.2 Bq L-1.

The present-day steady state of natural 3H levels in precipitation allow to research their inter-annual variability as driven by cosmogenic input, with particular emphasis on neutron flux intensity governed by the 11-year sunspot cycles. With almost two full solar cycle’s worth of observed 3H data in Vienna’s precipitation and other GNIP stations in the northern hemisphere, we discuss the impact of the neutron flux (as exemplified by the Oulu Neutron Monitor) in modulating the inter-annual variability. Our findings showed that while 52% of the interannual variability was explained by changes in the cosmogenic flux, an additional 31% of the variability resulted from the seasonal distribution of the amount of precipitation, a finding prominent in the previous solar cycle valley, particularly in the year 2015, that coincided with abnormally high winter precipitation.

While the regular oscillations of the neutron flux seem to constitute the main driver of the observed interannual changes of 3H contents in precipitation, atmospheric circulation processes were of varying importance in 15 GNIP stations. In spite of the relative data paucity (i.e. absence of sufficiently long records at even spatial distribution), we hypothesize that changes in precipitation seasonality, due to climate change impacts on global or regional atmospheric circulation patterns, may drive fluctuations in the natural steady-stage 3H levels in precipitation used to investigate atmospheric and hydrological processes. Hence, we stress the importance of spatially and temporally adequate observational baselines on a global level.

How to cite: Terzer-Wassmuth, S., Araguas-Araguas, L. J., Copia, L., and Miller, J. A.: Space or climate? Disentangling cosmogenic and climatic drivers of present-day tritium (3H) in global precipitation, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-11230,, 2022.

EGU22-12334 | Presentations | GI6.4

Rail-based cosmic ray neutron sensing (CRNS): pushing the boundaries towards expanding footprints and temporal resolutions 

Daniel Altdorff, Sascha Oswald, Steffen Zacharias, Carmen Zengerle, Hannes Mollenhauer, Peter Dietrich, Sabine Attinger, and Martin Schrön

Cosmic ray neutron sensing (CRNS) has become an established method for deriving the soil water content (SWC), based on the inverse relationship of neutron counting and the SWC of the surrounding area. The provided footprint, lateral up to 200m and vertical of several decimeter, qualifies CRNS to bridge the information gap between classical hydrogeophysical approaches and remote sensing. While stationary CRNS offers continuous long-term SWC measurements at high temporal resolution, the covered area remains fixed and predefined. Car-borne CRNS roving on the other hand, allows to expand the mapped area. However, the method requires active operation and is limited to snap shot information only. As an alternative, the operation of a permanent mobile CRNS platform on trains promises to combine the advantages from stationary and car-borne CRNS measurements, as recently suggested by Schrön et al. (2021), while also its technical implementation, data processing and interpretation raises new challenges and complexity.

In this study we introduce a fully automatic CRNS railway system, installed in a conventional locomotive of a freight train, as first and novel of its kind. Results of the first phase of operation will be presented. The measurements along an experimental rail track were supported by local SWC measurements, gravimetric and dielectric records (Mobile Wireless Ad-hoc Sensor Network), at three areas along the railway, and by a newly installed weather station. Additionally, car-borne CRNS data were recorded on two days close to the railway track.

Preliminary results of data collected between September and December 2021 showed very stable spatial pattern in relation to the segments crossed by the train, which have been confirmed by the car-borne dataset. Temporal variations within hours were also evident as direct or indirect response to local rain and snow events.  Based on the first results, we are confident, that rail-based CRNS offers the chance to play a prominent role in addressing soil hydrology at landscape scale in the future.

Schrön, M., Oswald, S. E., Zacharias, S., Kasner, M., Dietrich, P., & Attinger, S. (2021). Neutrons on rails: Transregional monitoring of soil moisture and snow water equivalent. Geophysical Research Letters, 48, e2021GL093924


How to cite: Altdorff, D., Oswald, S., Zacharias, S., Zengerle, C., Mollenhauer, H., Dietrich, P., Attinger, S., and Schrön, M.: Rail-based cosmic ray neutron sensing (CRNS): pushing the boundaries towards expanding footprints and temporal resolutions, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-12334,, 2022.

Cosmogenic isotopes are mostly produced in the stratosphere and troposphere, and the corresponding fractions depend on solar activity and tropopause altitude. Solar-cycle variability of cosmogenic isotope production is the strongest at high latitudes due to the lack of geomagnetic shielding. However, the exact zonal distribution of the production in troposphere and stratosphere regions, that is needed for the precise modelling of their transport and deposition, is not clear. In this work, we provide numerical estimates of cosmogenic isotopes production in the atmosphere for different conditions. Using the SOCOL-AER2-BE Chemical Climatic model (CCM), we present simulations of the production of cosmogenic isotopes ($^{14}$C, $^{36}$Cl, $^{10}$Be, and $^{7}$Be) and provide zonal distributions (tropical, subtropical, and polar regions) in the stratosphere and troposphere. The model is driven by four solar activity scenarios: 1) solar minimum year with solar modulation function - phi = 400 MeV and 2) solar maxima year with phi = 1100MeV. In these cases, the production is modulated by Galactic Cosmic Rays (GCR). Two other scenarios are 3) ground-level enhancement (GLE) event number 5 with hard spectrum on February 23, 1956 and 4) GLE event number 24 with soft spectrum on August 04, 1972. The production rates were calculated using a combination of the SOCOL and CRAC models.

How to cite: Golubenko, K.: Zonal distribution of cosmogenic isotopes in stratosphere and troposphere via CCM SOCOL, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-12618,, 2022.

EGU22-12812 | Presentations | GI6.4

Signal contribution of remote areas to cosmic-ray neutron sensors based on distance and sensitivity 

Martin Schrön, Markus Köhli, and Steffen Zacharias

Cosmic-Ray Neutron Sensing (CRNS) is an established measurement technique for water content in soils and snow. The high integration depth and the large measurement footprint is an important advantage compared to conventional point-scale sensors. However, the radial-symmetrical footprint definition based on the 86% quantile of detected neutrons is often not helpful to explain the influence of certain areas in complex fields. Many natural sites are highly heterogeneous and thus knowledge of the contribution of distant areas to the measurement signal would be very useful, e.g. to support calibration sampling, sensor location design, data interpretation, and uncertainty assessment. Here, CRNS calibration and validation remains a challenge, since the influence of the different fields and structures to the signal is usually not known.

In this presentation, we proposes a generalized analytical procedure to estimate the contribution of patches or fields in the footprint of a cosmic-ray neutron detector to its signal using the radial intensity functions. The proposed method could greatly support calibration sampling, sensor location design, and uncertainty assessment, e.g. in complex or vegetated terrain, without the need of computationally expensive neutron modeling. Furthermore, a new concept for a more practical definition of the sensor footprint is proposed, which represents the maximal distance to a field such that its soil moisture change is still sensible in terms of measurement precision. 

How to cite: Schrön, M., Köhli, M., and Zacharias, S.: Signal contribution of remote areas to cosmic-ray neutron sensors based on distance and sensitivity, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-12812,, 2022.

EGU22-486 | Presentations | ITS2.1/PS1.2

Enhancing planetary imagery with the holistic attention network algorithm 

Denis Maxheimer, Ioannis Markonis, Masner Jan, Curin Vojtech, Pavlik Jan, and Solomonidou Anezina

The recent developments in computer vision research in the field of Single Image Super Resolution (SISR)

can help improve the satellite imagery data quality and, thus, find application in planetary exploration.

The aim of this study is to enhance planetary surface imagery, in planetary bodies that there are

available data but in a low resolution. Here, we have applied the holistic attention network (HAN)

algorithm to a set of images of Saturn’s moon Titan from the Titan Radar Mapper instrument in its

Synthetic Aperture Radar (SAR) mode, which was on board the Cassini spacecraft. HAN can find

correlations among hierarchical layers, channels of each layer, and all positions of each channel, which

can be interpreted as an application and intersection of previously known models. The algorithm used

in our case-study was trained on 5000 grayscale images from HydroSHED Earth surface imagery dataset

resampled over different resolutions. Our experimental setup was to generate High Resolution (HR)

imagery from eight times lower resolution (x8 scale). We followed the standard workflow for this

purpose, which is to first train the network enhancing x2 scale to HR, then x4 scale to x2 scale, and

finally x8 scale to x4 scale, using subsequently the results of the previous training. The promising results

open a path for further applications of the trained model to improve the imagery data quality, and aid

in the detection and analysis of planetary surface features.

How to cite: Maxheimer, D., Markonis, I., Jan, M., Vojtech, C., Jan, P., and Anezina, S.: Enhancing planetary imagery with the holistic attention network algorithm, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-486,, 2022.

EGU22-692 | Presentations | ITS2.1/PS1.2

Autonomous lineament detection in Galileo images of Europa 

Caroline Haslebacher and Nicolas Thomas

Lineaments are prominent features on the surface of Jupiter's moon Europa. Analysing these linear features thoroughly leads to insights on their formation mechanisms and the interactions between the subsurface ocean and the surface. The orientation and position of lineaments is also important for determining the stress field on Europa. The Europa Clipper mission is planned to launch in 2024 and will fly by Europa more than 40 times. In the light of this, an autonomous lineament detection and segmentation tool would prove useful for processing the vast amount of expected images efficiently and would help to identify processes affecting the ice sheet. 

We have trained a convolutional neural network to detect, classify and segment lineaments in images of Europa returned by the Galileo mission. The Galileo images that make up the training set are segmented manually, following a dedicated guideline. For better performance, we make use of synthetically generated data to pre-train the network. The current status of the work will be described.

How to cite: Haslebacher, C. and Thomas, N.: Autonomous lineament detection in Galileo images of Europa, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-692,, 2022.

EGU22-1014 | Presentations | ITS2.1/PS1.2

Automatic detection of the electron density from the WHISPER instrument onboard CLUSTER II 

Emmanuel De Leon, Nicolas Gilet, Xavier Vallières, Luca Bucciantini, Pierre Henri, and Jean-Louis Rauch

The Waves of HIgh frequency and Sounder for Probing Electron density by Relaxation
(WHISPER) instrument, is part of the Wave Experiment Consortium (WEC) of the CLUSTER II
mission. The instrument consists of a receiver, a transmitter, and a wave spectrum
analyzer. It delivers active (when in sounding mode) and natural electric field spectra. The
characteristic signature of waves indicates the nature of the ambient plasma regime and, combined
with the spacecraft position, reveals the different magnetosphere boundaries and regions. The
thermal electron density can be deduced from the characteristics of natural waves in natural mode
and from the resonances triggered in sounding mode, giving access to a key parameter of scientific
interest and major driver for the calibration of particles instrument.
Until recently, the electron density derivation required a manual time/frequency domain
initialization of the search algorithms, based upon visual inspection of WHISPER active and natural
spectrograms and other datasets from different instruments onboard CLUSTER.
To automate this process, knowledge of the region (plasma regime) is highly desirable. A Multi-
Layer Perceptron model has been implemented for this purpose. For each detected region, a GRU,
recurrent network model combined with an ad-hoc algorithm is then used to determine the electron
density from WHISPER active spectra. These models have been trained using the electron density
previously derived from various semi-automatic algorithms and manually validated, resulting in an
accuracy up to 98% in some plasma regions. A production pipeline based on these models has been
implemented to routinely derive electron density, reducing human intervention up to 10 times. Work
is currently ongoing to create some models to process natural measurements where the data volume
is much higher and the validation process more complex. These models of electron density
automated determination will be useful for future other space missions.

How to cite: De Leon, E., Gilet, N., Vallières, X., Bucciantini, L., Henri, P., and Rauch, J.-L.: Automatic detection of the electron density from the WHISPER instrument onboard CLUSTER II, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-1014,, 2022.

EGU22-2765 | Presentations | ITS2.1/PS1.2

Extrapolation of CRISM based spectral feature maps using CaSSIS four-band images with machine learning techniques 

Michael Fernandes, Nicolas Thomas, Benedikt Elser, Angelo Pio Rossi, Alexander Pletl, and Gabriele Cremonese

Spectroscopy provides important information on the surface composition of Mars. Spectral data can support studies such as the evaluation of potential (manned) landing sites as well as supporting determination of past surface processes. The CRISM instrument on NASA’s Mars Reconnaissance Orbiter is a high spectral resolution visible infrared mapping spectrometer currently in orbit around Mars. It records 2D spatially resolved spectra over a wavelength range of 362 nm to 3920 nm. At present data collected covers less than 2% of the planet. Lifetime issues with the cryo-coolers prevents limits further data acquisition in the infrared band. In order to extend areal coverage for spectroscopic analysis in regions of major importance to the history of liquid water on Mars (e.g. Valles Marineris, Noachis Terra), we investigate whether data from other instruments can be fused to extrapolate spectral features in CRISMto these non-spectral imaged areas. The present work will use data from the CaSSIS instrument which is a high spatial resolution colour and stereo imager onboard the European Space Agency’s ExoMars Trace Gas Orbiter (TGO). CaSSIS returns images at 4.5 m/px from the nominal 400 km altitude orbit in four colours. Its filters were selected to provide mineral diagnostics in the visible wavelength range (400 – 1100 nm). It has so far imaged around 2% of the planet with an estimated overlap of ≲0.01% of CRISM data. This study introduces a two-step pixel based reconstruction approach using CaSSIS four band images. In the first step advanced unsupervised techniques are applied on CRISM hyperspectral datacubes to reduce dimensionality and establish clusters of spectral features. Given that these clusters contain reasonable information about the surface composition, in a second step, it is feasible to map CaSSIS four band images to the spectral clusters by training a machine learning classifier (for the cluster labels) using only CaSSIS datasets. In this way the system can extrapolate spectral features to areas unmapped by CRISM. To assess the performance of this proposed methodology we analyzed actual and artificially generated CaSSIS images and benchmarked results against traditional correlation based methods. Qualitative and quantitative analyses indicate that by this novel procedure spectral features of in non-spectral imaged areas can be predicted to an extent that can be evaluated quantitatively, especially in highly feature-rich landscapes.

How to cite: Fernandes, M., Thomas, N., Elser, B., Rossi, A. P., Pletl, A., and Cremonese, G.: Extrapolation of CRISM based spectral feature maps using CaSSIS four-band images with machine learning techniques, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-2765,, 2022.

EGU22-2994 | Presentations | ITS2.1/PS1.2

Interpretable Solar Flare Prediction with Deep Learning 

Robert Jarolim, Astrid Veronig, Tatiana Podladchikova, Julia Thalmann, Dominik Narnhofer, Markus Hofinger, and Thomas Pock

Solar flares and coronal mass ejections (CMEs) are the main drivers for severe space weather disturbances on Earth and other planets. While the geo-effects of CMEs give us a lead time of about 1 to 4 days, the effects of flares and flare-accelerated solar energetic particles (SEPs) are very immediate, 8 minutes for the enhanced radiation and as short as about 20 minutes for the highest energy SEPs arriving at Earth. Thus, predictions of solar flare occurrence at least several hours ahead are of high importance for the mitigation of severe space weather effects.

Observations and simulations of solar flares suggest that the structure and evolution of the active region’s magnetic field is a key component for energetic eruptions. The recent advances in deep learning provide tools to directly learn complex relations from multi-dimensional data. Here, we present a novel deep learning method for short-term solar flare prediction. The algorithm is based on the HMI photospheric line-of-sight magnetic field and its temporal evolution together with the coronal evolution as observed by multi-wavelengths EUV filtergrams from the AIA instrument onboard the Solar Dynamics Observatory. We train a neural network to independently identify features in the imaging data based on the dynamic evolution of the coronal structure and the photospheric magnetic field evolution, which may hint at flare occurrence in the near future.

We show that our method  can reliably predict flares six hours ahead, with 73% correct flaring predictions (89% when considering only M- and X-class flares), and 83% correct quiet active region predictions.

In order to overcome the “black box problem” of machine-learning algorithms, and thus to allow for physical interpretation of the network findings, we employ a spatio-temporal attention mechanism. This allows us to extract the emphasized regions, which reveal the neural network interpretation of the flare onset conditions. Our comparison shows that predicted precursors are associated with the position of flare occurrence, respond to dynamic changes, and align with characteristics within the active region.

How to cite: Jarolim, R., Veronig, A., Podladchikova, T., Thalmann, J., Narnhofer, D., Hofinger, M., and Pock, T.: Interpretable Solar Flare Prediction with Deep Learning, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-2994,, 2022.

EGU22-5721 | Presentations | ITS2.1/PS1.2

Magnetopause and bow shock models with machine learning 

Ambre Ghisalberti, Nicolas Aunai, and Bayane Michotte de Welle

The magnetopause (MP) and the bow shock (BS) are the two boundaries bounding the magnetosheath, the region between the magnetosphere and the solar wind. Their position and shape depend on the upstream solar wind and interplanetary magnetic field conditions.

Predicting their shape and position is the starting point of many subsequent studies of processes controlling the coupling between the Earth’s magnetosphere and its interplanetary environment. We now have at our disposal an important amount of data from a multitude of spacecraft missions allowing for good spatial coverage, as well as algorithms based on statistical learning to automatically detect the two boundaries. From the data of 9 satellites over 20 years, we identified around 19000 crossings of the BS and 36000 crossings of the MP. They were used, together with their associated upstream conditions, to train a regression model to predict the shape and position of the boundaries. 

Preliminary results indicate that the obtained models outperform analytical models without making simplifying assumptions on the geometry and the dependency over control parameters.

How to cite: Ghisalberti, A., Aunai, N., and Michotte de Welle, B.: Magnetopause and bow shock models with machine learning, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-5721,, 2022.

EGU22-5739 | Presentations | ITS2.1/PS1.2

Deep learning for surrogate modeling of two-dimensional mantle convection 

Siddhant Agarwal, Nicola Tosi, Pan Kessel, Doris Breuer, and Grégoire Montavon

Mantle convection plays a fundamental role in the long-term thermal evolution of terrestrial planets like Earth, Mars, Mercury and Venus. The buoyancy-driven creeping flow of silicate rocks in the mantle is modeled as a highly viscous fluid over geological time scales and quantified using partial differential equations (PDEs) for conservation of mass, momentum and energy. Yet, key parameters and initial conditions to these PDEs are poorly constrained and often require a large sampling of the parameter space to find constraints from observational data. Since it is not computationally feasible to solve hundreds of thousands of forward models in 2D or 3D, some alternatives have been proposed. 

The traditional alternative to high-fidelity simulations has been to use 1D models based on scaling laws. While computationally efficient, these are limited in the amount of physics they can model (e.g., depth-dependent material properties) and predict only mean quantities such as the mean mantle temperature. Hence, there has been a growing interest in machine learning techniques to come up with more advanced surrogate models. For example, Agarwal et al. (2020) used feedforward neural networks (FNNs) to reliably predict the evolution of entire 1D laterally averaged temperature profile in time from five parameters: reference viscosity, enrichment factor for the crust in heat producing elements, initial mantle temperature, activation energy and activation volume of the diffusion creep. 

We extend that study to predict the full 2D temperature field, which contains more information in the form of convection structures such as hot plumes and cold downwellings. This is achieved by training deep learning algorithms on a data set of 10,525 2D simulations of the thermal evolution of the mantle of a Mars-like planet. First, we use convolutional autoencoders to compress the size of each temperature field by a factor of 142. Second,  we compare the use of two algorithms for predicting the compressed (latent) temperature fields: FNNs and long-short-term memory networks (LSTMs).  On the one hand, the FNN predictions are slightly more accurate with respect to unseen simulations (99.30%  vs. 99.22% for the LSTM). On the other hand, Proper orthogonal decomposition (POD) of the LSTM and FNN predictions shows that despite a lower mean relative accuracy, LSTMs capture the flow dynamics better than FNNs. The POD coefficients from FNN predictions sum up to 96.51% relative to the coefficients of the original simulations, while for LSTMs this metric increases to 97.66%. 

We conclude the talk by stating some strengths and weaknesses of this approach, as well as highlighting some ongoing research in the broader field of fluid dynamics that could help increase the accuracy and efficiency of such parameterized surrogate models.

How to cite: Agarwal, S., Tosi, N., Kessel, P., Breuer, D., and Montavon, G.: Deep learning for surrogate modeling of two-dimensional mantle convection, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-5739,, 2022.

EGU22-6371 | Presentations | ITS2.1/PS1.2

STIX solar flare image reconstruction and classification using machine learning 

Hualin Xiao, Säm Krucker, Daniel Ryan, Andrea Battaglia, Erica Lastufka, Etesi László, Ewan Dickson, and Wen Wang

The Spectrometer Telescope for Imaging X-rays (STIX) is an instrument onboard Solar Orbiter. It measures X-rays emitted during solar flares in the energy range from 4 to 150 keV and takes X-ray images by using an indirect imaging technique, based on the Moiré effect. STIX instrument
consists of 32 pairs of tungsten grids and 32 pixelated CdTe detector units. Flare Images can be reconstructed on the ground using algorithms such as back-projection, forward-fit, and maximum-entropy after full pixel data are downloaded. Here we report a new image reconstruction and
classification model based on machine learning. Results will be discussed and compared with those from the traditional algorithms.

How to cite: Xiao, H., Krucker, S., Ryan, D., Battaglia, A., Lastufka, E., László, E., Dickson, E., and Wang, W.: STIX solar flare image reconstruction and classification using machine learning, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-6371,, 2022.

EGU22-8940 | Presentations | ITS2.1/PS1.2

Mars events polyphonic detection, segmentation and classification with a hybrid recurrent scattering neural network using InSight mission data 

Salma Barkaoui, Angel Bueno Rodriguez, Philippe Lognonné, Maarten De Hoop, Grégory Sainton, Mathieu Plasman, and Taichi kawamura

Since deployed on the Martian surface, the seismometer SEIS (Seismic Experiment for Interior Structure) and the APSS (Auxiliary Payload Sensors Suite) of the InSight (Interior Exploration using Seismic Investigations, Geodesy and Heat Transport) mission have been recorded the daily Martian respectively ground acceleration and pressure. These data are essential to investigate the geophysical and atmospheric features of the red planet. So far, the InSight team were able to detect multiple Martian events. We distinguish two types: the artificial events like the lander modes or the micro-tilts known as glitches or the natural events like the pressure drops which are important to estimate the Martian subsurface and the seismic events used to study the interior structure of Mars. Despite the data complexity, the InSight team was able to catalog these events (Clinton et al 2020 for the seismic event catalog, Banfield et al., 2018, 2020 for the pressure drops catalog and Scholz et al. (2020) for the glitches catalog). However, despite all this effort, we are still in front of multiple challenges. In fact,  the seismic events' detection is limited  due to the SEIS sensitivity, which is the origin of a high noise level that may contaminate the seismic events. Thus, we can miss some of them, especially in the noisy period. Besides, their detection is very challenging and require multiple preprocessing task which is time-consuming. For the pressure drops, the detection method used in Banfield et al.  2020 is limited by a threshold equal to 0.3 Pa. Thus, the rest of pressure drops are not included. Plus, due to lack of energy, the pressure sensor was off for several days. As a result, many pressure drops were missed. As a result, being able to detect them directly on the SEIS data which are, in contrast,  provided continuously, is very important.

In this regard, the aim of this study is to overcome these challenges and thus improve the Martian events detection and provide an updated catalog automatically. For that, we were inspired of one of the main technics used today in data processing and analysis in a complete automatic way: it is the Machine Learning and particularly in our case is the Deep Learning. The architecture used for that is the “Hybrid Recurrent Scattering Neural Network” (Bueno et al 2021)  adapted for Mars

How to cite: Barkaoui, S., Bueno Rodriguez, A., Lognonné, P., De Hoop, M., Sainton, G., Plasman, M., and kawamura, T.: Mars events polyphonic detection, segmentation and classification with a hybrid recurrent scattering neural network using InSight mission data, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-8940,, 2022.

EGU22-9077 | Presentations | ITS2.1/PS1.2

Automatic Detection of Interplanetary Coronal Mass Ejections 

Hannah Ruedisser, Andreas Windisch, Ute V. Amerstorfer, Tanja Amerstorfer, Christian Möstl, Martin A. Reiss, and Rachel L. Bailey

Interplanetary coronal mass ejections (ICMEs) are one of the main drivers for space weather disturbances. In the past,
different machine learning approaches have been used to automatically detect events in existing time series resulting from
solar wind in situ data. However, classification, early detection and ultimately forecasting still remain challenges when facing
the large amount of data from different instruments. We propose a pipeline using a Network similar to the ResUNet++ (Jha et al. (2019)), for the automatic detection of ICMEs. Comparing it to an existing method, we find that while achieving similar results, our model outperforms the baseline regarding GPU usage, training time and robustness to missing features, thus making it more usable for other datasets.
The method has been tested on in situ data from WIND. Additionally, it produced reasonable results on STEREO A and STEREO B datasets
with less input parameters. The relatively fast training allows straightforward tuning of hyperparameters and could therefore easily be used to detect other structures and phenomena in solar wind data, such as corotating interaction regions.

How to cite: Ruedisser, H., Windisch, A., Amerstorfer, U. V., Amerstorfer, T., Möstl, C., Reiss, M. A., and Bailey, R. L.: Automatic Detection of Interplanetary Coronal Mass Ejections, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-9077,, 2022.

EGU22-9621 | Presentations | ITS2.1/PS1.2

Machine Learning Techniques for Automated ULF Wave Recognition in Swarm Time Series 

Georgios Balasis, Alexandra Antonopoulou, Constantinos Papadimitriou, Adamantia Zoe Boutsi, Omiros Giannakis, and Ioannis A. Daglis

Machine learning (ML) techniques have been successfully introduced in the fields of Space Physics and Space Weather, yielding highly promising results in modeling and predicting many disparate aspects of the geospace. Magnetospheric ultra-low frequency (ULF) waves play a key role in the dynamics of the near-Earth electromagnetic environment and, therefore, their importance in Space Weather studies is indisputable. Magnetic field measurements from recent multi-satellite missions are currently advancing our knowledge on the physics of ULF waves. In particular, Swarm satellites have contributed to the expansion of data availability in the topside ionosphere, stimulating much recent progress in this area. Coupled with the new successful developments in artificial intelligence, we are now able to use more robust approaches for automated ULF wave identification and classification. Here, we present results employing various neural networks (NNs) methods (e.g. Fuzzy Artificial Neural Networks, Convolutional Neural Networks) in order to detect ULF waves in the time series of low-Earth orbit (LEO) satellites. The outputs of the methods are compared against other ML classifiers (e.g. k-Nearest Neighbors (kNN), Support Vector Machines (SVM)), showing a clear dominance of the NNs in successfully classifying wave events.

How to cite: Balasis, G., Antonopoulou, A., Papadimitriou, C., Boutsi, A. Z., Giannakis, O., and Daglis, I. A.: Machine Learning Techniques for Automated ULF Wave Recognition in Swarm Time Series, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-9621,, 2022.

The solar wind and its variability is well understood at Earth. However, at distances larger than 1AU the is less clear, mostly due to the lack of in-situ measurements. In this study we use transfer learning principles to infer solar wind conditions at Mars in periods where no measurements are available, with the aim of better illuminating the interaction between the partially magnetised Martian plasma environment and the upstream solar wind. Initially, a convolutional neural network (CNN) model for forecasting measurements of the interplanetary magnetic field, solar wind velocity, density and dynamic pressure is trained on terrestrial solar wind data. Afterwards, knowledge from this model is incorporated into a secondary CNN model which is used for predicting solar wind conditions upstream of Mars up to 5 hours in the future. We present the results of this study as well as the opportunities to expand this method for use at other planets.

How to cite: Durward, S.: Forecasting solar wind conditions at Mars using transfer learning, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-10105,, 2022.

EGU22-11501 | Presentations | ITS2.1/PS1.2

Automatic detection of solar magnetic tornadoes based on computer vision methods. 

Dmitrii Vorobev, Mark Blumenau, Mikhail Fridman, Olga Khabarova, and Vladimir Obridko

We propose a new method for automatic detection of solar magnetic tornadoes based on computer vision methods. Magnetic tornadoes are magneto-plasma structures with a swirling magnetic field in the solar corona, and there is also evidence for the rotation of plasma in them. A theoretical description and numerical modeling of these objects are very difficult due to the three-dimensionality of the structures and peculiarities of their spatial and temporal dynamics [Wedemeyer-Böhm et al, 2012, Nature]. Typical sizes of magnetic tornadoes vary from 102 km up to 106 km, and their lifetime is from several minutes to many hours. So far, quite a few works are devoted to their study, and there are no accepted algorithms for detecting solar magnetic tornadoes by machine methods. An insufficient number of identified structures is one of many problems that do not allow studying physics of magnetic tornadoes and the processes associated with them. In particular, the filamentous rotating structures are well delectable only at the limb, while one can only make suppositions about their presence at the solar disk.
Our method is based on analyzing SDO/AIA images at wavelengths 171 Å, 193 Å, 211 Å and 304 Å, to which several different algorithms are applied, namely, the convolution with filters, convolutional neural network, and gradient boosting. The new technique is a combination of several approaches (transfer learning & stacking) that are widely used in various fields of data analysis. Such an approach allows detecting the structures in a short time with sufficient accuracy. As test objects, we used magnetic tornadoes previously described in the literature [e.g., Wedemeyer et al 2013, ApJ; Mghebrishvili et al. 2015 ApJ]. Our method made it possible to detect those structures, as well as to reveal previously unknown magnetic tornadoes.

How to cite: Vorobev, D., Blumenau, M., Fridman, M., Khabarova, O., and Obridko, V.: Automatic detection of solar magnetic tornadoes based on computer vision methods., EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-11501,, 2022.

EGU22-12480 | Presentations | ITS2.1/PS1.2

A versatile exploration method for simulated data based on Self Organizing Maps 

Maria Elena Innocenti, Sophia Köhne, Simon Hornisch, Rainer Grauer, Jorge Amaya, Jimmy Raeder, Banafsheh Ferdousi, James "Andy" Edmond, and Giovanni Lapenta

The large amount of data produced by measurements and simulations of space plasmas has made it fertile ground for the application of classification methods, that can support the scientist in preliminary data analysis. Among the different classification methods available, Self Organizing Maps, SOMs [Kohonen, 1982] offer the distinct advantage of producing an ordered, lower-dimensional representation of the input data that preserves their topographical relations. The 2D map obtained after training can then be explored to gather knowledge on the data it represents. The distance between nodes reflects the distance between the input data: one can then further cluster the map nodes to identify large scale regions in the data where plasma properties are expected to be similar.

In this work, we train SOMs using data from different simulations of different aspects of the heliospheric environment: a global magnetospheric simulation done with the OpenGGCM-CTIM-RCM code, a Particle In Cell simulation of plasmoid instability done with the semi-implicit code ECSIM, a fully kinetic simulation of single X point reconnection done with the Vlasov code implemented in MuPhy2.

We examine the SOM feature maps, unified distance matrix and SOM node weights to unlock information on the input data. We then classify the nodes of the different SOMs into a lower and automatically selected number of clusters, and we obtain, in all three cases, clusters that map well to our a priori knowledge on the three systems. Results for the magnetospheric simulations are described in Innocenti et al, 2021. 

This classification strategy then emerges as a useful, relatively cheap and versatile technique for the analysis of simulation, and possibly observational, plasma physics data.

Innocenti, M. E., Amaya, J., Raeder, J., Dupuis, R., Ferdousi, B., & Lapenta, G. (2021). Unsupervised classification of simulated magnetospheric regions. Annales Geophysicae Discussions, 1-28.

How to cite: Innocenti, M. E., Köhne, S., Hornisch, S., Grauer, R., Amaya, J., Raeder, J., Ferdousi, B., Edmond, J. "., and Lapenta, G.: A versatile exploration method for simulated data based on Self Organizing Maps, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-12480,, 2022.

EGU22-12830 | Presentations | ITS2.1/PS1.2

Re-implementing and Extending the NURD Algorithm to the Full Duration of the Van Allen Probes Mission 

Matyas Szabo-Roberts, Karolina Kume, Artem Smirnov, Irina Zhelavskaya, and Yuri Shprits

Generating reliable databases of electron density measurements over a wide range of geomagnetic conditions is essential for improving empirical models of electron density. The Neural-network-based Upper hybrid Resonance Determination (NURD) algorithm has been developed for automated extraction of electron density from Van Allen Probes electric field measurements, and has been shown to be in good agreement with existing semi-automated methods and empirical models. The extracted electron density data has since then been used to develop the PINE (Plasma density in the Inner magnetosphere Neural network-based Empirical) model, an empirical model for reconstructing the global dynamics of the cold plasma density distribution based only on solar wind data and geomagnetic indices.
In this study we re-implement the NURD algorithm in both Python and Matlab, and compare the performance of these implementations to each other and previous NURD results. We take advantage of a labeled training data set now being available for the full duration of the Van Allen Probes mission to train the network and generate an electron density data set for a significantly longer time period. We perform detailed comparisons between this output, electron density produced from Van Allen Probes electric field measurements using the AURA semi-automated algorithm, and electron density obtained from existing empirical models. We also present preliminary results from the PINE plasmasphere model trained on this extended NURD electron density data set.

How to cite: Szabo-Roberts, M., Kume, K., Smirnov, A., Zhelavskaya, I., and Shprits, Y.: Re-implementing and Extending the NURD Algorithm to the Full Duration of the Van Allen Probes Mission, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-12830,, 2022.

EGU22-565 | Presentations | ST1.5

Linking In-situ Magnetic and Density Structures in the Low Latitude Slow Solar Wind to Solar Origins 

Thomas Woolley, Lorenzo Matteini, Timothy S. Horbury, Stuart D. Bale, Ronan Laker, Lloyd D. Woodham, and Julia E. Stawarz

To date, Parker Solar Probe has completed ten solar encounters and measured a wealth of in-situ data down to heliocentric distances of ~13 solar radii. This data provides a novel opportunity to investigate the near-Sun environment and understand the young slow solar wind. Typically, the slow solar wind observed in the inner heliosphere is split into an Alfvenic and a non-Alfvenic component. The Alfvenic slow wind is thought to originate from overexpanded coronal hole field lines, whereas the non-Alfvenic slow wind could originate from active regions, transient events, or reconnection at the tips of helmet streamers. In this work, we find structures associated with non-Alfvenic slow wind in the low latitude wind measured by Parker Solar Probe. We identify at least two distinct types of structure using magnetic field magnitude, electron pitch angle distributions, and electron number density. After statistically analysing these structures, with a focus on their plasma properties, shape, and location with respect to the heliospheric current sheet, we link them to solar origins. We find structures that are consistent with the plasma blobs seen previously in remote sensing observations.

How to cite: Woolley, T., Matteini, L., Horbury, T. S., Bale, S. D., Laker, R., Woodham, L. D., and Stawarz, J. E.: Linking In-situ Magnetic and Density Structures in the Low Latitude Slow Solar Wind to Solar Origins, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-565,, 2022.

EGU22-1353 | Presentations | ST1.5

Suprathermal Ion Observations Associated with the Heliospheric Current Sheet Crossings during Parker Solar Probe Encounters 7, 8, and 9 

Mihir Desai and the the Parker Solar Probe ISOIS, FIELDS & SWEAP Teams

We report observations of <100 keV/nucleon suprathermal (ST) H, He, O, and Fe ions in association with three separate crossings of the heliospheric current sheet that occurred near perhelia during PSP encounters 7, 8, and 9. In particular, we compare and contrast the ST ion time-intensity profiles, velocity dispersion, pitch-angle distributions, spectral forms, and maximum energies during the three HCS crossings. We find that these unique ST observations are remarkably different in each case, with those during E07 posing the most serious challenges for existing models of ST ion production in the inner heliosphere. In contrast, the ISOIS observations during E08 appear to be consistent with a scenario in which ST ions escape out of the reconnection exhausts into the separatrix layers after getting accelerated up to ~50-100 keV/nucleon by HCS-associated magnetic reconnection-driven processes. Finally, ST ions during the E09 HCS crossing have properties that are somewhat similar to those seen during both E07 and E08 crossings, with ion intensities being higher outside the exhausts and the separatrices, but significant intensity increases are also observed inside the reconnection exhausts. We discuss these new observations in terms of local versus remote acceleration sources as well as in terms of expectations of existing ST ion production and propagation, including reconnection-driven and diffusive acceleration in the inner heliosphere.

How to cite: Desai, M. and the the Parker Solar Probe ISOIS, FIELDS & SWEAP Teams: Suprathermal Ion Observations Associated with the Heliospheric Current Sheet Crossings during Parker Solar Probe Encounters 7, 8, and 9, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-1353,, 2022.

EGU22-1357 | Presentations | ST1.5

Evolution of power spectral density of magnetic field fluctuations in the inner heliosphere 

Jana Safrankova, Zdenek Nemecek, Frantisek Nemec, Tereza Durovcova, Luca Franci, and Alexander Pitna

The solar wind is a unique laboratory to study the turbulent processes occurring in a collisionless plasma with high Reynolds numbers. The paper analyzes power spectra of magnetic field fluctuations that are computed in the frequency range around the break between inertial and kinetic scales. We use observations made during first nine Parker Solar Probe encounters and compare them with observations of the spacecraft moving in other distances from the Sun (e.g., Solar Orbiter) and closer to 1 AU (Wind). A preliminary analysis of magnetic field fluctuations based on PSP and Wind measurements from the MHD to kinetic scales has shown that a relative level of compressive fluctuations increases until 0.25 AU and remains constant till 1 AU whereas a relative level of perpendicular fluctuations does not change with the distance from the Sun. We can conclude that the B slope is controlled by different process(es) close to the Sun against 1 AU and that, in spite of expectations, the critical distance for turbulence evolution is as large as 0.25 AU. We also discuss the role of important physical parameters (e.g., ion beta, temperature anisotropy, collisional age, magnetic field fluctuation amplitude) determining the properties of the turbulent cascade in different heliospheric locations.

How to cite: Safrankova, J., Nemecek, Z., Nemec, F., Durovcova, T., Franci, L., and Pitna, A.: Evolution of power spectral density of magnetic field fluctuations in the inner heliosphere, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-1357,, 2022.

EGU22-1809 | Presentations | ST1.5

Correlated Magnetic Field Dropouts and Plasma Wave Activity at Switchback Boundaries as Observed by Parker Solar Probe 

Anthony Rasca, William Farrell, Robert MacDowall, Stuart Bale, and Justin Kasper

The first solar encounters by the Parker Solar Probe revealed the magnetic field to be dominated by short field reversals in the radial direction referred to as "switchbacks."  While radial velocity and proton temperature were shown to increase inside the switchbacks, B exhibits very brief dropouts only at the switchback boundaries.  Brief intensifications in spectral density measurements near the electron plasma frequency, fpe, have also been observed at these boundaries, indicating the presence of plasma waves triggered by electron beams.  We perform a correlative study using observations from the Parker FIELDS Radio Frequency Spectrometer (RFS) and Fluxgate Magnetometer (MAG) to compare occurrences of spectral density intensifications at the electron plasma frequency (fpe intensifications) and B dropouts at switchback boundaries during Parker's first and second solar encounters.  We find that only a small fraction of minor B dropouts are associated with fpe intensifications.  This fraction increases with B dropout size until all dropouts are associated with fpe intensifications.  This suggests that in the presence of strong B dropouts, electron currents that create the perturbation in B along the boundaries are also stimulating plasma waves such as Langmuir waves.

How to cite: Rasca, A., Farrell, W., MacDowall, R., Bale, S., and Kasper, J.: Correlated Magnetic Field Dropouts and Plasma Wave Activity at Switchback Boundaries as Observed by Parker Solar Probe, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-1809,, 2022.

EGU22-2273 | Presentations | ST1.5

Sheath characteristics of interplanetary coronal mass ejections derived from Helios and PSP data 

Manuela Temmer and Volker Bothmer

Helios 1 and 2 data, covering the distance range from 0.3-1au, have been analysed to derive the characteristics of various substructures of interplanetary coronal mass ejections (ICMEs). We have investigated a data sample of 40 events observed by the Helios 1/2 spacecraft during the time period 1974-1981 with respect to the characteristics of different ICME features, such as sheath regions, leading edges and the magnetic ejecta (ME) themselves. For comparison and to investigate events at distances even closer to the Sun, we add a sample of 5 ICMEs observed with Parker Solar Probe during 2018-2021. We study the sheath density variations over distance and relate those to the ambient solar wind speed. The results show that the sheath region is moderately anti-correlated with the solar wind speed ahead of the disturbance. We further find that the sheath density becomes dominant over the ME density beyond about 0.2au and that its spatial extent constantly increases with distance. The results are important for better understanding the CME mass evolution due to sheath enlargements. Based on these analyses we derive an empirical relation between the sheath density and the local solar wind plasma speed upstream of the ICME shock. The empirical results can be used to model the sheath structure and help improve our understanding about CME propagation in the inner heliosphere.

How to cite: Temmer, M. and Bothmer, V.: Sheath characteristics of interplanetary coronal mass ejections derived from Helios and PSP data, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-2273,, 2022.

EGU22-2638 | Presentations | ST1.5

Alfvénicity of velocity and magnetic field Increments Observed by Parker Solar Probe 

Panisara Thepthong, Peera Pongkitiwanichakul, and David Ruffolo

Alfvénicity is a well-known property of the solar wind. The magnetic and velocity fluctuations are highly correlated over much of the region between the Sun and the Earth. The Parker Solar Probe (PSP) spacecraft enables us to probe Alfvénicity closer to the Sun than before. One previous finding is that the Alfvén ratio increases as the scales become smaller, as shown in some works by using Fourier analysis. This work measures Alfvénicity from PSP observations using 2nd-order structure functions and other quantities derived from magnetic and velocity increments as functions of time lag. We introduce a method to subtract noise from the velocity structure function. We also provide the relation between a time lag in our work and the frequency that contributes most to the 2nd-order structure function for such a time lag. This relation allows a direct comparison between functions of time lag and Fourier spectra. This research has been supported in part by grant RGNS 63-045 from Office of the Permanent Secretary, Ministry of Higher Education, Science, Research and Innovation, and by grant RTA6280002 from Thailand Science Research and Innovation.

How to cite: Thepthong, P., Pongkitiwanichakul, P., and Ruffolo, D.: Alfvénicity of velocity and magnetic field Increments Observed by Parker Solar Probe, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-2638,, 2022.

The high Reynolds number solar wind flow provides a natural laboratory for the study of turbulence in-situ. Parker Solar Probe has to date executed nine sampling distances between 0.2 AU to 1 AU, providing an opportunity to study how turbulence evolves in the expanding solar wind. We focus on data from the PSP/FIELDS[1] and the PSP/SWEAP[2] experiments which provide magnetic field and plasma observations respectively at sub-second cadence. We have identified multiple intervals of uniform solar wind turbulence, selected to exclude coherent structures such as pressure pulses and current sheets, and in which the primary proton population velocity varies by less than 20% of its mean value. We focus on events of multiple-hour duration, which span the spectral scales from the approximately 1/f range at low frequencies, through the magnetohydrodynamic (MHD) inertial range of turbulence and into the kinetic range, below the ion gyrofrequency. We perform a Haar wavelet decomposition[3] which provides accurate estimations of the exponents of these power-law ranges of the spectra, and of higher-order moments. This allows us to study how the spectral exponents may vary with distance from the sun and with solar wind conditions such as the plasma beta. We perform this analysis for both the magnetic field components and magnitude, which track Alfvenic and compressive turbulent fluctuations, respectively. At 1 AU, compressive fluctuations are known to exhibit scaling properties distinct from that of the individual magnetic field components.[4] Here we will investigate this behaviour at different distances from the Sun, plasma beta, and proton density.

We acknowledge the NASA Parker Solar Probe Mission and the SWEAP team led by J. Kasper and the FIELDS team led by S. D. Bale for use of data.

[1] Bale, S.D., Goetz, K., Harvey, P.R. et al. The FIELDS Instrument Suite for Solar Probe Plus.Space Sci Rev 204, 49–82 (2016).

[2] Kasper, J.C., Abiad, R., Austin, G. et al. Solar Wind Electrons Alphas and Protons (SWEAP) Investigation: Design of the Solar Wind and Coronal Plasma Instrument Suite for Solar Probe Plus. Space Sci Rev 204, 131–186 (2016).

[3] Kiyani, K.H., Chapman, S.C., Sahraoui, F., Hnat, B., Fauvarque, O., Khotyaintsev, Y.V.: Enhanced Magnetic Compressibility and Isotropic Scale Invariance at Sub-ion Larmor Scales in Solar Wind Turbulence. The Astrophysical Journal 763(1), 10 (2012).

[4] Hnat, B., Chapman, S., Gogoberidze, G., Wicks, R.: Scale-free texture of the fast solar wind. Physical review. E, Statistical, nonlinear, and soft matter physics 84, 065401 (2011).

How to cite: Wang, X., Chapman, S., Dendy, R., and Hnat, B.: Wavelet analysis of scaling in plasma fluctuations in the magnetohydrodynamic range of solar wind turbulence seen by Parker Solar Probe, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-2945,, 2022.

EGU22-3195 | Presentations | ST1.5

The emission of near-fce harmonic waves at small radial distances from the Sun 

Sabrina F. Tigik, Andris Vaivads, and David M. Malaspina

Near-fce harmonic waves are prevalent in the high-frequency electric field spectrum during Parker Solar Probe's (PSP) close encounters with the Sun. These waves are electrostatic and tend to occur in regions of a relatively stable magnetic field with low broadband magnetic fluctuation levels. We show that the emissions of near-fce harmonic waves are strongly connected to the magnetic field direction. We express the magnetic field direction in terms of spherical angles, where θB is the elevation angle and φB is the azimuthal angle. Then, we show that near-fce harmonics emissions occur when the magnetic field points in a narrow angular range, bounded by 80° ≤ θB ≤ 100° and 10° ≤ φB ≤ 30°, in most of the cases. We also show that the influence of magnetic field direction on near-fce harmonic waves goes down to the shortest time scales the FIELDS instrument can access. These results suggest that cross-scale interaction might play an essential role in the dynamics of the near-fce harmonics measured by PSP at small radial distances from the Sun. It may also provide important clues about the origin of these waves and their role at the early stages of solar wind evolution.

How to cite: F. Tigik, S., Vaivads, A., and M. Malaspina, D.: The emission of near-fce harmonic waves at small radial distances from the Sun, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-3195,, 2022.

EGU22-3859 | Presentations | ST1.5

Polarization of Langmuir waves observed by RPW-TDS instrument on Solar Orbiter during Type III radio bursts 

Jan Soucek, Sampath Bandara, David Pisa, Ondrej Santolik, Milan Maksimovic, Thomas Chust, Yuri Khotyaintsev, Antonio Vecchio, Matthieu Kretzschmar, Robert Wimmer-Schweingruber, Lars Berger, Javier Rodriquez-Pacheco, and Raul Gomez-Herrero

We investigate the polarization of Langmuir waves observed by the Time Domain Sampler (TDS) module of the Radio and Plasma Waves instrument on Solar Orbiter during several extensive Type III burst events. During its two-year-long cruise phase, Solar orbiter often crossed the source region of the Type III radio emission and observed the Langmuir waves generated by solar energetic electrons. The waves are known to exhibit complex modulation and often non-trivial elliptical polarization which sometimes rapidly changes on the timescales of tens of milliseconds. We show that the observed waveforms are typically composed of multiple sub-packets with a relatively short coherence length. We investigate the correlation between the polarization of the waves, simultaneously observed energetic electrons beams and other plasma properties.

How to cite: Soucek, J., Bandara, S., Pisa, D., Santolik, O., Maksimovic, M., Chust, T., Khotyaintsev, Y., Vecchio, A., Kretzschmar, M., Wimmer-Schweingruber, R., Berger, L., Rodriquez-Pacheco, J., and Gomez-Herrero, R.: Polarization of Langmuir waves observed by RPW-TDS instrument on Solar Orbiter during Type III radio bursts, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-3859,, 2022.

EGU22-3917 | Presentations | ST1.5

On the Deflections of Switchbacks 

Ronan Laker, Tim Horbury, Lorenzo Matteini, Thomas Woolley, Julia Stawarz, and Stuart Bale

Following their presence during Parker Solar Probe’s (PSP) first encounter, switchbacks have become an active area of research with several proposed mechanisms for their formation. Many of these theories have testable predictions, although it is not trivial to compare simulation results with in-situ data from PSP. For example, there is some debate regarding the deflection direction of switchbacks, with some theories predicting a consistent magnetic deflection in the +T direction in the RTN coordinate system. Such arguments are largely focussed on the first two PSP encounters, as these are the most studied encounters in the literature. We examine the deflection direction of switchbacks for the first eight PSP encounters, with the aim to clear up any ambiguity regarding this property of switchbacks. Much like the earlier results of Horbury et al. 2020 (during the first encounter) we find that switchbacks tend to deflect in the same direction for hours at a time. Although there is some consistency in deflection direction within an individual encounter, crucially we find that there is no preferred deflection direction across all the encounters. We speculate about the cause of these results and what implications they may have for switchback formation theories.

How to cite: Laker, R., Horbury, T., Matteini, L., Woolley, T., Stawarz, J., and Bale, S.: On the Deflections of Switchbacks, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-3917,, 2022.

EGU22-4130 | Presentations | ST1.5

Observations of the Time Domain Sampler receiver from the Radio and Plasma Wave instrument during the Solar Orbiter Earth flyby 

David Pisa, Jan Soucek, Ondrej Santolik, Miroslav Hanzelka, Milan Maksimovic, Antonio Vecchio, Yuri Khotyaintsev, Thomas Chust, Matthieu Kretzschmar, Lorenzo Matteini, and Timothy Horbury

On November 27, 2021, Solar Orbiter completed its only flyby of Earth on its way to the following Sun’s encounter in March 2022. Although this fast flyby was performed primarily to decrease the spacecraft’s velocity and change orbit to get closer to the Sun, the Radio and Plasma Wave (RPW) instrument had the opportunity to perform high cadence measurements in the Earth’s magnetosphere. We review the main observation of the Time Domain Sampler (TDS) receiver, a part of the RPW instrument, made during this flyby at frequencies below 200 kHz. The TDS receiver operated in a high cadence mode providing us with the regular waveform snapshot with 62 ms length every ten seconds for two electric components. Besides the regular captures, we have got more than five hundred onboard classified snapshots and the statistical products with a sixteen-second cadence. Before entering the terrestrial magnetosphere around 02:30UT, the spacecraft wandered through the foreshock region, registering intense bursts of Langmuir waves. After the bowshock crossing, Solar Orbiter was for more than two hours in the morning sector of the magnetosphere, recording various plasma wave modes. The closest approach was reached at 04:30UT above North Africa at an altitude of 460 km. Then the spacecraft continued into the Earth’s tail and entered the magnetosheath around 13:00UT. After 15:00UT, the Solar Orbiter crossed the bowshock, and bursts of Langmuir waves were detected again pointing out to the deep downstream foreshock region. Further from the Earth, intense Auroral Kilometric Radiation (AKR) at frequencies above 100 kHz was also detected.

How to cite: Pisa, D., Soucek, J., Santolik, O., Hanzelka, M., Maksimovic, M., Vecchio, A., Khotyaintsev, Y., Chust, T., Kretzschmar, M., Matteini, L., and Horbury, T.: Observations of the Time Domain Sampler receiver from the Radio and Plasma Wave instrument during the Solar Orbiter Earth flyby, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-4130,, 2022.

EGU22-8923 | Presentations | ST1.5

Origin of the solar wind observed by PSP 

Jasmina Magdalenic, Senthamizh Pavai Valliappan, and Luciano Rodriguez

The recently available novel Parker Solar Probe (PSP) observations allow mapping of the solar wind characteristics in the low solar corona, at only a few tens of solar radii, and study the characteristics of the solar wind and its transients close to the Sun. The first few perihelion passes of the PSP revealed the highly variable structure of the solar wind. These observations provided the unique possibility to study the solar wind originating from small coronal holes. Such studies were previously not possible as the majority of the in situ observations were taken at distances of about 1 AU, and the association of the solar wind originating from the small coronal holes and its sources on the Sun was extremely unreliable. We employ a magnetic connectivity tool (developed by ESA’s MADAWG group) to associate the solar wind parcels observed by the PSP with their source regions on the Sun.  Our study encompasses the first eight PSP perihelion passes, using a time window of about three weeks around each perihelion. The first results indicate that we can well distinguish and identify the solar wind observed by the PSP originating not only from the big, but also from the small coronal holes.

How to cite: Magdalenic, J., Valliappan, S. P., and Rodriguez, L.: Origin of the solar wind observed by PSP, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-8923,, 2022.

EGU22-8932 | Presentations | ST1.5

Solar wind modelling with EUHFORIA and comparison with the PSP observations 

Senthamizh Pavai Valliappan, Jasmina Magdalenic, Luciano Rodriguez, Stefaan Poedts, and Evangelia Samara

During recent decades, numerous studies were devoted to improve our understanding of the origin and the propagation of the fast solar wind, and its accurate modelling. The solar wind characteristics were poorly mapped up to now by the in situ observations, available mostly only at about 1AU. The novel Parker Solar Probe (PSP) observations, that are employed in this study, allow us to map the solar wind characteristics in the low solar corona at only few tens of solar radii, and to study the solar wind characteristics. These observations are also very important for understanding how accurately we can model the solar wind characteristics employing models such as EUHFORIA (European heliospheric forecasting information asset, Pomoell & Poedts, 2018), at distances rather close to the Sun.
In this study, we inspect the solar wind characteristics during the first eight closest PSP approaches to the Sun. The solar wind plasma characteristics observed by PSP are compared with the modelling results using the default set-up of EUHFORIA. We also calibrate the inner boundary (0.1 AU) conditions in EUHFORIA, but without changing the Wang-Sheeley-Arge formula which describes the solar wind characteristics at the inner boundary. The aim is to better model the solar wind at distances close to the Sun. We also use a magnetic connectivity tool (developed by ESA’s MADAWG group) to better associate the fast solar wind with its source region on the Sun.

How to cite: Valliappan, S. P., Magdalenic, J., Rodriguez, L., Poedts, S., and Samara, E.: Solar wind modelling with EUHFORIA and comparison with the PSP observations, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-8932,, 2022.

EGU22-9673 | Presentations | ST1.5

From Magnetic Reconnection at Chromospheric Network Boundaries to Switchbacks in the Inner Heliosphere 

Chuanpeng Hou, Jiansen He, Die Duan, Huichao Li, and Yajie Chen

During its perihelion encounter, Parker Solar Probe (PSP) has observed abundant kink and switchback patterns of magnetic field lines in the young solar wind. Such switchback structures have gained widespread attention due to their manifestation of Alfvenic and compressional properties, such as velocity enhancement, plasma density, and temperature variation. In particular, the origin mechanism has been a hot topic and waits to be observationally confirmed. Here we use a two-step ballistic backmapping method (tracing along the Parker Spiral and the PFSS-solution extrapolated from the GONG synoptic magnetogram) to determine the foot-points of the magnetic lines during switchback events measured by PSP on 24th-27th Jan. 2020. We identify ten jets corresponding to the PSP-switchbacks and find relevance between jets and switchback patches. We find that jets excite at the height of around low corona and position of chromospheric network boundaries. About 70% of jets accompany a magnetic cancelation, while 30% of jets are related to magnetic emergence. The variation in magnetic flux corresponding to magnetic cancelation and magnetic emergence is quantitatively equal to that of radial magnetic flux associated with switchback patches. These features suggest that switchbacks may originate from an interchange magnetic reconnection at chromospheric network boundaries, which provides direct evidence of switchbacks' solar origin. 

How to cite: Hou, C., He, J., Duan, D., Li, H., and Chen, Y.: From Magnetic Reconnection at Chromospheric Network Boundaries to Switchbacks in the Inner Heliosphere, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-9673,, 2022.

EGU22-12654 | Presentations | ST1.5

Non-thermal features in proton and alpha velocity distribution functions in the solar wind in the inner heliosphere 

Raffaella DAmicis, Roberto Bruno, Rossana De Marco, Marco Velli, Daniele Telloni, Denise Perrone, Luca Sorriso-Valvo, and Olga Panasenco

The solar wind is a highly variable, weakly collisional plasma originating from the Sun. The recent launches of PSP and Solar Orbiter have opened the way for the exploration of the innermost regions of our solar system and will greatly advance our understanding of several plasma phenomena occurring in the near-Sun environment, such as the heating and the acceleration of the solar wind. 

Plasma waves and wave-particle interactions play a relevant role in such phenomena, determining significant deviations of the Velocity Distribution Function (VDF) from the Local Thermodynamic Equilibrium. These deviations retain information of the interaction of particles with the turbulent electromagnetic fields and can be identified as thermal anisotropy, or non-thermal ion beams and heavy ion differential streaming in the ion component of the solar wind. This study will cover these topics with particular reference to new in-situ data from Solar Orbiter, PSP and with observations at L1  (e.g. Wind), with a focus on the central role Alfvénic fluctuations play in the evolution of the VDF features mentioned above.

How to cite: DAmicis, R., Bruno, R., De Marco, R., Velli, M., Telloni, D., Perrone, D., Sorriso-Valvo, L., and Panasenco, O.: Non-thermal features in proton and alpha velocity distribution functions in the solar wind in the inner heliosphere, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-12654,, 2022.

EGU22-1221 | Presentations | ST1.6

Diurnal anisotropy of polar neutron monitors, Dome C looks poleward 

Agnieszka Gil, Alexander Mishev, Stepan Poluianov, and Ilya Usoskin

Galactic cosmic rays (GCR) show a small local anisotropy detected as a diurnal variability of neutron monitor (NM) count rates. As the asymptotic directions of different NMs are diverse, the capability of the GCR diurnal variation observation is also various. Here we present that the Dome C (DOMC) NM is barely sensitive to the diurnal variation. Its amplitude is very small, 0.03%, in comparison to other polar NMs, for which the diurnal variability amplitudes vary from 0.16 to 0.4%. This fact is associated to the narrow asymptotic-direction cone of DOMC NM looking almost to the South pole with geographic latitude above 75o. Thus, DOMC NM is the only existing NM accepting cosmic-ray particles from the off-equatorial region, which makes this station a unique detector.

How to cite: Gil, A., Mishev, A., Poluianov, S., and Usoskin, I.: Diurnal anisotropy of polar neutron monitors, Dome C looks poleward, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-1221,, 2022.

The first solar proton event of solar cycle 25 was detected on 28 October 2021 by several neutron monitors (NMs) in the polar regions, the strongest signal was registered by the DOMC/DOMB monitors located at the Antarctic plateau at Concordia French-Italian research station. It is identified as the GLE (ground-level enhancement) #73 in the International GLE database. Here, we report the observations of the GLE by the global NM network and present the derived angular and spectral features of solar energetic protons with their dynamical evolution throughout the event employing a state-of-the-art model based on analysis of the neutron monitor data. Using the derived spectra we computed the related terrestrial effects, namely the cosmic rate induced ionization at several altitudes on a global map and discuss possible implications.

How to cite: Mishev, A., Larsen, N., and Usoskin, I.: The first GLE (# 73 – 28-Oct-2021) of solar cycle 25: a study of the related terrestrial effects using neutron monitor data, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-1243,, 2022.

EGU22-4215 | Presentations | ST1.6

Cosmic Ray Flux Correlation between McMurdo and Jang Bogo Neutron Monitor Stations vs. Time Lag 

Ekawit Kittiya, Waraporn Nuntiyakul, Achara Seripienlert, Paul Evenson, Alejandro Saiz, David Ruffolo, and Sueyon Oh

A neutron monitor is a large ground-based detector responding to the flux of cosmic ray particles in space by measuring atmospheric secondary neutrons. Any ground-based detector is sensitive to cosmic rays from a specific range of directions in space. In particular, a particle arriving from a specific sky direction with a specific rigidity (momentum per unit charge) was necessarily moving from a certain direction in space, called the asymptotic direction. McMurdo and Jang Bogo neutron monitor stations are Antarctic stations with similar geomagnetic latitudes but slightly different longitudes. From December 17, 2015 to January 9, 2017, six of the eighteen neutron counters from McMurdo had been transferred to Jang Bogo (with full transfer to Jang Bogo completed in December 2017). We present an analysis of the correlation of the cosmic ray flux between the McMurdo and Jang Bogo stations, during the time when both were operating, with ten-second time resolution. Although highly correlated, there are significant differences, including a systematic time lag of approximately 16 minutes between the data from the two stations. Although McMurdo observes with similar asymptotic directions to Jang Bogo, the response-weighted average directions still have a substantial difference of 21.9 degrees in geographic longitude, so with Earth’s rotation, time-independent anisotropy effects should induce a lag of 88 minutes.  Because the observed lag of 16 minutes is intermediate between 0 and 88 minutes, the joint observations reveal structure in the interplanetary cosmic-ray density that is consistent with a combination of simultaneous temporal variations and non-simultaneous variations with direction (i.e., anisotropy). The research is supported in part by a TA/RA scholarship (active recruitment) of Chiang Mai University and Thailand Science Research and Innovation via Research Team Promotion Grant RTA6280002.

How to cite: Kittiya, E., Nuntiyakul, W., Seripienlert, A., Evenson, P., Saiz, A., Ruffolo, D., and Oh, S.: Cosmic Ray Flux Correlation between McMurdo and Jang Bogo Neutron Monitor Stations vs. Time Lag, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-4215,, 2022.

EGU22-4352 | Presentations | ST1.6

Diurnal Dips and Forbush Decreases in Galactic Cosmic Rays Observed at High Cutoff Rigidity 

Chanoknan Banglieng, David Ruffolo, Alejandro Sáiz, Warit Mitthumsiri, Tanin Nutaro, Marc L. Duldig, and John E. Humble

A Forbush decrease (FD) is the decrease in Galactic cosmic ray (GCR) flux, e.g., as observed by a neutron monitor count rate, in association with a coronal mass ejection (CME) and/or its shock. The FD amplitude is known to decrease at higher cutoff rigidity. During Solar Cycle 24, the Mawson neutron monitor in Antarctica, with a low (atmosphere-limited) cutoff rigidity of ~1 GV, observed numerous FDs, while the Princess Sirindhorn Neutron Monitor (PSNM) located near the Earth's equator at Doi Inthanon, Thailand, with the world’s highest geomagnetic cutoff rigidity (17 GV), observed only a fraction of these as FDs.  Instead, we find that the shock arrival is often followed by repeated dips in the PSNM count rate at only certain times of day, while the GCR flux from other directions remains near the pre-shock level.  We refer to decreases of this type as diurnal dips.  In this work, we have surveyed FDs and diurnal dip events observed by PSNM and the dependence of their maximum amplitude on the solar wind speed, ICME speed, and interplanetary magnetic field.  We acknowledge logistical support from Australia's Antarctic Program for operating the Mawson NM and support from the National Astronomical Research Institute of Thailand and grant RTA6280002 from Thailand Science Research and Innovation.

How to cite: Banglieng, C., Ruffolo, D., Sáiz, A., Mitthumsiri, W., Nutaro, T., Duldig, M. L., and Humble, J. E.: Diurnal Dips and Forbush Decreases in Galactic Cosmic Rays Observed at High Cutoff Rigidity, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-4352,, 2022.

EGU22-5376 | Presentations | ST1.6

Boron-coated straw detector technology as an alternative to helium-3 and boron trifluoride based proportional counters for ground level neutron monitoring: a design study. 

Michael Aspinall, Cory Binnersley, Steve Bradnam, Stephen Croft, Malcolm Joyce, Lee Packer, and Jim Wild

The global network of neutron monitors comprises predominantly of the monitor standardised by Carmichael in 1964, the NM-64.  The design of these existing monitors and their instrumentation have changed very little over the last sixty years.  For example, their neutron detectors rely on gas filled proportional counters that are either filled with highly toxic boron trifluoride (BF3) or helium-3 (3He).  Almost the entire global supply of 3He is derived from a waste product of nuclear weapons programmes and, with the termination of such programmes and reducing nuclear weapons stockpile, the supply of 3He has become limited.  Consequently, 3He supply became strictly controlled in 2008 and its price has fluctuated since.  In some cases, new neutron monitors have reverted to BF3 filled counter tubes when the price of 3He has been at a premium.  Helium-3 filled proportional counters are also used extensively in radiation portal monitors deployed for homeland security and non-proliferation; objectives which have increased significantly over the last two decades.  The reduced production and increased demand for 3He has led to concerns over its supply and provided the research motivation for alternative neutron detection methods which are viable in terms of sensitivity, stability and gamma-rejection for certain applications.  One of these alternative technologies is based on boron-coated straws (BCS) manufactured and supplied by Proportional Technologies, Inc (PTI).  The technology is built on a patented low-cost technology that enables long copper tubes, known as ‘straws’, to be coated on the inside with a thin layer of 10B-enriched boron carbide (10B4C).  Thermal neutrons captured in the 10B are converted into secondary particles, through the 10B(n, α) reaction.  The straws can be of various diameter (circa 4 mm to 15 mm), length (up to 2 m) and shape (round, star or pie) to increase the surface area of 10B.  Multiple straws can be packed inside a 1” diameter aluminium tube acting as a single drop-in replacement for traditional 3He detectors or individually distributed directly throughout the moderating medium, thus increasing efficiency by detecting the thermal neutrons at the point that they are created.  BCS-based detectors are widely used in systems for homeland security, safeguards and neutron imaging in direct exchange for 3He tubes.  This study aims to design a neutron monitor utilising BCS technology that is cheaper, more compact and produces comparable results to the existing network of NM-64 monitors.  Monte Carlo simulations using the MCNP radiation transport code to model several BCS-based solutions and an NM-64 computational benchmark are reported.  These models are validated experimentally using a standard PTI portal monitor (PTI-110-NDME) to determine its efficiency, dieaway, deadtime and gamma rejection using a combination of bare 252Cf, AmLi and 137Cs sources.  The PTI-110-NDME consists of a 12” x 5” x 1 m high density polyethylene (HDPE) slab with thirty ~15-mm diameter straws, 93 cm active length, embedded uniformly throughout the moderator.  Funded by UK Research & Innovation (UKRI), this research is part of the Space Weather Instrumentation, Measurement, Modelling and Risk (SWIMMR) programme.

How to cite: Aspinall, M., Binnersley, C., Bradnam, S., Croft, S., Joyce, M., Packer, L., and Wild, J.: Boron-coated straw detector technology as an alternative to helium-3 and boron trifluoride based proportional counters for ground level neutron monitoring: a design study., EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-5376,, 2022.

EGU22-5504 | Presentations | ST1.6

Characteristics of the First Ground Level Enhancement (GLE) of Solar Cycle 25 on 28 October 2021 

Athanasios Papaioannou, Athanasios Kouloumvakos, Alexander Mishev, Rami Vainio, Ilya Usoskin, Konstantin Herbst, Alexis P. Rouillard, Anastasios Anastasiadis, Jan Giesler, Robert Wimmer-Schweingruber, and Patrick Kuhl

We present an overview of the first ground-level enhancement (GLE) event of solar cycle 25, recorded on 28 October 2021 (GLE73), based on the available neutron monitor (NM) network observations and on data from near-Earth spacecraft (GOES, SOHO, SolO). The maximum increase was ~7.3% for DOMC (Dome C NM at Concordia station) and 5.4% for SOPO (South Pole) conventional NMs located on the Antarctic plateau. Bare (lead-free) NMs at the same sites detected a higher response (14.0% for DOMB and 6.6% for SOPB). The Fort Smith (FSMT) NM shows the earliest increase among the high-latitude NMs, indicating a moderate anisotropy in the first phase of the GLE event. The maximum rigidity of accelerated protons did not exceed 2.4 GV. We estimated the solar release time (SRT) of ≥1 GV protons into open magnetic field lines at ~15:40 UT.  In-situ proton observations from near-Earth spacecraft were combined with the detection of a solar flare in soft X-rays (SXRs), a coronal mass ejection (CME), radio bursts and extreme ultraviolet (EUV) observations to identify the solar origin of the GLE. Around the ≥1 GV proton SRT the CME-driven shock was located at a height of ~2.33 Rs. The timing of the EUV wave evolution towards the field lines magnetically connected to Earth seem to be in good agreement with the inferred release time of ≥1 GV protons.

How to cite: Papaioannou, A., Kouloumvakos, A., Mishev, A., Vainio, R., Usoskin, I., Herbst, K., Rouillard, A. P., Anastasiadis, A., Giesler, J., Wimmer-Schweingruber, R., and Kuhl, P.: Characteristics of the First Ground Level Enhancement (GLE) of Solar Cycle 25 on 28 October 2021, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-5504,, 2022.

EGU22-5650 | Presentations | ST1.6

Coincident signals on multiple counters in a neutron monitor: Sparse vs. dense atmospheric secondary particles from cosmic ray showers 

Kullapha Chaiwongkhot, David Ruffolo, Poompong Chaiwongkhot, Alejandro Sáiz, and Chanoknan Banglieng

Neutron monitors were designed to measure atmospheric secondary neutrons from cosmic ray showers in order to track the cosmic ray flux vs. time.  Furthermore, at the Princess Sirindhorn Neutron Monitor (PSNM), an 18-counter NM64 detector at 2560-m altitude at Doi Inthanon, Thailand, the leader fraction (inverse multiplicity) inferred from time delay distribution between successive neutron events in the same counter has been used to track spectral variations.  More recent measurements of time delays between neutron events in different counters, as a function of counter separation, have confirmed that 1) the product neutrons from the interaction of a single atmospheric secondary neutron can spread among neighboring counters, with a cross-counter leader fraction that depends on whether the first counter is an end or middle counter, and 2) coincident counts between distant counters can be produced by multiple atmospheric secondary neutrons from the same primary cosmic ray, with a leader fraction that depends on whether the second counter is an end or middle counter.  Here we report on measurements of neutron signals in amplifier outputs at PSNM using a 4-channel oscilloscope, in order to further investigate these phenomena.  For a pair of neighboring counters located at the edge of the counter array, we find roughly equal event rates in either neighbor following a neutron trigger in one of them, implying that the difference in leader fraction relates to the base count rate of the first counter and is therefore lower if that is an end counter. In addition, an FPGA-based readout system was developed for more efficient collection of neutron events on two distant counters (Tubes 2 and 18) that were coincident within a 250-microsecond time window, while also monitoring Tube 10 in between.  The time distributions and neutron multiplicities indicate that a small fraction of the cosmic ray events that triggered both Tubes 2 and 18 also led to neutron events on Tube 10 with an enhanced rate of high multiplicity, indicating a few air shower events that densely “carpeted” the neutron monitor, while the majority of such coincidences apparently involved a sparse distribution of isolated secondary particles near Tube 2 and near Tube 18 and not near the intermediate Tube 10, which is consistent with a cross-counter leader fraction dependence on the count rate of the second counter.  This research was supported by the postdoctoral research sponsorship of Mahidol University, Thailand, by grant RTA6280002 from Thailand Science Research and Innovation, and by grant NSRF from the Program Management Unit for Human Resources & Institutional Development, Research and Innovation, NXPO, Thailand [grant number B05F640051]. 

How to cite: Chaiwongkhot, K., Ruffolo, D., Chaiwongkhot, P., Sáiz, A., and Banglieng, C.: Coincident signals on multiple counters in a neutron monitor: Sparse vs. dense atmospheric secondary particles from cosmic ray showers, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-5650,, 2022.

EGU22-6352 | Presentations | ST1.6

Multiple interactions in a Neutron Monitor 

Pierre-Simon Mangeard, John Clem, Paul Evenson, Waraporn Nuntiyakul, David Ruffolo, Alejandro Sáiz, Achara Seripienlert, and Surujhdeo Seunarine

The flux of low-energy (GeV-range) Galactic cosmic rays at Earth is modulated by the long term magnetic variations of the Sun (11-year sunspot cycle and 22-year magnetic solar cycle). This process known as Solar modulation is most pronounced at 1 GeV and below. However, it also operates at much higher energy, still exhibiting solar magnetic polarity dependence. For the last decades, ground-based neutron monitors provided valuable observations of the solar modulation up to a rigidity cutoff of about 17 GV. To extend the energy range of the neutron monitor observations, we recently upgraded the electronics of the Princess Sirindhorn Neutron Monitor in Thailand (PSNM, the operating neutron monitor at the highest geomagnetic rigidity cutoff) to record complex combinations of hits in multiple proportional counters. We present here the detection of multiple-hit events recorded at the PSNM. We discuss these observations with the help of a detailed Monte-Carlo simulation of energetic neutrons interacting in the detector. Finally, we estimate the nucleonic spectrum of the atmospheric secondary particles at the altitude of the detector.

How to cite: Mangeard, P.-S., Clem, J., Evenson, P., Nuntiyakul, W., Ruffolo, D., Sáiz, A., Seripienlert, A., and Seunarine, S.: Multiple interactions in a Neutron Monitor, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-6352,, 2022.

EGU22-11286 | Presentations | ST1.6

Integral fluences of GLEs: A new full reconstruction 

Ilya Usoskin, Sergey Koldobskiy, and Gennady Kovaltsov

For many reasons, from solar physics to terrestrial and engineering applications, it is important to estimate the total fluences of solar energetic particles (SEPs) especially for the strongest hard-spectrum events known as ground-level enhancements (GLEs). Here we present a revised reconstruction of the SEP spectral fluences using a recently developed probabilistic Monte-Carlo method, applied to major GLE events of the last decades. The method utilizes data from the ground-based neutron-monitor network in the higher-energy range and revised space-borne/ionospheric data for the lower-energy part. The fluences are reconstructed along with realistic uncertainty estimates which appear large for weak events and small for strong events.

How to cite: Usoskin, I., Koldobskiy, S., and Kovaltsov, G.: Integral fluences of GLEs: A new full reconstruction, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-11286,, 2022.

EGU22-11343 | Presentations | ST1.6

Characteristics of a long-lived CIR and analytical modelling of the corresponding depression in the GCR flux 

Mateja Dumbovic, Bojan Vrsnak, Bernd Heber, Manuela Temmer, Patrick Kuhl, and Anamarija Kirin

We observe a long-lived CIR recurring in 27 consecutive Carrington rotations 2057-2083 in the time period from June 2007 - May 2009. We characterize the in situ measurements of this long-lived CIR as well as the corresponding depression in the GCR count observed by SOHO/EPHIN, and analyze them throughout different rotations. We find that the inverted GCR count time-profile correlates well with that of the flow speed throughout different rotations. We perform a statistical analysis and find the GCR count amplitude correlated to the peak in the magnetic field and flow speed, as expected based on previous statistical studies. In order to characterize a generic CIR profile for modelling purposes, we perform the superposed epoch analysis using relative values of the key parameters. Based on the observed properties we propose a simple analytical model starting from the basic Fokker-Planck equation. We employ a convection-diffusion GCR propagation model and apply it to the solar wind and interplanetary magnetic field properties observed for the analyzed long-lived CIR. Our analysis demonstrates a very good match of the model results and observations.

How to cite: Dumbovic, M., Vrsnak, B., Heber, B., Temmer, M., Kuhl, P., and Kirin, A.: Characteristics of a long-lived CIR and analytical modelling of the corresponding depression in the GCR flux, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-11343,, 2022.

EGU22-13123 | Presentations | ST1.6

Galactic cosmic rays as signatures of interplanetary transients 

Mateja Dumbovic

Coronal mass ejections (CMEs), interplanetary shocks, and corotating interaction regions (CIRs) drive heliospheric variability, causing various interplanetary as well as planetary disturbances. One of their very common in-situ signatures are short-term reductions in the galactic cosmic ray (GCR) flux (i.e. Forbush decreases), which are measured by ground-based instruments at Earth and Mars, as well as various spacecraft throughout the heliosphere (most recently by Solar Orbiter). In general, interplanetary magnetic structures interact with GCRs producing depressions in the GCR flux. Therefore, different types of interplanetary magnetic structures cause different types of Forbush decreases, allowing us to distinguish between them. With new modelling efforts, as well as observational analysis we are one step closer in utilizing GCR measurements to provide information on interplanetary transients, especially where other measurements (e.g. plasma, magnetic field) are lacking.

How to cite: Dumbovic, M.: Galactic cosmic rays as signatures of interplanetary transients, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-13123,, 2022.

EGU22-13288 | Presentations | ST1.6

Measuring neutron monitor multiplicities at SANAE 

Du Toit Strauss

We present new results of the neutron monitor (NM) multiplicity, as measured by the SANAE NM with updated electronics, down to 10 microseconds. We identify the high-multiplicity component, formed when high-energy particles interact with the NM and produce multiple neutrons in the lead producer. This component is absent in the lead-free monitors and is absent when testing the leaded NM with a low-energy neutron source. We study the pressure dependence of both the high- and low-multiplicity components, as well as the ratio thereof. We show how this ratio, as a proxy for the energy spectrum of atmospheric particles incident on the NM, changes during a relatively small Forbush decrease observed in November 2021.

How to cite: Strauss, D. T.: Measuring neutron monitor multiplicities at SANAE, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-13288,, 2022.

EGU22-1608 | Presentations | ST1.7

High-frequency waves driven by pickup ion ring-beam distributions in the outer heliosheath 

Kaijun Liu, Ameneh Mousavi, and Sina Sadeghzadeh

Scattering of pickup ion ring-beam distributions in the outer heliosheath is a fundamental element in the spatial retention scenario of the energetic neutral atom (ENA) ribbon observed by the Interstellar Boundary EXplorer (IBEX). According to our earlier linear instability analysis, pickup ion ring-beam distributions trigger magnetic field-aligned, right-hand polarized unstable waves in two separate frequency ranges which are near and far above the proton cyclotron frequency, respectively. We have performed hybrid simulations to study the unstable waves near the proton cyclotron frequency. However, the high-frequency waves well above the proton cyclotron frequency are beyond the reach of hybrid simulations. In the present study, particle-in-cell simulations are carried out to investigate the parallel- and anti-parallel-propagating high-frequency waves excited by the outer heliosheath pickup ions at different pickup angles as well as the scattering of the pickup ions by the waves excited. In the early stages of the simulations, the results confirm the excitation of the parallel-propagating, right-hand polarized high-frequency waves as predicted by the earlier linear analysis. Later in the simulations, enhanced anti-parallel-propagating modes also emerge. Furthermore, the evolution of the pickup ion ring-beam distributions of the selected pickup angles reveals that the high-frequency waves do not significantly contribute to the pickup ion scattering. These results are favorable regarding the plausibility of the spatial retention scenario of the IBEX ENA ribbon.

How to cite: Liu, K., Mousavi, A., and Sadeghzadeh, S.: High-frequency waves driven by pickup ion ring-beam distributions in the outer heliosheath, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-1608,, 2022.

EGU22-3325 | Presentations | ST1.7

The Dispersive Nature of the Heliospheric Termination Shock 

Bertalan Zieger, Joe Giacalone, Marc Swisdak, Gary Zank, and Merav Opher

The Voyager spacecraft are the first man-made objects to cross the termination shock (TS), where the solar wind becomes sub-fast magnetosonic due to the interaction with the local interstellar medium. Voyager 2 observations revealed that classical single-fluid magnetohydrodynamic (MHD) or multispecies single-fluid MHD models are not sufficient to describe the microstructure of the TS and the observed nonlinear waves downstream the TS. Consequently, more sophisticated physical models, like multifluid, hybrid or fully kinetic solar wind models, are needed to capture nonlinear waves, dispersive shock waves, and ion-ion instabilities, where each ion species (and electrons) can move independently with their own bulk velocities, and the fluctuating parts of the ion velocities are often comparable to the mean velocity of the collective plasma fluid. The multifluid simulation of the TS by Zieger et al. [2015] shows a remarkable agreement with high-resolution Voyager 2 observations, reproducing not only the microstructure of the third TS crossing (TS3) but also the energy partitioning among thermal ions, pickup ions (PUI), and electrons across the shock. It was demonstrated that TS3 is a subcritical dispersive shock wave with low fast magnetosonic Mach number and high plasma ß. Here we present multifluid, hybrid, and particle-in-cell (PIC) simulations of the second TS crossing (TS2) by Voyager 2, which was somewhat stronger than TS3, with an observed compression ratio of 2.2. All three types of simulations confirm the dispersive nature of the TS in agreement with Voyager 2 observations. We conclude that TS2, just as TS3, is a subcritical dispersive shock wave with a soliton (overshoot) at the leading edge of the shock and a quasi-stationary nonlinear wave train downstream of the shock front. We compared the cross-shock electric field in the multifluid, hybrid, and PIC simulations and found a reasonable agreement. We show that the Hall electric field is dominating over the convective and ambipolar electric fields, which indicates that electrons play an important role in the shock transition. Finally, we demonstrate that the microstructure of the termination shock is controlled by dispersion rather than ion reflection, and only slightly affected by reflected solar wind ions in the hybrid and PIC simulations, which validates the multifluid model on fluid scale. The dispersive nature of the termination shock has important implications for the transition and acceleration of PUIs across the termination shock, which is revealed in the PUI distributions in our hybrid [Giacalone et al., 2021] and PIC simulations.

How to cite: Zieger, B., Giacalone, J., Swisdak, M., Zank, G., and Opher, M.: The Dispersive Nature of the Heliospheric Termination Shock, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-3325,, 2022.

The Voyager 1 and Voyager 2 (V1 & V2) crossings of the termination shock (TS; ~94 and ~84 AU, respectively), led to the first measurements of ions and electrons that constitute the heliosheath (HS). Their crossings of the heliopause (HP; ~122 AU and ~119 AU), pinpointed the extent of the upwind heliosphere's expansion into the Very Local Interstellar medium (VLISM). The Cassini/INCA >5.2 keV ENA images of the celestial sphere, have placed the local V1&2/LECP measurements in a global context and have led to the discovery of a high intensity and wide ENA region that encircles the celestial sphere, called “Belt” and corresponds to a “reservoir” of particles that exist within the HS. The heliosphere forms a time-dependent, roughly symmetric obstacle to the inward interstellar flow, responding within ~2-3 yrs, in both the nose and anti-nose directions to the outward propagating solar wind changes through the solar cycle. The shape of the ion energy spectra plays a critical role in determining the pressure balance and acceleration mechanisms inside the HS. Energy spectra from ~10 eV to 344 MeV show that the PUIs dominate the total pressure distribution inside the HS, but suprathermal ions provide a significant contribution that cannot be neglected, revealing that >5.2 keV ENAs serve as important indicators of the acceleration processes that the parent H+ population undergoes inside the HS, thus imposing a key constraint on any future interpretation concerning the HS dynamics. The combination of ENAs and ions in the HS show that the plasma beta is >>1, the magnetic field upstream at the heliopause required to balance the pressure from the HS is >0.5 nT (V1 direction) and ~0.67 nT (V2 direction) and that the neutral Hydrogen density is ~0.12/cm3. These inferred values are consistent with measurements from both V1 and V2 spacecraft. Energetic ion measurements from V1/LECP in and beyond the HP show an average radial inflow of 40-139 keV ions for ~10 AU inside the HS and an average radial outflow over a spatial range of ~28 AU past the HP. These particles correspond to an ion population leaking from the HS into interstellar space, most likely due to the flux tube interchange instability at the boundary and provide a direct observation of the communication between the HS and the VLISM. They may also provide an important constraint for future models that aim to explain the <6 keV ENA ribbon fluxes (measured from the IBEX mission), which are likely formed from the neutralization of energetic pickup ions gyrating in the IS magnetic field outside the HP, reflecting (in part) those ion distributions that are responsible for the formation of this unexpected structure.

How to cite: Dialynas, K., Krimigis, S., Decker, R., and Hill, M.: Highlights of the combination of “Ground truth” >28 keV in-situ ions from the Voyagers and >5.2 keV ENAs from Cassini in the study of the global heliosphere, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-4316,, 2022.

EGU22-6402 | Presentations | ST1.7

Interstellar Neutral He Parameters from Crossing Parameter Tubes with the Interstellar Mapping and Acceleration Probe (IMAP) informed by 10 Years of Interstellar Boundary Explorer (IBEX) Observations 

Nathan Schwadron, Eberhard Moebius, David McComas, Jon Bower, Maciek Bzowski, Stephen Fuselier, David Heirtzler, Marzena Kubiak, Marty Lee, Fatemeh Rahmanifard, Justyna Sokol, Pawel Swaczyna, and Reka Winslow

The Sun’s motion relative to the surrounding interstellar medium leads to an interstellar neutral (ISN) wind through the heliosphere. For several species, including He, this wind is moderately depleted by ionization and can be analyzed in-situ with pickup ions and direct neutral atom imaging.  Since 2009, observations of the wind at 1 AU with the Interstellar Boundary Explorer (IBEX) have returned a precise 4-dimensional parameter tube for the flow vector (speed VISN, longitude λISN, and latitude βISN) and temperature TISN of interstellar He in the local cloud, which organizes VISN, λISN, and TISN as a function of λISN, and the local flow Mach number (Vth−ISN/VISN).  We refer to this functional dependence as the 4D IBEX parameter tube. On IBEX, the limitation of measuring the ISN flow observations to nearly perpendicular to the Earth-Sun line limits the range of observations in ecliptic longitude to ∼ 30º.  This limitation results in large uncertainties along the IBEX parameter tube and relatively small uncertainties across the parameter tube.   Over the past three years, IBEX operations were modified to let the spin axis pointing of IBEX drift to the maximum offset (7º) west of the Sun, which is the limit for the IBEX spacecraft. This expansion of the IBEX viewing helps break the degeneracy of the ISN parameters along the 4D IBEX parameter tube. It complements the full χ-square-minimization to obtain the ISNs parameters through comparison with detailed models of the ISN flow. The next generation IBEX-Lo sensor on IMAP will be mounted on a pivot platform, enabling IMAP-Lo to follow the ISN flow over almost the entire spacecraft orbit around the Sun.  A near-continuous set of 4D parameter tubes on IMAP will be observed for He, and for O, Ne, and H that cross at varying angles in the full ISN parameter space. This analysis substantially reduces the flow parameter uncertainties for these species and mitigating systematic uncertainties, such as those from ionization effects and the presence of secondary components. We discuss implications of these measurements for understanding our environment and its relationship to the structure of the local interstellar medium. Thus, we discuss how IMAP will probe the interstellar neutral gas flow in detail to derive the precise parameters of the interstellar flow and relate these conditions to understand our place within the interstellar medium. 

How to cite: Schwadron, N., Moebius, E., McComas, D., Bower, J., Bzowski, M., Fuselier, S., Heirtzler, D., Kubiak, M., Lee, M., Rahmanifard, F., Sokol, J., Swaczyna, P., and Winslow, R.: Interstellar Neutral He Parameters from Crossing Parameter Tubes with the Interstellar Mapping and Acceleration Probe (IMAP) informed by 10 Years of Interstellar Boundary Explorer (IBEX) Observations, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-6402,, 2022.

EGU22-6502 | Presentations | ST1.7

Filtration and Scattering of Interstellar Neutral Helium beyond the Heliopause 

Paweł Swaczyna, Maciej Bzowski, Stephen Fuselier, André Galli, Jacob Heerikhuisen, Marzena Kubiak, David McComas, Eberhard Möbius, Fatemeh Rahmanifard, Nathan Schwadron, and Eric Zirnstein

The pristine very local interstellar medium (VLISM) is not available for in situ observations even with the Voyager spacecraft because the influence of the heliosphere on the VLISM plasma extends up to several hundred au from the Sun. Therefore, observations of interstellar neutral (ISN) helium, the species least modified at the heliospheric boundaries, are used to determine the pristine VLISM flow speed, direction, and temperature. For more than one solar cycle, the Interstellar Boundary Explorer (IBEX) has sampled ISN helium atoms at 1 au, significantly reducing the statistical uncertainties of the ISN helium flow parameters. Launching in 2025, the Interstellar Mapping and Acceleration Probe (IMAP) will further lower these uncertainties thanks to the utilization of a pivot platform, which provides a range of viewing orientations and reduces the parameter degeneracy seen from IBEX data. Even though interstellar helium is the least modified ISN species, recent studies show that the ISN helium flux is affected by charge exchange and elastic collisions beyond the heliopause. Charge exchange collisions outside the heliopause filter the primary ISN helium and produce a secondary population from perturbed He+ ions in the interstellar plasma. The secondary helium population, originally called the Warm Breeze, was discovered from IBEX observations. Moreover, the distribution function of the primary ISN helium population is modified by elastic collisions with slowed down and heated plasma ahead of the heliopause. Consequently, the combined primary and secondary ISN helium populations at 1 au are complex and cannot be separated. Furthermore, the modifications of the population properties are larger than the statistical uncertainty of IBEX observations. We use global heliosphere models to estimate the magnitude of the filtration and scattering caused by charge exchange and elastic collisions. Together with other sources of information about the VLISM, these estimates allow us to assess the pristine VLISM conditions outside of heliosphere influences. 

How to cite: Swaczyna, P., Bzowski, M., Fuselier, S., Galli, A., Heerikhuisen, J., Kubiak, M., McComas, D., Möbius, E., Rahmanifard, F., Schwadron, N., and Zirnstein, E.: Filtration and Scattering of Interstellar Neutral Helium beyond the Heliopause, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-6502,, 2022.

EGU22-11016 | Presentations | ST1.7

Recent results from neutral beam source calibration by means of the novel Absolute Beam Monitor 

Jonathan Gasser, André Galli, and Peter Wurz

The IMAP mission by NASA is dedicated to extending the physical understanding of our heliosphere and its interaction with the interstellar medium by enhancing and refining the results obtained from IBEX. The neutral atom analysis instrument IMAP-Lo will observe and map fluxes of low-energy heliospheric neutral atoms (ENAs) and interstellar neutral (ISN) H, D, He, O and Ne with energies as low as 10 eV up to 1000 eV.

The instrument testing and calibration with a neutral atom beam is foreseen in the MEFISTO test facility for ion and neutral particle instruments at the University of Bern. MEFISTO is equipped with an electron-cyclotron resonance ion source that provides ion beams at a beam energy 3keV/q up to 100 keV/q. The beam fed into the test chamber is decelerated to 10 eV/q – 3 keV/q and effectively neutralized in a removable neutralization stage via surface reflection on a highly polished single crystal tungsten surface. The relative neutral beam intensity is permanently monitored via the neutralizing surface current. The neutralization process induces a considerable reduction of particle kinetic energy and conical widening of the neutral beam.

Thus, one key improvement for the calibration of a neutral atom instrument such as IMAP-Lo is to be able to measure the absolute neutral particle flux and beam energy into the instrument in the test chamber. To achieve this goal, the Absolute Beam Monitor (ABM) was developed recently.

The ABM is a dedicated laboratory device for absolute neutral particle flux measurements below 3 keV and coarse kinetic energy determination. Neutrals entering the ABM aperture strike a single crystal W conversion surface at grazing angle and are reflected into a channeltron to generate a stop pulse. The simultaneous monitoring of secondary electrons released at the W surface as start signal, the stop signal and the coincidence event rate allows inferring the rate of neutral atoms into the ABM aperture. The average neutral beam energy is obtained from the start-stop time-of-flight spectrum. The ABM is the first and so far the only device to measure the absolute neutral atoms flux in this low energy range below a few 100 eV. It serves as a primary standard for gauging the MEFISTO neutral beam source.

We report on recent calibration results of neutral H, He, O, Ne beams in the 10 eV – 1 keV energy range with the ABM in MEFISTO.

How to cite: Gasser, J., Galli, A., and Wurz, P.: Recent results from neutral beam source calibration by means of the novel Absolute Beam Monitor, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-11016,, 2022.

The heliosphere surrounding our solar system is formed by the interaction between the solar wind and the local interstellar medium as the Sun moves through interstellar space. With dimensions on the order of hundreds to potentially thousands of au, it is extremely difficult to pinpoint the 3D structure of the heliosphere and its boundaries, and the properties of the plasma within it. However, we can remotely measure the properties of the heliosphere with energetic neutral atoms (ENAs) which are created as a product of charge exchange between interstellar neutrals and ions within the solar wind plasma. ENAs can propagate hundreds of au before ionizing, allowing us to remotely view the distant boundaries of the heliosphere.

The Interstellar Boundary Explorer (IBEX) mission, a NASA smaller explorer mission which has been measuring ENA fluxes at ~0.5-6 keV for more than a solar cycle, has revealed at least two separate sources of ENAs: the “ribbon” of enhanced ENA fluxes forming a narrow, circular band across the sky, and the “globally distributed flux” (GDF) that forms lower intensity, broad features near the nose and tail of the heliosphere. While it is believed that the ribbon is formed from secondary ENAs from outside the heliopause, and ENAs from the inner heliosheath inside the heliopause are a major contributor to the GDF, it is not clear exactly how much of the GDF may originate outside the heliopause, and how much of the flux observed in the Ribbon is from the GDF itself. To help solve this issue, we present recent developments in our understanding of the ribbon vs. GDF sources from both IBEX data analysis and modeling perspectives. Moreover, with the upcoming IMAP mission set to launch in 2025, we provide insight into how IMAP’s enhanced capabilities may improve our current understanding of the heliosphere.

How to cite: Zirnstein, E.: Differentiating the Ribbon and Globally Distributed ENA Flux in IBEX Observations, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-13186,, 2022.

As the Sun moves through its local galactic neighborhood, it disturbs the ionized component of interstellar matter, forming a complex bow-wave structure of a slowed and heated, magnetized plasma, flowing past the heliosphere. This region is called the outer heliosheath. The neutral component of interstellar matter is in equilibrium with the ionized component far ahead of the heliosphere, but within the outer heliosheath it decouples kinematically from the perturbed plasma. A complex interaction begins, mostly by charge exchange, between the ions and the neutral atoms, and as a result, a new population of neutral atoms is created, with kinematic parameters similar to that of the surrounding plasma. In addition, elastic collisions operate, additionally modifying both the original primary and the secondary population of neutral. Both the primary and the secondary populations penetrate inside the heliosphere and can be directly sampled at 1 au. We will review results of observations of these populations obtained from IBEX-Lo and results of modeling of the secondary populations of interstellar hydrogen, helium, and oxygen in the outer heliosheath. We will illustrate how these populations can be better observed owing to enhanced capabilities of the planned IMAP-Lo instrument and demonstrate how they can be leveraged to resolve the primary and secondary populations, and investigate the properties of local interstellar matter in the immediate neighborhood of the Sun.

How to cite: Bzowski, M.: The primary and secondary populations of interstellar neutral gas as seen now by IBEX and in the future by IMAP, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-13535,, 2022.

EGU22-13568 | Presentations | ST1.7

Pressure Balance at the Heliosphere Boundary and in the Local Interstellar Cloud 

Eberhard Möbius and Jeffey Linsky

The pressures exerted by the solar wind from the inside and the interstellar medium from the outside control the size and shape of the heliosphere. Magnetic fields and cosmic rays play important roles to varying extents on both sides. We will assess the pressure balance by assembling the relevant component pressures in the (inner) heliosheath inside the heliopause, in the Very Local Interstellar Medium (VLISM, interstellar medium affected by the presence of the heliosphere or outer heliosheath), and in the Local Interstellar Cloud (LIC) unaffected by the heliosphere. We take the cosmic ray pressure from Voyager observations, which don’t show substantial gradients beyond the heliopause. The remaining pressure on the heliopause from inside is due to thermal and suprathermal ions in the subsonic solar wind, obtained from IBEX and INCA ENA observations and Voyager in situ measurements at the higher energies.

Besides cosmic rays, the pressure in the undisturbed LIC is composed of magnetic field pressure taken from IBEX Ribbon observations and related modeling. Thermal and turbulent pressures are based on H and He neutral and ion densities from pickup ion and interstellar gas flow observations, combined with the temperature and turbulent speed from absorption-line observations. The total LIC pressure in its rest frame is almost 40% lower than the pressure inside the heliopause, whereas adding the full ram pressure based on the LIC velocity relative to the Sun exceeds that pressure substantially. We estimate the likely effective pressure on the heliopause by combining the compressed interstellar magnetic field, as measured by Voyager, and the compressed and heated interstellar plasma, resorting to results from global heliospheric modeling. An interesting result of these pressure comparisons is that the effective ram pressure on the heliopause is somewhat larger than the combined magnetic field, thermal, and turbulent pressure in the LIC, which points to the importance of the LIC ram pressure for the shape of the heliosphere. We also compare the LIC pressure with the gravitational pressure on the galactic disk at the location of the Sun.

How to cite: Möbius, E. and Linsky, J.: Pressure Balance at the Heliosphere Boundary and in the Local Interstellar Cloud, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-13568,, 2022.

Space plasmas reside in non-equilibrium stationary states described by kappa distributions. The high-energy asymptotic behavior of kappa distributions leads to a power-law relationship of the energy-flux spectra; this relationship, when observed, can be analyzed for determining the thermodynamic parameters of the plasma. In the presented analysis, we use Energetic Neutral Atom (ENA) observations from the IBEX-Hi sensor, converted to the corresponding proton plasma spectra of the inner heliosheath, and timestamped with the ENA creation time. We, then, model the proton spectra with kappa distributions and derive the sky maps of the (radially averaged) values of temperature, density, kappa, and other thermodynamic parameters of the proton plasma in the inner heliosheath. We examine the variations of the determined thermodynamics and whether a correlation exists with solar activity during the 24th solar cycle.

How to cite: Livadiotis, G.: Thermodynamics of the proton plasma in the inner heliosheath during the 24th solar cycle, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-13574,, 2022.

EGU22-564 | Presentations | ST1.8

Study on the association of solar gyroresonance emission sources with brightness temperature intensification at 17 GHz 

Abimael Amaro, Marlos Rockenbach da Silva, and Joaquim Eduardo Rezende Costa

Remarkable works done in the last decades by many authors on the solar gyroresonance mechanism have illuminated the way to establish the relationship between this form of emission and magnetic fields in the solar atmosphere and to know the magnetic nature of the middle and upper layers of the active regions. Despite all these advances, solar physics still needs a direct means (without magnetograms) of identifying the sources of gyroresonance emission.
In search of a solution to this problem, we used solar images at 17 GHz synthesized by the Nobeyama Radioheliograph (NoRH) to map the likely sources of gyroresonance. 
To achieve this result, we first hypothesized that gyroresonance and bremsstrahlung mechanisms can generate a large brightness temperature intensification due to the close relationship such mechanisms have with magnetic fields and because of the role of the magnetic field in controlling the brightness of the solar atmosphere in the radiofrequency range.
To test this hypothesis regarding the gyroresonance process, we selected 8 large active regions (ARs) among the HMI magnetograms generated by the Solar Dynamics Observatory (SDO) corresponding to the 1st half of the 24th solar cycle.  We then analyzed each AR through its magnetogram and its image at 17 GHz. Aiming to verify, in these radio maps, whether there is a correspondence between brightness bumps and parameters associated with gyroresonance emission, we constructed three categories of brightness maps for each active region, respectively presenting the field of: brightness temperature in BT maps, brightness temperature gradient in BTG maps, and brightness temperature gradient to brightness temperature ratio in BTG/BT maps. Such parameters are the characteristic circular polarization, whose modulus is greater than 30%, and characteristic magnetic field strengths, associated with the gyroresonance radiation at 17 GHz for 3rd and 4th harmonics. Such a step also aimed to verify which of these categories would best map the putative sources of gyroresonance emission.
On these maps, we then plotted the contours of the characteristic parameters. 
For each AR, we also obtained the degree of correlation between its brightness variables and the characteristic polarization. We then observed that the contours of characteristic magnetic field strengths are predominantly enveloped by the area of the brightness bumps, while the contours of the characteristic polarization fit well to such areas, being better fitted to the relative bumps (in the BTG/BT maps). In the statistical analysis, we observe that for each active region, there is a predominance of strong correlations between the brightness variables and the modulus of the characteristic polarization. Such correlation tends to be highest for the brightness temperature. For each brightness variable, the highest correlations tend to occur for the predominant polarization direction of the active region.
The data, therefore, indicate a high probability that the gyroresonance emission mechanism was at least one of the important causes of the radio bumps produced in the observed active regions.

How to cite: Amaro, A., Rockenbach da Silva, M., and Rezende Costa, J. E.: Study on the association of solar gyroresonance emission sources with brightness temperature intensification at 17 GHz, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-564,, 2022.

EGU22-730 | Presentations | ST1.8

Estimation of the solar wind extreme events. 

Carlos Larrodera, Lidia Nikitina, and Consuelo Cid

This research provides an analysis of extreme events in the solar wind and in the magnetosphere due to disturbances of the solar wind. Extreme value theory has been applied to a 20-year data set from the Advanced Composition Explorer spacecraft for the period 1998–2017. The solar proton speed, solar proton temperature, solar proton density, and magnetic field have been analyzed to characterize extreme events in the solar wind. The solar wind electric field, vBz has been analyzed to characterize the impact from extreme disturbances in the solar wind to the magnetosphere. These extreme values were estimated for 1-in-40- and 1-in-80-year events, which represent two and four times the range of the original data set. The estimated values were verified in comparison with measured values of extreme events recorded in previous years. Finally, our research also suggests the presence of an upper boundary in the magnitudes under study.

How to cite: Larrodera, C., Nikitina, L., and Cid, C.: Estimation of the solar wind extreme events., EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-730,, 2022.

EGU22-939 | Presentations | ST1.8

Covariation of solar granulation size and sunspot indices in activity cycle 24 

Alexey Sharov and Arnold Hanslmeier

New theoretical arguments and empirical evidence for correlated changes in sun granulation size and sunspot indices were obtained during this two-year study using blue continuum image data acquired by the Hinode solar optical telescope in 2006-2016 and simultaneous time series of daily sunspot numbers (SSN) and areas (SSA). An original set of simple scaling equations linking the relative variation in the average size of granular cells to the sum change in the total number of granular cells, SSA and surface gravity was written under the assumption that at the height of the blue continuum formation, the entire surface of the Sun is the sum of the areas of granular cells, pores and sunspots. The magnitudes of relative changes in the horizontal size of granular cells due to variations in global solar parameters were estimated and it was shown that variations in the SSA have the dominant influence on variations in the mean granular size. Periodic spurious changes in the mean granular size due to the sun's variable tilt with respect to the telescope and apparent changes in rotational speed and surface gravity were also mentioned.

Our empirical research focussed on the automatic identification and precise morphometric measurements of granular cells in high-resolution Hinode images using an efficient marker-controlled watershed segmentation algorithm. A total of seven image sequences with a 30-, 27- and 1-day cadence, all of which contained 840 images, were compiled, processed and corrected with regard to the variable sun-earth distance and heliographic coordinates. The resultant granulation parameters, including mean area, equivalent diameter, extent and contrast were compared to SSN and SSA data using temporal cross-correlation. An essential anti-correlation was measured between the mean size of granular cells and the daily SSN for the medium time intervals of four months, which were characterized with the optimal orbital conditions for the imaging of solar granulation. The decrease in the cellular scale by 3% with the increase in the average SSN index of 20%, was revealed for both ascending and descending phases of the 24th activity cycle, while the average contrast characterizing image quality remained almost unchanged. The approximately 14-day delay in the cause-effect relationship between the SSN and the granulation scale was revealed and a plausible explanation for this delay was given.

In the quiet sun at the disk centre, the mean equivalent diameter of granules was measured at 1.37 arcseconds, while the same parameter for granular cells was given as 2.06 arc seconds, which was in good agreement with the results of other researchers. In the close vicinity of sunspots, the mean size of granular cells decreased by a few percent, while the width of intergranular lanes decreased by 60%, indicating the existence of bright rings with higher temperatures around these sunspots. The significance of all these observations is that they confirm the results of the predecessors, positively support the 70-year-old hypothesis about the dependence of granulation properties on the sunspot cycle, and stimulate the use of the granulation scale as a cyclical proxy for solar activity over medium-term periods.

How to cite: Sharov, A. and Hanslmeier, A.: Covariation of solar granulation size and sunspot indices in activity cycle 24, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-939,, 2022.

EGU22-1364 | Presentations | ST1.8

Variations of alpha particle parameters across corotating rarefaction regions 

Tereza Durovcova, Jana Safrankova, and Zdenek Nemecek

A corotating rarefaction region (CRR) is formed between the slow solar wind stream and the fast stream in front of it. It is associated with a region of small longitudinal extent on the solar surface and plasma parameters would correspond to a smooth transition between the interacting streams. However, previous studies have shown that the stream interaction can change parameters of individual solar wind components in an unexpected way. In our study, we focused on behavior of the second most abundant ionic component, alpha particles. Preliminary analysis of measurements at 1 AU revealed large variations of the alpha relative abundance (AHe) that did not correlate with the solar wind speed and alpha-proton relative drift changes. To determine the global profile of alpha properties across CRR at 1 AU, we performed the superposed epoch analysis of identified CRRs. We found continuous but spatially separated transitions of the alpha-proton relative abundance, relative drift, and alpha-proton temperature ratio from values corresponding to fast streams to those typical for slow streams. Despite the expectation that AHe would decrease from the beginning of the rarefaction, it corresponds to values usually observed in the fast solar wind for a large part of the CRR. Moreover, AHe is often enhanced near the expected stream interface. Therefore, we propose that a major part of the CRR is filled by the wind from a coronal hole and discuss various scenarios that could explain the obtained profiles of alpha particle parameters across CRRs. Finally, we use CRR observations at different radial distances from the Sun (Solar Orbiter, Ulysses, etc.) to identify effects connected with a radial evolution of these large-scale structures.

How to cite: Durovcova, T., Safrankova, J., and Nemecek, Z.: Variations of alpha particle parameters across corotating rarefaction regions, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-1364,, 2022.

EGU22-1856 | Presentations | ST1.8

A Study of Small-Scale Brightenings using EUV Data from SPICE on board Solar Orbiter 

Jenny Marcela Rodriguez Gomez, Peter Young, and Therese Kucera

Recently small-scale solar dynamic features have been observed on the Sun by Solar Orbiter/EUI (Berghmans et al. 2021; Panesar et al. 2021). These events (initially referred to as “campfires”) are small-scale heating events observed mainly in the solar corona. The Spectral Imaging of Coronal Environment (SPICE; SPICE Consortium et al. 2020) is a high-resolution imaging spectrometer that observes the Sun in extreme ultraviolet (EUV) wavelengths. This study uses L2 SPICE data (calibrated images)  to characterize and follow the evolution of these small-scale brightenings. Additionally, we use HRIEUV 174 Å on board Solar Orbiter and Stereo A datasets to have a complementary view of small-scale brightenings during the study period.

How to cite: Rodriguez Gomez, J. M., Young, P., and Kucera, T.: A Study of Small-Scale Brightenings using EUV Data from SPICE on board Solar Orbiter, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-1856,, 2022.

Prominences and filaments are manifestations of magnetised, levitated plasma within the solar coronal atmosphere. Expanding on our previous 2.5D work presented in Jenkins & Keppens (2021), we will present a state-of-the-art magnetohydrodynamic simulation that yields the first fully 3D model to successfully unite the extreme-ultraviolet and Hydrogen-α prominence views that contain radial striations with the equivalent on-disk filaments comprised of finite width threads. Owed to the unprecedented resolution with which this simulation is carried out, we complete a full observational synthesis and provide predictions of exactly what the instruments associated with the upcoming Solar Orbiter and DKIST will observe. We then begin with an analysis of all hydromagnetic sources of the vorticity evolution and find the internal plasma dynamics to be consistent with the nonlinear development of the magnetic Rayleigh-Taylor instability. A further stability analysis that drops the strict, idealised mRTi initial conditions then enables us to tentatively characterise the preceding linear development as the general (quasi-) interchange gravitational instability. Our simulations and analyses show clearly how this universal interchange process operates, and how our results and conclusions finally unify the contradictory prominence/filament perspectives.

How to cite: Jenkins, J. and Keppens, R.: The Role of the Magnetic Rayleigh-Taylor Instability in Resolving the Solar Prominence/Filament Paradox, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-2438,, 2022.

EGU22-2663 | Presentations | ST1.8

Extrapolation of solar wind parameters in three-dimensions in the inner heliosphere 

Aniko Timar, Andrea Opitz, Gabor Facsko, Zsuzsanna Dalya, Gergely Koban, Nikolett Biro, Akos Madar, and Zoltan Nemeth

Solar wind parameters, such as the velocity, density or pressure of the solar wind, are one of the most important factors in space physics, and their knowledge at as many points in the heliosphere as possible contributes to a broader understanding of our solar system.

Solar wind parameters at various points in the inner heliosphere are estimated using extrapolation methods. Currently, all spacecraft measuring solar wind parameters are in the ecliptic plane, thus it is enough to extrapolate the data from space probes to other spacecraft or celestial bodies near the ecliptic. Solar Orbiter, on the other hand, will soon leave the ecliptic and reach heliocentric latitudes of 34 degrees by the end of the mission, opening a new perspective.

The ballistic method extrapolates solar wind parameters in one dimension from the point of measurement to a chosen heliospheric position. The simple ballistic model considers the average rotation period of the Sun for the extrapolation in longitude while assuming a constant solar wind velocity during radial propagation. Our improved solar wind propagation model takes into account the interaction of slow and fast solar wind by applying a pressure correction during the extrapolation.

Applying this pressure-corrected ballistic method to data from solar corona models, we determined the solar wind parameters in the heliosphere in three dimensions. The advantage of our pressure-corrected ballistic method is that it is simple, it requires little calculation and it can be easily applied to the data of solar corona models in order to obtain a fast and efficient prediction in three dimensions.

How to cite: Timar, A., Opitz, A., Facsko, G., Dalya, Z., Koban, G., Biro, N., Madar, A., and Nemeth, Z.: Extrapolation of solar wind parameters in three-dimensions in the inner heliosphere, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-2663,, 2022.

EGU22-3464 | Presentations | ST1.8

The role of coronal shocks for accelerating solar energetic electrons 

Nina Dresing, Athanasious Kouloumvakos, Rami Vainio, and Alexis Rouillard

We study the role of Coronal Mass Ejection (CME)-driven shocks in the acceleration of solar energetic electrons. Using observations by the two STEREO spacecraft, we correlate electron peak intensities of solar energetic particle events measured in situ with various parameters of the associated coronal shocks. These shock parameters were derived by combining 3D shock reconstructions with global modeling of the corona. This technique provides also shock properties in the specific shock regions that are magnetically connected to the two STEREO spacecraft. We find significant correlations between the peak intensities and the Mach number of the shock with correlation coefficients of about 0.7, which are similar for electrons at ∼ 1 MeV and protons at > 60 MeV. Lower energy electrons < 100 keV show a smaller correlation coefficient of 0.47. The causal relationship between electron intensities and the shock properties is supported by the vanishing correlations, when peak intensities at STEREO-A are related with the Alfvénic Mach numbers at the magnetic footpoint of STEREO-B and vice versa, which yields correlation coefficients of 0.03 and -0.13 for ∼ 1 MeV and < 100 keV electron peak intensities, respectively. We conclude that the high-energy electrons are accelerated mainly by the shock, while the low energy electrons are likely produced by a mixture of flare and shock-related acceleration processes.

How to cite: Dresing, N., Kouloumvakos, A., Vainio, R., and Rouillard, A.: The role of coronal shocks for accelerating solar energetic electrons, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-3464,, 2022.

EGU22-4121 | Presentations | ST1.8

Spatial Structure of CIRs 

Gergely Koban, Andrea Opitz, Zoltan Nemeth, Gabor Facsko, Akos Madar, Aniko Timar, Zsuzsanna Dalya, and Nikolett Biro

Co-rotating Interaction Regions are complex and fascianting structures in the Heliosphere that 
play an important role in space weather. They arise from the fast solar wind interacting with the 
slow solar wind streams. The interface between fast and slow solar wind is called the stream 
interface, and it is common for CIRs to produce forward shock at the leading edge and reverse 
shock at the trailing edge. CIRs often have considerable tilts in the north-south axis, owing to the magnetic 
conditions on the Sun.

Examination of the spatial structure of CIRs, – most importantly the aforementioned tilt – is not 
an easy task. We attempt a multi-spacecraft investigation in order to examine the spatial 
structure of CIRs on different distance scales. Using all available spacecraft data nearby, the tilt 
of the stream interface can be determined considering the time delays of the effects caused by 
the CIR recorded by each spacecraft. Our final aim is to improve solar wind propagation 
methods with these detailed CIR results.

How to cite: Koban, G., Opitz, A., Nemeth, Z., Facsko, G., Madar, A., Timar, A., Dalya, Z., and Biro, N.: Spatial Structure of CIRs, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-4121,, 2022.

EGU22-4155 | Presentations | ST1.8

Spatial variation of the background solar wind in the Inner Heliosphere 

Nikolett Biro, Andrea Opitz, Zoltan Nemeth, Aniko Timar, Akos Madar, Gergely Koban, and Zsuzsanna Dalya
The importance of background solar wind is unquestionable as it carries information on the solar surface conditions and has a major role in space weather events. The current solar minimum is a perfect time period for investigations regarding this field, with several space probes providing in-situ measurements.
Our aim is to determine the spatial variations in the background solar wind through multi-spacecraft data analysis, including recent missions, such as Parker Solar Probe and Solar Orbiter. We adjust for the radial and longitudinal time-lags between the different spacecraft, then compare their solar wind plasma measurements. The effects of latitudinal differences between the observations is then backmapped to coronagraph imagery. The results will be useful for further analysis of inner heliospheric structures, for the improvement of propagation models, and to support the analysis of out-of-ecliptic solar wind observations.

How to cite: Biro, N., Opitz, A., Nemeth, Z., Timar, A., Madar, A., Koban, G., and Dalya, Z.: Spatial variation of the background solar wind in the Inner Heliosphere, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-4155,, 2022.

EGU22-4178 | Presentations | ST1.8

Predicting the solar cycle amplitude with the new catalogue of hemispheric sunspot numbers 

Shantanu Jain, Tatiana Podladchikova, Astrid M. Veronig, Olga Sutyrina, Mateja Dumbović, Frédéric Clette, and Werner Pötzi

The sun’s magnetic field drives the 11-year solar cycle, and predicting its strength has practical importance for many space weather applications. Previous studies have shown that analysing the solar activity of the two hemispheres separately instead of the full sun can provide more detailed information on the activity evolution. However, the existing Hemispheric Sunspot Number data series (1945 onwards) was too short for meaningful solar cycle predictions. Based on a newly created hemispheric sunspot number catalogue for the time range 1874-2020 (Veronig et al. 2021, that is compatible with the International Sunspot Number from World Data Centre SILSO, we investigate the evolution of the solar cycle for the two hemispheres and develop a novel method for predicting the solar cycle amplitude. We demonstrate a steady relationship between the maximal growth rate of activity in the ascending phase of a cycle and its subsequent amplitude and form a 3rd order regression for the predictions. Testing this method for cycles 12-24, we show that the forecast made by the sum of the maximal growth rate from the North and South Hemispheric Sunspot number is more accurate than the same forecast from the Total Sunspot Number: The rms error of predictions is smaller by 27%, the correlation coefficient r is higher by 11% on average reaching values in the range r = 0.8-0.9 depending of the smoothing window of the monthly mean data. These findings demonstrate that empirical solar cycle prediction methods can be enhanced by investigating the solar cycle dynamics in terms of the hemispheric sunspot numbers, which is a strong argument supporting regular monitoring, recording, and analysing solar activity separately for the two hemispheres.

How to cite: Jain, S., Podladchikova, T., Veronig, A. M., Sutyrina, O., Dumbović, M., Clette, F., and Pötzi, W.: Predicting the solar cycle amplitude with the new catalogue of hemispheric sunspot numbers, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-4178,, 2022.

EGU22-4182 | Presentations | ST1.8

The large-scale structure of the solar wind: flux conservation, radial scalings, Mach numbers, and critical distances 

Daniel Verscharen, Stuart D. Bale, and Marco Velli

One of the key challenges in solar and heliospheric physics is to understand the acceleration of the solar wind. As a super-sonic, super-Alfvénic plasma flow, the solar wind carries mass, momentum, energy, and angular momentum from the Sun into interplanetary space. We present a framework  based on two-fluid magnetohydrodynamics to estimate the flux of these quantities based on spacecraft data independent of the heliocentric distance of the location of measurement.

Applying this method to the Ulysses dataset allows us to study the dependence of these fluxes on heliolatitude and solar cycle. The use of scaling laws provides us with the heliolatitudinal dependence and the solar-cycle dependence of the scaled Alfvénic and sonic Mach numbers as well as the Alfvén and sonic critical radii. Moreover, we estimate the distance at which the local thermal pressure and the local energy density in the magnetic field balance.

These results serve as predictions for observations with Parker Solar Probe, which currently explores the very inner heliosphere, and Solar Orbiter, which will measure the solar wind outside the plane of the ecliptic in the inner heliosphere during the course of the mission.

How to cite: Verscharen, D., Bale, S. D., and Velli, M.: The large-scale structure of the solar wind: flux conservation, radial scalings, Mach numbers, and critical distances, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-4182,, 2022.

EGU22-4327 | Presentations | ST1.8 | Highlight

First results from the SO/PHI instrument on Solar Orbiter 

Sami Khan Solanki, Julian Blanco, Fatima Kahil, Philipp Loeschl, Hanna Strecker, Johann Hirzberger, David Orozco Suárez, Jose Carlos Del Toro Iniesta, Joachim Woch, Achim Gandorfer, Alberto Alvarez-Herrero, Thierry Appourchaux, and Reiner Volkmer

The ESA/NASA Solar Orbiter mission, launched in February 2020, has just completed its cruise phase. During its nominal mission, it will explore the Sun and heliosphere from close up and from out of the ecliptic plane. It aims to address the overarching questions of how the Sun creates and controls the heliosphere, and why solar activity changes with time. Among the instruments onboard Solar Orbiter  is the Polarimetric and Helioseismic Imager (SO/PHI), which is the first magnetograph to leave the Sun-Earth line and to observe the Sun from different directions. Already during the cruise phase of Solar Orbiter, SO/PHI has provided a few glimpses of its capabilities, including the excellent quality of the data. In spite of the very limited amount of data gathered during cruise, a few interesting results have already been obtained. A selection of such results will be presented.


How to cite: Solanki, S. K., Blanco, J., Kahil, F., Loeschl, P., Strecker, H., Hirzberger, J., Orozco Suárez, D., Del Toro Iniesta, J. C., Woch, J., Gandorfer, A., Alvarez-Herrero, A., Appourchaux, T., and Volkmer, R.: First results from the SO/PHI instrument on Solar Orbiter, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-4327,, 2022.

EGU22-4987 | Presentations | ST1.8

Removal of false alarms from solar wind predictions 

Zsuzsanna Dalya and Andrea Opitz

Solar wind propagation models using in situ plasma observations as input can be improved by removing the signatures of Interplanetary Coronal Mass Ejections (ICMEs) from the input data. ICMEs are sporadic events that propagate in a given direction, hence their signatures in the plasma data should not be extrapolated to other heliocentric longitudes. As a result, in order to improve the prediction accuracy these should be filtered out. We create dedicated ICME lists providing the exact start and end times of ICMEs to numerous space probes such as ACE, STEREO A&B, SOHO, WIND, SolO, PSP, DSCOVER, VEX, MEX, Rosetta, BepiColombo, MAVEN and Messenger. We provide solar wind plasma and magnetic field predictions to any inner heliospheric position applying the ICME filter on the input dataset this way eliminating false alarms such as false ICME signature forecasts. Our corrected predictions contribute to the investigation of the background solar wind and related fields.

How to cite: Dalya, Z. and Opitz, A.: Removal of false alarms from solar wind predictions, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-4987,, 2022.

EGU22-5443 | Presentations | ST1.8

Large Amplitude Oscillations in Solar Filaments Caused by Solar Jets 

Reetika Joshi, Manuel Luna, Brigitte Schmieder, Fernando Moreno-Insertis, and Ramesh Chandra

Large amplitude oscillations (LAO) are often detected in filaments. Their origin has been associated with shock waves or to interaction with eruptions and jets. In this study we present  two examples of LAOs due to solar jets interaction with filaments on February 3-5 2015 and March 14 2015.  The filament eruption on March 14 was followed by a two step filament eruption along with a CME and become the strong geomagnetic storm of Solar Cycle 24 on 17 March 2015. These LAOs are  analysed by using time-distance diagnostics. The detected LAOs have periods of around 60 minutes and are damped after three oscillations. The observations are consistent with the results of a recent developed theoretical model of jet and filament interaction. The jets are associated with very weak flares which did not initiate any EUV wave. The role of waves as trigger of these oscillations can be discarded for these two events. 


How to cite: Joshi, R., Luna, M., Schmieder, B., Moreno-Insertis, F., and Chandra, R.: Large Amplitude Oscillations in Solar Filaments Caused by Solar Jets, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-5443,, 2022.

EGU22-5466 | Presentations | ST1.8

Temporal evolution of the background solar wind throughout the inner heliosphere 

Andrea Opitz, Aniko Timar, Zoltan Nemeth, Zsuzsanna Dalya, Gergely Koban, Nikolett Biro, and Akos Madar

The solar wind properties at a given point in the heliosphere depend strongly on the source surface characteristics, the dynamical effects during propagation and the transient events. We study the background solar wind structures after modelling their propagation throughout the 3-dimensional heliosphere. We remove the transient events from the observations, then apply the ballistic propagation method corrected for pressure gradients at stream interactions. A detailed multi-spacecraft investigation of the radial and latitudinal effects improves our model. These results are applied to study the temporal evolution of the solar wind by excluding the spatial effects through adjusting for the timelag calculated from the spacecraft separations.

How to cite: Opitz, A., Timar, A., Nemeth, Z., Dalya, Z., Koban, G., Biro, N., and Madar, A.: Temporal evolution of the background solar wind throughout the inner heliosphere, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-5466,, 2022.

EGU22-6088 | Presentations | ST1.8

ULF waves upstream and downstream of interplanetary shocks 

Xochitl Blanco-Cano, Primoz Kajdic, Diana Rojas-Castillo, Luis Preisser, Lan Jian, Christopher Russell, and Janet Luhmann

Interplanetary (IP) shocks can be driven in the solar wind by fast coronal mass ejections and by the interaction of fast solar wind with slow streams of plasma. These shocks can be preceded by extended wave and suprathermal ion foreshocks. We use STEREO data to study wave modes upstream and downstream of IP shocks. Understanding these waves is important because they contribute to shock acceleration processes and modify the solar wind as the shocks propagate in the heliosphere. We find that upstream regions can be permeated by whistler waves (f ~ 1 Hz) and/or ultra low frequency (ULF) right-handed waves (f~10-2–10-1 Hz). While whistlers appear to be generated at the shock, the origin of ULF waves is most probably associated with local kinetic ion instabilities. In contrast with planetary bow shocks, most IP shocks have a small Mach number (<4) and most of the upstream waves studied here are mainly transverse and no steepening occurs. Downstream of the shocks ion cyclotron and mirror mode waves can grow due to temperature anisotropy. The waves observed downstream of IP quasi-parallel shocks have larger amplitudes than waves in the regions downstream of quasi-perpendicular shocks. A variety of waves can be found in the sheath regions of IP shocks, even when IP shocks are weak, mostly for quasi-perpendicular shocks. These include ion cyclotron waves (ICW) with well defined peaks in frequency, broad spectra ICW, and mirror mode storms, which tend to occur for higher plasma beta.

How to cite: Blanco-Cano, X., Kajdic, P., Rojas-Castillo, D., Preisser, L., Jian, L., Russell, C., and Luhmann, J.: ULF waves upstream and downstream of interplanetary shocks, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-6088,, 2022.

EGU22-6620 | Presentations | ST1.8

Directional Discontinuities in the Inner Heliosphere 

Ákos Madár, Geza Erdos, Andrea Opitz, Zoltan Nemeth, Gabor Facsko, Aniko Timar, Nikolett Biro, Gergely Koban, and Zsuzsa Dalya

The Parker Solar Probe and Solar Orbiter spacecraft make whole new spatial and time scales available in the inner Heliosphere. With these new data, we study directional discontinuities that are common structures in the solar wind in this region. Their radial distribution can provide insight into the physical processes of this virtually collisionless plasma. Applying a method (Erdős & Balogh, 2008) based on minimum variance analysis to select directional discontinuities in magnetic field data, we determine their number as a function of distance from the Sun. How the directional discontinuity occurrence rate depends on the solar wind speed is also part of our analysis.

How to cite: Madár, Á., Erdos, G., Opitz, A., Nemeth, Z., Facsko, G., Timar, A., Biro, N., Koban, G., and Dalya, Z.: Directional Discontinuities in the Inner Heliosphere, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-6620,, 2022.

EGU22-6726 | Presentations | ST1.8

Multi-Spacecraft Observations of Gradual Solar Energetic Particle Events with Enhanced 3He Abundance 

Radoslav Bucik, Glenn M. Mason, Raúl Gómez-Herrero, Maher A. Dayeh, Mihir I. Desai, Samuel T. Hart, George C. Ho, David Lario, Vratislav Krupar, Robert F. Wimmer-Schweingruber, Javier Rodríguez-Pacheco, and Tilaye T. Asfaw

Suprathermal ions from coronal jets, characterized by enhanced 3He and heavy-ion abundances, are an essential component of the seed population accelerated by coronal mass ejection (CME)-driven shocks in gradual solar energetic particle (GSEP) events. However, the mechanisms through which CME-driven shocks gain access to these suprathermal ions and produce spectral and abundance variations in GSEP events remain largely unexplored. We study GSEP events simultaneously measured on at least two spacecraft, such as ACE, STEREO, and Solar Orbiter, where 3He finite mass peak is measured at least on one spacecraft. This presentation discusses the origin of vastly different abundances and spectral shapes in terms of variable remnant population from preceding impulsive SEP events. Furthermore, with the help of imaging observations from SDO and STEREO, we examine a possible direct contribution from parent active regions of GSEP events.

How to cite: Bucik, R., Mason, G. M., Gómez-Herrero, R., Dayeh, M. A., Desai, M. I., Hart, S. T., Ho, G. C., Lario, D., Krupar, V., Wimmer-Schweingruber, R. F., Rodríguez-Pacheco, J., and Asfaw, T. T.: Multi-Spacecraft Observations of Gradual Solar Energetic Particle Events with Enhanced 3He Abundance, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-6726,, 2022.

EGU22-6802 | Presentations | ST1.8 | Highlight

Interstellar Probe: A Mission to the Heliospheric Boundary and Interstellar Medium to Understand our Home in the Galaxy 

Pontus Brandt and the The Interstellar Probe Study Team

For the past 60, 000 years our Sun and its protective heliosphere have been plowing through the Local Interstellar Cloud (LIC), but is now in a historic transition region towards the G-cloud that could have dramatic consequences for the global heliospheric structure. An Interstellar Probe mission to the Very Local Interstellar Medium (VLISM) would bring new scientific discoveries of the mechanisms upholding our vast heliosphere and directly sample the Local Interstellar Clouds to allow us, not only to understand the current dynamics and shielding, but also how the heliosphere responded in the past and how it will respond in the new interstellar environment. An international team of scientists and experts have now completed a NASA-funded study led by The Johns Hopkins University Applied Physics Laboratory (APL) to develop pragmatic example mission concepts for an Interstellar Probe with a nominal design lifetime of 50 years. The team has analyzed dozens of launch configurations and demonstrated that asymptotic speeds in excess of 7.5 Astronomical Units (AU) per year can be achieved using existing or near-term propulsion stages with a powered or passive Jupiter Gravity Assist (JGA). These speeds are more than twice that of the fastest escaping man-made spacecraft to date, which is Voyager 1 currently at 3.59 AU/year. An Interstellar Probe would therefore reach the Termination Shock (TS) in less than 12 years and cross the Heliopause into the VLISM after about 16 years from launch.

In this presentation we provide an overview of the study, the science mission concept, discuss the compelling discoveries that await, and the associated example science payload, measurements and operations ensuring a historic data return that would push the boundaries of space exploration by going where no one has gone before.

How to cite: Brandt, P. and the The Interstellar Probe Study Team: Interstellar Probe: A Mission to the Heliospheric Boundary and Interstellar Medium to Understand our Home in the Galaxy, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-6802,, 2022.

EGU22-7831 | Presentations | ST1.8

Langmuir waves associated with magnetic holes in the solar wind 

Joan Jordi Boldu-O´Farrill Treviño, Daniel Graham, Michiko Morooka, Tomas Karlsson, Yuri Khotyaintsev, Mats André, Jan Souček, and Milan Maksimovic

Langmuir waves (electrostatic waves near the electron plasma frequency) are often observed in the solar wind, playing an important role in the energy dissipation of electrons. The largest amplitude waves are typically associated with type II and III solar radio bursts and planetary foreshocks. However, Langmuir waves not connected with radio bursts are also found in the solar wind. The causes of these Langmuir waves are not well understood. Langmuir waves are also found around magnetic holes, a localised depression of the magnetic field strength. This study aims to investigate the relationship between Langmuir waves and magnetic holes in the solar wind using electric and magnetic field measurements performed by the Solar Orbiter’s RPW and MAG instruments during 2020 and 2021. We identified a large set of Langmuir wave events from the RPW/TDS (Time Domain Sampler) waveform data using the plasma density estimated from the spacecraft’s potential obtained by RPW, showing that ~7% of them have been spotted inside magnetic holes. We will compare these events with local plasma conditions analysing the electron distribution functions, and discuss the mechanisms that may lead to the generation of Langmuir waves associated with magnetic holes in the solar wind.

How to cite: Boldu-O´Farrill Treviño, J. J., Graham, D., Morooka, M., Karlsson, T., Khotyaintsev, Y., André, M., Souček, J., and Maksimovic, M.: Langmuir waves associated with magnetic holes in the solar wind, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-7831,, 2022.

EGU22-7840 | Presentations | ST1.8

The exoplanetary magnetosphere extension in Sun-like stars based on the solar wind and solar UV emission 

Raffaele Reda, Luca Giovannelli, Tommaso Alberti, Francesco Berrilli, Luca Bertello, Dario Del Moro, Maria Pia Di Mauro, Piermarco Giobbi, and Valentina Penza

The solar activity in form of coronal mass ejections or solar wind disturbances, such as slow or high speed streams, affects the circumterrestrial electromagnetic environment, with a primary effect on the magnetosphere, compressing and perturbing it. Here, in order to connect the long-term solar activity variations to solar wind properties, we use measurement of a proxy for chromospheric activity, the Ca II K index, and solar wind OMNI data for the time interval 1965-2021, which almost entirely covers the last 5 solar cycles. By using both a cross correlation and a mutual information approach, a 3.6-year mean lag has been found between Ca II K index and solar wind dynamic pressure. This result allows us to obtain a relationship between the solar UV emission and the solar wind dynamic pressure, enabling us to derive the Earth’s magnetospheric extension over the last 5 solar cycles.
Moreover, the advantage of having used the Ca II K index proxy is that the relation found for the Sun can be easily extended to other stars with similar properties (i.e. Sun-like stars). To this scope, the model is then used to study the effect of stellar wind dynamic pressure on the magnetosphere of Earth-like planets orbiting at 1 AU around a sample of Sun-like stars.

How to cite: Reda, R., Giovannelli, L., Alberti, T., Berrilli, F., Bertello, L., Del Moro, D., Di Mauro, M. P., Giobbi, P., and Penza, V.: The exoplanetary magnetosphere extension in Sun-like stars based on the solar wind and solar UV emission, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-7840,, 2022.

EGU22-8435 | Presentations | ST1.8

Studies on Space Climate Made in the University of Extremadura (Spain) 

Víctor M.S. Carrasco, Alejandro J.P. Aparicio, José M. Nogales, Nieves Bravo-Paredes, Irene Tovar, Maricruz Gallego, and José M. Vaquero

Our research group has carried out several studies on space climate in the past. One of our lines of work is the recovery and analysis of the catalogues including historical sunspot observations. We have already published in digital version some sunspot catalogues made in observatories, for example, of the Iberian Peninsula (Aparicio et al. 2018). Recently, we have also digitized the catalogue from the sunspot observations made in the Stonyhurst College Observatory for 1921 – 1935 (Carrasco et al. 2021) and completed the first step for the publication of the Sacramento Peak Observatory sunspot catalogue (1947 – 2004) (Carrasco et al. 2020). Currently, we are analyzing the sunspot observations from drawings made in Boulder (U.S.A) for the period 1966 – 1992 to create a sunspot group catalogue including that information. Furthermore, we have also analyzed long-term observation series made by individual astronomers. Two examples of this kind of studies can be the analysis of the sunspot observations made by David Hadden during the period 1890 – 1931 (Carrasco et al. 2013) and those made by Eric Strach for 1969 – 2008 (Carrasco et al. 2019). Nowadays, in collaboration with other Italian research group, we are studying all the sunspot drawings made by Father Angelo Secchi in the second half of the 19th century. We are constructing the Wolf number, group number and the area series from these drawings. As future work, our objective is to publish a sunspot group catalogue from these observations.


Aparicio, A.J.P., Lefèvre, L., Gallego, M.C., Vaquero, J.M., Clette, F., Bravo-Paredes, N., Galaviz, P., Bautista, M.L.: 2018, A Sunspot Catalog for the Period 1952 – 1986 from Observations Made at the Madrid Astronomical Observatory, Solar Physics 293, 164.

Carrasco, V.M.S., Vaquero, J.M., Gallego, M.C., Trigo, R.M.: 2013, Forty two years counting spots: Solar observations by D.E. Hadden during 1890–1931 revisited. New Astronomy 25, 95.

Carrasco, V.M.S, Vaquero, J.M., Olmo-Mateos, V.M.: 2019, Eric Strach: Four Decades of Detailed Synoptic Solar Observations (1969‐2008), Space Weather 17, 796.

Carrasco, V.M.S., Pevtsov, A.A., Nogales, J.M., Vaquero, J.M.: 2020, The Sunspot Drawing Collection of the National Solar Observatory at Sacramento Peak (1947–2004), Solar Physics 296, 3.

Carrasco, V.M.S., Muñoz-Jaramillo, A., Nogales, J.M., Gallego, M.C., Vaquero, J.M.: 2021, Sunspot Catalog (1921–1935) and Area Series (1886–1940) from the Stonyhurst College Observatory, Astrophysical Journal Supplement Series 256, 38.

How to cite: Carrasco, V. M. S., Aparicio, A. J. P., Nogales, J. M., Bravo-Paredes, N., Tovar, I., Gallego, M., and Vaquero, J. M.: Studies on Space Climate Made in the University of Extremadura (Spain), EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-8435,, 2022.

With current observational methods it is not possible to determine the magnetic field in the solar corona accurately. Therefore, coronal magnetic field models have to rely on extrapolation methods using photospheric magnetograms as boundary conditions. In recent years, due to the increased resolution of observations and the need to resolve non-force-free lower regions of the solar atmosphere, there have been increased efforts to use magnetohydrostatic (MHS) field models instead of force-free extrapolation methods. Although numerical methods to calculate MHS solutions can deal with non-linear problems and hence provide more accurate models, analytical three-dimensional MHS equilibria can also be used as a numerically relatively “cheap” complementary method.


We discuss a family of analytical MHS equilibria that allows for a transition from a non-force-free region to a force-free region. The solution involves hypergeometric functions and while routines for the calculation of these are available, this can affect both the speed and the numerical accuracy of the calculations. We therefore look into the asymptotic behaviour of this solution in order to numerically approximate it through exponential functions aiming to improve the numerical efficiency. We present an illustrative example by comparing field line profiles, density and pressure differences between the exact solutions, the asymptotic solution and a hybrid model where the use of the hypergeometric function is restricted to an area around the transitional region between the non-force-free and the force-free domain.

How to cite: Nadol, L. and Neukirch, T.: Coronal Magnetic Field Extrapolation Using a Specific Family of Analytical 3D Magnetohydrostatic Equilibria — Practical Aspects, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-9430,, 2022.

EGU22-9622 | Presentations | ST1.8

Clarifying the relation between chromospheric emissions and photospheric magnetic fields using AIA 1600 Å and HMI data 

Ismo Tähtinen, Ilpo Virtanen, Alexei Pevtsov, and Kalevi Mursula

Intense photospheric magnetic fields manifest as enhanced emission in several spectral lines and parts of the continuum. Here we aim to improve the understanding of the relation between magnetic fields and radiative structures by using the seeing-free observations of the Atmospheric Imaging Assembly (AIA) and the Helioseismic Magnetic Imager (HMI), both on-board Solar Dynamics Observatory (SDO). We use AIA 1600 Å band, which captures the far-ultraviolet continuum emission originating from the temperature minimum of the solar chromosphere.

We developed a novel, objective method to define thresholds separating the brightest AIA 1600 Å pixels ("bright pixels") and the least bright pixels ("dark pixels") from the AIA 1600 Å brightness distribution. According to the method bright pixels are pixels whose standardized contrast (ratio of brightness to the center-to-limb variation) exceeds the level of I = 1.93. This threshold maximizes the average size of bright clusters (4-connected regions of bright pixels). Dark pixels are pixels whose standardized contrast is below I = 0.5. This corresponds to a threshold below which there are practically no pixels on quiet days. Comparing the AIA 1600 Å intensity and HMI magnetic field observations, we found that the AIA 1600 Å dark pixels correspond to the strongest magnetic field (B > 1325 G) pixels. These pixels are typically within sunspots. On the other hand, we found the AIA bright pixels correspond to moderate (55 G < B < 475 G) magnetic field intensity pixels of HMI.

We found that the percentage of AIA bright pixels on the solar surface almost entirely explains the observed variability of the total AIA 1600 Å emission, even in the presence of large sunspot groups. We developed a multi-linear regression model, which reliably predicts the magnitude of the disk-averaged unsigned magnetic field using the measured percentages of bright and dark pixels. We found that the large bright clusters have a constant mean unsigned magnetic field, similarly as earlier found for, e.g., Ca II K plages. However, the magnetic field strength of bright clusters is 246.5 ± 0.1 G, which is roughly 100 G larger than found earlier for Ca II K plages.

How to cite: Tähtinen, I., Virtanen, I., Pevtsov, A., and Mursula, K.: Clarifying the relation between chromospheric emissions and photospheric magnetic fields using AIA 1600 Å and HMI data, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-9622,, 2022.

EGU22-10969 | Presentations | ST1.8

Evolution and variation of the heliospheric current sheet during 1995-2009 

Kan Liou and Chin-Chun Wu

The heliospheric current sheet (HCS) is the largest known solar wind structure that exists persistently and continuously within the heliosphere. While it is discovered for more than half a century ago, owing to very limited and scattered in-situ solar wind observations in the heliosphere, the shape and evolution of the HCS are still little known and constitute the key subject in the study of global-scale solar wind. Currently the morphology of the HCS is derived largely from simple kinematic approaches that map the neutral sheet observed at 2.5 solar radii outward into the heliosphere. The dynamical effect of solar wind interactions on the evolution and global structure of the HCS is still poorly understood. Here we present results from a study of the HCS from 1995 to 2009 using time-dependent, three-dimensional, global magnetohydrodynamic (MHD) model simulations. We focus on the radial and solar cycle variations of the HCS tilt angle from 18 solar radii to 7 AU. We also compare our result with those obtained from ballistic outward projections of the neutral sheet observed at 2.5 solar radii. We present our analysis result and discuss possible implications.

How to cite: Liou, K. and Wu, C.-C.: Evolution and variation of the heliospheric current sheet during 1995-2009, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-10969,, 2022.

EGU22-11163 | Presentations | ST1.8

Deriving proper electron characteristics in the solar wind from Solar Orbiter observations 

Štěpán Štverák, Georgios Nicolaou, Christopher J. Owen, Milan Maksimovic, and Pavel M. Trávníček

Solar Orbiter, the latest ESA solar and heliospheric space mission, provides the most recent plasma measurements in the free streaming solar wind across a wide range of the radial distance from the Sun. Electron observations are enabled by the Electron Analyser System (EAS) being part of the Solar Wind Analyser (SWA) suit of instruments. Electron properties are measured in a form of velocity distribution functions (VDFs) with an unprecedent sampling rate of 10 s in the normal operational mode going down to even 8 Hz for limited burst mode snapshots. Regular EAS observations started in mid of 2020 and all data are being continuously made publicly available throughout the ESA's Solar Orbiter Archive (SOAR). Here we present some preliminary user-based methods and techniques of deriving proper electron characteristics from the Level 2 calibrated data set by calculating the required moments of measured VDFs. In particular, we focus on data correction with respect to possible instrument and/or spacecraft related effects. The obtained results are compared and validated by a cross calibration with respect to electron properties derived from the quasi-thermal noise measurements provided by the Radio and Plasma Waves (RPW) instrument, being also part of the in situ plasma payload of the Solar Orbiter mission.

How to cite: Štverák, Š., Nicolaou, G., Owen, C. J., Maksimovic, M., and Trávníček, P. M.: Deriving proper electron characteristics in the solar wind from Solar Orbiter observations, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-11163,, 2022.

EGU22-11326 | Presentations | ST1.8

Dynamics of helium abundance inside and around ICME 

Alexander Khokhlachev, Yuri Yermolaev, Maria Riazantseva, Liudmila Rakhmanova, and Irina Lodkina

The interplanetary manifestations of coronal mass ejections (ICME) typically characterized by an increased relative abundance of doubly ionized helium, which can be several time higher comparative to the slow solar wind streams and exceed 10%. At the same time the helium abundance can dynamically change as a result of local processes in solar wind plasma. However, the details of helium abundance dynamic have not yet been sufficiently studied.

The study is devoted to the changing of the helium abundance inside the ICME including both its types - Magnetic Clouds (MC) and EJECTA, and also in the vicinity of these large scale structures. We consider the helium abundance dynamic, and its relation with other solar wind parameters both at large (>106 km) and medium (104-105 km) scales, based on the statistic and correlation analysis of plasma and magnetic field parameters [Yermolaev et al., 2020; Khokhlachev et al., 2021]. Hourly average data from OMNI-2 database and 3 second WIND spacecraft measurements are used for the analysis of corresponding scales. Selection of the intervals related to ICME events (MC, EJECTA and SHEATH regions) was produced with the help of the catalog of large-scale events of Space Research Institute ( [Yermolaev et al., 2009]). It is shown that on large scales the helium abundance generally increases with a decrease in the plasma parameter β in the ICME. This relation can be predominantly explained by the strong positive correlation of helium abundance with magnetic pressure, while correlation with thermal pressure is ambiguous: weak negative for Magnetic Clouds and weak positive for EJECTA. This confirms the previously proposed hypothesis of the ion current flow enriched by helium ions inside the ICME [Yermolaev et al., 2020]. On medium scales, the trend of anticorrelation of the helium abundance with plasma β-parameter in the ICME is observed also. However, the dependences of the helium abundance on the plasma and magnetic field parameters can dynamically change at smaller scales. For example, several local medium-scale structures with a negative correlation of helium abundance and magnetic field magnitude can be observed on the background of positive correlation at large-scale structures.


Yermolaev Y.I. et al., Catalog of large-scale solar wind phenomena during 1976–2000, Cosmic Res., 2009, vol. 47, no. 2, pp. 81–94

Yermolaev Y.I. et al., Dynamics of large-scale solar-wind streams obtained by the double superposed epoch analysis. 4. Helium abundance, Journal of Geophysical Research, 2020, 125 (7)

Khokhlachev A. A. et al., Variations of Protons and Doubly Ionized Helium Ions in the Solar Wind, Cosmic Research, 2021, vol. 59, no. 6, pp. 415-426

How to cite: Khokhlachev, A., Yermolaev, Y., Riazantseva, M., Rakhmanova, L., and Lodkina, I.: Dynamics of helium abundance inside and around ICME, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-11326,, 2022.

EGU22-11910 | Presentations | ST1.8

Radial evolution of the key solar wind parameters according to the Parker Solar Probe and other missions 

Uliana Antsiferova and Olga Khabarova

Knowing the temporal and spatial behavior of the main plasma parameters in the solar wind is important because of both an academic interest and the necessity to build theoretical models. Meanwhile, the solar wind characteristics are available with good accuracy only from in situspacecraft observations. There are a very limited number of studies that analyzed the radial evolution of the interplanetary magnetic field, the speed, the density, and the temperature of the solar wind. Previously, data from missions carried out in the 1970s were used for these purposes. Meanwhile, none of them approached the Sun closer than ~0.3 AU, and the information available did not allow describing the processes occurring near the solar corona. The launch of the Parker Solar Probe mission has opened up vast opportunities for studying the solar wind, starting from distances of the order of the Alfven radius. At the moment, Parker Solar Probe data are available with approaches to the Sun up to 0.08 AU.

An analysis of data obtained from Parker Solar Probe and Helios 2 (up to 1 AU), and from IMP8 and Voyager 1 (further from the Earth's orbit and up to 7 AU) was performed. The dependencies of the interplanetary magnetic field, the temperature, the speed, and the density of the solar wind on heliocentric distance are revealed and their typical profiles are analyzed. The nature of deviations of observed values from theoretical expectations is discussed. A particular attention is paid to the comparison of the results from the Parker Solar Probe and Helios missions.

How to cite: Antsiferova, U. and Khabarova, O.: Radial evolution of the key solar wind parameters according to the Parker Solar Probe and other missions, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-11910,, 2022.