Diagnosing solar wind origins using in-situ measurements in the inner heliosphere

Robustly identifying the solar sources of individual packets of solar wind measured in interplanetary space remains an open problem. We set out to see if this problem is easier to tackle using solar wind measurements closer to the Sun than 1 AU, where the mixing and dynamical interaction of different solar wind streams is reduced. Using measurements from the Helios mission, we examined how the proton core temperature anisotropy and cross helicity varied with distance. At 0.3 AU there are two clearly separated anisotropic and isotropic populations of solar wind, that are not distinguishable at 1 AU. The anisotropic population is always Alfv\'enic and spans a wide range of speeds. In contrast the isotropic population has slow speeds, and contains a mix of Alfv\'enic wind with constant mass fluxes, and non-Alfv\'enic wind with large and highly varying mass fluxes. We split the in-situ measurements into three categories according these observations, and suggest that these categories correspond to wind that originated in the core of coronal holes, in or near active regions or the edges of coronal holes, and as small transients form streamers or pseudostreamers. Although our method by itself is simplistic, it provides a new tool that can be used in combination with other methods for identifying the sources of solar wind measured by Parker Solar Probe and Solar Orbiter.


I N T RO D U C T I O N
The solar wind is a continuous flow of plasma, travelling away from the Sun to fill interplanetary space. Although well studied by both in situ and remote sensing instruments, robustly determining the solar origin of all the solar wind measured in situ by spacecraft is still an open problem. Excluding large ejecta, the solar wind has traditionally been classified according to the average speed of protons, which make up ∼95-99 per cent of the wind by ion number density. It is well known that the solar source of fast solar wind (v > ∼ 500 km s −1 ) is open field lines rooted in coronal holes on the surface of the Sun (Krieger, Timothy & Roelof 1973;Sheeley, Harvey & Feldman 1976;Cranmer 2009), and that it is composed of a slowly varying bulk speed superposed with shorter time-scale Alfvénic velocity fluctuations above this background level (Belcher & Davis 1971;Thieme, Schwenn & Marsch 1989;Matteini et al. 2015). The Alfvénic fluctuations are believed to be E-mail: david.stansby14@imperial.ac.uk generated close to the surface of the Sun where the wind is sub-Alfvénic, and then propagate outwards into the heliosphere on longlived open magnetic field lines (Belcher & Davis 1971;Cranmer & van Ballegooijen 2005). In addition to the global dynamics, the local kinetic properties of fast solar wind are also well known: the protons in fast solar wind have large temperature anisotropies and a low plasma β in the inner heliosphere (Marsch et al. 1982b;Marsch, Ao & Tu 2004;Matteini et al. 2007). In addition, alpha particles, which constitute 1-5 per cent of the solar wind by ion number density, exhibit large magnetic field aligned drifts relative to the protons in the fast solar wind (Marsch et al. 1982a;Steinberg et al. 1996).
In contrast to the fast wind, the plasma properties of the slow solar wind (v < ∼ 500 km s −1 ) are much more variable (Schwenn 2007), and although it must have a number of different solar sources (Abbo et al. 2016), it is not clear how these sources contribute to the different parts of the slow solar wind measured in situ. In general there are two possible generation mechanisms for the slow wind: it can flow continuously on magnetic field lines that maintain a connection from the base of the corona to the heliosphere, in a similar manner to the fast solar wind (Wang & Sheeley 1990;Cranmer, van Ballegooijen & Edgar 2007;Wang 2010), or can be released transiently from closed magnetic field lines undergoing interchange reconnection with the open magnetic field lines that connect to the heliosphere (Sheeley et al. 1997;Einaudi et al. 2001;Rouillard et al. 2010;Higginson et al. 2017).
During initial analysis of Helios data, Marsch et al. (1981) discovered a portion of slow solar wind measured at 0.3 au that, apart from its speed, had the same properties of fast solar wind: large protonalpha drift speeds, large proton core temperature anisotropies, and highly Alfvénic wave activity. In addition, Roberts et al. (1987) described an 80-d interval in the Helios data where the purest Alfvénic fluctuations were in slow solar wind. The strong Alfvénic wave activity during these periods implies that the wind was released on open magnetic field lines, allowing the Alfvén waves to freely propagate outwards from the corona to the point of measurement. More recently the presence of an 'Alfvénic slow wind' has been studied statistically using data taken at 1 au, and independent of solar activity the slow solar wind is composed of both non-Alfvénic and Alfvénic components (D'Amicis & Bruno 2015). These results hint that some slow solar wind has exactly the same source (and therefore properties) as fast solar wind, but just happens to be released at a slower speed.
In addition to protons and alphas, much less abundant heavy ions are measured in the solar wind, which can be used as a more direct proxy for solar source than the proton speed. As the solar wind travels away from the Sun it effectively becomes collisionless within a few solar radii (Hundhausen, Gilbert & Bame 1968). This means that ions are no longer able to gain or lose electrons through electron-ion collisions, and the fraction of different charge states become frozen. Ion charge state ratios therefore act as a tracer of the plasma properties at the freezing point. The most commonly used ratios are O 7 + /O 6 + and C 6 + /C 4 + , which are positively correlated with the electron temperature at the freezing point (Hundhausen et al. 1968;Bochsler 2007;Landi et al. 2012). Low charge state ratios are present in wind that originates in coronal holes, which have relatively low electron temperatures, and high charge state ratios are present in streamer belt plasma that has relatively high electron temperatures. This information can therefore be used to distinguish between coronal hole and non-coronal hole solar wind (Geiss, Gloeckler & Von Steiger 1995;Zhao, Zurbuchen & Fisk 2009). As an example, the Alfvénic slow wind identified by D' Amicis & Bruno (2015) had similar low charge state ratios to the fast solar wind, indicating it had similar solar origins (D'Amicis, Bruno & Matteini 2016), further reinforcing the need to go beyond classifying solar wind based solely on the average proton speed.
Although the Helios mission was equipped with an instrument for measuring heavy ions in the solar wind , it is believed that the data from this instrument have been lost. However, the case studies of D' Amicis &Bruno (2015) andD'Amicis et al. (2016) hint that when heavy ion measurements are not available the Alfvénicity of the solar wind fluctuations may be a more reliable proxy for solar wind origin than speed. The problem with using Alfvénicity as a categorization variable is that the solar wind becomes systematically less Alfvénic with distance (Roberts et al. 1987;Bruno et al. 2007;Iovieno et al. 2016) due to both small-scale turbulent evolution and large-scale velocity shears and interaction regions (Bruno et al. 2006). This means not all solar wind that started off Alfvénic near the Sun is still Alfvénic when it is measured at 1 au. In this paper, we mitigate this problem by using the unique Helios data, with measurements of solar wind plasma from 0.3 to 1 au, to link properties measured in situ, that are only observable at distances <1 au, to solar sources. The Helios data are described in Section 2, and statistical results are presented in Section 3. Based on these observations we construct three different categories of solar wind, and in Section 4 we place possible solar sources to each of these categories. In Section 5, we compare our categorization scheme with other schemes, and then conclude and present a set of predictions for the upcoming Solar Orbiter and Parker Solar Probe missions in Section 6.

DATA
The data used here were measured by the twin Helios spacecraft, which were operational during the late 1970s and early 1980s. Both spacecraft had an electrostatic analyser for measuring the ion distribution function at 40.5 s cadence (Schwenn, Rosenbauer & Miggenrieder 1975), and two different fluxgate magnetometers for measuring the magnetic field (Musmann et al. 1975;Scearce et al. 1975). Here we use a reanalysis of the ion distribution functions that fitted a bi-Maxwellian to the proton core population present in the experimentally measured ion distribution functions. This data set provides the proton core number density (n p ), velocity (v p ), temperatures parallel (T p ), and perpendicular (T p⊥ ) to the magnetic field, and corresponding magnetic field values (B) at a maximum cadence of 40.5 s. For more details on the fitting procedure and access to the data see Stansby et al. (2018a). Note that all the data presented in this paper are parameters of the proton core population, and not numerical moments of the overall ion distribution. The total temperature was calculated as T p = (2T p⊥ + T p )/3, the temperature anisotropy as T p⊥ /T p , the parallel plasma β as β p = 2μ 0 n p k B T p / |B| 2 , and the Alfvén speed as v A = |B| / √ n p m p μ 0 .
To avoid contamination of very large transients all of the intervals listed as coronal mass ejections by Liu, Richardson & Belcher (2005) were removed from the data set before further analysis. The state of the Sun undergoes an 11-yr solar cycle, oscillating between solar minimum and solar maximum. Because the highest quality Helios data were taken early in the mission, only data taken in the years 1974-1978 inclusive (during the solar minimum between cycles 20 and 21) were used.

Alfvénicity
In order to classify the solar wind as Alfvénic or non-Alfvénic, the cross-helicity was calculated in every 20-min interval where at least 10 velocity and magnetic field data points were available. The cross-helicity is defined as where v = v p − v p0 are the proton velocity fluctuations in the wave frame, v p0 is the local Alfvén wave phase velocity, b = v A (B/ |B|) is the magnetic field in velocity units, and indicates a time average over all points in a 20-min interval (Bruno & Carbone 2013). v p0 was calculated using the method given by Sonnerup et al. (1987), which finds the local de-Hoffman Teller frame of reference in which |v × b| 2 is minimized; by construction, this is the value of v p0 for which the absolute value of σ c is maximized. Although a plasma with Alfvén waves propagating in opposite directions can have low values of σ c , in this paper 'Alfvénic' is specifically reserved to denote a plasma where Alfvén waves propagate predominantly in only one direction. The magnitude of σ c indicates whether the fluctuations in the plasma are predominantly unidirectional Alfvén waves (|σ c | ≈ 1) or not (|σ c | < 1). For Alfvénic periods, the sign of σ c determines the direction of wave propagation with respect to the local magnetic field. Because Alfvén waves in the solar wind almost always travel away from the Sun (Gosling et al. 2009), the sign of σ c is a reliable proxy for the magnetic polarity of the solar wind.

Entropy
Heavy ion charge state data measured at 1 au is commonly used to diagnose the solar origin of solar wind. Unfortunately there is no heavy ion data available from the Helios mission, so instead proton specific entropy was used as a proxy. The specific entropy argument is easily calculated from the proton distribution parameters and given by where α is the polytropic index of the fluid. Here α was taken to be 1.5, the value used by Pagel et al. (2004) and Stakhiv et al. (2016) who studied correlations between entropy and composition, and whose results are used in Section 4.2 to make an indirect link between proton temperature anisotropy and heavy ion charge states using S p as an intermediate variable.

Statistics
Fig . 1 shows the distribution of the solar wind in the β p -T p⊥ /T p plane measured by Helios at heliocentric distances of 0.3-0.4 and 0.9-1.0 au, and split into slow and fast wind using a simple cut in speed. The distribution at 0.3 au has previously been presented for an individual high speed stream by Matteini et al. (2007). The top right-hand panel of Fig. 1 shows that the ∼100-h continuous high As the wind propagates outwards the protons become more isotropic and the plasma β increases, primarily due to adiabatic evolution (Chew, Goldberger & Low 1956;Matteini et al. 2011), although by 1 au the fast wind protons have not yet reached the equilibrium state of T p⊥ /T p = 1. In contrast, at 1 au the slow solar wind is distributed around T p⊥ /T p = 1, where it is thought to be maintained during transit by a combination of collisions and kinetic instabilities that are active when β p ≥ 1 (e.g. Kasper, Lazarus & Gary 2002;Hellinger et al. 2006;Bale et al. 2009;Yoon 2016). At 0.3 au the majority of the slow solar wind is also spread around T p⊥ /T p = 1, but there is a significant fraction with T p⊥ /T p > 2 and β p 1 (Fig. 1, top left). This is the same region in the parameter space that fast solar wind occupies at 0.3 au, which implies that there is a portion of slow solar wind that has the same kinetic properties as fast solar wind. Instead of partitioning the data by speed, Fig. 2 shows the joint distribution of T p⊥ and T p for all measurements between 0.3 and 0.4 au. In this parameter space two populations are clearly distinct: one centred around T p⊥ /T p = 1 and other centred around T p⊥ /T p = 3. All fast solar wind occupies the anisotropic population, but the slow wind is split between the two populations (as shown in Fig. 1).
In addition to being anisotropic, the fast solar wind is filled by antisunward propagating Alfvén waves. This is a statement about the global dynamics of the plasma, in contrast to the local kinetic properties given by the parallel and perpendicular temperatures, and is a consequence of the wind being released on long-lived open field lines. To investigate whether the anisotropic wind as a whole has the same Alfvénic property as the fast solar wind, Fig. 3 shows the joint Figure 3. 2D histograms of temperature anisotropy against absolute cross-helicity (main panels) with adjoining 1D histograms of temperature anisotropy (right-hand panels) and cross-helicity (top panels). 2D bin counts are normalized to the maximum bin count. 1D histograms are linearly scaled. The horizontal and vertical dashed lines show the partitioning of the parameter space into three distinct regions; see text in Section 3.1 for more details. Percentages indicate the fraction of data points in each of the three regions. probability distribution of temperature anisotropy and cross-helicity at both 0.3 and 1 au. The distribution of temperature anisotropy is clearly bimodal at 0.3 au with a minimum at T p⊥ /T p = 1.7, a feature that is not observable at 1 au due to the 0.3 au anisotropic component becoming more isotropic with radial distance (as a result of adiabatic evolution). The distribution stops being clearly bimodal at radial distances greater than around 0.8 au (not shown). In the rest of this paper data taken at 0.3-0.4 au are presented, but the qualitative properties discussed are present at all radial distances from 0.3 to 0.8 au. It is clear from Fig. 3 that the fraction of solar wind that is Alfvénic is much higher at 0.3 au (∼80 per cent) as compared with 1 au (∼50 per cent), which agrees well with previous studies (Roberts et al. 1987;Bruno et al. 2007). At 0.3 au the anisotropic population is almost always Alfvénic (i.e. |σ c | > 0.8), but this is not the case for the isotropic wind. We therefore propose splitting the solar wind at 0.3 au into three populations based on the observable boundaries in this parameter space: (i) an anisotropic, Alfvénic population; (ii) an isotropic, Alfvénic population; (iii) an isotropic, non-Alfvénic population.
The split in anisotropy was chosen to be the saddle in between the two populations at T p⊥ /T p = 1.7, and the split in cross-helicity was chosen at the edge of the Alfvénic population at |σ c | = 0.8. These boundaries are shown in Fig. 3. At 0.3-0.4 au, ∼80 per cent of the wind was Alfvénic, split equally between isotropic and anisotropic, and the remaining 20 per cent was non-Alfvénic.
With this classification in mind, the top panel of Fig. 4 shows the radial velocity distribution of the solar wind in each category. Both isotropic populations consist primarily of solar wind with speeds less than 500 km s −1 , whereas the anisotropic population spans a wide range of speeds from 300 to 700 km s −1 . In fact, at 0.3-0.4 au Helios measured slightly more anisotropic solar wind below 500 km s −1 than above. This reinforces the idea that the concept of fast and slow solar winds breaks down at intermediate speeds where wind can have properties similar to either the very slow or very fast wind, and again suggests that some slow wind may have the same origin as fast wind. Another known property of the fast solar wind is that the radial mass flux does not depend on speed (Feldman et al. 1978;Wang 2010). To investigate whether this is true for the anisotropic wind as a whole, the main panel of Fig. 4 shows radial flux as a function of radial velocity for each of the three categories. The anisotropic wind has a constant flux that does not depend on speed; when evaluated at 1 au this average flux is ∼1.8 × 10 8 cm −2 s −1 , which agrees well with independent measurements made at 1 au and beyond (Phillips et al. 1995;Goldstein et al. 1996;Wang 2010). The isotropic Alfvénic wind has a slightly higher flux, whereas the non-Alfvénic wind has widely ranging fluxes varying up to 1-4 times the base value of the anisotropic wind, suggesting a distinct physical origin.
To investigate the link between the three categories and their source regions on the Sun, we also looked at the dependence of temperature anisotropy on proton specific entropy (used later as a proxy for composition). Fig. 5 shows the joint distribution of temperature anisotropy, and solar wind speed or proton specific entropy. At 0.3 au, low-speed wind (200-300 km s −1 ) is all isotropic, and at high speeds (450-700 km s −1 ) all of the wind is anisotropic, however, at intermediate speeds (300-450 km s −1 ) the wind is spread between the two different states of T ⊥ /T . In contrast, the variation between temperature anisotropy and specific entropy is slightly smoother; isotropic wind corresponds exclusively to low entropy and anisotropic wind exclusively to high entropy, with a continuous variation in between. This result is used only as a correlation, and we do not claim that there is any causal relationship between entropy and temperature anisotropy. In Section 4.2, we discuss how this correlation can be used as an intermediate step to infer the compositional properties of our three different categories.

Spatial distribution of the three solar wind populations
We have shown that at 0.3 au it is possible to distinguish between three types of solar wind based on the statistics of proton temperature anisotropy and Alfvénicity. To understand the spatial distribution of each population within the solar wind, Fig. 6 shows the time series measurements made by Helios 1 inside 0.5 au during its first perihelion pass.
The bimodal nature of proton temperature anisotropy is clear, even within the unaveraged and noisy 40.5-s cadence measurements, and the wind remained in one anisotropy state for days at a time. In contrast, within the isotropic category the Alfvénic and non-Alfvénic subcategories are well mixed and interspersed within each other. During some isotropic periods (e.g. around 1975-03-09) the non-Alfvénic wind is subdominant and appears to be embedded in the Alfvénic wind, whereas at other times (e.g. around 1975-03-30) there appears to be an approximately even mix of Alfvénic and non-Alfvénic wind.
The transition from isotropic to anisotropic wind was sharp, and always occurred at the leading edge of high-speed streams. It is known that composition and entropy undergo sharp changes at the leading edge of high-speed streams (Wimmer-Schweingruber, von Steiger & Paerli 1997;Lazarus et al. 2003;Crooker & McPherron 2012), but Fig. 6 demonstrates that a coincident temperature anisotropy boundary is also present. The sharp increases temperature anisotropy were driven by increases in T ⊥ , whilst T stayed roughly constant across the boundaries (not shown). The transition from anisotropic back to isotropic was also sharp, caused by sharp decreases in T ⊥ , and occurred inside the rarefaction edge of high-speed streams. The sudden drop in T ⊥ , which caused a coincident drop in total temperature and therefore a drop in specific entropy, was not correlated with changes in the cross-helicity. The only other observable change in the plasma and magnetic field data are magnetic field fluctuations that look qualitatively different on either side of the boundary (not shown). The time series gives a clear visual demonstration that it is impossible to cut the velocity in a single place to separate different types of wind, but clear bimodality makes a cut in temperature anisotropy easy. We reiterate that performing this separation is only possible at heliocentric distances <0.8 au, as at large distances the anisotropic wind becomes more isotropic, such that the two populations are no longer separable.

L I N K I N G IN SITU M E A S U R E M E N T S TO S O L A R S O U R C E S
We now use the observations made in Section 3 to link our three categories of solar wind to their solar sources. A summary of the conclusions drawn in this section is given in Table 1.

Known properties of coronal hole wind
It has long been known that wind originating on open field lines rooted inside large coronal holes forms the fast solar wind (Krieger et al. 1973;Sheeley et al. 1976;Cranmer 2009). Remote sensing measurements also show pronounced temperature anisotropies present above coronal holes, whilst the solar wind is still near to the Sun (Kohl et al. 1997;Cranmer, Panasyuk & Kohl 2008). Because the anisotropic category is the only one with high speeds (Fig. 4), we infer that wind produced in the core of coronal holes belongs to our anisotropic category. The spatial distribution of anisotropic wind, with slower speeds always occurring in the rarefaction edges of high-speed streams (Fig. 6), shows that rarefaction during transit is responsible for the relatively low speed of some anisotropic wind (Pizzo 1991). The reason slower speeds are not observable at the leading edge of high-speed streams is because by 0.3 au they have already been accelerated by the faster wind to form a corotating interaction region (Burlaga 1974;Pizzo 1991;McGregor et al. 2011;Richardson 2018).
Note that we have chosen to distinguish between the edges and the core of coronal holes; at the edge of coronal holes the magnetic field lines typically undergo large separations as a function of height in the corona, which has the effect of reducing both the wind speed (Levine, Altschuler & Harvey 1977;Wang & Sheeley 1991;Cranmer et al. 2007;Pinto, Brun & Rouillard 2016) and charge state ratios (Wang, Ko & Grappin 2009). In the next two sections further evidence is used to predict which one of our three categories coronal hole edge wind is part of. . The line at T p⊥ /T p = 1.7 shows the boundary between anisotropic (points above) and isotropic (points below) wind. Bottom panel shows absolute values of 20-min cross-helicity, coloured by sign. The line at |σ c | = 0.8 shows the threshold separating Alfvénic (points above) and non-Alfvénic (points below) wind.

Correlation of anisotropy, entropy, and composition
At distances beyond 1 au the proton specific entropy is anticorrelated with the O 7 + /O 6 + charge state ratio (Pagel et al. 2004). In addition, observations at 1 au show that specific entropy is anticorrelated with the C 6 + /C 4 + charge state ratio (Stakhiv et al. 2016). We have shown in Fig. 5 that entropy has a monotonic dependence on proton temperature anisotropy (but note again that this is not necessarily a causal relationship). Linking this newly observed relationship to the inferred relationship between entropy and charge state ratios suggests that anisotropic wind has low O 7 + /O 6 + and C 6 + /C 4 + charge state ratios, and isotropic wind has high charge state ratios. In addition to being related statistically, the sharp boundaries between anisotropic wind and isotropic wind mimic the locations of sharp composition boundaries found in other studies (see Section 3.2 for a discussion). This backs up the statistical inference derived between proton temperature anisotropy and heavy ion charge states. Using specific entropy as a bridge between anisotropy and composition therefore corroborates our previous conclusion that anisotropic wind originates in the core of large coronal holes, which are known to emit wind with low charge state ratios (Geiss et al. 1995;Wang et al. 2009). This allows us to infer that both categories of isotropic wind do not originate in the core of coronal holes, but may originate at coronal hole edges or outside coronal holes. This again agrees with remote sensing measurements that show reduced temperature anisotropies near the edges of coronal holes when compared to the core of coronal holes (Susino et al. 2008). We therefore suggest that the 'isotropic' wind forms what is commonly thought of as the 'slow solar wind'. There are a number of theories as to the origin of the slow solar wind (Abbo et al. 2016); in the next section we assign possible theories to either the Alfvénic or non-Alfvénic categories of the isotropic wind.

Alfvénicity and mass flux variability
Solar wind with a high Alfvénicity implies a steady state release of plasma on open field lines. This hypothesis is supported by the relatively constant mass flux in the Alfvénic isotropic wind (see Fig. 4). This means any Alfvénic wind must have been released on field lines that remained open for at least the 20-min resolution of the cross-helicity calculated from in situ data. Areas of long lasting open field on the Sun can be split into the core of coronal holes (already categorized), edges of coronal holes, and active regions. Remote sensing measurements have shown that active region outflows have high coronal electron temperatures (Neugebauer et al. 2002;Brooks & Warren 2012), contain open field lines allowing plasma to escape into the heliosphere (Slemzin et al. 2013), and can supply mass fluxes similar to those measured in situ (Brooks, Ugarte-Urra & Warren 2015). We therefore conclude active region wind is most consistent with the isotropic Alfvénic category, along with the wind from the edges of coronal holes that also contain long lasting open magnetic fields and has similar properties.
Finally, we predict that the non-Alfvénic wind is consistent with the final known type of slow solar wind, typically called 'number density structures' or 'blobs', which have been detected both remotely (e.g. Sheeley et al. 1997;Viall & Vourlidas 2015;DeForest et al. 2018) and in situ (e.g. Kepko & Spence 2003;Sanchez-Diaz et al. 2017;Stansby & Horbury 2018). These are non-steady state transient structures with a high density but similar speed as the surrounding slow wind, and therefore have enhanced mass fluxes relative to the background wind. This property is exactly what we have measured for the non-Alfvénic wind (Fig. 4), backing up our final categorization.
A summary of our mapping of possible solar sources to in situ solar wind categories is given in Table 1.

C O M PA R I S O N W I T H OT H E R C AT E G O R I Z AT I O N S C H E M E S
Recently several authors have also attempted to categorize sources of the solar wind using in situ observations, choosing the categories of coronal mass ejection wind, coronal hole wind, and interstream wind (for a summary, see Neugebauer, Reisenfeld & Richardson 2016). In this paper, we have deliberately removed coronal mass ejection wind from our data set, and have used proton temperature anisotropy as the only variable separating coronal hole wind (anisotropic) and interstream wind (isotropic). Zhao et al. (2009) used only the O 7 + /O 6 + ratio and solar wind speed measured at 1 au. This method is limited by the slow cadence (1 h) of charge state ratio measurements available, but has the advantage that the O 7 + /O 6 + ratio is known to be directly related to plasma properties near the Sun. Because the distribution of heavy ion charge states is only clearly bimodal at solar minimum , it is not clear if this method works well during solar maximum conditions. Because our method also assumes a bimodal distribution of charge state ratios, it is not clear if it is still applicable during solar maximum conditions either. Xu & Borovsky (2015) picked in situ measurements, manually categorized specific intervals of the measurements, and then tried to find boundaries in a multidimensional parameter space that reliably split the data into the assumed categories. These boundaries could then be applied to other intervals where the categorization is unknown. This method is practical and pragmatic for rapidly categorizing solar wind sources, but the boundaries between categories are somewhat arbitrary and do not necessarily directly relate to the different physics of solar wind formation at each solar wind source. The advantage of the Xu & Borovsky (2015) method is that it only uses single-point measurements of solar wind protons and magnetic fields, so the cadence at which it can be applied is limited only by that of the in situ measurements. In contrast our method is limited to a 20-min cadence, which is in practice limited by the number of 40.5-s cadence of proton measurements needed to reliably calculate σ c .
Other authors have backmapped solar wind measured at 1 au to try and determine the exact location on the Sun from which it originated (e.g. Neugebauer et al. 1998;Fu et al. 2015;Fazakerley, Harra & van Driel-Gesztelyi 2016;Peleikis et al. 2017;Zhao et al. 2017). This method assumes that the solar wind travels along magnetic field lines between the Sun's surface and magnetic field source surface at 2.5r s (solar radii), which are computed using a potential field source surface model, and then travels radially and at a constant speed to the in situ observer. This method has the advantage of drawing a direct link by trying to predict the exact solar wind source location of in situ measurements. Although it is successful in identifying sources on very large time-scales of ∼days, it is currently not possible to probe smaller scales, and does not take into account dynamical processing that occurs during transit between the Sun and the in situ observer.

C O N C L U S I O N S A N D P R E D I C T I O N S F O R F U T U R E M I S S I O N S
We have presented an attempt to map in situ measurements of solar wind to their sources, using properties of the solar wind observable at 0.3 au that are unobservable at 1 au due to dynamical interactions during transit. We find that the solar wind can be split into three categories (summarized in Table 1), based on in situ measurements of proton temperature anisotropy and Alfvénicity (Section 3), and sort possible solar origins of the solar wind into these three categories (Section 4). Although many other methods have been developed to attempt the same goal of solar source categorization (Section 5), the lack of in situ composition and remote sensing data available during the Helios era (1974)(1975)(1976)(1977)(1978)(1979)(1980)(1981)(1982)(1983)(1984) restricted our ability to use these more modern techniques. However, in the near future we will have access to simultaneous in situ measurements of protons in the inner heliosphere, in situ measurements of solar wind composition, and a wide range of remote sensing data. We finish by describing how our new categorization scheme can be applied to data from upcoming missions to the inner heliosphere.
Parker Solar Probe (PSP; Fox et al. 2016) will make in situ measurements of the solar wind at heliocentric distances inside 0.3 au, and the first comprehensive in situ solar wind measurements inside 1 au since Helios. Proton and magnetic field measurements made by PSP will allow us to perform the categorization scheme outlined in this paper. Advances in modelling and remote sensing since the Helios era mean that it will then be possible to backmap the solar wind measured by PSP to a predicted source location on the Sun. If our categorization is correct, the three categories of in situ solar wind will backmap to their respective inferred solar sources. Solar Orbiter (SO; Müller et al. 2013) will provide the first solar wind composition measurements between 0.3 and 1 au. This will allow us to directly test the correlation between temperature anisotropy and charge state ratios, without having to bridge the gap by using proton specific entropy as an intermediate variable. If our categorization is correct, isotropic wind will clearly correspond to high charge state ratios, and anisotropic wind will clearly correspond to low charge state ratios. In addition SO will carry on board remote sensing instruments that are designed to target the predicted solar sources of in situ measurements, making backmapping wind to its source even more accurate than using remote sensing instruments at 1 au.