Calculation of the barometric coefficients at the start of the 24th solar activity cycle for particle detectors of Aragats Space Environmental Center

After the major modernization of the data acquisition electronics of the particle detectors operated at Aragats Space Environmental Center (ASEC) calculations of the barometric coefficients of all the monitors were performed in the beginning of the 24th solar activity cycle. The barometric coefficients of particle detectors located at altitudes of 1000 m, 2000 m and 3200 m a.s.l. measuring various secondary cosmic ray fluxes were compared with theoretical expectations and monitors operated on different longitudes and latitudes. The barometric coefficients were also calculated for the several neutron monitors of recently established Eurasian database (NMDB) and SEVAN particle detector networks. The latitude and altitude dependencies of the barometric coefficients were investigated, as well as the dependence of coefficients on energy of the primary particles.     

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.     

The solar cosmic ray ground-level enhancements on 20 January 2005 and 13 December 2006

Close to the current solar activity minimum, two large solar cosmic ray ground-level enhancements (GLE) were recorded by the worldwide network of neutron monitors (NM). The enormous GLE on 20 January 2005 is the largest increase observed since the famous GLE in 1956, and the solar cosmic-ray event recorded on 13 December 2006 is among the largest in solar cycle 23. From the recordings of the NMs during the two GLEs, we determined the characteristics of the solar particle flux near Earth.     

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.     

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.     

Statistical analysis of the Thunderstorm Ground Enhancements (TGEs) detected on Mt. Aragats

Starting from 2008 experimental facilities of the Aragats Space Environmental Center (ASEC) routinely measure time series of secondary cosmic ray fluxes. At these years of the minimum of solar activity we analyze the new high-energy phenomena in the terrestrial atmosphere. Namely, Thunderstorm Ground Enhancements (TGEs) and Extensive Cloud Showers (ECSs). Several new particle detectors were designed and fabricated having lower energy threshold to detect particle fluxes from the thunderclouds; some of them have possibility to distinguish charged and neutral fluxes. During 2008–2012 years ASEC detectors located at Aragats, Nor Amberd and Yerevan were detected ∼300 TGE enhancements. Amplitude of majority of them is less than 5%; however, 13 TGEs have amplitude exceeding 20%. The maximal value of observed enhancement was 271% (September 19, 2009). The paper summarizes five-years study of the TGEs on Aragats. The statistical analysis revealing the month and day-of-time distributions of TGE events, as well as the amplitude and event duration diagrams are presented.     

Statistical study of the detection of solar protons of highest energies at 20 January 2005

On January 20, 2005, 7:02–7:04 UT the Aragats Multichannel Muon Monitor (AMMM) registered enhancement of the high energy secondary muon flux (energy threshold ∼5 GeV). The enhancement, lasting 3 min, has statistical significance of ∼4σ and is related to the X7.1 flare seen by the GOES satellite and the ground level enhancement detected by the world-wide network of neutron monitors and by muon detectors. The most probable proton energy corresponding to the measured 5 GeV muon flux is within 23–30 GeV. Due to upmost importance of the detection of solar particles of highest energies in presented paper we perform detailed statistical analysis of the detected peak. The statistical technique introduced in the paper is also appropriate for the searches of sources of ultra-high energy cosmic rays.     

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.     

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.     

Is there an instrumental drift in the counting rate of some high latitude neutron monitors?

For the last six decades the neutron monitors have provided a continuous string of very reliable data to the heliophysics community. Although neutron monitors are not the primary source of data for the galactic cosmic rays, these data serve as a baseline reference for the data collected by the detectors on board the satellites and deep space probes, far away from earth orbit. The pressure corrected hourly data are available from the World Data Centers. These data have been used to derive deep insights pertaining to the electromagnetic states of the heliosphere and the modes of transport of energetic charged particles in the tangled interplanetary magnetic fields. We present evidence that some of the high latitude neutron monitors are undergoing long-term drifts in their baselines. In particular, we argue that there is no physical basis to justify the observed long-term downward trend in the baseline of the South Pole neutron monitor. The real reason may have to do with its maintenance at a distant location with challenging logistics and an improper normalization of its data after the 26 months break in the 1970s.     

On the production of highest energy solar protons at 20 January 2005

On January 20, 2005, 7:02–7:05 UT the Aragats Multidirectional Muon Monitor (AMMM) located at 3200 m a.s.l. registered enhancement of the high energy secondary muon flux (threshold ∼5 GeV). The enhancement, lasting for 3 min, has statistical significance of ∼4σ and is related to the X7.1 flare seen by the GOES, and very fast (>2500 km/s) CME seen by SOHO, and the Ground Level Enhancements (GLE) #69 detected by the world-wide network of neutron monitors and muon detectors. The energetic and temporal characteristics of the muon signal from the AMMM are compared with the characteristics of other monitors located at the Aragats Space-Environmental Center (ASEC) and with other neutron and muon detectors. Since secondary muons with energies >5 GeV are corresponding to solar proton primaries with energies 20–30 GeV we conclude that in the episode of the particle acceleration at 7:02–7:05 UT 20 January 2005 solar protons were accelerated up to energies in excess of 20 GeV.     

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.     

Solar proton events recorded in the stratosphere during cosmic ray balloon observations in 1957–2008

Long-term balloon observations have been performed by the Lebedev Physical Institute since 1957 up to the present time. The observations are taken several times a week at the polar and mid latitudes and allow us to study dynamics of galactic and solar cosmic ray as well as secondary particle fluxes in the atmosphere and in the near-Earth space. Solar energetic particles (120) – mostly protons – (SEP) events with >100 MeV proton intensity above 1 cm⁻² s⁻¹ s⁻¹ were recorded during 1958–2006. Before the advent of the SEP monitoring on spacecraft these results constituted the only homogeneous series of >100 MeV SEP events. The SEP intensities and energy spectra inferred from the Lebedev Physical Institute observations are consistent with the results taken in the adjacent energy intervals by the spacecraft and neutron monitors. Joint consideration of the SEP events series recorded by balloons and by neutron monitors during solar cycles 20–23 makes it possible to restore the probable number of events in solar cycle 19, which was not properly covered by observations. Some correlation was found between duration of SEP event production in a solar cycle and sunspot cycle 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.     

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.     

Lunar cosmic ray radiation environments during Luna and Lunar Reconnaissance Orbiter missions

The RV-2N-series instruments onboard Luna missions and the Cosmic Ray Telescope for the Effects of Radiation (CRaTER) instrument onboard Lunar Reconnaissance Orbiter (LRO) were designed to characterize the global lunar radiation environment and its biological impacts by measuring cosmic ray (CR) intensity. In this study, we have shown that the RV-2N-series instruments onboard of Russian Luna missions and the CRaTER reliably detect both background CRs and solar proton events (SPEs) in the lunar radiation environment using the proton intensity measured by the RV-2N-series onboard Luna missions out of the Russian Luna program for the exploration of the Moon (November 1970–August 1975) and the CR intensity on the Moon observed by the CRaTER (June 2009–March 2011). Those were compared with the CR intensities observed by neutron monitors (McMurdo, Thule, Oulu) on the Earth. The sunspot number is used as the index of solar activity (NOAA National Geophysical Data Center). As a result, the background CR intensities on the Moon turned out to have a good anti-correlation with the solar activity. We have also identified the proton intensity increasing events on the Moon which have the similar profiles to those observed by neutron monitors on the Earth. Most of these events show the significant increase of proton intensities in the lunar radiation environment when the SPEs associated with solar eruptions are verified. Therefore, most of the proton intensity increasing events are associated with the energetic solar particles in the lunar environment.     

Median filtering algorithms for multichannel detectors

Particle detectors of worldwide networks are continuously measuring various secondary particle fluxes incident on Earth surface. At the Aragats Space Environmental Center (ASEC), the data of 12 cosmic ray particle detectors with a total of ∼280 measuring channels (count rates of electrons, muons and neutrons channels) are sent each minute via wireless bridges to a MySQL database. These time series are used for the different tasks of off-line physical analysis and for online forewarning services. Usually long time series contain several types of errors (gaps due to failures of high or low voltage power supply, spurious spikes due to radio interferences, abrupt changes of mean values of several channels or/and slowly trends in mean values due to aging of electronics components, etc.). To avoid erroneous physical inference and false alarms of alerting systems we introduce offline and online filters to “purify” multiple time-series. In the presented paper we classify possible mistakes in time series and introduce median filtering algorithms for online and off-line “purification” of multiple time-series.     

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.     

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.