The termination shock in a striped pulsar wind

I present observational and theoretical evidence that most of the pulsar spin-down energy is transferred away as a striped pulsar wind and that this energy is released by annihilation of the alternating magnetic field at the pulsar wind termination shock. One-dimensional particle-in-cells (PIC) simulations show that the alternating fields do annihilate at the termination shock in a striped wind. The particle acceleration should be studied in multidimensional simulations. As a first step, I simulated driven annihilation of alternating fields undergoing compression by an external force. It is shown that that in the course of this process, a particle distribution function is formed, which resembles that observed in plerions.     

Challenges to cosmic ray modeling: From beyond the solar wind termination shock

Voyager 1 crossed the solar wind termination shock on December 16, 2004 at a distance of 94 AU from the Sun, to become the first spacecraft to explore the termination shock region and to enter the heliosheath, the final heliospheric frontier. By the end of 2006, Voyager 1 will be at ∼101 AU, with Voyager 2 at ∼81 AU and still approaching the termination shock. Both spacecraft have been observing the modulation of galactic and anomalous cosmic rays since their launch in 1977. The recent observations close to or inside the heliosheath have provided several interesting ‘surprises’ with subsequent theoretical and modeling challenges. Examples are: what does the modulation of galactic cosmic rays amount to in this region?; how do the anomalous cosmic rays get accelerated and modulated?; why are there ‘breaks’ in the power-law slopes of the spectra of accelerated particles? Several numerical models have been applied to most of these topics over the years and comprehensive global predictions have been made the past decade, thought to be based on reasonable assumptions about the termination shock and the heliosheath. Examples of these predictions and assumptions are concisely discussed within the context of the main observed features of cosmic rays in the vicinity of the termination shock, ending with a discussion of some of the issues and challenges to cosmic ray modeling in particular.     

The role of radial perpendicular diffusion and latitude dependent acceleration along the solar wind termination shock

A numerical model, based on Parker’s transport equation, describing the modulation of anomalous cosmic rays and containing diffusive shock acceleration is applied. The role of radial perpendicular diffusion at the solar wind termination shock, and as the dominant diffusion coefficient in the outer heliosphere, is studied, in particular the role it plays in the effectiveness of the acceleration of anomalous protons and helium when its latitude dependence is changed. It is found that the latitudinal enhancement of radial perpendicular diffusion towards the heliospheric poles and along the termination shock has a prominent effect on the acceleration of these particles. It results in a ‘break’ in the energy spectrum for anomalous protons at ∼6.0 MeV, causing the spectral index to change from E^−1.38 to E^−2.23, but for anomalous helium at ∼3.0 MeV, changing the spectral index from E^−1.38 to E^−2.30. When approaching the simulated TS, the changes in the modulated spectra as they unfold to a ‘steady’ power law shape at energies below 50 MeV are much less prominent as a function of radial distances when radial perpendicular diffusion is increased with heliolatitude.     

The heliospheric modulation of cosmic rays: Effects of a latitude dependent solar wind termination shock

Observations made with the two Voyager spacecraft confirmed that the solar wind decelerates to form the heliospheric termination shock and that it has begun its merger with the local interstellar medium. The compression ratio of this shock affects galactic cosmic rays when they enter the heliosphere. Hydrodynamic (HD) models show that the compression ratio can have a significant latitude dependence; with the largest value in the nose direction of the heliosphere, becoming significantly less towards the polar regions. The modulation effects of such large latitude dependence are studied, using a well-established numerical drift and shock modulation model. We focus on computing the modulated spectra for galactic protons with emphasis on the radial and polar gradients in the equatorial plane and at a polar angle of θ = 55°, corresponding to the heliolatitude of Voyager 1. Two sets of solutions are computed and compared each time; with and without a latitude dependence for the compression ratio. All computations are done for the two magnetic field polarity cycles assuming solar minimum conditions. Including the termination shock in the model allows the study of the re-acceleration of galactic protons in the outer heliosphere. We find that for the A < 0 polarity cycle the intensity between ∼200 MeV and ∼1 GeV in the vicinity of the shock in the heliospheric equatorial plane may exceed the local interstellar value specified at the heliopause. Unfortunately, at θ = 55°, the effect is reduced. This seems not possible during an A > 0 cycle because significant modulation is then predicted between the heliopause and the termination shock, depending on how strong global gradient and curvature drifts are in the heliosheath. The overall effect of the shock on galactic protons in the equatorial plane is to reduce the total modulation as a function of radial distance with respect to the interstellar spectrum. Making the compression ratio latitude dependent enhances these effects at energies E < 200 MeV in the equatorial plane. At larger heliolatitudes these effects are even more significant. The differences in the modulation between the two drift cycles are compelling when the compression ratio is made latitude dependent but at Earth this effect is insignificant. A general result is that the computed radial gradient changes for galactic protons at and close to the TS and that these changes are polarity dependent. In line with previous work, large polarity dependent effects are predicted for the inner heliosphere and also close to the shock’s position in the equatorial plane. In contrast, at θ = 55°, the largest polarity effect occurs in the middle heliosphere (50 AU), enhanced by the latitude dependence of the compression ratio. At this latitude, the amount of proton modulation between the heliopause and the termination shock is much reduced. If galactic cosmic rays were to experience some diffusive shock acceleration over the 100–1000 MeV range at the shock, the radial gradient should change its sign in the vicinity of the shock, how large, depends on the compression ratio and the amount of drifts taking place in the outer heliosphere. The effective polar gradient shows a strong polarity dependence at Earth but this dissipates at θ = 55°, especially with increasing radial distance. This tendency is enhanced by making the compression ratio latitude dependent.     

Cosmic ray reacceleration on the galactic wind termination shock

Reacceleration of cosmic rays produced by galactic sources on the galactic wind termination shock is considered. The problem of the cosmic ray spectrum continuity is investigated. Numeric results are presented and discussed. We found that a smooth spectral transition from the galactic cosmic rays to the cosmic rays reaccelerated at the galactic wind termination shock is difficult to produce, if the maximum energy of accelerated particles is the same throughout the surface of the termination shock. The possible solution of this problem is the non-spherical termination shock with different maximum energies at different places of the shock.     

Modulation of cosmic rays: Perpendicular diffusion and drifts in a heliosphere with a solar wind termination shock

In this study the roles of polar perpendicular diffusion and drifts are illustrated in a model containing a heliosheath and diffusive shock acceleration as applied to the solar wind termination shock. Of particular interest is the relation of polar perpendicular diffusion to particle drifts and how the effectiveness of the termination shock acceleration of galactic and anomalous protons is influenced by this relation. We found that drifts have a more prominent effect than the polar enhancement of perpendicular diffusion so that its omission from termination shock models would produce unrealistically large shock acceleration and consequently also larger modulation effects throughout the heliosphere. The computed spectra at a heliolatitude of 35° are almost similar for the two polarity epochs indicating that the two Voyager spacecraft might not observe differences between the two cycles in future.     

Effects of the position of the solar wind termination shock and the heliopause on the heliospheric modulation of cosmic rays

The effects of changing the position of the solar wind termination shock and the position of the heliopause, and therefore the extent of the heliosheath, on the modulation of cosmic ray protons are illustrated. An improved numerical model with diffusive termination shock acceleration, a heliosheath and drifts is used. The modulation is computed in the equatorial plane and at 35 heliolatitude using recently derived diffusion coefficients applicable to a number of cosmic ray species during both magnetic polarity cycles of the Sun. It was found that qualitatively the modulation results for the different heliopause positions are similar although they differ quantitatively, e.g., clearly different radial gradients are predicted for the regions beyond the termination shock compared to inside the shock. The difference between the modulation for the two solar polarity cycles are less significant at a heliolatitude of 35° than in the equatorial plane. We found that moving the termination shock from 90 to 100 AU, with the heliopause fixed at 120 AU, caused only quantitative differences so that the exact position of the TS in the outer heliosphere seems not crucially important to global modulation. Moving the heliopause outwards, to represent the modulation in the tail region of the heliosphere, causes overall decreases in the cosmic ray intensities but not linearly as a function of energy, e.g., at 1 GeV the effect is insignificant. We conclude from this modelling that the modulation of protons in the heliospheric nose and tail regions are qualitatively similar although, clear quantitative and interesting differences occur.     

The nature of the distant solar wind

This review emphasizes the large scale, global aspects of the solar wind flow at large heliocentric distances and how that structure derives from coronal and interplanetary phenomena. Particular attention is paid to the solar cycle variation, and speculation on its three-dimensional manifestations in the outer heliosphere is offered. Solar cycle variations in the global solar wind structure may lead to systematic variations in the dynamical properties of the interplanetary medium which may have a significant effect upon galactic cosmic ray modulation and the three-dimensional shape of the heliospheric termination shock.     

Sharp boundaries of solar wind plasma structures and their relationship to solar wind turbulence

Sharp (<10 min) and large (>20%) solar wind ion flux changes are common phenomena in turbulent solar wind plasma. These changes are the boundaries of small- and middle-scale solar wind plasma structures which can have a significant influence on Earth’s magnetosphere. These solar wind ion flux changes are typically accompanied by only a small change in the bulk solar wind velocity, hence, the flux changes are driven mainly by plasma density variations. We show that these events occur more frequently in high-density solar wind. A characteristic of solar wind turbulence, intermittency, is determined for time periods with and without these flux changes. The probability distribution functions (PDF) of solar wind ion flux variations for different time scales are calculated for each of these periods and compared. For large time scales, the PDFs are Gaussian for both data sets. For small time scales, the PDFs from both data set are more flat than Gaussian, but the degree of flatness is much larger for the data near the sharp flux change boundaries.     

Acceleration at the termination shock: Anisotropies and spectra

Low-energy termination shock particle populations observed by the Voyagers upstream of the shock exhibited strong field-aligned beaming with anisotropies of the order of unity. The Parker transport equation is valid only for nearly isotropic phase space distributions and is inapplicable to these highly beamed populations. The usual approach is to revert to the more general focused transport equation retaining pitch-angle information. We developed a complimentary technique employing a three-moment expansion of the Skilling equation using Legendre polynomials. We investigate the effects of adiabatic focusing and reflection on the diffusive acceleration process at oblique shock waves. It is shown that low-energy particle intensities are discontinuous and sharply peaked at the shock, consistent with the observations. Particle spectra are not only harder than the power laws predicted from diffusive transport theory, but also exhibit spectral gaps near the low-energy acceleration threshold due to more efficient acceleration by scattering and mirroring. Our model also predicts upstream anisotropies as high as 100% for highly oblique shocks whereas downstream distributions are nearly isotropic.     

Correlations between neutral and ionized solar wind

We report results of a statistical study correlating ionized solar wind (ISW) fluxes observed by ACE during late 2000 and throughout 2001 with neutral solar wind (NSW) fluxes observed by IMAGE/LENA over the same period. The average correlation coefficient between the neutral and ionized solar wind is 0.66 with correlations greater than 0.80 occurring about 29% of the time. Correlations appear to be driven by high solar wind flux variability, similar to results obtained by in situ multi-spacecraft correlation studies. In this study, however, IMAGE remains inside the magnetosphere on over 95% of its orbits. As a function of day of year, or equivalently ecliptic longitude, the slope of the relationship between the neutral solar wind flux and the ionized solar wind flux shows an enhancement near the upstream direction, but the symmetry point appears shifted toward higher ecliptic longitudes than the interstellar neutral (ISN) flow direction by about 20°. The estimated peak interstellar neutral upstream density inside of 1 AU is about 7 × 10⁻³ cm⁻³.     

Atmospheric electrodynamics modulated by the solar wind

Variations of stratospheric temperature are connected with changes of the solar wind dynamic pressure. This effect could be explained in the framework of the global electric circuit concept. The energy of the solar wind modulates the energy balance of the global electric circuit where the stratosphere could be one of its other elements. The conductivity of the stratosphere in the polar region is equal to and sometimes more than the conductivity of the ground surface covered by ice or permafrost. Re-distribution of the global electric circuit currents between the stratosphere and the ground surface determines a different relation between solar wind dynamics and variations of the stratospheric temperature during different seasons.     

Chaos and multifractals in the solar wind

By using the false-nearest-neighbours method, we have argued that the deterministic component of solar wind plasma dynamics should be low-dimensional. In fact, the results we have obtained using the method of topological embedding indicate that the behaviour of the solar wind can be approximately described by a low-dimensional chaotic attractor in the inertial manifold, which is a subspace of system phase space. We have also shown that the multifractal spectrum of the solar wind attractor is consistent with that for the multifractal measure of the self-similar generalized weighted Cantor set with two different scaling parameters and one probability measure parameter responsible for nonuniform compression in phase space and multifractality. The values of the parameters fitted also demonstrate that the complex solar wind system could only be weakly non-conservative (small dissipation) and quantify nonlinear dynamics; some parts of the attractor in phase space are visited much more frequently than other parts. In addition, to quantify the multifractality of space plasma intermittent turbulence, we consider that generalized Cantor set also in the context of scaling properties of solar wind turbulence. We investigate the resulting multifractal spectrum of a one-dimensional phenomenological model of turbulence cascade depending on its parameters, especially for asymmetric scaling. In particular, we have shown that intermittent pulses are stronger for the cascade model with two different scaling parameters. Even thought solar wind turbulence appears to be rather space filling, a better agreement with the data is obtained, especially for the negative index of generalized dimensions. Therefore we argue that there is a need to use a two-scale asymmetric cascade model. We hope that this generalized multifractal model will be a useful tool for analysis of intermittent turbulence in space plasmas. We thus believe that fractal analysis of chaotic systems could lead us to a deeper understanding of their nature, and maybe even to predict their seemingly unpredictable behaviour.     

A new approach to solar wind monitoring

The monitoring of solar wind parameters is a key problem of the space weather program. We are presenting a new solution of plasma parameter determination suitable for small and fast solar wind monitors. The first version will be launched during the SPECTR-R project into a highly elongated orbit with apogee ∼350,000 km. The method is based on simultaneous measurements of the total ion flux and ion integral energy spectrum by six identical Faraday cups. Three of them are dedicated to determination of the ion flow direction, the other three (equipped with control grids supplied by a retarding potential) are used for determination of the density, temperature, and speed of the plasma flow. The version under development is primarily designed for the measurements in the solar wind and tail magnetosheath, thus for velocities range from 270 to 750 km/s, temperatures from 1 to 30 eV, and densities up to 200 cm⁻³. However, the instrument design can be simply modified for measurements in other regions with a substantial portion of low-energy plasma as a subsolar magnetosheath, cusp or low-latitude boundary layer. Testing of the engineering model shows that the proposed method can provide reliable plasma parameters with a high time resolution (up to 8 Hz). The paper presents not only the method and its technical realization but it documents all advantages and peculiarities of the suggested approach.     

Geomagnetic activity driven by solar wind turbulence

Within the framework of the solar wind—magnetosphere coupled system, intense perturbations in the solar wind, causing geomagnetic storms and substorms, have been widely studied by means of the so-called coupling parameters. However, remarkable variations in the geomagnetic field occur even in absence of such perturbations. In those conditions, solar wind MHD turbulence might have a role. Recent results have shown that solar wind turbulence can be described not only as a mixture of inward and outward stochastic Alfvénic fluctuations, but includes also advected structures, dominated by an excess of magnetic energy.     

Self-similar stochastic processes in solar wind turbulence

Solar wind data is used to estimate the autocorrelation function for the stochastic process x(τ) = y(t + τ) − y(t), considered as a function of τ, where y(t) is any one of the quantities B²(t), n_p(t)V²(t), or n_p(t). This process has stationary increments and a variance that increases like a power law τ^2γ where γ is the scaling exponent. For the kinetic energy density and the proton density the scaling exponent is close to the Kolmogorov value γ = 1/3, for the magnetic energy density it is slightly larger. In all three cases, it is shown that the autocorrelation function estimated from the data agrees with the theoretical autocorrelation function for a self-similar stochastic process with stationary increments and finite variance. This is far from proof, but it suggests that these stochastic processes may be self-similar for time scales in the small scale inertial range of the turbulence, that is, from approximately 10 to 10³ s.     

Flows and obstacles in the solar wind

Understanding the physics of the various disturbances in the solar wind is critical to successful forecasts of space weather. The STEREO mission promises to bring us new and deeper understanding of these disturbances. As we stand on the threshold of the first results from this mission, it is appropriate to review what we know about solar wind disturbances. Because of their complementary nature we discuss both the disturbances that arise within the solar wind due to the stream structure and coronal mass ejecta and the disturbances that arise when the solar wind collides with planetary obstacles, such as magnetospheres.     

A polytropic model for the solar wind

The solar wind fills the heliosphere and is the background medium in which coronal mass ejections propagate. A realistic modelling of the solar wind is therefore essential for space weather research and for reliable predictions. Although the solar wind is highly anisotropic, magnetohydrodynamic (MHD) models are able to reproduce the global, average solar wind characteristics rather well. The modern computer power makes it possible to perform full three dimensional (3D) simulations in domains extending beyond the Earth’s orbit, to include observationally driven boundary conditions, and to implement even more realistic physics in the equations. In general, MHD models for the solar wind often make use of additional source and sink terms in order to mimic the observed solar wind parameters and/or they hide the not-explicitly modelled physical processes in a reduced or variable adiabatic index. Even the models that try to take as much as possible physics into account, still need additional source terms and fine tuning of the parameters in order to produce realistic results. In this paper we present a new and simple polytropic model for the solar wind, incorporating data from the ACE spacecraft to set the model parameters. This approach allows to reproduce the different types of solar wind, where the simulated plasma variables are in good correspondence with the observed solar wind plasma near 1 AU.