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.     

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.     

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

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.     

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.     

Mean energy of response to galactic cosmic ray spectrum: IMP 8 and Climax neutron monitor

Solar modulations of galactic cosmic ray (GCR) intensity contain a wealth of information about their transport in the heliosphere. To extract this information from the data one studies the dependence of the observed modulations on the mean energy of response of detectors providing data for the analyses. There is a great deal of confusion about the detector energy response to GCR spectrum in the literature. We present a preliminary report on the computations of the mean energy of response for the Climax neutron monitor (CL/NM) and IMP 8 cosmic ray nuclear composition instrument to GCR protons for 1973–1998, covering the solar cycles 21 and 22. We find that for penetrating proton channel on IMP 8 the mean energy changes by a factor of over two whereas for the neutron monitor the change is only 21%. However, the corresponding change for the computed modulation function is a factor of about 3.5.     

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.     

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.     

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.     

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.     

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.     

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.     

宇宙线福布什下降的数值模拟

本文从银河宇宙线的太阳调制方程出发,认为激波对宇宙线的影响是由一扰动区产生的。在此扰动区中太阳风速度增加,扩散系数下降。由此进行了数值模拟,模拟结果表明:扩散系数的下降在产生福布什下降时要比太阳风速度的增加更为有效;福布什下降在近日球层内向外传播时,其幅度随径向距离的增大而衰减;两个无相互作用的激波同时存在时所产生的福布什下降为每个激波单独存在时的福布什下降的简单线性迭加。      

Rigidity dependence of Forbush decreases in the energy region exceeding the sensitivity of neutron monitors

Applicability of our present setup for solar modulation studies in a shallow underground laboratory is tested on four prominent examples of Forbush decrease during solar cycle 24. Forbush decreases are of interest in space weather application and study of energy-dependent solar modulation, and they have been studied extensively. The characteristics of these events, as recorded by various neutron monitors and our detectors, were compared, and rigidity spectrum was found. Linear regression was performed to find power indices that correspond to each event. As expected, a steeper spectrum during more intense extreme solar events with strong X-flares shows a greater modulation of galactic cosmic rays. Presented comparative analysis illustrates the applicability of our setup for studies of solar modulation in the energy region exceeding the sensitivity of neutron monitors.     

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.     

Cosmic ray anisotropy based on Yakutsk station in real time

The method has been developed to calculate galactic cosmic ray anisotropy parameters by using on-line data of the neutron monitor 24-NM-64 and muon telescope at the Yakutsk station. The preliminary analysis shows that characteristic changes in the anisotropy parameters caused by the first spherical harmonics of cosmic ray angular distribution are observed 1–2 days before the onset of the most part of large-scale geophysical disturbances on the Earth. There is reason to believe that the attraction of data of geophysical observations of other kinds will allow to develop the forecast methods for the arrival of large-scale interplanetary disturbances at the Earth.     

Energy dependence of the rigidity spectrum of Forbush decrease of galactic cosmic ray intensity

We show that rigidity spectrum of Forbush decrease (Fd) of galactic cosmic ray (GCR) intensity in September 9–23, 2005 clearly depends on energy. We calculated rigidity spectrum of the Fd based on the neutron monitors and Nagoya muon telescope channels’ data divided in three groups according to their cut off rigidities. We found that temporal changes of rigidity spectrum exponent γ are approximately similar for all cut off rigidity groups, but γ values are the larger the higher are cut off rigidities. We conclude that rigidity spectrum of Fd is hard for lower energy range and is soft for the higher energy range. We believe that an energy dependence of the power law rigidity spectrum of Fd is observed owing to the preferential convection–diffusion mechanism during Fd in September 9–23, 2005. It is a reflection of an influence of the temporal changes of the structure of the interplanetary magnetic field (IMF) turbulence in different range of frequency f during Fd. Particularly, a decisive role in formation of the character of the rigidity spectrum belongs to the changes of the exponent ν of the power spectral density (PSD) of the IMF turbulence (PSD ∝ f^−ν). The exponent ν is greater for high frequency region of the IMF turbulence (responsible for scattering of low rigidity particles of GCR), than for low frequency region of the IMF turbulence (being responsible for scattering of higher rigidity particles). Also, we challenge to estimate an existence of slab/2D structure of solar wind turbulence during the Fd in September 9–23, 2005 based on the distribution of average turbulence energy among the IMF’s components.