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.     

Ion production and ionization effect in the atmosphere during the Bastille day GLE 59 due to high energy SEPs

The influence of high energy particles, specifically cosmic rays, on atmospheric physics and chemistry is highly discussed. In most of the proposed models the role of ionization in the atmosphere due to cosmic rays is not negligible. Moreover, effect(s) on minor constituents and aerosols are recently observed, specifically over the polar regions during strong solar particle events. According to the recent findings for such effects it is necessary an essential increase of ion production, specifically during the winter period. The galactic cosmic rays are the main source of ionization in the Earth’s stratosphere and troposphere. Occasionally, the atmospheric ionization is significantly enhanced during strong solar energetic particles events, specifically over the polar caps. During the solar cycle 23 several strong ground level enhancements were observed. One of the strongest was the Bastille day event occurred on 14 July 2000. Using a full Monte Carlo 3-D model, we compute the atmospheric ionization, considering explicitly the contribution of cosmic rays with galactic and solar origin, focusing on high energy particles. The model is based on atmospheric cascade simulation with the PLANETOCOSMICS code. The ion production rate is computed as a function of the altitude above the sea level. The ion production rate is computed on a step ranging from 10 to 30?min throughout the event, considering explicitly the spectral and angular characteristics of the high energy part of solar protons as well as their time evolution. The corresponding event averaged ionization effect relative to the average due to galactic cosmic rays is computed in lower stratosphere and upper troposphere at various altitudes, namely 20?km, 15?km, 12?km and 8?km above the sea level in a sub-polar and polar regions. The 24^h and the weekly ionization effects are also computed in the troposphere and low stratosphere. Several applications are discussed.     

Periodic variations of ionization rate in the low and middle atmosphere

The cosmic ray ionization source functions which were obtained using a simplified extensive air shower model are used to calculate the eleven year cycle, seasonal and diurnal variations of ionization rate in the low and middle atmosphere. The ionization source function, as a function of the penetrating depth and the energy of cosmic ray particles, is the ionization rate per unit depth for a unit flux of incoming cosmic ray particles with certain energy.The calculation of the eleven year cycle variation of ionization rate in the low and middle atmosphere due to the modulation of galactic cosmic ray intensity by solar activity shows that the amplitude is larger at a higher magnetic latitude and is generally larger at higher altitudes. The relative amplitude of fluctuation of the ionization peak value (at altitudes near 15 km) is up to 45% in the magnetic polar region. The ionization rate, due to the seasonal variation of the atmospheric density, varies from several per cent below the ionization peak to several tens per cent above the peak. This seasonal variation of ionization rate reaches 35% at 70 km. The diurnal variation of atmospheric densities caused by atmospheric tidal oscillation can produce a diurnal variation of the ionization rate to an amplitude of several per cent at altitudes above 40 km. The diurnal oscillation is less than 1% below 35 km.     

Primary mass discrimination of high energy cosmic rays using PNN and k-NN methods

Probabilistic neural network (PNN) and k-Nearest Neighbors (k-NN) methods are widely used data classification techniques. In this paper, these two methods have been used to classify the Extensive Air Shower (EAS) data sets which were simulated using the CORSIKA code for three primary cosmic rays. The primaries are proton, oxygen and iron nuclei at energies of 100?TeV–10?PeV. This study is performed in the following of the investigations into the primary cosmic ray mass sensitive observables. We propose a new approach for measuring the mass sensitive observables of EAS in order to improve the primary mass separation. In this work, the EAS observables measurement has performed locally instead of total measurements. Also the relationships between the included number of observables in the classification methods and the prediction accuracy have been investigated. We have shown that the local measurements and inclusion of more mass sensitive observables in the classification processes can improve the classifying quality and also we have shown that muons and electrons energy density can be considered as primary mass sensitive observables in primary mass classification. Also it must be noted that this study is performed for Tehran observation level without considering the details of any certain EAS detection array.     

Analytical approach to cosmic ray ionization by nuclei with charge Z in the middle atmosphere – Distribution of galactic CR effects

The effects of galactic and solar cosmic rays (CR) in the middle atmosphere are considered in this work. A new analytical approach for CR ionization by protons and nuclei with charge Z in the lower ionosphere and middle atmosphere is developed in this paper. For this purpose the ionization losses (dE/dh) according to the Bohr–Bethe–Bloch formula for the energetic charged particles are approximated in three different energy intervals. More accurate expressions for energy decrease E(h) and electron production rate profiles q(h) are derived. The obtained formulas allow comparatively easy computer programming. The integrand in q(h) gives the possibility for application of adequate numerical methods – such as Romberg method or Gauss quadrature, for the solution of the mathematical problem. On this way the process of interaction of cosmic ray particles with the upper, middle and lower atmosphere will be described much more realistically. Computations for cosmic ray ionization in the middle atmosphere are made. The full CR composition is taken into account: protons, Helium (α-particles), light L, medium M, heavy H and very heavy VH group of nuclei.     

The extension of the INCL model for simulation of shielding in space

Radiation hazard for space missions is mainly due to cosmic ray protons, helium nuclei and light ions, whose energy spectrum is maximum around 1 GeV per nucleon but remains non-negligible for energies up to 15 GeV per nucleon. Nuclear reactions induced by high energy protons are often described by intranuclear cascade plus evaporation models. The attention is focused here on the Liège Intranuclear Cascade model (INCL), which has been shown to reproduce fairly well a great deal of experimental data for nucleon-induced reactions in the 200 MeV to 2 GeV range, when coupled with the ABLA evaporation-fission code. In order to extend the model to other conditions relevant for space radiation, three improvements of INCL are under development. They are reported on here. First, the reaction model has been extended to nucleon–nucleus reactions at incident energies up to 15 GeV, mainly by the inclusion of additional pion production channels in nucleon–nucleon collisions during the cascade. Second, a coalescence mechanism for the emission of light charged particles has been implemented recently. Finally, the model has been modified in order to accommodate light ions as projectiles. First results are shown and compared with illustrative experimental data. Implications for issues concerning radiation protection in space are discussed.     

Galactic cosmic ray composition and energy spectra.

Galactic cosmic ray nuclei represent a significant risk to long-duration spaceflight outside the magnetosphere. We review briefly existing measurements of the composition and energy spectra of heavy cosmic ray nuclei, pointing out which species and energy ranges are most critical to assessing cosmic ray risks for spaceflight. Key data sets are identified and a table of cosmic ray abundances is presented for elements from H to Ni (Z = 1 to 28). Because of the 22-year nature of the solar modulation cycle, data from the approaching 1998 solar minimum is especially important to reducing uncertainties in the cosmic ray radiation hazard. It is recommended that efforts to model this hazard take advantage of approaches that have been developed to model the astrophysical aspects of cosmic rays.     

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.     

Preliminary validation of computational procedures for a new atmospheric ionizing radiation (AIR) model.

A new computational procedure to determine particle fluxes in the Earth's atmosphere is presented. The primary cosmic ray spectrum has been modeled through an analysis of simultaneous proton and helium measurements made on high altitude balloon flights and spacecraft. An improved global fit to the data was achieved through applying a unique technique utilizing the Fokker-Plank equation with a non-linear rigidity-dependent diffusion coefficient. The propagation of primary particles through the Earth's atmosphere is calculated with a three-dimensional Monte Carlo transport program called FLUKA. Primary protons and helium nuclei (alphas) are generated within the rigidity range of 0.5 GV-20 TV uniform in cos2 theta. For a given location, primaries above the effective cutoff rigidity are transported through the atmosphere. Alpha particles are initially transported with a separate package called HEAVY to simulate fragmentation. This package interfaces with FLUKA to provide interaction starting points for each nucleon originating from a helium nucleus. Results from this calculation are presented and compared to measurements.     

太阳宇宙线在电离层D层中的电离

本文根据带电粒子对D层大气电离的理论,导出了太阳宇宙线在D层的电子产生率Q(h)的表达式,并计算了不同级别的太阳宇宙线事件、不同能谱参数下,Q(h)在极区随高度的分布。结果表明,不同级别、不同能谱的太阳宇宙线事件在极区产生的电离有显著的差别。同一级别,能谱指数γ越大,在较高的高度上电子产生率越大;能谱指数越小,在较低的高度上电子产生率越大。电子产生率的分布曲线出现明显的双峰,一个峰位于60公里左右,另一个峰位于85公里左右。前一个峰主要由太阳宇宙线中质子产生的,后一个峰主要是z≥2的重粒子成分产生的。本文所得结果明显好于Velinov等人的结果。      

Propagation of cosmic ray nuclei above 50 GeV/amu

It is now well-known that the relative abundance of secondary cosmic ray nuclei which are produced by spallation in interstellar space depends on energies. The mechanism responsible for this energy dependence are not well understood because of the scarcity of data on the cosmic ray composition at energy above 50 GeV/amu. The University of Chicago experiment on the Spacelab-2 flown by the Space Shuttle in 1985 is aimed at the measurement of the elemental composition in this unexplored energy regime. Our latest results on the energy dependence of the secondary to primary abundance ratio using the pairs of elements B,C and N,O, are discussed.     

Heavy particle fluxes in Salyut space station orbit

Results of a detailed analysis of heavy ion fluences measured on the Salyut and MIR orbital stations from 1978 to 1990 are presented. The analysis has made use of new models that describe the cosmic ray fluxes and their transition through the magnetosphere. The penetration of solar cosmic ray particles to the orbit, the increase of the (ScCr)/Fe flux ratio in the orbit, and the occurrences of 200–500 MeV/nucl heavy nuclei in <30° latitudes have been analyzed.     

The brief history of experimental research of cosmic ray variations in Yakutia

Some unknown historical facts of cosmic ray studies in the north-east of the former Soviet Union related to the Yakutsk scientific group are reported for the benefit of the international scientific community. It focuses on the founders of Yu.G. Shafer Institute of Cosmophysical Research and Aeronomy of Siberian Branch of Russian Academy of Sciences. A chronology of measurements of cosmic ray intensity variations since 1949 in Yakutia (Sakha Republic; NE Siberia) is given. In particular, for the first time the data of the first solar cosmic ray event registered at Yakutsk (GLE04), with a small ionization chamber S-2 (volume: 20 L) are presented. Moreover, the data of the large ionization chamber ASK-1 (volume: 950 L) for the 1953–2003 period useful for specialists in the field of cosmic ray variations are also shown.     

识别朝阳1号聚碳酸酯记录的宇宙线重粒子

本文介绍了根据一定REL区间所对应的粒子在固体核径迹探测器中的穿行距离与粒子原子序数的关系识别重粒子的方法,并用此法得出在我国卫星上放置的朝阳1号聚碳酸酯记录的宇宙线重粒子的原子序数和比能.      

Neutron spectra in the atmosphere from interactions of primary cosmic rays

Spectra of neutrons from interactions of primary cosmic rays in the earth's atmosphere are calculated with the Monte Carlo model fluka for various depths down to sea level. We discuss the environmental models describing the primary cosmic ray spectrum and details of the calculations. Neutron energy spectra are presented for different depths in the atmosphere and for different geographical locations. By comparing results of calculations to measurements on neutron spectra it is shown that fluka may serve as an important tool for the estimation of the radiation environment in the atmosphere.     

Preliminary results from the second flight of CREAM

Launched from McMurdo (Antarctica) in December 2005, the balloon experiment CREAM (cosmic ray energetics and mass) collected about 15 million triggers during its second flight of 28 days. Redundant charge identification, by two pixelated silicon arrays and a time resolved pulse shaping technique from a scintillator system, allowed a clear signature of the primary nuclei. The energy was measured with a tungsten/SciFi calorimeter preceded by a graphite target. Preliminary results from the analysis of the data of the second flight are presented.     

A program to measure new energetic particle nuclear interaction cross sections.

The Transport Collaboration, consisting of researchers from institutions in France, Germany, Italy and the USA, has established a program to make new measurements of nuclear interaction cross sections for heavy projectiles (Z > or = 2) in targets of liquid H2, He and heavier materials. Such cross sections directly affect calculations of galactic and solar cosmic ray transport through matter and are needed for accurate radiation hazard assessment. To date, the collaboration has obtained data using the LBL Bevalac HISS facility with 20 projectiles from 4He to 58Ni in the energy range 393-910 MeV/nucleon. Preliminary results from the analysis of these data are presented here and compared to other measurements and to cross section prediction formulae.     

超高能宇宙线电子在星际介质中的传播

本文研究能量高达1012eV以上宇宙线电子在星际介质中的传播特征, 得到了一些有趣的结果。作为一种自然的猜测, 宇宙线电子高能成份的区域性特征, 可能是导致银河系γ射线辐射的非均匀成团状分布结构的直接原因。