Open Access
Issue
A&A
Volume 613, May 2018
Article Number A62
Number of page(s) 7
Section The Sun
DOI https://doi.org/10.1051/0004-6361/201732567
Published online 01 June 2018

© ESO 2018

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://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

The solar wind is traditionally viewed as bimodal, consisting of fast (> ~500 km s−1) and slow (< ~500 km s−1) components. It is well established that fast solar wind originates from coronal holes on the Sun (Sheeley et al. 1976), however the origins of the slow solar wind are still unclear (Abbo et al. 2016). There are two types of theory that seek to explain the origin of the slow solar wind. The first type states that slow wind originates in open magnetic field regions, where the field at the base of the corona is permanently connected to the heliosphere. These regions can either be the edge of coronal holes (Levine et al. 1977; Wang & Sheeley 1991) or the edge of active regions (Harra et al. 2008, 2017; He et al. 2010). To produce slow wind speeds they must be areas where the expansion of the magnetic field between the base of the corona and the heliosphere is largest (Leer & Holzer 1980; Cranmer et al. 2007; Wang 2010) or where field line deviations from the radial direction at the coronal base are largest (Pinto et al. 2016). This theory predicts that the release of slow solar wind plasma into the heliosphere is a continuous steady state process.

The second set of theories are based on slow wind that starts in closed field regions, where the coronal magnetic field is not connected to the heliosphere. In these theories, interchange magnetic reconnection between adjacent open and closed field lines allows plasma that starts on closed field lines to access the open field lines and escape into the heliosphere (Crooker et al. 2002, 2004). This can either happen throughout the closed field regions (Fisk & Schwadron 2001; Fisk 2003; Fisk & Zhao 2008) or only at the open-closed boundary (Antiochos et al. 2011, 2012; Higginson et al. 2017). These theories predict that the release of plasma is a transient process, occurring only when interchange reconnection takes place.

One transient component of the slow solar wind that has been observed both remotely and in-situ are discrete structures with enhanced number density. They were first discovered in remote sensing measurements by Sheeley et al. (1997) using white light images, formed when sunlight Thompson scatters off electrons in the solar wind. When viewed in white light the structures are associated with the top of streamers, which are loops of closed magnetic field lines that separate open field regions of opposite polarity (Rouillard et al. 2010a, 2011; Viall & Vourlidas 2015). The two dimensional extension of the polarity inversion into the heliosphere forms the heliospheric current sheet (HCS), which starts above the tips of streamers, and is surrounded by dense material called the heliospheric plasma sheet (Winterhalter et al. 1994). Further away from the Sun dense structures form part of the material that surrounds the HCS (Wang et al. 1998). Once in the heliosphere their trajectories match what would be expected if they are periodically ejected from a fixed point co-rotating with the surface of the Sun (Sheeley & Rouillard 2010) and they are frozen in to and travel at the same speed as the surrounding solar wind (Viall et al. 2010). All of these features agree qualitatively with interchange reconnection taking place near the edge of streamers, releasing dense plasma to form the structures which subsequently surround the HCS as they travel outwards into the heliosphere (Sanchez-Diaz et al. 2017).

Remote observations of structures in the solar wind have a number of limitations. To be detectable the structures must have a large enough density relative to the background solar wind to scatter enough sunlight, have a large enough size to be resolvable by a finite resolution imaging instrument, and images must be exposed on short enough timescales so any measured intensity enhancements do not “wash out” due to smearing across multiple pixels. Some success has been had tracking individual structures up to 1 AU, however this often relies on the structures being compressed to increase their number density above a theoretically expected 1∕r2 decrease (Sheeley & Rouillard 2010; DeForest et al. 2011; Howard et al. 2012).

In-situ measurements do not suffer from these restrictions but instead are limited by the location of the observing spacecraft. Studies using data taken at 1 AU have shown individual cases of proton number density enhancements at the HCS that were accompanied by similar enhancements in Helium, Oxygen, and Carbon number densities (Kepko et al. 2016). The C6+ /C5+ ratio was used to infer the temperature of the corona where each packet of solar wind originated. The derived coronal temperatures were very variable on timescales of the same order as the structure sizes, indicating that the release of slow solar wind was not a steady state process, but there was no clear correlation between the location of enhanced proton number density and inferredcoronal temperatures. This meant the solar origin of the structures could not be unambiguously identified.

Previous studies, both remote sensing and in-situ, give an overview of structures generated near coronal streamers which form part of the material that surrounds the HCS in the heliosphere. In this paper we present the first in-situ statistical study of number density structures located in the slow solar wind, but away from in-situ crossings of the HCS. When observing the solar wind in-situ at 1 AU it is difficult to determine whether observations are truly a reflection of source properties, or are as a result of dynamical interactions as the solar wind travels radially outwards. For this reason we use data from the Helios mission taken between 0.3 AU and 0.5 AU, described in Sect. 2. In Sect. 3.1 we present observations of a single representative structure, and in Sect. 3.2 a survey of 140 identified structures summarises the plasma properties consistent across all structures. Finally, in Sect. 4 we discuss the implications of these observations for the sources and formation of the slow solar wind.

2 Data

2.1 Ion distribution function processing

This study uses data from the Helios mission (Porsche 1977). Helios 1 and Helios 2 travel in elliptical heliocentric orbits between 0.3 AU and 1 AU, and were operational between 1974–1985 and 1976–1980 respectively. Each spacecraft has an electrostatic analyser for measuring ion distribution functions, which measured a full 3D distribution every 40.5 s (Schwenn et al. 1975). The 3D distribution was then transmitted to the Earth; due to telemetry constraints the data is often not available at the full 40.5 s cadence. We have fitted the original ion distribution functions with a single bi-Maxwellian to extract the number density (np ), bulk velocity (v p ), and magnetic field parallel (Tp) and perpendicular (Tp) temperatures of the proton core population. The fitting was done to each distribution using a least squares method that minimised the cost function C= [ f(np,vp,Tp,Tp)fdata ]2\begin{equation*} C = \sum \left [f\left (n_{p}, \mathbf{v}_{p}, T_{p\parallel}, T_{p\perp} \right ) - f_{\mathrm{data}} \right ]^{2} \end{equation*}(1)

where f(np,vp,Tp,Tp)$f\left (n_{p}, \mathbf{v}_{p}, T_{p\parallel}, T_{p\perp} \right )$ is the fitted bi-Maxwellian distribution function, fdata is the measured distribution function, and the sum is taken over all points in fdata . Empirically the fitting process did a good job of fitting to only the proton core population, and was not affected by either the faster proton beam or alpha particle population that are sometimes present in the solar wind. An example fit can be viewed in Fig. 1 of the guide to the dataset (Stansby 2017). Along with the proton parameters, simultaneous magnetic field values were calculated as the average magnetic field during the time it took for each ion distribution function to be measured. Where available 4 samples/s data from the E2 fluxgate magnetometer were used (Musmann et al. 1975). When this was not the case 6 s cadence data from the E3 fluxgate magnetometer were used (Scearce et al. 1975).

The fitting process was applied to all the original ion distribution functions available from both Helios 1 and Helios 2. The code used to generate the dataset is openly available, making the fitted parameters easily reproducible. The source code and dataset, along with a more detailed description of the fitting procedure is available online (Stansby 2017).

From the fitted parameters the total proton temperature was calculated as Tp=(2Tp+Tp)/3$ T_{p} = \left (2T_{p\perp} &#x002B; T_{p\parallel} \right ) / 3 $, from which the proton thermal pressure is given by pp = npkBTp. The proton plasma beta is defined as βp = pppmag and the magnetic pressure given by pmag=| B |2/2μ0$ p_{\mathrm{mag}} = \left | \mathbf{B} \right |^{2} / 2\mu_{0} $.

As well as protons, the other major constituents of the solar wind are alpha particles and electrons. Because corresponding electron and alpha number densities, velocities and tempereatures are not available it was not possible to calculate their contributions to the thermal pressure on a point by point basis. In slow solar wind alpha particle temperatures are between 0.2–1 MK (Marsch et al. 1982), the alpha to proton number density ratio between 0.01–0.05 (Kasper et al. 2007), and electron temperatures lie in the range 0.2–0.4 MK (Štverák et al. 2015). Assuming plasma neutrality the proton, alpha, and electron number densities are related by np + 2nα = ne. This relation along with the ranges of possible alpha and electron parameters allowed lower and upper bounds to be placed on the total pressure, defined as ptot = (pp + pα + pe) + pmag. In the results section values of ptot are presented as bounds or error bars between the minimum and maximum possible values.

2.2 Number density structure selection

From the entire Helios plasma data set a search was performed to find all the occurrences of number density structures. The search was limited to heliocentric distances ranging from 0.3 AU to 0.5 AU to minimise the time for stream-stream interactions to alter the local plasma properties. Structures were identified by eye in time series data as a clear sharp increase and subsequent decrease in proton number density accompanied by a similar signature in βp . An example is shown in the bottom set of panels in Fig. 1. To distinguish structures from background noise they were required to last for more than 3 data points (~2 min).

Simultaneous increases in np and βp are also properties of the heliospheric plasma sheet, which surrounds the HCS (Winterhalter et al. 1994). In order to clearly distinguish between the heliospheric plasma sheet and transient structures, any occurrences of the plasma sheet were avoided by ignoring enhanced number densities that surround the HCS. The HCS was identified by a change in sign of the radial magnetic field component (Br ) and coincident change in sign of the magnetic field-velocity correlation. The magnetic field-velocity correlation was used to distinguish between global polarity reversals and magnetic field foldbacks that only locally change the sign of Br . This is similar to the common method of using solar wind electrons as a diagnostic (e.g. Owens et al. 2017), but replaces the direction of electron strahl propagation with the direction of Alfvén wave propagation. This method assumes that Alfvén waves are always travelling along magnetic field lines away from the end connected to the Sun; this a safe assumption to make as very few examples of waves travelling towards the Sun have been reported (Gosling et al. 2009; He et al. 2015). Reconnection exhausts in the solar wind also exhibit the same increases in np and βp (Gosling et al. 2006). To avoid selecting such exhausts any events which contained a smooth reversal in the magnetic field polarity were also discarded.

In total 140 structures were identified. The observing spacecraft and times of each structure are available online (Stansby 2018b).

3 Results

3.1 An example density structure

Figure 1 shows an example density structure, observed at 0.32 AU by Helios 1 on 28 March 1976. The top three panels show a 4 day overview of 10 minute averaged data, with the structure situated between the vertical black lines. The structure was located in slow solar wind (vr ≈ 350 km s−1) with a clearly defined negative magnetic polarity, around 18 h before the in-situ HCS crossing.

The bottom set of panels shows a shorter time interval, with the structure beginning shortly after 22:15 and lasting around 20 min. It is visible as a proton number density enhancement, from a background level of 75 cm−3 to a peak of 125 cm−3. The background number density had a similar value before and after the structure. The relative density increase along with similar before and after density values are comparable to previous in-situ observations of structures at 1 AU (Howard et al. 2012; Kepko et al. 2016). Outside the structure βp ≪ 1 as is typical at 0.3 AU, whereas during the structure βp rose to 0.5.

Both of the proton temperatures were higher during the structure; the parallel temperature increased by a factor of 2 from 0.07 MK to 0.15 MK, and the perpendicular temperature underwent a smaller increase from 0.05 MK to 0.08 MK. This equated to a doubling in total temperature from ~0.06 MK to ~0.12 MK. The large increase in thermal pressure (caused by both increased number density and temperature) was offset by a similar decrease in magnetic pressure. Within uncertainty in electron and alpha particle pressures this resulted in no variation in the total pressure, meaning the structure was pressure balanced with its surroundings. The proton velocity was slightly higher inside the structure, rising from 330 km s−1 to 360 km s−1 (not shown). Because the structure was pressure balanced and travelled at a similar speed to the surrounding plasma, it is not expected to have undergone a significant interaction with the surrounding solar wind as it travelled away from the Sun.

3.2 Statistical study

For all structures identified using the method described in Sect. 2.2 average densities, velocities, temperatures, magnetic field strengths, and pressures inside each structure were computed, along with the same quantities averaged outside each structure. For quantities outside the structures an average of 3 data points (≈2 min) before and 3 data points after the structure were taken.

Figure 2 shows scatter plots of quantities inside the structures against average quantities outside. Each point in the plots represents a single structure. The top left panel shows that number densities were always higher in the structures, as originally selected for. The other five panels show quantities that were not visible during manual selection of the structures. Temperatures were generally larger inside the structures, with an upper limit of 2 times hotter than the surrounding plasma. This increase in temperature was not consistently partitioned between either the parallel or perpendicular temperatures, but the temperature anisotropy was often significantly different inside the structures. The magnetic field strength was lower inside the structures compared to outside. This indicates that thermal pressure increases caused by number density and temperature increases were accompanied by coincident magnetic pressure decreases. Thermal pressures were higher inside the structures, magnetic pressures lower, but the total pressure of the surrounding solar wind was similar to the total pressure inside the structures. This is consistent with the structures being pressure balanced with their surroundings. The magnetic strength decreases inside the structure were often accompanied by large changes in the individual field components. Variations in the components did not consistently happen in one direction, but the field magnitude and direction was consistently the same before and after each structure (not shown). Finally, the radial velocity of the structures was very similar to the surrounding solar wind.

thumbnail Fig. 1

An example number density structure. Top three panels: a 4 day overview of 10 min averaged radial proton bulk speed (vr, top axis), radial magnetic field (Br, middle axis), and proton number density (np, bottom axis). The HCS crossing is marked with a red dashed line between days 87 and 88. Bottom four panels: a shorter 1.5 h interval, with the structure present between the grey dashed lines. The top axis shows proton number density (np ), and the second axis proton plasma beta (βp). The third axis shows proton parallel temperature (orange, Tp) and proton perpendicular temperature (blue, Tp). The fourth axis shows the proton thermal pressure (red, pp), magnetic pressure (blue, pmag), and total pressure (grey, ptot). The range of values in the total pressure is due to uncertainty in the number density and temperature of electrons and alpha particles.

thumbnail Fig. 2

Scatter plots of quantities inside the structures (y-axis) against quantities outside the structures (x-axis). One-to-one line are drawn in black on each plot. Points lying above these lines indicate a quantity larger inside a structure relative to outside and vice versa. Top left: proton core number density (np ) and top right total proton temperature (Tp). Middle left: magnetic field magnitude (| B |$\left | \mathbf{B} \right |$) and middle right magnetic (blue) proton thermal (red) and total (black) pressures. The total pressure crosses are error-bars in both directions due to uncertainty in the electron and alpha thermal pressure contribution to the total pressure. Bottom left: radial proton velocity (vr ) and bottom right proton temperature anisotropy (TpTp).

3.3 Viewing the structures in white light images

Structures can be observed either in-situ or using remote sensing. In this section we link these two measurement techniques by using in-situ observations to predict which individual structures would be visible in white light images. Thetwo structure properties that directly determine visibility in white light images are the fractional increase in electron number density and line of sight structure size. The fractional increase in number density was calculated as δnn0, where n0 is the average proton number density outside the structure and δn= n n0$\delta n = \left \langle n \right \rangle - n_{0}$ is the average number density increase inside the structure. The available in-situ measurements are of proton number density, which does not necessarily equal the electron number density. The other major positive ion species in the solar wind is a population of alpha particles, therefore the electron density ratio δnn0 is only equal to the proton density ratio under the assumption that the proton to alpha particle ratio does not change across the structure. This may not be the case and in the slow solar wind values of the alpha particle abundance range between 1 and 5% (Kasper et al. 2007). This range of possible abundances translates to a maximum uncertainty of ± 8% in the electron number density ratio.

When structures are viewed remotely the line of sight structure size is perpendicular to the radial direction of propagation. Only the flow parallel size is available from single spacecraft in-situ measurements, so this was used as a proxy for the line of sight size. This assumes that the structures have an aspect ratio of 1:1. The size of each structure was calculated as l=t| vp |$l = t \left | \left \langle \mathbf{v}_{p} \right \rangle \right |$ where l is the total advected size, t is the structure duration, and vp $\left \langle \mathbf{v}_{p} \right \rangle$ is the average proton velocity inside the structure.

The visibility of the structures also depends on their location relative to the observing instrument. The strongest white light signal is achieved when a structure is located on the Thompson sphere, which has its centre half way along the observer-Sun line and passes through the observer and the Sun (Howard & Tappin 2009). In Appendix A it is shown that for a structure located on the Thompson sphere the fractional increase in intensity when viewing the structure is given by δII00.64lr0δnn0\begin{equation*} \frac{\delta I}{I_{0}} \approx 0.64 \frac{l}{r_{0}} \frac{\delta n}{n_{0}}\end{equation*}(2)

where δI is the increase in white light intensity received when looking at a structure, I0 is the background intensity when the structure is not present, and r0 is the distance from the centre of the Sun to the structure. δII0 is used directly when detecting number density structures in white light images.

Figure 3 shows a scatter plot of fractional number density enhancement against structure size, with constant contours of δII0 from Eq. (2) in the background. The structures occur on scales ranging from 50 Mm to 2000 Mm, with no one characteristic scale size. The lower bound of 50 Mm is an instrumental limit caused by the 40.5 s cadence of plasma measurements available from Helios, along with the detection requirement that the structure duration spans at least three consecutive data points. The density contrast of the structures ranges from 10% above background levels to 200% above background levels with no single characteristic increase. Again the lower bound of 10% is likely to be an artificial limit set by the minimum increase needed to manually distinguish a structure from noise in the proton number density measurements.

Equation 2 assumes that a structure stays in a single pixel location during the time the image is exposed. In reality this is not the case, and δII0 is reduced by a factor of the number of pixels the structure travels through during a single exposure. This is known as “velocity smear”. For typical slow solar wind speeds of 300 km s−1 the instrumental velocity smear on Stereo HI-A reduces visibility by a factor of 10 (Howard et al. 2008; Viall et al. 2010). Typical values for previously observed structures were in the range δII0 = 2 × 10−4–10−3 (Viall et al. 2010), where the lower bound is due to the instrumental noise limit. Taking into account the velocity smear, values calculated using Eq. (2) must satisfy δII0 > 2 × 10−3 to be visible. 70% of structures in Fig. 3 do not meet this criteria, and thus would not be observable in Stereo HI-A images.

thumbnail Fig. 3

Scatter plot of electron density contrast (y-axis) against advected structure size (x-axis) . Uncertainties in electron density contrast are due to not knowing the variation of the alpha particle abundance within the structures. Background shading shows contours of Eq. (2) with r0 = 0.3 AU, which gives the fractional change in intensity a white-light instrument at 1 AU would see when viewing the blobs. The fractional intensity change values assume zero velocity smear across the observing instrument.

4 Discussion

In the previous section 140 discrete areas of enhanced number density coinciding with enhanced proton beta were identified (for an example see Fig. 1). It was then shown that compared to the surrounding plasma these regions were hotter, had lower magnetic pressure and higher thermal pressure, but remained pressure balanced, and travelled at similar speeds to the surrounding solar wind (Fig. 2). The structures here had no characteristic size or density increase, and 70% of them would not be visible in remote images taken at 1 AU (Fig. 3).

These observed features agree with previous individual in-situ measurements of pressure balanced structures in the inner heliosphere and at 1 AU (e.g. Burlaga & Ogilvie 1970; Thieme et al. 1990; Yao et al. 2012) and this work has expanded on these studies to provide a statistical survey of individual structures on a wide range of scales. Note that because the magnetic pressure decreased and thermal pressure increased within the structures, all of them had a high plasmabeta compared to their surroundings; this is in contrast to other pressure balanced structures with a low plasma beta which are also commonly observed in the solar wind, and often contain magnetic flux ropes (Moldwin et al. 2000; Cartwright & Moldwin 2008; Kilpua et al. 2012; Yu et al. 2014; Feng et al. 2015). Simulations have recently shown that pressure balanced structures can develop in a turbulent plasma such as the solar wind (Yang et al. 2017). Dense structures have been tracked all the way from the surface of the Sun to 0.5 AU in white light images (Sheeley & Rouillard 2010), which suggests that the most visible structures analysed here have a solar origin as opposed to developing as part of the solar wind turbulence. This does not rule out the possibility that the 70% of structures that would not have been visible in white light images were formed in-transit as part of the solar wind turbulence.

Pressure balance across the structures means that no net force was acting on the them to alter their shape or size at the time of observation. We therefore predict that the only way the structures changed shape as they propagated away from the Sun was by undergoing a simple 1∕r expansion in the two directions perpendicular to the radial direction of flow. In addition this validates the naming of the areas as “structures” which simply travel with the bulk solar wind flow, and agrees with remote sensing observations that came to a similar conclusion (Viall et al. 2010). In contrast to remote sensing observations, the range of observed structure sizes here (50 Mm to 2000 Mm) contains scales much lower than the previous remote sensing lower limit of 1000 Mm (Viall et al. 2010). We have shown quantitively that this is because the smallest structures are too small to be resolved by white light images taken from 1 AU.

The transient nature of the structures suggests that their release into the heliosphere was not continuous, and that (unlike coronal hole plasma) they did not originate on field lines that have a continuous connection from the corona to the solar wind. This conclusion leads to the interpretation that the structures were produced by interchange reconnection allowing hot, dense plasma originating on initially closed field lines to escape on to open field lines that map to the heliosphere (the mechanism sketched in Crooker et al. 2004, Fig. 5). The wide range of structure sizes implies that interchange reconnection was happening on a large range of physical length scales. Additionally all of the structures were observed in slow solar wind (|v p | < 500 km s−1).

In addition to streamers, closed field regions can also take the form of pseudo-streamers which are pairs of closed loops that separate open field lines of the same polarity, and therefore do not have an associated current sheet (Wang et al. 2007). The structures observed here were all located away from in-situ crossings of HCS (by selection), which means that if interchange reconnection produced these structures it may have happened near pseudo-streamers as opposed to streamers. However, with only single spacecraft in-situ measurements and no remote sensing observations it is hard to tell how close to the HCS the spacecraft was at any given time. Simulations of interchange reconnection at the edges of streamers have produced blobs with high number density and high-β which surround the HCS (Higginson & Lynch 2017). Further simulations could investigate whether structures can also be produced via interchange reconnection at pseudo-streamers, and any similarities and differences that could be used to distinguish between streamer and pseudo-streamer produced structures using in-situ measurements.

Finally, despite a search of the entire Helios mission only 140 examples of density structures away from crossings of the HCS were identified, constituting 40 h of the 3 917 h of slow solar wind (|v p | < 500 km s−1) Helios measured at distances < 0.5 AU. The structures presented here are rare, and did not make up the bulk of the slow solar wind. Similar blobs are frequently observed above streamers and surrounding the HCS, which we have not investigated here (Sheeley & Rouillard 2010; Rouillard et al. 2010a,b, 2011; Viall & Vourlidas 2015), and Sanchez-Diaz et al. (2017) estimate that this material forms 15% of the slow solar wind. This means that even if all of the structures presented in this paper were generated at the edges of pseudo streamers, structures generated near pseudo-streamers are much rarer than those generated near streamers.

5 Conclusion

We have presented the first in-situ statistical study of transient number density structures in the inner heliosphere. The structures had a high plasma beta, were hotter than their surroundings, pressure balanced, and occurred in the slow solar wind but away from in-situ crossings of the HCS. These results suggest that away from the heliospheric plasma sheet, a small fraction of slow solar wind is produced by interchange reconnection releasing hot, dense plasma from closed field lines on to open field lines and out into the heliosphere.

In the near future Solar Orbiter will provide the first solar wind composition measurements in the inner heliosphere, allowing diagnosis of the coronal source of the structures. In addition Parker Solar Probe will provide the first in-situ observations of the solar wind at 10 solar radii, and will measure the structures shortly after they are released from the Sun. In combination with remote imaging of the structures, these new measurements will unambiguously allow us to identify the solar source of the structures.

Acknowledgements

The authors acknowledge useful discussions with Chadi Salem, Lorenzo Matteini, and Alexis Roulliard, and thank the anonymous reviewer for comments that significantly improved discussion of the results. D. Stansby is supported by STFC studentship ST/N504336/1. T. S. Horbury is supported by STFC grant ST/N000692/1. The Helios plasma dataset used here and described in Sect. 2 is available online, along with the source code used to generate the dataset (Stansby 2017). Data was retrieved using HelioPy v0.4.2 (Stansby 2018a), processed using astropy v3.0 (Astropy Collaboration et al. 2018), and figures were produced using Matplotlib v2.1.1 (Hunter 2007; Droettboom et al. 2017).

Appendix A: White light intensity signal of a single structure

The total intensity received by a white light detector is (Howard & Tappin 2009) I=0z2 ne(z)G(z)dz\begin{equation*} I = \int_{0}^{\infty} z^{2} n_{e} \left ( z \right ) G \left ( z \right ) \mathrm{d}z \end{equation*}(A.1)

where the integral is taken along the line of sight, ne is the electron number density, and G(z)$G \left ( z \right )$ is the scattering expression. I0 is defined as the value of the integral when ne = n0 ∝ 1∕r2, where r is radial distance from the centre of the Sun (i.e. an unperturbed number density profile for the solar wind). Because the integral is linear in ne , the total intensity when a structure is added to the background profile is I=I0+z0l/2z0+l/2z2 [ ne(z)n0(z) ]G(z)dz\begin{equation*} I = I_{0} &#x002B; \int_{z_{0} - l / 2}^{z_{0} &#x002B; l / 2} z^{2} \left [ n_{e}\left ( z \right ) - n_{0}\left ( z \right ) \right ] G \left ( z \right ) \mathrm{d}z \end{equation*}(A.2)

where the structure is centred at z0 and has a line-of-sight length l. Under the assumption that the value of the integrand is constant over the structure this further reduces to I=I0+lz02[ ne(z0)n0(z0) ]G(z0)\begin{equation*} I = I_{0} &#x002B; l z_{0}^{2}\left [ n_{e}\left ( z_{0} \right ) - n_{0}\left ( z_{0} \right ) \right ] G \left ( z_{0} \right ) \end{equation*}(A.3)

The change in intensity relative to the background intensity is then given by δII0=[ z02G(z0)n0(z0)I0 ]lδnn0\begin{equation*} \frac{\delta I}{I_{0}}= \left [\frac{z_{0}^{2} G \left ( z_{0} \right ) n_{0} \left ( z_{0} \right )}{I_{0}} \right ] l\ \frac{\delta n}{n_{0}} \end{equation*}(A.4)

where δI = II0 and δn = nen0. The term insquare brackets has units of inverse distance. It depends explicitly on the structure location (z0 ) and implicitly on the observer location and look direction. It takes its largest value when z0 is located on the Thompson sphere (Howard & Tappin 2009). For an observer situated at 1 AU, small observation angles with respect to the Sun, and a structure located on the Thompson sphere, the expression in square brackets evaluates numerically to ≈ 0.64∕r0, where r0 is the distancefrom the centre of the Sun to the number density structure. This gives the final expression δII00.64lr0δnn0.\begin{equation*} \frac{\delta I}{I_{0}} \approx 0.64 \frac{l}{r_{0}} \frac{\delta n}{n_{0}} .\end{equation*}(A.5)

This equation gives the change in intensity relative to the background intensity if a structure is located at the most favourable position for viewing. Constant contours of δII0 calculated using Eq. (A.5) are plotted on Fig. 3 for the range of structures sizes and density increases observed in this paper.

References

  1. Abbo, L., Ofman, L., Antiochos, S. K., et al. 2016, Space Sci. Rev., 201, 55 [NASA ADS] [CrossRef] [Google Scholar]
  2. Antiochos, S. K., Mikić, Z., Titov, V. S., Lionello, R., & Linker, J. A. 2011, ApJ, 731, 112 [Google Scholar]
  3. Antiochos, S. K., Linker, J. A., Lionello, R., et al. 2012, Space Sci. Rev., 172, 169 [NASA ADS] [CrossRef] [Google Scholar]
  4. Astropy Collaboration (Price-Whelan, A., et al.) 2018, ArXiv e-prints [arXiv:1801.02634] [Google Scholar]
  5. Burlaga, L. F., & Ogilvie, K. W. 1970, Sol. Phys., 15, 61 [NASA ADS] [CrossRef] [Google Scholar]
  6. Cartwright, M. L., & Moldwin, M. B. 2008, J. Geophys. Res.: Space Phys., 113, A09105 [NASA ADS] [CrossRef] [Google Scholar]
  7. Cranmer, S. R., van Ballegooijen, A. A., & Edgar, R. J. 2007, ApJS, 171, 520 [Google Scholar]
  8. Crooker, N. U., Gosling, J. T., & Kahler, S. W. 2002, J. Geophys. Res.: Space Phys., 107, 1028 [NASA ADS] [CrossRef] [Google Scholar]
  9. Crooker, N. U., Huang, C.-L., Lamassa, S. M., et al. 2004, J. Geophys. Res., 109, A03107 [NASA ADS] [CrossRef] [Google Scholar]
  10. DeForest, C. E., Howard, T. A., & Tappin, S. J. 2011, ApJ, 738, 103 [Google Scholar]
  11. Droettboom, M., Caswell, T. A., Hunter, J., et al. 2017, https://zenodo.org/record/1098480 [Google Scholar]
  12. Feng, H. Q., Zhao, G. Q., & Wang, J. M. 2015, J. Geophys. Res.: Space Phys., 120, 175 [Google Scholar]
  13. Fisk, L. A. 2003, J. Geophys. Res., 108, 1157 [Google Scholar]
  14. Fisk, L. A., & Schwadron, N. A. 2001, ApJ, 560, 425 [NASA ADS] [CrossRef] [Google Scholar]
  15. Fisk, L. A., & Zhao, L. 2008, Proc. Int. Astron. Union, 4, 109 [CrossRef] [Google Scholar]
  16. Gosling,J. T., Eriksson, S., & Schwenn, R. 2006, J. Geophys. Res., 111, A10102 [Google Scholar]
  17. Gosling, J. T., McComas, D. J., Roberts, D. A., & Skoug, R. M. 2009, ApJ, 695, L213 [NASA ADS] [CrossRef] [Google Scholar]
  18. Harra, L. K., Sakao, T., Mandrini, C. H., et al. 2008, ApJ, 676, L147 [NASA ADS] [CrossRef] [Google Scholar]
  19. Harra, L. K., Ugarte-Urra, I., De Rosa, M., et al. 2017, PASJ, 69, 47 [NASA ADS] [CrossRef] [Google Scholar]
  20. He, J., Pei, Z., Wang, L., et al. 2015, ApJ, 805, 176 [NASA ADS] [CrossRef] [Google Scholar]
  21. He, J.-S., Marsch, E., Tu, C.-Y., Guo, L.-J., & Tian, H. 2010, A&A, 516, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Higginson, A. K., & Lynch, B. J. 2017, ApJ, 859, 6 [NASA ADS] [CrossRef] [Google Scholar]
  23. Higginson, A. K., Antiochos, S. K., DeVore, C. R., Wyper, P. F., & Zurbuchen, T. H. 2017, ApJ, 837, 113 [NASA ADS] [CrossRef] [Google Scholar]
  24. Howard, T. A., & Tappin, S. J. 2009, Space Sci. Rev., 147, 31 [Google Scholar]
  25. Howard, R. A., Moses, J. D., Vourlidas, A., et al. 2008, Space Sci. Rev., 136, 67 [NASA ADS] [CrossRef] [Google Scholar]
  26. Howard, T. A., DeForest, C. E., & Reinard, A. A. 2012, ApJ, 754, 102 [NASA ADS] [CrossRef] [Google Scholar]
  27. Hunter, J. D. 2007, Comput. Sci. Eng., 9, 90 [Google Scholar]
  28. Kasper, J. C., Stevens, M. L., Lazarus, A. J., Steinberg, J. T., & Ogilvie, K. W. 2007, ApJ, 660, 901 [NASA ADS] [CrossRef] [Google Scholar]
  29. Kepko, L., Viall, N. M., Antiochos, S. K., et al. 2016, Geophys. Res. Lett., 43, 4089 [Google Scholar]
  30. Kilpua, E. K. J., Jian, L. K., Li, Y., Luhmann, J. G., & Russell, C. T. 2012, Sol. Phys., 281, 391 [NASA ADS] [Google Scholar]
  31. Leer, E., & Holzer, T. E. 1980, J. Geophys. Res.: Space Phys., 85, 4681 [NASA ADS] [CrossRef] [Google Scholar]
  32. Levine, R. H., Altschuler, M. D., & Harvey, J. W. 1977, J. Geophys. Res., 82, 1061 [NASA ADS] [CrossRef] [Google Scholar]
  33. Marsch, E., Mühlhäuser, K.-H., Rosenbauer, H., Schwenn, R., & Neubauer, F. M. 1982, J. Geophys. Res., 87, 35 [Google Scholar]
  34. Moldwin, M. B., Ford, S., Lepping, R., Slavin, J., & Szabo, A. 2000, Geophys. Res. Lett., 27, 57 [NASA ADS] [CrossRef] [Google Scholar]
  35. Musmann, G., Neubauer, F. M., Maier, A., & Lammers, E. 1975, Raumfahrtforschung, 19, 232 [NASA ADS] [Google Scholar]
  36. Owens, M. J., Lockwood, M., Riley, P., & Linker, J. 2017, J. Geophys. Res.: Space Phys. 122, 10980 [Google Scholar]
  37. Pinto, R. F., Brun, A. S., & Rouillard, A. P. 2016, A&A, 592, A65 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Porsche, H. 1977, J. Geophys., 42, 551 [Google Scholar]
  39. Rouillard, A. P., Davies, J. A., Lavraud, B., et al. 2010a, J. Geophys. Res.: Space Phys., 115, A04103 [NASA ADS] [Google Scholar]
  40. Rouillard, A. P., Lavraud, B., Davies, J. A., et al. 2010b, J. Geophys. Res.: Space Phys., 115, A04104 [NASA ADS] [Google Scholar]
  41. Rouillard, A. P., Sheeley, N. R., Cooper, T. J., et al. 2011, ApJ, 734, 7 [Google Scholar]
  42. Sanchez-Diaz, E., Rouillard, A. P., Davies, J. A., et al. 2017, ApJ, 851, 32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Scearce, C., Cantarano, S., Ness, N., et al. 1975, Raumfahrtforschung, 19, 237 [NASA ADS] [Google Scholar]
  44. Schwenn, R., Rosenbauer, H., & Miggenrieder, H. 1975, Raumfahrtforschung, 19, 226 [NASA ADS] [Google Scholar]
  45. Sheeley, N. R., & Rouillard, A. P. 2010, ApJ, 715, 300 [NASA ADS] [CrossRef] [Google Scholar]
  46. Sheeley, N. R., Harvey, J. W., & Feldman, W. C. 1976, Sol. Phys., 49, 271 [NASA ADS] [CrossRef] [Google Scholar]
  47. Sheeley, N. R., Wang, Y., Hawley, S. H., et al. 1997, ApJ, 484, 472 [NASA ADS] [CrossRef] [Google Scholar]
  48. Stansby, D. 2017, DOI: 10.5281/zenodo.1009506 [Google Scholar]
  49. Stansby, D. 2018a, DOI: 10.5281/zenodo.1166053 [Google Scholar]
  50. Stansby, D. 2018b, DOI: 10.5281/zenodo.1186586 [Google Scholar]
  51. Štverák, Š., Trávníček, P. M., & Hellinger, P. 2015, J. Geophys. Res.: Space Phys., 120, 8177 [NASA ADS] [CrossRef] [Google Scholar]
  52. Thieme, K., Marsch, E., & Schwenn, R. 1990, Ann. Geophys., 8, 713 [NASA ADS] [Google Scholar]
  53. Viall, N. M., & Vourlidas, A. 2015, ApJ, 807, 176 [NASA ADS] [CrossRef] [Google Scholar]
  54. Viall, N. M., Spence, H. E., Vourlidas, A., & Howard, R. 2010, Sol Phys., 267, 175 [Google Scholar]
  55. Wang, Y., Sheeley, N. R., & Rich, N. B. 2007, ApJ, 658, 1340 [NASA ADS] [CrossRef] [Google Scholar]
  56. Wang, Y.-M. 2010, ApJ, 715, L121 [Google Scholar]
  57. Wang, Y.-M., & Sheeley, N. R. 1991, ApJ, 372, L45 [Google Scholar]
  58. Wang, Y.-M., Sheeley, N. R., Walters, J. H., et al. 1998, ApJ, 498, L165 [NASA ADS] [CrossRef] [Google Scholar]
  59. Winterhalter, D., Smith, E. J., Burton, M. E., Murphy, N., & McComas, D. J. 1994, J. Geophys. Res., 99, 6667 [NASA ADS] [CrossRef] [Google Scholar]
  60. Yang, L., He, J., Tu, C., et al. 2017, ApJ, 836, 69 [NASA ADS] [CrossRef] [Google Scholar]
  61. Yao, S., Tu, C., He, J., Wei, W., & Meng, X. 2012, Chin. Sci. Bull., 57, 1421 [CrossRef] [Google Scholar]
  62. Yu, W., Farrugia, C. J., Lugaz, N., et al. 2014, J. Geophys. Res.: Space Phys., 119, 689 [Google Scholar]

All Figures

thumbnail Fig. 1

An example number density structure. Top three panels: a 4 day overview of 10 min averaged radial proton bulk speed (vr, top axis), radial magnetic field (Br, middle axis), and proton number density (np, bottom axis). The HCS crossing is marked with a red dashed line between days 87 and 88. Bottom four panels: a shorter 1.5 h interval, with the structure present between the grey dashed lines. The top axis shows proton number density (np ), and the second axis proton plasma beta (βp). The third axis shows proton parallel temperature (orange, Tp) and proton perpendicular temperature (blue, Tp). The fourth axis shows the proton thermal pressure (red, pp), magnetic pressure (blue, pmag), and total pressure (grey, ptot). The range of values in the total pressure is due to uncertainty in the number density and temperature of electrons and alpha particles.

In the text
thumbnail Fig. 2

Scatter plots of quantities inside the structures (y-axis) against quantities outside the structures (x-axis). One-to-one line are drawn in black on each plot. Points lying above these lines indicate a quantity larger inside a structure relative to outside and vice versa. Top left: proton core number density (np ) and top right total proton temperature (Tp). Middle left: magnetic field magnitude (| B |$\left | \mathbf{B} \right |$) and middle right magnetic (blue) proton thermal (red) and total (black) pressures. The total pressure crosses are error-bars in both directions due to uncertainty in the electron and alpha thermal pressure contribution to the total pressure. Bottom left: radial proton velocity (vr ) and bottom right proton temperature anisotropy (TpTp).

In the text
thumbnail Fig. 3

Scatter plot of electron density contrast (y-axis) against advected structure size (x-axis) . Uncertainties in electron density contrast are due to not knowing the variation of the alpha particle abundance within the structures. Background shading shows contours of Eq. (2) with r0 = 0.3 AU, which gives the fractional change in intensity a white-light instrument at 1 AU would see when viewing the blobs. The fractional intensity change values assume zero velocity smear across the observing instrument.

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.