Issue
A&A
Volume 650, June 2021
Parker Solar Probe: Ushering a new frontier in space exploration
Article Number A25
Number of page(s) 13
Section Catalogs and data
DOI https://doi.org/10.1051/0004-6361/202039833
Published online 02 June 2021

© R. C. Allen et al. 2021

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

Stream interaction regions (SIRs) form in interplanetary space where faster solar wind streams, originating from coronal holes, overtake the preceding slower solar wind streams (e.g., Belcher et al. 1971; Pizzo 1978; Richardson 2018). Where the two streams meet, a density “pileup” of compressed plasma is formed upstream of the stream interface, which is typically identified as the peak in total pressure between the two solar wind populations, and is followed by a rarefaction region in the fast solar wind (see Richardson 2018, and references therein). As the SIR propagates away from the Sun, the compression can lead to the formation of an SIR-associated forward and reverse shock pair at and often beyond 1 au, which efficiently accelerates particles (see Jian et al. 2006, 2008; Richardson 2018, and references therein). As such, SIRs have been found to be a major source of energetic particles in interplanetary space (e.g., Van Hollebeke et al. 1981; Tsurutani et al. 1982; Richardson 2018; Kollmann et al. 2019). As the main source of the high-speed stream, the coronal hole rotates around the solar rotation axis, and the SIR structure rotates with it. After a complete solar rotation, the SIR is classified as a corotating interaction region (CIR). SIRs and CIRs can also be highly geoeffective, triggering geomagnetic storms (e.g., Tsurutani & Gonzalez 1997; Turner et al. 2006) and affecting the ionosphere and thermosphere of Earth (e.g., Chen et al. 2014). Additionally, these structures, and their associated shocks, can modulate the intensity of galactic cosmic rays from the inner heliosphere out to the very local interstellar medium (e.g., Hill et al. 2020a).

The properties of SIRs and CIRs vary with radial distance and latitude, as well as over the lifetime of the SIR and CIR as the footpoint of the structure rotates around the solar rotation axis. Prior to the launch of Helios 1 and 2 in 1974, Hunhausen (1973) constructed a hydrodynamic model of the formation of SIRs in the inner heliosphere. The model predicted that as the fast streams flow away from the Sun, the velocity transition from slow-to-fast wind would steepen, the density pileup upstream of the SIR would grow, and the pressure peak and temperature increase near the interface would both become relatively stronger and sharper (Hunhausen 1973). These modeling results were later confirmed by superposed epoch investigations of SIRs where events between 0.3–0.4 au at 0.9–1 au were compared (Richter & Luttrell 1986). Later comparisons of SIR observations at Pioneer Venus Orbiter (at 0.7 au) and those at 1 au expanded upon these studies, finding that the rate of SIRs with sharp stream interfaces increased from 9% at 0.7 au to 21% at 1 au (Jian 2008).

Regarding the intensity of energetic particles associated with SIRs at different radial distances, Van Hollebeke et al. (1978) compared SIR and CIR events between Helios 1 and 2 (spanning 0.3–1 au), IMP 7 (at 1 au), Pioneer 10 (spanning 9–10 au), and Pioneer 11 (at 3.8 au). These comparisons revealed that the SIR- and CIR-associated suprathermal flux increased in intensity with radial distance from the Sun to ~3 au, where the intensity peaked, before decreasing again at larger radial distances (Van Hollebeke et al. 1978). This motivated the current understanding that most SIR- and CIR-associated shocks form around 3 au (see Richardson 2018, and references therein). Particles accelerated at these shock fronts can propagate back into the inner heliosphere along magnetic field lines, although scattering processes and magnetic cooling decrease the intensity and change the spectral form of the particles as they approach the Sun (e.g., Mewaldt et al. 1978; Fisk & Lee 1980; Mason et al. 1999; Zhao et al. 2016). In addition to radial variations in the ion spectra, the suprathermal ion spectra and intensity can also be different between spacecraft closely separated in longitude, depending on where along the shock surface at greater heliospheric distances the particles are magnetically connected (e.g., Zhao et al. 2015).

Along with particles accelerated further out in the heliosphere, local, non-shock-related acceleration has been proposed to occur in the inner heliosphere associated with the stream interface (e.g., Schwadron et al. 1996; Giacalone et al. 2002; Ebert et al. 2012; Filwett et al. 2017, 2019; Allen et al. 2020a). These studies have suggested that in the inner heliosphere, the lower-energy component of SIR- and CIR-associated suprathermal ions may be a result of local acceleration, while the higher-energy component is a result of distant shock acceleration (e.g., Ebert et al. 2012; Filwett et al. 2019). This has been supported by observations of acceleration at SIR compressions within 1 au, where shocks have not yet formed (Chotoo et al. 2000; Chen et al. 2015). Recent observations by Parker Solar Probe (PSP) at 0.3 au (McComas et al. 2019) have further suggested the importance of compression-related acceleration at lower heliospheric distances than where SIR- and CIR-associated shocks typically form (Joyce et al. 2021). PSP observations have also observed weaker modulation of the energetic particles from transport than expected, suggesting sub-Parker field lines may allow a more direct access to shock regions (Joyce et al. 2021; Schwadron et al. 2021).

In addition to the radial evolution of SIRs and CIRs as well as their associated suprathermal particles, these structures also vary temporally as they corotate with the Sun. The bulk plasma properties and suprathermal ion intensity can vary from rotation to rotation (e.g., Mason et al. 2009; Allen et al. 2021) as the initial stream forms into a well-developed CIR with associated shocks and continues to evolve. More generally, SIR and CIR properties as well as occurrences have a solar cycle dependence (Jian et al. 2011, 2019), as does the SIR and CIR suprathermal ion composition (e.g., Mason et al. 2008, 2012; Allen et al. 2019). The combination of temporal and radial evolution along with longitudinal, latitudinal, and solar cycle variations underscore the importance of multi-point and successive observations between various missions to better study the dynamics and evolution of SIRs and CIRs.

With the addition of PSP to the Heliophysics System Observatory (HSO), we have the ability to investigate SIRs and CIRs within the inner heliosphere for the first time since the Helios era (see McComas et al. 2019; Allen et al. 2020a, 2021; Cohen et al. 2020; Desai et al. 2020; Joyce et al. 2020, 2021). To better understand the variability and evolution of SIRs and CIRs, we present here an initial catalog of SIRs and CIRs observed by PSP during the first five orbits of the mission, along with corresponding observations at the Solar Terrestrial Relations Observatory-Ahead (STEREO-A), Advanced Composition Explorer (ACE), and Wind spacecraft. Future iterations of this catalog will also include corresponding SIR and CIR observations from Solar Orbiter (Müller et al. 2020) as the data become available. Updates to this catalog will correspond to each new public release of PSP data and will be publicly available through the PSP Science Gateway1. This will allow for researchers to not only find fortuitous events, but also facilitate new statistical studies needed to further our understanding of the spatiotemporal evolution of SIRs and CIRs as well as their associated energetic particle populations. Section 2 discusses the different missions and datasets currently used for this catalog and describes the criteria used to identify SIRs and CIRs for inclusion in this list. Section 3 illustrates the PSP observed SIRs and CIRs over the first five orbits along with corresponding observations by STEREO-A and L1 missions (ACE and Wind). A detailed look at two of these identified SIR and CIR events and preliminary statistical analysis of the radial evolution of SIR and CIR properties are also discussed. Finally, a summary of this catalog, along with a discussion of how it could be used for future investigations, is given in Sect. 4.

2 Methodology

2.1 Parker Solar Probe

The PSP mission (Fox et al. 2016) launched into orbit around the Sun on 2018 August 12 and is providing key measurements of the solar wind and interplanetary magnetic field in the inner heliosphere. Over the lifetime of the mission, the perihelion of PSP has gradually lowered through the use of Venus gravity assists. The observations shown here are 1-min moment averages of bulk solar wind measurements from the Solar Probe Cup (SPC, Case et al. 2020), which is part of the Solar Wind Electrons Alpha and Protons(SWEAP) instrument suite (Kasper et al. 2016), along with 1-min averages of magnetic field measurements from the FIELDS suite (Bale et al. 2016). To identify SIRs and CIRs, the solar wind tangential flow deflection is computed using the equation αtan = tan−1(vTvR), where vT and vR are the tangential and radial velocity components, respectively, in the radial-tangential-normal (RTN) coordinate system. TheSPC proton thermal velocity data product is used to obtain a scalar proton temperature through (1)

where T is the temperature, m is the mass of a proton, kB is Boltzmann’s constant, and vtherm is the thermal velocity. The specific entropy argument (S) of the solar wind is computed from the SPC proton density (n) and the derived temperature using the formulation (e.g., Jian et al. 2006), (2)

Additionally, this study uses a combined proton plasma and magnetic pressure (P) derived from the FIELDS magnetic field magnitude (B) and the SPC proton density and derived temperature, (3)

To investigate the suprathermal ion abundance during the SIR and CIR events, 30-min averaged measurements from the Energetic Particle Instrument-Lo (EPI-Lo, Hill et al. 2017, 2020b), which is part of the Integrated Science Investigation of the Sun (IS⊙IS, McComas et al. 2016), are used. This paper presents only one EPI-Lo data product, the “IonToF” product, which is a time-of-flight only measurement that allows for observations of ions from ~30–586 keV nuc−1, but does not distinguish between separate ion species.

thumbnail Fig. 1

Overview of the first five orbits of PSP showing: (a) heliocentric distance, (b) radial H+ velocity, (c) H+ density, (d) H+ temperature, (e) specific entropy argument, (f) sum of H+ plasma and magnetic pressures, (g) magnetic field magnitude, and (h) suprathermal ion count rate. Blue and red dashed lines indicate the times of clean and ambiguous SIR and CIR observations, respectively.

2.2 STEREO-A

Since the aim of this catalog is to highlight corresponding observations between PSP and other spacecraft in the HSO, observationsfrom the STEREO mission (Kaiser et al. 2008), a set of two spacecraft that were launched on 2006 October 25, are used. However, only observations from STEREO-A (STA) are available in the PSP era. Similar to PSP, 1-min averaged magnetic field measurements from the magnetometer (Acuña et al. 2008) and suprathermal particles from the Solar Electron and Proton Telescope (SEPT, Müller-Mellin et al. 2008) instruments on the In situ Measurements of Particles and CME Transients (IMPACT) investigation (Luhmann et al. 2008) are used to identify SIRs and CIRs. SEPT measures ions in the energy range of 60–7000 keV, but it cannot differentiate the ion species. Additionally, 1-min bulk solar wind properties (i.e., velocity, temperature, and density) from the PLAsma and SupraThermal Ion Composition investigation (PLASTIC, Galvin et al. 2008) are used. Calculations of T, S, and P following the same procedures as for PSP are also performed.

2.3 L1 Missions: wind and ACE

In addition to STA, this catalog links SIR and CIR observations at PSP to those from the Wind mission (Acuña et al. 1995). Wind, which launched in 1994, has been stationed near the Sun-Earth Lagrange point 1 (L1) since2004. The SIR and CIR search and catalog entries are based on 1-min resolution bulk solar wind plasma properties from the Solar Wind Experiment (SWE, Ogilvie et al. 1995) and magnetic field measurements from the Magnetic Fields Investigation (MFI, Lepping et al. 1995). As is done with STA measurements, T, S, and P are calculated following the same procedures as for PSP.

The last observatory currently included in the catalog is the ACE mission (Stone et al. 1998). ACE, which launched in 1997, is stationed near L1 and is used in tandem with Wind since their observations are nearly identical. The combination of Wind and ACE allows for a more continuous dataset. Solar wind bulk properties (velocity, density, and temperature) are provided by the Solar Wind Electron Proton Alpha Monitor (SWEPAM, McComas et al. 1998). ACE magnetic field measurements are provided by the Magnetic Field Experiment (MAG, Smith et al. 1998), and suprathermal ion observations are from the Electron, Proton, and Alpha Monitor-Low Energy Magnetic Spectrometer (EPAM-LEMS, Gold et al. 1998) instrument. Our analysis uses 64 s averages of SWEPAM observations, 1-min averages of MAG observations, and 5-min ion intensity averages from EPAM/LEMS120.

Table 1

PSP SIR and CIR events.

Table 2

Corresponding STA SIR and CIR events.

thumbnail Fig. 2

Overview of STA observations corresponding to SIRs and CIRs observed at PSP showing: (a) radial H+ velocity, (b) H+ density, (c) H+ temperature, (d) specific entropy argument, (e) sum of H+ plasma and magnetic pressures, (f) magnetic magnitude, and (g) suprathermal ion intensities for energies between 37.5 and 3000 keV. Vertical blue dashed lines illustrate the times of SIR and CIR observations that were also observed by PSP.

2.4 Catalog criteria

Following the same criteria as Jian et al. (2006, 2013), PSP-observed SIR and CIR events were selected that matched the majority of the following criteria given the available data: (1) velocity increase, (2) total pressure peak observed during the velocity increase, (3) enhancement of density ahead of the interface, (4) local maxima in the magnetic field magnitude during the velocity increase, (5) temperature and entropy increase from the slow to fast solar wind, and (6) a tangential flow deflection around the interface. These events are then categorized into one of the following three subsets identified with quality flags: clean events without data gaps with at least five of the six criteria satisfied (denoted by a value of zero); clean events(i.e., at least five of six criteria satisfied) with data available through the interface region but with gaps in the slow and/or fast solar wind near the velocity increase (value of one); and events missing some desired data with larger data gaps occurring during the velocity increase (value of two), which still satisfy criteria one, three, and five. Examples of these are shown in the following section. It should be noted that these SIR and CIR events are identified with only observations from SWEAP and FIELDS, and as such data are only available for a subset of the orbit. There are intervals, particularly when PSP is further out in its orbit closer to aphelion, where IS⊙IS data are available without corresponding observations from FIELDS and SWEAP (see events in Cohen et al. 2020; Desai et al. 2020). As a result, it is likely that SIR and CIR observations within the IS⊙IS dataset are missing from this catalog. Additionally, to focus on clean SIR and CIR observations, events that were determined to occur near interplanetary coronal mass ejections (ICMEs) are not included in this list. As this catalog uses more stringent selection criteria, some of the high-speed streams from Allen et al. (2020a) are not included in this list.

After identifying an SIR and CIR observed by PSP, the structure is tracked using the approach of Allen et al. (2020a). First, the structure is identified in a global 3D magnetohydrodynamic heliospheric Wang-Sheeley-Arge (WSA)-ENLIL model (Arge et al. 2004; Odstrcil et al. 2004) using synoptic solar magnetic field maps derived from magnetograms. For the purpose of building this catalog, zero-point-corrected quick-reduce magnetogram synoptic maps from the Global Oscillation Network Group (GONG, Harvey et al. 1996) were used. After identifying the potential SIR and CIR structure in the WSA-ENLIL model2, the structure was followed both forward and backward in time until it would reach either STA or L1. This coarse timing approximation is then used as a guide in searching the STA or L1 datasets for observations of the structure. If the observational characteristics suggest that an SIR and CIR, which was observed by PSP, was seen by either STA or L1, then the 1 au observation time is included in the list. This process is then repeated until the SIR or CIR is not observed, allowing for the creation of the catalog of SIR and CIR observations that correspond to those seen by PSP. We note that there are other SIR and CIR events that have been listed by others (e.g., Jian et al. 2019) which were not observed by PSP, and as such are not included here. Future work will expand on this catalog of SIRs and CIRs as the PSP mission progresses. Additionally, this same process will be expanded to include Solar Orbiter as those data become available.

thumbnail Fig. 3

Overview of L1 observations (ACE in red and Wind in black for panels a–f) in the same format as Fig. 2. Panel g: suprathermal ion intensities for the energy range of 56–3020 keV from ACE-EPAM. Vertical blue dashed lines illustrate the times of SIR and CIR observations that were also observed by PSP.

3 Results

3.1 PSP SIRs and CIRs

In situ observations from PSP during the first five orbits are shown in Fig. 1. During these orbits, 11 SIRs and CIRs were observed that met our above criteria, and they are denoted by color coded vertical lines (Fig. 1). Data gaps often occur when PSP is at larger heliospheric distances, especially for SPC, as well as during Ka-band downlinking intervals, and so this limits the time periods that are searchable for SIRs and CIRs. Following the methodology listed in Sect. 2.4, no SIRs or CIRs were observed during orbit four. This is likely due to PSP traveling fast enough around the Sun so that it likely remained between high-speed streams emanating from coronal holes during this time period (not shown). Apparent velocity enhancements that are not flagged as SIR and CIR events in Fig. 1 are found to correspond to ICMEs, sector boundary crossings, or slow velocity increases without the formation of a stream interface.

Table 1 lists the eleven SIRs and CIRs observed at PSP during the first five orbits. The event numbers listed in the first column are used for later reference with the identified corresponding events. The SIR and CIR start and end times are identified by eye based on the plasma and fields measurements. Events that occur near data gaps may not have fully accurate start and/or stop times identified. The interface time is defined as the time of peak pressure during this velocity ramp. Additional information (radial distance at the start time, max P, max B, max n, min V, and max V) are also provided. The P and B maxima are defined as the peak value within the velocity increase, while the maximum of n is defined as the peak density within two days of the beginning of the velocity increase and two days following the velocity increase. All quantities are validated by eye to ensure the values picked are not outliers and that the peak in density corresponds to the density pileup region. The min V and max V are the minimum and maximum velocity in between the SIR and CIR start and end times identified. When data gaps are present during the velocity increase, the interface time and max P, B, and n are all set to “-”, since these quantities cannot be reliably identified. An “ambiguous” flag is assigned using the criteria discussed in Sect. 2.4. Table 1 indicates if each event has a corresponding observation at STA and/or L1, as well as if the structure was a CIR or not, that is to say if it had been observed to have persisted for a complete solar rotation. The determination of whether the structure is an SIR or CIR is done using corresponding observations. If the structure can be successfully tracked back in time by at least one solar rotation then it is considered a CIR, and, ifnot, it is deemed an SIR. Finally, the rightmost columns of Table 1 list whether or not the structure was observed with a shock and if it was observed near a heliospheric current sheet (HCS) crossing

Table 3

Corresponding L1 SIR and CIR events.

3.2 Corresponding SIRs and CIRs

To probe the evolution of SIRs and CIRs in both space (i.e., radial evolution) and time, each PSP SIR and CIR observation is tracked both backward and forward in time in a search for a corresponding observation at STA and L1 (ACE and Wind). Figure 2 illustrates the STA observations over the same time period as the PSP observationsin Fig. 1, with vertical lines denoting the times of SIR and CIR observations corresponding to SIRs andCIRs identified at PSP. While the speed at which CIRs corotate has some variability (e.g., Allen et al. 2020b), the near-circular orbit of STA results in CIR observations to appear with a repeatable cadence approximately every 25 days. It is also clear that the elliptical orbit of PSP, along with its super-corotating speeds near perigee, can lead PSP to miss some SIR and CIR structures that are observed at STA. As such, the list of corresponding STA observations should not be taken as an exhaustive list of STA-observed 1 au SIRs and CIRs as has been compiled by others (e.g., Jian et al. 2019). Table 2 lists the corresponding SIRs and CIRs identified in the STA data. As in Table 1, the first 11 columns mark the event times and characteristics. The following columns denote the PSP event number each STA observation corresponds to, whether or not there were also corresponding observations at L1, if the event at STA had been observed for at least one solar rotation prior to this observation, if the SIR and CIR had an associated shock, and if it was observed near a HCS crossing. If a CIR persists long enough, such as for STA events 1, 3, 5, and 7, then there can be multiple STA observations corresponding to multiple observations at PSP (PSP events 1 and 2 in this example).

Observations from ACE and Wind are shown in Fig. 3 in the same fashion as Fig. 2 (with ACE measurements in red and Wind measurements in black). As was the case for the STA events, the list of corresponding SIRs and CIRs identified in the ACE and Wind measurements should not be constructed as an exhaustive list at L1 as done by others (e.g., Jian et al. 2006, 2011). Table 3 provides the catalog of corresponding SIRs and CIRs observations at L1 in the same fashion as Table 2.

thumbnail Fig. 4

Example of a clean SIR and CIR event observed by PSP showing: (a) radial H+ velocity, (b) tangential H+ velocity deflection, (c) H+ density, (d) H+ temperature, (e) specific entropy argument, (f) sum of H+ plasma and magnetic pressure, (g) magnetic field magnitude, (h) magnetic ΘRTN, (i) magnetic ΦRTN, and (j) 30–586 keV suprathermal ion count rate. The blue and orange vertical dashed lines indicate where the radial velocity transitions from slow to fast solar wind. The red dotted line denotes the pressure peak within the velocity increase.

thumbnail Fig. 5

Nearest corresponding observations at L1 (left) and STEREO-A (right) of the SIR and CIR observed by PSP in Fig. 4. Panels show (a and j) radial H+ velocity, (b and k) H+ density, (c and l) H+ temperature, (d and m) specific entropy argument, (e and n) sum of H+ plasma and magnetic pressure, (f and o) magnetic field magnitude, (g and p) magnetic ΘRTN, (h and q) magnetic ΦRTN, and (i and r) suprathermal ion intensities. Vertical lines indicate the start of the velocity increase (blue dashed), peak pressure during the velocity increase (red dotted), and the end of the velocity increase (orange dashed).

3.3 Example of a clean event

An example of a “clean” SIR and CIR event observed by PSP is shown in Fig. 4. This event exhibited all the classical CIR characteristics and is discussed in more detail in other studies (Allen et al. 2020a; Joyce et al. 2021; Schwadron et al. 2021). The velocity increase is observed near 00:00 UT on 2018 November 15. The verticalblue and orange dashed lines denote the “start” and “stop” times listed in Table 1, respectively,while the vertical dotted red line indicates the interface (max P between the start and stop times). A clear tangential flow deflection occurred during the velocity ramp (Fig. 4b). Additionally, the solar wind density is observed to have a pileup region upstream of the interface, followed by a lower density region downstreamin the fast-solar wind (Fig. 4c). The temperature and entropy are both higher in the high-speed stream as compared to the slow solar wind (Figs. 4d and e). The pressure and magnetic field magnitude both peak at the interface (Figs. 4f and g). For this event at PSP, the suprathermal ion intensity has an isolated peak at the SIR and CIR interface (Fig. 4j), as well as a broader enhancement on 2018 November 16; these suprathermal ion enhancements are further discussed in Allen et al. (2020a), and the ion spectra are investigated in detail by Joyce et al. (2021). Since all of the criteria for SIR and CIR selection are satisfied, and there are no data gaps for this event, this SIR and CIR is deemed a clean event (i.e., ambiguous flag of 0). Additionally, this structure was also observed by both STA and L1 (see Tables 2 and 3, respectively), as well as by PSP at an earlier time (event 1 in Table 1). From the preceding observations at 1 au and at PSP, this event is identified as a CIR since the first observation of this structure was 98 days prior to the PSP observation on 2018 November 14–15 (i.e., more than one solar rotation).

Figure 5 illustrates the two temporally closest corresponding observations from L1 and STA. For this solar rotation, the CIR structure was first observed at L1 on 2018 November 4–5, and last observed at STA on 2018 November 22–23. Both L1 and STA detect a velocity enhancement, density pileup, temperature and entropy increase, and peak of P and B at the stream interface. However, notably, the suprathermal ion measurements at L1 do not show a large enhancement, while STA detects a clear enhancement following the interface and stretching into the high-speed stream. The differences between the three observations (L1, PSP, and STA) illustrate the importance of multipoint measurements to better understand the evolution of large-scale structures in the heliosphere. As discussed in Allen et al. (2020a), the differences in the suprathermal ion populations at and near the stream interface may be due to the relative contribution of locally-accelerated populations versus the inward transport of shock-accelerated ions (see also Schwadron et al. 1996; Filwett et al. 2019). Through investigating the PSP SIR and CIR structures in the context of 1 au measurements, insight into both the radial and temporal evolution can be gathered.

thumbnail Fig. 6

Example of an ambiguous SIR and CIR event observed by PSP in the same format as Fig. 4.

thumbnail Fig. 7

Closest corresponding observations of the SIR and CIR observed by PSP in Fig. 6 at STEREO-A (left) and L1 (right) shown in the same format as Fig. 5.

3.4 Example of an ambiguous event

An example of an ambiguous SIR and CIR event observed by PSP is shown in Fig. 6. This event has previously been studied in depth by Allen et al. (2021). Due to data gaps caused by Ka-band downlinking from the spacecraft, only fragments of this interval were collected. From the intervals with data, however, the beginning of the velocity increase is observed along with a tangential flow deflection, density pileup, and enhancement in temperature and entropy from slow to fast wind. As most of the velocity increase is missing, it is not possible to identify the interface, which would have occurred sometime between 18:00 UT on 2019 September 19 to 12:00 UT on 2019 September 20 (Fig. 6a). However, the suprathermal ions did exhibit an enhancement during the start of the high-speed stream (Fig. 6j). Given the available data, this event is categorized as an SIR and CIR, and so it is included in our catalog of events (Table 1). However, due to data gaps occurring during the velocity enhancement, this event is flagged as ambiguous (ambiguous counter of 2). It should also be noted that this structure persisted for multiple solar rotations as observed at 1 au (Tables 2 and 3) and as such is identified as a CIR. While these observations from PSP are not able to provide a full picture of the CIR, such as information of the peak pressure at the interface, these types of observations can be nonetheless scientifically valuable. As an example, for the event shown in Fig. 6 (Event 10 in Table 1), the CIR was observed when PSP and STA were nearly radially aligned, enabling an examination of the radial variation in the CIR properties that were measured by PSP.

Figure 7 shows the most recent L1 and STA observations corresponding to the PSP event shown in Fig. 6 (all observations of this CIR are discussed in Allen et al. 2021). The CIR was observed at STA approximately 1.77 days after the PSP observations and at L1 ~7.5 days after the PSP observation. All of the signatures of the CIR are seen to be very similar between the three observations, indicating that, in contrast to the event in Figs. 4 and 5, the structure observed in late September 2019 is not evolving much between observations by the different spacecraft. This static nature of the mature CIR shown in Figs. 6 and 7 allows for investigations of the radial variations of CIRs and their associated energetic particles (see Allen et al. 2021).

3.5 Initial statistical trends

To investigate the evolution of SIR and CIR properties as a function of heliocentric distance, the maximum density, peak pressure, and peak magnetic field magnitude for the events are scaled by radial distance to 1 au. Following Perrone et al. (2019), the density is scaled as nR−1.96 and magnetic field magnitude as BR−1.59. The scaled pressure was computed by first scaling the temperature using TR−0.93 (Perrone et al. 2019) as well as the scaled density and magnetic field magnitude to recompute the pressure using Eq. (3). Due to data gaps in ambiguous events potentially leading to a misidentification of the max density and peak magnetic field magnitude and pressure, only clean PSP events are used for this initial statistical study. Scaling the observations by radial distance allows for an investigation of radial variations due to the interaction of the fast and slow speed streams (i.e., this scaling removes the underlying systematic radial variations due to processes such as volumetric expansion that would otherwise dominate the radial trends).

Figure 8 shows the radial evolution of the scaled max density, magnetic field magnitude, and pressure. While there is considerable variation in the SIR and CIR properties at 1 au, they are often higher than the six clean PSP events.These results are consistent with an earlier interpretation of Helios observations (e.g., Richter & Luttrell 1986), where the SIR and CIR structures develop the density pile-up and compression of the magnetic field, resulting in a pressure enhancement that grows with radial distance. However, because of the limited number of clean events in the inner heliosphere, as well as the spread in values at 1 au and known temporal variations in the SIRs and CIRs, additional observations are required before statistically meaningful relations can be extrapolated from these observations. These early statistical trends do indicate that future investigations using PSP, STEREO-A, L1 observatories, and eventually Solar Orbiter will allow for further investigations in the radial evolution of SIRs and CIRs.

thumbnail Fig. 8

Radial profile of (a) adjusted max density, (b) adjusted peak magnetic field magnitude, and (c) adjusted peak pressure. Values were adjusted to remove bulk radial variations following the radial relations in Perrone et al. (2019).

4 Summary

SIRs and CIRs are large-scale, and often long-lived structures that are an important aspect of the solar wind and dynamics in interplanetary space. Understanding how these structures evolve in the inner heliosphere as they form and persist for multiple solar rotations requires a robust database of observations at small heliospheric distances and complementary measurements from other observatories at different locations and successive points in time. Such a living catalog of SIR and CIR events observed by PSP with corresponding observations at STA and L1 (ACE and Wind), and the criteria for inclusion in the catalog, are presented here. Future iterations of the catalog will also include the identification of corresponding events at Solar Orbiter. These catalogs of SIRs and CIRs will be updated and published on the PSP Science Gateway3 concurrent to PSP public data releases.

While adhering to stringent identification criteria, the catalog of PSP events includes intervals where measurements are not fully available for complete characterization of the stream interface and parts of the upstream or downstream regions. However, as demonstrated by the two different PSP events highlighted here, certain features of the SIR and CIR structures can still be analyzed in depth. As also illustrated by the two examples of SIR and CIR events shown here, corresponding observations, even for incomplete measurement intervals, allow for studies of both temporal (longitudinal) and radial evolution of SIRs and CIRs between observatories in a way that has not been possible since the Helios era. This is further demonstrated by the initial analysis of the radial variation of certain SIR and CIR properties also included here, which shows evidence of systematic increases in bulk characteristics at the structure interface with increasing heliocentric distance. Additionally, the advanced instrumentation of modern missions, along with latitudinal separation provided by Solar Orbiter, will allow for more detailed investigations of SIR and CIR structures and associated energetic plasma populations.

Acknowledgements

This work was supported by NASA’s Parker Solar Probe Mission, contract NNN06AA01C and by NASA’s ACE grant NNX17AC05G. C.M. thanks the Austrian Science Fund (FWF): P31659-N27. The ISOIS data and visualization tools are availableto the community at: https://spacephysics.princeton.edu/missions-instruments/isois; data are also available via the NASA Space Physics Data Facility (https://spdf.gsfc.nasa.gov/) and from https://sppgway.jhuapl.edu/event_list. The STEREO SEPT data are available at http://www2.physik.uni-kiel.de/stereo/data/sept/, and STEREO magnetic field and plasma datacan be found at the STEREO Science Center at https://stereo-ssc.nascom.nasa.gov. Wind data are available at https://cdaweb.gsfc.nasa.gov/. The ACE data are publicly available at the ACE mission website at http://www.srl.caltech.edu/ACE/. WSA-ENLIL movies of PSP orbits are available at http://helioweather.net/missions/psp/. Living versions of the catalog presented in this paper are available and kept up to date on both the Parker Solar Probe Science Gateway (https://sppgway.jhuapl.edu/event_list/).

References

  1. Acuña, M. H., Ogilvie, K. W., Baker, D. N., et al. 1995, Space Sci. Rev., 71, 5 [Google Scholar]
  2. Acuña, M. H., Curtis, D., Scheifele, J. L., et al. 2008, Space Sci. Rev., 136, 203 [Google Scholar]
  3. Allen, R. C., Ho, G. C., & Mason, G. M. 2019, ApJ, 883, L10 [Google Scholar]
  4. Allen, R. C., Lario, D., Odstrcil, D., et al. 2020a, ApJS, 246, 36 [Google Scholar]
  5. Allen, R. C., Ho, G. C., Jian, L. K., et al. 2020b, Space Weather, 18, 4 [Google Scholar]
  6. Allen, R. C., Ho, G. C., Mason, G. M., et al. 2021, Geophys. Res. Lett., 48, e2020GL091376 [Google Scholar]
  7. Arge, C. N., Luhmann, J. G., Odstrcil, D., Schrijver, C. J., & Li, Y. 2004, J. Atmos. Sol. Terr. Phys., 66, 1295 [Google Scholar]
  8. Bale, S. D., Goetz, K., Harvey, P. R., et al. 2016, Space Sci. Rev., 204, 49 [Google Scholar]
  9. Belcher, J. W., & Davis, L. 1971, J. Geophys. Res., 76, 3534 [Google Scholar]
  10. Case, A. W., Kasper, J. C., Stevens, M. L., et al. 2020, ApJS, 246, 43 [Google Scholar]
  11. Chen, G.-M., Xu, J., Wang, W., & Burns, A. G. 2014, J. Geophys. Res. Space Phys., 119, 7928 [Google Scholar]
  12. Chen, J. H., Schwadron, N. A., Möbius, E., & Gorby, M. 2015, J. Geophys. Res. Space Phys., 120, 9269 [Google Scholar]
  13. Chotoo, K., Schwadron, N. A., Mason, G. M., et al. 2000, J. Geophys. Res., 105, 23107 [Google Scholar]
  14. Cohen, C. M. S., Christian, E. R., Cummings, A. C., et al. 2020, ApJS, 246, 20 [Google Scholar]
  15. Desai, M. I., Mitchell, D. G., Szalay, J. R., et al. 2020, ApJS, 246, 56 [Google Scholar]
  16. Ebert, R. W., Dayeh, M. A., Desai, M. I., & Mason, G. M. 2012, ApJ, 749, 73 [Google Scholar]
  17. Filwett, R. J., Desai, M. I., Dayeh, M. A., & Broiles, T. W. 2017, ApJ, 838, 23 [Google Scholar]
  18. Filwett, R. J., Desai, M. I., Ebert, R. W., & Dayeh, M. A. 2019, ApJ, 876, 88 [Google Scholar]
  19. Fisk, L. A., & Lee, M. A. 1980, ApJ, 237, 620 [Google Scholar]
  20. Fox, N. J., Velli, M. C., Bale, S. D., et al. 2016, Space Sci. Rev., 204, 7 [NASA ADS] [CrossRef] [Google Scholar]
  21. Galvin, A. B., Kistler, L. M., Popecki, M. A., et al. 2008, Space Sci. Rev., 136, 437 [Google Scholar]
  22. Giacalone, J., Jokipii, J. R., & Kóta, J. 2002, ApJ, 573, 845 [Google Scholar]
  23. Gold, R. E., Krimigis, S. M., Hawkins, S. E. III, et al. 1998, Space Sci Rev., 86, 541 [Google Scholar]
  24. Harvey, J. W., Hill, F., Hubbard, R. P., et al. 1996, Science, 272, 1284 [Google Scholar]
  25. Hill, M. E., Mitchell, D. G., Andrews, G. B., et al. 2017, J. Geophys. Res. Space Phys., 122, 1513 [Google Scholar]
  26. Hill, M. E., Allen, R. C., Kollmann, P., et al. 2020a, ApJ, 905 69 [Google Scholar]
  27. Hill, M. E., Mitchell, D. G., Allen, R. C., et al. 2020b, ApJS, 246, 65 [Google Scholar]
  28. Hundhausen, A. J. 1973, J. Geophys. Res., 78, 1528 [CrossRef] [Google Scholar]
  29. Jian, L. 2008, PhD thesis, University of California, Los Angeles, USA [Google Scholar]
  30. Jian, L., Russell, C. T., Luhmann, J. G., & Skoug, R. M. 2006, Sol. Phys., 239, 337 [NASA ADS] [CrossRef] [Google Scholar]
  31. Jian, L. K., Russell, C. T., Luhmann, J. G., Skoug, R. M., & Steinberg, J. T. 2008, Sol. Phys., 250, 375 [Google Scholar]
  32. Jian, L. K., Russell, C. T., & Luhmann, J. G. 2011, Sol. Phys., 274, 321 [Google Scholar]
  33. Jian, L. K., Russell, C. T., Luhmann, J. G., Galvin, A. B., & Simunac, K. D. 2013, AIP Conf. Proc., 1539, 191 [Google Scholar]
  34. Jian, L. K., Luhmann, J. G., Russell, C. T., & Galvin, A. B. 2019, Sol. Phys., 294, 31 [Google Scholar]
  35. Joyce, C. J., McComas, D. J., Christian, E. R., et al. 2020, ApJS, 246, 41 [Google Scholar]
  36. Joyce, C. J., McComas, D. J., Schwadron, N. A., et al. 2021, A&A, 650, L5 (PSP SI) [EDP Sciences] [Google Scholar]
  37. Kaiser, M. K, Kucera, T. A., Davila, J. M., et al. 2008, Space Sci. Rev., 136, 5 [Google Scholar]
  38. Kasper, J. C., Abiad, R., Austin, G., et al. 2016, Space Sci. Rev., 204, 131 [NASA ADS] [CrossRef] [Google Scholar]
  39. Kollmann, P., Hill, M. E., McNutt, R. L., Jr., et al. 2019, ApJ, 876, 46 [Google Scholar]
  40. Lepping, R. P., Acuna, M. H., Burlaga, L. F., et al. 1995, Space Sci. Rev., 71, 207 [Google Scholar]
  41. Luhmann, J. G., Curtis, D. W., Schroeder, P., et al. 2008, Space Sci. Rev., 136, 117 [Google Scholar]
  42. Mason, G. M., Von Steiger, R., Decker, R. B., et al. 1999, Corotating Interaction Regions, eds. A. Balogh, J. T. Gosling, J. R. Jokipii, R. Kallenbach, & H. Kunow (Dordrecht: Springer Netherlands), 327 [Google Scholar]
  43. Mason, G. M., Leske, R. A., Desai, M. I., et al. 2008, ApJ, 678, 1458 [Google Scholar]
  44. Mason, G. M., Desai, M. I., Mall, U., et al. 2009, Sol. Phys., 256, 393 [Google Scholar]
  45. Mason, G. M., Desai, M. I., & Li, G. 2012, ApJ, 748, L31 [Google Scholar]
  46. McComas, D. J., Bame, S. J., Barker, P., et al. 1998, Space Sci Rev., 86, 563 [Google Scholar]
  47. McComas, D. J., Alexander, N., Angold, N., et al. 2016, Space Sci. Rev., 204, 187 [Google Scholar]
  48. McComas, D. J., Christian, E. R., Cohen, C. M. S., et al. 2019, Nature, 576, 223 [Google Scholar]
  49. Mewaldt, R. A., Stone, E. C., & Vogt, R. E. 1978, Geophys. Res. Lett., 5, 965 [Google Scholar]
  50. Müller, D., St Cyr, O. C., Zouganelis, I., et al. 2020, A&A, 642, A1 [CrossRef] [EDP Sciences] [Google Scholar]
  51. Müller-Mellin, R., Böttcher, S., Falenski, J., et al. 2008, Space Sci. Rev., 136, 363 [Google Scholar]
  52. Odstrcil, D., Riley, P., & Zhao, X. P. 2004, J. Geophys. Res., 109, A02116 [Google Scholar]
  53. Ogilvie, K. W., Chornay, D. J., Frizenreiter, R. J., et al. 1995, Space Sci. Rev., 71, 55 [Google Scholar]
  54. Perrone, D., Stansby, D., Horbury, T. S., & Matteini, L. 2019, MNRAS, 483, 3730 [Google Scholar]
  55. Pizzo, V. 1978, J. Geophys. Res., 83, 5563 [Google Scholar]
  56. Richardson, I. G. 2018, Liv. Rev. Sol. Phys., 15, 1 [Google Scholar]
  57. Richter, A. K., & Luttrell, A. H. 1986, J. Geophys. Res., 91, 5873 [NASA ADS] [CrossRef] [Google Scholar]
  58. Schwadron, N. A., Fisk, L. A., & Gloackler, G. 1996, Geophys. Res. Lett., 23, 2871 [Google Scholar]
  59. Schwadron, N. A., Joyce, C. J., Aly, A., et al. 2021, A&A, 650, A24 (PSP SI) [EDP Sciences] [Google Scholar]
  60. Smith, C. W., L’Heureux, J., Ness, N. F., et al. 1998, Space Sci. Rev., 86, 613 [Google Scholar]
  61. Stone, E. C., Frandsen, A. M., Mewaldt, R. A., et al. 1998, Space Sci. Rev., 86, 1 [Google Scholar]
  62. Tsurutani, B. T., & Gonzalez, W. D. 1997, Geophys. Monogr. Ser., 98, 77 [Google Scholar]
  63. Tsurutani, B. T., Smith, E. J., Pyle, K. R., & Simpson, J. A. 1982, J. Geophys. Res., 87, 7389 [Google Scholar]
  64. Turner, N. E., Mitchell, E. J., Knipp, D. J., & Emery, B. A. 2006, Geophys. Monogr., 167, 113 [Google Scholar]
  65. Van Hollebeke, M. A. I., McDonald, F. B., Trainor, J. H., & von Rosenvinge, T. T. 1978, J. Geophys. Res., 83, 4723 [Google Scholar]
  66. Van Hollebeke, M. A. I, McDonald, F. B., Trainor, J. H., & von Rosenvinge, T. T. 1981, Solar Wind 4, Proceedings of the Conference held in August 18-September 1, 1978 in Burghausen, FDR, ed. H. Rosenbauer. MPAE-W-100-81-31. Garching, FDR: Max-Planck-Institute für Aeronomie, 497 [Google Scholar]
  67. Zhao, L., Li, G., Ebert, R. W., et al. 2015, J. Geophys. Res. Space Phys., 121, 77 [Google Scholar]
  68. Zhao, L., Zhang, M., & Rassoul, H. K. 2016, ApJ, 821, 62 [Google Scholar]

All Tables

Table 1

PSP SIR and CIR events.

Table 2

Corresponding STA SIR and CIR events.

Table 3

Corresponding L1 SIR and CIR events.

All Figures

thumbnail Fig. 1

Overview of the first five orbits of PSP showing: (a) heliocentric distance, (b) radial H+ velocity, (c) H+ density, (d) H+ temperature, (e) specific entropy argument, (f) sum of H+ plasma and magnetic pressures, (g) magnetic field magnitude, and (h) suprathermal ion count rate. Blue and red dashed lines indicate the times of clean and ambiguous SIR and CIR observations, respectively.

In the text
thumbnail Fig. 2

Overview of STA observations corresponding to SIRs and CIRs observed at PSP showing: (a) radial H+ velocity, (b) H+ density, (c) H+ temperature, (d) specific entropy argument, (e) sum of H+ plasma and magnetic pressures, (f) magnetic magnitude, and (g) suprathermal ion intensities for energies between 37.5 and 3000 keV. Vertical blue dashed lines illustrate the times of SIR and CIR observations that were also observed by PSP.

In the text
thumbnail Fig. 3

Overview of L1 observations (ACE in red and Wind in black for panels a–f) in the same format as Fig. 2. Panel g: suprathermal ion intensities for the energy range of 56–3020 keV from ACE-EPAM. Vertical blue dashed lines illustrate the times of SIR and CIR observations that were also observed by PSP.

In the text
thumbnail Fig. 4

Example of a clean SIR and CIR event observed by PSP showing: (a) radial H+ velocity, (b) tangential H+ velocity deflection, (c) H+ density, (d) H+ temperature, (e) specific entropy argument, (f) sum of H+ plasma and magnetic pressure, (g) magnetic field magnitude, (h) magnetic ΘRTN, (i) magnetic ΦRTN, and (j) 30–586 keV suprathermal ion count rate. The blue and orange vertical dashed lines indicate where the radial velocity transitions from slow to fast solar wind. The red dotted line denotes the pressure peak within the velocity increase.

In the text
thumbnail Fig. 5

Nearest corresponding observations at L1 (left) and STEREO-A (right) of the SIR and CIR observed by PSP in Fig. 4. Panels show (a and j) radial H+ velocity, (b and k) H+ density, (c and l) H+ temperature, (d and m) specific entropy argument, (e and n) sum of H+ plasma and magnetic pressure, (f and o) magnetic field magnitude, (g and p) magnetic ΘRTN, (h and q) magnetic ΦRTN, and (i and r) suprathermal ion intensities. Vertical lines indicate the start of the velocity increase (blue dashed), peak pressure during the velocity increase (red dotted), and the end of the velocity increase (orange dashed).

In the text
thumbnail Fig. 6

Example of an ambiguous SIR and CIR event observed by PSP in the same format as Fig. 4.

In the text
thumbnail Fig. 7

Closest corresponding observations of the SIR and CIR observed by PSP in Fig. 6 at STEREO-A (left) and L1 (right) shown in the same format as Fig. 5.

In the text
thumbnail Fig. 8

Radial profile of (a) adjusted max density, (b) adjusted peak magnetic field magnitude, and (c) adjusted peak pressure. Values were adjusted to remove bulk radial variations following the radial relations in Perrone et al. (2019).

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.