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.     

27-day variations of the galactic cosmic ray intensity and anisotropy in the minimum of the 23rd solar activity cycle

We study the 27-day variations of the solar wind velocity, galactic cosmic ray (GCR) intensity and anisotropy in the last minimum epoch of solar activity (2007–2009, A < 0). The average amplitude of the 27-day variation of the galactic cosmic ray anisotropy (A27A) in the current minimum epoch of solar activity (2007–2009, A < 0) is lesser than in previous positive polarity period as it is expected from the drift theory. So, polarity dependence rule for the 27-day variation of the GCR anisotropy is fully kept. It is a universal principle for the amplitudes of the 27-day variation of the GCR anisotropy. At the same time, the average amplitude of the 27-day variation of the GCR intensity (A27I) remains at the same level as for previous minimum epoch 1995–1997 (A > 0) showing by the same token an violation of its polarity dependence rule established earlier. We assume that this phenomenon could be generally related with the well established 27-day variation of the solar wind velocity being in anti-correlation with the similar changes of the 27-day variation of the GCR intensity. Generally, a character of the heliolongitudinal asymmetry of spatial large-scale structure of the solar wind velocity (SWV) established in the recent minimum epoch, preferentially pronounces in the behavior of the 27-day variation of the GCR intensity than anisotropy. The formation of the 27-day variation of the GCR anisotropy preferentially takes place in a restricted disk like local vicinity in the helioequatorial region, whilst the 27-day variation of the GCR intensity is formed in the global three dimensional vicinity of the heliosphere.     

The role of drift and convection–diffusion mechanisms in small energy global cosmic ray modulation: Application to the hysteresis phenomenon of proton and alpha particle satellite data

The hysteresis effect for small energies of galactic cosmic rays is due to two effects. The first is the same as for neutron monitor energies – the delay of the interplanetary processes responsible for cosmic ray modulation with respect to the initiating solar processes, according to the effective velocity of solar wind and shock waves propagation. Then, the observed cosmic ray intensity is connected to the solar activity variations during many months before the time of cosmic ray measurement. The second is caused by the time delay of small energy cosmic ray diffusion from the boundary of modulation region to the Earth’s orbit. The model describing the connection between solar activity variation and cosmic ray convection–diffusion global modulation for neutron monitor energies is here developed by taking into account also the time-lag of the small energy particle diffusion in the Heliosphere. We use theoretical results on drifts and analytically approximate the dependences of drifts from tilt angle, and take into account the dependence from the sign of primary particles, and from the sign of polar magnetic field (A > 0 or A < 0). The obtained results are applied on proton and alpha-particle satellite data. We analyze satellite 5-min data of proton fluxes with energies >1 MeV, >2 MeV, >5 MeV, >10 MeV, >30 MeV, >50 MeV, >60 MeV, >100 MeV, and in intervals 10–30 MeV, 30–60 MeV, and 60–100 MeV during January 1986–December 1999. We exclude periods with great cosmic ray increases caused by particle acceleration in solar flare events. Then, we determine monthly averaged fluxes, as well as 5-month and 11-month smoothed data. We analyze also satellite 5-min data on alpha-particle fluxes in the energy intervals 60-160 MeV, 160–260 MeV and 330–500 MeV during January 1986–May 2000. We correct observation data for drifts and then compare with what is expected according to the convection–diffusion mechanism. We assume different dimensions of the modulation region (by the time propagation X₀ of solar wind from the Sun to the boundary of modulation region), for X₀ values from 1 to 60 average months, by one-month steps. For each value of X₀ we determine the correlation coefficient between variations of expected and observed cosmic ray intensities (the estimation of cosmic ray intensities values is given in Section 3 by Eq. (9), and the determination of correlation and regression coefficients in Section 3 by Eq. (8)). The dimension of modulation region is determined by the value of X_0 max, for which the correlation coefficient reaches the maximum value. Then the effective radial diffusion coefficient and residual modulation in small energy region can be estimated.     

Dependence of the 27-day variation of cosmic rays on the global magnetic field of the Sun

We show that the higher range of the heliolongitudinal asymmetry of the solar wind speed in the positive polarity period (A > 0) than in the negative polarity period (A < 0) is one of the important reasons of the larger amplitudes of the 27-day variation of the galactic cosmic ray (GCR) intensity in the period of 1995–1997 (A > 0) than in 1985–1987 (A < 0). Subsequently, different ranges of the heliolongitudinal asymmetry of the solar wind speed jointly with equally important corresponding drift effect are general causes of the polarity dependence of the amplitudes of the 27-day variation of the GCR intensity. At the same time, we show that the polarity dependence is feeble for the last unusual minimum epoch of solar activity 2007–2009 (A < 0); the amplitude of the 27-day variation of the GCR intensity shows only a tendency of the polarity dependence. We present a three dimensional (3-D) model of the 27-day variation of GCR based on the Parker’s transport equation. In the 3-D model is implemented a longitudinal variation of the solar wind speed reproducing in situ measurements and corresponding divergence-free interplanetary magnetic field (IMF) derived from the Maxwell’s equations. We show that results of the proposed 3-D modeling of the 27-day variation of GCR intensity for different polarities of the solar magnetic cycle are in good agreement with the neutron monitors experimental data. To reach a compatibility of the theoretical modeling with observations for the last minimum epoch of solar activity 2007–2009 (A < 0) a parallel diffusion coefficient was increased by ∼40%.     

The spatial distribution of galactic and anomalous cosmic rays in the heliosphere at solar minimum

The spatial distributions of galactic and anomalous cosmic rays in the heliosphere at the solar minima of Cycles 20/22 (qA > 0) and of Cycle 21 (qA < 0) are studied, using data from IMP 8, Voyagers 1/2 and Pioneer 10. It is found that the radial dependences of intensities J can be approximated by a power of radial distance r as J ∝ r^α with a different value of a constant in the inner and outer heliosphere with a transition at a radial distance of 10–15 AU. To study the physical meaning of these radial intensity profiles we examined the rigidity dependences of the intensity gradients by determining the particle mean free paths, using a simple one-dimensional modulation model. The particle mean free path λ was assumed to be a separable function of distance of the form r^γ and rigidity R of R^δ over the range of 0.5–3.0 GV in the inner and outer heliosphere. It was shown that λ of rigidity dependence of R^1.6 determined for Cycle 20/22 (qA > 0) with anomalous He is about 4–5 times larger than that of Cycle 21 (qA < 0) with R^0.9 at around 1 GV in the outer heliosphere, and that the radial dependences are r^1.4 and r^1.1, respectively, for Cycles 20/22 and for Cycle 21.     

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.     

Features of the 11-year variation of galactic cosmic rays in different periods of solar magnetic cycles

Variations of galactic cosmic ray intensity have been studied based on the neutron monitors and interplanetary magnetic field experimental data for different ascending and descending epochs of solar activity. The dependence of the diffusion coefficient on the cosmic ray particles rigidity R is stronger in the maxima epoch than in the minima epoch of solar activity. For the period of 1977–1981 (qA > 0) the diffusion coefficient for the minimum epoch is, χ_min ∝ R^{0.7 ± 0.04} and for the maximum χ_max ∝ R^{1.3 ± 0.05}; for the period of 1987–1990 (qA < 0), χ_min ∝ R^{0.8 ± 0.05} and χ_max ∝ R^{1.1 ± 0.04}. The exponents ν_y and ν_z of the power spectral density of the B_y and B_z components of the IMF in the region of the frequencies (10⁻⁶– 4 × 10⁻⁶) Hz are larger for the minimum epoch of 1987 (ν_y ≈ 2.0 and ν_z ≈ 1.93) than for the maximum epoch of 1990 (ν_y ≈ 1.43 and ν_z ≈ 1.27).     

α–ω Effect and peculiarities of the 27-day variations of the galactic cosmic rays intensity,solar wind and solar activity parameters

We show that the amplitudes of the 27-day variations of galactic cosmic ray (GCR) intensity, solar wind and solar activity parameters have a periodicity with duration of three to four Carrington rotation periods (3–4 CRP). We assume that the general reason for this phenomenon may be related to similar cyclicity of topological structure of the solar magnetic field lines created owing to the asymmetry of turbulent solar dynamo and solar differential rotation transforming the Sun’s poloidal magnetic field to the toroidal (α–ω effect), and vice versa.     

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.     

The Compton–Getting effect of energetic particles with an anisotropic pitch-angle distribution

This paper presents a simulation of anisotropy measurements by the low-energy charged particle (LECP) experiment on Voyager 1 for cases when the particle pitch-angle distribution function in the solar wind plasma reference frame is not isotropic. The model includes both the Compton–Getting anisotropy and perpendicular diffusion anisotropy that possibly exists in the upstream region of the termination shock. The results show that the Voyager 1 data cannot rule out either the model with zero solar wind speed or the one with a finite speed on qualitative basis. The determination of solar wind speed using the Compton–Getting effect is affected by the assumption of the magnetic field direction and perpendicular diffusion anisotropy. Because the pitch-angle distribution anisotropy is so large, a small uncertainty in the magnetic field direction can produce very different solar wind speeds ranging from ∼0 to >400 km/s. In fact, if the magnetic field is chosen to be in the Parker spiral direction, which is consistent with the magnetometer measurement on Voyager 1, the derived solar wind speed is still close to the supersonic value. Only the two lowest-energy channels of the LECP instrument may give a definitive answer to the solar wind speed. However, because these channels contain a very high level of cosmic ray background, an uncertainty of just a few percent in the background can entirely hamper the estimate of solar wind speed.     

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.     

F2 layer characteristics and electrojet strength over an equatorial station

The data presented in this work describes the diurnal and seasonal variation in hmF2, NmF2, and the electrojet current strength over an African equatorial station during a period of low solar activity. The F2 region horizontal magnetic element H revealed that the Solar quiet Sq(H) daily variation rises from early morning period to maximum around local noon and falls to lower values towards evening. The F2 ionospheric current responsible for the magnetic field variations is inferred to build up at the early morning hours, attaining maximum strength around 1200 LT. The Sq variation across the entire months was higher during the daytime than nighttime. This is ascribed to the variability of the ionospheric parameters like conductivity and winds structure in this region. Seasonal daytime electrojet (EEJ) current strength for June solstice, March and September equinoxes, respectively had peak values ranging within 27–35 nT (at 1400 LT) , 30–40 nT (at 1200 LT) and 35–45 nT (at 1500 LT). The different peak periods of the EEJ strength were attributed to the combined effects of the peak electron density and electric field. Lastly, the EEJ strength was observed to be higher during the equinoxes than the solstice period.     

Heliospheric implications of structure in the interstellar medium

Equilibrium models of diffuse interstellar material (ISM) near the Sun show a range of cloud densities, ionization, and temperatures which are consistent with data, although the local ISM must be inhomogeneous over ∼2 pc scales. The ISM close to the Sun has properties that are consistent with the sheetlike warm neutral (and partially ionized) gas detected in the Arecibo Millennium Survey. Local interstellar magnetic fields are poorly understood, but data showing weak polarization for nearby stars indicate dust may be trapped in fields or currents in the heliosheath nose region. Implications of this dust capture are widespread, and may impact the interpretation of the cosmic microwave background data. Observations of interstellar H⁰ inside of the solar system between 1975 and 2000 do not suggest any variation in the properties or structure of local interstellar H⁰ over distance scales of ∼750 AU to within the uncertainties.     

Evaluation of solar proton spectra using balloon cosmic ray observations and Monte Carlo simulation results

We have developed a method to evaluate the spectrum of solar energetic protons at the top of the Earth’s atmosphere from the measurements of our balloon cosmic ray experiment. By using the Monte Carlo PLANETOCOSMICS code based on Geant4 we compute the interaction of solar protons [10 MeV–10 GeV] with the Earth’s atmosphere. We obtain the angular and energy distributions of secondary particles (p, e⁻, e⁺, photons, muons) at different atmospheric levels as a function of primary proton spectra. By comparing the calculated depth dependence of the particle flux with the data obtained by our balloon experiment we can deduce the parameters of the solar proton spectrum that best fit the observations. In this paper we discuss our solar proton spectrum estimation method, and present results of its application to selected solar proton events from 2001 to 2005.     

Prediction of expected global climate change by forecasting of galactic cosmic ray intensity time variation in near future based on solar magnetic field data

A method of prediction of expected part of global climate change caused by cosmic ray (CR) by forecasting of galactic cosmic ray intensity time variation in near future based on solar activity data prediction and determined parameters of convection-diffusion and drift mechanisms is presented. This gave possibility to make prediction of expected part of global climate change, caused by long-term cosmic ray intensity variation. In this paper, we use the model of cosmic ray modulation in the Heliosphere, which considers a relation between long-term cosmic ray variations with parameters of the solar magnetic field. The later now can be predicted with good accuracy. By using this prediction, the expected cosmic ray variations in the near Earth space also can be estimated with a good accuracy. It is shown that there are two possibilities: (1) to predict cosmic ray intensity for 1–6 months by using a delay of long-term cosmic ray variations relatively to effects of the solar activity and (2) to predict cosmic ray intensity for the next solar cycle. For the second case, the prediction of the global solar magnetic field characteristics is crucial. For both cases, reliable long-term cosmic ray and solar activity data as well as solar magnetic field are necessary. For solar magnetic field, we used results of two magnetographs (from Stanford and Kitt Peak Observatories). The obtained forecasting of long-term cosmic ray intensity variation we use for estimation of the part of global climate change caused by cosmic ray intensity changing (influenced on global cloudiness covering).     

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.     

A search for the signature of microquasars in the cosmic ray iron spectrum measured by TIGER

In a recent paper Heinz and Sunyaev suggest that relativistic jets observed in microquasars might result in narrow features in the energy spectra of heavy cosmic rays with ≈1 to ≈10 GeV/nuc. They further argue that such features might be observable if there has been one or more microquasars nearby within the last few million years. We report preliminary results of a search for evidence of such features using data from a 32-day balloon flight of the Trans-Iron Galactic Element Recorder (TIGER). Although this flight took place near solar maximum, calculations of the broadening effects of solar modulation indicate that a narrow feature of sufficient intensity should still be observable. An energy spectrum for iron with high statistical significance has been derived from ≈100,000 Fe events in the energy range from about 2.5 to 10 GeV/nuc. Although our preliminary results do not reveal any obvious features, we will discuss the possibility of observing such features with TIGER and other instruments.     

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.     

Modeling and experimental study of the 27-day variation of galactic cosmic-ray intensity for a solar wind velocity depending on heliolongitude

We develop a three-dimensional (3-D) model of the 27-day variation of galactic cosmic-ray (GCR) intensity with a spatial variation of the solar wind velocity. A consistent, divergence-free interplanetary magnetic field is derived by solving the corresponding Maxwell equations with a variable solar wind speed, which reproduces in situ observed experimental data for the time interval to be analyzed (24 August 2007–28 February 2008). We perform model calculations for the GCR intensity using the variable solar wind and the corresponding magnetic field. Results are compatible with experimental data; the correlation coefficient between our model predictions and observed 27-day GCR variation is 0.80 ± 0.05.