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.     

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.     

Relationships between neutron fluxes and rain flows

Registration of secondary cosmic ray neutrons is a convenient tool for investigation of primary cosmic ray variations and meteorological effects as well. At present a large network of neutron monitors exists, providing the studies of cosmic ray variations related to the interplanetary conditions and geomagnetic activity. At the same time cosmic ray variations may be caused by some atmospheric processes. In this connection, using the data from standard and lead-free neutron monitors, and gamma and muon detectors, we studied relations between rain flows and neutron, gamma and ionization component behavior. To explain observable results the calculations of neutron and gamma absorption and albedo neutron spectra have been performed on the basis of universal software package FLUKA-2006. In this study we used hourly data on the neutron flux, corrected for barometric pressure and data from local meteorological stations. It was shown that secondary neutron radiation, recorded by lead-free NM, and gamma radiation as well are strongly effected by meteorological factors. The neutron component behavior depends on the moisture content in the soil, and above its surface.     

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.     

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.     

Galactic cosmic ray flux simulation and prediction.

A dynamic galactic cosmic ray model is proposed to quantitatively describe the z=1-28 ions and electrons of E=10-10(5) MeV/nucleon and their particle flux variations around the Earth's orbit and beyond the Earth's magnetosphere due to diverse large-scale variations of solar activity factors. The variations of large-scale heliospheric magnetic fields and the galactic cosmic ray flux variation time delays relative to solar activity variations are simulated. The lag characteristics and sunspot number predictions having been determined in detail, the model can be used to predict galactic cosmic ray flux levels.     

多方向μ介子通量一致性特征研究

通过地面μ介子望远镜可以探测不同方向到达的宇宙线通量,得到从外空间入射的宇宙线受到不均匀结构的调制情况,判断CME的特征.通过分析,发现大地磁暴前Nagoya台站东向和南向的探测数据存在固定的2h时间差,认为这是由两个入射方向的宇宙线粒子先后穿越CME结构引起的.分别计算了两个方向相同时间和南向相位后移2h后通量探测数据的相关系数,以及两种情形下通量差的变化幅度,定量描述了CME接近地球过程中两个方向通量的相关特征.通过比较发现,CME接近地球过程中,经过相位变换的两个方向的相关系数明显高于未经变换的情况,经过相位变换的两个方向的通量差幅度明显小于未经变换的情况;CME到达地球后,两种情形的相关系数和通量差幅度则趋于相同.对2003-2005年K_p=9的地磁暴事件的分析均发现了这种现象.对2006年12月14日大地磁暴前的μ介子通量特征进行了分析,也完全符合上述特征.      

Calculated cosmic rays dynamics before the geomagnetic storms onset caused by the interplanetary shocks

Dynamics of the galactic cosmic ray intensity caused by their interactions with a shock front in the September 9, 1992 event has been determined. Corresponding variations of the cosmic ray intensity have been calculated for different stations of the world network of neutron monitors and muon telescopes of stations Nagoya and Sakashita. Comparison calculated results with observational data shows, in general, satisfactory consensus both on amplitude and in time. The developed method can be used for investigation dynamics of the solar wind disturbances precursors in the cosmic rays.     

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.     

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 local time dependence of the anisotropic solar cosmic ray flux.

The distribution of the solar cosmic radiation flux over the earth is not uniform, but the result of complex phenomena involving the interplanetary magnetic field, the geomagnetic field and latitude and longitude of locations on the earth. The latitude effect relates to the geomagnetic shield; the longitude effect relates to local time. For anisotropic solar cosmic ray events the maximum particle flux is always along the interplanetary magnetic field direction, sometimes called the Archimedean spiral path from the sun to the earth. During anisotropic solar cosmic ray event, the locations on the earth viewing "sunward" into the interplanetary magnetic field direction will observe the largest flux (when adjustments are made for the magnetic latitude effect). To relate this phenomena to aircraft routes, for anisotropic solar cosmic ray events that occur during "normal quiescent" conditions, the maximum solar cosmic ray flux (and corresponding solar particle radiation dose) will be observed in the dawn quadrant, ideally at about 06 hours local time.     

History of cosmic ray research in Finland

The history of cosmic ray research in Finland can be traced back to the end of 1950s, when first ground-based cosmic ray measurements started in Turku. The first cosmic ray station was founded in Oulu in 1964 performing measurements of cosmic rays by a muon telescope, which was later complemented by a neutron monitor. Since the 1990s, several research centers and universities, such as The Finnish Meteorological Institute, Helsinki University of Technology, University of Oulu, University of Turku and University of Helsinki have been involved in space science projects, such as SOHO, AMS, Cluster, Cassini, BepiColombo, etc. At the same time, ground-based cosmic ray measurements have reached a new level, including a fully automatic on-line database in Oulu and a new muon measuring underground site in Pyhäsalmi. Research groups in Helsinki, Oulu and Turku have also extensive experience in theoretical investigations of different aspects of cosmic ray physics. Cosmic ray research has a 50-year long history in Finland, covering a wide range from basic long-running ground-based observations to high-technology space-borne instrumentation and sophisticated theoretical studies. Several generations of researchers have been involved in the study ensuring transfer of experience and building the recognized Finnish research school of cosmic ray studies.     

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.     

Estimation of long-term cosmic ray intensity variation in near future and prediction of their contribution in expected global climate change

On the basis of results obtained in our paper [Dorman, L.I. Long-term cosmic ray intensity variation and part of global climate change, controlled by solar activity through cosmic rays, Paper D2.1/C2.2/E3.1-0097-04. Adv. Space Res., 2004 (accepted)], we determine: the dimension of the Heliosphere (modulation region), radial diffusion coefficient and other parameters of convection–diffusion; drift mechanisms of long-term variations of cosmic ray (CR) dependence on particle energy; level of solar activity (SA); and generally, the solar magnetic field. We obtain this important information on the basis of CR and SA data in the past, taking into account the theory of convection–diffusion and global drift modulation of galactic CR in the Heliosphere. By using these results and other regularly published predictions of expected SA variation in the near future, as well as predictions of the next SA cycle, we may make predictions of long-term cosmic ray intensity variation expected in the near future (up to 10–12 years). In [Dorman, L.I. Long-term cosmic ray intensity variation and part of global climate change, controlled by solar activity through cosmic rays, Paper D2.1/C2.2/E3.1-0097-04. Adv. Space Res., 2004 (accepted)], properties of connections between long-term variation in CR intensity and some part of a global climate change were estimated, controlled by solar activity through CR. We show that in this way we may make predictions of some part of a global climate change expected in the near future (up to 10–12 years and maybe more, depending upon the period during which definite predictions of SA can be made), controlled by solar activity through CR. In this case, estimations of expected long-term changes in the planetary distribution of cutoff rigidities, which also influence CR intensity, as well as CR-influenced effects on global climate variation, become important.     

Short-term periodicities in the downward longwave radiation and their associations with cosmic ray and solar interplanetary data

In this study downward longwave (LW) atmospheric radiation data for the period of 2014–2020 were used to search for short-term periodicities using fast Fourier transform (FFT). Several local peaks in the power spectrum density were found and established. The time series exhibits a series of significant peaks (exceeding the 95% confidence limit), such as at 273 days, 227 days, 200 days, 178 days, 157 days, 110 days, 120 days, 87 days, 73 days, 53–56 days, 35–30 days, 25–27 days, 21 days, 13 days, and 9–10 days.Moreover, cosmic ray data from KACST muon detector and the Oulu neutron monitor, as well as the data for the solar radio flux at 10.7 cm (F10.7 cm), Dst index, and solar wind speed for the same period as the LW data, were used to look for common cyclic variations and periodicities matching those found in the LW radiation. This was done to investigate the possible effect of the solar activity parameters on LW radiation. Several common periodicities were observed in the spectra of all the variables considered, such as 227 days, 154–157 days, 25–27 days, and 21 days. Some of the periodicities found in the LW radiation spectrum can be attributed to the modulation of the cosmic ray intensity by solar activity. Others are attributed to the disturbances in the interplanetary magnetic field. Based on the spectral results, we suggest that the solar signals may directly or indirectly affect the variations of the downward longwave radiation, which in turn may affect climate change.     

CME-geometry and cosmic-ray anisotropy observed by a prototype muon detector network

We analyze the cosmic-ray anisotropy observed by a prototype network of muon detectors during geomagnetic storms associated with coronal mass ejections (CMEs). The network currently consists of multidirectional surface muon detectors at Nagoya (Japan) and Hobart (Australia), together with a prototype detector at São Martinho (Brazil) which has been in operation since March, 2001. In this report, we analyze the anisotropy recorded in both the muon detector and neutron monitor (the Spaceship Earth) networks and find significant enhancements of cosmic-ray anisotropy during geomagnetic storms. Following the analysis by Bieber and Evenson [Bieber, J.W., Evenson, P. CME geometry in relation to cosmic ray anisotropy. Geophys. Res. Lett. 25 (1998) 2955–2958] for the neutron monitor data at 10 GeV, we also derive cosmic-ray density gradients from muon data at higher-energy (50 GeV), possibly reflecting the larger-scale geometry of CMEs causing geomagnetic storms. We particularly find in some events the anisotropy enhancement clearly starting prior to the storm onset in both the muon and neutron data. This is the first result of the CME-geometry derived from simultaneous observations of the anisotropy with networks of multidirectional muon detectors and neutron monitors.     

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.     

Test alert service against very large SEP Events

The Aragats Solar Environment Center provides real time monitoring of different components of secondary cosmic ray fluxes. We plan to use this information to establish an early warning alert system against extreme, very large solar particle events with hard spectra, dangerous for satellite electronics and for the crew of the Space Station. Neutron monitors operating at altitude 2000 and 3200 m are continuously gathering data to detect possible abrupt variations of the particle count rates. Additional high precision detectors measuring muon and electron fluxes, along with directional information are under construction on Mt. Aragats. Registered ground level enhancements, in neutron and muon fluxes along with correlations between different species of secondary cosmic rays are analyzed to reveal possible correlations with expected times of arrival of dangerous solar energetic particles.     

南极中山站宇宙线观测

2014年由两个具有0.5m×0.5m闪烁体面积的探测器组成的宇宙线μ子望远镜在南极中山站建成,观测数据通过卫星链路传回并发布.通过对中山站μ子望远镜数据的气象效应分析,发现中山站的气压变化对宇宙线计数率影响显著并呈负相关.计算得到宇宙线高能质子在中山站的垂直截止刚度R_C为0.076GV,对比中国西藏羊八井和北京高能质子的垂直截止刚度,中山站非常有利于对太阳高能质子事件的观测.      

Possible influence of cosmic rays on climate through thunderstorm clouds

With decreasing of cosmic ray (CR) intensity caused by increasing of solar activity (SA) or in some short periods of Forbush-decreases, the intensity of secondary CR relativistic electrons decreases and the probability of formation of thunderstorm clouds and discharges between clouds or between clouds and ground is also expected to decrease. This will influence on weather and climate. In this case is very important to have more detail information on the atmospheric electric field distribution in the atmosphere, additional to information what gave now electric field sensors (EFS) only in about one point near the ground. We show that CR not only influenced on atmospheric electric field phenomenon, but can give practically continuous information on the atmospheric electric field distribution in the atmosphere. We extend our theory of CR atmospheric electric field effect on electron–photon, muon and neutron component including different multiplicities. We take into account that about 0.07 of neutron monitor counting rate caused by negative soft muons captured by lead nucleons and formed mesoatoms with generation of several MeV energy neutrons from lead. In this case the neutron monitor or neutron super-monitor works as analyzer that detects muons of only one, negative sign. It is very important because the atmospheric electric field effect have opposite signs for positive and negative muons that main part of this effect in the muon telescope or in ionization chamber is compensated and we can observe only small part of total effect of one sign muons. On the basis of our general theory of CR atmospheric electric field effects with taking into account of negative soft muon acceleration and deceleration in the Earth atmosphere (in dependence of direction and intensity of electric field) we discuss the possibility of existing this effect in CR neutron monitor counting rate and in different multiplicities and calculate the expected effects in dependence of atmospheric electric field distribution in the atmosphere. We show that the comparison of observed effects with theoretically expected will give important information on the value of atmospheric electric field and its distribution in the atmosphere. We consider also the possible influence of secondary relativistic electrons of CR and relativistic electrons precipitated from the Earth’s radiation belts on thunderstorms and lightnings, and through this – on climate change.