Issue 
J. Space Weather Space Clim.
Volume 14, 2024
Topical Issue  CMEs, ICMEs, SEPs: Observational, Modelling, and Forecasting Advances



Article Number  13  
Number of page(s)  17  
DOI  https://doi.org/10.1051/swsc/2024011  
Published online  08 May 2024 
Research Article
The spheroid CME model in EUHFORIA
^{1}
SolarTerrestrial Centre of Excellence – SIDC, Royal Observatory of Belgium, Av. Circulaire 3, 1180 Brussels, Belgium
^{2}
Space Science Center, University of New Hampshire, 8 College Road, Durham, NH 03824, USA
^{3}
Predictive Science Inc., 9990 Mesa Rim Rd, Suite 170, San Diego, CA 92121, USA
^{*} Corresponding author: camilla.scolini@oma.be
Received:
16
September
2023
Accepted:
4
April
2024
Predictions of coronal mass ejection (CME) propagation and impact in the heliosphere, in either research or operational settings, are usually performed by employing magnetohydrodynamic (MHD) models. Within such simulations, the CME ejecta is often described as a hydrodynamic pulse that lacks an internal magnetic field and is characterized by a spherical geometry – leading to the socalled cone CME model. Whitelight observations of CMEs in the corona, however, reveal that the morphology of these structures resembles more closely that of a croissant, i.e., exhibiting an elongated crosssection of their front. It follows that, in space weather forecasts, the assumption of a spherical geometry may result in erroneous predictions of CME impacts in the heliosphere in terms of hit/miss and arrival time/speed, especially in the case of flank encounters. A spheroid CME model is expected to provide a more accurate description of the elongated morphology that is often observed in CMEs. In this paper, we describe the implementation and initial validation of the spheroid CME model within the MHD EUropean Heliospheric FORecasting Information Asset (EUHFORIA) code. We perform EUHFORIA simulations of an idealized CME as well as a “real” event to compare the spheroidal model with the traditional cone one. We show how the initial ejecta geometry can lead to substantially different estimates in terms of CME impact, arrival time/speed, and geoeffectiveness, especially with increasing distance to the CME nose.
Key words: Coronal mass ejections / Space weather predictions / Magnetohydrodynamics
© C. Scolini & E. Palmerio, Published by EDP Sciences 2024
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Coronal mass ejections (CMEs) are among the most severe drivers of geomagnetic effects on Earth and, as such, are routinely monitored, analyzed, and modeled by space weather forecasting agencies (e.g., Pizzo et al., 2011). Once a CME erupts, its geometric and kinematic properties are usually estimated from whitelight coronagraph imagery, and the resulting set of parameters is employed as input for MagnetoHydroDynamic (MHD) models of CME propagation in the inner heliosphere. A variety of models have been developed to describe CME structures in the context of heliospheric MHD models, including magnetized fluxrope models that account for various geometrical and magnetic field descriptions of CMEs (see e.g. Kataoka et al., 2009; Zhou et al., 2014; Zhang et al., 2019; Singh et al., 2020; Shen et al., 2021a,2021b; Maharana et al., 2022). However, for the purpose of operational space weather forecasting, CMEs in MHD models are typically treated without an internal magnetic field description, i.e., the injected structure is simplified as a hydrodynamic cloud or pulse (e.g., Odstrcil, 2003; Pomoell & Poedts, 2018). This approach – despite yielding not fully realistic CME interactions with the ambient solar wind and/or other solar transients – not only allows for more rapid predictions in terms of computational time but also avoids the complex description of the magnetic configuration of the embedded flux rope, which is more difficult to characterize in realtime (e.g., Kilpua et al., 2019). The threedimensional (3D) deprojected geometric and kinematic properties of CMEs in the corona are usually derived via triangulation techniques, such as the Space Weather Prediction Center CME Analysis Tool (SWPCCAT; Millward et al., 2013, commonly employed in operational settings) and the Graduated Cylindrical Shell (GCS; Thernisien, 2011, widely used in the research community) model.
When observed in the solar corona, CMEs often exhibit an elongated crosssection of their front, leading to a socalled “croissant” morphology that is often associated with an internal flux rope structure (e.g., Cremades & Bothmer, 2004; Vourlidas et al., 2013). To take this aspect into account, many forward models employed in CME research – such as the GCS technique mentioned above, or the coronagraph fitting tool implemented within the Flux Rope in 3D (FRi3D; Isavnin, 2016) architecture – are described by a number of free parameters that allow for a toroidal axis of the structure and two wellseparated legs that connect back to the Sun. However, these models are only meant to reproduce the morphology of CMEs in the corona, and thus consist of “hollow” parameterized shells, with no information on the internal magnetic configuration – especially since coronal fields cannot be measured routinely at present. On the other hand, these reconstruction techniques permit relatively quick estimates of the geometric and kinematic properties of CMEs using one or more viewpoints, and their outputs can be directly employed in CME propagation models, either on their own (leading to the injection of a hydrodynamic pulse as described above) or combined with indirect information on the magnetic field of the corresponding flux rope (e.g., from solar disk observations; Palmerio et al., 2017).
The CME geometry assumed by most CME propagation models themselves (especially those used in forecasting), however, is much simpler, and is typically described as a spherical shape. In other words, these models are significantly more compatible with the geometry of the socalled “icecream cone” (Fisher & Munro, 1984) description of CMEs than with that of the croissant – which, on the other hand, is generally regarded as more realistic. An example of a CME reconstructed using whitelight data employing both the cone and croissant models is shown in Figure 1. When adapting fitting results based on coronagraph imagery to CME propagation models that assume a spherical morphology, there are two possible approaches. The first is to reconstruct the CME by applying directly the cone model (e.g., Palmerio et al., 2019; Scolini et al., 2019), while the second is to first use the croissant model and then derive a circular front by considering, e.g., the maximum elongation or the average radius of the resulting structure (e.g., Scolini et al., 2020; Asvestari et al., 2021). Either way, not accounting for the observed elongation in the CME frontal crosssection can result in wrong predictions of a CME’s impact at a certain location in terms of its hit/miss outcome, as well as of its arrival time and/or speed, especially in the case of glancing encounters. A prominent example is that of the 7 January 2014 event (featured in Fig. 1), whose interplanetary propagation was studied in detail by by Mays et al. (2015) and Möstl et al. (2015). This CME, despite originating from the vicinity of the central meridian of the Earthfacing Sun and being very fast and energetic, experienced a strong westward deflection in the solar corona and resulted only in a flank encounter at Earth. Mays et al. (2015) showed that, in this case, hindcasts that assumed a spherical CME geometry performed significantly worse than those that employed a tilted ellipsoidal CME morphology.
Figure 1 Example of a CME (from 7 January 2014) fitted in the solar corona using the cone and croissant models. The top row (a–c) shows whitelight data from three viewpoints (Earth/SOHO, STEREOA, and STEREOB), the middle row (d–f) shows the same set of data with the cone wireframe (in green) overlaid, and the bottom row (g–i) shows the same set of data with the croissant wireframe (in magenta) overlaid. Both sets of reconstructions were performed using the GCS model, with the difference that for the cone case the socalled halfwidth parameter (quantifying the halfangular separation between the CME legs) was set to zero. 
In this work, we focus on the MHD EUropean Heliospheric FORecasting Information Asset (EUHFORIA; Pomoell & Poedts, 2018) model. Specifically, to improve current EUHFORIA space weather research and forecasting capabilities, we have increased the flexibility of the (default) cone model by introducing the modeling of CME geometry as a spheroid. In this paper, we report the details of such a new capability and present a first validation, which is achieved through comparison with the standard cone model and with insitu signatures (typical CME or wellstudied event). The paper is structured as follows. Section 2 presents a technical description of the model implementation. Section 3 presents a first validation of the model capabilities for a hypothetical (idealized) CME event, while Section 4 presents the improvement in the case of a real CME event observed on 7 January 2014 and previously studied by Mays et al. (2015). Section 5 summarizes our conclusions and future perspectives.
2 Spheroid CME model implementation
The spheroid CME model is a spheroidshaped hydrodynamical CME model, similar to the “traditional” cone model described by Pomoell & Poedts (2018). The CME is assumed to propagate radially outwards from the source location with a uniform speed, density, and temperature throughout. Similarly to the traditional cone model, the spheroidal model requires as input the following geometric and kinematic parameters at 0.1 au (i.e., the heliospheric inner boundary of EUHFORIA): the latitude and longitude of the CME center (θ_{CME} and ϕ_{CME}, in HEEQ coordinates), the speed of the CME (v_{CME}, in km s^{−1}), and the time at which the CME front first intersects the model inner boundary (t_{CME}). However, while the cone model traditionally assumes a spherical CME shape, the newly implemented spheroidal model allows for a more accurate description of the elongated morphology that is often observed in CMEs. Thus, while the cone model only requires the halfangular width as an input parameter (ω_{CME}/2, in degrees), the spheroidal model introduced here requires three parameters to characterize a given CME’s geometry: the semimajor axis (r_{maj}, in solar radii, R_{s}), the semiminor axis(r_{min}, in R_{s}), and the tilt angle (γ_{CME}, in degrees and computed anticlockwise from the solar west direction as seen from the equatorial plane). In addition to the aforementioned set of geometrical and kinematic parameters, the spheroidal CME model also assumes a homogeneous density (ρ_{CME}) and temperature (T_{CME}) as does the traditional cone model.
We note that in the traditional cone model first introduced in its EUHFORIA implementation by Pomoell & Poedts (2018), the conversion from the CME angular width to the CME radius in the code is implemented based on the tangent relation (Eq. 1 in Scolini et al., 2018a). For the purpose of directly comparing with the traditional cone model, here we decide to keep this relation to operate the translation between the angular width and the radius of the CME at 0.1 au in the spheroidal model as well, although a more appropriate alternative is provided by Equation 2 in Scolini et al. (2018a).
In detail, for each (θ, ϕ) cell part of the inner boundary of the model domain, the polar angle with respect to the CME center is computed as
$${\theta}^{\text{'}}\left(\theta ,\varphi \right)=\mathrm{arctan}\left(\frac{\theta {\theta}_{\mathrm{CME}}}{\varphi {\varphi}_{\mathrm{CME}}}\right).$$(1)
Then, a static mask (in units of R_{s}) is introduced to identify which cells on the inner boundary are going to be part of the CME frontal crosssection during the insertion period, as follows:
$$\mathrm{mas}{\mathrm{k}}_{{R}_{\mathrm{s}}}\left(\theta ,\varphi \right)=\frac{{r}_{\mathrm{min}}{r}_{\mathrm{maj}}}{\sqrt{\left({r}_{\mathrm{min}}^{2}{r}_{\mathrm{maj}}^{2}\right){\mathrm{cos}}^{2}\left({\theta}^{\text{'}}\left(\theta ,\varphi \right){\gamma}_{\mathrm{CME}}\right)+{r}_{\mathrm{maj}}^{2}}}.$$(2)
The normalization to radians is performed as
$$\mathrm{mas}{\mathrm{k}}_{\mathrm{rad}}\left(\theta ,\varphi \right)=\mathrm{arctan}\left(\frac{\mathrm{mas}{\mathrm{k}}_{{R}_{\mathrm{s}}}\left(\theta ,\varphi \right)}{21.5{R}_{\mathrm{s}}}\right).$$(3)
The time dependence of the mask is obtained by multiplying the static mask by a timedependent “normalized opening angle” (running between 0 and 1), similar to what is done for the traditional cone model in EUHFORIA (Pomoell & Poedts, 2018; Scolini et al., 2018a):
$$\mathrm{mas}{\mathrm{k}}_{\mathrm{rad}}\left(\theta ,\varphi ,t\right)=\mathrm{mas}{\mathrm{k}}_{\mathrm{rad}}\left(\theta ,\varphi \right)\cdot \zeta \left(t\right)=\mathrm{mas}{\mathrm{k}}_{\mathrm{rad}}\left(\theta ,\varphi \right)\cdot \mathrm{sin}\left(\frac{\pi}{2}\frac{t}{{t}_{1/2}}\right),$$(4)
where t_{1/2} = r_{min}/v_{CME} is the time at which the CME geometrical center crosses the inner boundary. We impose that at any given time t during the CME insertion period, all (θ, ϕ) points on the inner boundary that satisfies the condition
$$(\theta {\theta}_{\mathrm{CME}}{)}^{2}+(\varphi {\varphi}_{\mathrm{CME}}{)}^{2}<\mathrm{mas}{\mathrm{k}}_{\mathrm{rad}}^{2}(\theta ,\varphi ,t)$$(5)
are inside the CME, while all that do not satisfy this relation are outside of the CME. At each point within the CME, the speed, density, and temperature are initialized to the v_{CME}, ρ_{CME}, and T_{CME}specified by the user. For each point outside of the CME, the values are retained to those of the ambient solar wind. Because this is an unmagnetized CME model, for all points the magnetic field conditions are retained to those of the ambient solar wind.
As a result, in the spheroidal model CMEs have a spherical crosssection only in the r–θ plane, while they have elliptical crosssections when sliced in θ–ϕ surfaces of constant r. Figure 2 provides a 3D visualization of the geometry of a spherical and spheroidal CME shape, as well as a projection of their crosssections onto the various planes (for a Cartesian coordinate system x′, y′, and z′ centered at the CME center).
Figure 2 3D visualization of the (a) spherical and (b) spheroidal CME shapes. x′, y′, and z′ constitute a local coordinate system centered at the center of the CME. Dimensions are normalized to 1 for convenience. The spherical CME has circular crosssections in the directions perpendicular to x′, y′, and z′, as shown by the 2D shadow projections. The spheroidal CME has a circular crosssections only in the direction perpendicular to the major axis (parallel to y′), while it has elliptical crosssections in the directions perpendicular to the minor axes (parallel to x′ and y′), as shown by the 2D shadow projections. 
3 Spheroid CME model validation: fictitious CME event
To exemplify the similarities between the cone and spheroidal CME models, as well as to highlight the advantages associated with the use of a spheroidal model over the cone model under certain circumstances, we simulate hypothetical CME events propagating through a synthetic solar wind background. This solar wind background is the same as previously employed by Scolini et al. (2021, 2023), and it includes a heliospheric current and plasma sheet with sinusoidal behavior in longitude which reaches up to ±15° in latitude. The solar wind has a baseline value of 450 km s^{−1} (intermediate between slow and fast solar wind) everywhere except for regions around the plasma sheet, where the speed reaches values as low as 300 km s^{−1}. Full details are provided by Scolini et al. (2021, 2023), while a 2D full view of the solar wind mass density and magnetic field boundary conditions is provided in Figure 1 by Scolini et al. (2021). Characterizing the behavior of this newlyimplemented spheroidal CME model in EUHFORIA in such an idealized solar wind condition represents a first but necessary step to later be able to conduct a meaningful interpretation of the results obtained for CME structures propagating through different (incl. more realistic, see Sect. 4) solar wind conditions.
We run the same fictitious CME event with both the cone and spheroidal models, using input parameters as similar as possible. The complete set of input CME parameters at 21.5 R_{s} used to initialize the EUHFORIA runs is listed in Table 1. For the spheroidal model, we choose the major axis to be aligned with the equatorial plane (corresponding to a tilt angle of γ_{CME} = 0°), so that the advantage of using a spheroidal CME geometry is maximal in the region around the ecliptic plane (i.e., where most planets, as well as most past and currently operational spacecraft probing the solar wind conditions are located).
CME input parameters used for the spherical (cone) and spheroidal models to initialize the fictitious CME cases.
In addition to the input parameters required by EUHFORIA, the bottom portion of Table 1 lists some derived parameters, namely the CME total volumes (V_{CME}), masses (m_{CME}), and kinetic energies (E_{kin,CME}). As a result of the choices of input parameters, these differ by less than 2% in the two models, and they are representative of the typical population of fast CMEs observed in Solar Cycles 23 and 24 (Gopalswamy, 2006; Jang et al., 2016; Lamy et al., 2019). Despite the idealized numerical setup employed in this fictitious CME experiment, we therefore expect the results to apply to a potentially large set of real CME events. In the spheroidal model, we set r_{maj} to be twice as large as r_{min}, so that the CME aspect ratio is 2:1 in the longitudinal versus latitudinal directions.
3.1 Comparison with the spherical (cone) model at 0.1 au
Figure 3 shows the radial speed and particle number density at the inner boundary (r = 0.1 au) during the CME insertion period for the cone and spheroidal simulations. The area where the CME is being inserted (v_{r} = 800 km s^{−1}, n = 600 cm^{−3}) is visible as the red (orange) patch on the left (right). The CME crosssection appears visibly different between the two runs, especially in the longitudinal direction where the CME has a radius of 10.02 R_{s} in the cone model run, and of 16.20 R_{s} in the spheroidal model run. On the other hand, in the latitudinal direction the two CME present a more similar extension, as dictated by their radius of 10.02 R_{s} in the cone model run and of 7.82 R_{s} in the spheroidal model run.
Figure 3 3D views of the radial velocity (v_{r}, in units of km s^{−1}) and number density (n, in units of cm^{−3}) at the heliospheric inner boundary during the CME insertion, for the two simulations performed. (a) v_{r} for the cone (spherical) CME model, (b) n for the cone (spherical) CME model, (c) v_{r} for the spheroidal CME model, (d) n for the spheroidal CME model. The area where the CME is being inserted (v_{r} = 800 km s^{−3}, n = 600 cm^{−3}) is visible as the red (orange) patch on the left (right). 
As a result, even if the injected CMEs are both centered on the fictitious Sun–Earth line, the spheroid one is expected to be better suited to reproduce insitu arrivals across a wider range of helio longitudes – particularly amenable for flank encounters and multispacecraft validations – while maintaining a similar latitudinal extent.
3.2 Comparison with the spherical (cone) model in the heliosphere
Next, we explore how the differences at the model inner boundary of 0.1 au translate in the heliospheric domain. Figure 4 shows the radial speed in the ecliptic and meridional planes for the sphere (cone) and spheroid CMEs. The results in the two runs appear rather similar in the meridional plane, where the two CMEs are initiated with a difference of only 5° (∼2 R_{s}) in their angular extent. On the other hand, the spheroid CME is visibly wider in the equatorial plane, which coincides with the direction of its major axis. Although the CME nose is located at a similar heliocentric distance for both simulations at the time of the snapshot shown in Figure 4 (suggesting a similar arrival time at 1 au, see Section 3.3), it is clear that the spheroid CME appears approximately 20° wider than the cone one – note that the green lines in Figures 4a and 4c mark increments of 10° in longitude from the fictitious Sun–Earth line, located at (θ, ϕ) = (0°, 0°).
Figure 4 Modeling results for the fictitious CME with the sphere and spheroidal CME models. The panels show the radial speed v_{r} in the equatorial (left) and meridional (right) plane on 20230101 05:57 UT when the CME nose was close to 1 au in both simulations. Representative scaled number density nr^{2} values are indicated by the gray contour lines. (a, b) spherical (cone) CME model, (c, d) spheroidal CME model. The green lines marks directions to representative virtual spacecraft located at θ = 0° and ϕ between −90° and 90°, with 10° increments. 
3.3 Comparison with the spherical (cone) model at 1 au
Figure 5 shows the spherical (cone) and spheroidal run time series for the radial speed and number density at 1 au at the CME nose (θ, ϕ) = (0°, 0°) and as a function of the impact parameter (from the center to the flanks, with virtual spacecraft placed at 10° separation). Locations are shown at longitudes of 0° (a), 10° (b), 20° (c), 30° (d), 40° (e), and 50° (f) on the western CME flank.
Figure 5 Radial speed and number density recorded by virtual spacecraft at 1 au for the cone (teal) and spheroid (violet) CME models. As a reference, the ambient solar wind values are plotted in gray. Locations are shown at ϕ values of (a) 0°, (b) 10°, (c) 20°, (d) 30°, (e) 40°, and (f) 50°. Vertical dashed lines mark the CME arrival time in the two runs. Dots mark the CME impact (peak) speed and density. 
We identify the CME arrival time at each location and for each run by comparing the CME time series with a time series obtained from running the sole solar wind without any CME inserted, and by identifying the first time the two time series diverge in terms of their speed and/or density parameters. The resulting CME arrival times at each location are marked by vertical lines in Figure 5. Comparable results are found at the spacecraft locations on the eastern CME flank (not shown).
Even from a visual inspection of the results, we can appreciate how the difference in the CME arrival time increases with the longitudinal offset from the CME longitudinal direction of propagation (initiated at ϕ = 0°).
In Figure 6, we further quantify this longitudinal effect for the CME arrival time, impact speed, and impact density.
Figure 6 Difference between the spheroid and sphere (cone) CME impact properties at 1 au, as a function of the spacecraft longitude. (a) Difference in the CME arrival time, (b) Difference in the CME impact (peak) speed, (c) Difference in the CME impact (peak) density. The orange lines are calculated by subtracting the sphere to the spheroid properties. The teal (violet) lines show the subtracted values between the sphere (spheroid) and the ambient solar wind. 
Panel (a) reports the time difference in the CME arrival time between the two simulations, calculated as Δt = t_{spheroid} − t_{sphere}. We observe that while at ϕ = ϕ_{CME} = 0° the CMEs in the two runs arrive within ∼1 h from each other, the larger the longitudinal offset, the more the CME in the cone model lags behind the spheroidal CME. This delay reaches a maximum of ∼5–6 h at ϕ = ±50°, and is consistent with the broader front of the spheroidal CME (due to its larger width) in the equatorial plane compared to the cone CME. Conversely, the cone CME is characterized by a more convex front due to its spherical, narrower geometry. Furthermore, spacecraft located at ϕ = ±60° only detect the passage of the CME with spheroidal geometry (corresponding to the missing points in Fig. 6a). Spacecraft located at larger longitudinal separations from the CME direction of propagation do not detect the CME in either of the two simulations (corresponding to the zero values in Fig. 6a). It is important to note that although the values for Δt reported for this specific case are smaller than the current prediction capability for the CME arrival time (estimated to be around ±10 h; Riley et al., 2018; Vourlidas et al., 2019), our results demonstrate how effects related to the CME geometry are contributing to the overall arrival time prediction uncertainty unless the CME geometry is not properly accounted for. Furthermore, Δt may be larger, possibly exceeding ±10 h, for other combinations of CME and solar wind parameters: in such a case, these differences would become comparable to (if not greater than) prediction uncertainties, thus directly affecting the prediction capabilities for given CME events. These results also demonstrate how the choice of CME geometry not only affects predictions of the CME arrival time, but also predictions of the CME hit/miss at a given spacecraft location.
Figures 6b and 6c show the results for the difference in the impact (i.e., peak) CME speed and density, calculated as Δv_{r} = v_{r,spheroid} − v_{r,sphere} and Δn = n_{spheroid} − n_{sphere}. Spacecraft within ϕ = ±10° correspond to the locations where the cone CME arrives slightly earlier (Δt > 0) and faster (Δv°<°0) than the spheroid CME. In this region, the cone CME is also denser (Δn°<°0). Conversely, spacecraft at 20° ≤ ϕ ≤ 50° show the spheroidal CME arriving earlier (Δt < 0), faster (Δv > 0), and denser (Δn > 0) than the cone one, consistently with its broader front. Overall, Δv ranges from −20 km s^{−1} at ϕ = 0°, to 32 (20) km s^{−1} at ϕ = −40° (+40°), while Δn goes from −2 cm^{−3} at ϕ = 0°, to 3 (2) cm^{−3} at ϕ = −40° (+40°). While these Δv and Δn values are somehow moderate when considered individually, they are expected to have a cumulative effect in the level of perturbation induced on geospace. We explore this aspect in Section 3.4 by estimating the level of CME geoeffectiveness hypothetically expected at each synthetic spacecraft location (each representing the location of a fictitious Earth) in terms of induced magnetospheric compression and geomagnetic storm strength.
3.4 Comparison with the spherical (cone) model: prediction of CME geoeffectiveness
To conclude our analysis of fictitious CME events, we want quantify the potential impact of different CME shapes on the space weather predictions at 1 au. To do so, we follow the approach by Scolini et al. (2018a) and use EUHFORIA simulation outputs at various virtual spacecraft at 1 au to compute: (1) the minimum magnetopause standoff distance (dso) in the subsolar direction, based on the Shue model (Shue et al., 1997); and (2) the maximum Kp index, proxy of the CMEdriven geomagnetic activity, by means of the coupling function proposed by Newell et al. (2008). The choice of these two metrics is based on the consideration that these are mainly susceptible to the CME dynamic pressure at impact, rather than to its internal magnetic structure (which is not realistically modeled in the simulations considered here). Furthemore, both official forecasts (e.g., from NOAA Space Weather Prediction Center, https://www.swpc.noaa.gov/products/3dayforecast) and researchoriented forecasts (e.g., in the NASA Community Coordinated Modeling Center DONKI catalog, https://ccmc.gsfc.nasa.gov/tools/DONKI/) are based on predictions from solar wind–Kp coupling functions, which makes the predicted Kp a standard metric for this sort of validation efforts. Despite the limitations common to such coupling functions (e.g., Luo et al., 2017; Lockwood, 2022), the use of the specific solar wind–Kp coupling function by Newell et al. (2008) also has the advantage to enable backcomparison with previous works on the validation of the EUHFORIA cone model (e.g., Scolini et al., 2018a,b). We calculate the predicted minimum dso and maximum Kp index based on EUHFORIA in situ time series for the full array of virtual spacecraft located at θ = 0° and evenly spanned from ϕ = −90° to ϕ = 90°.
Figure 7 compares the predictions obtained from EUHFORIA time series at different virtual spacecraft for the two CME models. Both the minimum magnetopause standoff distance and the maximum Kp index significantly depend on the impact location and CME model used. Figure 7a shows that the minimum magnetopause standoff distance is moderately affected by the different CME models tested. The predicted minimum magnetopause standoff distance is expected at ϕ = 0° in both models, where it reaches altitudes of 6.1 earth radii (R_{e}) and 6.3 R_{e} for the cone and spheroidal models, respectively. Figure 7a also shows that the choice of CME model influences whether and at which locations the magnetopause would be expected to cross the geosynchronous orbits (located at an altitude of 6.6 R_{e}), thereby significantly affecting the prediction of the expected CME impact on technological systems in space. In our simulations, the spacecraft locations where predictions are most uncertain with respect to whether the magnetopause will cross or not the geosynchronous orbits (meaning the cone and spheroidal model do not agree in their prediction) are those located at ϕ values between 20° and 30°. All locations at ϕ ≤ 20° are predicted to have the magnetopause pushed below geosynchronous orbits, while all locations at ϕ ≥ 30° are predicted to have the magnetopause remaining above geosynchronous orbits. The maximum difference in the predicted minimum magnetopause standoff distance between the cone and spheroidal model is of 0.7 R_{e} at ϕ = −40°, confirming predictions related to the impact of CME flanks are those most susceptible to the particular CME model used.
Figure 7 CME geoeffectiveness prediction as a function of the spacecraft longitude. (a) Estimate of the minimum magnetopause standoff distance calculated using the model by Shue et al. (1997). The black horizontal line marks the geosynchronous orbits at 6.6 R_{e}. (b) Estimate of the maximum Kp index calculated using the relation by Newell et al. (2008). 
Figure 7b shows the results for the predicted maximum Kp index as function of the impact location and CME model used. The maximum Kp values predicted are 5.3 and 5.0 (corresponding to a minor storm level according to the National Oceanic and Atmospheric Administration (NOAA) Geomagnetic Storm Scale classification^{1}) in the cone and spheroidal model, respectively. Such values are expected at ϕ = 0°. At larger ϕ values the predicted maximum Kp ranges between 4.9 and 1.9, corresponding to quiettime levels. The maximum difference in the predicted maximum Kp between the two models is of 0.74 and is achieved at ϕ = −30°, showing that, just as for the magnetopause standoff distance, predictions at locations impacted by the CME flank are those most affected by the specific CME model used.
We conclude by noting that, as both models employed in this work model CMEs as nonmagnetized plasma clouds inserted in the heliosphere without an internal magnetic field structure, the magnetic field is generally significantly underestimated compared to observations. Therefore, such estimates may provide just a lower limit for the CME geoeffectiveness, and the impact that a realistic magnetic field (e.g., obtained using a flux rope CME model) would have on the variation of the Kp index predictions cannot be assessed. However, prediction centers (e.g., the NOAA/Space Weather Prediction Centre and the UK MET office) still employ cone CME models to produce daily predictions of the space weather condition near Earth. Also, Kp forecasts formulated by employing cone CME models and empirical relations – similar to the one used in this work – are still widely used to validate model performances for both scientific and operational purposes (see, e.g., Mays et al., 2015). Therefore, we regard the discussion presented here as relevant given the current operational status and recent research publications.
4 Case study CME on 7 January 2014
In the following, we test the spheroidal CME model in EUHFORIA against a real CME event. As a case study, we choose the CME observed on 7 January 2014 and previously analyzed and modeled by Mays et al. (2015) using a tilted ellipsoid CME shape in the WSA–ENLIL+Cone model (e.g. Odstrcil et al., 2004).
Figure 1 shows the CME in the solar corona as seen by SOHO (representing the Earth viewpoint), STEREOA, and STEREOB (top row), and fitted with the cone (middle row) and croissant (bottom row) models. Both sets of reconstructions were performed using the GCS model and fitting tool, with the difference that for the cone case, the socalled halfwidth parameter (quantifying the halfangular separation between the CME legs) was set to zero.
We use the results from the GCS reconstructions shown in Figure 1 to construct a set of input parameters for the CME in both the cone and spheroid CME models as described in Section 2. The resulting set of CME input parameters used in EUHFORIA is provided in Table 2. While the values for latitude, longitude, and tilt are direct outputs of the GCS reconstruction, the angular width parameters are obtained from the GCSspecific halfwidth and aspect ratio (i.e., the ratio of the CME size at two orthogonal directions) using the formulation outlined by Palmerio et al. (2023, Appendix A). The CME speed is estimated by fitting the CME with the GCS model in two successive time frames (on 7 January 2014 at ∼19:10 and ∼19:40 UT) and by considering the corresponding heliocentric distance traveled by the CME apex (or nose). The CME time of crossing through the model inner boundary is determined by linearly extrapolating the obtained CME kinematics to the 21.5 R_{s} (0.1 au) boundary height. The ratio of the CME major radius to the minor radius corresponds to a spheroid aspect ratio of 2. This set of GCSbased CME input parameters for EUHFORIA are consistent with those reconstructed by Mays et al. (2015) for the same event. Similarly to Table 1, the bottom portion of Table 2 lists derived parameters such as the CME total volume (V_{CME}), mass (m_{CME}), and kinetic energy (E_{kin,CME}) in the two cases. In contrast to the fictitious event discussed in Section 3, however, the values for the spheroidal model in this event are 57% smaller than for the spherical (cone) model. Since each of these parameters likely has an influence on the global propagation of the CME front, we expect them to also impact the arrival time/speed of the CME and the resulting geoeffectiveness at a specific target location such as Earth (the latter one mainly due to the different dynamic pressure in the two EUHFORIA runs). We further discuss these aspects below.
Input parameters used to model the 7 January 2014 CME with the cone and spheroidal models.
To model the ambient solar wind, we drive the EUHFORIA coronal model using a single National Solar Observatory Global/Oscillation Network Group (NSO/GONG; Harvey et al., 1996) dailyupdated standard synoptic magnetogram taken on 7 January 2014 at 01:04 UT. The coronal model set up used is the same as in Pomoell & Poedts (2018). Simulations are performed over a computational domain extending from 0.1 au to 2 au in the radial (r) direction, and covering ±60° in the latitudinal (θ) direction and ±180° in the longitudinal (ϕ) direction. We employ a uniform grid with 256 cells in the radial direction and a resolution of 4° in the latitudinal and longitudinal directions.
It is important to note that both uncertainties in the solar wind and CME parameters can have significant effects on the predicted propagation and impact of CMEs at given target locations. For example, Riley & BenNun (2021) explored the sources of uncertainty in the predicted CME arrival time at Earth using a set of numerical MHD simulations of cone CMEs in different ambient solar wind backgrounds. They found that uncertainties in each of the CME initial parameters, such as longitude, latitude, width, and speed, introduce between 2.5 h and 7.5 h of uncertainty into the predicted CME arrival time at Earth. Furthermore, they concluded that the ambient solar wind structure was the largest source of such an uncertainty, and that without better constraints on the initial conditions of heliospheric simulations (i.e., from coronal models), it is likely that the CME arrival time error will remain close to ±10 h (Riley et al., 2018). Even more so, Ledvina et al. (2023) found that input photospheric (i.e., magnetogram maps) and coronal (i.e., magnetic field reconstructions) conditions in modeling CME propagation can have a profound effect on the structures predicted in situ even in the case of a simplified, uniform background solar wind. At the same time, the choice of CME input parameters has a direct effect on the resulting propagation and estimated impact(s). It has been shown that, even when employing a single propagation model and keeping the background solar wind conditions unchanged, CME properties derived from coronagraphic reconstructions performed by different users, or even by the same user assuming different CME morphologies, can result in CME predictions that vary significantly in timing and severity (Palmerio et al., 2022). Overall, it has been demonstrated that CME reconstructions performed independently may be expected to differ by 4.0° in the latitude, 8.0° in the longitude, 24.0° in the tilt, 9.3° in the angular width, and 115 km s^{−1} in the speed (Kay & Palmerio, 2024). In the context of the GCS model employed here, a sensitivity analysis performed by Thernisien et al. (2009) and a multiuser investigation conducted by Verbeke et al. (2023) reported typical uncertainties of a few degrees in the latitude and longitude (the former being slightly better constrained), ∼20° in the tilt, ∼10° in the halfwidth, and ∼0.5 R_{S} in the nose height. Such uncertainties are also compatible with those reported in Table 1 by Mays et al. (2015) for the specific CME considered here (i.e., the 7 January 2014 event), based on the GCS as well as other reconstruction methods. For this specific event, a typical ∼10° uncertainty in the halfwidth derived from the GCS reconstruction also implies the two shapes employed in the cone (spherical) and spheroidal EUHFORIA runs are clearly distinct from each other given their halfwidths differ for more than the typical uncertainty, at least along the CME minor axis direction. More generally, Kay & NievesChinchilla (2021) studied the sensitivity of CME input parameters in the framework of analytical modeling, and found that different properties tend to be more important in the cases of CMEs with different speeds, the most critical parameters being generally the mass, magnetic field strength, adiabatic index, and axial field profile. Hence, we caution against drawing strong conclusions from the results presented in this section, especially in light of the relatively large uncertainties involved in modeling the different domains of the Sun–Earth (or Sun–heliosphere) chain.
4.1 Comparison of CME models at 0.1 au
Figure 8 shows the modeling results on the spherical surface at r = 0.1 au on 7 January 2014 around 22:00 UT, just after the CME insertion in the heliospheric domain. Figure 8 demonstrates the difference in the CME front shape between the two runs, with the cone CME being more extended than the spheroid one. In the phases immediately after their insertion in the heliosphere, both CMEs are propagating tangentially to the Sun–Earth direction.
Figure 8 3D views of the radial velocity (v_{r}, in units of km s^{−1}) at r = 0.1 au on 7 January 2014 around 22:00 UT (during the CME insertion), for the two simulations performed. (a) Cone, (b) Spheroid. The green lines mark the central meridian and the solar equator. The area where the CME is being inserted is visible as a red patch is the southwestern quadrant. 
We note that, in contrast to the fictitious CME case (see Fig. 3), for this event the crosssection of the spheroidal CME intersects the solar equatorial plane significantly less than the spherical one. Even more so, while the cone CME crosses Earth’s projected position with its northeastern edge, the spheroid does not feature a component along the Sun–Earth line at 0.1 au. Nevertheless, an encounter at Earth may still take place due to processes such as CME expansion in the heliosphere and interactions with the structured solar wind (see, e.g., the reviews of Manchester et al., 2017; Luhmann et al., 2020, on the interplanetary evolution of CMEs).
4.2 Comparison of CME models in the heliosphere
Figure 9 shows the modeling results in the equatorial and meridional planes around 07:00 UT on 9 January 2014, when the CME nose was close to 1 au in the two simulations. Figure 9 shows that the cone CME is generally more extended and faster than the spheroid CME in both the ecliptic plane (Figs. 9a and 9c) and in the meridional one (Figs. 9b and 9d) at this stage of propagation. While close to the CME nose this is most likely caused by the higher mass of the cone CME compared to the spheroidal CME, the differences near the CME flanks are primarily due to the different geometry employed in the two runs. Furthermore, contrary to the situation depicted in Figure 8 at the model inner boundary, both runs show that the CME has a propagation component along the Sun–Earth direction, most likely due to the expansion of its shock front during the initial propagation phases in the heliospheric domain.
Figure 9 Modeling results for the 7 January 2014 CME with the cone and spheroidal CME models on 9 January 2014 around 07:00 UT. The format is the same as Figure 4. The green lines mark the Sun–Earth direction. 
4.3 Comparison of CME models at Earth
Figure 10 compares the EUHFORIA time series at Earth in the two simulations with real observations from the Wind (Ogilvie & Desch, 1997) spacecraft, specifically from the Solar Wind Experiment (SWE; Ogilvie et al., 1995) instrument.
Figure 10 EUHFORIA predictions at Earth for the 7 January 2014 CME. (a) Radial speed, (b) Proton number density. Cone and spheroidal model results are shown in teal and violet, respectively. As a reference, ambient solar wind model results are plotted in gray. Real measurements from Wind are shown in black. Vertical dashed lines mark the CME arrival time in the two EUHFORIA runs and in Wind observations. The CME impact speed and density in EUHFORIA runs is marked by the dots. 
The CME is observed to arrive at Wind on 9 January at 19:40 UT, when its interplanetary shock was detected (from the IPshocks database^{2}; Kilpua et al., 2015). The shock was followed by a sheath region and by the CME ejecta, which started on January 10 at about 07:00 UT and ended on January 11 at about 03:00 UT as indicated by magnetic field rotations and the presence of bidirectional suprathermal electron flux (not shown). As shown in Figure 10, EUHFORIA modeled the CME to arrive at Earth on 9 January at 07:41 UT using the spherical (cone) model, and on 9 January at 17:03 UT using the spheroidal model. Using the input parameters described above, the spheroidal model thus performs significantly better than the cone in predicting the CME arrival at Earth (Δt_{err} = t_{predicted} − t_{observed} = −11.97 h for the cone, compared to Δt_{err} = −2.6 h for the spheroid). Furthermore, we note that the relative difference in the arrival time between the two runs at Earth is of almost 10 h, even though the two CMEs were launched with the same initial speed in the two simulations. This result clearly demonstrates how the choice of CME model can have an effect on CME arrival time predictions that is at least comparable to current prediction capabilities (estimated to be around ±10 h; Riley et al., 2018; Vourlidas et al., 2019).
In Wind observations, the CME had an impact peak speed of about 480 km s^{−1} and a peak proton density of about 12 cm^{−3}. Figure 10 shows that the spherical (cone) model predicted a CME peak speed of 617 km s^{−1}, and a peak proton density of 22 cm^{−3}. These values are significantly larger than the observed values (about 140 km s^{−1} and 10 cm^{−3} for the impact speed and density, respectively). On the other hand, the spheroidal model was able to predict the CME impact speed and density remarkably well with respect to the observed values. In this case, the modeled CME peak speed was 494 km s^{−1}s, and the peak proton density was 16 cm^{−3}, which are only about 15 km s^{−1} and 4 cm^{−3} higher than observed. For comparison, Mays et al. (2015) report that the uncertainties in the CME initial parameters affect the CME arrival time at Earth by up to +9/−7 h compared to the mean arrival time predicted by the WSA–ENLIL+Cone model. However, in their simulations the mean predicted arrival time was shifted by about 18 h compared to the observed arrival time at Earth. Additionally, the CME impact speed and density were overpredicted by 100–700 km s^{−1} and by a factor 3–5 compared to the observed values, respectively. With respect to these metrics, the EUHFORIA spheroidal CME simulation presented in this work therefore performs significantly better. It is however important to emphasize that the spheroidal CME simulation performed in this work only represents a single realization: while we can expect similar uncertainties in the predictions at Earth as those reported by Mays et al. (2015), more extensive ensemble simulations should be performed in order to determine the actual prediction performances for a broad set of ambient solar wind realizations and CME initial parameters.
Although it remains challenging to disentangle the effects of different CME geometries from that of different CME masses and kinetic energies, the latter are most likely dominant close to the CME nose, while geometrical effects are expected to be predominant close to the flanks (i.e., including at Earth location). We caution that future targeted validation efforts are required to precisely pinpoint the magnitude and localization of these effects across CME fronts, for different combinations of CME and solar wind parameters.
4.4 Comparison of CME models: prediction of CME geoeffectiveness
As a final step, we compute the predicted CME geoeffectiveness (in terms of magnetopause compression and Kp index, as in Sect. 3.4) based on EUHFORIA time series at Earth for both the cone (spherical) and spheroidal simulations, and we compare them with predictions obtained from Wind time series, as well as with actual measurements of the Kp index (from the GFZ – Helmholtz Center Potsdam^{3}) at 3hour cadence. Results are shown in Figure 11.
Figure 11 Modeled and observed CME geoeffectiveness as a function of time. (a) Magnetopause standoff distance calculated from EUHFORIA and Wind time series using the model by Shue et al. (1997). The black horizontal line marks the geosynchronous orbits at 6.6 R_{e}. Cone and spheroidal model results are shown in teal and violet, respectively. As a reference, ambient solar wind model results are plotted in gray. Estimates from Wind measurements are shown in black. (b) Kp index calculated from EUHFORIA and Wind time series using the relation by Newell et al. (2008). Threehour Kp measurements are shown in red. 
In terms of the CMEinduced compression of the Earth’s magnetopause (Fig. 11a), the cone simulation predicts the dso to reach a minimum value of about 5.6 R_{e}, while the spheroid simulation predicts it to only reach a minimum of 6.3 R_{e}. Both predictions are below geosynchronous orbits, thus in this case, the two EUHFORIA simulations agree in predicting that satellites in geosynchronous orbits may become exposed to solar wind conditions during the CME impact. However, the time evolution of the dso predicted by the two simulations is quite different, as expected from the different CME arrival times and impact parameters predicted. The spheroidal model, in particular, is the one that matches better with dso predictions based on Wind data. This is true both for the dso temporal evolution, as well as for the predicted minimum dso (6.6 R_{e} based on Wind data, compared to 6.3 R_{e} from the spheroidal EUHFORIA simulation).
The two EUHFORIA simulations also predict significantly different space weather impacts in terms of geomagnetic storm strength based on the Kp index (Fig. 11b). The cone simulation predicts a maximum Kp of 6.3 (corresponding to a moderate storm according to NOAA levels), while the spheroidal simulation predicts a maximum Kp of only 3.3 (corresponding to a quiet state). For comparison, the Wind time series predicts a maximum Kp of around 3.9, and the observed maximum Kp associated with the event was classified as 3 (corresponding to quiet conditions) by NOAA. As predictions from the cone model result in a dramatic overestimation of the CME geoeffectiveness for this particular event, this case clearly demonstrates how an appropriate choice of CME model in heliospheric simulations can provide a significant improvement to the predicted CME geoeffectiveness based on the Kp index. Finally, we note that the measured Kp for this specific event is very well reproduced by applying the solar wind–coupling function by Newell et al. (2008) empirical relation to Wind time series. This gives us additional confidence that this empirical relation, despite dating back more than a decade, is doing a rather good job in capturing the main coupling between the solar wind and geomagnetic response in terms of the Kp for this specific event.
Ultimately disentangling the impact of different CME geometries from that of different CME masses and kinetic energies on the resulting geoeffectiveness remains challenging without performing adhoc, extensive validation studies. However, we note here as well that differences in the CME masses and kinetic energies are most likely dominant close to the CME nose, while geometrical effects are expected to be predominant close to the flanks (i.e., including Earth’s location).
5 Conclusions
In this paper, we have described the implementation and initial validation of the spheroid CME model in the 3D MHD EUHFORIA code. After presenting the mathematical description of the spheroidal structure, we have shown a comparison of EUHFORIA simulation runs employing the traditional cone (spherical) ejecta and the newlyimplemented spheroidal one for an idealized CME, as well as for an event that erupted on 7 January 2014 and that was later observed at Earth as a flank encounter. We have demonstrated that, at least under the specific conditions considered here, the new spheroidal implementation provided an improvement compared to the traditional cone model.
In the first validation experiment (i.e., the one with fictitious CMEs), we have found that differences in CME arrival time and speed between the two models become more prominent with increasing distance from the CME nose. As a result, the expected CME geoeffectiveness (evaluated in terms of the Kp index and magnetopause standoff distance) estimated from the two models’ results at 1 au also becomes more different as the distance from the CME nose increases. Furthermore, we have demonstrated the choice of CME geometry not only affects predictions of the CME arrival time and impact properties, but can also impact the predictions of the CME hit/miss at spacecraft impacted by the CME flank extremities. Since the fictitious CME modeled represents an average fast CME in terms of speed and angular width, despite the idealised numerical setup employed, one can expect the results to apply to a potentially large set of real CME events.
In the second validation test (i.e., the one with the real CME), we have found that a spheroidal CME structure is able to reproduce a better arrival time and speed of the observed flank encounter at Earth over the cone ejecta. Specifically, we noted that the CME arrival time and speed differed by ∼10 h and ∼120 km s^{−1}, respectively, for the two CMEs that were launched with the same velocity, mass density, and temperature, and only differed for the morphology of their crosssection. The predicted geoeffectiveness was also highly affected by the choice of the CME model. Particularly the prediction of the induced geomagnetic storm strength varied from a quiet state (Kp of 3 for the spheroid CME model, in remarkable agreement with actual Kp measurements at Earth) to a moderate storm (Kp above 6 for the cone CME model). This event thus clearly demonstrated how an appropriate description of the CME flanks in heliospheric simulations can provide a significant improvement to the predicted CME geoeffectiveness, particularly in the case flank encounters.
Although further investigations are needed in the future to quantify exactly in how many cases and to which extent the spheroidal model performs better than the cone model, the results presented in this work clearly demonstrate the potential of the new spheroid approach in providing improved space weather forecasts, particularly with respect to flank CME encounters. Employing a spheroidal CME morphology in 3D MHD heliospheric simulations can constitute a significant improvement over the traditional cone model, while maintaining the physics and assumptions of the overall simulation – and, thus, the computational time – fairly unchanged. This is especially advantageous in the context of realtime space weather forecasts, where magnetized CMEs are currently not (yet) employed. We note that the WSA–ENLIL+Cone model, currently used for operations by several forecasting agencies, has also implemented a socalled ellipsoidal geometry for CMEs (as the one showcased in the work of Mays et al., 2015), which is however only employed in research applications at the time of writing (Odstrcil, 2023]. One source of difficulty may be related to the adopted coronal reconstruction technique to derive CME input parameters (e.g., the “simpler” SWPCCAT versus the “more complex” GCS). Nevertheless, it has been shown that CMEs can be relatively easily be fitted with ellipses in coronagraph imagery even when a single viewpoint is available (e.g., Yurchyshyn et al., 2007), and the method could be extended to two or more simultaneous observers without recurring to geometrical descriptions characterized by many free parameters. We also note that this work presents a first validation of the performance of the spheroid CME model versus the traditional cone CME model without considering the effect of the uncertainties affecting the determination of the initial CME parameters. Thus, while our results suggest this new CME implementation may perform better than the traditional cone model under specific conditions, and while a general improvement compared to the traditional cone model is also expected, whether this is indeed achieved remain to be proven in future studies.
In conclusion, a larger validation study to evaluate the advantages of employing a spheroidal versus a spherical CME model will provide more robust benchmarking, but the results presented in this work show that the spheroid approach is a promising avenue to pursue toward the improvement of current space weather forecasting capabilities.
Acknowledgments
EUHFORIA is developed as a joint effort between the University of Helsinki and KU Leuven. The simulations were carried out at the VSC – Flemish Supercomputer Center, funded by the Hercules foundation and the Flemish Government – Department EWI. EP acknowledges support from NASA’s HTMS (grant no. 80NSSC20K1274) and LWSSC (grant no. 80NSSC22K0893) programs as well as NSF’s PREEVENTS (grant no. ICER1854790) program. The editor thanks two anonymous reviewers for their assistance in evaluating this paper.
References
 Asvestari E, Pomoell J, Kilpua E, Good S, Chatzistergos T, Temmer M, Palmerio E, Poedts S, Magdalenic J. 2021. Modelling a multispacecraft coronal mass ejection encounter with EUHFORIA. A&A 652: A27. https://doi.org/10.1051/00046361/202140315. [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Cremades H, Bothmer V. 2004. On the threedimensional configuration of coronal mass ejections. A&A 422: 307–322. https://doi.org/10.1051/00046361:20035776. [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fisher RR, Munro RH. 1984. Coronal transient geometry. I – The flareassociated event of 1981 March 25. Astrophys J 280: 428–439. https://doi.org/10.1086/162009. [CrossRef] [Google Scholar]
 Gopalswamy N. 2006. Coronal mass ejections of solar cycle 23. J Astrophys Astron 27: 243–254. https://doi.org/10.1007/BF02702527. [CrossRef] [Google Scholar]
 Harvey JW, Hill F, Hubbard RP, Kennedy JR, Leibacher JW, et al. 1996. The Global Oscillation Network Group (GONG) project. Science 272(5266): 1284–1286. https://doi.org/10.1126/science.272.5266.1284. [CrossRef] [Google Scholar]
 Isavnin A. 2016. FRiED: A novel threedimensional model of coronal mass ejections. Astrophys J 833(2): 267. https://doi.org/10.3847/15384357/833/2/267. [CrossRef] [Google Scholar]
 Jang S, Moon YJ, Kim RS, Lee H, Cho KS. 2016. Comparison between 2D and 3D parameters of 306 frontside halo CMEs from 2009 to 2013. Astrophys J 821(2): 95. https://doi.org/10.3847/0004637X/821/2/95. [CrossRef] [Google Scholar]
 Kataoka R, Ebisuzaki T, Kusano K, Shiota D, Inoue S, Yamamoto TT, Tokumaru M. 2009. Three dimensional MHD modeling of the solar wind structures associated with 13 December 2006 coronal mass ejection. J Geophys Res 114(A10): A10102. https://doi.org/10.1029/2009JA014167. [Google Scholar]
 Kay C, NievesChinchilla T. 2021. Modeling interplanetary expansion and deformation of CMEs with ANTEATR PARADE: relative contribution of different forces. J Geophys Res Space Phys 126(5): 2020JA028911. https://doi.org/10.1029/2020JA028911. [CrossRef] [Google Scholar]
 Kay C, Palmerio E. 2024. Collection, collation, and comparison of 3d coronal cme reconstructions. Space. Weather 22(1): e2023SW003796. https://doi.org/10.1029/2023SW003796. [CrossRef] [Google Scholar]
 Kilpua EKJ, Lugaz N, Mays ML, Temmer M. 2019. Forecasting the structure and orientation of earthbound coronal mass ejections. Space Weather 17(4): 498–526. https://doi.org/10.1029/2018SW001944. [NASA ADS] [CrossRef] [Google Scholar]
 Kilpua EKJ, Lumme E, Andreeova K, Isavnin A, Koskinen HEJ. 2015. Properties and drivers of fast interplanetary shocks near the orbit of the Earth (1995–2013). J Geophys Res Space Phys 120(6): 4112–4125. https://doi.org/10.1002/2015JA021138. [CrossRef] [Google Scholar]
 Lamy PL, Floyd O, Boclet B, Wojak J, Gilardy H, Barlyaeva T. 2019. Coronal mass ejections over solar cycles 23 and 24. Space Sci Rev 215(5): 39. https://doi.org/10.1007/s112140190605y. [CrossRef] [Google Scholar]
 Ledvina VE, Palmerio E, Kay C, AlHaddad N, Riley P. 2023. Modeling CME encounters at parker solar probe with OSPREI: dependence on photospheric and coronal conditions. A&A 673: A96. https://doi.org/10.1051/00046361/202245445. [CrossRef] [EDP Sciences] [Google Scholar]
 Lockwood M. 2022. Solar wind – magnetosphere coupling functions: pitfalls, limitations, and applications. Space Weather 20(2): e2021SW002989. https://doi.org/10.1029/2021SW002989. [CrossRef] [Google Scholar]
 Luhmann JG, Gopalswamy N, Jian LK, Lugaz N. 2020. ICME evolution in the inner heliosphere. Sol Phys 295(4): 61. https://doi.org/10.1007/s11207020016240. [CrossRef] [Google Scholar]
 Luo B, Liu S, Gong J. 2017. Two empirical models for shortterm forecast of Kp. Space Weather 15(3): 503–516. https://doi.org/10.1002/2016SW001585. [CrossRef] [Google Scholar]
 Maharana A, Isavnin A, Scolini C, Wijsen N, Rodriguez L, Mierla M, Magdalenić J, Poedts S. 2022. Implementation and validation of the FRi3D flux rope model in EUHFORIA. Adv Space Res 70(6): 1641–1662. https://doi.org/10.1016/j.asr.2022.05.056. [CrossRef] [Google Scholar]
 Manchester W, Kilpua EKJ, Liu YD, Lugaz N, Riley P, Török T, Vršnak B. 2017. The physical processes of CME/ICME evolution. Space Sci Rev 212(3–4): 1159–1219. https://doi.org/10.1007/s1121401703940. [CrossRef] [Google Scholar]
 Mays ML, Thompson BJ, Jian LK, Colaninno RC, Odstrcil D, et al. 2015. Propagation of the 7 January 2014 CME and resulting geomagnetic nonevent. Astrophys J 812(2): 145. https://doi.org/10.1088/0004637X/812/2/145. [CrossRef] [Google Scholar]
 Millward G, Biesecker D, Pizzo V, de Koning CA. 2013. An operational software tool for the analysis of coronagraph images: determining CME parameters for input into the WSAEnlil heliospheric model. Space Weather 11(2): 57–68. https://doi.org/10.1002/swe.20024. [NASA ADS] [CrossRef] [Google Scholar]
 Möstl C, Rollett T, Frahm RA, Liu YD, Long DM, et al. 2015. Strong coronal channelling and interplanetary evolution of a solar storm up to Earth and Mars. Nat. Commun. 6: 7135. https://doi.org/10.1038/ncomms8135. [CrossRef] [Google Scholar]
 Newell PT, Sotirelis T, Liou K, Rich FJ. 2008. Pairs of solar windmagnetosphere coupling functions: Combining a merging term with a viscous term works best. J Geophys Res 113(A4): A04218. https://doi.org/10.1029/2007JA012825. [Google Scholar]
 Odstrcil D. 2003. Modeling 3D solar wind structure. Adv Space Res 32(4): 497–506. https://doi.org/10.1016/S02731177(03)003326. [CrossRef] [Google Scholar]
 Odstrcil D. 2023. Heliospheric 3D MHD ENLIL simulations of multiCME and multispacecraft events. Front. Astron. Space Sci. 10: 1226992. https://doi.org/10.3389/fspas.2023.1226992. [CrossRef] [Google Scholar]
 Odstrcil D, Riley P, Zhao XP. 2004. Numerical simulation of the 12 May 1997 interplanetary CME event. J Geophys Res 109(A2): A02116. https://doi.org/10.1029/2003JA010135. [Google Scholar]
 Ogilvie KW, Chornay DJ, Fritzenreiter RJ, Hunsaker F, Keller J, et al. 1995. SWE, a comprehensive plasma instrument for the wind spacecraft. Space Sci Rev 71(1–4): 55–77. https://doi.org/10.1007/BF00751326. [CrossRef] [Google Scholar]
 Ogilvie KW, Desch MD. 1997. The wind spacecraft and its early scientific results. Adv Space Res 20: 559–568. https://doi.org/10.1016/S02731177(97)004390. [CrossRef] [Google Scholar]
 Palmerio E, Kilpua EKJ, James AW, Green LM, Pomoell J, Isavnin A, Valori G. 2017. Determining the intrinsic CME flux rope type using remotesensing solar disk observations. Sol Phys 292(2): 39. https://doi.org/10.1007/s112070171063x. [CrossRef] [Google Scholar]
 Palmerio E, Lee CO, Mays ML, Luhmann JG, Lario D, et al. 2022. CMEs and SEPs during November–December 2020: a challenge for realtime space weather forecasting. Space Weather 20(5): e2021SW002993. https://doi.org/10.1029/2021SW002993. [NASA ADS] [Google Scholar]
 Palmerio E, Maharana A, Lynch BJ, Scolini C, Good SW, Pomoell J, Isavnin A, Kilpua EKJ. 2023. Modeling a coronal mass ejection from an extended filament channel. II. Interplanetary propagation to 1 au. Astrophys J 958(1): 91. https://doi.org/10.3847/15384357/ad0229. [CrossRef] [Google Scholar]
 Palmerio E, Scolini C, Barnes D, Magdalenić J, West MJ, et al. 2019. Multipoint study of successive coronal mass ejections driving moderate disturbances at 1 au. Astrophys J 878(1): 37. https://doi.org/10.3847/15384357/ab1850. [CrossRef] [Google Scholar]
 Pizzo V, Millward G, Parsons A, Biesecker D, Hill S, Odstrcil D. 2011. WangSheeleyArgeEnlil cone model transitions to operations. Space Weather 9(3): 03004. https://doi.org/10.1029/2011SW000663. [Google Scholar]
 Pomoell J, Poedts S. 2018. EUHFORIA: European heliospheric forecasting information asset. J Space Weather Space Clim 8: A35. https://doi.org/10.1051/swsc/2018020. [Google Scholar]
 Riley P, BenNun M. 2021. On the sources and sizes of uncertainty in predicting the arrival time of interplanetary coronal mass ejections using global MHD models. Space Weather 19(6): e02775. https://doi.org/10.1029/2021SW002775. [CrossRef] [Google Scholar]
 Riley P, Mays ML, Andries J, Amerstorfer T, Biesecker D, et al. 2018. Forecasting the arrival time of coronal mass ejections: analysis of the CCMC CME scoreboard. Space Weather 16(9): 1245–1260. https://doi.org/10.1029/2018SW001962. [CrossRef] [Google Scholar]
 Scolini C, Chané E, Temmer M, Kilpua EKJ, Dissauer K, et al. 2020. CME–CME interactions as sources of CME geoeffectiveness: the formation of the complex ejecta and intense geomagnetic storm in 2017 early September. Astrophys J Suppl 247(1): 21. https://doi.org/10.3847/15384365/ab6216. [CrossRef] [Google Scholar]
 Scolini C, Verbeke C, Poedts S, Chané E, Pomoell J, Zuccarello FP. 2018a. Effect of the initial shape of coronal mass ejections on 3D MHD simulations and geoeffectiveness predictions. Space Weather 16(6): 754–771. https://doi.org/10.1029/2018SW001806. [CrossRef] [Google Scholar]
 Scolini C, Messerotti M, Poedts S, Rodriguez L. 2018b. Halo coronal mass ejections during solar cycle 24: reconstruction of the global scenario and geoeffectiveness. J Space Weather Space Clim 8: A9. https://doi.org/10.1051/swsc/2017046. [Google Scholar]
 Scolini C, Rodriguez L, Mierla M, Pomoell J, Poedts S. 2019. Observationbased modelling of magnetised coronal mass ejections with EUHFORIA. A&A 626: A122. https://doi.org/10.1051/00046361/201935053. [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Scolini C, Winslow RM, Lugaz N, Poedts S. 2021. Evolution of interplanetary coronal mass ejection complexity: a numerical study through a swarm of simulated spacecraft. Astrophys J Lett. 916(2): L15. https://doi.org/10.3847/20418213/ac0d58. [CrossRef] [Google Scholar]
 Scolini C, Winslow RM, Lugaz N, Poedts S. 2023. Characteristic scales of complexity and coherence within interplanetary coronal mass ejections: insights from spacecraft swarms in global heliospheric simulations. Astrophys J 944(1): 46. https://doi.org/10.3847/15384357/aca893. [CrossRef] [Google Scholar]
 Shen F, Liu Y, Yang Y. 2021a. Numerical research on the effect of the initial parameters of a CME fluxrope model on simulation results. Astrophys J Suppl 253(1): 12. https://doi.org/10.3847/15384365/abd4d2. [CrossRef] [Google Scholar]
 Shen F, Liu Y, Yang Y. 2021b. Numerical research on the effect of the initial parameters of a CME fluxrope model on simulation results. II. Different locations of observers. Astrophys J 915(1): 30. https://doi.org/10.3847/15384357/ac004e. [CrossRef] [Google Scholar]
 Shue JH, Chao JK, Fu HC, Russell CT, Song P, Khurana KK, Singer HJ. 1997. A new functional form to study the solar wind control of the magnetopause size and shape. J Geophys Res 102(A5): 9497–9512. https://doi.org/10.1029/97JA00196. [CrossRef] [Google Scholar]
 Singh T, Yalim MS, Pogorelov NV, Gopalswamy N. 2020. A modified spheromak model suitable for coronal mass ejection simulations. Astrophys J 894(1): 49. https://doi.org/10.3847/15384357/ab845f. [CrossRef] [Google Scholar]
 Thernisien A. 2011. Implementation of the graduated cylindrical shell model for the threedimensional reconstruction of coronal mass ejections. Astrophys J Suppl. 194(2): 33. https://doi.org/10.1088/00670049/194/2/33. [CrossRef] [Google Scholar]
 Thernisien A, Vourlidas A, Howard RA. 2009. Forward modeling of coronal mass ejections using STEREO/SECCHI data. Sol Phys 256(1–2): 111–130. https://doi.org/10.1007/s1120700993465. [CrossRef] [Google Scholar]
 Verbeke C, Mays ML, Kay C, Riley P, Palmerio E, et al. 2023. Quantifying errors in 3D CME parameters derived from synthetic data using whitelight reconstruction techniques. Adv Space Res 72(12): 5243–5262. https://doi.org/10.1016/j.asr.2022.08.056. [CrossRef] [Google Scholar]
 Vourlidas A, Lynch BJ, Howard RA, Li Y. 2013. How many CMEs have flux ropes? Deciphering the signatures of shocks, flux ropes, and prominences in coronagraph observations of CMEs. Sol Phys 284(1): 179–201. https://doi.org/10.1007/s1120701200848. [Google Scholar]
 Vourlidas A, Patsourakos S, Savani NP. 2019. Predicting the geoeffective properties of coronal mass ejections: current status, open issues and path forward. Phil Trans R Soc London Ser A 377(2148): 20180096. https://doi.org/10.1098/rsta.2018.0096. [Google Scholar]
 Yurchyshyn V, Hu Q, Lepping RP, Lynch BJ, Krall J. 2007. Orientations of LASCO Halo CMEs and their connection to the flux rope structure of interplanetary CMEs. Adv Space Res 40(12): 1821–1826. https://doi.org/10.1016/j.asr.2007.01.059. [CrossRef] [Google Scholar]
 Zhang M, Feng XS, Yang LP. 2019. Threedimensional MHD simulation of the 2008 December 12 coronal mass ejection: from the Sun to Interplanetary space. J Space Weather Space Clim 9: A33. https://doi.org/10.1051/swsc/2019034. [CrossRef] [EDP Sciences] [Google Scholar]
 Zhou Y, Feng X, Zhao X. 2014. Using a 3D MHD simulation to interpret propagation and evolution of a coronal mass ejection observed by multiple spacecraft: the 3 April 2010 event. J Geophys Res Space Phys 119(12): 9321–9333. [CrossRef] [Google Scholar]
Cite this article as: Scolini C & Palmerio E. 2024. The spheroid CME model in EUHFORIA. J. Space Weather Space Clim. 14, 13. https://doi.org/10.1051/swsc/2024011.
All Tables
CME input parameters used for the spherical (cone) and spheroidal models to initialize the fictitious CME cases.
Input parameters used to model the 7 January 2014 CME with the cone and spheroidal models.
All Figures
Figure 1 Example of a CME (from 7 January 2014) fitted in the solar corona using the cone and croissant models. The top row (a–c) shows whitelight data from three viewpoints (Earth/SOHO, STEREOA, and STEREOB), the middle row (d–f) shows the same set of data with the cone wireframe (in green) overlaid, and the bottom row (g–i) shows the same set of data with the croissant wireframe (in magenta) overlaid. Both sets of reconstructions were performed using the GCS model, with the difference that for the cone case the socalled halfwidth parameter (quantifying the halfangular separation between the CME legs) was set to zero. 

In the text 
Figure 2 3D visualization of the (a) spherical and (b) spheroidal CME shapes. x′, y′, and z′ constitute a local coordinate system centered at the center of the CME. Dimensions are normalized to 1 for convenience. The spherical CME has circular crosssections in the directions perpendicular to x′, y′, and z′, as shown by the 2D shadow projections. The spheroidal CME has a circular crosssections only in the direction perpendicular to the major axis (parallel to y′), while it has elliptical crosssections in the directions perpendicular to the minor axes (parallel to x′ and y′), as shown by the 2D shadow projections. 

In the text 
Figure 3 3D views of the radial velocity (v_{r}, in units of km s^{−1}) and number density (n, in units of cm^{−3}) at the heliospheric inner boundary during the CME insertion, for the two simulations performed. (a) v_{r} for the cone (spherical) CME model, (b) n for the cone (spherical) CME model, (c) v_{r} for the spheroidal CME model, (d) n for the spheroidal CME model. The area where the CME is being inserted (v_{r} = 800 km s^{−3}, n = 600 cm^{−3}) is visible as the red (orange) patch on the left (right). 

In the text 
Figure 4 Modeling results for the fictitious CME with the sphere and spheroidal CME models. The panels show the radial speed v_{r} in the equatorial (left) and meridional (right) plane on 20230101 05:57 UT when the CME nose was close to 1 au in both simulations. Representative scaled number density nr^{2} values are indicated by the gray contour lines. (a, b) spherical (cone) CME model, (c, d) spheroidal CME model. The green lines marks directions to representative virtual spacecraft located at θ = 0° and ϕ between −90° and 90°, with 10° increments. 

In the text 
Figure 5 Radial speed and number density recorded by virtual spacecraft at 1 au for the cone (teal) and spheroid (violet) CME models. As a reference, the ambient solar wind values are plotted in gray. Locations are shown at ϕ values of (a) 0°, (b) 10°, (c) 20°, (d) 30°, (e) 40°, and (f) 50°. Vertical dashed lines mark the CME arrival time in the two runs. Dots mark the CME impact (peak) speed and density. 

In the text 
Figure 6 Difference between the spheroid and sphere (cone) CME impact properties at 1 au, as a function of the spacecraft longitude. (a) Difference in the CME arrival time, (b) Difference in the CME impact (peak) speed, (c) Difference in the CME impact (peak) density. The orange lines are calculated by subtracting the sphere to the spheroid properties. The teal (violet) lines show the subtracted values between the sphere (spheroid) and the ambient solar wind. 

In the text 
Figure 7 CME geoeffectiveness prediction as a function of the spacecraft longitude. (a) Estimate of the minimum magnetopause standoff distance calculated using the model by Shue et al. (1997). The black horizontal line marks the geosynchronous orbits at 6.6 R_{e}. (b) Estimate of the maximum Kp index calculated using the relation by Newell et al. (2008). 

In the text 
Figure 8 3D views of the radial velocity (v_{r}, in units of km s^{−1}) at r = 0.1 au on 7 January 2014 around 22:00 UT (during the CME insertion), for the two simulations performed. (a) Cone, (b) Spheroid. The green lines mark the central meridian and the solar equator. The area where the CME is being inserted is visible as a red patch is the southwestern quadrant. 

In the text 
Figure 9 Modeling results for the 7 January 2014 CME with the cone and spheroidal CME models on 9 January 2014 around 07:00 UT. The format is the same as Figure 4. The green lines mark the Sun–Earth direction. 

In the text 
Figure 10 EUHFORIA predictions at Earth for the 7 January 2014 CME. (a) Radial speed, (b) Proton number density. Cone and spheroidal model results are shown in teal and violet, respectively. As a reference, ambient solar wind model results are plotted in gray. Real measurements from Wind are shown in black. Vertical dashed lines mark the CME arrival time in the two EUHFORIA runs and in Wind observations. The CME impact speed and density in EUHFORIA runs is marked by the dots. 

In the text 
Figure 11 Modeled and observed CME geoeffectiveness as a function of time. (a) Magnetopause standoff distance calculated from EUHFORIA and Wind time series using the model by Shue et al. (1997). The black horizontal line marks the geosynchronous orbits at 6.6 R_{e}. Cone and spheroidal model results are shown in teal and violet, respectively. As a reference, ambient solar wind model results are plotted in gray. Estimates from Wind measurements are shown in black. (b) Kp index calculated from EUHFORIA and Wind time series using the relation by Newell et al. (2008). Threehour Kp measurements are shown in red. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.