The IGY and beyond: A brief history of ground-based cosmic-ray detectors

The state of art of ground-based cosmic-ray research from its discovery to present is reviewed. After discovery of cosmic rays by Hess in 1912, the nature of the primary and secondary radiation was established from recordings by a variety of instruments, sensitive to various components of cosmic rays and operated at different latitudes, longitudes and altitudes, including instruments carried by balloons. The IGY formalized international co-operation and coordinated study of cosmic rays, which is vital for meaningful interpretation of cosmic-ray data. Data collected at different geographic locations require an effective cutoff rigidity as a data ordering parameter. This parameter is obtained from tracing trajectories of primary cosmic rays in the Earth’s magnetic field. After 50 years the world’s neutron monitor network remains still the backbone for studying intensity variations of primary cosmic rays in the rigidity ranges between 1 and 15 GV, associated with transport and with transient events. Also the penetrating muon and neutrino components of secondary cosmic rays have a long history of recording and fundamental problem investigations. Valuable data about composition and spectrum of primary cosmic rays in ever increasing high-energy regions have been obtained during the years of investigations with various configurations and types of extensive air shower detectors. The culture of personal involvement of the physicist in carrying out experiments and data acquisition characterized the continued vitality of cosmic-ray investigations ranging from its atmospheric, geomagnetic and heliospheric transport through to its solar and astrophysical origins.     

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.     

Non-linear diffusion of cosmic rays

The propagation of cosmic rays in the interstellar medium after their release from the sources – supernova remnants – can be attended by the development of streaming instability. The instability creates MHD turbulence that changes the conditions of particle transport and leads to a non-linear diffusion of cosmic rays. We present a self-similar solution of the equation of non-linear diffusion for particles ejected from a SNR and discuss how obtained results may change the physical picture of cosmic ray propagation in the Galaxy.     

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.     

Some additional results in the regular balloon monitoring of cosmic rays in the Earth’s atmosphere

The different types of the data recorded in the experiment of the regular balloon monitoring of cosmic rays (carried out since 1957 by Lebedev Physical Institute, Moscow, Russia, in several locations) are described. So called detailed information (the form of each pulse detected by the ground-based receiver) recorded during the last 12 years is discussed in more details. The use of these data both for getting and correcting the standard results of the experiment and for obtaining some additional information on the cosmic rays in the Earth’s atmosphere is considered.     

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.     

Research by the Tasmanian cosmic ray group during the International Geophysical Year

Systematic recording of the cosmic radiation commenced in Hobart in 1946 and at Mawson in Antarctica in 1955, making these two of the longest running cosmic ray observatories in the world. For the IGY, observations were also made at a sub-Antarctic island and near the equator, and an airborne survey of the nucleonic component was made from Geomagnetic Latitude −60°, south of Australia, to Japan and back. At Hobart there were neutron monitors, vertical and inclined muon telescopes, an ionization chamber, and two muon telescopes at ∼40 m of water equivalent underground. The research based on these and other observations determined the energy dependence of the Forbush and 11-year variations and concentrated, in particular, on understanding the anisotropic nature of galactic cosmic rays up to 150 GeV; the anisotropies in the onset phase of Forbush decreases; and the anisotropies in solar cosmic ray events. An investigation was initiated to calculate the trajectories and cutoff rigidities of cosmic rays in a high order simulation of the geomagnetic field. This was completed in 1959–60.     

Study of the 27-day variations of the galactic cosmic ray intensity and anisotropy

We demonstrate that the general features of the radial and azimuthal components of the anisotropy of galactic cosmic rays can be studied by the harmonic analysis method using data from an individual neutron monitor with cut off rigidity <5 GV. In particular, we study the characteristics of the 27-day (solar rotation period) variations of the galactic cosmic ray intensity and anisotropy, solar wind velocity, interplanetary magnetic field strength and sunspot number. The amplitudes of the 27-day variations of the galactic cosmic ray anisotropy are greater, and the phases more clearly established, in A > 0 polarity periods than in A < 0 polarity periods at times of minimum solar activity. The phases of the 27-day variations of the galactic cosmic rays intensity and anisotropy are opposite with respect to the similar changes of the solar wind velocity in A > 0 polarity periods. No significant dependence of the amplitude of the 27-day variation of the galactic cosmic ray anisotropy on the tilt angle of the heliospheric neutral sheet is found. Daily epicyclegrams obtained by Chree’s method show that the 27-day variations of the galactic cosmic ray anisotropy during A > 0 polarity periods follow elliptical paths with the major axes oriented approximately along the interplanetary magnetic field. The paths are more irregular during A < 0 polarity periods.     

Radial distribution of galactic cosmic rays at solar maximum

In this paper we analyze the spatial distribution of galactic cosmic rays during periods of maximum solar activity of the cycles 21, 22 and 23. We have used a two dimensional model to solve the cosmic ray transport equation. This model includes all relevant physical processes: diffusion, convection, drift and shock effects on cosmic ray propagation inside the heliosphere. We focused on the study of the radial distribution of galactic cosmic rays, and compare our results with the spacecraft observations for two energies (175 MeV H and 265 MeV/n He). Although the radial intensities of galactic cosmic rays can be explained qualitatively with all three local interstellar spectra (LISs) used in this work, we applied a reduced chi-squared analysis to investigate the best LIS that could fit the data.     

Fluctuations of cosmic rays and IMF in the vicinity of interplanetary shocks

Fluctuations of cosmic rays and interplanetary magnetic field upstream of interplanetary shocks are studied using data of ground-based polar neutron monitors as well as measurements of energetic particles and solar wind plasma parameters aboard the ACE spacecraft. It is shown that coherent cosmic ray fluctuations in the energy range from 10 keV to 1 GeV are often observed at the Earth’s orbit before the arrival of interplanetary shocks. This corresponds to an increase of solar wind turbulence level by more than the order of magnitude upstream of the shock. We suggest a scenario where the cosmic ray fluctuation spectrum is modulated by fast magnetosonic waves generated by flux of low-energy cosmic rays which are reflected and/or accelerated by an interplanetary shock.     

On possible genesis of fractal dimensions in the turbulent pulsations of cosmic plasma – galactic cosmic rays – turbulent pulsation in planetary atmosphere system

The spectrum of turbulent pulsations induced in the atmosphere by the galactic cosmic rays is defined. A possible manifestation of genesis of fractal dimensions in the system of “spectrum of turbulent pulsations of cosmic plasma – galactic cosmic rays’ spectrum – spectrum of atmospheric turbulent pulsations” is analyzed.     

Ion production and ionization effect in the atmosphere during the Bastille day GLE 59 due to high energy SEPs

The influence of high energy particles, specifically cosmic rays, on atmospheric physics and chemistry is highly discussed. In most of the proposed models the role of ionization in the atmosphere due to cosmic rays is not negligible. Moreover, effect(s) on minor constituents and aerosols are recently observed, specifically over the polar regions during strong solar particle events. According to the recent findings for such effects it is necessary an essential increase of ion production, specifically during the winter period. The galactic cosmic rays are the main source of ionization in the Earth’s stratosphere and troposphere. Occasionally, the atmospheric ionization is significantly enhanced during strong solar energetic particles events, specifically over the polar caps. During the solar cycle 23 several strong ground level enhancements were observed. One of the strongest was the Bastille day event occurred on 14 July 2000. Using a full Monte Carlo 3-D model, we compute the atmospheric ionization, considering explicitly the contribution of cosmic rays with galactic and solar origin, focusing on high energy particles. The model is based on atmospheric cascade simulation with the PLANETOCOSMICS code. The ion production rate is computed as a function of the altitude above the sea level. The ion production rate is computed on a step ranging from 10 to 30?min throughout the event, considering explicitly the spectral and angular characteristics of the high energy part of solar protons as well as their time evolution. The corresponding event averaged ionization effect relative to the average due to galactic cosmic rays is computed in lower stratosphere and upper troposphere at various altitudes, namely 20?km, 15?km, 12?km and 8?km above the sea level in a sub-polar and polar regions. The 24^h and the weekly ionization effects are also computed in the troposphere and low stratosphere. Several applications are discussed.     

Empirical model of cosmic ray spectrum in energy interval 1 MeV–100 GeV during 11-year solar cycle

The galactic cosmic rays (GCR) are the main ionization source at altitude of ∼3–35 km in the atmosphere. For high latitude anomalous cosmic ray (ACR) component has also a significant influence on the atmospheric ionization. We propose an empirical model for differential spectra D(E) of galactic and anomalous cosmic rays in energy interval 1 MeV–100 GeV during solar cycle. In the model data are used which cover three solar cycles: 20, 22 and 23. The LEAP87, IMAX92, CAPRICE94, AMS98 and BESS experimental spectra for protons and alpha particles are fitted to the proposed empirical model. The modulated GCR differential spectra are compared with force-field approximation to the one-dimensional transport equation and with solutions of two-dimensional cosmic ray transport equation. For experimental spectra, the calculation of the model parameters is performed by Levenberg–Marquardt algorithm, applied to the special case of least squares. Algorithm that combines the rapid local convergence of Newton–Raphson method with globally convergent method for non-linear systems of equations is applied for theoretically obtained differential spectra. The described programmes are realized in algorithmic language C++. The proposed model gives practical possibility for investigation of experimental data from measurements of galactic cosmic rays and their anomalous component.     

The physics of cosmic-ray modulation

The theory of the modulation of galactic cosmic rays by the solar wind is reviewed. The basic transport equation is presented, interpreted and then applied to cosmic-ray transport in a model heliosphere immersed in a constant uniform bath of galactic cosmic rays. The results of numerical modelling are presented and the dominant physical effects analyzed. A variety of observational tests of the model which were reported over the last several years are summarized and shown, generally, to support a model in which particle drifts play an important role. Recent measurements which show that the latitudinal gradient of cosmic rays changed sign in the recent sunspot minimum (relative the last sunspot minimum) are shown to provide additional, strong, support for the model. A new picture of the interplanetary magnetic field is presented, which gives promise of improving considerably the agreement between the theory and observations in the few remaining problem areas.     

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.     

Effects of solar modulation on the cosmic ray positron fraction

We implemented a 2D Monte Carlo model to simulate the solar modulation of galactic cosmic rays. The model is based on the Parker’s transport equation which contains diffusion, convection, particle drift and energy loss. Following the evolution in time of the solar activity, we are able to modulate a local interstellar spectrum (LIS), that we assumed isotropic beyond the termination shock, down to the Earth position inside the heliosphere. In this work we focused our attention to the cosmic ray positron fraction at energy below ∼10 GeV, showing how the particle drift processes could explain different results for AMS-01 and PAMELA. We compare our modulated spectra with observations at Earth, and then make a prediction of the cosmic ray positron fraction for the AMS-02 experiment.     

Measurement of cosmic ray elemental composition from the CAKE balloon experiment

CAKE (Cosmic Abundances below Knee Energies) was a prototype balloon experiment for the determination of the charge spectra and abundances of the primary cosmic rays (CR) with Z > 10. It was a passive instrument made of layers of CR39® and Lexan®/Makrofol® nuclear track detectors; it had a geometric acceptance of ∼0.7 m² sr for Fe nuclei. Here, the scanning and analysis strategies, the algorithms used for the off-line filtering and for the tracking in automated mode of the primary cosmic rays are presented, together with the resulting CR charge distribution and their abundances.     

Cosmic-ray modulation and the heliospheric magnetic field

The heliospheric magnetic field plays a key role in any model for the modulation of cosmic rays. It enters into all diffusion coefficients, and its magnitude, spatial gradient and direction determine drifts patterns of cosmic rays in the heliosphere. While the first axisymmetric model of E.N. Parker proved quite successful to explain in situ measurements in the ecliptic plane, new insight into the origin and the nature of the field, especially at high heliographic latitudes, has led to the development of complex fully-three-dimensional, time-dependent models. In this review, we discuss a selection of models for the heliospheric magnetic field, and discuss how some of the more recent Fisk-type models affect the modulation of cosmic rays.     

The scaler mode in the Pierre Auger Observatory to study heliospheric modulation of cosmic rays

The impact of the solar activity on the heliosphere has a strong influence on the modulation of the flux of low energy galactic cosmic rays arriving at Earth. Different instruments, such as neutron monitors or muon detectors, have been recording the variability of the cosmic ray flux at ground level for several decades. Although the Pierre Auger Observatory was designed to observe cosmic rays at the highest energies, it also records the count rates of low energy secondary particles (the scaler mode) for the self-calibration of its surface detector array. From observations using the scaler mode at the Pierre Auger Observatory, modulation of galactic cosmic rays due to solar transient activity has been observed (e.g., Forbush decreases). Due to the high total count rate coming from the combined area of its detectors, the Pierre Auger Observatory (its detectors have a total area greater than 16,000 m²) detects a flux of secondary particles of the order of ∼10⁸ counts per minute. Time variations of the cosmic ray flux related to the activity of the heliosphere can be determined with high accuracy. In this paper we briefly describe the scaler mode and analyze a Forbush decrease together with the interplanetary coronal mass ejection that originated it. The Auger scaler data are now publicly available.