Be-7 nuclei produced by galactic cosmic rays and solar energetic particles in the earth''s atmosphere

Be-7 radioactive nuclei with a half-life of 53.3 days result from spallation reactions of galactic cosmic rays(GCR) and solar energetic particles (SEP) with N and O nuclei in the Earth's atmosphere. We calculate the average global production of Be-7 in the atmosphere by GCR and SEP The result indicates that an intense SEP event produces a large amount of Be-7 in the polar stratosphere and part of them could be transported to the surface at lower latitudes. The ground-level measurement of Be-7 in Japan exhibits the possibility of enhancement in the Be-7 radioactivity associated with the intense SEP event on July 14, 2000. In addition, the present experiment shows seasonal variations in the surface Be-7 concentration which peaks in spring and autumn. We discuss the possible air mass mixing between the stratosphere and troposphere to explain the measured seasonal variations. The surface concentration of Pb-210 nuclei indicates a similar trend to that of Be-7 and we suggest two possible explanations.     

Regional and temporal variability of solar activity and galactic cosmic ray effects on the lower atmosphere circulation

In this work we studied the spatial and temporal structure of long-term effects of solar activity (SA) and galactic cosmic ray (GCR) variations on the lower atmosphere circulation as well as possible reasons for the peculiarities of this structure. The study revealed a strong latitudinal and regional dependence of SA/GCR effects on pressure variations in the lower troposphere which seems to be determined by specific features of baric systems formed in different regions. The temporal structure of SA/GCR effects on the troposphere circulation at high and middle latitudes is characterized by a roughly 60-year periodicity which is apparently due to the epochs of the large-scale atmospheric circulation. It is suggested that a possible mechanism of long-term effects of solar activity and cosmic ray variations on the troposphere circulation involves changes in the evolution of the polar vortex in the stratosphere of high latitudes, as well as planetary frontal zones.     

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.     

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.     

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.     

Stratospheric polar vortex as a possible reason for temporal variations of solar activity and galactic cosmic ray effects on the lower atmosphere circulation

Possible reasons for the temporal instability of long-term effects of solar activity (SA) and galactic cosmic ray (GCR) variations on the lower atmosphere circulation were studied. It was shown that the detected earlier ∼60-year oscillations of the amplitude and sign of SA/GCR effects on the troposphere pressure at high and middle latitudes (Veretenenko and Ogurtsov, Adv.Space Res., 2012) are closely related to the state of a cyclonic vortex forming in the polar stratosphere. The intensity of the vortex was found to reveal a roughly 60-year periodicity affecting the evolution of the large-scale atmospheric circulation and the character of SA/GCR effects. An intensification of both Arctic anticyclones and mid-latitudinal cyclones associated with an increase of GCR fluxes at minima of the 11-year solar cycles is observed in the epochs of a strong polar vortex. In the epochs of a weak polar vortex SA/GCR effects on the development of baric systems at middle and high latitudes were found to change the sign. The results obtained provide evidence that the mechanism of solar activity and cosmic ray influences on the lower atmosphere circulation involves changes in the evolution of the stratospheric polar vortex.     

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.     

Estimation of the solar proton spectrum in the GLE70 event

The GLE event on December 13, 2006 as observed by network station neutron monitor data is investigated (is considered). The GLE energetic spectrum suggested for this event is estimated taking into account the primary differential spectrum of galactic cosmic rays, coupling coefficients and integral multiplicities of concrete detectors at different latitudes and levels of observation. It is noted that the additional increase of solar protons is also manifested in the ionization chamber ASK-1 data at the Yakutsk station.     

Model for induced ionization by galactic cosmic rays in the Earth atmosphere and ionosphere

The ionization profiles produced by galactic cosmic rays in the Earth atmosphere and ionosphere are obtained on the basis of Monte Carlo simulations. Cascade processes in the atmosphere are simulated using CORSIKA 6.52 code with FLUKA 2006 and QGSJET II hadronic interaction subroutines. Proton induced showers are considered using a realistic atmospheric model (US Standard Atmosphere). The energy deposit from different components is taken into account, namely electromagnetic, hadron and muon components. The curvature of the atmosphere is considered in the computer code. On the basis of the computed ionization yield function the ion pair production rate in the atmosphere is obtained for different conditions and locations. The model is applicable to the entire atmosphere, from ground level to upper atmosphere. Several applications of the obtained results are discussed. The Monte Carlo simulation model considers nuclear interactions below the altitude of 35 km. It is compared with analytical–numerical electron production rate model. The latter model which takes into account the electromagnetic interactions above altitudes of 35 km has two main regions of application: above 50 km (thin target model) and between 35 and 50 km (intermediate target model). A good agreement between the CORSIKA results and analytical–numerical model results is found above altitude of 35 km.     

The atmospheric cosmic- and solar energetic particle radiation environment at aircraft altitudes

Galactic cosmic rays interact with the solar wind, the earth's magnetic field and its atmosphere to produce hadron, lepton and photon fields at aircraft altitudes. In addition to cosmic rays, energetic particles generated by solar activity bombard the earth from time to time. These particles, while less energetic than cosmic rays, also produce radiation fields at aircraft altitudes which have qualitatively the same properties as atmospheric cosmic rays. We have used a code based on transport theory to calculate atmospheric cosmic-ray quantities and compared them with experimental data. Agreement with these data is seen to be good. We have then used this code to calculate equivalent doses to aircraft crews. We have also used the code to calculate radiation doses from several large solar energetic particle events which took place in 1989, including the very large event that occurred on September 29^th and 30^th of that year. The spectra incident on the atmosphere were determined assuming diffusive shock theory.     

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.     

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.     

The relationship between plasma effects and cosmic radiation studied with TriTel-LMP measurements during the ESEO mission

The main point of the paper is to use the simultaneous measurements of the energetic particle flux by TriTel and those of electron density by a Langmuir probe to study the question of to what extent solar electromagnetic and corpuscular radiation (galactic cosmic rays, particle precipitation from the radiation belts) are responsible for the ionization of the atmosphere. The electron density measured by the Langmuir probe is the sum of the ionization produced by the solar electromagnetic radiation and that due to the corpuscular radiation. The ionization produced by the solar electromagnetic radiation may be computed. The flux of energetic particles in an energy range may be determined by taking the difference between the threshold energy of the TriTel telescopes and the energy corresponding to the local cut-off rigidity. As the ESEO satellite will have a quasi-polar and circular orbit, the cut-off rigidity will change from low to high latitudes, thus enabling the assignment of different energy bands for the telescopes. Thus, it will be possible to determine which energy bands of particle produce ionization at different latitudes.     

大气中子实时计算模式

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

Challenges to cosmic ray modeling: From beyond the solar wind termination shock

Voyager 1 crossed the solar wind termination shock on December 16, 2004 at a distance of 94 AU from the Sun, to become the first spacecraft to explore the termination shock region and to enter the heliosheath, the final heliospheric frontier. By the end of 2006, Voyager 1 will be at ∼101 AU, with Voyager 2 at ∼81 AU and still approaching the termination shock. Both spacecraft have been observing the modulation of galactic and anomalous cosmic rays since their launch in 1977. The recent observations close to or inside the heliosheath have provided several interesting ‘surprises’ with subsequent theoretical and modeling challenges. Examples are: what does the modulation of galactic cosmic rays amount to in this region?; how do the anomalous cosmic rays get accelerated and modulated?; why are there ‘breaks’ in the power-law slopes of the spectra of accelerated particles? Several numerical models have been applied to most of these topics over the years and comprehensive global predictions have been made the past decade, thought to be based on reasonable assumptions about the termination shock and the heliosheath. Examples of these predictions and assumptions are concisely discussed within the context of the main observed features of cosmic rays in the vicinity of the termination shock, ending with a discussion of some of the issues and challenges to cosmic ray modeling in particular.     

Radiation factors in space and a system for their monitoring.

The radiation environment is of special concern when the spaceship flies in deep space. The annual fluence of the galactic cosmic rays is approximately 10(8) cm-2 and the absorbed dose of the solar cosmic rays can reach 10 Gy per event behind the shielding thickness of 3-5 g cm-2 Al. For the radiation environment monitoring it is planned to place a measuring complex on the space probes "Mars" and "Spectr" flying outside the magnetosphere. This complex is to measure: cosmic rays composition, particle flux, dose equivalent, energy and LET spectra, solar X-rays spectrum. On line data transmission by the space probes permits to obtain the radiation environment data in space.     

Cosmic ray reacceleration on the galactic wind termination shock

Reacceleration of cosmic rays produced by galactic sources on the galactic wind termination shock is considered. The problem of the cosmic ray spectrum continuity is investigated. Numeric results are presented and discussed. We found that a smooth spectral transition from the galactic cosmic rays to the cosmic rays reaccelerated at the galactic wind termination shock is difficult to produce, if the maximum energy of accelerated particles is the same throughout the surface of the termination shock. The possible solution of this problem is the non-spherical termination shock with different maximum energies at different places of the shock.     

Normalized ionization yield function for various nuclei obtained with full Monte Carlo simulations

Several recent results important for production of ion pairs in the Earth atmosphere by various primary cosmic ray nuclei are presented. The direct ionization by various primary cosmic ray nuclei is explicitly obtained. The longitudinal profile of atmospheric cascades is sensitive to the energy and mass (charge) of the primary particle. In this study different cosmic ray nuclei are considered as primaries, namely Helium, Oxygen and Iron nuclei. The cosmic ray induced ionization is obtained on the basis of CORSIKA 6.52 code simulations using FLUKA 2006 and QGSJET II hadronic interaction models. The energy of the primary particles is normalized to GeV per nucleon. In addition, the ionization yield function Y is normalized as ion pair production per nucleon. The obtained ionization yield functions Y for various primaries are compared. The presented results and their application are discussed.     

The atmospheric radiation response to solar-particle-events.

High-energy solar particles, produced in association with solar flares and coronal mass ejections, occasionally bombard the earth's atmosphere. resulting in radiation intensities additional to the background cosmic radiation. Access of these particles to the earth's vicinity during times of geomagnetic disturbances are not adequately described by using static geomagnetic field models. These solar fluxes are also often distributed non uniformly in space, so that fluxes measured by satellites obtained at great distances from the earth and which sample large volumes of space around the earth cannot be used to predict fluxes locally at the earth's surface. We present here a method which uses the ground-level neutron monitor counting rates as adjoint sources of the flux in the atmosphere immediately above them to obtain solar-particle effective dose rates as a function of position over the earth's surface. We have applied this approach to the large September 29-30, 1989 ground-level event (designated GLE 42) to obtain the magnitude and distribution of the solar-particle effective dose rate from an atypically large event. The results of these calculations clearly show the effect of the softer particle spectra associated with solar particle events, as compared with galactic cosmic rays, results in a greater sensitivity to the geomagnetic field, and, unlike cosmic rays, the near-absence of a "knee" near 60 degrees geomagnetic latitude.     

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

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