Analytical approach to cosmic ray ionization by nuclei with charge Z in the middle atmosphere – Distribution of galactic CR effects

The effects of galactic and solar cosmic rays (CR) in the middle atmosphere are considered in this work. A new analytical approach for CR ionization by protons and nuclei with charge Z in the lower ionosphere and middle atmosphere is developed in this paper. For this purpose the ionization losses (dE/dh) according to the Bohr–Bethe–Bloch formula for the energetic charged particles are approximated in three different energy intervals. More accurate expressions for energy decrease E(h) and electron production rate profiles q(h) are derived. The obtained formulas allow comparatively easy computer programming. The integrand in q(h) gives the possibility for application of adequate numerical methods – such as Romberg method or Gauss quadrature, for the solution of the mathematical problem. On this way the process of interaction of cosmic ray particles with the upper, middle and lower atmosphere will be described much more realistically. Computations for cosmic ray ionization in the middle atmosphere are made. The full CR composition is taken into account: protons, Helium (α-particles), light L, medium M, heavy H and very heavy VH group of nuclei.     

Observation of enhanced sub-iron (Sc---Cr) to iron ration in low energy cosmic rays of 50–100 MeV/N in spacelab-3

Relative abundances of sub-iron (Sc-Cr) to iron nuclei in low energy (50–100 MeV/N) galactic cosmic rays have been determined from an analysis of about 100 events of heavy ions (Z = 10−28) recorded in a detector assembly flown in the Anuradha cosmic ray experiment in the Spacelab-3 on a six day mission in April–May 1985. The measured abundance ratio of (Sc-Cr)/Fe nuclei in 50–100 MeV/N energy range is 1.1 ± 0.3, and the present result of enhanced ratio of sub-iron to iron nuclei is in agreement with other experimental results in 200–800 MeV/N range. The over-abundance of iron secondaries at these low energies cannot be explained in the conventional models for propagation of cosmic rays. Available experimental data indicate a very different time history for the low energy iron-group, as compared to those of lighter nuclei in galactic cosmic rays.     

Propagation and origin of ultra high-energy cosmic rays

Propagation of UHE protons through CMB radiation leaves the imprint on energy spectrum in the form of Greisen–Zatsepin–Kuzmin (GZK) cutoff, bump (pile-up protons) and dip. The dip is a feature in energy range 1 × 10¹⁸–4 × 10¹⁹ eV, caused by electron–positron pair production on CMB photons. Calculated for power-law generation spectrum with index γ_g = 2.7, the shape of the dip is confirmed with high accuracy by data of Akeno-AGASA, HiRes, Yakutsk and Fly’s Eye detectors. The predicted shape of the dip is robust: it is valid for the rectilinear and diffusive propagation, for different discreteness in the source distribution, for local source overdensity, deficit, etc. This property of the dip allows us to use it for energy calibration of the detectors. The energy shift λ for each detector is determined by minimum χ² in comparison of observed and calculated dip. After this energy calibration the absolute fluxes, measured by AGASA, HiRes and Yakutsk detectors remarkably coincide in energy region 1 × 10¹⁸–1 × 10²⁰ eV. Below the characteristic energy E_c ≈ 1 × 10¹⁸ eV the spectrum of the dip flattens for both diffusive and rectilinear propagation, and more steep galactic spectrum becomes dominant at E < E_c. The energy of transition E_tr < E_c approximately coincides with the position of the second knee E_2kn, observed in the cosmic ray spectrum. The dip-induced transition from galactic to extragalactic cosmic rays at the second knee is compared with traditional model of transition at ankle, the feature observed at energy 1 × 10¹⁹ eV.     

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.     

The composition of cosmic rays at high energies

Measurements of the composition of the cosmic rays at high energies, and of the energy spectra of the individual components provide the basis for the understanding of the sources, of the acceleration mechanism, and of the galactic containment of these particles. We briefly review the presently available information, and we describe a recent measurement on the Space Shuttle that we performed in order to substantially extend the range of energies in which the elemental composition is known. We present and discuss the results, and we also summarize and discuss recent data on the electron component of cosmic rays. The body of data now available contains several features that are difficult to explain within current models of galactic shock acceleration and “leaky box” containment. We emphasize the need for further measurements, and we briefly discuss possible opportunities for future work.     

IONOSAT – Ionospheric satellite cluster

The IONOSAT project (from IONOspheric SATellites) is proposed by National Space Agency of Ukraine for First European Space Program as a part of Space Weather (SW) Program. As it is commonly accepted, Space Weather means the changes of the conditions on the Sun, in solar wind, magnetosphere and ionosphere which may affect the operation and reliability of on-board and ground technological systems and threaten human health. In this chain ionosphere is specific and integral part of SW formation. Moreover, namely in the ionosphere main part of the energy absorption of Sun-activated sporadic corpuscular and radiation fluxes takes places. The excitation of ionosphere by falling fluxes produces its “luminescence” in wide frequency band – from ULF waves till ultraviolet – and by this ionosphere works as an efficient “screen” or SW indicator.A goal of the proposed project is long-term spatial–temporal monitoring of main field and plasma parameters of ionosphere with aim to further develop fundamental conceptions of solar-terrestrial connections physics, nowcasting and forecast of SW, and diagnostics of natural and technogenic hazards with the help of scientific payload installed on-board a cluster of 3 low-Earth orbit (LEO) microsatellites (tentative launch date – 2012 year).The state of the project proposal and realization plans are discussed.     

Peculiarities of galactic cosmic ray Forbush effects during October–November 2003

Features of two successive Forbush effects of the galactic cosmic ray intensity in October–November 2003 have been studied based on the neutron monitors data. The rigidity spectrum of the galactic cosmic ray intensity in the course of the first Forbush effect (22–27 October) is gradually hardening, while the rigidity spectrum of the second Forbush effect (28 October–10 November) from the starting moment is very hard. As far, the energy range of the turbulence of the interplanetary magnetic field is in general responsible for the diffusion of galactic cosmic ray particles of the energy 5–50 GeV (to which neutron monitors are sensitive), we postulate that the gradually hardening (from day to day) of the rigidity spectrum of the first Forbush effect is associated with the enhancement of the power spectral density in the energy range of the interplanetary magnetic field turbulence caused by the large scale irregularities generated due to the interaction of the extending high speed disturbances with the background solar wind. The very hard rigidity spectrum (from the starting moment) of the second Forbush effect is generally associated with the well established new structure of the energy range of the interplanetary magnetic field turbulence enriched by the already created large scale irregularities. The gradually softening of the rigidity spectrum during the recovery phase of the second Forbush effect confirms that the disturbed interplanetary magnetic field turbulence step by step returns to the initial state.     

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.     

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.     

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.     

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 structure of the solar cycle maximum phase in the galactic cosmic ray 11-year variation

Two phenomena connected with the maximum phase of the 11-year solar cycle in the galactic cosmic ray intensity – the change in the energy dependence of the intensity variations and the double-peak structure in the intensity modulation time profile – are considered for the last five solar cycles (Nos. 19–23). The distinct 22-year cycle in the magnitude of the so called energy hysteresis is observed.The periods of the solar cycle maximum phase in the galactic cosmic ray intensity, characterized by the specific energy dependence of the intensity, are estimated. It is found that the double-peak structures belonging to the solar cycle maximum phase and those around it are very similar both in the amplitude and in its energy dependence.     

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.     

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.     

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.     

A review of recent results on cosmic gamma-ray bursts

The results obtained on cosmic gamma-ray bursts over the last several years are reviewed and compared with the older “historical” results. Fine time resolution measurements of burster light curves continue to reveal structure at the millisecond and sub-millisecond level, suggesting a compact object origin. Similarly, the evolution of the low energy X-ray spectra of bursts towards shapes consistent with 1–2 keV blackbodies may be interpreted in terms of a neutron star origin, as can the continuing detection of absorption and emission features. The statistical evidence, however, argues strongly for an isotropic distribution which has been completely sampled. To reconcile this with galactic neutron stars requires the assumption that they are Population II objects. Counterpart searches have evolved to the point where they may be carried out within days of an event, and a soft X-ray source has now been detected in the error box of one recent burst.     

On ultra-high energy cosmic rays: Origin in AGN jets and transport in expanding universe

The cosmic ray source spectrum produced by AGN (active galactic nucleus) jets is calculated. A distinctive feature of these calculations is the account for the jet distribution on kinetic energy. The expected cosmic ray spectrum at the Earth is determined with the use of a simple numerical code which takes into account interactions of ultra-high energy protons and nuclei with the background radiation in an expanding universe.     

Rapid variability of 10–140 keV X-rays from Cygnus X-1

On five occasions in 1977 and 1978, Cygnus X-1 was observed using the Low-Energy Detectors of the UCSD/MIT Hard X-Ray and Low-Energy Gamma-Ray Experiment on the HEAO-1 satellite. Rapid (0.08 s ≤ t ≤ 1000 s) variability was found in the 10 – 140 keV band. The power spectrum was “white” for 10⁻³ Hz < f ≤ 5 × 10⁻² Hz and was proportional to f⁻¹ for 5 × 10⁻² Hz ≤ f < 3 Hz, indicating correlations on all time scales < 20 s. If the emission is produced by Comptonization of a soft photon flux in a hot cloud, the heating of the cloud cannot be constant; it must vary on time scales up to 20 seconds. A variable accretion rate could cause the observed effects.     

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.     

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.