-
How coronal mass ejections are influenced by the morphology and toroidal flux of their source magnetic flux ropes?
Authors:
J. H. Guo,
L. Linan,
S. Poedts,
Y. Guo,
B. Schmieder,
A. Lani,
Y. W. Ni,
M. Brchnelova,
B. Perri,
T. Baratashvili,
S. T. Li,
P. F. Chen
Abstract:
Coronal mass ejections (CMEs) stand as intense eruptions of magnetized plasma from the Sun, playing a pivotal role in driving significant changes of the heliospheric environment. Deducing the properties of CMEs from their progenitors in solar source regions is crucial for space weather forecasting. Deducing the properties of CMEs from their progenitors in solar source regions is crucial for space…
▽ More
Coronal mass ejections (CMEs) stand as intense eruptions of magnetized plasma from the Sun, playing a pivotal role in driving significant changes of the heliospheric environment. Deducing the properties of CMEs from their progenitors in solar source regions is crucial for space weather forecasting. Deducing the properties of CMEs from their progenitors in solar source regions is crucial for space weather forecasting. The primary objective of this paper is to establish a connection between CMEs and their progenitors in solar source regions, enabling us to infer the magnetic structures of CMEs before their full development. To this end, we create a dataset comprising a magnetic flux rope series with varying projection shapes, sizes and toroidal fluxes, using the Regularized Biot-Savart Laws (RBSL). Thereafter, we simulate the propagation of these flux ropes from the solar surface to a distance of 25$R_{\odot}$ with our global coronal MHD model which is named COCONUT. Our parametric survey reveals significant impacts of source flux ropes on the consequent CMEs. We find that the projection shape can influence the magnetic structures of CMEs at 20$R_{\odot}$, albeit with minimal impacts on the propagation speed. However, these impacts diminish as source flux ropes become fat. In terms of toroidal flux, our simulation results demonstrate a pronounced correlation with the propagation speed of CMEs, as well as the successfulness in erupting. This work builds the bridge between the CMEs in the outer corona and their progenitors in solar source regions. Our parametric survey suggests that the projection shape, cross-section radius and toroidal flux of source flux ropes are crucial parameters in predicting magnetic structures and propagation speed of CMEs, providing valuable insights for space weather prediction.
△ Less
Submitted 12 July, 2024;
originally announced July 2024.
-
Modelling the propagation of coronal mass ejections with COCONUT: implementation of the Regularized Biot-Savart Laws flux rope model
Authors:
Jinhan Guo,
L. Linan,
S. Poedts,
Y. Guo,
A. Lani,
B. Schmieder,
M. Brchnelova,
B. Perri,
T. Baratashvili,
Y. W. Ni,
P. F. Chen
Abstract:
Context: Coronal mass ejections (CMEs) are rapid eruptions of magnetized plasma that occur on the Sun, which are known as the main drivers of adverse space weather. Accurately tracking their evolution in the heliosphere in numerical models is of utmost importance for space weather forecasting. Aims: The main objective of this paper is to implement the Regularized Biot-Savart Laws (RBSL) method in…
▽ More
Context: Coronal mass ejections (CMEs) are rapid eruptions of magnetized plasma that occur on the Sun, which are known as the main drivers of adverse space weather. Accurately tracking their evolution in the heliosphere in numerical models is of utmost importance for space weather forecasting. Aims: The main objective of this paper is to implement the Regularized Biot-Savart Laws (RBSL) method in a new global corona model COCONUT. This approach has the capability to construct the magnetic flux rope with an axis of arbitrary shape. Methods: We present the implementation process of the RBSL flux rope model in COCONUT, which is superposed onto a realistic solar wind reconstructed from the observed magnetogram around the minimum of solar activity. Based on this, we simulate the propagation of an S-shaped flux rope from the solar surface to a distance of 25 solar radii. Results: Our simulation successfully reproduces the birth process of a CME originating from a sigmoid in a self-consistent way. The model effectively captures various physical processes and retrieves the prominent features of the CMEs in observations. In addition, the simulation results indicate that the magnetic topology of the CME flux rope at around 20 solar radii deviates from a coherent structure, and manifests as a mix of open and closed field lines with diverse footpoints. Conclusions: This work demonstrates the potential of the RBSL flux rope model in reproducing CME events that are more consistent with observations. Moreover, our findings strongly suggest that magnetic reconnection during the CME propagation plays a critical role in destroying the coherent characteristic of a CME flux rope.
△ Less
Submitted 22 November, 2023;
originally announced November 2023.
-
CME Propagation Through the Heliosphere: Status and Future of Observations and Model Development
Authors:
M. Temmer,
C. Scolini,
I. G. Richardson,
S. G. Heinemann,
E. Paouris,
A. Vourlidas,
M. M. Bisi,
writing teams,
:,
N. Al-Haddad,
T. Amerstorfer,
L. Barnard,
D. Buresova,
S. J. Hofmeister,
K. Iwai,
B. V. Jackson,
R. Jarolim,
L. K. Jian,
J. A. Linker,
N. Lugaz,
P. K. Manoharan,
M. L. Mays,
W. Mishra,
M. J. Owens,
E. Palmerio
, et al. (9 additional authors not shown)
Abstract:
The ISWAT clusters H1+H2 have a focus on interplanetary space and its characteristics, especially on the large-scale co-rotating and transient structures impacting Earth. SIRs, generated by the interaction between high-speed solar wind originating in large-scale open coronal magnetic fields and slower solar wind from closed magnetic fields, are regions of compressed plasma and magnetic field follo…
▽ More
The ISWAT clusters H1+H2 have a focus on interplanetary space and its characteristics, especially on the large-scale co-rotating and transient structures impacting Earth. SIRs, generated by the interaction between high-speed solar wind originating in large-scale open coronal magnetic fields and slower solar wind from closed magnetic fields, are regions of compressed plasma and magnetic field followed by high-speed streams that recur at the ca. 27 day solar rotation period. Short-term reconfigurations of the lower coronal magnetic field generate flare emissions and provide the energy to accelerate enormous amounts of magnetised plasma and particles in the form of CMEs into interplanetary space. The dynamic interplay between these phenomena changes the configuration of interplanetary space on various temporal and spatial scales which in turn influences the propagation of individual structures. While considerable efforts have been made to model the solar wind, we outline the limitations arising from the rather large uncertainties in parameters inferred from observations that make reliable predictions of the structures impacting Earth difficult. Moreover, the increased complexity of interplanetary space as solar activity rises in cycle 25 is likely to pose a challenge to these models. Combining observational and modeling expertise will extend our knowledge of the relationship between these different phenomena and the underlying physical processes, leading to improved models and scientific understanding and more-reliable space-weather forecasting. The current paper summarizes the efforts and progress achieved in recent years, identifies open questions, and gives an outlook for the next 5-10 years. It acts as basis for updating the existing COSPAR roadmap by Schrijver+ (2015), as well as providing a useful and practical guide for peer-users and the next generation of space weather scientists.
△ Less
Submitted 9 August, 2023;
originally announced August 2023.
-
Impact of the solar activity on the propagation of ICMEs: Simulations of hydro, magnetic and median ICMEs at minimum and maximum of activity
Authors:
Barbara Perri,
Brigitte Schmieder,
Pascal Démoulin,
Stefaan Poedts,
Florian Regnault
Abstract:
The propagation of Interplanetary Coronal Mass Ejections (ICMEs) in the heliosphere is influenced by many physical phenomena, related to the internal structure of the ICME and its interaction with the ambient solar wind and magnetic field. As the solar magnetic field is modulated by the 11-year dynamo cycle, our goal is to perform a theoretical exploratory study to assess the difference of propaga…
▽ More
The propagation of Interplanetary Coronal Mass Ejections (ICMEs) in the heliosphere is influenced by many physical phenomena, related to the internal structure of the ICME and its interaction with the ambient solar wind and magnetic field. As the solar magnetic field is modulated by the 11-year dynamo cycle, our goal is to perform a theoretical exploratory study to assess the difference of propagation of an ICME in typical minimum and maximum activity backgrounds. We define a median representative CME at 0.1~au, using both observations and numerical simulations, and describe it using a spheromak model. We use the heliospheric propagator European Heliospheric FORecasting Information Asset (EUHFORIA) to inject the same ICME in two different background wind environments. We then study how the environment and the internal CME structure impact the propagation of the ICME towards Earth, by comparison with an unmagnetized CME. At minimum of activity, the structure of the heliosphere around the ecliptic causes the ICME to slow down, creating a delay with the polar parts of the ejecta. This delay is more important if the ICME is faster. At maximum of activity, a southern coronal hole causes a northward deflection. For these cases, we always find that the ICME at maximum of activity arrives first, while the ICME at minimum of activity is actually more geo-effective. The helicity sign of the ICME is also a crucial parameter but at minimum of activity only, since it affects the magnetic profile and the arrival time of up to 8 hours.
△ Less
Submitted 27 June, 2023;
originally announced June 2023.
-
COCONUT, a novel fast-converging MHD model for solar corona simulations: II. Assessing the impact of the input magnetic map on space-weather forecasting at minimum of activity
Authors:
Barbara Perri,
Blazej Kuzma,
Michaela Brchnelova,
Tinatin Baratashvili,
Fan Zhang,
Peter Leitner,
Andrea Lani,
Stefaan Poedts
Abstract:
This paper is dedicated to the new implicit unstructured coronal code COCONUT, which aims at providing fast and accurate inputs for space weather forecast as an alternative to empirical models. We use all 20 available magnetic maps of the solar photosphere covering the date of the 2nd of July 2019 which corresponds to a solar eclipse on Earth. We use the same standard pre-processing on all maps, t…
▽ More
This paper is dedicated to the new implicit unstructured coronal code COCONUT, which aims at providing fast and accurate inputs for space weather forecast as an alternative to empirical models. We use all 20 available magnetic maps of the solar photosphere covering the date of the 2nd of July 2019 which corresponds to a solar eclipse on Earth. We use the same standard pre-processing on all maps, then perform coronal MHD simulations with the same numerical and physical parameters. In the end, we quantify the performance for each map using three indicators from remote-sensing observations: white-light total solar eclipse images for the streamers' edges, EUV synoptic maps for coronal holes and white-light coronagraph images for the heliospheric current sheet. We discuss the performance for space weather forecasts and we show that the choice of the input magnetic map has a strong impact. We find performances between 24% to 85% for the streamers' edges, 24% to 88% for the coronal hole boundaries and a mean deviation between 4 to 12 degrees for the heliospheric current sheet position. We find that the HMI runs are globally performing better on all indicators, with the GONG-ADAPT being the second-best choice. HMI runs perform better for the streamers' edges, GONG-ADAPT for polar coronal holes, HMI synchronic for equatorial coronal holes and for the streamer belt. We especially showcase the importance of the filling of the poles. This demonstrates that the solar poles have to be taken into account even for ecliptic plane previsions.
△ Less
Submitted 12 October, 2022;
originally announced October 2022.
-
To E or not to E: Numerical Nuances of Global Coronal Models
Authors:
Michaela Brchnelova,
Błażej Kuźma,
Barbara Perri,
Andrea Lani,
Stefaan Poedts
Abstract:
In the recent years, global coronal models have experienced an ongoing increase in popularity as tools for forecasting solar weather. Within the domain of up to 21.5Rsun, magnetohydrodynamics (MHD) is used to resolve the coronal structure using magnetograms as inputs at the solar surface. Ideally, these computations would be repeated with every update of the solar magnetogram so that they could be…
▽ More
In the recent years, global coronal models have experienced an ongoing increase in popularity as tools for forecasting solar weather. Within the domain of up to 21.5Rsun, magnetohydrodynamics (MHD) is used to resolve the coronal structure using magnetograms as inputs at the solar surface. Ideally, these computations would be repeated with every update of the solar magnetogram so that they could be used in the ESA Modelling and Data Analysis Working Group (MADAWG) magnetic connectivity tool (http://connect-tool.irap.omp.eu/). Thus, it is crucial that these results are both accurate and efficient. While much work has been published showing the results of these models in comparison with observations, not many of it discusses the intricate numerical adjustments required to achieve these results. These range from details of boundary condition formulations to adjustments as large as enforcing parallelism between the magnetic field and velocity. By omitting the electric field in ideal-MHD, the description of the physics can be insufficient and may lead to excessive diffusion and incorrect profiles. We formulate inner boundary conditions which, along with other techniques, reduce artificial electric field generation. Moreover, we investigate how different outer boundary condition formulations and grid design affect the results and convergence, with special focus on the density and the radial component of the B-field. The significant improvement in accuracy of real magnetic map-driven simulations is illustrated for an example of the 2008 eclipse.
△ Less
Submitted 9 September, 2022;
originally announced September 2022.
-
Effects of Mesh Topology on MHD Solution Features in Coronal Simulations
Authors:
Michaela Brchnelova,
Fan Zhang,
Peter Leitner,
Barbara Perri,
Andrea Lani,
Stefaan Poedts
Abstract:
Magnetohydrodynamic (MHD) simulations of the solar corona have become more popular with the increased availability of computational power. Modern computational plasma codes, relying upon Computational Fluid Dynamics (CFD) methods, allow for resolving the coronal features using solar surface magnetograms as inputs. These computations are carried out in a full 3D domain and thus selection of the rig…
▽ More
Magnetohydrodynamic (MHD) simulations of the solar corona have become more popular with the increased availability of computational power. Modern computational plasma codes, relying upon Computational Fluid Dynamics (CFD) methods, allow for resolving the coronal features using solar surface magnetograms as inputs. These computations are carried out in a full 3D domain and thus selection of the right mesh configuration is essential to save computational resources and enable/speed up convergence. In addition, it has been observed that for MHD simulations close to the hydrostatic equilibrium, spurious numerical artefacts might appear in the solution following the mesh structure, which makes the selection of the grid also a concern for accuracy. The purpose of this paper is to discuss and trade off two main mesh topologies when applied to global solar corona simulations using the unstructured ideal MHD solver from the COOLFluiD platform. The first topology is based on the geodesic polyhedron and the second on UV mapping. Focus will be placed on aspects such as mesh adaptability, resolution distribution, resulting spurious numerical fluxes and convergence performance. For this purpose, firstly a rotating dipole case is investigated, followed by two simulations using real magnetograms from the solar minima (1995) and solar maxima (1999). It is concluded that the most appropriate mesh topology for the simulation depends on several factors, such as the accuracy requirements, the presence of features near the polar regions and/or strong features in the flow field in general. If convergence is of concern and the simulation contains strong dynamics, then grids which are based on the geodesic polyhedron are recommended compared to more conventionally used UV-mapped meshes.
△ Less
Submitted 28 February, 2022;
originally announced February 2022.
-
Unifying the Validation of Ambient Solar Wind Models
Authors:
Martin A. Reiss,
Karin Muglach,
Richard Mullinix,
Maria M. Kuznetsova,
Chiu Wiegand,
Manuela Temmer,
Charles N. Arge,
Sergio Dasso,
Shing F. Fung,
Jose Juan Gonzalez Aviles,
Siegfried Gonzi,
Lan Jian,
Peter MacNeice,
Christian Möstl,
Mathew Owens,
Barbara Perri,
Rui F. Pinto,
Lutz Rastätter,
Pete Riley,
Evangelia Samara,
ISWAT H1-01 Team Members
Abstract:
Progress in space weather research and awareness needs community-wide strategies and procedures to evaluate our modeling assets. Here we present the activities of the Ambient Solar Wind Validation Team embedded in the COSPAR ISWAT initiative. We aim to bridge the gap between model developers and end-users to provide the community with an assessment of the state-of-the-art in solar wind forecasting…
▽ More
Progress in space weather research and awareness needs community-wide strategies and procedures to evaluate our modeling assets. Here we present the activities of the Ambient Solar Wind Validation Team embedded in the COSPAR ISWAT initiative. We aim to bridge the gap between model developers and end-users to provide the community with an assessment of the state-of-the-art in solar wind forecasting. To this end, we develop an open online platform for validating solar wind models by comparing their solutions with in situ spacecraft measurements. The online platform will allow the space weather community to test the quality of state-of-the-art solar wind models with unified metrics providing an unbiased assessment of progress over time. In this study, we propose a metadata architecture and recommend community-wide forecasting goals and validation metrics. We conclude with a status update of the online platform and outline future perspectives.
△ Less
Submitted 16 May, 2022; v1 submitted 31 January, 2022;
originally announced January 2022.
-
Powering Stellar Magnetism: Energy Transfers in Cyclic Dynamos of Sun-like Stars
Authors:
Allan Sacha Brun,
Antoine Strugarek,
Quentin Noraz,
Barbara Perri,
Jacobo Varela,
Kyle Augsutson,
Paul Charbonneau,
Juri Toomre
Abstract:
We use the ASH code to model the convective dynamo of solar-type stars. Based on a series of 15 3-D MHD simulations spanning 4 bins in rotation and mass, we show what mechanisms are at work in these stellar dynamos with and without magnetic cycles and how global stellar parameters affect the outcome. We also derive scaling laws for the differential rotation and magnetic field based on these simula…
▽ More
We use the ASH code to model the convective dynamo of solar-type stars. Based on a series of 15 3-D MHD simulations spanning 4 bins in rotation and mass, we show what mechanisms are at work in these stellar dynamos with and without magnetic cycles and how global stellar parameters affect the outcome. We also derive scaling laws for the differential rotation and magnetic field based on these simulations. We find a weaker trend between differential rotation and stellar rotation rate, ($ΔΩ\sim (|Ω|/Ω_{\odot})^{0.46}$) in the MHD solutions than in their HD counterpart $(|Ω|/Ω_{\odot})^{0.66})$, yielding a better agreement with the observational trends based on power laws. We find that for a fluid Rossby number between $0.15 \lesssim Ro_f \lesssim 0.65$ the solutions possess long magnetic cycle, if $Ro_f \lesssim 0.42$ a short cycle and if $Ro_f \gtrsim 1$ (anti-solar-like differential rotation) a statistically steady state. We show that short-cycle dynamos follow the classical Parker-Yoshimura rule whereas the long-cycle period ones do not. We further demonstrate that the Rossby number dependency of the large-scale surface magnetic field in the simulation ($B_{L,surf} \sim Ro_{f}^{-1.26}$) agrees better with observations ($B_{V} \sim Ro_{s}^{-1.4 \pm 0.1}$) and differs from dynamo scaling based on the global magnetic energy ($B_{bulk} \sim Ro_{f}^{-0.5}$). We also show that up to few percents of the stellar luminosity can be channelled into the star's magnetism, hence providing a large energy reservoir for possible surface eruptive events.
△ Less
Submitted 31 January, 2022;
originally announced January 2022.
-
The two-step Forbush decrease: a tale of two substructures modulating galactic cosmic rays within coronal mass ejections
Authors:
Miho Janvier,
Pascal Démoulin,
Jingnan Guo,
Sergio Dasso,
Florian Regnault,
Sofia Topsi-Moutesidou,
Christian Gutierrez,
Barbara Perri
Abstract:
Interplanetary Coronal Mass Ejections (ICMEs) are known to modify the structure of the solar wind as well as interact with the space environment of planetary systems. Their large magnetic structures have been shown to interact with galactic cosmic rays, leading to the Forbush decrease (FD) phenomenon. We revisit in the present article the 17 years of Advanced Composition Explorer spacecraft ICME d…
▽ More
Interplanetary Coronal Mass Ejections (ICMEs) are known to modify the structure of the solar wind as well as interact with the space environment of planetary systems. Their large magnetic structures have been shown to interact with galactic cosmic rays, leading to the Forbush decrease (FD) phenomenon. We revisit in the present article the 17 years of Advanced Composition Explorer spacecraft ICME detection along with two neutron monitors (McMurdo and Oulu) with a superposed epoch analysis to further analyze the role of the magnetic ejecta in driving FDs. We investigate in the following the role of the sheath and the magnetic ejecta in driving FDs, and we further show that for ICMEs without a sheath, a magnetic ejecta only is able to drive significant FDs of comparable intensities. Furthermore, a comparison of samples with and without a sheath with similar speed profiles enable us to show that the magnetic field intensity, rather than its fluctuations, is the main driver for the FD. Finally, the recovery phase of the FD for isolated magnetic ejecta shows an anisotropy in the level of the GCRs. We relate this finding at 1 au to the gradient of the GCR flux found at different heliospheric distances from several interplanetary missions.
△ Less
Submitted 29 September, 2021;
originally announced September 2021.
-
Modeling Solar Wind Variations over an 11-yr Cycle with Alfvén Wave Dissipation: a Parameter Study
Authors:
Soumitra Hazra,
Victor Réville,
Barbara Perri,
Antoine Strugarek,
Allan Sacha Brun,
Eric Buchlin
Abstract:
We study the behaviour and properties of the solar wind using a 2.5D Alfvén wave driven wind model. We first systematically compare the results of an Alfvén wave (AW) driven wind model with a polytropic approach. Polytropic magnetohydrodynamic wind models are thermally driven, while Alfvén waves act as additional acceleration and heating mechanisms in the Alfvén wave driven model. We confirm that…
▽ More
We study the behaviour and properties of the solar wind using a 2.5D Alfvén wave driven wind model. We first systematically compare the results of an Alfvén wave (AW) driven wind model with a polytropic approach. Polytropic magnetohydrodynamic wind models are thermally driven, while Alfvén waves act as additional acceleration and heating mechanisms in the Alfvén wave driven model. We confirm that an AW-driven model is required to reproduce the observed bimodality of slow and fast solar winds. We are also able to reproduce the observed anti-correlation between the terminal wind velocity and the coronal source temperature with the AW-driven wind model. We also show that the wind properties along an eleven-year cycle differ significantly from one model to the other. The AW-driven model again shows the best agreement with observational data. Indeed, solar surface magnetic field topology plays an important role in the Alfvén wave driven wind model, as it enters directly into the input energy sources via the Poynting flux. On the other hand, the polytropic wind model is driven by an assumed pressure gradient; thus it is relatively less sensitive to the surface magnetic field topology. Finally, we note that the net torque spinning down the Sun exhibits the same trends in the two models, showing that the polytropic approach still captures correctly the essence of stellar winds.
△ Less
Submitted 27 January, 2021;
originally announced January 2021.
-
Impact of solar magnetic field amplitude and geometry on cosmic rays diffusion coefficients in the inner heliosphere
Authors:
Barbara Perri,
Allan Sacha Brun,
Antoine Strugarek,
Victor Réville
Abstract:
Cosmic rays (CRs) are tracers of solar events when they are associated with solar flares, but also galactic events when they come from outside our solar system. SEPs are correlated with the 11-year solar cycle while GCRs are anti-correlated due to their interaction with the heliospheric magnetic field and the solar wind. Our aim is to quantify separately the impact of the amplitude and the geometr…
▽ More
Cosmic rays (CRs) are tracers of solar events when they are associated with solar flares, but also galactic events when they come from outside our solar system. SEPs are correlated with the 11-year solar cycle while GCRs are anti-correlated due to their interaction with the heliospheric magnetic field and the solar wind. Our aim is to quantify separately the impact of the amplitude and the geometry of the magnetic field on the propagation of CRs of various energies in the inner heliosphere. We focus especially on the diffusion caused by the magnetic field along and across the field lines. To do so, we use the results of 3D MHD wind simulations running from the lower corona up to 1 AU. The wind is modeled using a polytropic approximation, and fits and power laws are used to account for the turbulence. Using these results, we compute the parallel and perpendicular diffusion coefficients of the Parker CR transport equation, yielding 3D maps of the diffusion of CRs in the inner heliosphere. By varying the amplitude of the magnetic field, we change the amplitude of the diffusion by the same factor, and the radial gradients by changing the spread of the current sheet. By varying the geometry of the magnetic field, we change the latitudinal gradients of diffusion by changing the position of the current sheets. By varying the energy, we show that the distribution of values for SEPs is more peaked than GCRs. For realistic solar configurations, we show that diffusion is highly non-axisymmetric due to the configuration of the current sheets, and that the distribution varies a lot with the distance to the Sun with a drift of the peak value. This study shows that numerical simulations and theory can help quantify better the influence of the various magnetic field parameters on the propagation of CRs. This study is a first step towards generating synthetic CR rates from numerical simulations.
△ Less
Submitted 5 October, 2020;
originally announced October 2020.