- Institute for Physical Science and Technology, University of Maryland, College Park, MD, United States
The earliest evidence of spatial distributions of solar energetic particles (SEPs) compared events from many different source longitudes on the Sun, but the early Pioneers provided the first evidence of the large areas of equal SEP intensities across the magnetically confined “reservoirs” late in the events. More detailed measurements of the importance of self-generated waves and trapping structures around the shock waves that accelerate SEPs were obtained from the Helios mission plus IMP 8, especially during the year when the two Voyager spacecraft also happened by. The extent of the dozen widest SEP events in a solar cycle, which effectively wrap around the Sun, was revealed by the widely separated STEREO spacecraft with three-point intensities fit to Gaussians. Element abundances of the broadest SEP events favor average coronal element abundances with little evidence of heavy-element-enhanced “impulsive suprathermal” ions that often dominate the seed population of the shocks, even in extremely energetic local events. However, it is hard to define a distribution with two or three points. Advancing the physics of SEPs may require a return to the closer spacing of the Helios era with coverage mapped by a half-dozen spacecraft to help disentangle the distribution of the SEPs from the underlying structure of the magnetic field and the accelerating shock.
1 Introduction
The spatial distribution of solar energetic particles (SEPs), and their variation with time, particle species, and energy, is fundamental to an understanding of the physics of particle acceleration and transport in the heliosphere. How much of the variation we see at a single spacecraft is a true time variation and how much is spatial variation being convected past? Is the SEP source itself broadly extended in space or do SEPs somehow diffuse out of a limited source? Does a shock source sample different abundances of seed ions from different places or, for example, does Fe simply scatter less than O, to produce early enhancements and later suppressions in Fe/O? Multispacecraft comparisons can be a key to distinguishing the physical effects dependent upon space and time.
1.1 SEP history and context
Multispacecraft measurements and the perceived need for them generally followed the study of SEPs on or near Earth. The SEP events observed first (Forbush, 1946) were rare, large, energetic “ground-level enhancements” (GLEs), where GeV protons produce a nuclear cascade through the atmosphere to the ground level that enhances the continuous signal produced similarly by the galactic cosmic rays (GCRs). Since solar flares were found to accompany these early SEP events, flares were considered a possible source. However, the spatial span of these flares, stretching from the east on the Sun to behind the western limb, raised a significant problem: how could the SEPs cross magnetic field lines radiating out from the Sun to find their way to Earth?
Meanwhile, solar radio astronomers, also using ground-based instruments, had identified different sources triggered by energetic solar electrons. Radio emission, excited at the local plasma frequency, depends upon the square root of the local electron density, which decreases with distance from the Sun. Wild et al. (1963) described radio type III bursts where frequencies decreased rapidly, excited by 10–100 keV electrons that streamed out from a source near the Sun. There are also type II bursts where the source moved at the slower speed of a ∼1,000 km s−1 shock wave. Wild et al. (1963) suggested two types of SEP sources, namely, point sources of mostly electrons near the Sun and fast shock waves that could accelerate energetic protons, like those that produce GLEs. Even though shock acceleration was well-known in other contexts, like the supernova sources of GCRs, Wild et al. (1963) were far ahead of their time in solar physics. After 20 years, the clear 96% association of large SEP events with shocks driven by fast, wide coronal mass ejections (CMEs) was established by Kahler et al. (1984). Yet, 10 years later, Gosling 1993) and Gosling (1994) still needed to point out the error of the “solar flare myth”.
Parker (1965) explained particle transport in terms of pitch-angle scattering as SEPs followed magnetic field lines out from the Sun. Diffusion theory is an important tool when there is actually a physical mechanism, like pitch-angle scattering, to produce the random walk. There is also a random walk of the magnetic field footpoints (Jokipii and Parker, 1969; Li and Bian, 2023) prior to events, which can produce an effective random walk perpendicular to the mean magnetic field, but it is independent of SEP-event parameters and is completely inadequate to explain a huge spread of SEPs far from a presumed source longitude near a flare. Other schemes such as the “birdcage” model (Newkirk and Wenzel, 1978) were also invented to spread SEPs from a flare source (see review, Sec 2.3 in Reames, 2021a). Reinhard and Wibberenz (1974) envisioned a mysterious “fast propagation region” extending 60⁰ from the flare to spread the SEPs prior to their slower interplanetary journey. Could this region actually match the surface of a shock? After decades of resistance, a flare source has mainly been abandoned for the largest SEP events that are now generally attributed to spatially extensive CME-driven shock waves (Mason et al., 1984; Gosling, 1993; Reames, 1995b; 1999; 2013; 2021b; Zank et al., 2000; 2007; Lee et al., 2012; Desai and Giacalone, 2016), especially for GLEs (Tylka and Dietrich, 2009; Gopalswamy et al., 2012; 2013; Mewaldt et al., 2012; Raukunen et al., 2018). Observations (e.g., Kahler et al., 1984) did replace the “fast propagation region” with the surface of a CME-driven shock wave that actually accelerates the particles, beginning at 2–3 solar radii (Tylka et al., 2003; Cliver et al., 2004; Reames, 2009a; Reames, 2009b) and continuing far out into the heliosphere. We will soon study from STEREO observations (e.g., Figure 5) that shocks easily wrap around the Sun, expanding widely across magnetic field lines where SEPs alone cannot pass.
As the element and isotope abundances in SEPs began to be measured, they would present new evidence for two different physical sources of SEPs. The earliest measurements, during large SEP events with nuclear emulsions on sounding rockets, extended element abundances up to S (Fichtel and Guss, 1961) and then to Fe (Bertsch et al., 1969). Later studies showed that average SEP element abundances in large events were a measure of coronal abundances that differed from photospheric abundances as a simple function of the elements’ first ionization potential (FIP; Meyer, 1985; Reames, 1995; Reames, 2014; Reames, 2021a; Reames, 2021b). The FIP dependence of SEPs differs fundamentally from that of the solar wind (Reames, 2018a; Laming et al., 2019), probing the physics of the formation of the corona itself.
However, early measurements in space soon identified a completely new type of event, distinguished by extremely high abundances of 3He in some events, such as 3He/4He = 1.52 ± 0.1 (e.g., Serlemitsos and Balasubrahmanyan, 1975), vs. a solar value of ∼5 × 10−4. Production of 3He from 3He by nuclear fragmentation was ruled out by the lack of 2H and Li, Be, and B fragments from C and O, e.g., Be/O and B/O < 2 × 10−4 (McGuire et al., 1979; Cook et al., 1984). The huge enhancements of 3He were produced by new physics, involving a wave–particle resonance (e.g., Fisk, 1978; Temerin and Roth, 1993) with complex spectra (Liu et al., 2004; Liu et al., 2006; Mason, 2007). The 3He-rich events were associated with the beamed non-relativistic electrons (Reames et al., 1985) and their type III radio bursts (Reames and Stone, 1986) studied by Wild et al. (1965). Element abundances of these “impulsive” 3He-rich events were also different from those of the large “gradual” shock-associated events (Mason et al., 1986; Reames, 1988); here, the enhancements increased as a power of the ion mass-to-charge ratio A/Q (Reames et al., 1994), which was especially clear when it became possible to measure elements above Fe with element groups resolved as high as Au and Pb (Reames, 2000; Mason et al., 2004; Reames and Ng, 2004). Average enhancements varied as (A/Q)3.6 with the Q value determined at a temperature of ∼3 MK (Reames et al., 2014a; Reames et al., 2014b). This power law can be used in a best-fit method to determine source plasma temperatures (Reames, 2016; Reames, 2018b). Impulsive SEP events have been associated with magnetic reconnection (Drake et al., 2009) on open field lines in solar jets (Bučík, 2020), which also eject CMEs that sometimes drive shocks (Kahler et al., 2001; Nitta et al., 2006; Nitta et al., 2015; Wang et al., 2006; Bučík et al., 2018a; Bučík et al., 2018b; Bučík et al., 2021) fast enough to reaccelerate the enhanced 3He and heavy ions along with ambient protons and ions (Reames, 2019; Reames, 2022b). Derived abundances from γ-ray lines suggest that flares involve similar physics (Murphy et al., 1991; Mandzhavidze et al., 1999; Murphy et al., 2016), but those accelerated particles are trapped on loops, losing their energy to γ-rays, electron bremstrahlung, and hot, bright plasma. The opening magnetic reconnections that drive jets also must close neighboring fields to form flares (see, e.g., Figure 4 in Reames, 2021b).
Shock waves can also accelerate residual suprathermal impulsive ions that can pool to provide a seed population (Mason et al., 1999; Tylka et al., 2001; Desai et al., 2003; Tylka et al., 2005; Tylka and Lee, 2006). The combination of two fundamental physical acceleration mechanisms, magnetic reconnection, and shock acceleration, as well as two distinctive element abundance patterns, led Reames (2020) and Reames (2022b) to suggest four distinguishable SEP event abundance pathways:
SEP1: “Pure” impulsive magnetic reconnection in solar jets with no fast shock.
SEP2: Jets with fast, narrow CMEs driving shocks that reaccelerate SEP1 ions plus ambient coronal plasma. Pre-enhanced SEP1 ions dominate high Z, and ambient protons dominate low Z.
SEP3: Fast, wide CME-driven shocks accelerate SEP1 residues from active-region pools from many jets, plus ambient plasma. Again, the SEP1 seed ions dominate high Z.
SEP4: Fast, wide CME-driven shocks accelerate ions where ambient plasma completely dominates.
Persistent pools of SEP1-residual seed ions available for reacceleration in SEP3 events have now been widely observed and reported (Richardson et al., 1990; Desai et al., 2003; Bučík et al., 2014; Bučík et al., 2015; Chen et al., 2015; Reames, 2022a), and may be fed by numerous impulsive events too small to be distinguished individually.
Large gradual SEP events are frequently accompanied initially by type III bursts. In principle, these impulsive jets could inject a SEP1 seed population. However, any such injection of SEP1 ions seems to be swamped by coronal seed ions in the many large SEP4 events. The extremely abundant type III electrons can be distinguished from shock-accelerated electrons by their proximity to the event source; the latter emerge only in poorly connected events (Cliver and Ling, 2007).
The spatial extent of these features and the underlying physics are of considerable interest. Can we find spatial differences in the seed populations sampled by shocks? Unfortunately, each point on a shock moving radially across Parker spiral fields will have spread particles over 50⁰—60⁰ by the time it reaches 1 AU (Reames, 2022a). This tends to blur any initial spatial variations in abundances.
2 Near-earth observations
Before we could monitor a single SEP event at multiple locations, we could study multiple events at Earth from different source longitudes on the Sun. Using data from IMP 4, 5, 7, 8, and ISEE 3 over nearly 20 years, Cane et al. (1988) ordered the 20-MeV proton profiles of 235 large SEP events as a function of the solar source longitude. This study gave the typical extent and time profiles of an “average” SEP event, sliced as a function of longitude. Recently, Reames (2023) revisited this study, adding examples that better illustrate the evolving role of the shock source as it propagates outward from the corona. A collection of sample time distributions of protons vs. longitude from this work is shown in Figure 1. Although the SEP events shown in Figure 1 are all different, we expect that slices of a single event at different longitudes could show similar features.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1254266/fspas-10-1254266-HTML/image_m/fspas-10-1254266-g001.jpg)
FIGURE 1. IMP 8/GME proton intensities of the listed energies are shown vs. time in each panels (a–g) around a central map showing their nominal distribution around a fixed CME-driven shock source. Event onset times in each panel are flagged with the source longitude and the CME speed (when available), and the times of shock passage are noted with the shock speed and θBn, the angle between the field B and the shock normal (when available). All panels have the same intensity scale. Dashed lines in the central map illustrate that a shock must first encounter any field line at its footpoint on the east flank but may first strike it far from the Sun on the west (Reames, 2023).
Historically, the intensity peaks at the shock seen on several of the profiles in Figure 1 were called “energetic storm particle” (ESP) events. During their acceleration, particles are trapped near the shock by self-generated Alfvén waves (Bell, 1978a; Bell, 1978b; Lee, 1983; Lee, 2005), creating an autonomous structure that can propagate onto new field lines where earlier accelerated particles may be absent, as shown in Figure 1E (Reames, 2023).
3 Multispacecraft observations
3.1 Pioneer
Some of the early Pioneer spacecraft were launched into Earth-like solar orbits, although coverage was only hours per day. When there was no spacecraft near Earth, this led to the awkward comparison of 15-MeV proton data from Pioneer 6 and 7 with GeV ground-level neutron monitor measurements at Earth (Bukata et al., 1969). Fortunately, McKibben (1972) was able to include data from IMP 4 near Earth. Although most of these spatial distributions were analyzed in terms of adjustable coefficients in the fashionable perpendicular diffusion, rather than shock acceleration, McKibben also noted that late in SEP events, proton intensities could be identical over large spans of longitude. Later named “reservoirs” by Roelof et al. (1992), these regions involved particles quasi-trapped magnetically behind the shock; as the volume of this magnetic bottle expands, adiabatic deceleration (e.g., Kecskeméty et al., 2009) decreases all intensities, preserving spectral shapes (Reames et al., 1997; Reames, 2013). Extreme SEP scattering, once used to explain this slow decay, is actually found to be negligible in reservoirs since particles from new 3He-rich events travel scatter-free across them (Mason et al., 1989; Reames, 2021a).
3.2 Ulysses
Whenever there was a dependable spacecraft stationed near Earth, like IMP 8, any traveling spacecraft, such as the solar-polar-orbiting Ulysses, allowed two-spacecraft comparisons. The reservoir comparisons by Roelof et al. (1992) found uniform intensities behind the shock, extending radially over 2.5 AU from IMP 8 to Ulysses. Ulysses observed reservoirs at heliolatitudes up to >70⁰, N and S (Lario, 2010), and in other electron observations (Daibog et al., 2003).
3.3 Helios, IMP, and Voyager
The two Helios spacecraft followed neighboring solar orbits from 0.3 to 1.0 AU, beginning in 1974. Beeck et al. (1987) used data from these spacecraft to study the radial and energy dependence of diffusive scattering of protons, while, in a larger study, Lario et al. (2006) fit the peak intensities and the fluence of events to the form R-n exp[- k(ϕ - ϕ0)2], where n and k are constants, R is the observer’s radius in AU, ϕ is the solar longitude of the observer, and ϕ0 is the solar longitude of the event centroid. For a Gaussian distribution, k = (2σ2)−1. Lario et al. (2006, 2007) also considered the time variation of the point where the observer’s field line intercepts the expanding shock source, an important feature of shock acceleration defined by Heras et al. (1995).
Helios provided an excellent opportunity to study spatial distributions of shock-accelerated particles and their reservoirs, especially during 1978 when the Voyager spacecraft were nearby (Reames et al., 1996; Reames et al., 1997; Reames et al., 2013; Reames, 2023). An especially interesting event is shown in Figure 2 where the intensities of SEPs at Voyager 2 (Figures 2A or 2E) do not begin to increase until after the shock passes S1 (Figure 2B), presumably because this is where it first intercepts the field line to Voyager. Proton intensities then slowly rise as Voyager becomes connected to stronger and stronger regions of the approaching shock until the ESP structure finally arrives at S4 on 6 January. It should be noticed in Figure 2A that the peak intensities near the shock are similar at all four spacecraft. The ESP structure forms as protons streaming away from the shock generate resonate waves of wave number k ≈ B/μP, where B is the field strength, P is the proton rigidity, and μ is the cosine of its pitch angle (Lee, 1983; Lee, 2005; Ng and Reames, 2008; Reames, 2023). These waves trap particles in ESP structures. As the structure moves out to lower B, the resonance shifts so that high energies preferentially leak away early, as also seen in Figure 1E. Voyager sees the pure “naked” ESP event with few of the streaming protons that created it. The shock that eventually arrives at Voyager generated the intense early streaming protons seen early by Helios 1 and an intermediate structure as it passed Helios 2 and IMP 8.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1254266/fspas-10-1254266-HTML/image_m/fspas-10-1254266-g002.jpg)
FIGURE 2. In (A), intensities of 6–11-MeV protons are compared for Helios 1, Helios 2, IMP 8, and Voyager 2 during the 1 January 1978 SEP event, while (B) shows the spatial configuration of the spacecraft on their initial field lines, and stages in the expansion of the CME-driven shock at S1, S2, etc., are sketched. Intensity–time profiles for a full list of energy intervals are shown for (C) Helios 1, (D) IMP 8, and (E) Voyager 2. MC is a magnetic cloud from the original CME (Burlaga et al., 1981). The onset time of the event is flagged by E6 and the shock passage at each spacecraft by H1, H2, I8, and V2 (Reames, 2023).
The spacecraft distribution for the event shown in Figure 2 is unique; in this spacecraft sample, the SEPs are produced as the source shock moves radially. The width of this CME is limited, but the spacecraft are positioned to follow the evolution of the event: 1) well-connected Helios 1 samples its central production near the Sun, 2) Helios 2 and IMP 8 sample its production near 1 AU, and 3) Voyager 2 scans its western flank and then samples its central strength as it passes 2 AU. A fortuitous occurrence occurred: NASA did not design Voyager as a complement to Helios.
Another fortuitous observation with these spacecraft occurred in September 1978 and is shown in Figure 3. In Figure 3A, well-connected IMP 8 shows a fast increase in intensities and a peak near the time of shock passage, while Helios 1 and 2, far around the west flank, increase more slowly and then join IMP 8 in a reservoir behind the shock on 25 September. Distant Voyagers show a slow SEP increase to a plateau near the time S2, where IMP 8 joins the same intensity once it is no longer constrained by the east flank of the shock. Then, at S3, the western flank of the shock strikes field lines that send particles sunward to IMP 8 and outward to Voyager 2, and then to Voyager 1. This second SEP peak is clearly seen at energies above 40 MeV in Figure 3C and probably even in Figure 3D. A “left” and then a “right” from the same shock are observed; at Voyager, the two peaks are comparable in size.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1254266/fspas-10-1254266-HTML/image_m/fspas-10-1254266-g003.jpg)
FIGURE 3. In (A), intensities of 3–6-MeV protons are compared for IMP 8 (blue), Helios 1 (green), Helios 2 (yellow), Voyager 1 (red), and Voyager 2 (violet) in the 23 September 1978 SEP event, while (B) maps the configuration of the spacecraft on their initial field lines and shows the expansion of a CME-driven shock at S1, S2, and S3. The onset time of the event is flagged as W50, and the shock passage at each spacecraft as H1, H2, I8, V1, and V2. (B) Western flank of shock S3 intercepts the blue and red fields, where arrows direct particles accelerated sunward to IMP 8 and then outward to Voyager 2, respectively, and later to Voyager 1. Intensity–time profiles for a full list of energy intervals are shown for (C) IMP 8 and (D) Voyager 1 (Reames, 2023).
It should be noted that the second peak in the IMP 8 data in Figure 3C shows significant velocity dispersion corresponding to the ∼6-AU path inward from the new source (actually a larger delay than in the first peak), while none is seen at Voyager 2 (Figure 3D) near this source. Furthermore, the intensities at IMP 8 and Voyager are quite similar in the second peak; any new injection from the Sun would have produced a huge difference in intensities like that seen in the first peak because of the great difference in radial distances since IMP 8 and the Voyagers seem to be on similar field lines. At the second peak, the shock has filled these field lines, forming a reservoir with similar intensities at IMP and Voyager.
Another parameter that can depend upon solar longitude is the solar particle release (SPR) time derived from velocity dispersion (Tylka et al., 2003; Reames, 2009a; Reames, 2009b). If the first particles of each energy to arrive at the spacecraft have been released at nearly a single time, the SPR time, and have scattered little, their travel time will be dt = L/v, where L is the field line length from the source to observer and v is the particle velocity. Plotting the onset times vs. v−1 will yield a linear fit with slope L and intercept at the SPR time. Figure 4B shows such a plot from Reames and Lal (2010) for the spacecraft distributed, as shown in Figure 4A, during the GLE of 22 November 1977. The parabolic fit for the height shown in Figure 4D, not from this event, is a fit of 26 GLEs observed from Earth (Reames, 2009b); presumably, it is a first-order correction for weaker, more slowly evolving, shock flanks. Gopalswamy et al. (2013) fit the SPR height for GLEs at Earth directly to the source longitude (uncorrected for foot-point motion). Type II bursts begin near ≈1.3 RS, but they certainly need not correspond to the same field line longitude or shock physics as the measured SEPs. Non-relativistic electrons that produce type II bursts do not resonate with Alfvén waves like ions and hence are more likely to be accelerated by quasi-perpendicular regions of the shocks.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1254266/fspas-10-1254266-HTML/image_m/fspas-10-1254266-g004.jpg)
FIGURE 4. (A) Distribution of spacecraft during the 22 November 1977 SEP event at W40; (B) plots and least-square fits of the particle onset times vs. v−1 at available energies for several spacecraft; (C) time delay; (D) corresponding shock height of the SPR time vs. solar source longitude minus the footpoint longitude. Data from Helios 1 are limited (Reames and Lal, 2010). The parabola shown in (D), from Reames (2009b), may just be a first-order correction for the occurrence of the shock on its flanks.
Thus, Helios allowed spatial comparisons of the beginnings of events, their evolution and ESP events at the middle, and their reservoirs at the end (Reames et al., 1997), noted previously.
3.4 STEREO, with Wind, ACE, and SOHO
The launch of Wind in November 1994 began a new era in SEP coverage from Earth, and it was later joined by SOHO and ACE, adding different capabilities. This coverage still exists in 2023. STEREO ahead (A) and behind (B) were launched together in late 2006 at the beginning of solar minimum. By September 2012, they had reached the ±120⁰ longitude, providing three equally spaced observing points around the Sun and optimal coverage for finding the most extensive SEP events.
With STEREO, it has become even clearer that SEP events can be quite extensive. For >25-MeV proton events, Richardson et al. (2014) found that 17% spanned three spacecraft; 34%, two spacecraft; and 36%, one spacecraft, with 13% unclear. Studying the abundances of H, He, O, and Fe at 0.3, 1, and 10 MeV amu−1, Cohen et al. (2017) found only 10 three-spacecraft events, out of 41.
Most of the studies of particle distributions with STEREO have involved Gaussian fits of the intensity peaks or fluences at three longitudes. Three points determine a parabola, and a Gaussian is a parabola in logarithmic space, so Gaussians tend to fit the data very well. Xie et al. (2019) studied 19–30-MeV protons in 28 events finding σ = 39⁰ ± 6.8⁰. Paassilta et al. (2018) compiled a list of 46 wide-longitude events above 55 MeV, of which seven were suitable, averaging σ = 43.6 ± 8.3⁰, and for 14 events with E > 80 MeV, de Nolfo et al. (2019) found an average σ ≈ 41⁰. Lario et al. (2013) studied the 15–40 and 25–53-MeV proton peak intensities and found σ ≈ 45 ± 2⁰ for both proton energy ranges. Cohen et al. (2017) found average three-spacecraft events centered at ϕθ ≈22 ± 4° west of the flare site and ≈43° ± 1° wide; they found no dependence of the width on the charge-to-mass ratio Q/A of the elements. This lack of dependence on A/Q seems to argue against lateral diffusive transport. Kahler et al. (2023) fit the data on 20-MeV protons to hourly Gaussians, showing substantial broadening of the distribution with time, especially on the western flank where the shock expanded across new spiral field lines. However, all of these similar widths and parameters only apply to the largest ≈20% of events that span three spacecraft. We have no widths from the one-spacecraft third of the events. For electron events, Klassen et al. (2016) examined events for closely spaced (<72⁰) STEREO spacecraft and found they were not well-fit by Gaussians, while Dresing et al. (2018) pointed out that peak intensities used for the Gaussian fits may not represent the real spatio-temporal intensity distribution as the intensity peaks may have been measured at different times.
Some of the space–time coupling in SEP distributions is illustrated by the event shown in Figure 5. Here, the shock itself is seen at each spacecraft, yet the slight onset delay at Wind suggests it misses the base of that field line. Both STEREO spacecraft show fast intensity increases, followed by early declines, but the intensity at Wind peaks well after the shock passage. How do we distinguish variations in space and time? When does a Gaussian spatial distribution apply? Is the Gaussian formed by the peaks, or by the fluences; is it spatial or also temporal? Even hourly Gaussians miss the true behavior late in an event.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1254266/fspas-10-1254266-HTML/image_m/fspas-10-1254266-g005.jpg)
FIGURE 5. Schematic representation of the time variations in ∼20-MeV protons in the large three-spacecraft SEP event of 25 February 2014 on the left and the map of spacecraft configuration on the right. Vertical time flags mark the event onset from S12E77 and the shock passage times at each spacecraft. In the map, circles follow expansion of a spherical shock wave, and initial Parker spirals connect the Sun with each spacecraft; these field lines would be distorted as the shock passes. Pink field lines measure solar rotation (see text).
The two pink field lines labeled 1 and 2 in the map shown in Figure 5 will be carried, from their initial positions shown, to the red location of Wind by the time the shock reaches 1 or 2 solar radii, respectively, because of solar rotation. Hence, the time profile we see at Wind is greatly enhanced behind the shock by SEPs swept in on field lines that were much more centrally located initially. Meanwhile, STEREO A is rapidly acquiring field lines with decreased intensities of SEPs that have rotated from behind the Sun, and STEREO B eventually joins it in a suppressed reservoir region late on 2 March. Can SEP models follow this mixture of space–time variations? Time variations longer than about a day need to accommodate the 13⁰ day−1 solar rotation. This event evolves over a week.
It should also be noted in Figure 5 that the accelerating shock wave is seen at all three spacecraft, as is often the case. These shock waves are able to wrap around the Sun, crossing field lines that the particles alone cannot cross. Thus, STEREO shows us the way the physics of shock acceleration can easily replace the early confusion of the “birdcage model,” “coronal diffusion,” and the “fast propagation region”—diffusion from a point source—that once held sway.
3.5 SEP element abundances
Spatial distributions of fluences of element abundances were studied extensively by Cohen et al. (2017) using STEREO and ACE data. The Gaussian distributions of H, He, O, and Fe were found to be similar in 10 three-spacecraft events, suggesting that the widths are independent of rigidity or transport. They found the average three-spacecraft Gaussian distributions to be 43° ± 1° wide, although with significant variations. None of these large three-spacecraft SEP events showed any of the enhancements of heavy elements, e.g., Fe/O, typically found in shock-reaccelerated impulsive suprathermal ions of SEP3 events. An analysis of the A/Q-dependence of the element abundances in a typical large SEP4 event is shown in Figure 6, where power-law fits of flat or suppressed heavy elements extend to include protons (Reames, 2020; Reames, 2022b).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1254266/fspas-10-1254266-HTML/image_m/fspas-10-1254266-g006.jpg)
FIGURE 6. (A) Map of the spacecraft distributions during the SEP event of 23 January 2012; (B) corresponding intensity–time profiles of 20-MeV protons at each spacecraft. The abundance enhancements, relative to SEP-average (coronal) abundances, for elements noted by Z, are shown vs. A/Q for the listed time intervals for (C) STEREO B, (D) Wind, and (E) STEREO A. Best fit power laws vs. A/Q for elements with Z > 2 are shown extended down to protons at A/Q = 1.
For each time period, the observed abundances of Z ≥ 6 ions are divided by the corresponding coronal abundances and fit to power-laws. Most SEP4 events show flat fits, i.e., abundances the same as coronal, or declining power laws, as shown in Figure 6. The declining power laws may result from reduced scattering of heavier ions that allows them to leak more easily from the acceleration region SEP3 events have both enhanced protons and power-law heavy-element enhancements that increase with A/Q, reflecting shock acceleration of seed populations of normal ambient coronal ions and impulsive suprathermal ions with their characteristic high-Z enhancement. SEP3 events can be very large. In solar cycle 23, about half of the GLEs were SEP3 events (Reames, 2022a) and half SEP4 events. In weaker cycle 24, when STEREO was available, we found only one of the two-spacecraft events listed by Cohen et al. (2017) that was an SEP3 event (4 August 2011); this is an SEP3 event at Wind (figure 10 in Reames, 2020) and shows similar enhancements at STEREO A, although with poorer statistics. Why are there so few wide SEP3 events? Is it the weak solar cycle or are SEP3 events inherently narrow, perhaps because the primary pools of seed particles are confined?
3.6 Impulsive SEP events
A distinction of impulsive SEP events is that their longitude spread is much more limited than that of gradual events (Reames, 1999), but the width of that distribution depends upon the sensitivity of the instruments observing them and has increased somewhat with time (e.g., Reames et al., 2014a). A significant factor in this width is that variations in the solar wind speed vary with the longitude of the footpoint of the observer’s field line, but there is also variation due to the random walk of the footpoints of the field lines (Jokipii and Parker, 1969; Li and Bian, 2023) prior to the event.
A search for 3He-rich events between ISEE 3 and Helios (Reames et al., 1991) found several corresponding events. The well-studied event of 17 May 1979 (e.g., Reames et al., 1985) with 3He/4He ≥ 10 was seen with similar 3He enhancement by Helios 1 near 0.3 AU, as shown in Figure 7 as a sharp spike of very short (∼1 h) duration, presumably resulting from reduced scattering. Associated electron trajectories were tracked spatially, using the direction and frequency of the radio type III burst as measured from ISEE 3 (Reames et al., 1991).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1254266/fspas-10-1254266-HTML/image_m/fspas-10-1254266-g007.jpg)
FIGURE 7. Extremely 3He-rich event of 17 May 1979 is shown by the track of its radio type III burst in the upper panel, and time histories of electrons and ions are indicated at Helios 1 and ISEE 3 in the following panels. Indicated energies are in MeV for electrons and protons, and MeV amu−1 for He. Arrows above the central panel mark two times the event sources (Reames et al., 1991).
The time history for the Kiel Instrument on Helios 1 in Figure 7 shows composite He; isotope ratios were tabulated separately using pulse-height data. Despite differences in He energies at the two spacecraft, both show the event at 0550 UT with 3He/4He ≈ 10 and that at 1700 UT with 3He/4He ≈ 1.
In the early STEREO era, before the spacecraft were widely separated, Wiedenbeck et al. (2013) found electrons and ions ±20⁰ ahead and behind Earth, during an event that occurred during a solar quiet time. The ions showed strong intensity gradients. The event was associated with a weak CME. Moreover, Klassen et al. (2015) studied an electron beam event early in the STEREO mission.
The discovery of fast CMEs associated with impulsive SEP events (Kahler et al., 2001) suggested that impulsive SEPs could be spread laterally when nearly radial shocks in these SEP2 events distributed particles across spiral field lines. Particles from SEP1 jets without fast shocks would be expected to follow any open field lines from the jet, presumably less widely distributed. Suggestions of greater spreads for SEP2 impulsive events with shocks, based upon current understanding, have not been explored, particularly because the available spacecraft are too widely separated.
Opportunities for multiple measurements of impulsive SEP events are rare. Perhaps the best opportunity is for measurements of the radial variation in SEP scattering between Parker Solar Probe (PSP) or Solar Orbiter and spacecraft near Earth. Are the time profiles of SEPs at PSP near the Sun similar to those of X-ray profiles of the event? To what extent are the jets that release SEPs part of more-extensive flaring systems? Magnetic reconnection that opens some field lines must also close others (e.g., see Figure 4 in Reames, 2021b), but reconnecting closed field lines with other closed lines would allow no escape.
4 Discussion
STEREO has given us an improved sense of the widths of some of the most extensive SEP events, but are the SEP distributions really characterized by nice smooth Gaussians? Are the SPR times or heights really parabolic? During the Helios era, when the Voyagers happened by, a few events showed us greater complexity (e.g., Figures 2, 3). Such events allow us to explore the underlying physics. STEREO also allows 3D modeling of the CME and its shock (e.g., Rouillard et al., 2011; Rouillard et al., 2012; Kouloumvakos et al., 2019; Zhang et al., 2023), but this is only beginning to be extended to a 3D modeling of the SEP distribution this shock produces. Correlations between SEP peaks and shock properties at the base of the observer’s field line (e.g., Koulomvakos et al., 2019) may have reached their limits. SEP peak intensities and times are controlled, and often limited, by transport and are not determined by any single point on the evolving shock; transport and shock contributions both vary with time in complex ways. No single parameter represents an entire SEP event very well.
Kahler’s (1982) “big flare syndrome,” while expressed for a specific case, should be a general warning that correlations do not imply causality. Increasingly energetic magnetic reconnection events at the Sun can spawn bigger flares, faster CMEs, and larger SEP events. They are all correlated, yet Hα flares do not cause CMEs, or GLEs; ultimately, they are all consequences of magnetic reconnection. Later, Kahler (1992) asked “how did we form such a fundamentally incorrect view of the effects of flares after so much observational and theoretical work?” We need a clearer resolution of the underlying physics that does connect a cause and its effects.
Gopalswamy et al. (2013) sought to understand possible differences between the GLE of 17 May 2012 and six other large SEP events with similar CME speeds that were not GLEs. This event had a smaller flare (M5.1) than any of the others. A large effect was produced by the small latitudinal distance of the shock-nose from the ecliptic, i.e., the GLE was better connected to Earth. The GeV protons in this GLE were all produced within ∼8 min. This suggests that a highly localized region of production might occur soon after shock formation (Ng and Reames, 2008) and may differ from the global properties of the CME, discounting correlations. Incidentally, the 17 May 2012 event was the only GLE included in the study by Cohen et al. (2017), and it had very few heavy ions to be studied at three spacecraft; it was not one of the 10 three-spacecraft SEP events.
It has also been observed that GLEs are more likely when shocks pass through solar streamers where the higher densities and lower Alfvén speeds produce higher Alfvénic Mach numbers (Liu et al., 2023) and regions of higher θBn (e.g., Kong et al., 2017; 2019) that can enhance acceleration.
Intensities of GeV protons only exceed GCRs for short periods and are too weak to be a practical radiation hazard. However, ∼100-MeV protons are much more numerous and persistent and hence a significant hazard to astronauts outside the Earth’s magnetic fields. When an event is near central meridian, intensities of high-energy SEPs can be extended and increased by the ESP event when the nose of the shock passes over us. Some historical examples are shown in Figure 8. For western sources, we see the greatest effect of the shock nose early, but the shock then weakens toward our longitude and weakens in strength out to 1 AU. The early SEP intensities at 1 AU are constrained by wave growth, which can establish the “streaming limit” early (Reames and Ng, 1998; Ng et al., 1999; Ng et al., 2003; Reames and Ng, 2010; Ng et al., 2012; Reames and Ng, 2014), while ESP intensities are unbounded. Thus, the spatial distribution of SEP intensities at the shock, i.e., the ESP event, is of both fundamental and practical importance. The SEP events of October 1989 shown in Figure 8 are the basis of SEP “storm shelter” radiation-shielding requirements for astronauts in missions beyond a low Earth orbit (Townsend et al., 2018).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1254266/fspas-10-1254266-HTML/image_m/fspas-10-1254266-g008.jpg)
FIGURE 8. High-energy proton intensities during historic large SEP events near the central meridian can be enhanced and extended in time by their ESP events, even at energies >100 MeV (a significant radiation hazard). Early intensities are bounded by the “streaming limit,” but the ESP peak is not. All of these events are GLEs, which often occur near the central meridian, allowing the ESP to surface.
The ESP structure is formed very early in an event when protons streaming away from the shock amplify resonant Alfvén waves (Melrose, 1980; Stix, 1992). Waves trap ions of a given rigidity, scattering them back and forth across the shock so they gain velocity on each transit, subsequently amplifying waves of lower k (longer wavelength) which trap ions of higher P, etc., creating the ESP structure of energetic ions trapped around the shock (Lee, 1983; Lee, 2005). At the streaming limit, higher SEP intensities simply grow more waves, trapping more ions back near the shock. In gradual SEP events, the ESP structure always exists. When the shock is near the Sun, the ESP is initially hidden among the same SEPs that form it as they stream away. The shock we see at 1 AU was much stronger when it began near the Sun. The shock and the ESP event can explain all of the SEPs in the gradual event, early and late. It is a key to the physics of SEP acceleration in these events, whether it emerges at your particular longitude or not. The basic ESP structure will only be exposed at some longitudes in some events.
When the flank of a shock crosses to new field lines, where the early emission is absent, the “naked” ESP event emerges, as seen at Voyager in Figure 2E or at IMP 8 in Figure 1E. The ESP structure persists as the shock moves outward, but the CME speed decreases, and decreasing B shifts the resonance so the highest energy ions preferentially begin to leak away. Earlier peaks of ∼30–100-MeV protons are seen along with the naked ESP events in these figures because they leaked out after the ESP events have crossed to the new field lines.
Generally, the time of the peak intensity is a contest: the initial intensities are bounded at the streaming limit, while the ESP peak is unbounded, yet the highest energies in the ESP begin to leak away as the shock slows and moves out to lower B. Understanding the high-energy particle acceleration of SEPs is about understanding shock acceleration, wave trapping, and the formation, spatial extent, and persistence of ESP events. The large events shown in Figure 8 are relatively rare, and we tend to assume that the same physics can be studied at lower energies. However, the physics changes for those energies where resonant waves dominate—a function of space and time. Following this, few general questions do arise:
1) Comparing Figures 1C, E, how important is θBn in maintaining the ESP event?
2) The increasing intensity ramp at Voyager in Figure 2A involves SEPs leaked from the approaching ESP event. Does it show strengthening of the shock with longitude toward the east or mainly increasing proximity?
3) Are the SEP3 events with enhanced heavy ions limited in longitude? Surely, seed-particle populations in an event can change with longitude; where are the events that change from SEP3 to SEP4? Why are STEREO three-spacecraft events all SEP4 events? Half of GLEs are SEP3s.
4) How does the longitude extent of the SEPs compare with that of the shock itself? How does each vary with time?
5) Are GLEs or other high-energy events more limited in longitude? Would they be one- or two-spacecraft events for STEREO? Are broader shocks weaker or less efficient?
We should expect that shocks can accelerate whatever seed population or populations they encounter, wherever they go. We cannot exclude mixtures; we can only distinguish which one dominates at high Z and thus label an event SEP3 or SEP4.
We are constantly hampered by the correlated mixture of space and time. The footpoint of a field line from Earth lies 50–60⁰ to our west. As our connection point on the shock scans to the east, this shock also weakens with time. How much is the change with longitude; how much is the change with time? The only resolution is to measure spatial variations with a scale substantially less than 50⁰. By the time the shock arrives at 1 AU, it has mixed SEPs from ≈50⁰ of solar longitude; in longer time intervals, solar rotation causes greater mixing.
5 What is needed next?
SEP evolution in space and time is complicated. STEREO spacecraft separation from Earth of ∼120⁰ was much too coarse to resolve the SEP-shock evolution we happened to observe in a Helios-IMP-Voyager period. What kind of observations would help? Multiple points with a better spatial resolution would help. Consider two primary spacecraft, equipped as STEREO was, with coronagraphs and in situ instruments, which are fixed in Earth’s solar orbit at ±60⁰ from Earth, as shown in Figure 9. On each side, between each primary spacecraft and Earth are two much smaller spin-stabilized spacecraft (similar to Wind), 20⁰ apart, each capable of measuring SEPs, magnetic fields, and solar-wind plasma to map SEP events, shocks, and interplanetary CME structures. This configuration presumes that capabilities similar to a primary spacecraft are preexisting near Earth.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1254266/fspas-10-1254266-HTML/image_m/fspas-10-1254266-g009.jpg)
FIGURE 9. A possible multispacecraft configuration of the future might include two primary spacecraft, similar to STEREO, and four small spacecraft measuring SEPs, plasma, and fields, to study in situ spatial distributions of SEPs, interplanetary CMEs, and shock waves. The spacecraft would maintain their positions relative to Earth for years, spanning at least one solar maximum.
The ≈20⁰ spacing would provide at least two measurements between the longitude of a spacecraft and that of its typical solar magnetic footpoint. It would allow meaningful coverage of many spatially small- and moderate-sized events with sampling of variations along the shock for the larger-sized events. “Smaller” SEPs do not mean weaker. Of course, it would be better to have more complete coverage, but the configuration in Figure 9 represents a major improvement in spatial resolution, and reasonable (33%) coverage, at a modest increase in cost and complexity over the dual-spacecraft STEREO mission. A spacing of ≤20⁰ seems essential and provides a reasonable tradeoff between spacing and coverage. This would provide a map of the shock strength, direction of propagation, and θBn at up to seven points at 1 AU that could be compared with the coronagraph mapping near the Sun and could give seven SEP/ESP profiles with differing onset times and intensities. What does the shock front really look like? Is the nose of the shock a hot spot that could produce a localized GLE; for what longitudes does the peak intensity arrive with the shock?
Are the highest energies in SEP events limited to short time periods and small spatial intervals as Gopalswamy et al. (2013) conclude? If so, what physics defines those intervals, and why do they differ from the apparent STEREO finding of similarly broad Gaussians for most energy bands? Up to some energy, we would expect high SEP intensities to generate enough waves to extend the spatial trapping and the duration of an event (e.g., Figure 8), although not at the absolute peak energy where there are yet few resonant waves. These peaks could be highly localized in well-connected regions that become poorly correlated with the average properties of the CME. Moreover, in some events, particle trapping is increased by the presence of multiple shock waves, as in the 14 July 2000 “Bastille Day” GLE (Lepping et al., 2001). To what extent can the SEP energy profile and the peak energy be predicted from an early coronagraph map of the shock strength? We would have an extra day to predict the strength of the ESP event.
Author contributions
DR: conceptualization, writing-original manuscript, and writing-review and editing.
Acknowledgments
The author thanks Ed Cliver for helpful discussions.
Conflict of interest
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Beeck, J., Mason, G. M., Hamilton, D. C., Wibberenz, G., Kunow, H., Hovestadt, D., et al. (1987). A multispacecraft study of the injection and transport of solar energetic particles. Astrophys. J. 322, 1052. doi:10.1086/165800
Bell, A. R. (1978a). The acceleration of cosmic rays in shock fronts. I. Mon. Not. Roy. Astron. Soc. 182, 147–156. doi:10.1093/mnras/182.2.147
Bell, A. R. (1978b). The acceleration of cosmic rays in shock fronts. II. Mon. Not. Roy. Astron. Soc. 182, 443–455. doi:10.1093/mnras/182.3.443
Bertsch, D. L., Fichtel, C. E., and Reames, D. V. (1969). Relative abundance of iron-group nuclei in solar cosmic rays. Astrophys. J. Lett. 157, L53. doi:10.1086/180383
Bučík, R. (2020). 3He-rich solar energetic particles: solar sources. Space Sci. Rev. 216, 24. doi:10.1007/s11214-020-00650-5
Bučík, R., Innes, D. E., Chen, N. H., Mason, G. M., Gómez-Herrero, R., and Wiedenbeck, M. E. (2015). Long-lived energetic particle source regions on the Sun. J. Phys. Conf. Ser. 642, 012002. doi:10.1088/1742-6596/642/1/012002
Bučík, R., Innes, D. E., Mall, U., Korth, A., Mason, G. M., and Gómez-Herrero, R. (2014). Multi-spacecraft observations of recurrent 3He-rich solar energetic particles. Astrophys. J. 786, 71. doi:10.1088/0004-637X/786/1/71
Bučík, R., Innes, D. E., Mason, G. M., Wiedenbeck, M. E., Gómez-Herrero, R., and Nitta, N. V. (2018a). 3He-rich solar energetic particles in helical jets on the Sun. Astrophys. J. 852, 76. doi:10.3847/1538-4357/aa9d8f
Bučík, R., Mulay, S. M., Mason, G. M., Nitta, N. V., Desai, M. I., and Dayeh, M. A. (2021). Temperature in solar sources of 3He-rich solar energetic particles and relation to ion abundances. Astrophys. J. 908, 243. doi:10.3847/1538-4357/abd62d
Bučík, R., Wiedenbeck, M. E., Mason, G. M., Gómez-Herrero, R., Nitta, N. V., and Wang, L. (2018b). 3He-rich solar energetic particles from sunspot jets. Astrophys. J. Lett. 869, L21. doi:10.3847/2041-8213/aaf37f
Bukata, R. P., GronstalPalmeira, P. T. R. A. R., McCracken, K. G., and Rao, U. R. (1969). Neutron monitor and Pioneer 6 and 7 studies of the January 28, 1967 solar flare event. Sol. Phys. 10, 198–211. doi:10.1007/BF00146170
Burlaga, L., Sittler, E., Mariani, F., and Schwenn, R. (1981). Magnetic loop behind an interplanetary shock: voyager, helios, and IMP 8 observations. J. Geophys. Res. 86, 6673. doi:10.1029/JA086iA08p06673
Cane, H. V., Reames, D. V., and von Rosenvinge, T. T. (1988). The role of interplanetary shocks in the longitude distribution of solar energetic particles. J. Geophys. Res. 93, 9555. doi:10.1029/JA093iA09p09555
Chen, N. H., Bučík, R., Innes, D. E., and Mason, G. M. (2015). Case studies of multi-day 3He-rich solar energetic particle periods. Astron. Astrophys. 580, 16. doi:10.1051/0004-6361/201525618
Cliver, E. W., Kahler, S. W., and Reames, D. V. (2004). Coronal shocks and solar energetic proton events. Astrophys. J. 605, 902–910. doi:10.1086/382651
Cliver, E. W., and Ling, A. G. (2007). Electrons and protons in solar energetic particle events. Astrophys. J. 658, 1349–1356. doi:10.1086/511737
Cohen, C. M. S., Mason, G. M., and Mewaldt, R. A. (2017). Characteristics of solar energetic ions as a function of longitude. Astrophys. J. 843, 132. doi:10.3847/1538-4357/aa7513
Cook, W. R., Stone, E. C., and Vogt, R. E. (1984). Elemental composition of solar energetic particles. Astrophys. J. 279, 827. doi:10.1086/161953
Daibog, E. I., Stolpovskii, V. G., and Kahler, S. W. (2003). Invariance of charged particle time profiles at late stages of SCR events from the data of multisatellite observations. Cosm. Res. 41, 128.
de Nolfo, G. A., Bruno, A., Ryan, J. M., Dalla, S., Giacalone, J., Richardson, I. G., et al. (2019). Comparing long-duration gamma-ray flares and high-energy solar energetic particles. Astrophys. J. 879 (90), 90. doi:10.3847/1538-4357/ab258f
Desai, M. I., and Giacalone, J. (2016). Large gradual solar energetic particle events. Living Rev. Sol. Phys. 13, 3. doi:10.1007/s41116-016-0002-5
Desai, M. I., Mason, G. M., Dwyer, J. R., Mazur, J. E., Gold, R. E., Krimigis, S. M., et al. (2003). Evidence for a suprathermal seed population of heavy ions accelerated by interplanetary shocks near 1 AU. Astrophys. J. 588, 1149–1162. doi:10.1086/374310
Drake, J. F., Cassak, P. A., Shay, M. A., Swisdak, M., and Quataert, E. (2009). A magnetic reconnection mechanism for ion acceleration and abundance enhancements in impulsive flares. Astrophys. J. Lett. 700, L16–L20. doi:10.1088/0004-637X/700/1/L16
Dresing, N., Gómez-Herrero, R., Heber, B., Klassen, A., Temmer, M., and Veronig, A. (2018). Long-lasting injection of solar energetic electrons into the heliosphere. Astron. Astrophys. 613, A21. doi:10.1051/0004-6361/201731573
Fichtel, C. E., and Guss, D. E. (1961). Heavy nuclei in solar cosmic rays. Phys. Rev. Lett. 6, 495–497. doi:10.1103/PhysRevLett.6.495
Fisk, L. A. (1978). He-3-rich flares - a possible explanation. Astrophys. J. 224, 1048. doi:10.1086/156456
Forbush, S. E. (1946). Three unusual cosmic ray increases possibly due to charged particles from the Sun. Phys. Rev. 70, 771–772. doi:10.1103/physrev.70.771
Gopalswamy, N., Xie, H., Akiyama, S., Yashiro, S., Usoskin, I. G., and Davila, J. M. (2013). The first ground level enhancement event of solar cycle 24: direct observation of shock formation and particle release heights. Astrophys. J. Lett. 765, L30. doi:10.1088/2041-8205/765/2/L30
Gopalswamy, N., Xie, H., Yashiro, S., Akiyama, S., Mäkelä, P., and Usoskin, I. G. (2012). Properties of Ground level enhancement events and the associated solar eruptions during solar cycle 23. Space Sci. Rev. 171, 23–60. doi:10.1007/s11214-012-9890-4
Gosling, J. T. (1994). Correction to “The solar flare myth”. J. Geophys. Res. 99, 4259. doi:10.1029/94JA00015
Gosling, J. T. (1993). The solar flare myth. J. Geophys. Res. 98, 18937–18949. doi:10.1029/93JA01896
Heras, A. M., Sanahuja, B., Lario, D., Smith, Z. K., Detman, T., and Dryer, M. (1995). Three low-energy particle events: modeling the influence of the parent interplanetary shock. Astrophys. J. 445, 497. doi:10.1086/175714
Jokipii, J. R., and Parker, E. N. (1969). Stochastic aspects of magnetic lines of force with application to cosmic ray propagation. Astrophys. J. 155, 777. doi:10.1086/149909
Kahler, S. W., Ling, A. G., and Reames, D. V. (2023). Spatial evolution of 20-MeV solar energetic proton events. Astrophys. J. 942, 68. doi:10.3847/1538-4357/aca7c0
Kahler, S. W., Reames, D. V., and Sheeley, N. R. (2001). Coronal mass ejections associated with impulsive solar energetic particle events. Astrophys. J. 562, 558–565. doi:10.1086/323847
Kahler, S. W., Sheeley, N. R., Howard, R. A., Koomen, M. J., Michels, D. J., McGuire, R. E., et al. (1984). Associations between coronal mass ejections and solar energetic proton events. J. Geophys. Res. 89, 9683. doi:10.1029/JA089iA11p09683
Kahler, S. W. (1992). Solar flares and coronal mass ejections. Ann. Rev. Astron. Astrophys. 30, 113–141. doi:10.1146/annurev.aa.30.090192.000553
Kahler, S. W. (1982). The role of the big flare syndrome in correlations of solar energetic proton fluxes and associated microwave burst parameters. J. Geophys. Res. 87, 3439. doi:10.1029/JA087iA05p03439
Kecskeméty, K., Daibog, E. I., Logachev, Y. I., and Kóta, J. (2009). The decay phase of solar energetic particle events. J. Geophys. Res. A 114, 6102. doi:10.1029/2008JA013730
Klassen, A., Dresing, N., Gómez-Herrero, R., and Heber, B. (2015). First simultaneous observations of a near-relativistic electron spike event by both STEREO spacecraft. Astron. Astrophys. 580, A115. doi:10.1051/0004-6361/201525700
Klassen, A., Dresing, N., Gómez-Herrero, R., Heber, B., and Müller-Mellin, R. (2016). Unexpected spatial intensity distributions and onset timing of solar electron events observed by closely spaced STEREO spacecraft. Astron. Astrophys. 593, A31. doi:10.1051/0004-6361/201628734
Kong, X., Guo, F., Chen, Y., and Giacalone, J. (2019). The acceleration of energetic particles at coronal shocks and emergence of a double power-law feature in particle energy spectra. Astrophys. J. 883, 49. doi:10.3847/1538-4357/ab3848
Kong, X., Guo, F., Giacalone, J., Li, H., and Chen, Y. (2017). The acceleration of high-energy protons at coronal shocks: the effect of large-scale streamer-like magnetic field structures. Astrophys. J. 851, 38. doi:10.3847/1538-4357/aa97d7
Kouloumvakos, A., Rouillard, A. P., Wu, Y., Vainio, R., Vourlidas, A., Plotnikov, I., et al. (2019). Connecting the properties of coronal shock waves with those of solar energetic particles. Astrophys. J. 876, 80. doi:10.3847/1538-4357/ab15d7
Laming, J. M., Vourlidas, A., Korendyke, C., Chua, D., Cranmer, S. R., Ko, Y. K., et al. (2019). Element abundances: A new diagnostic for the solar wind. Astrophys. J. 879, 124. doi:10.3847/1538-4357/ab23f1
Lario, D., Aran, A., Agueda, N., and Sanahuja, B. (2007). Radial dependence of proton peak intensities and fluences in SEP events: influence of the energetic particle transport parameters. AdSpR 40, 289–294. doi:10.1016/j.asr.2007.01.057
Lario, D., Aran, A., Gómez-Herrero, R., Dresing, N., Heber, B., Ho, G. C., et al. (2013). Longitudinal and radial dependence of solar energetic particle peak intensities: STEREO, ACE, SOHO, GOES, and MESSENGER observations. Astrophys. J. 767, 41. doi:10.1088/0004-637X/767/1/41
Lario, D., Kallenrode, M.-B., Decker, R. B., Roelof, E. C., Krimigis, S. M., Aran, A., et al. (2006). Radial and longitudinal dependence of solar 4-13 MeV and 27-37 MeV proton peak intensities and fluences: helios and IMP 8 observations. Astrophys. J. 653, 1531–1544. doi:10.1086/508982
Lario, D., Maksimovic, M., Issautier, K., Meyer-Vernet, N., Moncuquet, M., and Pantellini, F. (2010). Heliospheric energetic particle reservoirs: ulysses and ACE 175-315 keV electron observations. Proc. 12th Sol. Wind Conf. AIP Conf. Proc. 1216, 625. doi:10.1063/1.3395944
Lee, M. A. (2005). Coupled hydromagnetic wave excitation and ion acceleration at an evolving coronal/interplanetary shock. Astrophys. J. Suppl. 158, 38–67. doi:10.1086/428753
Lee, M. A. (1983). Coupled hydromagnetic wave excitation and ion acceleration at interplanetary traveling shocks. J. Geophys. Res. 88, 6109. doi:10.1029/JA088iA08p06109
Lee, M. A., Mewaldt, R. A., and Giacalone, J. (2012). Shock acceleration of ions in the heliosphere. Space Sci. Rev. 173, 247–281. doi:10.1007/s11214-012-9932-y
Lepping, R. P., Berdichevsky, D. B., Burlaga, L. F., Lazarus, A. J., Kasper, J., Desch, M. D., et al. (2001). The bastille day magnetic clouds and upstream shocks: near earth interplanetary observations. Sol. Phys. 204, 285–303. doi:10.1023/A:1014264327855
Li, G., and Bian, N. H. (2023). Lagrangian stochastic model for the motions of magnetic footpoints on the solar wind source surface and the path lengths of boundary-driven interplanetary magnetic field lines. Astrophys. J. 945, 150. doi:10.3847/1538-4357/acbd43
Liu, S., Petrosian, V., and Mason, G. M. (2004). Stochastic acceleration of 3He and 4He in solar flares by parallel-propagating plasma waves. Astrophys. J. Lett. 613, L81–L84. doi:10.1086/425070
Liu, S., Petrosian, V., and Mason, G. M. (2006). Stochastic acceleration of 3He and 4He in solar flares by parallel-propagating plasma waves: general results. Astrophys. J. 636, 462–474. doi:10.1086/497883
Liu, W., Kong, X., Guo, F., Zhao, L., Feng, S., Yu, F., et al. (2023). Effects of coronal magnetic field configuration on particle acceleration and release during the ground level enhancement events in solar cycle 24. Available at: https://arxiv.org/abs/2307.12191.
Mandzhavidze, N., Ramaty, R., and Kozlovsky, B. (1999). Determination of the abundances of subcoronal 4He and of solar flare-accelerated 3He and 4He from gamma-ray spectroscopy. Astrophys. J. 518, 918–925. doi:10.1086/307321
Mason, G. M. (2007). 3He-rich solar energetic particle events. Space Sci. Rev. 130, 231–242. doi:10.1007/s11214-007-9156-8
Mason, G. M., Gloeckler, G., and Hovestadt, D. (1984). Temporal variations of nucleonic abundances in solar flare energetic particle events. II - evidence for large-scale shock acceleration. Astrophys. J. 280, 902. doi:10.1086/162066
Mason, G. M., Mazur, J. E., Dwyer, J. R., Jokippi, J. R., Gold, R. E., and Krimigis, S. M. (2004). Abundances of heavy and ultraheavy ions in 3He-rich solar flares. Astrophys. J. 606, 555–564. doi:10.1086/382864
Mason, G. M., Mazur, J. E., and Dwyer, J. R. (1999). [TSUP]3[/TSUP]H[CLC]e[/CLC] Enhancements in Large Solar Energetic Particle Events. Astrophys. J. Lett. 525, L133–L136. doi:10.1086/312349
Mason, G. M., Ng, C. K., Klecker, B., and Green, G. (1989). Impulsive acceleration and scatter-free transport of about 1 MeV per nucleon ions in 3He-rich solar particle events. Astrophys. J. 339, 529. doi:10.1086/167315
Mason, G. M., Reames, D. V., Klecker, B., Hovestadt, D., and von Rosenvinge, T. T. (1986). The heavy-ion compositional signature in He-3-rich solar particle events. Astrophys. J. 303, 849. doi:10.1086/164133
McGuire, R. E., von Rosenvinge, T. T., and McDonald, F. B. (1979). A survey of solar cosmic ray composition. Proc. 16thInt. Cosm. Ray Conf. 5, 61.
McKibben, R. B. (1972). Azimuthal propagation of low-energy solar-flare protons as observed from spacecraft very widely separated in solar azimuth. J. Geophys. Res. 77, 3957–3984. doi:10.1029/JA077i022p03957
Mewaldt, R. A., Looper, M. D., Cohen, C. M. S., Haggerty, D. K., Labrador, A. W., Leske, R. A., et al. (2012). Energy spectra, composition, and other properties of ground-level events during solar cycle 23. Space Sci. Rev. 171, 97–120. doi:10.1007/s11214-012-9884-2
Meyer, J. P. (1985). The baseline composition of solar energetic particles. Astrophys. J. Suppl. 57, 151. doi:10.1086/191000
Murphy, R. J., Kozlovsky, B., and Share, G. H. (2016). Evidence for enhanced 3He in flare-accelerated particles based on new calculations of the gamma-ray line spectrum. Astrophys. J. 833, 196. doi:10.3847/1538-4357/833/2/196
Murphy, R. J., Ramaty, R., Kozlovsky, B., and Reames, D. V. (1991). Solar abundances from gamma-ray spectroscopy: comparisons with energetic particle, photospheric, and coronal abundances. Astrophys. J. 371, 793. doi:10.1086/169944
Newkirk, G., and Wenzel, D. G. (1978). Rigidity-independent propagation of cosmic rays in the solar corona. J. Geophys. Res. 83, 2009. doi:10.1029/JA083iA05p02009
Ng, C. K., Reames, D. V., and Tylka, A. J. (1999). Effect of proton-amplified waves on the evolution of solar energetic particle composition in gradual events. Geophys. Res. Lett. 26, 2145–2148. doi:10.1029/1999GL900459
Ng, C. K., and Reames, D. V. (2008). Shock acceleration of solar energetic protons: the first 10 minutes. Astrophys. J. Lett. 686, L123–L126. doi:10.1086/592996
Ng, C. K., Reames, D. V., and Tylka, A. J. (2003). Modeling shock-accelerated solar energetic particles coupled to interplanetary Alfvén waves. Astrophys. J. 591, 461–485. doi:10.1086/375293
Ng, C. K., Reames, D. V., and Tylka, A. J. (2012). Solar energetic particles: shock acceleration and transport through self-amplified waves. AIP Conf. Proc. 1436, 212. doi:10.1063/1.4723610
Nitta, N. V., Mason, G. M., Wang, L., Cohen, C. M. S., and Wiedenbeck, M. E. (2015). Solar sources of 3He-rich solar energetic particle events in Solar Cycle 24. Astrophys. J. 806, 235. doi:10.1088/0004-637X/806/2/235
Nitta, N. V., Reames, D. V., DeRosa, M. L., Yashiro, S., and Gopalswamy, N. (2006). Solar sources of impulsive solar energetic particle events and their magnetic field connection to the earth. Astrophys. J. 650, 438–450. doi:10.1086/507442
Paassilta, M., Papaioannou, A., Dresing, N., Vainio, R., Valtonen, E., and Heber, B. (2018). Catalogue of > 55 MeV wide-longitude solar proton events observed by SOHO, ACE, and the STEREOs at ≈ 1 AU during 2009 – 2016. Sol. Phys. 293, 70. doi:10.1007/s11207-018-1284-7
Parker, E. N. (1965). The passage of energetic charged particles through interplanetary space. Space Sci. 13, 9–49. doi:10.1016/0032-0633(65)90131-5
Raukunen, O., Vainio, R., Tylka, A. J., Dietrich, W. F., Jiggens, P., Heynderickx, D., et al. (2018). Two solar proton fluence models based on ground level enhancement observations. J. Spa. Wea. Spa. Clim. 8, A04. doi:10.1051/swsc/2017031
Reames, D. V. (2000). Abundances of trans-iron elements in solar energetic particle events. Astrophys. J. Lett. 540, L111–L114. doi:10.1086/312886
Reames, D. V. (2018b). Abundances, ionization states, temperatures, and FIP in solar energetic particles. Space Sci. Rev. 214, 61. doi:10.1007/s11214-018-0495-4
Reames, D. V., Barbier, L. M., and Ng, C. K. (1996). The spatial distribution of particles accelerated by coronal mass ejection-driven shocks. Astrophys. J. 466, 473. doi:10.1086/177525
Reames, D. V. (1988). Bimodal abundances in the energetic particles of solar and interplanetary origin. J. Lett. 330, L71. doi:10.1086/185207
Reames, D. V., Cliver, E. W., and Kahler, S. W. (2014a). Abundance enhancements in impulsive solar energetic-particle events with associated coronal mass ejections. Sol. Phys. 289, 3817–3841. doi:10.1007/s11207-014-0547-1
Reames, D. V., Cliver, E. W., and Kahler, S. W. (2014b). Variations in abundance enhancements in impulsive solar energetic-particle events and related CMEs and flares. Sol. Phys. 289, 4675–4689. doi:10.1007/s11207-014-0589-4
Reames, D. V. (1995a). Coronal Abundances determined from energetic particles. Adv. Space Res. 15 (7), 41–51. doi:10.1016/0273-1177(94)00018-v
Reames, D. V. (2014). Element abundances in solar energetic particles and the solar corona. Sol. Phys. 289, 977–993. doi:10.1007/s11207-013-0350-4
Reames, D. V. (2022a). Energy spectra vs. element abundances in solar energetic particles and the roles of magnetic reconnection and shock acceleration. Sol. Phys. 297, 32. doi:10.1007/s11207-022-01961-2
Reames, D. V. (2020). Four distinct pathways to the element abundances in solar energetic particles. Space Sci. Rev. 216, 20. doi:10.1007/s11214-020-0643-5
Reames, D. V. (2023). How do shock waves define the space-time structure of gradual solar energetic-particle events? Space Sci. Rev. 219, 14. doi:10.1007/s11214-023-00959-x
Reames, D. V. (2019). Hydrogen and the abundances of elements in impulsive solar energetic-particle events. Sol. Phys. 294, 37. doi:10.1007/s11207-019-1427-5
Reames, D. V., Kahler, S. W., and Ng, C. K. (1997). Spatial and temporal invariance in the spectra of energetic particles in gradual solar events. Astrophys. J. 491, 414–420. doi:10.1086/304939
Reames, D. V., Kallenrode, M.-B., and Stone, R. G. (1991). Multi-spacecraft observations of solar 3He-rich events. Astrophys. J. 380, 287. doi:10.1086/170585
Reames, D. V., and Lal, N. (2010). A multi-spacecraft view of solar-energetic-particle onsets in the 1977 November 22 event. Astrophys. J. 723, 550–554. doi:10.1088/0004-637X/723/1/550
Reames, D. V., Meyer, J. P., and von Rosenvinge, T. T. (1994). Energetic-particle abundances in impulsive solar flare events. Astrophys. J. Suppl. 90, 649. doi:10.1086/191887
Reames, D. V., and Ng, C. K. (2004). Heavy-element abundances in solar energetic particle events. Astrophys. J. 610, 510–522. doi:10.1086/421518
Reames, D. V., and Ng, C. K. (1998). Streaming-limited intensities of solar energetic particles. Astrophys. J. 504, 1002–1005. doi:10.1086/306124
Reames, D. V., and Ng, C. K. (2010). Streaming-limited intensities of solar energetic particles on the intensity plateau. Astrophys. J. 723, 1286–1293. doi:10.1088/0004-637X/723/2/1286
Reames, D. V., and Ng, C. K. (2014). The streaming limit of solar energetic-particle intensities. Living a Star Workshop Extreme Space Weather Events. Available at: https://arxiv.org/abs/1412.2279.
Reames, D. V., Ng, C. K., and Tylka, A. J. (2013). Spatial distribution of solar energetic particles in the inner heliosphere. Sol. Phys. 285, 233–250. doi:10.1007/s11207-012-0038-1
Reames, D. V. (1999). Particle acceleration at the Sun and in the heliosphere. Space Sci. Rev. 90, 413–491. doi:10.1023/A:1005105831781
Reames, D. V. (2021b). Sixty years of element abundance measurements in solar energetic particles. Space Sci. Rev. 217, 72. doi:10.1007/s11214-021-00845-4
Reames, D. V. (2021a). “Solar energetic particles,” in Lec. Notes phys. 978 (Cham, Switzerland: Springer Nature). open access. doi:10.1007/978-3-030-66402-2
Reames, D. V. (1995b). Solar energetic particles: A paradigm shift. Revs. Geophys. Suppl. 33, 585. doi:10.1029/95RG00188
Reames, D. V. (2022b). Solar energetic particles: spatial extent and implications of the H and He abundances. Space Sci. Rev. 218, 48. doi:10.1007/s11214-022-00917-z
Reames, D. V. (2009b). Solar energetic-particle release times in historic ground-level events. Astrophys. J. 706, 844–850. doi:10.1088/0004-637X/706/1/844
Reames, D. V. (2009a). Solar release times of energetic particles in ground-level events. Astrophys. J. 693, 812–821. doi:10.1088/0004-637X/693/1/812
Reames, D. V., and Stone, R. G. (1986). The identification of solar 3He-rich events and the study of particle acceleration at the sun. Astrophys. J. 308, 902. doi:10.1086/164560
Reames, D. V. (2016). Temperature of the source plasma in gradual solar energetic particle events. Sol. Phys. 291, 911–930. doi:10.1007/s11207-016-0854-9
Reames, D. V. (2018a). The “FIP effect” and the origins of solar energetic particles and of the solar wind. Sol. Phys. 293, 47. doi:10.1007/s11207-018-1267-8
Reames, D. V. (2013). The two sources of solar energetic particles. Space Sci. Rev. 175, 53–92. doi:10.1007/s11214-013-9958-9
Reames, D. V., von Rosenvinge, T. T., and Lin, R. P. (1985). Solar 3He-rich events and nonrelativistic electron events - a new association. Astrophys. J. 292, 716. doi:10.1086/163203
Reinhard, R., and Wibberenz, G. (1974). Propagation of flare protons in the solar atmosphere. Sol. Phys. 36, 473–494. doi:10.1007/BF00151216
Richardson, I. G., Reames, D. V., Wenzel, K. P., and Rodriguez-Pacheco, J. (1990). Quiet-time properties of low-energy (less than 10 MeV per nucleon) interplanetary ions during solar maximum and solar minimum. Astrophys. J. Lett. 363, L9. doi:10.1086/185853
Richardson, I. G., von Rosenvinge, T. T., Cane, H. V., Christian, E. R., Cohen, C. M. S., Labrador, A. W., et al. (2014). 25 MeV proton events observed by the high energy telescopes on the stereo a and b spacecraft and/or at Earth during the first ∼ seven years of the STEREO mission, Sol. Phys. 289 3059–3107. doi:10.1007/s11207-014-0524-8
Roelof, E. C., Gold, R. E., Simnett, G. M., Tappin, S. J., Armstrong, T. P., and Lanzerotti, L. J. (1992). Low-energy solar electrons and ions observed at ULYSSES February-April, 1991 - the inner heliosphere as a particle reservoir. Geophys. Res. Lett. 19, 1243–1246. doi:10.1029/92GL01312
Rouillard, A. C., Odstrčil, D., Sheeley, N. R., Tylka, A. J., Vourlidas, A., Mason, G., et al. (2011). Interpreting the properties of solar energetic particle events by using combined imaging and modeling of interplanetary shocks. Astrophys. J. 735, 7. doi:10.1088/0004-637X/735/1/7
Rouillard, A., Sheeley, N. R., Tylka, A., Vourlidas, A., Ng, C. K., Rakowski, C., et al. (2012). The longitudinal properties of a solar energetic particle event investigated using modern solar imaging. Astrophys. J. 752, 44. doi:10.1088/0004-637X/752/1/44
Serlemitsos, A. T., and Balasubrahmanyan, V. K. (1975). Solar particle events with anomalously large relative abundance of 3He. Astrophys. J. 198, 195. doi:10.1086/153592
Temerin, M., and Roth, I. (1992). The production of 3He and heavy ion enrichment in 3He-rich flares by electromagnetic hydrogen cyclotron waves. Astrophys. J. Lett. 391, L105. doi:10.1086/186408
Townsend, L. W., Adams, J. H., Blattnig, S. R., Clowdsley, M. S., Fry, D. J., Jun, I., et al. (2018). Solar particle event storm shelter requirements for missions beyond low Earth orbit. Life Sci. Space Res. 17, 32–39. doi:10.1016/j.lssr.2018.02.002
Tylka, A. J., Cohen, C. M. S., Dietrich, W. F., Lee, M. A., Maclennan, C. G., Mewaldt, R. A., et al. (2005). Shock geometry, seed populations, and the origin of variable elemental composition at high energies in large gradual solar particle events. Astrophys. J. 625, 474–495. doi:10.1086/429384
Tylka, A. J., Cohen, C. M. S., Dietrich, W. F., Krucker, S., McGuire, R. E., Mewaldt, R. A., et al. (2003). “Inspirehep,” in Proceedingof the 28th International Cosmic Ray Conference (ICRC 2003) - INSPIRE, Tsukuba, Japan, August 2003, 3305.
Tylka, A. J., Cohen, C. M. S., Dietrich, W. F., Maclennan, C. G., McGuire, R. E., Ng, C. K., et al. (2001). Evidence for remnant flare suprathermals in the source population of solar energetic particles in the 2000 bastille day event. Astrophys. J. Lett. 558, L59–L63. doi:10.1086/323344
Tylka, A. J., and Dietrich, W. F. (2009). “A new and comprehensive analysis of proton spectra in ground-level enhanced (GLE) solar particle events,” in Proceedings of 31st International Cosmic Ray Conference Lodz, Lodz, Poland, July 2009.
Tylka, A. J., and Lee, M. A. (2006). A model for spectral and compositional variability at high energies in large, gradual solar particle events. Astrophys. J. 646, 1319–1334. doi:10.1086/505106
Wang, Y.-M., Pick, M., and Mason, G. M. (2006). Coronal holes, jets, and the origin of 3He-rich particle events. Astrophys. J. 639, 495–509. doi:10.1086/499355
Wiedenbeck, M. E., Mason, G. M., Cohen, C. M. S., Nitta, N. V., Gómez-Herrero, R., and Haggerty, D. K. (2013). Observations of solar energetic particles from 3He-rich events over a wide range of heliographic longitude. Astrophys. J. 762, 54. doi:10.1088/0004-637X/762/1/54
Wild, J. P., Smerd, S. F., and Weiss, A. A. (1963). Solar bursts. Annu. Rev. Astron. Astrophys. 1, 291–366. doi:10.1146/annurev.aa.01.090163.001451
Xie, H., Cyr, O. C., Mäkelä, P., and Gopalswamy, N. (2019). Statistical study on multispacecraft widespread solar energetic particle events during solar cycle 24. J. Geophys. Res. 124, 6384–6402. doi:10.1029/2019JA026832
Zank, G. P., Li, G., and Verkhoglyadova, O. (2007). Particle acceleration at interplanetary shocks. Space Sci. Rev. 130, 255–272. doi:10.1007/s11214-007-9214-2
Zank, G. P., Rice, W. K. M., and Wu, C. C. (2000). Particle acceleration and coronal mass ejection driven shocks: A theoretical model. J. Geophys. Res. 105, 25079–25095. doi:10.1029/1999JA000455
Zhang, M., Cheng, L., Zhang, J., Riley, P., Kwon, R. Y., Lario, D., et al. (2023). A data-driven, physics-based transport model of solar energetic particles accelerated by coronal mass ejection shocks propagating through the solar coronal and heliospheric magnetic fields. Astrophys. J. Suppl. 266, 35. doi:10.3847/1538-4365/accb8e
Keywords: solar energetic particles, shock waves, coronal mass ejections, solar jets, solar system abundances, multispacecraft missions, heliosphere
Citation: Reames DV (2023) Review and outlook of solar energetic particle measurements on multispacecraft missions. Front. Astron. Space Sci. 10:1254266. doi: 10.3389/fspas.2023.1254266
Received: 06 July 2023; Accepted: 11 August 2023;
Published: 31 August 2023.
Edited by:
David Gary Sibeck, National Aeronautics and Space Administration, United StatesReviewed by:
Gang Li, University of Alabama in Huntsville, United StatesNina Dresing, University of Turku, Finland
Copyright © 2023 Reames. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Donald V. Reames, dvreames@gmail.com