Model for induced ionization by galactic cosmic rays in the Earth atmosphere and ionosphere

The ionization profiles produced by galactic cosmic rays in the Earth atmosphere and ionosphere are obtained on the basis of Monte Carlo simulations. Cascade processes in the atmosphere are simulated using CORSIKA 6.52 code with FLUKA 2006 and QGSJET II hadronic interaction subroutines. Proton induced showers are considered using a realistic atmospheric model (US Standard Atmosphere). The energy deposit from different components is taken into account, namely electromagnetic, hadron and muon components. The curvature of the atmosphere is considered in the computer code. On the basis of the computed ionization yield function the ion pair production rate in the atmosphere is obtained for different conditions and locations. The model is applicable to the entire atmosphere, from ground level to upper atmosphere. Several applications of the obtained results are discussed. The Monte Carlo simulation model considers nuclear interactions below the altitude of 35 km. It is compared with analytical–numerical electron production rate model. The latter model which takes into account the electromagnetic interactions above altitudes of 35 km has two main regions of application: above 50 km (thin target model) and between 35 and 50 km (intermediate target model). A good agreement between the CORSIKA results and analytical–numerical model results is found above altitude of 35 km.     

Time-dependent cosmic ray modulation

Time-dependent cosmic ray modulation is calculated over multiple solar cycles using our well established two-dimensional time-dependent modulation model. Results are compared to Voyager 1, Ulysses and IMP cosmic ray observations to establish compatibility. A time-dependence in the diffusion and drift coefficients, implicitly contained in recent expressions derived by , ,  and , is incorporated into the cosmic ray modulation model. This results in calculations which are compatible with spacecraft observations on a global scale over consecutive solar cycles. This approach compares well to the successful compound approach of Ferreira and Potgieter (2004). For both these approaches the magnetic field magnitude, variance of the field and current sheet tilt angle values observed at Earth are transported time-dependently into the outer heliosphere. However, when results are compared to observations for extreme solar maximum, the computed step-like modulation is not as pronounced as observed. This indicates that some additional merging of these structures into more pronounced modulation barriers along the way is needed.     

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.     

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.     

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.     

Solar, geomagnetic and cosmic ray intensity changes, preceding the cyclone appearances around Mexico

Recently it has been suggested that there exist specific changes in the cosmic ray intensity and some solar and geomagnetic parameters during the days, preceding the hurricane appearances over the North Atlantic Ocean. To understand better these phenomena, data for all hurricanes born not only over the Atlantic but also over the Pacific waters in the last 55 years that hit the Mexican borders were elaborated. As basic hurricane parameters the maximum rotational velocity and the estimated total energy were used. To avoid any interference all hurricanes, overlapping the preceding ones with more than 20 days were not included. Then the behavior of the cosmic ray (CR) intensity, the sunspot (SS) numbers, and the geomagnetic parameters (AP) and (KP) in 35 days prior and 20 days after the cyclone start were investigated. The CR, SS, AP and KP showed much more intensive disturbances in the periods preceding and following the hurricane appearance. For SS this disturbance gradually increase with the hurricane strength. A characteristic peak in the CR intensity appears before the hurricane start. But its place varies between 5 and 20 days before that start. Specific changes were observed in the SS. For major hurricanes they begins sometimes more than 20 days in advance. The AP and the KP show series of bursts, spread over the whole period of 30 preceding days. The obtained results from the performed correlational analysis are enough interesting to motivate a further statistical analysis with more precise techniques: in particular a common periodicity of 30 years found in the number of tropical storms landing into Mexico, the averaged rotational wind velocity and the ACE must be studied in connection with the solar Hale cycle. Using coherence wavelet spectral analysis we present a comparative study between one terrestrial and one cosmophysical phenomena that presumable influence hurricanes development: African dust outbreaks versus cosmic rays for all North Atlantic tropical cyclones. It is shown that the cosmophysical influence cannot be considered as a negligible effect.     

Modeling ground and space based cosmic ray observations

After entering our local astrosphere (called the heliosphere), galactic cosmic rays, as charged particles, are affected by the Sun’s turbulent magnetic field. This causes their intensities to decrease towards the inner heliosphere, a process referred to as modulation. Over the years, cosmic ray modulation has been studied extensively at Earth, utilizing both ground and space based observations. Moreover, modelling cosmic ray modulation and comparing results with observations, insight can be gained into the transport of these particles, as well as offering explanations for observed features. We review some of the most prominent cosmic ray observations made near Earth, how these observations can be modelled and what main insights are gained from this modelling approach. Furthermore, a discussion on drifts, as one of the main modulation processes, are given as well as how drift effects manifest in near Earth observations. We conclude by discussing the contemporary challenges, fuelled by observations, which are presently being investigated. A main challenge is explaining observations made during the past unusual solar minimum.     

Relation between different theories for cosmic ray cross field diffusion

In the past few year several theories have been presented for describing cosmic ray scattering across the mean magnetic field. It is the purpose of the present article to discuss the relation between these different theories in order to improve our understanding of cosmic ray perpendicular scattering and to explore the parameter regimes for which these different theories are valid.     

Galactic cosmic ray modulation for sunspot cycle 23

We present a study of the galactic cosmic ray modulation for sunspot cycle 23. We use the monthly and the annual mean hourly, pressure corrected, data from neutron monitors of the global network (monthly rate is calculated as the average of the hourly pressure corrected values). We draw attention to an asymmetry in the galactic cosmic ray (GCR) recovery during odd and even cycles for the monthly mean hourly rate data. For over half a century of observations, we find that the recovery for the odd cycles is to a higher level than for the even cycles. Qualitatively the effect is ascribed to charged particle drifts in inhomogeneous interplanetary magnetic field. Even so it has not been possible to arrive at a quantitative, self-consistent, explanation in terms of drifts at higher and lower GCR rigidities. We also study the rigidity dependence of the amplitude of 11-year modulation over a wide range (1–200 GV) of GCR spectrum; it is a power law in rigidity with an exponent −1.22. We discuss the implication of these findings on quasi-linear diffusion theories of modulation. We reflect on GCR recovery pattern for 2006–2009.     

The brief history of experimental research of cosmic ray variations in Yakutia

Some unknown historical facts of cosmic ray studies in the north-east of the former Soviet Union related to the Yakutsk scientific group are reported for the benefit of the international scientific community. It focuses on the founders of Yu.G. Shafer Institute of Cosmophysical Research and Aeronomy of Siberian Branch of Russian Academy of Sciences. A chronology of measurements of cosmic ray intensity variations since 1949 in Yakutia (Sakha Republic; NE Siberia) is given. In particular, for the first time the data of the first solar cosmic ray event registered at Yakutsk (GLE04), with a small ionization chamber S-2 (volume: 20 L) are presented. Moreover, the data of the large ionization chamber ASK-1 (volume: 950 L) for the 1953–2003 period useful for specialists in the field of cosmic ray variations are also shown.     

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.     

Rigidity spectrum of the 27-day variation of the galactic cosmic ray intensity in different epochs of solar activity

We study the temporal evolution of the power rigidity spectrum of the first (27 days) and the second (14 days) harmonics of the 27-day variation of the galactic cosmic ray intensity measured by neutron monitors in the period of 1965–2002. The rigidity spectrum of these variations can be approximated by a power law. We show the rigidity spectra of the first and the second harmonics of the 27-day variation of the galactic cosmic ray intensity have similar time profiles. These spectra are hard (γ ≈ 0.5 ± 0.1) and soft (γ ≈ 1.1 ± 0.2) during solar maximum and minimum activity, respectively. We ascribe this to the alternation of the sizes of the modulation regions responsible for the 27-day variation of the galactic cosmic ray intensity in different epochs of solar activity. Especially, the average radial sizes of the modulation regions of the 27-day variation and the heliolatitudinal extension of the heliolongitudinal asymmetry are smaller during solar minimum than during solar maximum. We show also, that the temporal changes of the power rigidity spectra of the first and the second harmonics of the 27-day variation of the galactic cosmic ray intensity are in a negative correlation with the changes of the rigidity spectrum of the corresponding 11-year variation.     

Early history of cosmic rays and solar wind – Some personal remembrances

The early history of solar wind was replete with prejudices and strong opposition to Parker’s formulation. It was only after conclusive evidence from satellite data was obtained that the idea of solar wind was accepted. Some personal experiences of mine during my stay at the University of Chicago in 1953–1954, including the encounter of Dr. Simpson with Dr. Biermann and the inconclusive discussion between them about a possible perpetual solar outflow of particles are presented and further developments when Parker came to Chicago in 1956 and formulated his idea of solar wind, as narrated to me later by Dr. Simpson, are described.     

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.     

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.     

Evolution of magnetic fields and cosmic ray acceleration in supernova remnants

Observations show that the magnetic field in young supernova remnants (SNRs) is significantly stronger than can be expected from the compression of the circumstellar medium (CSM) by a factor of four expected for strong blast waves. Additionally, the polarization is mainly radial, which is also contrary to expectation from compression of the CSM magnetic field. Cosmic rays (CRs) may help to explain these two observed features. They can increase the compression ratio to factors well over those of regular strong shocks by adding a relativistic plasma component to the pressure, and by draining the shock of energy when CRs escape from the region. The higher compression ratio will also allow for the contact discontinuity, which is subject to the Rayleigh–Taylor (R–T) instability, to reach much further out to the forward shock. This could create a preferred radial polarization of the magnetic field. With an Adaptive Mesh Refinement MHD code (AMRVAC), we simulate the evolution of SNRs with three different configurations of the initial CSM magnetic field, and look at two different equations of state in order to look at the possible influence of a CR plasma component. The spectrum of CRs can be simulated using test particles, of which we also show some preliminary results that agree well with available analytical solutions.     

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.     

On the effects of dynamical turbulence on the perpendicular diffusion of low-energy cosmic ray electrons

Diffusion perpendicular to the heliospheric magnetic field plays an integral role in the transport of charged particles in the heliosphere. In this study the perpendicular diffusion coefficient of low-energy cosmic ray electrons is calculated, using an equation derived from the random ballistic decorrelation interpretation of nonlinear guiding centre theory. An observationally motivated 2D turbulence power spectrum is assumed and the effects of various turbulence inputs on the resulting perpendicular diffusion coefficient are investigated. The perpendicular diffusion coefficients are first determined at 1 AU, for both magnetostatic and dynamical turbulence conditions. These solutions are also evaluated for radial distances of 0.1 AU to 10 AU to further investigate the values of the perpendicular diffusion coefficients in the very inner heliosphere. The results of this study show that the dissipation range of the turbulence power spectrum provides a negligible contribution towards the perpendicular diffusion coefficient, and that solutions derived using only the energy containing range serve as good approximations for solutions derived assuming the full 2D turbulence power spectrum. Finally, it is shown that the effects of dynamical turbulence, as considered in the present study, do not affect the perpendicular diffusion coefficients derived from the scattering theory considered here.     

First measurements of periodicities and anisotropies of cosmic ray flux observed with a water-Cherenkov detector at the Marambio Antarctic base

A new water-Cherenkov radiation detector, located at the Argentine Marambio Antarctic Base (64.24S-56.62 W), has been monitoring the variability of galactic cosmic ray (GCR) flux since 2019. One of the main aims is to provide experimental data necessary to study interplanetary transport of GCRs during transient events at different space/time scales. In this paper we present the detector and analyze observations made during one full year. After the analysis and correction of the GCR flux variability due to the atmospheric conditions (pressure and temperature), a study of the periodicities is performed in order to analyze modulations due to heliospheric phenomena. We can observe two periods: (a) 1 day, associated with the Earth’s rotation combined with the spatial anisotropy of the GCR flux; and (b) $～$ 30 days due to solar impact of stable solar structures combined with the rotation of the Sun. From a superposed epoch analysis, and considering the geomagnetic effects, the mean diurnal amplitude is $～$ 0.08% and the maximum flux is observed in $～$ 15 h local time (LT) direction in the interplanetary space. In such a way, we determine the capability of Neurus to observe anisotropies and other interplanetary modulations on the GCR flux arriving at the Earth.