Modeling the acceleration and modulation of anomalous cosmic ray oxygen

After the solar wind termination shock crossings of the Voyager spacecraft, the acceleration of anomalous cosmic rays has become a very contentious subject. In this paper we examine several topics pertinent to anomalous cosmic ray oxygen acceleration and transport using a numerical cosmic ray modulation model. These include the effects of drifts on a purely Fermi I accelerated spectra, the effects of introducing higher charge states of oxygen into the modulation model, examining the viability of momentum diffusion as a re-acceleration process in the heliosheath and examining energy spectra, and intensity gradients, in the inner heliosphere during consecutive drift cycles.     

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.     

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.     

Acceleration of interplanetary pick-up ions and anomalous cosmic rays

It may not be doubted anymore that anomalous cosmic rays (ACRs) are produced in the heliosphere from interplanetary pick-up ions through their acceleration at the solar wind termination shock. However, there is no general agreement in the community of heliospheric researchers concerning the mechanism of injection of the pick-up ions into the shock acceleration. We discuss here three possible ways for pick-up ions to be involved into the acceleration process at the termination shock: (1) preacceleration of pick-up ions in the whole region from the Sun up to the termination shock by solar wind turbulences and interplanetary shock waves, (2) local preacceleration of pick-up ions in a vicinity of the termination shock by shock surfing, and (3) formation of high-velocity tails in pick-up ion spectra consisting of secondary pick-up ions which are produced in the supersonic solar wind due to ionization of energetic neutral atoms entering from the inner heliosheath.     

Modulation of galactic cosmic rays in a north–south asymmetrical heliosphere

Observations made with the two Voyager spacecraft confirmed that the solar wind decelerates to form the heliospheric termination shock. Voyager 1 crossed this termination shock at ∼94 AU in 2004, while Voyager 2 crossed it in 2007 at a different heliolatitude, about 10 AU closer to the Sun. These different positions of the termination shock confirm the dynamic and cyclic nature of the shock’s position. Observations from the two Voyager spacecraft inside the heliosheath indicate significant differences between them, suggesting that apart from the dynamic nature caused by changing solar activity there also may exist a global asymmetry in the north–south (polar) dimensions of the heliosphere, in addition to the expected nose–tail asymmetry. This relates to the direction in which the heliosphere is moving in interstellar space and its orientation with respect to the interstellar magnetic field. In this paper we focus on illustrating the effects of this north–south asymmetry on modulation of galactic cosmic ray Carbon, between polar angles of 55° and 125°, using a numerical model which includes all four major modulation processes, the termination shock and the heliosheath. This asymmetry is incorporated in the model by assuming a significant dependence on heliolatitude of the thickness of the heliosheath. When comparing the computed spectra between the two polar angles, we find that at energies E < ∼1.0 GeV the effects of the assumed asymmetry on the modulated spectra are insignificant up to 60 AU from the Sun but become increasingly more significant with larger radial distances to reach a maximum inside the heliosheath. In contrast, with E > ∼1.0 GeV, these effects remain insignificant throughout the heliosphere even very close to the heliopause. Furthermore, we find that a higher local interstellar spectrum for Carbon enhances the effects of asymmetric modulation between the two polar angles at lower energies (E < ∼300 MeV). In conclusion, it is found that north–south asymmetrical effects on the modulation of cosmic ray Carbon depend strongly on the extent of the geometrical asymmetry of the heliosheath together with the assumed value of the local interstellar spectrum.     

Modelling of galactic Carbon in an asymmetrical heliosphere: Effects of asymmetrical modulation conditions

Observations of galactic cosmic rays (GCRs) from the two Voyager spacecraft inside the heliosheath indicate significant differences between them, suggesting that in addition to a possible global asymmetry in the north–south dimensions (meridional plane) of the heliosphere, it is also possible that different modulation (turbulence) conditions could exist between the two hemispheres of the heliosphere. We focus on illustrating the effects on GCR Carbon of asymmetrical modulation conditions combined with a heliosheath thickness that has a significant dependence on heliolatitude. To reflect different modulation conditions between the two heliospheric hemispheres in our numerical model, the enhancement of both polar and radial perpendicular diffusion off the ecliptic plane is assumed to differ from heliographic pole to pole. The computed radial GCR intensities at polar angles of 55° (approximating the Voyager 1 direction) and 125° (approximating the Voyager 2 direction) are compared at different energies and for both particle drift cycles. This is done in the context of illustrating how different values of the enhancement of both polar and radial perpendicular diffusion between the two hemispheres contribute to causing differences in radial intensities during solar minimum and moderate maximum conditions. We find that in the A > 0 cycle these differences between 55° and 125° change both quantitatively and qualitatively for the assumed asymmetrical modulation condition as reflected by polar diffusion, while in the A < 0 cycle, minute quantitative differences are obtained. However, when both polar and radial perpendicular diffusion have significant latitude dependences, major differences in radial intensities between the two polar angles are obtained in both polarity cycles. Furthermore, significant differences in radial intensity gradients obtained in the heliosheath at lower energies may suggest that the solar wind turbulence at and beyond the solar wind termination shock must have a larger latitudinal dependence.     

Cosmic ray anisotropies in the outer heliosphere

The observation of the directional distribution of energetic and cosmic ray particles has been done with the Voyager spacecraft over a long period. Since 2002, when the first flux enhancements of charged particles associated with the approach of Voyager 1 to the solar wind termination shock were observed, these anisotropy measurements have become of special interest. They play an important role to understand the magnetic field and shock structure and the basics of the modulation of cosmic ray and anomalous particles at and beyond the termination shock. They also serve as motivation to study the spatial behavior of galactic and anomalous cosmic ray anisotropies with numerical modulation models in order to illustrate how the radial anisotropy, at different energies, change from upstream to downstream of the termination shock. Observations made by Voyager 1 indicate that the termination shock is a complicated region than previously thought, hence the effects of the latitude dependence of the termination shock’s compression ratio and injection efficiency on the radial anisotropies of galactic and anomalous protons will be illustrated. We find that the magnitude and direction of the radial anisotropy strongly depends on the position in the heliosphere and the energy of particles. The effect of the TS on the radial anisotropy is to abruptly increase its value in the heliosheath especially in the A > 0 cycle for galactic protons and in both polarity cycles for anomalous protons. Furthermore, the global effect of the latitude dependence of the shock’s compression ratio is to increase the radial anisotropy for galactic protons throughout the heliosphere, while when combined with the latitude dependence of the injection efficiency this increase depends on modulation factors for anomalous protons and can even alter the direction of the radial anisotropy.     

Modelling of low-energy galactic electrons in the heliosheath

The modulation of cosmic ray electrons in the heliosphere plays an important role in improving our understanding and assessment of the processes applicable to low-energy galactic electrons. A full three-dimensional numerical model based on Parker’s transport equation is used to study the modulation of 10 MeV galactic electrons, in particular inside the heliosheath. The emphasis is placed on the role that perpendicular diffusion plays in causing the extraordinary large increase in the observed intensities of these electrons in the heliosheath. The modelling is compared with observations of 6–14 MeV electrons from the Voyager 1 mission. Results are shown for the radial intensity profiles of these electrons, as well as the modulation effects of varying the extent of the heliosheath by changing the location of the termination shock and the heliopause and the value of the local interstellar spectrum. We confirm that the heliosheath acts as a modulation ‘barrier’ for low-energy galactic electrons. The significance of this result depends on how wide the inner heliosheath is; on how high the very local interstellar spectrum is at these low energies (E < 100 MeV) and on how small perpendicular diffusion is inside the inner heliosheath.     

Cosmic ray modulation with a Fisk-type heliospheric magnetic field and a latitude-dependent solar wind speed

The effect of a latitude-dependent solar wind speed on a Fisk heliospheric magnetic field [Fisk, L. A. Motion of the footpoints of heliospheric magnetic field lines at the Sun: implications for recurrent energetic particle events at high heliographic latitudes. J. Geophys. Res. 101, 15547–15553, 1996] was first discussed by Schwadron and Schwadron and McComas [Schwadron, N.A. An explanation for strongly underwound magnetic field in co-rotating rarefaction regions and its relationship to footpoint motion on the the sun. Geophys. Res. Lett. 29, 1–8, 2002. and Schwadron, N.A., McComas, D.J. Heliospheric “FALTS”: favored acceleration locations at the termination shock. Geophys. Res. Lett. 30, 41–1, 2003]. Burger and Sello [Burger, R.A., Sello, P.C. The effect on cosmic ray modulation of a Parker field modified by a latitudinal-dependent solar wind speed. Adv. Space Res. 35, 643–646, 2005] found a significant effect for a simplified 2D version of a latitude-dependent Fisk-type field while Miyake and Yanagita [Miyake, S., Yanagita, S. The effect of a modified Parker field on the modulation of the galactic cosmic rays. In: Proceedings of 30th International Cosmic Ray Conference. Merida, Mexico, vol. 1, 445–448, 2007] found a smaller effect. The current report improves on a previous attempt Hitge and Burger [Hitge, M., Burger, R.A. The effect of a latitude-dependent solar wind speed on cosmic-ray modulation in a Fisk-type heliospheric magnetic field. In: Proceedings of 30th International Cosmic Ray Conference. Merida, Mexico, vol. 1, pp. 449–450, 2007] where the global change in the solar wind speed and not the local speed gradient was emphasized. The sheared Fisk field of Schwadron and McComas [Schwadron, N.A., McComas, D.J. Heliospheric “FALTS”: Favored acceleration locations at the termination shock. Geophys. Res. Lett. 30, 41–1, 2003.) is similar to the current Schwadron–Parker hybrid field. Little difference is found between the effects of a Parker field and a Schwadron–Parker hybrid field on cosmic-ray modulation, in contrast to the results of Burger and Sello and Miyake and Yanagita [Burger, R.A., Sello, P.C. The effect on cosmic ray modulation of a Parker field modified by a latitudinal-dependent solar wind speed. Adv. Space Res. 35, 643–646, 2005 and Miyake, S., Yanagita, S. The effect of a modified Parker field on the modulation of the galactic cosmic rays. In: Proceedings of 30th International Cosmic Ray Conference. Merida, Mexico, vol. 1, pp. 445–448, 2007]. The two-dimensional approximation used by these authors is therefore inadequate to model the complexities of the actual three-dimensional field. We also show that a Fisk-type field with a latitude-dependent solar wind speed (Schwadron–Parker hybrid field) decreases both the relative amplitude of recurrent cosmic ray intensity variations and latitude gradients and yields similar constants of proportionality for these quantities as for the constant solar wind speed case.     

Cosmic ray modulation in the outer heliosphere: Predictions for cosmic ray intensities up to the heliopause along Voyager 1 and 2 trajectories

Time dependent cosmic ray modulation in the outer heliosphere is calculated and results are compared to Voyager 1 and 2 observations using a two-dimensional time-dependent cosmic ray transport model. We predict possible future 133–242 MeV proton observations along the Voyager 1 and 2 spacecraft trajectories. Recent theoretical advances in cosmic ray transport parameters are introduced in order to provide a time-dependence for the assumed transport parameters used in the model. This leads to results that are in general compatible with the spacecraft observations in the inner and outer heliosphere over multiple solar cycles. However, for the outer heliosphere, we find that the Voyager 1 and 2 spacecraft observations cannot be fitted with an identical set of parameters along both trajectories. This indicates a possible asymmetric heliosphere or a symmetric heliosphere but with different diffusion parameters in the northern and southern hemispheres, respectively. Furthermore, results indicate that Voyager 2 observations are still under the influence of solar cycle related changes because of the large modulation volume between the heliopause and spacecraft location in contrast to Voyager 1 which shows a steady increase in cosmic ray intensities.     

The onset of sunspot cycle 24 and galactic cosmic ray modulation

The annual mean sunspot number (SSN) has a minimum value in 2008, while the monthly mean SSN has a value of zero in August 2009. The galactic cosmic ray modulation for cycle 24 began at earth orbit in January 2010. We study the onset characteristics of the new modulation cycle using data from the global network of neutron monitors. They respond to time variations in different segments of the galactic cosmic ray rigidity spectrum. The corresponding temporal variations in the interplanetary magnetic field intensity (B) and solar wind velocity (V) as well as the tilt angle of the heliospheric current sheet are also studied. There is a lag of 3 months between a large, sharp increase of the tilt angle of the heliospheric current sheet and the onset of modulation. Some neutron monitors are undergoing long-term drifts of unknown origin.     

The dynamic heliosphere,solar activity,and cosmic rays

This brief review addresses the relation between solar activity, cosmic ray variations and the dynamics of the heliosphere. The global features of the heliosphere influence what happens inside its boundaries on a variety of time-scales. Galactic and anomalous cosmic rays are the messengers that convey vital information on global heliospheric changes in the manner that they respond to these changes. By observing cosmic rays over a large range of energies at Earth, and with various space detectors, a better understanding is gained about space weather and climate. The causes of the cosmic ray variability are reviewed, with emphasis on the 11-year and 22-year cycles, step modulation, charge-sign dependent modulation and particle drifts. Advances in this field are selectively discussed in the context of what still are some of the important uncertainties and outstanding issues.     

North–South asymmetry in cosmic ray fluxes as measured in the stratosphere and in selected solar wind parameters in the near-Earth space

North–South asymmetry in the cosmic ray fluxes as resulted from the long-term balloon measurements in the northern and southern polar stratosphere does not agree with that found from the neutron monitor data. In order to reveal possible sources of the observed asymmetry, selected interplanetary parameters were examined. North–South asymmetry relative to the heliospheric neutral sheet was considered for solar wind velocity, plasma density and some other solar plasma parameters. It is shown that North–South asymmetry of the solar wind velocity and plasma density depends on the Earth’s heliolatitude and the phase of the 11-year solar activity cycle. This may be relevant to the results of cosmic ray measurements in the stratosphere.     

宇宙线福布什下降的数值模拟

本文从银河宇宙线的太阳调制方程出发,认为激波对宇宙线的影响是由一扰动区产生的。在此扰动区中太阳风速度增加,扩散系数下降。由此进行了数值模拟,模拟结果表明:扩散系数的下降在产生福布什下降时要比太阳风速度的增加更为有效;福布什下降在近日球层内向外传播时,其幅度随径向距离的增大而衰减;两个无相互作用的激波同时存在时所产生的福布什下降为每个激波单独存在时的福布什下降的简单线性迭加。      

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.     

The charge-sign dependent effect in the solar modulation of cosmic rays

The centennial anniversary of the discovery of cosmic rays was in 2012. Since this discovery considerable progress has been made on several aspects related to galactic cosmic rays in the heliosphere. It is known that they encounter a turbulent solar wind with an imbedded heliospheric magnetic field when entering the Sun’s domain. This leads to significant global and temporal changes in their intensity inside the heliosphere, a process known as the solar modulation of cosmic rays. The prediction of a charge-sign dependent effect in solar modulation in the late 1970s and the confirmatory observational discoveries can also be considered as a milestone. A short review is given of these predictions based on theoretical and numerical modelling work, the observational confirmation and related issues.     

MHD modeling of the outer heliosphere: Achievements and challenges

An overview is presented of magnetic-field-related effects in the solar wind (SW) interaction with the local interstellar medium (LISM) and the different theoretical approaches used in their investigation. We discuss the possibility that the interstellar magnetic field (ISMF) introduces north–south and east–west asymmetries of the heliosphere, which might explain observational data obtained by the Voyager 1 and Voyager 2 spacecraft. The SW–LISM interaction parameters that are responsible for the deflection of the interstellar neutral hydrogen flow from the direction of propagation of neutral helium in the inner heliosheath are outlined. The possibility of a strong ISMF, which increases the heliospheric asymmetry and the H–He flow deflection, is discussed. The effect of the combination of a slow-fast solar wind during solar minimum over the Sun’s 11-year activity cycle is illustrated. The consequences of a tilt between the Sun’s magnetic and rotational axes are analyzed. Band-like areas of an increased magnetic field distribution in the outer heliosheath are sought in order to discover regions of possible 2–3 kHz radio emission.     

On the transport of cosmic rays in the distant heliosheath

This paper discusses the transport of energetic charged particles through a sectored magnetic field in distant regions of the inner heliosheath. As the plasma flow slows down on approach to the stagnation point on the heliopause, the distance between the folds of the current sheet decreases to the point where it becomes comparable to the cyclotron radius of a cosmic ray particle. Under these conditions a particle can effectively drift across the stack of magnetic sectors with a speed comparable with the particle’s velocity. For a random distribution of current sheet separation distances, a diffusive transport across the stack of sectors occurs instead. The proposed mechanism could have contributed to unusually high intensities of galactic cosmic rays measured by Voyager 1 in the heliosheath during 2009–2010.     

Observed solar cycle modulation of galactic cosmic rays over a range of rigidities

We have studied the long-term, steady-state, solar cycle modulation of galactic cosmic ray intensity for seven cycles (17–23). Our analysis is based on the data obtained with a variety of detectors on earth (neutron monitors of the global network and muon detectors) as well as telescopes flown on high altitude balloons and on-board near-earth satellites. The median rigidity of response for these detectors to galactic cosmic ray spectrum lies in the range 1–70 GV. We correlate cosmic ray data to sunspot numbers, Ap, solar wind bulk speed (V), magnetic field (B), as well as to the cycle maximum (M), minimum (m), and the epochs of the solar polar field reversals. This enables us to derive the rigidity dependence of observations, and helps us to define the characteristics of the modulation function in the heliosphere.