Search

In mid-2012, a GMIR observed by Voyager 2 crossed through the heliosheath and collided with the heliopause, generating a pressure pulse that propagated into the very local interstellar medium. The effects of the transmitted wave were seen by Voyager 1 just 93 days after its own heliopause crossing. The passage of the transient was accompanied by long-lasting decreases in galactic cosmic ray intensities that occurred from ~2012.55 to ~2013.35 and ~2012.91 to ~2013.70 at Voyager 2 and Voyager 1, respectively. Omnidirectional (>20 MeV) proton-dominated measurements from each spacecraft's Cosmic Ray Subsystem reveal a remarkable similarity between these causally-related events, with a correlation coefficient of 91.2% and a time-lag of 130 days. Knowing the locations of the two spacecraft, we use the observed time-delay to calculate the GMIR's average speed through the heliosheath (inside the heliopause) as a function of temperature in the very local interstellar medium. This, combined with particle, field, and plasma observations enables us to infer previously unmeasured properties of the heliosheath, including a range of sound speeds and total effective pressures. For a nominal temperature of ~20,000 K just outside the heliopause, we find a sound speed of 314 (+/-) 32 km/s and total effective pressure of 267 (+/-) 55 fPa inside the heliopause. We compare these results with the Interstellar Boundary Explorer's data-driven models of heliosheath pressures derived from energetic neutral atom fluxes (the globally distributed flux) and present them as additional evidence that the heliosheath's dynamics are driven by suprathermal energetic processes.

It is currently believed that the turbulent fluctuations pervade the outermost heliosphere. Turbulence, magnetic reconnection, and their link may be responsible for magnetic energy conversion in these regions. The governing mechanisms of such anisotropic and compressible magnetic turbulence in the inner heliosheath (IHS) and in the local interstellar medium (LISM) still lack a thorough description. The present literature mainly concerns large scales which are not representative of the inertial-cascade dynamics of turbulence. Moreover, lack of broadband spectral analysis makes the IHS dynamics critically understudied. Our recent study shows that 48 s magnetic-field data from the Voyager mission are appropriate for a spectral analysis over a frequency range of six decades, from 5 x 10-8 Hz to 10-2 Hz. Here, focusing on the Voyager 2 observation interval from 2013.824 to 2016.0, we describe the structure of turbulence in a sector zone of the IHS. A spectral break around 7 x 10-7 Hz (magnetic structures with size l~1.3 Astronomical Units) separates the energy-injection regime from the inertial-cascade regime of turbulence. A second scale is observed around 6 x 10-5 Hz (l~ 0.017 AU) and corresponds to a peak of compressibility and intermittency of fluctuations., Comment: 10

The understanding of inertial-scale dynamics in the heliosheath is not yet thorough. Magnetic field fluctuations across the inner heliosheath (IHS) and the local interstellar medium (LISM) are here considered to provide accurate and highly resolved statistics over different plasma conditions between 88 and 136 au. By using the unique in situ 48 s measurements from the Voyager Interstellar Mission, we investigate different fluctuation regimes at the magnetohydrodynamic (MHD) scales, down to the MHD-to-kinetic transition. We focus on a range of scales exceeding five frequency decades (5 x 10-8 < f < 10-2 Hz), which is unprecedented in literature analysis. A set of magnetic field data for eight intervals in the IHS, in both unipolar and sector regions, and four intervals in the LISM is used for the analysis. Results are set forth in terms of the power spectral density, spectral compressibility, structure functions, and intermittency of magnetic field increments. In the heliosheath, we identify the energy-injection regime displaying a ~1/f energy decay, and the inertial-cascade regime. Here, the power spectrum is anisotropic and dominated by compressive modes, with intermittency that can reach kurtosis values of up to 10. In the interstellar medium the structure of turbulence is anisotropic as well, with transverse fluctuations clearly prevailing after 2015 May. Here, we show that intermittent features occur only at scales smaller than 10-6 Hz., Comment: 32 pages, 12 figures

The Voyager plasma observations show that the physics of the heliosheath is rather complex, and especially that the temperature derived from observation differs from expectations. To explain this fact the temperature in the heliosheath should be based on $\kappa$ distributions instead of Maxwellians because the former allows for much higher temperature. Here we show an easy way to calculate the $\kappa$ temperatures when those estimated from the data are given as Maxwellian temperatures. We use the moments of the Maxwellian and $\kappa$ distributions to estimate the $\kappa$ temperature. Moreover, we show that the pressure (temperature) given by a truncated $\kappa$ distribution is similar to that given by a Maxwellian and only starts to increase for higher truncation velocities. We deduce a simple formula to convert the Maxwellian to $\kappa$ pressure or temperature. We apply this result to the Voyager-2 observations in the heliosheath., Comment: 9 pages 7 Figures

Directly tracking an interplanetary coronal mass ejection (ICME) by widely separated spacecrafts is a great challenge. However, such an event could provide us a good opportunity to study the evolution of embedded Alfv\'enic fluctuations (AFs) inside ICME and their contributions to local plasma heating directly. In this study, an ICME observed by Wind at 1.0 au on March 4-6 1998 is tracked to the location of Ulysess at 5.4 au. AFs are commonly found inside the ICME at 1.0 au, with an occurrence rate of 21.7% and at broadband frequencies from 4$\times 10^{-4}$ to 5$\times 10^{-2}$ Hz. When the ICME propagates to 5.4 au, the Aflv\'enicity decreases significantly, and AFs are rare and only found at few localized frequencies with the occurrence rate decreasing to 3.0%. At the same time, the magnetic field intensity at the AF-rich region has an extra magnetic dissipation except ICME expansion effect. The energetics of the ICME at different radial distance is also investigated here. Under similar magnetic field intensity situations at 1.0 au, the turbulence cascade rate at the AF-rich region is much larger than the value at the AF-lack region. Moreover, it can maintain as the decrease of magnetic field intensity if there is lack of AFs. However, when there exists many AFs, it reduces significantly as the AFs disappear. The turbulence cascade dissipation rate within the ICME is inferred to be 1622.3 $J\cdot kg^{-1}\cdot s^{-1}$, which satisfies the requirement of local ICME plasma heating rate, 1653.2 $J\cdot kg^{-1}\cdot s^{-1}$. We suggest that AF dissipation is responsible for extra magnetic dissipation and local plasma heating inside ICME., Comment: 15 pages, 3 figures, 1 table

Nonlinear cascade of low-frequency Alfvenic fluctuations (AFs) is regarded as one candidate of the energy sources to heat plasma during the non-adiabatic expansion of interplanetary coronal mass ejections (ICMEs). However, AFs inside ICMEs were seldom reported in the literature. In this study, we investigate AFs inside ICMEs using observations from Voyager 2 between 1 and 6 au. It is found that AFs with high degree of Alfvenicity frequently occurred inside ICMEs, for almost all the identified ICMEs (30 out of 33 ICMEs), and 12.6% of ICME time interval. As ICMEs expand and move outward, the percentage of AF duration decays linearly in general. The occurrence rate of AFs inside ICMEs is much less than that in ambient solar wind, especially within 4 au. AFs inside ICMEs are more frequently presented in the center and at the boundaries of ICMEs. In addition, the proton temperature inside ICME has a similar distribution. These findings suggest significant contribution of AFs on local plasma heating inside ICMEs.

As a follow-up study on Sun-to-Earth propagation of fast coronal mass ejections (CMEs), we examine the Sun-to-Earth characteristics of slow CMEs combining heliospheric imaging and in situ observations. Three events of particular interest, the 2010 June 16, 2011 March 25 and 2012 September 25 CMEs, are selected for this study. We compare slow CMEs with fast and intermediate-speed events, and obtain key results complementing the attempt of \citet{liu13} to create a general picture of CME Sun-to-Earth propagation: (1) the Sun-to-Earth propagation of a typical slow CME can be approximately described by two phases, a gradual acceleration out to about 20-30 solar radii, followed by a nearly invariant speed around the average solar wind level, (2) comparison between different types of CMEs indicates that faster CMEs tend to accelerate and decelerate more rapidly and have shorter cessation distances for the acceleration and deceleration, (3) both intermediate-speed and slow CMEs would have a speed comparable to the average solar wind level before reaching 1 AU, (4) slow CMEs have a high potential to interact with other solar wind structures in the Sun-Earth space due to their slow motion, providing critical ingredients to enhance space weather, and (5) the slow CMEs studied here lack strong magnetic fields at the Earth but tend to preserve a flux-rope structure with axis generally perpendicular to the radial direction from the Sun. We also suggest a "best" strategy for the application of a triangulation concept in determining CME Sun-to-Earth kinematics, which helps to clarify confusions about CME geometry assumptions in the triangulation and to improve CME analysis and observations., Comment: 37 pages, 13 figures, accepted for publication in ApJ Supplement

The Voyager probes are the furthest, still active, spacecraft ever launched from Earth. During their 38-year trip, they have collected data regarding solar wind properties (such as the plasma velocity and magnetic field intensity). Unfortunately, a complete time evolution of the measured physical quantities is not available. The time series contains many gaps which increase in frequency and duration at larger distances. The aim of this work is to perform a spectral and statistical analysis of the solar wind plasma velocity and magnetic field using Voyager 2 data measured in 1979, when the gaps/signal ratio is of order of unity. This analysis is achieved using four different data reconstruction techniques: averages on linearly interpolated subsets, correlation of linearly interpolated data, compressed sensing spectral estimation, and maximum likelihood data reconstruction. With five frequency decades, the spectra we obtained have the largest frequency range ever computed at 5 astronomical units from the Sun; spectral exponents have been determined for all the components of the velocity and magnetic field fluctuations. Void analysis is also useful in recovering other spectral properties such as integral scales (see for instance Table 4) and, if the confidence level of the measurements is sufficiently high, the decay variation in the small scale range due, for instance, to dissipative effects., Comment: 11 pages, 7 figures

The largest geomagnetic storms of solar cycle 24 so far occurred on 2015 March 17 and June 22 with $D_{\rm st}$ minima of $-223$ and $-195$ nT, respectively. Both of the geomagnetic storms show a multi-step development. We examine the plasma and magnetic field characteristics of the driving coronal mass ejections (CMEs) in connection with the development of the geomagnetic storms. A particular effort is to reconstruct the in situ structure using a Grad-Shafranov technique and compare the reconstruction results with solar observations, which gives a larger spatial perspective of the source conditions than one-dimensional in situ measurements. Key results are obtained concerning how the plasma and magnetic field characteristics of CMEs control the geomagnetic storm intensity and variability: (1) a sheath-ejecta-ejecta mechanism and a sheath-sheath-ejecta scenario are proposed for the multi-step development of the 2015 March 17 and June 22 geomagnetic storms, respectively; (2) two contrasting cases of how the CME flux-rope characteristics generate intense geomagnetic storms are found, which indicates that a southward flux-rope orientation is not a necessity for a strong geomagnetic storm; and (3) the unexpected 2015 March 17 intense geomagnetic storm resulted from the interaction between two successive CMEs plus the compression by a high-speed stream from behind, which is essentially the "perfect storm" scenario proposed by \citet[][i.e., a combination of circumstances results in an event of unusual magnitude]{liu14a}, so the "perfect storm" scenario may not be as rare as the phrase implies., Comment: 15 pages, 5 figures, accepted for publication in ApJ Letters

The transition between the supersonic solar wind and the subsonic heliosheath, the termination shock (TS), was observed by Voyager 2 (V2) on 2007 August 31-September 1 at a distance of 84 AU from the Sun. The data reveal multiple crossings of a complex, quasi-perpendicular supercritical shock. These experimental data are the starting point for a more sophisticated analysis that includes computer modeling of a shock in the presence of pickup ions (PUIs). here, we present two-dimensional (2-D) particle-in-cell (PIC) simulations of the TS including PUIs self-consistently. We also report the ion velocity distribution across the TS using the Faraday cup data from V2. A relatively complete plasma and magnetic field data set from V2 gives us the opportunity to do a full comparison between the experimental data and PIC simulation results. Our results show that: (1) The nonstationarity of the shock front is mainly caused by the ripples along the shock front and these ripples from even if the percentage of PUIs is high. (2) PUIs play a key role in the energy dissipation of the TS, and most of the incident ion dynamic energy is transferred to the thermal energy of PUIs instead of solar wind ions (SWIs). (3) The simulated composite heliosheath ion velocity distribution function is a superposition of a cold core formed by transmitted SWIs, the shoulders contributed by the hot reflected SWIs and directly transmitted PUIs, and the wings of the distribution dominated by the very hot reflected PUIs. (4) The V2 Faraday cups observed the cool core of the distribution, so they saw only a tip of the iceberg. For the evolution of the cool core distribution function across the TS, the computed results agree reasonably well with the V2experimental results., Comment: The paper has been accepted by ApJ

In the majority of the literature on plasma shock waves, electrons play the role of "ghost particles," since their contribution to mass and momentum flows is negligible, and they have been treated as only taking care of the electric plasma neutrality. In some more recent papers, however, electrons play a new important role in the shock dynamics and thermodynamics, especially at the solar-wind termination shock. They react on the shock electric field in a very specific way, leading to suprathermal nonequilibrium distributions of the downstream electrons, which can be represented by a kappa distribution function. In this paper, we discuss why this anticipated hot electron population has not been seen by the plasma detectors of the Voyager spacecraft downstream of the solar-wind termination shock. We show that hot nonequilibrium electrons induce a strong negative electric charge-up of any spacecraft cruising through this downstream plasma environment. This charge reduces electron fluxes at the spacecraft detectors to nondetectable intensities. Furthermore, we show that the Debye length $\lambda _{\mathrm D}^{\kappa}$ grows to values of about $\lambda _{\mathrm D}^{\kappa}/\lambda _{\mathrm D}\simeq 10^{6}$ compared to the classical value $\lambda _{\mathrm D}$ in this hot-electron environment. This unusual condition allows for the propagation of a certain type of electrostatic plasma waves that, at very large wavelengths, allow us to determine the effective temperature of the suprathermal electrons directly by means of the phase velocity of these waves. At moderate wavelengths, the electron-acoustic dispersion relation leads to nonpropagating oscillations with the ion-plasma frequency $\omega _{\mathrm p}$ , instead of the traditional electron plasma frequency., Comment: 6 pages, 2 figures

On 2012 September 30 - October 1 the Earth underwent a two-step geomagnetic storm. We examine the Sun-to-Earth characteristics of the coronal mass ejections (CMEs) responsible for the geomagnetic storm with combined heliospheric imaging and in situ observations. The first CME, which occurred on 2012 September 25, is a slow event and shows an acceleration followed by a nearly invariant speed in the whole Sun-Earth space. The second event, launched from the Sun on 2012 September 27, exhibits a quick acceleration, then a rapid deceleration and finally a nearly constant speed, a typical Sun-to-Earth propagation profile for fast CMEs \citep{liu13}. These two CMEs interacted near 1 AU as predicted by the heliospheric imaging observations and formed a complex ejecta observed at Wind, with a shock inside that enhanced the pre-existing southward magnetic field. Reconstruction of the complex ejecta with the in situ data indicates an overall left-handed flux rope-like configuration, with an embedded concave-outward shock front, a maximum magnetic field strength deviating from the flux rope axis and convex-outward field lines ahead of the shock. While the reconstruction results are consistent with the picture of CME-CME interactions, a magnetic cloud-like structure without clear signs of CME interactions \citep{lugaz14} is anticipated when the merging process is finished., Comment: 15 pages, 5 figures. Accepted for publication in ApJ Letters

In 2012 March the Sun exhibited extraordinary activities. In particular, the active region NOAA AR 11429 emitted a series of large coronal mass ejections (CMEs) which were imaged by STEREO as it rotated with the Sun from the east to west. These sustained eruptions are expected to generate a global shell of disturbed material sweeping through the heliosphere. A cluster of shocks and interplanetary CMEs (ICMEs) were observed near the Earth, and are propagated outward from 1 AU using an MHD model. The transient streams interact with each other, which erases memory of the source and results in a large merged interaction region (MIR) with a preceding shock. The MHD model predicts that the shock and MIR would reach 120 AU around 2013 April 22, which agrees well with the period of radio emissions and the time of a transient disturbance in galactic cosmic rays detected by Voyager 1. These results are important for understanding the "fate" of CMEs in the outer heliosphere and provide confidence that the heliopause is located around 120 AU from the Sun., Comment: 13 pages, 5 figures, accepted for publication in ApJ Letters