New index of long-term variations of galactic cosmic ray intensity

We find that the soft rigidity spectrum of the Galactic Cosmic Ray (GCR) intensity variations for the maximum epoch and the hard rigidity spectrum for the minimum epoch calculated based on the neutron monitors experimental data (1960–2002) are related with the various dependence of the diffusion coefficient on the GCR particle’s rigidity for different epoch of solar activity. This dependence is stronger in the maximum epoch than in the minimum epoch of solar activity, and is provided by the essential temporal rearrangements of the structure of the Interplanetary Magnetic Field (IMF) turbulence from the maxima to minima epoch of solar activity. We also show that the rigidity spectrum of GCR intensity variations is harder for the effective rigidities ∼(10–15) GV (by neutron monitors data), than for the effective rigidities ∼(25–30) GV (by neutron monitors and muon telescopes data). A general scenario of GCR modulation versus solar activity is settled on the essential temporal rearrangements of the structure of the IMF turbulence. Therefore, the temporal changes of the power law rigidity spectrum exponent can be considered as a vital (new) index to explain the 11-year variations of the GCR intensity. We assume that ∼(70–80)% of the changes of the amplitudes of the 11-year variations of GCR intensity is related with the changes of the IMF turbulence versus solar activity.     

Effects of geomagnetic secular variations on cosmic ray access to the terrestrial environment

Taking advantage of the cutoff computations performed for more than a hundred locations from 1955 to 1995 [every 5 years; Shea, M.A., Smart, D.F. Vertical cutoff rigidities for cosmic ray stations since 1955, in: Proceedings of the ICRC 2001, Hamburg, Copernicus Gesellschaft, vol. 10, pp. 4063–4066, 2001], we carefully checked the relationship between the vertical cutoff rigidity and the McIlwain parameter introduced by Shea et al. [Shea, M.A., Smart, D.F., Gentile, L.C. Estimating cosmic ray vertical cutoff rigidities as a function of the McIlwain L-parameter for different epochs of the geomagnetic field. Phys. Earth Planet. Int., 48, 200–205, 1987]. We derived an updated algorithm that can be used outside the polar and equatorial regions, avoiding time consuming computations. Results for the European area and 1990 epoch suggest that the fast evaluation is accurate within 0.1 GV in 26 out of the 30 considered locations.     

The cosmic ray ground level enhancement of 6 November 1997.

The relativistic solar proton event of 6 November 1997 resulted in the first ground-level enhancement (GLE) of solar cycle 23. The earliest onset was around 1215 UT but was up to 15 minutes later at some neutron monitor locations. The time of maximum intensity also varied significantly over the world-wide neutron monitor network. The modeled particle distributions and spectra are presented. The apparent particle arrival direction is found to be largely consistent with propagation outward from the sun along interplanetary magnetic field lines.     

Evaluation of the dynamic cutoff rigidity model using dosimetry data from the STS-28 flight.

We have a developed a dynamic cutoff rigidity model based on computed world grids of vertical cutoff rigidities derived from employing the Tsyganenko magnetospheric model. The dynamic range of this model covers all magnetic activity levels specified by integer values of the Kp magnetic index. We present comparisons of the measured dose observed on the space shuttle during the August 1989 solar proton event with the dose computed from solar particles predicted to be allowed through the magnetosphere to the space shuttle position. We find a one-to-one correspondence between the portion of the orbit predicted to be subjected to solar protons and the portion of the orbit where solar particle dose measurements were obtained.     

The space-developed dynamic vertical cutoff rigidity model and its applicability to aircraft radiation dose.

We have developed a dynamic geomagnetic vertical cutoff rigidity model that predicts the energetic charged particle transmission through the magnetosphere. Initially developed for space applications, we demonstrate the applicability of this library of cutoff rigidity models for computing aircraft radiation dose. The world grids of vertical cutoff rigidities were obtained by particle trajectory tracing in a magnetospheric model. This reference set of world grids of vertical cutoff rigidities calculated for satellite altitudes covers all magnetic activity levels from super quiet to extremely disturbed (i.e., Kp indices ranging from 0 to 9+) for every three hours in universal time. We utilize the McIlwain "L" parameter as the basis of the interpolation technique to reduce these initial satellite altitude vertical cutoff rigidities to cutoff rigidity values at aircraft altitudes.     

Cosmic radiation exposure of aircraft occupants on simulated high-latitude flights during solar proton events from 1 January 1986 through 1 January 2008

From 1 January 1986 through 1 January 2008, GOES satellites recorded 170 solar proton events. For 169 of these events, we estimated effective and equivalent dose rates and doses of galactic cosmic radiation (GCR) and solar cosmic radiation (SCR), received by aircraft occupants on simulated high-latitude flights. Dose rate and dose estimates that follow are for altitudes 30, 40, 50, and 60 kft, in that order.     

Variations in cosmic ray cutoff rigidities during the great geomagnetic storm of November 2004

A very strong interplanetary and magnetospheric disturbance observed on 7–13 November 2004 can be regarded as one of the strongest events during the entire period of space observations. In this paper we report on the studies of cosmic ray cutoff rigidity variations during 7–13 November 2004 showing how storm conditions can affect the direct cosmic ray access to the inner magnetosphere. Effective cutoff rigidities have been calculated for selected points on the ground by tracing trajectories of cosmic ray particles through the magnetospheric magnetic field of the “storm-oriented” Tsyganenko 2003 model. Cutoff rigidity variations have also been determined by the spectrographic global survey method on the basis of experimental data of the neutron monitor network. Relations between the calculated and experimental cutoff rigidities and the geomagnetic Dst-index and interplanetary parameters have been investigated. Correlation coefficients between the cutoff rigidities obtained by the trajectory tracing method and the spectrographic global survey method have been found to be in the limits 0.76–0.89 for all stations except the low-latitude station Tokyo (0.35). The most pronounced correlation has been revealed between the cutoff rigidities that exhibited a very large variation of ∼1–1.5 GV during the magnetic storm and the Dst index.     

南极中山站宇宙线观测

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

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.     

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.     

The IGY and beyond: A brief history of ground-based cosmic-ray detectors

The state of art of ground-based cosmic-ray research from its discovery to present is reviewed. After discovery of cosmic rays by Hess in 1912, the nature of the primary and secondary radiation was established from recordings by a variety of instruments, sensitive to various components of cosmic rays and operated at different latitudes, longitudes and altitudes, including instruments carried by balloons. The IGY formalized international co-operation and coordinated study of cosmic rays, which is vital for meaningful interpretation of cosmic-ray data. Data collected at different geographic locations require an effective cutoff rigidity as a data ordering parameter. This parameter is obtained from tracing trajectories of primary cosmic rays in the Earth’s magnetic field. After 50 years the world’s neutron monitor network remains still the backbone for studying intensity variations of primary cosmic rays in the rigidity ranges between 1 and 15 GV, associated with transport and with transient events. Also the penetrating muon and neutrino components of secondary cosmic rays have a long history of recording and fundamental problem investigations. Valuable data about composition and spectrum of primary cosmic rays in ever increasing high-energy regions have been obtained during the years of investigations with various configurations and types of extensive air shower detectors. The culture of personal involvement of the physicist in carrying out experiments and data acquisition characterized the continued vitality of cosmic-ray investigations ranging from its atmospheric, geomagnetic and heliospheric transport through to its solar and astrophysical origins.     

Effective non-vertical and apparent cutoff rigidities for a cosmic ray latitude survey from Antarctica to Italy in minimum of solar activity

In this paper we will report the results of the computation of cutoff rigidities of vertical and non-vertical incident cosmic ray particles. Non-vertical effective cutoff rigidities have been computed by tracing particle trajectories through the “real” geomagnetic magnetic field comprising the International Geomagnetic Reference Field model (IGRF95, IAGA Division 5 Working Group 8, 1996: Sabaka, T.J., Langel, R.A., Baldwin, R.T., Conrad, J.A. The geomagnetic field, 1900–1995, including the large scale fields from magnetospheric sources and NASA candidate models for the 1995 IGRF revision. J. Geomag. Geoelect. 49, 157–206, 1997.) and the Tsyganenko [Tsyganenko, N.A. A magnetospheric magnetic field model with a warped tail current sheet. Planet. Space Sci. 37, 5–20, 1989.] magnetosphere model. The computation have been done for the backward route (from Antarctica to Italy) of the Italian Antarctic ship survey 1996–1997, for geographic points corresponding to the daily average coordinates of the ship; for zenith angles 15°, 30°, 45° and 60°, and azimuth angles from 0° to 360° in steps of 45°. By means of the obtained non-vertical cutoffs the apparent cutoff rigidities have been calculated. The information on integral multiplicities of secondary neutrons detected by the neutron monitor in dependence of the zenith angle of incoming primary cosmic ray particles have also been used. This information is based on the theoretical calculations of meson-nuclear cascades of primary protons with different rigidities arriving to the Earth’s atmosphere at the zenith angles of 0°, 15°, 30°, 45°, 60° and 75°. The difference between the computed apparent and vertical cutoff rigidities reaches ∼1 GV at rigidities >7–8 GV. At rigidities of 10–16 GV, the difference between the apparent and vertical cutoff rigidities is larger than that obtained earlier by Clem et al. [Clem, J.M., Bieber, J.W., Duldig, M., Evenson, P., Hall, D., Humble, J.E. Contribution of obliquely incident particles to neutron monitor counting rate. J. Geophys. Res. 102, 26919–26926, 1997.] and Dorman et al. [Dorman, L.I., Villoresi, G., Iucci, N., Parisi, M., Tyasto, M.I., Danilova, O.A., Ptitsyna, N.G. Cosmic ray survey to Antarctica and coupling functions for neutron component near solar minimum (1996–1997), 3. Geomagnetic effects and coupling functions. J. Geophys. Res. 105, 21047–21056, 2000.].     

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.     

1989年9月29日GLE事件的特征

1989年9月29日发生了30多年来最大的地面宇宙线增强事件(即GLE事件)。本文描述这次世界范围内的大事件的主要特征,包括强度时间变化、能谱、各向异性以及地理位置上的分布特征,并简述了这次事件对地球环境产生的影响。      

The limitations of using vertical cutoff rigidities determined from the IGRF magnetic field models for computing aircraft radiation dose.

Vertical cutoff rigidities derived from the International Geomagnetic Reference Fields (IGRF) are normally used to compute the radiation dose at a specific location and to organize the radiation dose measurements acquired at aircraft altitudes. This paper presents some of the usually ignored limits on the accuracy of the vertical cutoff rigidity models and describes some of the computational artifacts present in these models. It is noted that recent aircraft surveys of the radiation dose experienced along specific flight paths is sufficiently precise that the secular variation of the geomagnetic field is observable.     

Geomagnetic cutoffs: a review for space dosimetry applications.

The earth's magnetic field acts as a shield against charged particle radiation from interplanetary space, technically described as the geomagnetic cutoff. The cutoff rigidity problem (except for the dipole special case) has "no solution in closed form". The dipole case yields the Stormer equation which has been repeatedly applied to the earth in hopes of providing useful approximations of cutoff rigidities. Unfortunately the earth's magnetic field has significant deviations from dipole geometry, and the Stormer cutoffs are not adequate for most applications. By application of massive digital computer power it is possible to determine realistic geomagnetic cutoffs derived from high order simulation of the geomagnetic field. Using this technique, "world-grids" of directional cutoffs for the earth's surface and for a limited number of satellite altitudes have been derived. However, this approach is so expensive and time consuming it is impractical for most spacecraft orbits, and approximations must be used. The world grids of cutoff rigidities are extensively used as lookup tables, normalization points and interpolation aids to estimate the effective geomagnetic cutoff rigidity of a specific location in space. We review the various options for estimating the cutoff rigidity for earth-orbiting satellites.     

Effect of “Noisy” sun conditions on aircrew radiation exposure

In computer codes used to estimate the aircrew radiation exposure from galactic cosmic radiation, a quiet sun model is usually assumed. A revised computer code (PCAIRE ver. 8.0f) is used to calculate the impact of noisy sun conditions on aircrew radiation exposure. The revised code incorporates the effect of solar storm activity, which can perturb the geomagnetic field lines, altering cutoff rigidities and hence the shielding capability of the Earth’s magnetic field. The effect of typical solar storm conditions on aircrew radiation exposure is shown to be minimal justifying the usual assumptions.     

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.     

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.