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.     

Space storm measurements of the July 2005 solar extreme events from the low corona to the Earth

The Athens Neutron Monitor Data Processing (ANMODAP) Center recorded an unusual Forbush decrease with a sharp enhancement of cosmic ray intensity right after the main phase of the Forbush decrease on 16 July 2005, followed by a second decrease within less than 12 h. This exceptional event is neither a ground level enhancement nor a geomagnetic effect in cosmic rays. It rather appears as the effect of a special structure of interplanetary disturbances originating from a group of coronal mass ejections (CMEs) in the 13–14 July 2005 period. The initiation of the CMEs was accompanied by type IV radio bursts and intense solar flares (SFs) on the west solar limb (AR 786); this group of energetic phenomena appears under the label of Solar Extreme Events of July 2005. We study the characteristics of these events using combined data from Earth (the ARTEMIS IV radioheliograph, the Athens Neutron Monitor (ANMODAP)), space (WIND/WAVES) and data archives. We propose an interpretation of the unusual Forbush profile in terms of a magnetic structure and a succession of interplanetary shocks interacting with the magnetosphere.     

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.     

Evolution of Dst index,cosmic rays and global temperature during solar cycles 20–23

We have studied conditions in interplanetary space, which can have an influence on galactic cosmic ray (CR) and climate change. In this connection the solar wind and interplanetary magnetic field parameters and cosmic ray variations have been compared with geomagnetic activity represented by the equatorial Dst index from the beginning 1965 to the end of 2012. Dst index is commonly used as the solar wind–magnetosphere–ionosphere interaction characteristic. The important drivers in interplanetary medium which have effect on cosmic rays as CMEs (coronal mass ejections) and CIRs (corotating interaction regions) undergo very strong changes during their propagation to the Earth. Because of this CMEs, coronal holes and the solar spot numbers (SSN) do not adequately reflect peculiarities concerned with the solar wind arrival to 1 AU. Therefore, the geomagnetic indices have some inestimable advantage as continuous series other the irregular solar wind measurements. We have compared the yearly average variations of Dst index and the solar wind parameters with cosmic ray data from Moscow, Climax, and Haleakala neutron monitors during the solar cycles 20–23. The descending phases of these solar cycles (CSs) had the long-lasting solar wind high speed streams occurred frequently and were the primary contributors to the recurrent Dst variations. They also had effects on cosmic rays variations. We show that long-term Dst variations in these solar cycles were correlated with the cosmic ray count rate and can be used for study of CR variations. Global temperature variations in connection with evolution of Dst index and CR variations is discussed.     

Sunspot activity and cosmic ray modulation at 1 a.u. for 1900–2013

The descent of sunspot cycle 23 to an unprecedented minimum of long duration in 2006–2009 led to a prolonged galactic cosmic ray (GCR) recovery to the highest level observed in the instrumental era for a variety of energetic charged particle species on Earth, over a wide range of rigidities. The remarkable GCR increase measured by several ground-based, balloon-borne, and detectors on a satellite is described and discussed. It is accompanied by a decrease in solar wind velocity and interplanetary magnetic field at 1 a.u., reaching the lowest values since measurements of the solar wind began in October 1963; the solar polar field strength (μT) measured at the Wilcox Solar Observatory (WSO) is also significantly reduced compared to prior cycles since the start of the program in 1976, the polar field in the northern hemisphere reversed in June 2012 and again in February 2014, that in the southern hemisphere reversed in July 2013. If updates of WSO data confirm the second reversal in northern solar hemisphere, it would pose a serious challenge to the Dynamo Theory. The long-term change in solar behavior may have begun in 1992, perhaps earlier. The physical underpinnings of these solar changes need to be understood and their effect on GCR modulation processes clarified. The study discusses the recent phenomena in the context of GCR modulation since 1900. These happenings affected our empirical predictions for the key parameters for the next two sunspot cycles (they may be progressively less active than sunspot cycle 24) but it enhanced support for our prediction that solar activity is descending into a Dalton-like grand minimum in the middle of the twentyfirst century, reducing the frequency of the coronal mass ejections; they determine the space weather affecting the quality of life on Earth, radiation dose for hardware and human activities in space as well as the frequency of large Forbush decreases at 1 a.u.     

Cosmic ray variations of solar origin in relation to human physiological state during the December 2006 solar extreme events

There is an increasing amount of evidence linking biological effects to solar and geomagnetic disturbances. A series of studies is published referring to the changes in human physiological responses at different levels of geomagnetic activity. In this study, the possible relation between the daily variations of cosmic ray intensity, measured by the Neutron Monitor at the Cosmic Ray Station of the University of Athens (http://cosray.phys.uoa.gr) and the average daily and hourly heart rate variations of persons, with no symptoms or hospital admission, monitored by Holter electrocardiogram, is considered. This work refers to a group of persons admitted to the cardiological clinic of the KAT Hospital in Athens during the time period from 4th to 24th December 2006 that is characterized by extreme solar and geomagnetic activity. A series of Forbush decreases started on 6th December and lasted until the end of the month and a great solar proton event causing a Ground Level Enhancement (GLE) of the cosmic ray intensity on 13th December occurred. A sudden decrease of the cosmic ray intensity on 15th December, when a geomagnetic storm was registered, was also recorded in Athens Neutron Monitor station (cut-off rigidity 8.53 GV) with amplitude of 4%. It is noticed that during geomagnetically quiet days the heart rate and the cosmic ray intensity variations are positively correlated. When intense cosmic ray variations, like Forbush decreases and relativistic proton events produced by strong solar phenomena occur, cosmic ray intensity and heart rate get minimum values and their variations, also, coincide. During these events the correlation coefficient of these two parameters changes and follows the behavior of the cosmic ray intensity variations. This is only a small part of an extended investigation, which has begun using data from the year 2002 and is still in progress.     

Ionization effect of solar protons in the Earth atmosphere – Case study of the 20 January 2005 SEP event

The cosmic ray ground level enhancement on January 20, 2005 is among the largest recorded events in the history of cosmic ray measurements. The solar protons of MeV energies cause an excess of ionization in the atmosphere, specifically over polar caps following major solar disturbances. The ionization effect in the Earth atmosphere is obtained for various latitudes on the basis of solar proton energy spectra, reconstructed from GOES 11 measurements and subsequent full Monte Carlo simulation of cosmic ray induced atmospheric cascade. The estimation of ionization rates is based on a numerical model for cosmic ray induced ionization. The evolution of atmospheric cascade is performed with the CORSIKA 6.52 code using FLUKA 2006b and QGSJET II hadron interaction models. The atmospheric ion rate ionization is explicitly obtained for 40°N, 60°N and 80°N latitudes. The time evolution of obtained ion rates is presented. It is demonstrated that ionization effect is negative for 40°N and small for 60°N, because of accompanying Forbush decrease. The ionization effect is significant only in sub-polar and polar atmosphere during the major ground level enhancement of 20 January 2005.     

Fluctuations of cosmic rays and IMF in the vicinity of interplanetary shocks

Fluctuations of cosmic rays and interplanetary magnetic field upstream of interplanetary shocks are studied using data of ground-based polar neutron monitors as well as measurements of energetic particles and solar wind plasma parameters aboard the ACE spacecraft. It is shown that coherent cosmic ray fluctuations in the energy range from 10 keV to 1 GeV are often observed at the Earth’s orbit before the arrival of interplanetary shocks. This corresponds to an increase of solar wind turbulence level by more than the order of magnitude upstream of the shock. We suggest a scenario where the cosmic ray fluctuation spectrum is modulated by fast magnetosonic waves generated by flux of low-energy cosmic rays which are reflected and/or accelerated by an interplanetary shock.     

Experimental data and analysis of the October 2003 Forbush decrease

On October 28, 2003 an Earthward-directed coronal mass ejection (CME) was observed from SOHO/LASCO imagery in conjunction with an X17 solar flare. The CME, traveling at nearly 2000 km/s, impacted the Earth on October 29, 2003 causing ground-based particle detectors to register a counting rate drop known as a Forbush decrease. The CME was not only responsible for affecting the rate of cosmic rays, but also caused anisotropies in their direction of incidence. Data from Project GRAND, an array of proportional wire chambers which detects secondary muons, are presented during the time of this Forbush decrease.     

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.     

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

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

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.     

Effect of Forbush decreases of cosmic ray flux on ozone at higher middle latitudes

The Forbush decreases of cosmic ray flux occur prevailingly together with geomagnetic storms, because these phenomena have a similar origin in solar/interplanetary processes. To study the effects of large Forbush decreases on total ozone at middle latitudes, we use the TOMS total ozone data along latitudinal circles 40°N and 50°N. The effects of Forbush decreases are found to occur or to be non-measurable under the same conditions as those of geomagnetic storms: certain effect occurs only at 50°N (not 40°N), in winter, under conditions of high solar activity and the east phase of the QBO. However, the effects of the analyzed Forbush decreases are weaker than the effects of strong geomagnetic storms.     

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.     

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.     

宇宙线强度变化与磁扰K类型

本文把1966—1983年期间发生的679个地磁暴进行了分类,利用统计方法分析了各类磁暴发生前后宇宙线强度的变化特征.突发急始脉冲发生后,宇宙线强度没有出现显著的变化;缓始型暴发生后,宇宙线强度出现Forbush下降,但下降幅度很小;急始型暴发生后,宇宙线强度出现十分明显的Forbush下降.当把急始型暴按K指数大小和持续时间分为5种类型,发现它们伴随的宇宙线Forbush下降是不一样的,其下降幅度随磁暴的增强而加大,下降的速率随磁暴的增强而加快,扰动的持续时间随磁暴的减弱而增加.      

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.     

Large geomagnetic storms of extreme solar event periods in solar cycle 23

During extreme solar events such as big flares or/and energetic coronal mass ejections (CMEs) high energy particles are accelerated by the shocks formed in front of fast interplanetary coronal mass ejections (ICMEs). The ICMEs (and their sheaths) also give rise to large geomagnetic storms which have significant effects on the Earth’s environment and human life. Around 14 solar cosmic ray ground level enhancement (GLE) events in solar cycle 23 we examined the cosmic ray variation, solar wind speed, ions density, interplanetary magnetic field, and geomagnetic disturbance storm time index (D_st). We found that all but one of GLEs are always followed by a geomagnetic storm with D_st  −50 nT within 1–5 days later. Most(10/14) geomagnetic storms have D_st index  −100  nT therefore generally belong to strong geomagnetic storms. This suggests that GLE event prediction of geomagnetic storms is 93% for moderate storms and 71% for large storms when geomagnetic storms preceded by GLEs. All D_st depressions are associated with cosmic ray decreases which occur nearly simultaneously with geomagnetic storms. We also investigated the interplanetary plasma features. Most geomagnetic storm correspond significant periods of southward B_z and in close to 80% of the cases that the B_z was first northward then turning southward after storm sudden commencement (SSC). Plasma flow speed, ion number density and interplanetary plasma temperature near 1 AU also have a peak at interplanetary shock arrival. Solar cause and energetic particle signatures of large geomagnetic storms and a possible prediction scheme are discussed.