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 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.     

Applications and usage of the real-time Neutron Monitor Database

A high-time resolution Neutron Monitor Database (NMDB) has started to be realized in the frame of the Seventh Framework Programme of the European Commission. This database will include cosmic ray data from at least 18 neutron monitors distributed around the world and operated in real-time. The implementation of the NMDB will provide the opportunity for several research applications most of which will be realized in real-time mode. An important one will be the establishment of an Alert signal when dangerous solar cosmic ray particles are heading to the Earth, resulting into ground level enhancements effects registered by neutron monitors. Furthermore, on the basis of these events analysis, the mapping of all ground level enhancement features in near real-time mode will provide an overall picture of these phenomena and will be used as an input for the calculation of the ionization of the atmosphere. The latter will be useful together with other contributions to radiation dose calculations within the atmosphere at several altitudes and will reveal the absorbed doses during flights. Moreover, special algorithms for anisotropy and pitch angle distribution of solar cosmic rays, which have been developed over the years, will also be set online offering the advantage to give information about the conditions of the interplanetary space. All of the applications will serve the needs of the modern world which relies at space environment and will use the extensive network of neutron monitors as a multi-directional spectrographic detector. On top of which, the decreases of the cosmic ray intensity – known as Forbush decreases – will also be analyzed and a number of important parameters such as galactic cosmic ray anisotropy will be made available to the users of NMDB. A part of the NMDB project is also dedicated to the creation of a public outreach website with the scope to inform about cosmic rays and their possible effects on humans, technological systems and space-terrestrial environment. Therefore, NMDB will also stand as an informative gate on space research through neutron monitor’s data usage.     

New index of long-term variations of galactic cosmic ray intensity

We find that the soft rigidity spectrum of the Galactic Cosmic Ray (GCR) intensity variations for the maximum epoch and the hard rigidity spectrum for the minimum epoch calculated based on the neutron monitors experimental data (1960–2002) are related with the various dependence of the diffusion coefficient on the GCR particle’s rigidity for different epoch of solar activity. This dependence is stronger in the maximum epoch than in the minimum epoch of solar activity, and is provided by the essential temporal rearrangements of the structure of the Interplanetary Magnetic Field (IMF) turbulence from the maxima to minima epoch of solar activity. We also show that the rigidity spectrum of GCR intensity variations is harder for the effective rigidities ∼(10–15) GV (by neutron monitors data), than for the effective rigidities ∼(25–30) GV (by neutron monitors and muon telescopes data). A general scenario of GCR modulation versus solar activity is settled on the essential temporal rearrangements of the structure of the IMF turbulence. Therefore, the temporal changes of the power law rigidity spectrum exponent can be considered as a vital (new) index to explain the 11-year variations of the GCR intensity. We assume that ∼(70–80)% of the changes of the amplitudes of the 11-year variations of GCR intensity is related with the changes of the IMF turbulence versus solar activity.     

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.     

Forbush decreases and solar events seen in the 10–20 GeV energy range by the Karlsruhe Muon Telescope

Since 1993, a muon telescope located at Forschungszentrum Karlsruhe (Karlsruhe Muon Telescope) has been recording the flux of single muons mostly originating from primary cosmic-ray protons with dominant energies in the 10–20 GeV range. The data are used to investigate the influence of solar effects on the flux of cosmic rays measured at Earth. Non-periodic events like Forbush decreases and ground level enhancements are detected in the registered muon flux. A selection of recent events will be presented and compared to data from the Jungfraujoch neutron monitor. The data of the Karlsruhe Muon Telescope help to extend the knowledge about Forbush decreases and ground level enhancements to energies beyond the neutron monitor regime.     

Modeling and data analysis of a Forbush decrease

We study the Forbush decrease of the galactic cosmic ray intensity observed in 9–25 September 2005 using the experimental data and a newly developed time-dependent three dimensional modeling. We analyze neutron monitors and muon telescopes, and the interplanetary magnetic field data. We demonstrate a clear relationship between the rigidity (R) spectrum exponent (γ) of the Forbush decrease and the exponent (ν) of the power spectral density of the components of the interplanetary magnetic field in the frequency range of ∼ 10⁻⁶–10 ⁻⁵ Hz. We confirm that an inclusion of the time-dependent changes of the exponent ν makes the newly developed nonstationary three dimensional model of the Forbush decrease compatible with the experimental data. Also, we show that the changes of the rigidity spectrum exponent γ does not depend on the level of convection of the galactic cosmic rays stream by solar wind; depending on the changes of the exponent ν, i.e. on the state of the turbulence of the interplanetary magnetic field.     

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.     

Relativistic solar protons in the ground level event of 23 February 1956: New study

In this study we applied again to the outstanding solar particle event of 23 February 1956, the largest one in the entire history of observations of solar cosmic rays. Due to significant improvement of the analysis/modeling techniques and new understanding of physical processes in the solar atmosphere and interplanetary space, a possibility arises to interpret the old data in the light of modern concept of multiple particle acceleration at/near the Sun. In our new analysis the data of available then neutron monitors and muon telescopes are used. The technique of the analysis includes: (a) calculation of asymptotic cones of ground-based detectors; (b) modeling of cosmic ray detector responses at variable parameters of the flux of solar relativistic protons; (c) determination of primary solar proton parameters outside magnetosphere by comparison of computed responses with observations. Certain evidence was obtained that the flux of relativistic solar protons consisted of two distinct components: prompt and delayed ones. The prompt component with exponential energy spectrum caused a giant impulse-like increase at a number of European cosmic ray stations. The delayed component had a power-law spectrum and was a cause of gradual increase at cosmic ray stations in the North American region. A numerical simulation of the proton acceleration in the vicinity of the magnetic reconnection region brings to the proton spectrum with exponential dependence on energy. This agrees with observational data for the prompt component. It is also shown that the huge increase in ∼5000% on neutron monitors was due to the prompt component only with the exponential proton spectrum. The power-law spectrum of comparable intensity gave considerably smaller effect.     

Two intense Forbush decreases of solar activity cycle 22

We study two intense Forbush decreases that occurred during two adjacent SOLTIP (Solar connection of Transient Interplanetary Processes) intervals; namely SOLTIP 1 (22–27 March 1991) and SOLTIP 2 (1–17 June 1991); galactic cosmic ray intensity at the depth of the second Forbush decrease was the lowest ever recorded since continuous monitoring by Climax neutron monitor began in 1951 (58% below the solar minimum value of 1954), indicating extreme conditions in the heliosphere that prevented galactic cosmic rays from reaching the Earth. These decreases were seen propagating in outer heliosphere by the deep space missions Voyagers 1, 2 and Pioneer 10, 11, with suitable time delays. We analyze hourly, pressure corrected, neutron monitor data from the global sites in both hemispheres, and muon telescopes located underground; they respond to 10–300 GV range of the galactic cosmic ray spectrum. This circumstance provides us an ideal opportunity to study the rigidity dependence of the amplitudes of the two Forbush decreases. In both cases the amplitude is found to be a power law in rigidity, with negative exponents.     

New insights from cross-correlation studies between solar activity indices and cosmic-ray flux during Forbush decrease events

Observed galactic cosmic ray intensity can be subjected to a transient decrease. These so-called Forbush decreases are driven by coronal mass ejection induced shockwaves in the heliosphere. By combining in situ measurements by space borne instruments with ground-based cosmic ray observations, we investigate the relationship between solar energetic particle flux, various solar activity indices, and intensity measurements of cosmic rays during such an event. We present cross-correlation study done using proton flux data from the SOHO/ERNE instrument, as well as data collected during some of the strongest Forbush decreases over the last two completed solar cycles by the network of neutron monitor detectors and different solar observatories. We have demonstrated connection between the shape of solar energetic particles fluence spectra and selected coronal mass ejection and Forbush decrease parameters, indicating that power exponents used to model these fluence spectra could be valuable new parameters in similar analysis of mentioned phenomena. They appear to be better predictor variables of Forbush decrease magnitude in interplanetary magnetic field than coronal mass ejection velocities.     

Case study of Forbush decreases: Energy dependence of the recovery

Case study is presented for three Forbush decreases in 2004–2005, using cosmic ray data from ground-based detectors – neutron monitors and a muon detector. One of them was a typical event (September 2005), while two other were quite unusual (November 2004 and January 2005). Two unusual features, not expected from the standard theory, are revealed: (1) the recovery time of a Forbush decrease can strongly depend on the energy; (2) an over-recovery is observed in the most energetic cosmic ray data (muon detector). A simple scenario is suggested for the observed phenomenon.     

Modeling the solar cosmic ray event of 13 December 2006 using ground level neutron monitor data

In order to understand the physics under extreme solar conditions such as those producing ground level enhancements of solar cosmic rays, it is important to use accurate and reliable models. The NM-BANGLE Model is a new cosmic ray model which couples primary solar cosmic rays at the top of the Earth’s atmosphere with the secondary ones detected at ground level by neutron monitors during GLEs. This model calculates the evolution of several GLE parameters such as the solar cosmic ray spectrum, anisotropy and particle flux distribution, revealing crucial information on the energetic particle propagation and distribution. The total output of the NM-BANGLE Model is a multi-dimensional GLE picture that gives an important contribution to revealing the characteristics of solar energetic particle events recorded at ground level. In this work, the results of the NM-BANGLE Model application to the recent GLE of 13 December 2006 are presented and discussed. Moreover, a comparison with the extreme event of 20 January 2005 (GLE69) has been realized.     

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.     

Investigation of diurnal variations of cosmic ray fluxes measured with using ASEC and NMDB monitors

A study of daily variations of secondary Cosmic Rays (CR) is performed using data on charged and neutral CR fluxes. Particle detectors of Aragats Space-Environmental Center (ASEC), Space Environmental Viewing and Analysis Network (SEVAN) and neutron monitors of the Neutron Monitor Database (NMDB) are used. ASEC detectors continuously register various species of secondary CR with different threshold energies and incident angles. NMDB joins data of 12 Eurasian neutron monitors. Data at the beginning of the 24th solar activity cycle are used to avoid biases due to solar transient events and to establish a benchmark for the monitoring of solar activity in the new started solar cycle.     

Unexpected burst of solar activity recorded by neutron monitors during October–November 2003

During the extreme burst of solar activity in October–November 2003, a series of outstanding events distinguished by their magnitude and peculiarities were recorded by the ground based neutron monitor network. The biggest and most productive in 23rd solar cycle active region 486 generated the most significant series of solar flares among of which the flare X28/3B on November 4, 2003 was the mostly powerful over the history of X-ray solar observations. The fastest arrival of the interplanetary disturbance from the Sun after the flare event in August 1972 and the highest solar wind velocity and IMF intensity were observed during these events. In one-week period three ground level enhancements (GLEs) of solar cosmic rays were recorded by neutron monitor network (28, 29 October and 2 November 2003). Maximum proton energy in these events seems to be ranged from 5 to 10 GeV. Joint analysis of data from ground level stations (neutron monitors) and satellite measurements allows the estimation of the particle path length, the onset time of the injection on the Sun and some other proton flux characteristics.     

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.     

23rd solar cycle: Solar activity parameters and associated Forbush decreases

We report Forbush decreases (FD) in cosmic ray intensity from January 1996 to December 2008, the whole Solar Cycle 23rd. Statistical analysis is done for only 152 events for which associated solar flare position, flare classes, and Coronal Mass Ejections (CME) speed are given. We applied FD parameters taken from the Forbush Effects and Interplanetary Disturbances databases maintained by the Pushkov Institute of Terrestrial Magnetism, Ionosphere, and Radiowave Propagation (IZMIRAN), obtained by processing the data of the worldwide neutron monitor network using the global survey method (GSM) (A. Belov et al., 2018). For the said number of events, we examine their effect on interplanetary space and the decrease of the galactic cosmic rays (GCR) near Earth. We found that the 11–20° latitudinal belt shows more FD- associated flare events than the other latitudinal belts, and on this belt, the Southern hemisphere is more active. The results reveal that FDs and solar flares are well correlated. Statistical analysis is carried out for the magnitude of the CR decrease with solar and geomagnetic parameters.     

Neutron monitor asymptotic directions of viewing during the event of 13 December 2006

During the recent ground level enhancement of 13 December 2006, also known as GLE70, solar cosmic ray particles of energy bigger that ∼500 MeV/nucleon propagated inside the Earth’s magnetosphere and finally accessed low-altitude satellites and ground level neutron monitors. The magnitude and the characteristics of this event registered at different neutron monitor stations of the worldwide network can be interpreted adequately on the basis of an estimation of the solar particle trajectories in the near Earth interplanetary space. In this work, an extended representation of the Earth’s magnetic field was realized applying the Tsyganenko 1989 model. Using a numerical back-tracing technique the solar proton trajectories inside the magnetospheric field of the Earth were calculated for a variety of particles, initializing their travel at different locations, covering a wide range of energies. In this way, the asymptotic directions of viewing were calculated for a significant number of neutron monitor stations, providing crucial information on the Earth’s “magnetospheric optics” for primary solar cosmic rays, on the top of the atmosphere, during the big solar event of December 2006. The neutron monitor network has been treated, therefore, as a multidimensional tool that gives insights into the arrival directions of solar cosmic ray particles as well as their spatial and energy distributions during extreme solar events.