The 17 April 2021 widespread solar energetic particle event

Context. A complex and long-lasting solar eruption on 17 April 2021 produced a widespread Solar Energetic Particle (SEP) event that was observed by ﬁve longitudinally well-separated observers in the inner heliosphere covering distances to the Sun from 0.42 to 1 au: BepiColombo, Parker Solar Probe, Solar Orbiter, STEREO A, and near-Earth spacecraft. The event was the second widespread SEP event dected in solar cycle 25 and produced relativistic electrons and protons. It was associated with a long-lasting solar hard X-ray ﬂare showing multiple hard X-ray peaks over a duration of one hour. The event was further accompanied by a medium fast Coronal Mass Ejection (CME) with a speed of 880 km s − 1 driving a shock, an EUV wave as well as long-lasting and complex radio burst activity showing four distinct type III burst groups over a period of 40 minutes. Aims. We aim at understanding the reason for the the wide spread of elevated SEP intensities in the inner heliosphere as well as identifying the underlying source regions of the observed energetic electrons and protons. Methods. A comprehensive multi-spacecraft analysis of remote-sensing observations and in-situ measurements of the energetic particles and interplanetary context is applied to attribute the SEP observations at the diﬀerent locations to the various potential source regions at the Sun. An ENLIL simulation is used to characterize the complex interplanetary state and its role for the energetic particle transport. The magnetic connection between each spacecraft and the Sun is determined using ballistic backmapping in combination with potential ﬁeld source surface extrapolations in the lower corona. In combination with a reconstruction of the coronal shock front we then determine the times when the shock establishes magnetic connections with the diﬀerent observers. Radio observations are used to characterize the directivity of the four main injection episodes, which are then employed in a 2D SEP transport simulation to test the importance of these diﬀerent injection episodes. Results. A comprehensive timing analysis of the inferred solar injection times of the SEPs observed at each spacecraft suggests diﬀerent source processes being important for the electron and the proton event. Comparison among the characteristics and timing of the potential particle sources, such as the CME-driven shock or the ﬂare, suggests a stronger shock contribution for the proton event and a more likely ﬂare-related source of the electron event. Conclusions. In contrast to earlier studies on widespread SEP events, we ﬁnd that in this event an important ingredient for the wide SEP spread was the wide longitudinal range of about 110 ◦ covered by distinct SEP injections, which is also supported by our SEP transport modeling.


Introduction
Solar energetic particle (SEP) events are characterized by a rich and complex set of physical processes responsible for the acceleration and propagation of the particles.Since the early observations of Forbush (1946), an enormous amount of knowledge has been built around SEPs, highlighting their importance for understanding the behavior of the outer layers of the Sun's atmosphere, as well as addressing fundamental questions related to energetic particle propagation in astrophysical environments (e.g., Reames 2021).Multipoint observations of SEP events at different heliospheric locations provide an invaluable opportunity to study the production and transport of energetic particles, with several recent studies addressing the problem from a variety of perspectives (e.g., Dresing et al. 2014; Gómez-Herrero locations in the heliographic equatorial plane together with nominal Parker field lines connecting each observer with the Sun depicted in the center of the plot.The black arrow marks the longitude of the associated flare (identified using Solar Orbiter STIX measurements as described in Sect.4.1), and the dashed black spiral denotes the nominal magnetic field line connecting to this location.BepiColombo (yellow) was the spacecraft with the best nominal connection to the flare site.Parker Solar Probe (purple) and Solar Orbiter (blue) were approximately equally separated in longitude from the flaring region, but on different sides.STEREO A (red) and Earth (black) were further separated to the west of the flare.Despite the large angular separation between all spacecraft, the SEP event was observed at these five locations as shown in Fig. 1 (right), which makes it a widespread event (e.g., Dresing et al. 2014).The top panel shows ∼1 MeV electron intensities and the bottom panel ∼25 MeV proton intensities.As expected due to its closest magnetic connection, BepiColombo observed the highest intensities.Parker Solar Probe, being situated closest at 0.42 au from the Sun to the Sun, observed significantly higher intensities than Solar Orbiter, although their total separation angles are comparable.STEREO A observed a weak proton event but no significant increase of MeV electron intensities.At Earth/L1 (SOHO), the location with the poorest nominal magnetic connection with the flare site, the electron event seems more intense and distinct when compared with STEREO A. While Parker Solar Probe and Solar Orbiter were situated close to the ecliptic plane at slightly northern latitudes, BepiColombo, STEREO A and Earth were situated at southern latitudes (see Table 1, which summarizes the observer locations and their magnetic connections to the Sun) with a maximum of −7.2 • in the case of BepiColombo.The 17 April 2021 event, therefore, shows not only a spatial asymmetry with respect to the flare longitude but also clear differences between the electron and proton distributions.
We investigate here the drivers for this wide SEP spread as well as the reasons for the observed asymmetries.The most common explanations for widespread events have been a large acceleration region, e.g. an extended coronal shock (e.g., Rouillard et al. 2012;Gómez-Herrero et al. 2015;Lario et al. 2016;Rodríguez-García et al. 2021;Kouloumvakos et al. 2022), and/or efficient perpendicular transport in the corona or interplanetary medium (e.g., Dresing et al. 2012;Dröge et al. 2016).
The paper is organized as follows.After discussing the instrumentation used in this study in Sect.2, we begin with a detailed analysis of the magnetic connectivity between the different observer locations and the Sun (Sect.3).Section 4.1 discusses the complex and long-lasting flare of the event, Sect.4.2 describes the analysis of the associated coronal mass ejection (CME) and CME-driven shock.Sect.4.3 presents observation and analysis of the associated extreme ultra-violet (EUV) wave.Observations and a reconstruction of the coronal CME-driven shock are presented in Sect. 4.4. Sect. 4.5 provides an analysis of the various type II and III radio bursts observed during the event.In Sect.5, we study the interplanetary context, involving also simulations of the state of the interplanetary medium using 3D magnetohydrodynamic (MHD) simulations, and present an overview of the multi-spacecraft SEP observations.A more detailed analysis of SEP onset times, velocity dispersion, and pitchangle distributions is presented in Sect.5.3.In Sect.6, we relate the timing of SEP arrivals with solar counterpart observations to infer the parent source regions of the SEPs.In Sect.7, interplanetary transport modeling results are presented assuming two different scenarios: the first being the standard scenario of a single SEP injection into interplanetary space from a single source region, while a second scenario assumes multiple SEP injections from different particle sources at different times.Sections 8 and 9 provide the discussion and conclusions of the study presented here, respectively.

BepiColombo
Several data sets from the cruise phase of BepiColombo (Benkhoff et al. 2021) en route to Mercury are used in this study, such as from the Solar Intensity X-Ray and Particle Spectrometer (SIXS; Huovelin et al. 2020) on board the Mercury Planetary Orbiter (MPO, the European spacecraft involved in the BepiColombo mission).SIXS provides measurements of high-energy electrons and protons with the SIXS-P particle detector.This instrument consists of five orthogonal 150 µm thick Si PIN detectors, also called 'Sides', and a 5 × 5 × 6.3 mm 3 CsI(Tl) scintillator with photodiode read-out.It detects electrons in the range 50 keV to 3 MeV and protons in the range 1 to 30 MeV with a total nominal geometric factor of about 0.19 cm 2 sr.We note that Sides 0 and 4 are partially and totally obstructed by the spacecraft cruise shield, respectively.
We also use data from the BepiColombo Environment Radiation Monitor (BERM; Pinto et al. 2022) on board MPO, which is a particle detector that consists of a single silicon stack telescope with a small particle entrance of 0.5 mm 2 and a 50 µm beryllium cutoff window.Particles are arranged into 5 channels for electrons (∼0.15-10MeV), 8 channels for protons (1.5-100 MeV), and 5 channels for heavy ions (1-50 MeV•mg -1 •cm -2 ).BERM is mounted behind the radiator panel and faces the anti-sunward direction.
In addition, data from the Solar Particle Monitor (SPM) on board the Japanese s/c Mercury Magnetospheric Orbiter (MMO, also known as Mio; Murakami et al. 2020) are employed.SPM is a particle detector that forms part of the housekeeping suite, and it consists of two silicon photodiodes (SPM1 and SPM2), each one with an effective area of 10 mm × 10 mm and a depletion layer thickness of 0.3 mm.Each sensor has four different deposited energy channels, which cover the energy ranges 70-1170 keV and 50-200 MeV respectively.A calibration of the sensors is currently being performed with Monte Carlo simulations based on Geant4 (Agostinelli et al. 2003).
Finally, we also use data from the BepiColombo MPO magnetometer (MPO-MAG; Heyner et al. 2021), which is composed of two tri-axial fluxgate magnetometers placed on a 2.9 m boom.MPO-MAG measures magnetic field up to 128 Hz in a ±2048 nT range.

Parker Solar Probe
Energetic particles at Parker Solar Probe (PSP; Fox et al. 2016) are measured by the Integrated Science Investigation of the Sun (IS IS; McComas et al. 2016) suite.Low-energy electrons and ions (∼20 keV to 20 MeV/nucleon over 2π stereoradians) are covered by the Energetic  Wiedenbeck et al. 2017), consisting of three telescopes of stacked solid-state detectors, using the standard dE/dx versus residual energy technique to measure ions from ∼1 to >100 MeV/nuc and electrons in the range ∼ 0.5-6 MeV.The first two Low Energy Telescopes (LETs) of EPI-Hi consist of a double-ended detector, providing oppositely viewing apertures (LETA and LETB) and one single-ended detector (LETC) with a viewing axis perpendicular to that of LETA.The third telescope (High Energy Telescope; HET) covers the highest energies and is doubleended with two apertures (HETA and HETB) providing roughly sunward and anti-sunward viewing directions along the nominal Parker spiral.
Observations of the magnetic field are obtained from the fluxgate magnetometer part of the FIELDS (Bale et al. 2016) suite, and solar wind measurements are provided by the Solar Probe Cup (SPC) instrument part of the Solar Wind Electrons Alphas and Protons (SWEAP; Kasper et al. 2016) investigation.
Radio observations are provided by the Radio Frequency Spectrometer (RFS; Pulupa et al. 2017) part of FIELDS, which is a dual-channel digital spectrometer designed for both remote-sensing of radio waves and in-situ measurements of electrostatic fluctuations between 10 kHz and 19.2 MHz.Here, we use the RFS data when input channels were set to the two pairs of crossed dipoles.Besides the radio flux density, it also allows us to retrieve the degree of circular polarization (Pulupa et al. 2020).

Solar Orbiter
Data from several instruments on board Solar Orbiter (Müller et al. 2020) are used.The Spectrometer/Telescope for Imaging X-rays (STIX; Krucker et al. 2020) provides imaging spectroscopy in the X-ray range (4-150 keV).It has a full-disk field of view (FOV) and sub-second time resolution.The Radio and Plasma Waves (RPW; Maksimovic et al. 2020Maksimovic et al. , 2021;;Vecchio et al. 2021) instrument on Solar Orbiter consists of several subsystems including the Thermal Noise and High Frequency Receiver (TNR-HFR or THR) with a dual channel sweeping receiver in the range from 4 kHz up to 16 MHz.In particular, THR provides measurements of the plasma quasi-thermal noise (QTN) in the range 4 kHz -1 MHz.When the QTN signal is quite strong, the spectral peak at the electron plasma frequency can be identified from which the in-situ absolute electron density can be derived (Meyer-Vernet et al. 2017;Khotyaintsev et al. 2021).The properties of energetic particles as measured by Solar Orbiter are studied using the Electron Proton Telescope (EPT) and the High Energy Telescope (HET) of the Energetic Particle Detector (EPD; Rodríguez-Pacheco et al. 2020) instrument suite.Both sensors consist of two double-ended telescopes.EPT and measures ions and electrons in the energy ranges 20 keV -15 MeV and 20-400 keV, respectively, and HET relativistic electrons between 0.3 and 30 MeV and protons between 7 and 107 MeV.
The Solar Orbiter Magnetometer (MAG; Horbury et al. 2020) is a fluxgate vector magnetometer, yielding in-situ measurements of the interplanetary magnetic field with 16 vectors/s (normal mode) and up to 128 vectors/sec (burst mode).
The lower-energy, thermal, and suprathermal particles are measured by the Solar Wind Analyzer (SWA; Owen et al. 2020) suite.In this work, measurements from the SWA Proton and Alphas Sensor (PAS), sampling 3D velocity distribution functions of protons and alpha particles in the 0.2-20 keV energy range with a 4 s time cadence, are used to address the in-situ plasma moments, such as the solar wind's bulk flow speed and density.

STEREO A
Observations from several instruments on board the Solar Terrestrial Relations Observatory (STEREO; Kaiser et al. 2008) are used in this study.As the STEREO B spacecraft is inactive since October 2014 due to multiple hardware anomalies, only data from instruments on board STEREO A are available for the period under consideration.
The STEREO/WAVES (S/WAVES; Bougeret et al. 2008) instrument provides comprehensive measurements of all components of the electric field fluctuations between 2.5 kHz and 16 MHz.It allows us to locate sources, and calculate the polarization state (including apparent source sizes) of radio emissions in a heliocentric distance range from 4 R to 1 au, while the flux density can be measured even down to 2 R (Krupar et al. 2014).Unfortunately, direction-finding data were not available for this event.
Interplanetary magnetic field measurements are provided by the Magnetic Field Experiment (MFE; Acuña et al. 2008), part of the In situ Measurements of Particles And CME Transients (IMPACT; Luhmann et al. 2008) instrument suite.MFE is a triaxial fluxgate magnetometer mounted on a telescopic boom at a distance of ∼3 m from the spacecraft body, reaching a maximum cadence of 32 vectors/s.
Energetic particle observations with 1-minute cadence are provided by several instruments, part of the IMPACT investigation.The Solar Electron Proton Telescope (SEPT; Müller-Mellin et al. 2008) consists of dual double-ended magnet/foil particle telescopes measuring 30-400 keV electrons and 60-7000 keV ions.Two separate units provide anisotropy information in four different looking directions: Sun, Asun (pointing sunward and anti-sunward along the nominal Parker spiral, respectively), North, and South (pointing towards the North and South ecliptic poles, respectively).Since July 2015, after the solar conjunction, the spacecraft was rolled 180 • about the spacecraft-Sun line and these nominal pointing directions changed, with the Sun and Asun telescopes pointing perpendicular to the nominal Parker spiral direction, North pointing southward and South pointing northward.The Low Energy Telescope (LET; Mewaldt et al. 2008) measures protons from ∼2 to ∼13 MeV and heavier ions from ∼2 to >40 MeV/nuc (species-dependent energy range).The field of view is divided into 16 different sectors, providing directional information.The High Energy Telescope (HET;von Rosenvinge et al. 2008) provides the highest energy measurements, including 0.7-4 MeV electrons and 13-100 MeV protons.
Physical properties of the solar wind plasma are obtained by the Plasma and Suprathermal Ion Composition (PLASTIC; Galvin et al. 2008) instrument, in particular by the Solar Wind Sector (SWS), sampling solar wind proton bulk parameters.
Remote-sensing observations from STEREO A are provided by several instruments that are part of the Sun Earth Connection Coronal and Heliospheric Investigation (SECCHI; Howard et al. 2008).This instrument suite includes an Extreme UltraViolet Imager (EUVI; Wuelser et al. 2004), two coronagraphs (COR1 and COR2) imaging the corona from 1.4 up to 15 R , and two Heliospheric Imager (HI; Eyles et al. 2009) cameras (HI1 and HI2).

Near-Earth spacecraft
From the Wind spacecraft (Ogilvie & Desch 1997), we use the Magnetic Field Investigation (MFI; Lepping et al. 1995) instrument, measuring at a cadence of 11 vectors/s.The Solar Wind Experiment (SWE; Ogilvie et al. 1995) and Three-Dimensional Plasma and Energetic Particle Investigation (3DP; Lin et al. 1995) instruments provide energetic particle measurements from which we use electron observations in the range of ∼40-600 keV.Finally, the Wind/WAVES (WAVES; Bougeret et al. 1995) instrument measures the electric field from 0.3 Hz up to 13 MHz using three dipolar antennas.
From the Solar and Heliospheric Observatory (SOHO; Domingo et al. 1995), we use energetic proton measurements of the Energetic and Relativistic Nuclei and Electron (ERNE; Torsti et al. 1995) covering energies of a few to a hundred MeV, and energetic electron measurements in the MeV range by the Electron Proton Helium Instrument (EPHIN), part of the Comprehensive Suprathermal and Energetic Particle Analyser (COSTEP; Müller-Mellin et al. 1995) suite, and coronagraph observations by the Large Angle and Spectrometric Coronagraph (LASCO; Brueckner et al. 1995).
Extreme ultraviolet (EUV) solar images were obtained by the Atmospheric Imaging Assembly (AIA; Lemen et al. 2012) on board the Solar Dynamics Observatory (SDO; Pesnell et al. 2012).

Mars
Observations at Mars have been obtained from two spacecraft in orbit around the planet.Energetic particle data come from the Solar Energetic Particle (SEP; Larson et al. 2015) instrument on board the Mars Atmosphere and Volatile Evolution (MAVEN; Jakosky et al. 2015) mission.MAVEN/SEP is a solid-state telescopic detector with two identical sensors (SEP1 and SEP2), each containing two oppositely arranged double-ended telescopes (A and B) and measuring ions in the energy range ∼20-6000 keV and electrons in the range of ∼20-200 keV.SEP1 and SEP2 are body-mounted onto the MAVEN spacecraft to provide orthogonal look directions, with each telescope providing a 42 • × 31 • FOV coverage.
Solar wind density and velocity come from the Analyzer of Space Plasmas and Energetic Atoms (ASPERA-3; Barabash et al. 2006) on board Mars Express (MEX; Chicarro 2004), and in particular from the Ion Mass Analyser (IMA) sensor that measures ions in the energy range 10 eV/q -30 keV/q.

Spacecraft constellation and magnetic connectivity
For all locations at which this event was observed, with the exception of Mars, we derive the instantaneous magnetic connectivity to the solar surface at 16:00 UT on 17 April 2021 using two different coronal models.The first model is a standard potential-field source surface (PFSS; Schatten et al. 1969;Altschuler & Newkirk 1969;Wang & Sheeley 1992) model out to 2.5 R , and the second is the Wang-Sheeley-Arge (WSA; Arge & Pizzo 2000;Arge et al. 2003Arge et al. , 2004;;McGregor et al. 2008) model.The latter makes use of the Schatten Current Sheet (SCS) model to extend the PFSS solution in this work to 5 R , providing a more realistic magnetic field topology of the upper corona.This height is appropriate for applications of deriving spacecraft connectivity since WSA is magnetostatic and is designed to be most accurate in low beta regimes, whereas when coupled with an MHD model as its inner boundary condition one would want to derive the coronal field to between 20-30 R to ensure that the solar wind is supersonic and super alfvénic (Arge et al. 2004).In both the WSA and traditional PFSS approach, the observed solar wind speed at the time of the event is used to backmap the event to the model-derived coronal field, assuming a Parker spiral.
In the case of BepiColombo, which did not measure the solar wind speed during cruise phase, a nominal value of 400 km s −1 is used, which coincides with the simulated solar wind given by the ENLIL model (Odstrcil et al. 2004), as discussed in Sect. 5.In this section, we present the results from both models, discuss any differences between the solutions, and elaborate on any uncertainties associated with magnetic connectivity.
Both the PFSS and WSA coronal solutions are obtained using an identical Air Force Data Assimilative Photospheric Flux Transport (ADAPT; Arge et al. 2010Arge et al. , 2011Arge et al. , 2013;;Hickmann et al. 2015) time-dependent photospheric field map, derived using data from Global Oscillation Network Group (GONG; Harvey et al. 1996) magnetograms.ADAPT uses flux-transport modeling (Worden & Harvey 2000) to account for solar time-dependent phenomena (e.g., differential rotation, meridional, and supergranulation flows) when observational data are not available.This is especially useful for studying events that originate on the solar far-side (i.e., the solar hemisphere not visible from Earth).Since ADAPT is an ensemble model, it provides 12 possible states (i.e., realizations) of the solar surface magnetic field, ideally representing the best estimate of the range of possible global photospheric flux distribution solutions at any given moment in time.When coupled with the WSA model, ADAPT-WSA derives an ensemble of 12 realizations representing the global coronal field and spacecraft magnetic connectivity to 1 R (i.e., the solar surface) for a given moment in time, providing an estimate of the uncertainty in the ADAPT-WSA solution.The best realization is then determined by comparing the model-derived and the in-situ observed radial interplanetary magnetic field (IMF) and solar wind speed.
Figure 2 shows the instantaneous connectivity derived with the PFSS coronal field solution, with the corresponding footpoint connectivity listed in Table 1 columns (6)-( 7).The plot shows the Sun in the center, the source surface (dashed circle), which is the outer boundary of potentialfield models, and the spacecraft constellation in the heliospheric Carrington coordinate system, where the unit of distance is the solar radius.The magnetic connectivity to different spacecraft is estimated as a nominal Parker spiral connecting to the source surface, from which magnetic field lines are tracked downwards to the photosphere using a PFSS extrapolation.Table 1 shows the magnetic connection points from the various spacecraft to the photosphere and the observed solar wind speed that is used to calculate the Parker spiral.Note that the scale in the plot is logarithmic above the source surface and linear below.
Figure 3 shows the WSA-derived instantaneous magnetic connectivity and thus the model-estimated magnetic footpoint for each of the five spacecraft on 17 April 2021 at 16:00 UT for all 12 realizations of the ADAPT-WSA output.The 12 realizations often produce similar results causing overlapping footpoints, which are overlaid onto the WSA-derived coronal holes shaded in red (positive/outwardly-directed field) and blue (negative/inwardly-directed field).Since this plot is a summary of the connectivity for all 12 realizations, the coronal hole shading represents any grid cell derived as open by any of the 12 realizations.The heliospheric current sheet (HCS) Fig. 2: Semi-logarithmic representation of the spacecraft constellation in Carrington coordinate system.The plot is linear inside the dashed black circle, which marks the distance of the potential field source surface at 2.5 R in this case), and the orange circle marks the Sun.Above 2.5 R , the plot is logarithmic in distance.Color-coded solid circles mark the various spacecraft of the constellation, and the lines connected to them represent the nominal Parker spiral solutions computed considering their heliocentric distances and the observed solar wind speeds.Inside the dashed black line, the magnetic connection is extrapolated with a PFSS solution, where the color of the lines corresponds to heliospheric latitude.The black arrow corresponds to the flare location.
is overplotted in yellow and is nearly parallel with the solar equatorial plane.The average and standard deviation over all realizations of footpoint connectivity are calculated and shown in Table 1 columns ( 9)-( 10), with the exception of STEREO A. For this spacecraft, the values in Table 1 represent eight of the 12 ADAPT-WSA realizations which derived the magnetic footpoint at the southern polar coronal hole boundary.The other four realizations derive the source region of this event at the northern polar coronal hole boundary, with an average and standard deviation of 308.8 • ±1.1 • Carrington longitude, 63.1 • ±0.3 • heliographic latitude.This is common when the spacecraft is near the HCS (discussed in more detail below).It is important to note that in several instances on this table, the standard deviation of the footpoint connectivity falls within the 2.0 • resolution limit of the WSA model.The standard deviation is only included to show the precision and range of variance among the 12 realizations.
For all spacecraft except STEREO A, the 12 ADAPT-WSA realizations produce very similar results for the model-determined magnetic footpoints.The largest standard deviation that was calculated was 8.4 • for Parker Solar Probe's longitudinal footpoint connectivity (purple).This is likely because Parker Solar Probe observed this event on the solar far-side as seen from Earth, where there are no observations of the photospheric field to update our solution.All other spacecraft have nominal standard deviations in footpoint latitude and Carrington longitude, giving us more confidence in our results.
When comparing the results from the PFSS model vs. WSA, both the model-derived polarities (Table 1 column  AR 12818 Fig. 3: ADAPT-WSA derived instantaneous magnetic connectivity on 17 April 2021 at 16:00:00 UT for five of the spacecraft observing the SEP event, overlain onto the corresponding ADAPT-GONG map used to derive the coronal field, and the WSA-derived coronal holes (red/positive, blue/negative).The footpoint connectivity for each spacecraft is labeled by 12 colored points, one for each ADAPT-WSA realization, and the WSA-derived HCS is overplotted in yellow.The locations of two ARs that emerge on the solar far-side are labeled with an "×", yet are not incorporated into the ADAPT map until several days after this event.AR 12818 associated with a solar flare (discussed in Sect.4) is labeled with an orange "×".Orange shading marks the portion of the Sun not observed by remote imagers on board STEREO A, Solar Orbiter, or spacecraft near Earth.
(8)) and the footpoint connectivity agree overall, with the exception of the magnetic footpoint derived for Parker Solar Probe.The PFSS model derives the source region of the Parker Solar Probe-observed event on the boundary of the northern polar coronal hole extension (at 28.6 • latitude, 146.7 • Carrington longitude), whereas the ADAPT-WSA derived source region is at the boundary of the northern polar coronal hole (at 65.3 • latitude, 131.1 • Carrington longitude).Differences in the two model solutions could arise for a few reasons, a primary one being that Parker Solar Probe observed this event on the solar far-side where we do not have recent observations of the photospheric field to drive coronal models.Nevertheless, both models derived the footpoint locations of this event for all five spacecraft as originating from the boundaries of coronal holes, with each spacecraft situated within 5 • of the HCS.When the solar wind originates from the HCS at locations where it is nearly parallel to the solar equatorial plane, there is increased uncertainty in the backmapped locations of the magnetic footpoints for observers at 1 R when using any coronal model.This is because a difference of a few degrees from the HCS (i.e., 1 -2 model grid cells) could connect the spacecraft to either side of the streamer belt.It is also common in this scenario for the spacecraft-measured polarity to fluctuate between inward and outward connectivity as the spacecraft never becomes sufficiently separated from the HCS.How-ever, for this event, both models accurately derive the solar wind magnetic field polarity, which is measured by each of the five spacecraft in situ (Table 1 column (8)), giving us more confidence in our results.
Selecting the best ADAPT input map to drive both models is particularly challenging for this event because the spacecraft were widely separated in longitude, whereas this type of modeling produces the most accurate results for spacecraft connected to the most recently added photospheric field observations (i.e., in this case Earth and STEREO A).Additionally, Parker Solar Probe and Bepi-Colombo observed this event on the solar far-side.Solar Orbiter was also located on the far-side; however, the spacecraft was connected to the near-side (i.e., the solar hemisphere visible from Earth) at 1 R .Lastly, there were two far-side active regions (ARs) that rotated onto the nearside on 19 and 22 April (labeled with an "×" in Fig. 3).Although they are not visible in the ADAPT map from 17 April, the locations of these ARs are labeled with crosses in Fig. 3. New far-side AR emergence is problematic for all coronal models (Wallace et al. 2022).To account for this evolution, we test various input maps from 17 -23 April into both the PFSS model and WSA.We find that the connectivity for each spacecraft does not change drastically by using any particular map in this date range.Therefore, we select a map from the time closest to the SEP event, 17 April 2021 at 16:00:00 UT.It is important to note that the two far-side ARs fall inside the visible hemispheres observed by Solar Orbiter and STEREO A during the time of the SEP event, making it possible to identify if any of these ARs are associated with a solar flare.One of the farside ARs located at −19.09 • latitude, 204.73 • Carrington longitude was associated with a solar flare (discussed in detail in Sect.4).Additionally, we can be confident that none of these far-side ARs produced flares within the longitudinal sector in which no remote observations of the solar corona are available (i.e., from ∼45 -125 • Carrington longitude, shaded in orange in Fig. 3) comprising the location of Parker Solar Probe, situated at 104 • Carrington longitude.

Observations of the associated flare
The SEP event was associated with a solar flare occurring in the active region that was assigned the NOAA AR number 12818 when it rotated onto the Earth-facing hemisphere three days later.While the flare was clearly visible in the field of view of the STEREO A EUVI instrument, it was initially occulted as seen from from Earth, but later phases of the eruption could be seen above the limb.Starting around 15:45 UT, the eruption of a flux rope was observed in EUV with SDO/AIA above the northeastern limb.Note that because we use observational assets at different locations, we have shifted all times pertaining to flare observations to UT at the Sun.From 16:03 UT on, the cusp-shaped top of a flaring arcade became visible in the 131 Å channel.At wavelengths corresponding to lower temperatures, flaring loops appeared only after 18:15 UT, consistent with a considerable occultation angle.The GOES soft X-ray flux started to increase at 16:15 UT and peaked at 17:10 UT as a GOES class of B9.7 (see top panel in Fig. 4).
In contrast to Earth-based assets, the whole flare was visible from Solar Orbiter, and observed in hard X-rays (HXR) with STIX.We show STIX count rates integrated over two energy bins in the bottom panel of Fig. 4. The counts are background-subtracted and normalized to the peak count rate in the two ranges.The thermal HXR emission at 4-10 keV (generated by the hot plasma) increased from 15:55 UT onward, peaked at 16:22 UT, and decayed to pre-event background levels only at around 19:00 UT, thus indicating a long-duration event.Again, these times refer to when events have happened on the Sun.At Solar Orbiter, they were observed 7 min later.Based on a statistical comparison of STIX and GOES/SXR X-ray fluxes for flares that were fully visible for both instruments, the true GOES class can be estimated to ∼C5 1 This is also shown by the fact that the GOES fluxes peak more than half an hour after the STIX thermal count rate, since GOES sees significant emission only when larger loops due to magnetic reconnection become filled by hot plasma later in the event.
Above 15 keV, the STIX light curves show the spiky behavior typical for the non-thermal HXR emission generated 1 See the STIX website for a description of the method: https://datacenter.stix.i4ds.net/wiki/index.php?title= GOES_Flux_vs_STIX_counts.The discrepancy between the B9.7 obtained from actual GOES observations and the class estimate from STIX is mainly due to occultation of the majority of hot flare plasma as seen from Earth.  .Bottom: Normalized background-subtracted STIX count rates integrated over two different energy bands.Note the gradual evolution of the thermal emission at 4-10 keV (blue) as opposed to the multiple non-thermal spikes seen at 15-25 keV (green) and at 25-50 keV (red; multiplied by 0.3 for clarity).For both GOES and STIX, times have been shifted so that they refer to UT at the Sun.by accelerated electrons that precipitate into the chromosphere.At least 13 non-thermal spikes are identified.While the non-thermal emission phase is usually restricted to a few minutes in typical C-class flares (cf.Veronig et al. 2002), this flare shows non-thermal emission over 50 min.In particular, there are two major peaks in the late phase that in contrast to the other show emission above 25 keV, indicative of a comparatively harder spectrum.In Appendix A, we provide a full spectral analysis of the event using the STIX data.
HXR images can be reconstructed from pixelated STIX science data.Figure 5 shows the HXR sources for the eight HXR peaks that had the largest number of non-thermal counts overplotted on STEREO A EUVI 304 Å images that have been rotated so that they correspond to the viewpoint of Solar Orbiter.In the EUVI frames, we mainly see the chromospheric flare ribbons, thus such a reprojection that assumes that all features are lying on the solar surface does not introduce significant projection artefacts.Red contours show the coronal thermal source (6-10 keV), and the blue contours show the chromospheric non-thermal footpoints (15-25 keV).All images are reconstructed with the Expectation Maximization algorithm (Massa et al. 2019).Normally, the precise source locations are provided by the STIX Aspect System (Warmuth et al. 2020).However, Solar Orbiter was at a heliocentric distance of 0.84 au, which is too far from the Sun to provide a reliable pointing solution.We therefore apply the average image displacement obtained from other events from the cruise phase where aspect information was available as implemented in the STIX imaging software.This method yields a mean position uncertainty of about ±10 .The flare position (plotted in heliocentric Cartesian coordinates in Fig. 5) is at the heliographic co-Table 1: Magnetic connectivity between spacecraft and the Sun.Columns (1)-( 4) present the respective observer and its location in Carrington coordinates (with the first row providing the flare location).Column (5) lists the measured solar wind speed, ( 6)-( 7) and ( 9)-(10) provide the backmapped magnetic footpoints of the observer at the solar surface using the simple PFSS and the ADAPT-WSA models, respectively.Column (8) presents the magnetic field polarity observed (O) in situ and modeled (M) by ADAPT-WSA and PFSS. (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) PFSS (a)  ADAPT-WSA    The coronal thermal source undergoes little evolution little evolution throughout the flare.One might expect to observe a pair of non-thermal sources consistent with the footpoints of the magnetic loops containing the hot plasma (cf.Fletcher et al. 2011).However, most HXR peaks only show a single footpoint at the eastern edge of the thermal source.The issue here is that while the individual nonthermal peaks are all very clearly defined, the total number of counts above 15 keV per peak is quite low, on the order of 1 000-2 000 counts.This is marginal for imaging, particularly in case there is more than one source present.Nevertheless, we find that all non-thermal peaks are originating from the same active region, and there is no evidence of a second remote source.While the presence of such a secondary source cannot be ruled out, it would have to be weaker than the main source by a factor of 5-10.We conclude that the footpoint brightness was very asymmetric in this event, with the eastern footpoint clearly dominat-  ing.This is consistent with the flare ribbons seen at 304 Å, where also the southwestern ribbon is the dominating one.
The different non-thermal peaks are not associated with changing footpoint locations.

CME observations
A CME erupted from the same active region as the associated flare, NOAA AR 12818, located in the southern solar hemisphere at heliographic coordinates E111S18 (203 • Carrington longitude) on the day of the event.The active region entered Earth's field of view on 20 April.The evolution of the CME was observed ∼30 • from the eastern limb towards the central meridian as seen from STEREO A. The bottom panels of Fig. 6 show EUV images taken by SDO/AIA (left) and STEREO A/EUVI (right) at ∼16:10 UT (all times refer to observation times at the spacecraft).At this time, we observe the first clear indication of the eruption, when the CME exhibits prominent signatures of overexpansion (e.g.Patsourakos & Vourlidas 2009), evidenced by the bubble-like appearance in EUV images.At the same time, the flare ribbons are activated along an arched path.
Between 16:12 UT and 16:25 UT, the CME continued to expand as it reached the edge of the field-of-view of the EUV instruments.The left panel of Fig. 7 shows the time (∼16:50 UT) at which the CME was clearly visible by the coronagraph imagery: COR1 (second column) and LAS-CO/C2 data (first column).As shown in the images, the CME morphology in white light images at these heights is consistent with classic flux-rope characteristics, namely featuring the presence of a bright outer rim followed by a cavity (e.g.Vourlidas et al. 2013).
The angular separation between STEREO A and Earth was ∼53 • , which still enables a reliable 3D reconstruction of the CME (e.g.Balmaceda et al. 2018;Verbeke et al. 2022).For this purpose, we used the graduated cylindrical shell (GCS; Thernisien et al. 2006) model to reproduce the CME appearance by fitting pairs of EUV (at distances below 1.5 R ) and white-light (from ∼2.5 to 22 R ) images.The model consists of a croissant-like structure fully described by six free parameters: three for location and orientation (latitude and longitude of the CME leading-edge, and tilt or inclination of the main axis of the CME with respect to the solar equator), and three for the geometry (height; aspect ratio, which sets the rate of expansion versus the height of the CME; and angular separation of the legs or half-angle).The sensitivity (deviations) in the parameters of the GCS analysis is given in Table 2 of Thernisien et al. (2009).It is worth noting that these parameters are sensitive to image quality and human interpretation (Verbeke et al. 2022).The routine used for the reconstruction is rtcloudwidget.pro,available as part of the scraytrace package in the SolarSoft IDL library2 .
The bottom panels of Fig. 6 and Fig. 7 show the GCS fit analysis, where the green mesh represents the flux rope structure.The 3D reconstruction shows that the CME follows a non-radial path towards the solar equator in the early evolution with the latitude varying from −14 • to −9 • from 16:12 to 17:23 UT.The longitude and the tilt angle, meanwhile, do not show deviations, staying at fixed values of −116 • and −70 • , respectively.The GCS parameters were chosen to better describe the portion of the CME oriented towards Solar Orbiter, as the croissant-like shape used for the fitting could not represent fully the CME due to its non-radial propagation and curved axis.The last term was introduced by Rodríguez-García et al. ( 2022) to refer to flux ropes that may deviate from the nominal semi-circular (croissant-like) shape and have instead an undulating axis.
The CME speed at the leading-edge estimated from the linear fit to the height-time measurements is 880 km s −1 .The width of the CME is estimated based on Dumbović et al. (2019), where the semi-angular extent in the equatorial plane is expressed by Then, the total angular extent of the CME is 46 • .The value of R maj (face-on CME half-width) is calculated by adding R min (edge-on CME half-width) to the half-angle, and R min was calculated as the arcsin(aspect ratio).The CME width deviation was derived from the mean half-angle error, estimated by Thernisien et al. (2009) as +13 • /−7 • .Thus, at the latest time of the 3D reconstruction at 19:23 UT, corresponding to a CME height of 15.5 R , the narrow CME (∼46 • ) is propagating in the direction E116S09 with a moderate speed (∼880 km s -1 ).Notes.a Height from the solar center of the intersection point between the reconstructed CME-driven shock and the field line connecting to the respective spacecraft.

EUV wave observations
Figure 8 together with the accompanying movie shows an overview of the EUV wave evolution in STEREO A/EUVI 195 Å running-difference images created with a lag of 150 s.
The prominent signatures of the EUV wave, which exhibits a quasi-circular propagation away from the eruptive center over the solar disk, are already clearly seen around 16:10 UT and can be followed for about 40 min in STEREO A. As follows from the derivation of EUV wave kinematics and perturbation characteristics in Appendix C, the EUV wave on the solar disk extends to about 680 Mm from the source region with a mean velocity of 223-327 km s −1 .Above the solar limb, in the northern direction, the EUV wave can be followed to a distance of about 740 Mm, propagating with speeds of 260-450 km s −1 (for heights increasing from 1.05 to 1.15 R ).At the same time, in the southern direction, the wave is seen only to about 350 Mm, propagating with speeds of 220-300 km s −1 .As seen from the movie, accompanying Fig. 8, the EUV wave reaches the backmapped magnetic footpoints of BepiColombo (yellow) at around 16:55 UT (point 3) and at around 17:00 UT (points 1 and 2).Each of the points 1, 2 and 3 correspond to the spacecraft's magnetic field footpoints obtained using either ADAPT-WSA (point 1), the PFSS model at 1 R (point 2) and PFSS at a height of 100 Mm above the photosphere (point 3).The footpoints of other spacecraft, which lie on the visible hemisphere as seen by STEREO A, are displayed in other colors as described in the figure legend.

CME-driven shock observations
In white-light images the signatures of a shock wave formed in front of the expanding flux rope are faint.We use calibrated, excess-mass images (i.e.pre-event image subtracted) and display them in Fig. 7.By 16:30 UT, when the CME front is visible in COR1 FOV, the EUV wave is still visible on the surface.The CME exhibits a diffuse front ahead the brighter rim, more clearly seen at the north flank in both COR1 and LASCO-C2 images (marked with white arrows in the left panels of Fig. 7).This typical "two-front" morphology is generally interpreted as evidence of a CMEdriven shock in white-light images (Ontiveros & Vourlidas The magnetic footpoints are determined using a combination of ballistic backmapping in the heliosphere and backmapping below the source surface using ADAPT-WSA to 1 Rs (points 1), a standard PFSS model to 1 R (points 2), and to a height of 100 Mm above the photosphere (points 3).As seen in the accompanying video, the EUV wave reaches the BepiColombo footpoints at around 16:55 UT (point 3) and at around 17:00 UT (points 1 and 2).Times refer to the observation time at STEREO A.
2009; Vourlidas et al. 2013).The EUV wave is visible on the disk until 16:50 UT.By this time, the CME reaches the edge of COR1 FOV.At larger distances, namely COR2 FOV, a diffuse arched-shaped feature (white arrow in right panel of Fig. 7) is also seen propagating on the northwest quadrant.This feature is best visible between 18:23 UT and 20:23 UT in COR2 FOV and may result from the compression of a relatively weak shock wave against the underlying coronal structures.We use these features to estimate the angular extension of the shock.For this, a spherical surface (Kwon & Vourlidas 2017) is used to model the 3D appearance of the shock (represented by the red mesh in Fig. 7).
From the 3D reconstruction, we estimate that the shock reaches a speed of ∼1500 km s −1 below 5 R and is propagating on the direction between Solar Orbiter and Bepi-Colombo, consistent with the direction estimated from the CME 3D modeling in Sect.4.2.Following Kwon & Vourlidas (2017), shown in their figure 2, we determine the angular width of the shock to be ∼180 • at 19:23 UT, corresponding to a height of the shock nose of ∼16.3 R .Table 3 shows the timing of the first intersection between the coronal shock reconstruction and the magnetic field lines obtained with the ADAPT-WSA model connecting to the different spacecraft analyzed in this study.All times refer to observation times at 1 au.

Radio observations
The radio emission associated with the eruption on 17 April 2021 was observed by ground-based and space-borne instruments and includes both type II and type III radio bursts.Type II (TII hereon) bursts are related to acceleration of energetic electrons at shock waves (Krasnoselskikh et al. 1985;Benz & Thejappa 1988, and references therein), while type III (TIII hereon) radio bursts are signatures of fast electron beams propagating via open (or quasi-open) magnetic field lines from the corona to interplanetary space (Zheleznyakov 1965;Jebaraj et al. 2023b).In Fig. 9, we present a dynamic radio spectrum using measurements from the ground-based e-Callisto instrument located at Swiss Landschlacht, providing observations in the 80 MHz to 10 MHz range.This spectrum shows a poorly-observed decametric TII radio burst starting at 16:26 UT and exhibiting both fundamental and second harmonic emission lanes.Harmonic emission is brighter compared to the fundamental due to the large angle between the source and Earth (directivity of harmonic emission has a wider angle than the fundamental; Zheleznyakov & Zaitsev 1970;Tkachenko et al. 2021).
Figure 10 presents a combined dynamic radio spectrum of all available hecto-kilometer observations from all the observing spacecraft, namely, Parker Solar Probe, Solar Orbiter, STEREO A, and Wind.The spectrum shows a number of different radio emissions including groups of TIII bursts and distinctly patchy TII emission.A list of the starting times of each radio burst as observed by different spacecraft is provided in Table 4.An interesting aspect is that most of the TIII radio burst groups were best or exclusively observed by PSP/FIELDS/RFS partly due to its enhanced resolution and sensitivity (Pulupa et al. 2017).But the radial distance of the spacecraft from the Sun and also the directivity of the emission at the source (Thejappa et al. 2012) also play a key role.Jebaraj et al. (2020) have suggested that the intensity of a radio burst is higher in the direction of the source propagation.Therefore, the intensity of the radio emission at different observers depends on both the position of the observing spacecraft and the intrinsic directivity of the radio source.This explains the emission intensity at Parker Solar Probe, which was the closest spacecraft to the Sun during the flare-CME event.As we show in Fig. 1 (left), the spacecraft were located at different longitudes and radial distances (Table 1).In the following, we use the Parker Solar Probe spectra (Fig. 10  The different groups of TIII radio bursts and multiple components of the TII bursts exhibit interesting characteristics as far as their spectral morphology is concerned.The first and third TIII groups (TIII(1) and TIII(3) hereon) are rather faint at the short-hectometer wavelengths and appear to be intense across all spacecraft observations.The second and fourth type III groups (TIII(2) and TIII(4) hereon) were observed almost only by Parker Solar Probe and consisted of a large number of individual TIII bursts that were better distinguishable in the short-hectometer wavelengths.This indicates that during the time when TIII(2) and TIII(4) were observed, there were multiple smaller episodes of electron acceleration and subsequent release into the open magnetic field lines in the direction of Parker Solar Probe.This is further corroborated by the polarization measurements made by PSP/FIELDS/RFS, which indicate that the energetic electron beams (the sources of type III bursts) were strongly directed towards Parker Solar Probe.Appendix B discusses the details of the polarization measurements extensively.The results indicate that TIII(2) and TIII(4) originated from a region of negative magnetic field polarity.The relatively high degree of polarization (Fig. 10 panel 2) of TIII(4) at its origin also hints at a region with high magnetic field strength (e.g., ramp of a quasiperpendicular shock wave).As for the magnetic connectivity, Parker Solar Probe also observed Langmuir waves (see; Ginzburg & Zhelezniakov 1958;Melrose 1985) on multiple occasions, close to the local plasma frequency.This indicates that the electron beams generating TIII(2) and TIII(4) were directly sampled by Parker Solar Probe.
The TII bursts associated with the event are distinctly patchy and complex in the hectometer wavelengths (see Fig. 10).It is likely that the different TII components are associated with the same shock wave but at different regions.All TII bursts also appear bursty in terms of intensity variations (marked in Fig. 10), suggesting an on-and-off emission process at the shock front (Mann & Classen 1995).On-and-off TII bursts are believed to be emitted from locations on the shock wave where the upstream plasma conditions induce rapid changes to its obliquity and other characteristics (e.g., Schmidt & Cairns 2014;Jebaraj et al. 2021;Kouloumvakos et al. 2021).Due to their patchy nature, it is somewhat difficult to distinguish between them, but we identify two main TII radio components (TII(1) and TII(2), which are marked in Fig. 10) observed in the short hectometer wavelengths (16-13 MHz) together with TIII(1) and TIII(3), and TIII(4), respectively.Overall, it seems these patchy TII bursts were observed from the start of the event and the beginning of TIII(1) and continued even after TIII(4).Furthermore, we note another interesting temporal and spectral phenomenon observed together with TIII(4), namely, the presence of TII herringbone-like features (TII(HB) hereon).The observation of such a feature may either indicate interaction of the shock wave and TIII(4), or that some of the electron beams generating TIII(4) may just be herringbones accelerated at the near-perpendicular shock front.Herringbone features with no clear backbone emission are often clear indicators of shock fronts with nearperpendicular geometry (θ Bn = 87 • −89.9 • ), which are able to accelerate electrons along either side of the magnetic field lines interacting with the said shock front (Mann & Klassen 2005).We also note that observations beyond TIII(4) of patchy TII bursts may be associated with both TII(1) and TII(2).
A similar mechanism may also contribute to TIII(2), which was also observed uniquely by Parker Solar Probe.The polarization analysis of TIII(2) presented in Appendix B suggests that if there would have been a herringbone-like feature at the origin of these type III bursts, it would be observed in the decamater wavelengths.However, due to the lack of meter-decameter observations, it is not possible to make such a conclusion.
Solar Orbiter was the second-closest radio observer radially and also the second-closest spacecraft to the flaring active region in terms of the magnetic connectivity (see Fig. 1 and Fig. 2).The spacecraft observed mainly TIII(1) and TIII(3) and also TIII(4) at lower frequencies.Due to the limited survey-mode observations during the initial phase of the mission (Maksimovic et al. 2020), the low resolution HFR observations from the Solar Orbiter/RPW instrument make it difficult to recognize the strongly patchy type II burst.However, the likely intensity variations (in the frequency range 16-5 MHz) of the type II bursts can be seen in Fig. 10 panel 3.
At the time of the event, STEREO A was located almost diametrically opposite (∼180 • ) from Parker Solar Probe.It observed well the TIII(1) and TIII(3) and also partially the TIII(2).TIII(4) was observed faintly at lower frequencies at this location.This indicates that for STEREO A the source region of the TIII(2) and TIII(4) may have been partially and fully occulted.A number of type II patches corresponding to the ones observed by Parker Solar Probe were also observed.TII(1) and TII(2), which are indicated by the red and yellow rectangles in Fig. 10, were observed to be nearly as intense as in Parker Solar Probe.However, TII(HB) was considerably weaker (marked by the orange rectangle).Such variations in intensity may indicate that the source directivity was in the direction of the spacecraft, which observed the brighter emission.In this case, the faintness of TII(HB) in STEREO A observations further supports that the source of the herringbones was likely located close to the line-of-sight of Parker Solar Probe and at the periphery of STEREO A.
Wind was the furthest spacecraft from the flare location and therefore only observed the low frequency parts of TIII(1) and TIII(3).Most other bursts were either too faint or not observed at all by the spacecraft.An interesting feature here is that Wind observed very faint signatures of both TII(1) and TII(2) as indicated by the rectangles in Fig. 10.Their fluxes however were an order of magnitude smaller than the ones observed by STEREO A. Considering that TII(1) was observed by ground based instrumentation, it is likely that the source of the emission was visible from Earth and therefore for Wind as well.
The multi-vantage point observations also introduce the phenomena of time delay (light travel time to spacecraft).By combining the time delay and the intensity variations between different spacecraft, it is possible to locate the spatial position of the source at a given frequency.We present a detailed analysis of the radio source propagation in Appendix B. Figure B.1 shows the radio source locations of TIII(1) and TIII(3) estimated using a directivity model.The results of the analysis suggest that TIII(1) propagated in the longitude −121.0 • ± 3.2 • (slightly east of the flare longitude), and the electron beam generating TIII(3) propagated in the longitude −98.3 • ± 4.1 • (slightly west of the flare longitude).This fits the by-eye analysis of the radio bursts based on their visibility to each observer.

Interplanetary context and SEP observations
The heliospheric conditions through which particles and CME-driven shocks propagate at the time of their release can significantly affect the SEP timing and intensity profiles (e.g.Laitinen et al. 2013;Dalla et al. 2020;Lario et al. 2022a).We use both multi-point solar wind and IMF observations and the WSA-ENLIL+Cone model (Odstrcil et al. 2004) to provide a comprehensive understanding of the geometry, not only of the interplanetary structures and their possible influence in the propagation of the SEPs, but also of the shocks and their role in forming the observed intensity-time profiles.In this section, we first describe the ENLIL simulation and then discuss in-situ plasma, magnetic field, and multi-spacecraft SEP observations.

The state of the interplanetary medium as derived with the ENLIL model
ENLIL is a global 3D MHD model3 that provides a time-dependent background characterization of the heliosphere outside 21.5 R .ENLIL uses time-dependent magnetograms as a background, into which spheroidal-shaped high-pressure structures without any internal magnetic field can be inserted to mimic observed CME-associated solar wind disturbances.ENLIL-modelled CMEs have an artificially higher thermal pressure to compensate for the lack of a strong magnetic field (Odstrcil et al. 2004, and references therein).To improve the characterization of the heliosphere, multi-point coronagraph observations are used to infer CME parameters, using the GCS model described in Sect.4.2.The inner boundary condition is given by the WSA V5.2 model, using inputs from the standard quickreduce zero-point corrected magnetograms from GONG (GONGZ), available on the National Solar Observatory website4 .In this case, the GONGZ magnetograms fit the in-situ solar wind speed and magnetic field polarity better (not shown).The reliability of the CME arrival predictions depends strongly on the initial CME input parameters, such as speed, direction, and width (Lee et al. 2013;Mays et al. 2015;Kay et al. 2020), but also on the errors that can arise in the ambient model parameters and on the accuracy of the solar wind background derived from the magnetograms (Lee et al. 2013).Based on Wold et al. (2018), the mean absolute arrival-time prediction error is 10.4 ± 0.9 hours, with a tendency to an early prediction of −4.0 hours.
The magnetic connectivity at the onset time of the SEP event is relevant to the understanding of the SEP observations, and considering the ENLIL-modelled varying solar wind conditions to calculate the IMF lines is an alternative to using the nominal Parker spirals.The preconditioning of the heliosphere and the interaction of the IP structures with the ambient solar wind that might be present at the SEP onset time can actively influence this connectivity (Masson et al. 2012;Palmerio et al. 2021;Lario et al. 2022b).Therefore, we choose an ENLIL simulation time from 10 to 24 April 2021 (i.e. from seven days before to seven days after the SEP event onset).This interval encompasses possible previous CMEs as well as subsequent CMEs propagating through the structured solar wind streams up to a distance tures present near Solar Orbiter, Earth, STEREO A, and Mars at the time of the onset of the particles that might modify the magnetic connectivity and SEP propagation conditions.There is a relatively small ICME reaching Solar Orbiter during the ongoing SEP event.According to ENLIL, this ICME does not extend to any other investigated spacecraft.Ahead of this ICME there is a clearly wider eruption covering about 140 • in longitude.Its western edge encloses STEREO A at the time of the SEP injection from the Sun and it is between the Sun and Mars.None of the CMEs inserted into ENLIL impacts Earth, but the leading edge of a stream interaction region (SIR) is reaching the planet at the time of the onset of the particle intensity increase seen at Earth.
The ENLIL simulation also shows that at the time of the initial SEP injection from the Sun (left panel of Fig. 11) the IP medium is relatively undisturbed between the Sun and BepiColombo as well as Parker Solar Probe.We note that the wide ICME discussed previously crossed BepiColombo and Solar Orbiter, but this was before the SEPs were injected at the Sun.Nevertheless, this ICME may still have an effect on the propagation conditions of SEPs.The simulated status of the heliosphere around the SEP onset time agrees overall with the in-situ plasma and magnetic field measurements as discussed further below.The right panel of Fig. 11 shows the heliosphere two days later, on 19 April 2021 at 21:00 UT.The ICME that was associated to the SEP event has then reached BepiColombo and Solar Orbiter.The simulation suggests that the ICME nose propagates between these two spacecraft and both of them cross the structure near the flanks.
The five bottom panels of Figs. 12, 13, and 14 present the in-situ plasma and magnetic field data over-plotted with the pink line showing the result of the ENLIL simulation from 17 April to mid 23 April.The whole set of panels in these figures present, from top to bottom, energetic electron intensities at different energies (1), proton/ion intensities at different energies (2), the magnetic field magnitude (3), the magnetic field latitudinal (4) and azimuthal (5) angles in spacecraft centered radial-tangential-normal (RTN) coordinates, namely θ B-RTN and φ B-RTN , the solar wind proton speed (6), and the solar wind proton density (7).As specified in the following section, ENLIL follows the general trend of the measured solar wind speed at the locations of Solar Orbiter and Mars, which were only separated by 9 • in longitude during the SEP event, while at the remaining locations there are some differences with in-situ measurements.Although ENLIL reproduced the overall features of high-speed streams present in the heliosphere during the period of study, the differences between the modeled and measured solar wind profiles could be explained by complex coronal holes that render the comparison between measurements and ENLIL results at Earth difficult, as well as preventing accurately resolving the glancing encounter of the SIR at STEREO A, as discussed below.
ENLIL successfully predicts the arrival of the several ICMEs observed in situ within the uncertainty of the model, as shown in the increase of the speed, density, or magnetic field in the pink profiles in Figs. 12, 13, and 14.Due to the absence of the internal magnetic field in the simulated CMEs, the magnetic field magnitude increase is, however, lower than what was measured in situ.In particular, based on ENLIL simulations and measured in situ as discussed in Sect.5.2, the ICME related to the SEP event is in-tercepted by BepiColombo and Solar Orbiter, while Mars might be only observing the associated IP shock, not the ICME ejecta.We relate the better simulation of the arrival time of this ICME at Solar Orbiter location in comparison with BepiColombo to the fact that we chose the CME parameters which better reproduced the portion of the CME oriented towards Solar Orbiter as ENLIL input, as discussed in Sect.4.2.The minimum longitudinal extent of the ICME related to the SEP event is ∼45 • , as shown in the right panel of Fig. 11.This value is in agreement with the angular extent of the CME along the equatorial plane (∼46 • ) estimated from the GCS reconstruction presented in Sect.4.2.

Multi-spacecraft in-situ plasma, magnetic field and SEP observations in context with the ENLIL simulation
As discussed in Sect.5.1, during the period of study there are several IP structures impacting the spacecraft under consideration, which may in turn affect the SEP particle profiles.In the following, we discuss the energetic particle observations and their relation with the interplanetary context.

BepiColombo
The two top panels of Fig. 12 (left) show the SEP event as observed by BepiColombo detected over a broad energy range by the MPO/SIXS, MPO/BERM, and Mio/SPM instruments with the time of the flare onset marked by the arrows at the top.Panel (1) shows the impulsive energetic electron event that reaches energies of at least 2 MeV.The proton intensity-time profile (panel 2) shows a more gradual increase.The Mio/SPM observations show that the event was observed even at proton energies > 200 MeV.There is no plasma information available, but at the time when particle intensities started to increase, the solar wind speed given by ENLIL is ∼400 km s -1 , as shown by the pink line in panel (6).
Commonly, the in-situ identification of the passage of ICMEs is based on a set of signatures typically observed in magnetic field and plasma data as well as some other proxies, such as bi-directional suprathermal electron (BDE) profiles (e.g., Zurbuchen & Richardson 2006;Kilpua et al. 2017).BepiColombo lacked plasma data, but the in-situ MPO-MAG magnetic field observations in panels (3)-( 5) do not show any evidence of typical ICME signatures (e.g., enhanced field, low field variety, coherent field rotation) during the onset and rising phase of the SEP event.This agrees with the previously discussed ENLIL simulation results and confirms that there was no large-scale solar wind structure at BepiColombo that could have directly influenced the SEP time profiles.
The increase in the magnitude of the magnetic field observed by BepiColombo on 19 April marks the arrival of the ICME related to the SEP event.The IP shock arrives at 11:40 UT (vertical solid line), while ENLIL simulates the ICME arrival time ∼6 hours earlier.Unfortunately, SIXS has a data gap at that time so that a potential low-energy particle intensity response to the shock passage (i.e., an energetic storm particle event) could not be studied in detail.However, BERM fluxes at ∼1.5-5.9MeV do not show a sig- nificant increase at the shock or a response to the passage of the ejecta, which follows the shock (gray shaded area).
The leading edge of the ejecta was observed at 13:57 UT on 19 April identified by a change in the magnetic field polarity along with the presence of coherent and organized magnetic field.Specifically, we observe a smooth and monotonic change of the magnetic field latitudinal and azimuthal angles shown in panels (4) and ( 5) that lasted until 20 April 00:04 UT.No other structures are observed until the end of the period shown in Fig. 12.

Parker Solar Probe
Panel (6) of Fig. 12 (right) shows that at the time of the SEP event onset the solar wind speed at Parker Solar Probe is ∼320 km s -1 .The SEP event has a very impulsive time profile both in the electrons shown in panel (1) and in the protons shown in panel (2).Compared with BepiColombo the event has a shorter duration, namely a faster decay.PSP/EPI-Hi/HET observes observes intensity increases at electron energies above 2 MeV and protons above 50 MeV.
Based on the plasma and magnetic field data given by the SWEAP and FIELDS instruments, no IP structures can be identified during the whole period shown in the right column of Fig. 12 location of Parker Solar Probe.This is in agreement with the ENLIL simulation results.

Solar Orbiter
Panels (1) and (2) of Fig. 13 (left) show the SEP event observed by Solar Orbiter.While the electron event is ob-served to reach energies up to ∼1 MeV, it is not as impulsive as the event observed by BepiColombo and Parker Solar Probe but shows more of a plateau-like profile.The intervening structures present at the time of the SEP onset, as suggested by ENLIL (Sect.5.1), might be associated with this behavior as they might hinder the SEP transport.This may also be the reason for the low anisotropy observed at the onset of the event as described in Sect.5.3.3.The energetic ion observations by Solar Orbiter/EPD/EPT allow us to discern the initial phase of the event only at energies 400 keV, reaching energies up to ∼60 MeV as observed by EPD/HET.
The solar wind speed at the time of the electron event onset is ∼380 km s -1 , as shown in panel (6), measured by the SWA instrument, which is well reproduced by ENLIL (pink line).In a later phase of the SEP event, Solar Orbiter observes several IP structures identified using MAG, SWA and the RPW instruments on board Solar Orbiter.A first ejecta (first gray shaded area in the left column of Fig. 13) arrives at 05:24 UT on 18 April.While it does not affect the energetic electron intensity time profiles or the highenergy ( 2 MeV) ion intensity time profiles, it seems to have acted as a particle barrier for low-energy ions.Only after its passage, at 14:18 UT on the same day, an increase in the 400 keV energy ions is observed.These particles are likely associated with the solar eruption on 17 April.The ICME-driven shock associated with this eruption arrives at Solar Orbiter on April 17 at 20:20 UT on 19 April.The shock simulated by ENLIL arrives a about ∼30 minutes later than the measured shock as shown by the pink line in panels (3), (6-7) in the left column of Figure 13.In-situ plasma and magnetic field observations.The panels present, from top to bottom, the magnetic field magnitude, the magnetic field latitudinal and azimuthal angles, θ B-RTN and φ B-RTN , the solar wind speed, and the proton density, where RTN stands for radial-tangential-normal coordinates (IP structures as described in top panel).The pink lines represent the ENLIL simulation results.
The shock obliquity (θ Bn , namely the angle between the shock normal and the upstream magnetic field) is estimated at Solar Orbiter using the magnetic coplanarity method (e.g., Paschmann & Schwartz 2000).A value of θ Bn ∼ (21 ± 5) • is computed, employing a systematic variation of the upstream and downstream averaging window lengths between 3 and 13 minutes, with the method described in Trotta et al. (2022).The lack of plasma data around the shock crossing limits further analyses.However, using the novel method introduced by Gedalin et al. ( 2021), an estimation for the Alfvénic Mach number using magnetic field only data yields M A ∼ 1.8, consistent with the fact that the shock passage has no significant influence on the energetic particle population at higher energies.
The low-energy ions keep rising until a peak is observed shortly before the shock passage, after which the intensities decrease.Right after the shock passage, another solar energetic electron event is observed on 20 April, which is not re-lated to the event under study but originated from an M1.1 flare in AR 12816 (at S24E25 as seen from Earth), peaking at 23:42 UT on 19 April.The ∼1 MeV/nucleon ions associated with this new injection showed a large enrichment of 3 He, with a 3 He/ 4 He ratio of, with a ratio of ∼5% (not shown).The second ejecta, which arrives at 06:51 UT on 20 April, corresponds to the ICME associated with the 17 April SEP event.It is marked by a smooth magnetic field and monotonic and coherent rotations in the magnetic field angles.While the second energetic electron event shows a depression of fluxes during the ejecta passage, low-energy ions show an enhancement inside the ejecta during its first half and a decrease during the second half.

STEREO A
Observations of the SEP event at STEREO A are shown in Fig. 13 (right).STEREO A observes a clear electron event, Fig. 13: In-situ SEP time profiles as well as plasma and magnetic field observations by Solar Orbiter (left) and STEREO A (right).Top: Energetic electron and proton temporal profiles observed from several energy channels.We use the sunward looking sectors of Solar Orbiters' EPD/EPT and HET.For STEREO A, as not all instruments provide sectored measurements, we use omni-directional data.The salmon shaded area indicates an SIR observed by STEREO A, while the stream interface is shown as a dashed line.Flare time and rest of IP structures indicated as in Fig. 12. Bottom: In-situ plasma and magnetic field observations.Solar wind densities for Solar Orbiter are obtained from RPW/QTN measurements.Panels as in Fig. 12. however only at near-relativistic ( 400 keV) energies.Similarly to Solar Orbiter, the proton event is only well observed at higher energies, namely 2 MeV.The maximum energy of the proton event seems to be lower than at the other spacecraft, barely reaching ∼25 MeV.At lower ion energies, SEPT observes an |bf enhanced pre-event intensity background, (most likely due to the SIR as described below), which might mask the SEP event.
As shown by the magnetic field and plasma data in panels (3)-( 7), the SEP event onset takes place during the passage of an SIR at STEREO A, which is indicated by the salmon-shaded vertical bar at the beginning of the time interval displayed in the right column of Fig. 13.The signatures indicate a glancing cross of the SIR structure, with only a very modest increase of the solar wind speed.The speed rises from ∼400 to ∼450 km s -1 ; sudden changes of the magnetic field polarity close to the stream interface (dashed vertical line), and drops in the magnetic field strength together with temperature increases (not shown) and proton density enhancement, which suggests that local reconnections are occurring.The ENLIL simulation also suggests an SIR arrival (not shown), but several hours earlier than observed, and infers a clearer intersection of the high-speed stream with the spacecraft.At STEREO A, no signatures of ICMEs are detected.This is in agreement with the ENLIL results.
The lowest-energy ion channels of SEPT show a clear variation in their intensities happening right after the stream interface.At the same time, the thermal proton density drops.The energetic electron increases observed after the data gap are associated to later SEP events that are not related to the event under study.passage, as measured by the MFI and SWE instruments on board Wind.The arrival at Earth of the high-speed stream simulated by ENLIL arrives a few hours later than that actually measured.The rear boundary of the SIR is difficult to define, as there is not a clear step-like speed increase and the dynamic pressure does not show any clear peak.The reason behind this behavior might be that the solar wind arriving at Earth originates from multiple and complex coronal holes.There is a big southern coronal hole extending to the equator and some low-latitude large patchier holes (not shown).The solar wind speed at the onset of the particle event is ∼600 km s -1 , as shown in panel (6).Panel (1) shows that only SOHO/EPHIN, which has a very low instrumental background, observes a clear but very gradual electron event at 0.25-0.7 MeV.The lower energies covered by Wind/3DP are showing an enhanced background that possibly masks the SEP event and may contain also ion contamination.This enhanced background, likely caused by the SIR, also dominates the low-energy ion observations by Wind/3DP.However, SOHO/EPHIN and ERNE show a proton event extending into the deca-MeV range, which is small, gradual, and clearly delayed with respect to the time of the flare.

Mars
On 17 April 2021 Mars was located at a heliocentric distance of 1.6 au 22 • west of the flaring active region at 225 • Carrington longitude.The top two panels of the right column of Figure 14 show ∼60-210 keV electron and ∼70-7000 keV proton intensities as measured in different energy channels of MAVEN/SEP.Panels (3)-( 5) show only the ENLIL simulations of the magnetic field, as no measurements are available.The solar wind speed (6) and density (7) measurements by MEX/ASPERA-3/IMA are rather sparse, however, they show an overall good agreement with the ENLIL simulation for the solar wind speed.In this case, we also show the dashed pink lines corresponding to the background solar wind simulation, without including any CME.The separation of the solid and dashed lines indicate the effects produced by the passage of the interplanetary structures, based only on ENLIL results.A first interplanetary shock is modeled to arrive at 09:00 UT on 18 April.Two pre-SEP event ICMEs arrive at 04:00 UT on 19 April and 10:00 UT on 21 April, as simulated by ENLIL, that might be the same ICMEs measured earlier by Solar Orbiter.Lastly, the interplanetary shock related to the SEP event under study impacts Mars at 10:00 UT on 22 April, based on both the simulation and on the increase in solar wind speed and density measured in situ.According to ENLIL, the shock is, however, not followed by an ejecta.The ICME flank might therefore have missed Mars.
It is difficult to associate the energetic particle increases observed by MAVEN with the 17 April SEP event.Nevertheless, an electron increase is observed in the higherenergy channels right after the flare (marked by an arrow).Although the onset times are very hard to determine and might suggest a too-early onset to account for the expected travel time of these electrons, a potential SEP contribution cannot be excluded.More likely, on the other hand, is that the CME-driven shock associated with the event has contributed to the electron and proton increases observed on 22 April because the peaks of the SEP increases agree well with the shock arrival time simulated by ENLIL.However, another possible source of this increase could be the same new SEP event as observed also by STEREO A on 22 April (at S24E25 as seen from Earth), which is magnetically well-connected with Mars during the period under study (see Fig. 1).

SEP pitch-angle distributions and first arriving particles
Figure 1 (right) combines the SEP observations as measured by the five inner-heliospheric spacecraft and shows how strongly the event characteristics such as intensitytime profiles, onset times, and peak intensities vary from observer to observer.Given the well-separated positions of these spacecraft shown in Fig. 1 (left) and their varying separations with respect to the parent flare location, this is not unexpected.It has been found in earlier studies that the longitudinal distribution of peak intensities usually shows a decrease with increasing longitudinal separation angle from the associated flare longitude (e.g., Lario et al. 2013;Dresing et al. 2014;Richardson et al. 2014).These authors described the longitudinal peak-intensity distributions with Gaussian functions.However, due to the limitation of only three well-separated observers, these analyses suffered from large uncertainties.The new, larger spacecraft fleet will allow us to analyze these longitudinal distributions in a better way, however, instrument inter-calibrations, especially of the new mission's payload, are still pending.Looking at the ∼20-25 MeV proton peak intensities observed by each spacecraft (Fig. 1) we find, however, a deviation from the expected ordering of peak intensities with absolute longitudinal separation angle.Parker Solar Probe, which is slightly worse connected (|∆Φ| = 72 • ) than Solar Orbiter (|∆Φ| = 65 • ), observes not only a significantly higher-intensity event, but also a clearly more impulsive time profile.While the higher peak intensities at Parker Solar Probe could be explained by its smaller radial distance from the Sun as compared to Solar Orbiter (0.42 au vs. 0.84 au), the significantly different time profiles rather suggest a different connectivity to the SEP injection region.
In this section, we analyze the SEP observations in more detail to determine the timing of the first arriving particles, which allows us to relate the SEPs with their solar counterpart observations (see Sect. 6).Furthermore, we analyze pitch-angle distributions (PADs) to characterize the degree of pitch-angle diffusion, namely the importance of transport effects.

BepiColombo
BepiColombo detects the most intense event out of all observers.This is expected based on both its closer radial distance from the Sun (0.63 au) and its fairly good connection to the associated flaring active region -1 • (79 • ) longitudinal (total) separation angle between the flare site and the spacecraft magnetic footpoint (cf.Table 1).BepiColombo also observes the earliest SEP onsets, e.g.16:30 UT for 71 keV electrons, and the corresponding inferred injection times are the earliest out of all observers (see Sect. 6).Surprisingly, BepiColombo/SIXS detects a 5-minute earlier onset time for 71 keV electrons than for 960 keV electrons (see Table 4).Although these onset times almost agree within the error bars, the difference between the inferred injection times is significant.The much longer travel time of the lower energy electrons yields an 11-minutes earlier injection time compared to the ∼1 MeV electrons.The first 25 MeV protons are detected at 17:00 UT±4 min, which corresponds to an inferred injection time situated between that of the low-energy electrons and the high-energy electrons.However, given the error bars, it would agree with both of the inferred electron injection times.Figure 15 shows sectored energetic particle measurements by SIXS (middle panel) as well as the pitch-angles covered by the center of the four different viewing directions (top panel) and the intensity-PAD in the bottom panel.The left-hand figure shows ∼100 keV electrons and the right-hand figure shows 8 MeV protons.Although the event is anisotropic, as can be seen by the different intensity levels observed by the different sides of the SIXS instrument, a velocity dispersion analysis did not yield meaningful results neither for electrons nor for protons.Therefore, we can only apply the time-shift analysis to infer the particle injection times at the Sun at specific energy bands.The proton anisotropy is stronger and the anisotropic phase is clearly longer for protons, lasting about two hours.Unfortunately, the electron onset falls into a period of poor pitch-angle coverage of the sector where particles streaming from the Sun along the outward magnetic field would enter (pitch angle 0).The electron anisotropy could therefore be underestimated during the onset phase and this could lead also to a determination of too late electron onset times.

Parker Solar Probe
Based on the ∼25 MeV proton observations (see Fig. 1, right), Parker Solar Probe observes the second most intense event after BepiColombo.Because Parker Solar Probe's electron observations are not yet available in units of intensity, it is not possible to compare their intensity level with that of other spacecraft.While EPI-Hi does not provide the necessary time resolution to discern velocity dispersion in these relativistic electrons, the time resolution of EPI-Lo would be sufficient to discern velocity dispersion in the near-relativistic electrons.However, the limited statistics at these energies makes it challenging to conclude whether EPI-Lo observed electron velocity dispersion or not.Nevertheless, a small but significant anisotropy is present in the 90 keV electron observations, as denoted by the higher intensity of the sunward-looking wedge W3 (black) compared to the antisunward viewing wedge W7 (red).Still during the rising phase of the electron event, around 17:00 UT, we observe a second step, marked by the blue dashed line in the third panel of Fig. 16 (left), which is observed by both EPI-Lo and EPI-Hi.It does not correlate with any changes in the magnetic field and therefore seems not to be caused by a local effect.Later, around 17:30-18:00 UT, we observe a phase of stronger anisotropy consistent between EPI-Lo and EPI-Hi that appears to be tied to changes in the magnetic field.
Figure 16 (right) shows proton observations at 10 MeV (top panel) and 25 MeV (middle panel) in different viewing directions as provided by EPI-Hi/LET.As for BepiColombo, the proton observations show a stronger anisotropy than the electrons, which also lasts longer (>6 hours).In both the 10 MeV and 25 MeV timeintensity profiles, the sunward-facing aperture (LETA) shows the fastest onset and highest intensity.
In contrast to the electron observations, the protons show a clear velocity dispersion.Figure 17 shows a velocity dispersion analysis (VDA) that results in a path length of L = 0.63 au traveled by the protons and an inferred proton injection time at 16:46 UT±10 min.Even considering the uncertainties, this injection time is significantly later than those determined for electrons through a time shift analysis (TSA), using the same path length, that result in 16:26 UT (16:30 UT) for 920 keV (90 keV).

Solar Orbiter
Solar Orbiter's magnetic footpoint at the Sun is similarly far separated in longitude from the flare location to that of Parker Solar Probe, but it is situated on the other side, namely west of the flare.Solar Orbiter observes significantly lower proton intensities than Parker Solar Probe (see Fig. 1, right).Furthermore, in contrast to BepiColombo and Parker Solar Probe, who observe an impulsive proton time profile, Solar Orbiter observes a gradual profile both in electrons and protons.While in case of CME-driven shock acceleration a more gradual time profile is expected for an observer situated to the east of the source region as compared to an observer situated to the west due to their different connections to the CME-driven shock front (e.g., Cane et al. 1988), the difference in their peak intensities is not expected to show such strong asymmetry (e.g., Richardson et al. 2014).However, Solar Orbiter's distance to the Sun, which is twice of Parker Solar Probe's distance, is expected to contribute to this intensity difference.
Energetic electron observations do not show any significant anisotropy, neither at lower energies as illustrated by the ∼100 keV electron PADs (Fig. 18, left) observed by Solar Orbiter/EPT, nor at MeV energies (not shown).The right-hand part of Fig. 18 shows the PAD of ∼8 MeV pro- tons as detected by Solar Orbiter/HET, which shows that the early phase of the MeV proton event is anisotropic for about seven hours, showing higher fluxes in the sunwardlooking telescope that corresponds to pitch angles near 180 • , consistent with the inward magnetic polarity (see also column 8 in Table 1).
To perform a VDA, we determine the onset time based on the proton time profiles in the HET sunward telescope.Therefore, we use the energy channels between 7 MeV and 45 MeV and reconstruct the energy bins by combining every three proton channels.We then apply the Poisson-CUSUM method (Huttunen-Heikinmaa et al. 2005) and derive the onset time in those new channels.The VDA for protons is based on these onset times and results in an inferred injection time at 17:14±12 min and a path length of L = 1.24±0.18au.We display the results in Fig. 19 where we overplot the resulting VDA fit on the dynamic proton spectrogram.
For electrons, no VDA was possible because the onset times of many different energy channels were basically the same.A possible reason could be the rather poor pitchangle coverage during the event onset not covering the direction along the magnetic field.We therefore determine inferred injection times for selected energy channels based on TSA assuming the same path length as derived from the proton VDA (see Table 4).

STEREO A
STEREO A is a far-separated observer with 129 • (109 • ) of longitudinal (total angular) separation between the flare location and the spacecraft's magnetic footpoint at the Sun computed with ADAPT-WSA.It is therefore not surprising that the SEP event at STEREO A is less intense than those explored so far, and that the intensity-time profiles are more gradual and isotropic (e.g., Dresing et al. 2014).Figure 20 (left) shows the electron PAD observed by STEREO A/SEPT at ∼100 keV, which shows no anisotropy except during the time of the maximum, where the intensity in the anti-sunward sector is slightly higher.We note that, since the spacecraft was put upside-down after the superior solar conjunction in 2015, the Sun and anti-Sun sectors no longer point along the nominal Parker spiral but perpendicular to it.Figure 20 (right, third panel from top) shows the 4-6 MeV proton intensities observed in the 16 sectors of STEREO A/LET.The top panel shows the color-coded intensity-PAD and the second panel shows the pitch-angles of the sector centers.The statistics in the single sectors are poor, which is why the bottom panel shows averaged intensities of the eight A and B sectors, respectively, and in black an average of all sectors.Interestingly, LET shows a double-peak time profile with the first peak, starting shortly after 18:00 UT on 17 April, being much more anisotropic than the second peak as almost no intensity is yet observed in the B-side sectors of LET.The depletion between the peaks at ∼6 UT on 18 April is not caused by poor pitch-angle coverage.Indeed, the pitch-angle coverage is better during this phase than during neighboring periods.As shown in Fig. 13 (right), there is no clear interplanetary structure that can be associated with this dip.We therefore argue that it is either caused by a change of the magnetic connection to the parent source region or a distinct new particle injection, which is also supported by the differently strong anisotropies during both peaks.
Due to the gradual nature of the event and rather poor statistics, it was not possible to apply a VDA, and in order to determine onset times we had to average the data, leading to significant uncertainties.We obtain an onset at 18:25±10 min for 85-125 keV electrons and 19:30±1h for 13.6-23.8MeV protons.Assuming a scatter-free propagation along a nominal Parker spiral with a length of 1.16 au, this would translate to inferred injection times of 18:08±10 min for the electrons and 18:20±1 h for the protons.The event is also less energetic at STEREO A. Different to all other inner spacecraft, STEREO A does not detect electrons in the MeV range and the event at 25 MeV protons is very weak (see Fig. 1, right).However, this could be also due to instrumental differences with STEREO A/HET being less sensitive.

Near-Earth Spacecraft
As already discussed in Sect.5.2.5, the SEP event at the Sun-Earth L1 point is only observed at high energies, both in electrons and protons.However, as the event is weak and very gradual, no VDA was possible and the determined onset times (see Table 4) suffer large uncertainties.Altogether, the observations of a gradual, delayed, and small event at Earth suggest that the event was only observed due to perpendicular particle diffusion (e.g., Dresing et al. 2012) since there was probably no direct magnetic connection to a source region.

Combined timing analysis and implications on the sources of the SEP event
Table 4 presents a timeline of the main features of the 17 April 2021 SEP event with all observations times at the different spacecraft (column 2) being shifted to the Sun (column 1).This means, remote-sensing observations are corrected for the varying light travel times based on the different spacecraft distances, and energetic particle onset times are used to infer the corresponding injection times at the Sun.Studying the SEP event as observed by the multiple spacecraft implies the use of a multitude of instruments that provide different energy ranges and channel widths, different instrumental backgrounds, varying signal-to-noise ratios based on their locations with respect to the SEP source region, as well as different local interplanetary conditions that can influence the SEP observations (see Sect. 5).All these factors make a comparison of the SEP observations at the different spacecraft challenging.For example, for many spacecraft locations or species, a VDA was not possible (see Sect. 5.3).In these cases, we apply a simple TSA to infer the SEP injection times and use an energy channel showing a clear onset time.This implies that we sometimes have to use different energy ranges to infer SEP injection times.
Table 4 (see also Fig. 21) shows that the first feature of the event was the start of the HRX flare, which happened already at 16:00 UT.The flare then continued for almost an hour showing multiple HXR peaks with the last and strongest one at 16:44:30 UT (see Sect. 4.1).The first radio type III burst onset (TIII(1)) was only observed 16 minutes after the start of the flare, followed by the first type II (TII(1)) burst observation at 16:18 UT (see also Sect. 4.5).Around this time (16:18±4min), the first SEPs were inferred to be injected towards BepiColombo as derived from the 71 keV electron onset.The protons of about 25 MeV were injected were injected later towards BepiColombo at 16:25±4min, temporally situated between TIII(1) and TIII(2), at the end of TII(1).BepiColombo's SPM instrument on board Mio was even able to detect > 200 MeV protons, which arrived, however, significantly delayed with an inferred injection time at 16:43±5min.Surprisingly, we determine a significant later injection time for ∼1 MeV electrons at 16:29±1min, as compared to the 71 keV electrons, which happened during TIII(2).Due to the impulsive, highintensity event observed by BepiColombo/SIXS the onset times are well-defined carrying only small uncertainties, which are assumed to be the same for the inferred injection times.This strongly suggests that not only the electrons and protons observed at BepiColombo are related to different injection episodes, but also that the near-relativistic and relativistic particles suggest different injection times, a feature which was also observed during the 9 October 2021 SEP event (Jebaraj et al. 2023a).
Figure 21 illustrates the inferred SEP injection times (vertical lines) in comparison with the HXR flare observations taken by Solar Orbiter/STIX and the radio observations by PSP/RFS.In contrast to Table 4, Fig. 21 only displays those injection times that were inferred to happen during the early phase of the event, namely during the radio active phase.Therefore, only times corresponding to Bepi-Colombo, Parker Solar Probe, and Solar Orbiter, the three best-connected spacecraft, are included.
In the case of Parker Solar Probe, we find later injection times compared to BepiColombo and a significantly earlier injection of electrons compared to protons.Both relativistic and near-relativistic electrons are inferred to be injected during TIII(2) at 16:26±5min (∼1 MeV) and 16:30±5min (∼90 keV).Because TIII(2) was found to be strongly directed towards Parker Solar Probe, this association is not surprising.We do not find evidence of SEPs related with TIII(1) to reach Parker Solar Probe's location.However, the inferred injection time of a step-like feature in the rising phase of Parker Solar Probe's electron event (see Sect. 5.3.2) shows a temporal correlation with TIII(4), the second type III burst, which shows a strong directivity towards Parker Solar Probe (cf.Sect.4.5).The clear velocity dispersion observed by Parker Solar Probe for deka-MeV protons yields an injection time at 16:46±10min, which is temporally situated between TIII(3) and TIII(4).This sug- gests, similar to BepiColombo observations, that electrons and protons were injected during different episodes and could possibly be related to different acceleration mechanisms and locations.
As discussed in Sect.5.2.3 the interplanetary conditions between Solar Orbiter and the Sun were disturbed by minor transient events, likely affecting the SEP transport and leading to the comparatively lower peak intensities and less well-defined onsets.This could also lead to delayed SEP onsets and consequently yield too-late-inferred injection times, especially when using TSA, which we do for the electrons observed by Solar Orbiter.Nevertheless, using the same path length as derived from proton VDA we find the same pattern of electrons inferred to be injected earlier at 16:41±15min (∼1.6 MeV) and 16:55±2min (∼100 keV) than protons for which we were able to perform a VDA yielding an injection time at 17:14±12min for protons between 7 and 45 MeV.Only the 106 keV electron injection time could potentially be related to a radio feature, that is, TIII(4).The inferred proton injection time is about 20 minutes later than the last type III burst (TIII(4)).
The onset times of SEPs at STEREO and Earth are so delayed and uncertain that we cannot infer a direct connection with any of the early activity phenomena of the event, as shown in Fig. 21.Furthermore, the injection times of the three best-connected spacecraft (BepiColombo, Parker Solar Probe, and Solar Orbiter) spread already over the whole radio active time period of about 40 min.This implies that all four type III radio bursts could mark distinct SEP injections that have contributed to the global multi-spacecraft SEP event.The different directions of these radio bursts (see Sect. 4.5) furthermore opens the possibility that the multiple injection episodes were differently important for the different observer locations.
The vertical shaded regions in Fig. 21 denote the times (including uncertainties) at which a magnetic connection with the CME-driven shock was established with each of the five inner-heliosphere spacecraft according to the analysis reported in Sect.4.4 and summarized in Table 3.Although we find the shock to potentially connect already early and at low heights with all the five spacecraft locations, several inferred injection times happened already before, suggesting that the shock was not the main accelerator of these For STEREO A and Earth, the shock connection times happened earlier than any inferred injection times, but because of the strongly delayed and uncertain onset times, it is not possible to pin down a clear role of the shock against the potentially involved transport effects.

Interplanetary transport modeling
In this section, we present simulations of the interplanetary transport of SEPs using the spatially 2D model of Strauss & Fichtner (2015).Simulations are performed for 150 keV electrons to qualitatively illustrate the transport concepts discussed in this work like the role of perpendicular diffusion vs. a direct magnetic connection to the source region.
As input to the model we implement the pitch-angle and perpendicular diffusion coefficients used by Strauss et al. (2017Strauss et al. ( , 2020) ) that are based on fundamental turbulent quantities and optimized to reproduce the catalog of widespread events from Dresing et al. (2014).The top left panel of Fig. 22 shows a contour plot of the omni-directional intensity of 150 keV electrons, calculated from this model when assuming a single SEP source, at five hours after particle injection.Here, the position of the different spacecraft are shown, along with their magnetic connectivity to the inner model boundary, assumed to be the Alfvén surface approximately located at r ∼ 10 R (0.05 au).The dashed magnetic field line connects to the maximum of the injected SEP distribution, which was assumed to be a Gaussian spatial distribution with a broadness of 5 • .Below this, we show the temporal profile of the omni-directional intensity at the different spacecraft, the assumed profile of the SEP source, and lastly the corresponding particle anisotropies.As in previous work, we assume a Reid-Axford profile for the SEP source with an acceleration time of 0.1 hours and an escape time of one hour.For the temporal profile of the differential intensity two sets of solutions are shown: The first, in dashed lines, are model solutions using the default set-up of Strauss et al. (2017), while the solid curves are solutions where the parallel mean-free-path is reduced, in an admittedly ad-hoc fashion, by a factor of 5.This is done to account for the possibly more disturbed nature of the inner heliosphere during this event as discussed in Sect.5, in contrast to the basic model that assumes quiet solar minimum conditions.Of course, a smaller parallel mean-freepath, namely more pitch-angle scattering, leads to a slower rise to maximum, a slower decay phase, and a smaller level of anisotropy.The exact levels of the transport parameters appropriate to reproduce this specific event will be the topic of a future, more detailed, modeling endeavor.Interestingly, and the main conclusion from this basic setup, is the fact that the model, while reproducing the general trends observed by most of the spacecraft, consistently under-estimates the SEP intensity at the position of the Parker Solar Probe spacecraft for a range of transport parameters.As a possible remedy for this discrepancy, the right panels of Fig. 22 show the modeling scenario of multiple SEP sources, releasing particles into the inner heliosphere at different positions.The dotted and dashed lines in the top panel show the position of these four injection sources where these are chosen to approximately correspond to the inferred position of the four observed radio bursts from Sect.4.5.The magnitude of the four injections are, how-ever, not well constrained and chosen here rather arbitrarily to roughly correspond to the measurements.The normalization of these injections are chosen such that the total fluence of the electrons introduced into the heliosphere is the same for the left and right panels.For the default model set-up (i.e., the dashed curves), the different injections are visible in the calculated temporal profiles of the magnetically well-connected spacecraft, while any such prominent peaks disappear for the case of more pitch-angle scattering.More importantly, the level of the simulated profile at Parker Solar Probe is also now more consistent with the observations.

Discussion
The SEP event associated with the flare-CME on 17 April 2021 was observed by five well-separated spacecraft in the inner heliosphere, with additional constrains provided by observations at Mars.The multi-vantage point observations portray a complex picture of the event, which involves significantly different characteristics of both the energetic electron and the proton/ion event and an asymmetry in the longitudinal distribution of their intensities.We find evidence that the reason for the wide SEP spread involves a number of different mechanisms with varying importance for different vantage points, which we will discuss in the following.
While the associated CME was relatively slow and narrow (speed: ∼880 km s −1 and width: ∼46 • , see Sect.4.2) as compared to other widespread SEP events with high-energy particles (e.g., Lario et al. 2017;Kouloumvakos et al. 2019), the solar flare emission in HXR was exceptionally longlasting (one hour) and complex (cf.Sect.4.1).The radio event was also equally long-lasting (Sect.4.5) and showed four distinct type III burst groups indicating particle injection episodes over a period of about 40 min.Several type II bursts were also observed, which suggests particle acceleration at the different flank regions of the associated shock.Although the event was observed by a fleet of five wellseparated spacecraft, no full optical coverage of the solar surface was available, leaving some sectors of potential flare locations unobserved.However, our comprehensive analysis of the available X-ray, EUV, white-light, and radio observations suggests that the event was caused by activity related to a single source active region at the Sun (cf.Sect.4).
The earliest SEP onsets and inferred injection times are found for BepiColombo, which was the best-connected spacecraft to the flaring active region.Significantly further, but similarly far separated from the active region in heliolongitude were Parker Solar Probe (east of the flare) and Solar Orbiter (west of the flare).However, both spacecraft observed dramatically different SEP characteristics suggesting a longitudinal asymmetry: Parker Solar Probe observed a more intense and impulsive event while Solar Orbiter observed a more gradual, less intense and delayed event.While the intensity difference could also be explained by the different radial distances of the spacecraft, the strongly different intensity-time profiles rather suggest a different magnetic connection to the source region, which could be different portions of the CME-driven shock front (e.g., Cane et al. 1988) and/or combination of differently directed SEP injections.
Our detailed radio analysis reveals that two of the four observed TIII radio burst episodes were directed signifi-Fig.22: Transport modeling results for 150 keV electrons.The left panels represent the standard case of a single SEP injection into interplanetary space, while the right panels are for multiple injections.The top panels are normalized contour plots of the SEP intensity at five hours after the initial injection, while the bottom panels show the resulting particle intensities, as a function of time, at a number of spacecraft positions, the different SEP injections, and the resulting particle anisotropies.More details are given in the main text.cantly more towards Parker Solar Probe as compared to the other radio bursts.These two injection episodes are likely the main source of the electron event at Parker Solar Probe, which is also supported by the inferred injection times (cf.Sect.6) and the results of the transport modeling (Sect.7).These injections may also have contributed to the Parker Solar Probe proton event as we infer an earlier injection of protons than the shock connection time (cf.Sect.6).This would also explain the comparatively high intensities observed at Parker Solar Probe when compared to Solar Orbiter.However, the long-lasting anisotropy of the proton event at Parker Solar Probe of about nine hours (see Fig. 16), which is not observed in the case of the electrons, suggests a long-lasting proton injection, most likely related to the CME-driven shock.A similar picture is presented by Solar Orbiter, at which a clear anisotropy, lasting about seven hours, was observed for protons but not for electrons (Fig. 18), suggesting that the shock played an important role in creating the proton event at Solar Orbiter.Also, BepiColombo observes a long-lasting anisotropy for protons (about seven hours), while electrons only show significant anisotropies during the rise phase of the event (Fig. 15), which is in agreement with a short, likely flare-related, injection.While an initial flare contribution to the proton event cannot be excluded for BepiColombo and Parker Solar Probe, a later shock contribution is, therefore, likely.
An even less intense and more delayed SEP event was observed by STEREO A and near-Earth spacecraft, which were far separated in heliolongitude (126 • for STEREO A and 144 • for Earth, see Table 2).At both locations, no significant anisotropies were observed for electrons and only at STEREO A a short anisotropic period during the early phase of the event was visible in LET proton measurements (Fig. 20).Although we determine potential shock connection times for both positions already before any inferred injection time of SEPs reaching STEREO A or Earth, the observed significantly lower intensities and missing anisotropies suggest that no direct connection with the shock nor flare-related injection was established but rather that perpendicular diffusion was involved in distributing the SEPs.However, the short anisotropic period observed by STEREO/LET for 4-6 MeV protons could be the trace of a shorter-lasting connection to the shock.The presence of interplanetary structures such as pre-event ICMEs and the SIRs at the spacecraft locations could further have modified the magnetic field topology and enhanced scattering conditions, which likely contributed to the unclear and delayed SEP onset times.However, a wider SEP injection region either provided by the extended shock front of 180 • or the four type III burst-related injection episodes, which covered a longitudinal angle of about 110 • , may have been a key ingredient in producing the widespread SEP event reaching also STEREO A and Earth.Especially TIII(2) and TIII(4), marking injections close to the longitudinal location angle of Parker Solar Probe (see Sect. 4.5), would provide another, significantly closer injection region for the Earth's location, facilitated over the western limb, as compared to the location of the flare.This could explain the comparatively early electron onset time detected by SOHO/EPHIN around 22:00 UT on 17 April.
At Mars, which was magnetically well-connected with STEREO A during the onset of the event (see Fig. 1, left), we did not observe a clear SEP increase associated with the early phase of the event.However, the associated CME-driven shock could have reached Mars on 22 April and an energetic particle increase was observed, which could have been related to the shock.

Summary and conclusions
The 17 April 2021 SEP event is the second widespread event of solar cycle 25 and the first one that was ever observed by five well-separated space missions in the inner heliosphere (within 1 au) constrained also by observations at Mars.It is an energetic event showing electrons up to the MeV range and 25 MeV protons reaching all inner heliospheric spacecraft positions, which span a longitudinal range of 210 • .BepiColombo observations by Mio/SPM even show the presence of >200 MeV protons.The closest observer to the Sun was Parker Solar Probe (r = 0.42 au) followed by BepiColombo (r = 0.63 au) and Solar Orbiter (r = 0.84 au).As outlined in Sect.8, the interplanetary SEP event was likely formed by a combination of different processes with varying importance at different spacecraft positions.For instance, the observations suggest a different origin of the electron and proton SEP event.This is most clear in the case of the three best-connected observers, BepiColombo, Parker Solar Probe, and Solar Orbiter: at Parker Solar Probe and Solar Orbiter we find significantly earlier inferred injection times for electrons (at all energies) than for ∼25 MeV protons.At BepiColombo only the nearrelativistic electron injection is found to be earlier than that of the protons.Also the much longer-lasting anisotropies observed in the proton event compared to electrons suggest an extended injection for protons only.Furthermore, different spacecraft were likely fed by different injections related to the various radio features with different injection directions as suggested by the radio directivity analyses (see Sect. 4.5).The timing analysis (see Sect. 4) shows that BepiColombo detected electrons injected already during the first episode (TIII(1)), while Parker Solar Probe likely detected only electrons from the later episodes, mainly from TIII(2) and TIII(4), which were directed towards Parker Solar Probe.
A possible alternative source of the two type III groups, namely TIII(2) and TIII(4) observed by Parker Solar Probe, could be the shock wave.In Sect.4.5, we show that TIII(2) and TIII(4) were strongly polarized.This high degree of polarization indicates that the source is a region with strong magnetic fields.A highly compressive shock wave may provide such conditions where the electron beams are accelerated via a shock drift acceleration mechanism (SDA; Ball & Melrose 2001;Mann et al. 2018).The energy gained by the electrons in these cases may also dependent on other factors such as the upstream electron distribution.If we were to assume that the thermal electrons (∼1% of the speed of light, Halekas et al. 2020) are being accelerated, then the maximum energy gain in a short period through SDA can be a factor of 14 which leads to a leads to 14% speed of light or ∼10 keV.However, a small portion of the tail electrons may be accelerated to higher energies.Our multi-spacecraft analysis further emphasizes that the location where the TIII(2) and TIII(4) originated was in the direction of Parker Solar Probe.There is a strong possibility that some of the type III bursts within TIII(4) were accelerated by the shock wave since they are observed to be emanating from TII(HB) (cf.Sect.4.5).However, the lack of meter-decameter wave measurements limits us from corroborating the generation of TIII(2) by the shock wave.In order to grasp the mechanisms of electron acceleration in the corona, full meter-decameter measurement would be necessary (cf.Jebaraj et al. 2023a).
Due to the not yet available Parker Solar Probe electron measurements in units of intensity, we cannot compare the electron intensity levels of Parker Solar Probe with other spacecraft.However, because the time profiles at Parker Solar Probe are much more impulsive and peaked for both electrons and protons compared to Solar Orbiter, situated at a comparable absolute longitudinal separation angle, it is plausible that the overall ordering of intensities observed at the different spacecraft is similar for electrons and protons.Based on this assumption we performed initial interplanetary transport modeling of this event for electrons (see Sect. 7) supporting the idea that SEPs were released from several different source longitudes.These show that without a SEP source near the magnetic footpoint of Parker Solar Probe, the measured intensity at that spacecraft cannot be reproduced by the model, independently of the adopted transport coefficients.While these simulations show promise, a more detailed modeling study is required, taking the disturbed nature of the interplanetary medium into consideration while sufficiently optimizing the transport coefficients used in the model.The detailed 3D modeling of the particle transport will be based on the treatment developed by Dröge et al. (2010), which employs focused particle transport along the large-scale heliospheric magnetic field as well as diffusion perpendicular to the field.We will also take into account the disturbances in the largescale magnetic field caused by preceding CMEs as comes out from the observations and predictions from the ENLIL model used in the current version of the paper.Such an extensive study is currently ongoing.
An important feature, which is not observed for the electron event, is the presence of long-lasting periods of proton anisotropies as observed by BepiColombo, Parker Solar Probe, and Solar Orbiter.This requires an extended, likely shock-associated proton injection.However, a flare contribution to the early phase of the proton event at these spacecraft cannot be excluded and is especially likely for Parker Solar Probe, for which the shock-connection time is determined to be established only after the inferred proton injection time.In the case when the scattering of protons, in particular through 90 • pitch angle, is reduced a long-lasting anisotropy can arise as well.
In the case of the two farthest separated observers, STEREO A at a longitudinal separation angle of 126 • , and Earth at 144 • , SEP intensities were significantly lower, showing a more gradual profile and significantly delayed onsets, which suggests that these observers did not establish a direct magnetic connection with any of the potential SEP source regions.Missing anisotropies together with the aforementioned characteristics suggest that perpendicular diffusion was involved in distributing the SEPs to these far separated longitudes.The presence of interplanetary structures such as the pre-event ICMEs and the SIRs may have contributed to modifying the magnetic field topology and enhancing scattering conditions leaving also room for a potential direct magnetic connection that was masked by a strongly disturbed parallel transport.However, even in the case of perpendicular transport being involved, we consider it likely that the widespread SEP observations were supported by an extended injection region.This could either have been provided by the extended shock front (∼ 180 • ) or by the different injection directions marked by the four radio type III burst episodes covering in total a longitudinal range of about 110 • .A likely evidence for an extended shock front is the presence of multiple type II radio bursts, namely TII(1), TII(2), and TII(HB) (see Sect. 4.5), which are emitted at different locations on the expanding shock front.Our analysis of the radio intensity and directivity suggests that the sources of TII(1) and TII(2) were directed towards Solar Orbiter and STEREO A, while that of TII(HB) was clearly directed towards Parker Solar Probe.
The study of the 17 April 2021 widespread SEP event allowed us to perform a comprehensive multi-spacecraft analysis combining remote-sensing and in-situ observations of six well-separated observer positions and taking full advantage of the various complementary data sets.The advanced spacecraft fleet enabled us to characterize signatures of a very complex SEP event, which would not have been possible with fewer observers.We were able to identify significant differences between the electron and proton SEP event, as observed by the different spacecraft, with a more likely flare association of the electron event and a more likely shock source for the proton event.However, a mixing of both cannot be excluded.Thanks to the position of Parker Solar Probe, we were able to observe otherwise hidden SEP features that highlight the role of significantly different injection directions of the four different injection episodes, which we consider a new scenario that has to be taken into account as a potential contributor to widespread events.
Future case studies of additional widespread events with the currently available spacecraft fleet will hopefully allow us to further characterize the necessary ingredients of widespread events and the different scenarios that are able to produce these rather rare events.backward exponential smoothing method (Brown 1963) in the left panels of Fig. C.1.The x-axis shows the distance (in Mm) from the eruptive center, which corresponds to the outer boundary of each ring, and the y-axis gives the resulting integral intensity in the considered angular sectors.The EUV wave is represented by the intensity increase in the perturbation profiles toward the peak amplitudes (big red markers), followed by a decay to background level.To obtain a robust estimate of the wavefront locations, we extract the outer distance at an intensity value corresponding to 30% of the peak value above the background (small orange markers).The right panels of Fig. C.1 show the obtained evolution of the wavefront locations as function of time.From the linear fits (solid lines) to the wavefront locations, which were followed up to about 600 Mm from the source region, we derive a mean velocity increasing gradually from sector 4 to sector 1 from 223 to 327 km s −1 .In sector 2, the wave propagation could be further seen up to about 680 Mm from the source region; however, we do not include it in our analysis because of the relatively weak signal at these distances.
To study the EUV wave kinematics above the solar limb and as function of height, we show in Fig. C.2 stack plots generated from STEREO A/EUVI 195 Å base-difference images in vertical slits at different heights, from 1.05 to 1.30 R as measured from the Sun center (cf., Veronig et al. 2018).To this aim, the emission is integrated over a vertical layer of a width of 18 EUVI pixels and stacked in time.In these stack plots, the EUV wave in the northern direction can be followed to a distance of about 740 Mm, propagating with velocities from 260 to 450 km s −1 (for heights increasing from 1.05 to 1.15 R ).As can be also seen in Fig. C.2, there is a wave reflection, which is best observed at a height of 1.15 R .In the southern direction, the velocities are smaller, increasing from 220-300 km s −1 with height, and the EUV wave propagation is seen only to about 350 Mm.This smaller extent of the wave propagation is most probably related to the southern coronal hole, which represents an obstacle to the wave propagation due to its high Alfvén speed (e.g., Piantschitsch et al. 2018;Downs et al. 2021).In the stack plots, one can also see that close to the eruption center, it is difficult to disentangle the CME flanks and the wave, but thereafter it is well seen that the CME flank expansion stops whereas the wave decouples from it and further propagates as a freely propagating fast MHD wave (e.g., Warmuth 2015).Such behavior has been observed in previous events and led to the interpretation that the EUV wave is initiated by the fast initial lateral expansion of the CME flanks (Veronig et al. 2008;Veronig et al. 2018;Kienreich et al. 2009;Patsourakos & Vourlidas 2009).

Fig. 1 :
Fig. 1: Longitudinal spacecraft constellation and magnetic connectivity at 16:00 UT on 17 April 2021 (left) together with multi-spacecraft SEP observations (right).The upper panel shows ∼1 MeV electron intensities and the lower panel ∼25 MeV proton intensities (respectively ions) observed by the spacecraft indicated by the same colors in the left figure.The observers' configuration plot (left panel) has been produced using the Solar MAgnetic Connection HAUS (Solar-MACH; Gieseler et al. 2022) tool.
Fig.4: X-ray observations of the associated flare.Top: Soft X-ray fluxes as recorded by GOES-16.Bottom: Normalized background-subtracted STIX count rates integrated over two different energy bands.Note the gradual evolution of the thermal emission at 4-10 keV (blue) as opposed to the multiple non-thermal spikes seen at 15-25 keV (green) and at 25-50 keV (red; multiplied by 0.3 for clarity).For both GOES and STIX, times have been shifted so that they refer to UT at the Sun.
b) Lat. (b) V obs Long. (b) Lat. (b) Polarity Long. (b,c)Lat.(b,c) Notes.(a)PFSS and ADAPT-WSA footpoints at 1 R ;(b) Longitude and latitude values are given in the Carrington coordinate system; (c) Average and standard deviation values calculated from all 12 realizations of the ADAPT-WSA.(d)T his is not an observational value but a nominal value consistent with ENLIL simulations 17-Apr-2021 15:58:23-16:01

Fig. 5 :
Fig.5: Flare evolution as seen in a series of STIX HXR images overlaid on STEREO A EUVI 304 Å images that have been rotated so that they correspond to the viewpoint of Solar Orbiter.Depicted are the coronal thermal source (red contours) and the chromospheric non-thermal footpoints (blue contours) reconstructed with the Expectation Maximization algorithm.The integration times (UT at the Sun) correspond to the eight non-thermal HXR peaks with the highest number of counts above 15 keV.Additionally, the observation times of the EUVI images are shown.For reference, a longitude-latitude grid (in Stonyhurst coordinates) with a spacing of 5 • is overplotted.

Fig. 6 :
Fig. 6: EUV observations by SDO/AIA (left) and STEREO A/EUVI (right) at the same instant of time.The green mesh corresponding to the 3D reconstruction of the CME is overlaid to base-difference images shown in the upper panels.

Fig. 7 :
Fig. 7: Base-difference images of the coronagraph observations by SOHO/LASCO/C2 and STEREO A/COR1 (left) and LASCO/C3 and STEREO A/COR2 (right) at different times.The green (red) mesh corresponding to the 3D reconstruction of the CME (CME-driven shock) is shown in the lower panels.The white arrows indicate the signatures of the CME-driven shock.

Fig. 8 :
Fig.8: EUV wave overview as observed in STEREO A/EUVI 195 Å running-difference images from 16:15 to 16:40 UT.We follow the EUV wave in four angular sectors 1-4.A movie accompanying the figure is available online (movie1).Markers show magnetic footpoints derived for STEREO A (red), BepiColombo (yellow), Solar Orbiter (blue), and ACE (green) spacecraft.The magnetic footpoints are determined using a combination of ballistic backmapping in the heliosphere and backmapping below the source surface using ADAPT-WSA to 1 Rs (points 1), a standard PFSS model to 1 R (points 2), and to a height of 100 Mm above the photosphere (points 3).As seen in the accompanying video, the EUV wave reaches the BepiColombo footpoints at around 16:55 UT (point 3) and at around 17:00 UT (points 1 and 2).Times refer to the observation time at STEREO A.

Fig. 9 :
Fig. 9: Decametric type II radio burst observed by the Swiss-Landschlacht e-Callisto receiver.Fundamental and second harmonic lanes are marked by F and H, respectively.
panel 1) to describe the spectral morphology of the TIII and TII bursts.The identification of the different TIII and TII bursts in Parker Solar Probe observations provide the foundation for the multi-spacecraft directivity analysis presented in Appendix B, where we combine the identification of the bursts with the cross-calibrated data from other spacecraft to locate the source in interplanetary space (Fig. B.1).

Fig. 10 :
Fig. 10: Radio spectrograms from all available space-borne observatories.Panels 1 & 2 show Stokes I and the net polarization (Stokes V/I) from Parker Solar Probe.Panels 3 -5 show the Stokes I measurements from Solar Orbiter, STEREO A, and Wind, respectively.The different bursts are indicated in panel 1.The TII bursts are marked in other panels by rectangular boxes of red (TII(1)), yellow (TII(2)), and orange (TII(HB)).

Fig. 11 :
Fig. 11: Radial velocity contour plot from the ENLIL simulation in the ecliptic plane.The black and white dashed lines represent the IMF lines, and the black contours track the ICME ejecta.The white lines correspond to the HCS, which separates the regions with opposite magnetic polarity, shown in blue (negative) or red (positive) on the outer edge of the simulation region.Left panel: magnetic connectivity of the different spacecraft around the particle solar release time.Right panel: SEP event-related ICME arrival to Solar Orbiter.

Fig. 12 :
Fig. 12: In-situ SEP time profiles as well as plasma and magnetic field observations by BepiColombo (left) and Parker Solar Probe (right).Top: Energetic electron and proton temporal profiles observed from several energy channels and instruments.For SIXS, we use fluxes detected in side 2 of the detector.The flare eruption time is represented by the arrow on the upper x-axes.The vertical solid line and gray shaded area, respectively, indicate IP shock and ejecta transit observed by BepiColombo.Bottom:In-situ plasma and magnetic field observations.The panels present, from top to bottom, the magnetic field magnitude, the magnetic field latitudinal and azimuthal angles, θ B-RTN and φ B-RTN , the solar wind speed, and the proton density, where RTN stands for radial-tangential-normal coordinates (IP structures as described in top panel).The pink lines represent the ENLIL simulation results.

Figure 14 (
Figure 14 (left) shows the SEP event observed at near-Earth spacecraft.Similar to STEREO A, Earth is embedded in the trailing portion of an SIR, after the stream interface

Figure 16 (
Figure16 (left)  shows the energetic electron observations at 920 keV (top panel) and 90 keV (third panel) in the different viewing directions as provided by EPI-Hi/LET and in two wedges of EPI-Lo, respectively.The determined onset times using the sunward-looking sectors and 5-min averaged data are marked by the red dashed lines.

Fig. 15 :
Fig. 15: Pitch-angle distribution of 106 keV electrons (left) and 8.02 MeV protons (right) measured by BepiColombo/SIXS.Top: pitch-angle coverage of sides 0-3, middle: intensities measured in sides 0-3, bottom: pitch-angle distribution with color-coded intensities normalized to the median of each time step.Gray pitch-angle bins mark no pitch-angle coverage, while white bins are zero-count periods.

Fig. 17 :
Fig.17: VDA of protons from PSP EPI-Hi/LETA from 1 to 30 MeV.Top panel: the red line is a 'by-eye' fit to the onset of the observed intensities as a function of 1/v and time.Bottom panel: the same data are plotted, but with the velocity dispersion removed.The legend provides the derived path length and injection times corresponding to the fit line.

Fig. 18 :
Fig. 18: Pitch-angle distribution of 86-130 keV electrons (left) and 7.4-9.2MeV protons (right) observed by Solar Orbiter/EPD-EPT and EPD-HET, respectively.Top: Pitch-angle coverage of the four different sensor apertures.Middle: Intensities observed by each field of view.Bottom: Pitch-angle distribution with color-coded intensities.

Fig. 19 :
Fig. 19: VDA of protons measured by Solar Orbiter/HET sun (red points) and EPT sun (blue points, not included in the VDA fit).The vertical red line and shade represent the derived injection time and uncertainty.

Fig. 20 :
Fig. 20: Pitch-angle distribution of 85-125 keV electrons (left) and 4-6 MeV protons (right) observed by STEREO A/SEPT and LET, respectively.Left plot shows from top to bottom: Pitch-angle coverage of the different sensor apertures, intensities observed by each field of view, and pitch-angle distribution with color-coded intensities.Right plot shows from top to bottom: Pitch-angle distribution with color-coded intensities, pitch-angle pointing of the different LET sectors, intensities observed by each sector, and average intensities measured by the eight sectors on each side of the instrument.

Fig. 21 :
Fig. 21: Inferred SEP injection times (vertical lines with temporal error bars on top) overplotted on the radio spectrogram as observed by PSP/RFS and the 15-25 keV X-ray observations by Solar Orbiter/STIX.All times have been shifted to the Sun by assuming the propagation time of the emission to the respective spacecraft.The shaded ranges mark the times when the spacecraft establish a magnetic connection with the CME-driven shock including an uncertainty of ±3 min.
the questions related to particle propagation in interplanetary space during the meeting of the team No. 459 (led by G. Li and L. Wang).B.S.-C.acknowledges support through UK-STFC Ernest Rutherford Fellowship ST/V004115/1 and STFC grant ST/V000209/1.The work of F.S. was supported by DLR grant No. 50 OT 1904.N.W. acknowledges support from NASA program NNH17ZDA001N-LWS and from the Research Foundation -Flanders (FWO-Vlaanderen, fellowship No. 1184319N).C.O.L. acknowledges support from NASA's LWS Program (Grant No. 80NSSC21K1325) and the MAVEN project funded through the NASA Mars Exploration Program.C.O.L. and C.M.S.C. acknowledge support from the IMPACT Investigation by the NASA Heliophysics Division through the STEREO Project Office at NASA GSFC (Grant No. 80NSSC18K1446).M.L. acknowledges support from the Italian Space Agency and the National Institute of Astrophysics through the ASI-INAF n. 2020-35-HH.0 agreement for the development of the ASPIS prototype of scientific data centre for Space Weather.ENLIL simulation results have been provided by the CCMC at NASA Goddard Space Flight Center (GSFC) through their public Runs on Request system (http://ccmc.gsfc.nasa.gov;run ID Laura_Rodriguez-Garcia_041322_SH_1).The WSA model was developed by N. Arge, currently at GSFC, and the ENLIL Model was developed by D. Odstrcil, currently at George Mason University.

Fig
Fig. C.1: Kinematics of the EUV wave observed on the solar disk.Left panels: Perturbation profiles demonstrating the EUV wave propagation in the four angular sectors considered.Big red markers indicate the peak amplitude of each profile, while small orange markers show the estimated wavefront locations extracted at 30% of the peak value, which we use for the assessment of the EUV wave kinematics presented in the right panels.A movie accompanying the figure is available online (movie2) and represents the dynamics of the intensity perturbation profiles.Right panels: Distance-time profiles of the EUV wave fronts for the four angular sectors.Dots indicate data points, and the solid lines show the linear fit.
Fig. C.2: Stack plots derived from STEREO A/EUVI 195 Å base-difference images in vertical slits at heights from 1.05 to 1.30 R above the solar surface from the Sun center.The EUV wave fronts are shown by black lines.

Table 2 :
Separation angles between location of the flare and spacecraft magnetic footpoints based on ADAPT-WSA values Spacecraft Lon.sep.Lat.sep.Total sep.

Table 3 :
First intersection between the coronal shock and magnetic field lines connecting to the spacecraft as determined with the ADAPT-WSA model (Sect.3).Times refer to observation times at 1 au.

Table 4 :
Timing of solar phenomena and inferred SEP injection times t inj .All times shifted to the Sun.In contrast, for Solar Orbiter a sole shock source could be justified as the shock connection time happens during the first inferred injection time (taking into account the large error bar), which is the one of the MeV electrons, and well before the inferred injections of the ∼100 keV electrons and that of the protons.