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

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

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

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.     

南极中山站宇宙线观测

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

强磁场扰动对宇宙线调制的统计研究

对1978─1982太阳活动高年时发生的激波、强磁场扰动及激波与强磁场扰动共存这三类事件引起的宇宙线变化进行了统计研究,得到如下结果:(1)激波与强磁场扰动共存时引起的宇宙线强度下降最为显着;只有激波或强磁场扰动时,宇宙线的强度变化相对较小;(2)标志速度间断的激波是产生宇宙线Forbush下降的重要因素;(3)速度间断在强磁场扰动对宇宙线的调制中可能起一个触发的作用。      

Prediction of expected global climate change by forecasting of galactic cosmic ray intensity time variation in near future based on solar magnetic field data

A method of prediction of expected part of global climate change caused by cosmic ray (CR) by forecasting of galactic cosmic ray intensity time variation in near future based on solar activity data prediction and determined parameters of convection-diffusion and drift mechanisms is presented. This gave possibility to make prediction of expected part of global climate change, caused by long-term cosmic ray intensity variation. In this paper, we use the model of cosmic ray modulation in the Heliosphere, which considers a relation between long-term cosmic ray variations with parameters of the solar magnetic field. The later now can be predicted with good accuracy. By using this prediction, the expected cosmic ray variations in the near Earth space also can be estimated with a good accuracy. It is shown that there are two possibilities: (1) to predict cosmic ray intensity for 1–6 months by using a delay of long-term cosmic ray variations relatively to effects of the solar activity and (2) to predict cosmic ray intensity for the next solar cycle. For the second case, the prediction of the global solar magnetic field characteristics is crucial. For both cases, reliable long-term cosmic ray and solar activity data as well as solar magnetic field are necessary. For solar magnetic field, we used results of two magnetographs (from Stanford and Kitt Peak Observatories). The obtained forecasting of long-term cosmic ray intensity variation we use for estimation of the part of global climate change caused by cosmic ray intensity changing (influenced on global cloudiness covering).     

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 local time dependence of the anisotropic solar cosmic ray flux.

The distribution of the solar cosmic radiation flux over the earth is not uniform, but the result of complex phenomena involving the interplanetary magnetic field, the geomagnetic field and latitude and longitude of locations on the earth. The latitude effect relates to the geomagnetic shield; the longitude effect relates to local time. For anisotropic solar cosmic ray events the maximum particle flux is always along the interplanetary magnetic field direction, sometimes called the Archimedean spiral path from the sun to the earth. During anisotropic solar cosmic ray event, the locations on the earth viewing "sunward" into the interplanetary magnetic field direction will observe the largest flux (when adjustments are made for the magnetic latitude effect). To relate this phenomena to aircraft routes, for anisotropic solar cosmic ray events that occur during "normal quiescent" conditions, the maximum solar cosmic ray flux (and corresponding solar particle radiation dose) will be observed in the dawn quadrant, ideally at about 06 hours local time.     

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

Estimation of long-term cosmic ray intensity variation in near future and prediction of their contribution in expected global climate change

On the basis of results obtained in our paper [Dorman, L.I. Long-term cosmic ray intensity variation and part of global climate change, controlled by solar activity through cosmic rays, Paper D2.1/C2.2/E3.1-0097-04. Adv. Space Res., 2004 (accepted)], we determine: the dimension of the Heliosphere (modulation region), radial diffusion coefficient and other parameters of convection–diffusion; drift mechanisms of long-term variations of cosmic ray (CR) dependence on particle energy; level of solar activity (SA); and generally, the solar magnetic field. We obtain this important information on the basis of CR and SA data in the past, taking into account the theory of convection–diffusion and global drift modulation of galactic CR in the Heliosphere. By using these results and other regularly published predictions of expected SA variation in the near future, as well as predictions of the next SA cycle, we may make predictions of long-term cosmic ray intensity variation expected in the near future (up to 10–12 years). In [Dorman, L.I. Long-term cosmic ray intensity variation and part of global climate change, controlled by solar activity through cosmic rays, Paper D2.1/C2.2/E3.1-0097-04. Adv. Space Res., 2004 (accepted)], properties of connections between long-term variation in CR intensity and some part of a global climate change were estimated, controlled by solar activity through CR. We show that in this way we may make predictions of some part of a global climate change expected in the near future (up to 10–12 years and maybe more, depending upon the period during which definite predictions of SA can be made), controlled by solar activity through CR. In this case, estimations of expected long-term changes in the planetary distribution of cutoff rigidities, which also influence CR intensity, as well as CR-influenced effects on global climate variation, become important.     

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.     

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.     

Peculiarities of galactic cosmic ray Forbush effects during October–November 2003

Features of two successive Forbush effects of the galactic cosmic ray intensity in October–November 2003 have been studied based on the neutron monitors data. The rigidity spectrum of the galactic cosmic ray intensity in the course of the first Forbush effect (22–27 October) is gradually hardening, while the rigidity spectrum of the second Forbush effect (28 October–10 November) from the starting moment is very hard. As far, the energy range of the turbulence of the interplanetary magnetic field is in general responsible for the diffusion of galactic cosmic ray particles of the energy 5–50 GeV (to which neutron monitors are sensitive), we postulate that the gradually hardening (from day to day) of the rigidity spectrum of the first Forbush effect is associated with the enhancement of the power spectral density in the energy range of the interplanetary magnetic field turbulence caused by the large scale irregularities generated due to the interaction of the extending high speed disturbances with the background solar wind. The very hard rigidity spectrum (from the starting moment) of the second Forbush effect is generally associated with the well established new structure of the energy range of the interplanetary magnetic field turbulence enriched by the already created large scale irregularities. The gradually softening of the rigidity spectrum during the recovery phase of the second Forbush effect confirms that the disturbed interplanetary magnetic field turbulence step by step returns to the initial state.     

Data reduction and correction algorithm for digital real-time processing of cosmic ray measurements: NM64 monitoring at Dourbes

We present a data correction algorithm for real-time data processing for the NM64 galactic cosmic ray neutron monitor at the Royal Meteorological Institute (RMI) in Dourbes, Belgium. The correction is based on three main tests: a continuity test, tube ratios test and a derivative test. The continuity test works as a high pass filter with a threshold based on the entire recorded dataset. Additionally, it monitors whether the logging takes place at regular intervals (continuously). The ratios test identifies noisy sections and the final derivative test criterion will identify single or double spikes by testing them against the median increase of the intensities. Using these criteria, all data from the cosmic ray station at Dourbes is corrected in real time. Test results have been compared with data from verified neutron monitor stations with a similar geomagnetic cutoff rigidity.     

The role of drift and convection–diffusion mechanisms in small energy global cosmic ray modulation: Application to the hysteresis phenomenon of proton and alpha particle satellite data

The hysteresis effect for small energies of galactic cosmic rays is due to two effects. The first is the same as for neutron monitor energies – the delay of the interplanetary processes responsible for cosmic ray modulation with respect to the initiating solar processes, according to the effective velocity of solar wind and shock waves propagation. Then, the observed cosmic ray intensity is connected to the solar activity variations during many months before the time of cosmic ray measurement. The second is caused by the time delay of small energy cosmic ray diffusion from the boundary of modulation region to the Earth’s orbit. The model describing the connection between solar activity variation and cosmic ray convection–diffusion global modulation for neutron monitor energies is here developed by taking into account also the time-lag of the small energy particle diffusion in the Heliosphere. We use theoretical results on drifts and analytically approximate the dependences of drifts from tilt angle, and take into account the dependence from the sign of primary particles, and from the sign of polar magnetic field (A > 0 or A < 0). The obtained results are applied on proton and alpha-particle satellite data. We analyze satellite 5-min data of proton fluxes with energies >1 MeV, >2 MeV, >5 MeV, >10 MeV, >30 MeV, >50 MeV, >60 MeV, >100 MeV, and in intervals 10–30 MeV, 30–60 MeV, and 60–100 MeV during January 1986–December 1999. We exclude periods with great cosmic ray increases caused by particle acceleration in solar flare events. Then, we determine monthly averaged fluxes, as well as 5-month and 11-month smoothed data. We analyze also satellite 5-min data on alpha-particle fluxes in the energy intervals 60-160 MeV, 160–260 MeV and 330–500 MeV during January 1986–May 2000. We correct observation data for drifts and then compare with what is expected according to the convection–diffusion mechanism. We assume different dimensions of the modulation region (by the time propagation X₀ of solar wind from the Sun to the boundary of modulation region), for X₀ values from 1 to 60 average months, by one-month steps. For each value of X₀ we determine the correlation coefficient between variations of expected and observed cosmic ray intensities (the estimation of cosmic ray intensities values is given in Section 3 by Eq. (9), and the determination of correlation and regression coefficients in Section 3 by Eq. (8)). The dimension of modulation region is determined by the value of X_0 max, for which the correlation coefficient reaches the maximum value. Then the effective radial diffusion coefficient and residual modulation in small energy region can be estimated.     

History of cosmic ray research in Finland

The history of cosmic ray research in Finland can be traced back to the end of 1950s, when first ground-based cosmic ray measurements started in Turku. The first cosmic ray station was founded in Oulu in 1964 performing measurements of cosmic rays by a muon telescope, which was later complemented by a neutron monitor. Since the 1990s, several research centers and universities, such as The Finnish Meteorological Institute, Helsinki University of Technology, University of Oulu, University of Turku and University of Helsinki have been involved in space science projects, such as SOHO, AMS, Cluster, Cassini, BepiColombo, etc. At the same time, ground-based cosmic ray measurements have reached a new level, including a fully automatic on-line database in Oulu and a new muon measuring underground site in Pyhäsalmi. Research groups in Helsinki, Oulu and Turku have also extensive experience in theoretical investigations of different aspects of cosmic ray physics. Cosmic ray research has a 50-year long history in Finland, covering a wide range from basic long-running ground-based observations to high-technology space-borne instrumentation and sophisticated theoretical studies. Several generations of researchers have been involved in the study ensuring transfer of experience and building the recognized Finnish research school of cosmic ray studies.     

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.     

The Tien-Shan mountain cosmic ray station of the Ionosphere Institute of Kazakhstan Republic

The mountain cosmic ray (CR) station of the Ionosphere Institute of Kazakhstan Republic (CR station Alma-Ata B, 43.1N latitude, 76.6E longitude, geomagnetic rigidity cutoff 6.69 GV) is a center for an experimental study of the non-stationary processes caused by cosmic rays of different origin in the interplanetary and near-Earth space, so as their influence on the state of the Earth’s magnetosphere and the upper atmosphere layers. This paper summarizes efforts performed over the years by the Almaty CR group till its present status.