On the possibility to check the magnetosphere’s model by CR: The strong geomagnetic storm in November 2003

We report on a study of cosmic ray cutoff rigidity variations during the strong geomagnetic storm of 18–24 November 2003. We employed the most recent Tsyganenko magnetospheric model to represent the very strong disturbed Magnetosphere. We used this magnetic field for the cosmic ray trajectory calculations to determine the geomagnetic cutoff rigidity throughout this period of severe geomagnetic disturbance. We determine the cutoff rigidity changes during this period by two methods, by trajectory calculations and by the spectrographic global survey method. The values of the change of cutoff rigidities obtained by two different methods are correlated with the D_st and interplanetary magnetic field and plasma parameters and result in correlation coefficients in the range 0.63–0.84 for the various cosmic ray stations. The result of this study indicates that the most significant contributions to the cutoff rigidity changes are due to D_st variation although the influence of solar wind density and B_z and B_y components of IMF variations is significant.     

World grid of cosmic ray vertical cut-off rigidity for the last decade

Cosmic ray cut-off rigidity tables and maps over the world concerning the epochs 2010, 2015 and the current one 2020 have been constructed. These maps display the effective cut-off rigidity in every five degrees in latitude and in longitude at the altitude of 20 km above the surface of the international reference ellipsoid. The values of the geomagnetic cut-off rigidity were calculated in every 5° in latitude and in every 15° in longitude applying the well-known method of particle trajectory calculations resulted from the theory of the particle motion in the Earth's magnetic field. The applied software employed the 12th Generation of the International Geomagnetic Reference Field (IGRF 12) and trajectories were calculated at 0.01 GV intervals in order to determine the vertical cut-off rigidity for each location. Beyond the use of the calculated cut-off rigidity values as a basic reference of charged particle access to different geographical locations during quiet and/or more intense geomagnetic periods, these results can be used for a long- term forecasting of the geomagnetic conditions variations.     

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.     

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.].     

质子地磁垂直截止刚度在2004年11月7-8日两个时刻处的数值模拟

利用单粒子轨道理论, 在T04和IGRF2000磁场模式建立的磁层模型基础上, 应用四阶龙格库塔方法, 模拟计算宇宙线带电质子在地球磁层中的运动以及沿天顶方向入射到达地球磁场内的某一特定位置, 得到了距地球表面450km高度处全球质子垂直截止刚度在2004年11月7-8日中两个时刻的计算值. 根据计算得到的684个不同位置处的截止刚度值, 分析了同一时刻地磁垂直截止刚度随磁纬和磁经的变化. 与此同时, 模拟计算了相同时刻下磁纬为30°, 磁经为0°, 45°, 90°, 135°, 180°, 225°, 270°, 315°处, 地磁有效垂直截止刚度随高度的变化情况. 结果表明, 在磁纬30°处, 其质子有效垂直截止刚度随距离地心高度的变化沿向阳处向磁尾处方向, 变化越来越缓慢.      

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.     

A method for operative calculation of charged particle penetration into the LEO

We develop a simple method for calculating the effective vertical cutoff rigidity of charged particles, taking into account the Kp-index and the local time, on the basis of generalization of the results of extensive trajectory calculations for trial particles moving in the geomagnetic field. The vertical cutoff rigidities, calculated by the Tsyganenko-89 model, are presented as an International Geomagnetic Reference Field (IGRF) model calculated and thereafter corrected in accordance with the geomagnetic disturbance and local time conditions. The fits from the proposed method agree with the results of cutoff rigidity measurements carried out by satellites. The method is intended for applications using cutoff calculations, such as evaluating particle penetration of spatial boundaries, calculating magnetospheric transmissions for low-orbital spacecrafts flights and interpreting the results of orbital experiments.     

宇宙线变化特征在强地磁暴预报中的应用

利用宇宙线中子探测数据定性分析了地面宇宙线多台站之间的相互联系以及大磁暴与宇宙线之间的响应关系. 以Irkutsk和Oulu宇宙线台站为例, 运用小波去噪技术提高数据的稳定性. 结果表明, 相同世界时条件下, 两站宇宙线通量相关性在事件发生时较高; 而相同地方时条件下, 相关性则在平静期较高. 进一步采用相同地方时条件对不同宇宙线台站的通量在平静期和扰动期的相对变化进行分析, 选取2004年7月强地磁暴典型事例进行直观分析, 发现大地磁暴前Irkutsk和Oulu台站的宇宙线相对通量发生明显差异, 可以尝试作为强地磁暴宇宙线先兆特征. 通过对2001年3月至2005年5月的强磁暴和中强磁暴进行统计, 得到与强地磁暴相关的适当宇宙线相对差异阈值. 将得到的阈值对2005年9月至2011年12月所有强磁暴及中强磁暴进行验证, 总成功率达到87.5%, 误报率为35.7%, 结果较好.      

宇宙线暴时增加及其特征

利用11个地面超中子堆的资料分析了6个不同地磁暴期间由于赤道环电流增强引起的暴时增加,并分析了它们的特征。结果表明,地磁暴愈强,暴时增加愈明显。暴时增加与地磁H分量的变化有很好的负相关,并有明显的纬度效应或与地磁截止刚度有依赖关系。北京地区是暴时增加相当明显的区域。暴时增加还有明显的地方时依赖关系。在白昼一侧的增加更大,表明环电流具有不对称性的特征。      

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.     

利用地面宇宙线强度变化预报地磁暴方法初步研究

分析了Nagoya宇宙线闪烁体望远镜探测数据的变化特点,定性地探讨了CME可能引起的地面宇宙线的变化特征,通过实例证实了地面宇宙线通量的异常波动是地磁暴发生的重要先兆特征,并且将8 h内宇宙线通量与该时间段内平均通量的偏差D8(t)参数应用到宇宙线数据分析中.通过数据分析与讨论,认为D8(t)参数达到一定阈值是地磁暴的重要先兆特征,但不是充分条件,虚假信号仍占多数;D8(t)参数与太阳质子事件探测结果相结合,对于大地磁暴的预报有较好的效果.      

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.     

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.     

由环电流的宇宙线效应验证二维对称环电流模型

本文利用不同纬度超中子堆记录的宇宙线强度, 对不同磁暴类型下环电流的宇宙线效应进行了分析。加强的西向环电流使宇宙线截止刚度下降, 在某些特大的磁暴中, 中低纬字宙线截止刚度下降可达1GV以上, 从而使宇宙线强度增加。利用二维对称环电流模式, 从宇宙线的角度计算了环电流的磁暴效应。与实际观测值的对比显示, 除非特大的磁暴过程, 利用该模式所得结果与实际观测符合得较好。      

Cosmic ray exposure factors for Shuttle altitudes derived from calculated cut-off rigidities.

The allowed cosmic radiation flux accessible to an earth-orbiting spacecraft is a complex function of the satellite position and the geomagnetic cutoff characteristics at each zenith and azimuth angle at each position. We have determined cosmic ray exposure factors for the galactic cosmic ray spectrum for typical shuttle altitudes and inclinations up to 50 degrees. We have utilized d world grid of trajectory-derived cutoff rigidity calculations at 400 km altitude to determine geomagnetic transmission functions that permit a simple and direct calculation of the allowed cosmic ray spectrum to a 400 km satellite orbit. If the interplanetary cosmic ray spectrum is multiplied by the orbit-averaged geomagnetic transmission function the result is the allowed cosmic ray spectrum at the spacecraft.     

Variations of geomagnetic cutoff rigidity in the southern hemisphere close to 70°W (South-Atlantic Anomaly and Antarctic zones) in the period 1975–2010

We report the existence of rapid variations in (effective) geomagnetic cutoff rigidity (Rc) between the equatorial and Antarctic zones adjacent to the Andes Mountains, revealed by the variation rate of geomagnetic cutoff rigidity (VRc) in the period 1975–2010. Our analysis is based on empirical records and theoretical models of the variations in cosmic rays and on the structure of geomagnetic fields. These have given us a different view of variations in Rc in time and space along the 70°W meridian, where secular variations in the geomagnetic field are strongly influenced by the proximity of the South Atlantic Magnetic Anomaly (SAMA), one of the most important characteristics of the terrestrial magnetic field that affects our planet, close from the equator to the 50°S parallel and from South America to South Africa. The VRc presents rapid changes in mid-latitudes where SAMA exerts its influence despite the existence of smooth changes in the geomagnetic field. This shows that these changes occur mainly in the spatial configuration, rather than in the temporal evolution of Rc. The analysis was performed using measurements from the Chilean Network of Cosmic Rays and Geomagnetism Observatories, equipped with BF-3 and latest generation He-3 neutron monitors, Fluxgate magnetometers, geomagnetic reference field (IGRF) and Tsyganenko 2001 model (just for completeness).     

Eccentric dipole approximation of the geomagnetic field: Application to cosmic ray computations

A comparison of the full IGRF model of the geomagnetic field with two simplified models, the truncated IGRF and the eccentric dipole model, is performed. The simplified models were found to provide a reasonable approximation for the large scale geomagnetic field distribution. In the application of the simplified geomagnetic models to the shielding of cosmic rays in the magnetosphere as quantified via the geomagnetic cut-off rigidity, the eccentric dipole and the truncated IGRF provide a good large scale view. The use of the simplified model does not introduce any additional systematic errors at the global scale but may be a source of moderate uncertainty at the regional scale in the tropical Atlantic region. This study quantitatively validates the use of such simplified geomagnetic models when describing the shielding of cosmic rays in the magnetosphere.     

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.     

Large geomagnetic storms of extreme solar event periods in solar cycle 23

During extreme solar events such as big flares or/and energetic coronal mass ejections (CMEs) high energy particles are accelerated by the shocks formed in front of fast interplanetary coronal mass ejections (ICMEs). The ICMEs (and their sheaths) also give rise to large geomagnetic storms which have significant effects on the Earth’s environment and human life. Around 14 solar cosmic ray ground level enhancement (GLE) events in solar cycle 23 we examined the cosmic ray variation, solar wind speed, ions density, interplanetary magnetic field, and geomagnetic disturbance storm time index (D_st). We found that all but one of GLEs are always followed by a geomagnetic storm with D_st  −50 nT within 1–5 days later. Most(10/14) geomagnetic storms have D_st index  −100  nT therefore generally belong to strong geomagnetic storms. This suggests that GLE event prediction of geomagnetic storms is 93% for moderate storms and 71% for large storms when geomagnetic storms preceded by GLEs. All D_st depressions are associated with cosmic ray decreases which occur nearly simultaneously with geomagnetic storms. We also investigated the interplanetary plasma features. Most geomagnetic storm correspond significant periods of southward B_z and in close to 80% of the cases that the B_z was first northward then turning southward after storm sudden commencement (SSC). Plasma flow speed, ion number density and interplanetary plasma temperature near 1 AU also have a peak at interplanetary shock arrival. Solar cause and energetic particle signatures of large geomagnetic storms and a possible prediction scheme are discussed.     

Cosmic ray variations of solar origin in relation to human physiological state during the December 2006 solar extreme events

There is an increasing amount of evidence linking biological effects to solar and geomagnetic disturbances. A series of studies is published referring to the changes in human physiological responses at different levels of geomagnetic activity. In this study, the possible relation between the daily variations of cosmic ray intensity, measured by the Neutron Monitor at the Cosmic Ray Station of the University of Athens (http://cosray.phys.uoa.gr) and the average daily and hourly heart rate variations of persons, with no symptoms or hospital admission, monitored by Holter electrocardiogram, is considered. This work refers to a group of persons admitted to the cardiological clinic of the KAT Hospital in Athens during the time period from 4th to 24th December 2006 that is characterized by extreme solar and geomagnetic activity. A series of Forbush decreases started on 6th December and lasted until the end of the month and a great solar proton event causing a Ground Level Enhancement (GLE) of the cosmic ray intensity on 13th December occurred. A sudden decrease of the cosmic ray intensity on 15th December, when a geomagnetic storm was registered, was also recorded in Athens Neutron Monitor station (cut-off rigidity 8.53 GV) with amplitude of 4%. It is noticed that during geomagnetically quiet days the heart rate and the cosmic ray intensity variations are positively correlated. When intense cosmic ray variations, like Forbush decreases and relativistic proton events produced by strong solar phenomena occur, cosmic ray intensity and heart rate get minimum values and their variations, also, coincide. During these events the correlation coefficient of these two parameters changes and follows the behavior of the cosmic ray intensity variations. This is only a small part of an extended investigation, which has begun using data from the year 2002 and is still in progress.