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.     

Real-time GLE alert in the ANMODAP Center for December 13, 2006

Within the last years, a real-time system to monitor high energy cosmic rays for space weather use has been operated at Athens cosmic ray station. Neutron monitors and satellite high resolution data in real time are used, making it possible to observe cosmic rays in dual energy range observations. In large solar energetic particle (SEP) events, ground level enhancement (GLE) can provide the earliest alert for the onset of the SEP event. This system watches for count rate increases recorded in real time by 23 neutron monitors, which triggers an alarm if a ground level enhancement (GLE) of cosmic ray intensity is detected.     

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.     

On transmissivity of low energy cosmic rays in disturbed magnetosphere

Access of low energy cosmic rays to any position on the Earth depends on the state of the magnetosphere. Anisotropy of cosmic rays, deduced from the neutron monitor network, must assume the variable transmissivity of the magnetosphere especially during the geomagnetic disturbances. We illustrate that computations based on different available models of geomagnetic field during selected strong geomagnetic disturbances in 2003 and 2004 imply different profiles of cut-off rigidities in time, different transmissivity functions and different asymptotic directions. Using of cosmic ray records by neutron monitors at middle and low latitudes during geomagnetically active periods, in addition to cosmic ray anisotropy in interplanetary space deduced from high and low energy cosmic ray ground based measurements, may be used for checking validity of geomagnetic field models.     

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.     

A study of the ground level enhancement of 23 February 1956

One of the greatest and most famous increase of solar cosmic rays over the neutron monitor epoch is the ground level enhancement in 1956. All future proton events are inevitable when compared with this one and therefore it is necessary to provide the efficiency of such a comparison derived from the existing data. In this paper, we return to the analysis of ground level observations on 23 February 1956 in order to model more precisely the solar cosmic ray behaviour. The extremely high magnitude of this effect allowed various spectral characteristics of solar cosmic rays, their anisotropy, differential and integral proton fluxes, and angular distribution of the source of solar particle anisotropy to be obtained with sufficient accuracy on the basis of available data from 13 neutron monitors. The most outstanding feature of this event was a narrow and extremely intensive beam of ultra relativistic particles arriving at Earth at the beginning of the event. This unique beam was not long and its width did not exceed 30–40°, thus, its contribution to solar particle density was not significant. Many features of this GLE are apparently explained by the peculiarity of particle interplanetary propagation from a remote (limb or behind of limb) source.     

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.     

Interactive database of cosmic ray anisotropy (DB-A10)

The worldwide neutron monitor network is a unique tool for obtaining with high accuracy the information on density variations, energy spectrum and anisotropy of comic rays at the Earth, outside its atmosphere and magnetosphere. These hourly averaged parameters were obtained over the whole period of cosmic ray monitoring by the ground level neutron monitor network (from 1957 till present) and are collected within the MySQL database. The Internet-project has developed for free access and supplying of cosmic ray density and anisotropy data in different formats.     

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.     

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.     

大气中子实时计算模式

为实时评估0~100km高度范围内的大气中子全球分布,对宇宙线在地磁场和大气中的传输过程进行了分析.利用蒙特卡罗方法工具包Geant4,预先计算不同能量的粒子在大气层中产生的次级粒子能谱分布,形成大气次级粒子数据库,并与相关模型进行对比,验证了该数据库的有效性和可靠性.以实测或预报的空间环境参数作为输入,计算同步轨道银河宇宙线和太阳质子事件能谱以及100km高度上的地磁垂直截止刚度,最终得到大气层顶上的粒子能谱.通过对大气次级粒子数据库的线性插值,实现1h分辨率的大气中子能谱和辐射剂量全球分布的实时计算.      

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.     

Empirical model of cosmic ray spectrum in energy interval 1 MeV–100 GeV during 11-year solar cycle

The galactic cosmic rays (GCR) are the main ionization source at altitude of ∼3–35 km in the atmosphere. For high latitude anomalous cosmic ray (ACR) component has also a significant influence on the atmospheric ionization. We propose an empirical model for differential spectra D(E) of galactic and anomalous cosmic rays in energy interval 1 MeV–100 GeV during solar cycle. In the model data are used which cover three solar cycles: 20, 22 and 23. The LEAP87, IMAX92, CAPRICE94, AMS98 and BESS experimental spectra for protons and alpha particles are fitted to the proposed empirical model. The modulated GCR differential spectra are compared with force-field approximation to the one-dimensional transport equation and with solutions of two-dimensional cosmic ray transport equation. For experimental spectra, the calculation of the model parameters is performed by Levenberg–Marquardt algorithm, applied to the special case of least squares. Algorithm that combines the rapid local convergence of Newton–Raphson method with globally convergent method for non-linear systems of equations is applied for theoretically obtained differential spectra. The described programmes are realized in algorithmic language C++. The proposed model gives practical possibility for investigation of experimental data from measurements of galactic cosmic rays and their anomalous component.     

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.     

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.     

Computation of radiation environment during ground level enhancements 65, 69 and 70 at equatorial region and flight altitudes

The radiation environment in the troposphere of the Earth is governed by cosmic rays of galactic and solar origin. During major solar energetic particles events the radiation environment changes dramatically. As a results the risk of biological effects due to exposure to ionizing radiation of aircrew increases. Here we present a numerical model for computation of absorbed dose in air due to cosmic rays of galactic and solar origin. It is applied for computation of radiation environment at flight altitude in the equatorial region during several major ground level enhancements, namely GLE65 on 28 October 2003, GLE69 on 20 January 2005 and GLE70 on 13 December 2006. The model is based on a full Monte Carlo simulation of cosmic ray induced atmospheric cascade. The cascade simulation is carried out with CORSIKA 6.990 code with corresponding hadron generators FLUKA 2011 and QGSJET II. The contribution of different cascade components, namely electromagnetic, hadron and muon is explicitly obtained. The spectra of arriving solar energetic particles are calculated from ground level measurements with neutron monitors and satellite data from GOES. The obtained results are discussed.     

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.     

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.     

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.