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

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.     

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.     

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 ray modulation with a Fisk-type heliospheric magnetic field and a latitude-dependent solar wind speed

The effect of a latitude-dependent solar wind speed on a Fisk heliospheric magnetic field [Fisk, L. A. Motion of the footpoints of heliospheric magnetic field lines at the Sun: implications for recurrent energetic particle events at high heliographic latitudes. J. Geophys. Res. 101, 15547–15553, 1996] was first discussed by Schwadron and Schwadron and McComas [Schwadron, N.A. An explanation for strongly underwound magnetic field in co-rotating rarefaction regions and its relationship to footpoint motion on the the sun. Geophys. Res. Lett. 29, 1–8, 2002. and Schwadron, N.A., McComas, D.J. Heliospheric “FALTS”: favored acceleration locations at the termination shock. Geophys. Res. Lett. 30, 41–1, 2003]. Burger and Sello [Burger, R.A., Sello, P.C. The effect on cosmic ray modulation of a Parker field modified by a latitudinal-dependent solar wind speed. Adv. Space Res. 35, 643–646, 2005] found a significant effect for a simplified 2D version of a latitude-dependent Fisk-type field while Miyake and Yanagita [Miyake, S., Yanagita, S. The effect of a modified Parker field on the modulation of the galactic cosmic rays. In: Proceedings of 30th International Cosmic Ray Conference. Merida, Mexico, vol. 1, 445–448, 2007] found a smaller effect. The current report improves on a previous attempt Hitge and Burger [Hitge, M., Burger, R.A. The effect of a latitude-dependent solar wind speed on cosmic-ray modulation in a Fisk-type heliospheric magnetic field. In: Proceedings of 30th International Cosmic Ray Conference. Merida, Mexico, vol. 1, pp. 449–450, 2007] where the global change in the solar wind speed and not the local speed gradient was emphasized. The sheared Fisk field of Schwadron and McComas [Schwadron, N.A., McComas, D.J. Heliospheric “FALTS”: Favored acceleration locations at the termination shock. Geophys. Res. Lett. 30, 41–1, 2003.) is similar to the current Schwadron–Parker hybrid field. Little difference is found between the effects of a Parker field and a Schwadron–Parker hybrid field on cosmic-ray modulation, in contrast to the results of Burger and Sello and Miyake and Yanagita [Burger, R.A., Sello, P.C. The effect on cosmic ray modulation of a Parker field modified by a latitudinal-dependent solar wind speed. Adv. Space Res. 35, 643–646, 2005 and Miyake, S., Yanagita, S. The effect of a modified Parker field on the modulation of the galactic cosmic rays. In: Proceedings of 30th International Cosmic Ray Conference. Merida, Mexico, vol. 1, pp. 445–448, 2007]. The two-dimensional approximation used by these authors is therefore inadequate to model the complexities of the actual three-dimensional field. We also show that a Fisk-type field with a latitude-dependent solar wind speed (Schwadron–Parker hybrid field) decreases both the relative amplitude of recurrent cosmic ray intensity variations and latitude gradients and yields similar constants of proportionality for these quantities as for the constant solar wind speed case.     

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

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

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.     

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.     

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.     

Midday bottomside electron density profiles during moderate solar activity and comparison with IRI-2001

Bottomside electron density (Ne-h) profiles during midday (10–14 h) are analyzed using modern digital ionosonde observations at a low-middle latitude station, New Delhi (28.6N, 77.2E, dip 42.4N), for the period from January 2003 to December 2003, pertaining to moderate solar activity (MSA). Each individual profile is normalized with respect to the peak height and density (hmF2, NmF2) of the F2-region. These profiles are compared with those obtained from the International Reference Ionosphere (IRI-2001) model. Bilitza [Bilitza, D. International Reference Ionosphere 2000. Radio Sci. 36 (2), 261–275, 2001] using both the options namely: Gulyaeva’s model [Gulyaeva, T.L. Progress in ionospheric informatics based on electron density profile analysis of ionograms. Adv. Space Res. 7 (6) 39–48, 1987] and B0 Tab. option [Bilitza, D., Radicella, S.M., Reinisch, B.W., Adeniyi, J.O., Mosert Gonzalez, M.E., Zhang, S.R., Obrou, O. New B0 and B1 models for IRI. Adv. Space Res. 25 (1), 89–95, 2000]. The study reveals that during summer and equinox, the IRI model with B0 Tab. option in general, produces better agreement with the observed median profiles, while the IRI predictions using Gulyaeva’s option, overestimate the electron density distribution at all the heights below the F2-peak. However, during winter, in general, the IRI model, using both the options, reveals shows fairly good agreement with the observations.     

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.     

Reminiscences of cosmic ray research in Mexico

Cosmic ray research in Mexico dates from the early 1930s with the work of the pioneering physicist, Manuel Sandoval Vallarta and his students from Mexico. Several experiments of international significance were carried out during that period in Mexico: they dealt with the geomagnetic latitude effect, the north–south and west–east asymmetry of cosmic ray intensity, and the sign of the charge of cosmic rays. The international cosmic ray community has met twice in Mexico for the International Cosmic Ray Conferences (ICRC): the fourth was held in Guanajuato in 1955, and the 30th took place in Mérida, in 2007. In addition, an international meeting on the Pierre Auger Collaboration was held in Morelia in 1999, and the International Workshop on Observing UHE Cosmic Rays took place in Metepec in 2000. A wide range of research topics has been developed, from low-energy Solar Energetic Particles (SEP) to the UHE. Instrumentation has evolved since the early 1950s, from a Simpson type neutron monitor installed in Mexico City (2300 m asl) to a solar neutron telescope and an EAS Cherenkov array, (within the framework of the Auger International Collaboration), both at present operating on Mt. Sierra La Negra in the state of Puebla (4580 m asl). Research collaboration has been undertaken with many countries; in particular, the long-term collaboration with Russian scientists has been very fruitful.     

Flight safety implications of the extreme solar proton event of 23 February 1956

There is considerable speculation about the effects at aircraft altitudes resulting from extreme solar proton events. The ground level event (GLE) of 23 February 1956 (GLE 5), remains the largest solar proton event of the neutron monitor era in terms of its influence on count rates at monitors near sea level. During this GLE the count rate was increased by as much as 4760% (15-min average) at the Leeds monitor relative to the count rate from galactic cosmic radiation (GCR). Two modern models of the event cumulative solar proton spectrum for this event, a 6-parameter fit in energy and a 4-parameter Band fit in rigidity, are compared with 1-h of GCR at solar minimum. While effective doses calculated with CARI-7A for both models at low geomagnetic cutoff rigidities are indeed high when compared with GCR and can exceed recommended exposure limits, both GLE spectra exhibit a much stronger dependence on cutoff rigidity than GCR, and a larger fraction of the dose from neutrons. At locations with cutoff rigidities above 4.2 and 6.4?GV, respectively, the GLE effective doses are smaller than the GCR hourly dose. At locations with cutoff rigidities above about 4?GV, GCR was the dominant source of exposure in 10?h or less at all altitudes examined. This suggests that if a similar event occurs in the future, low- and mid-latitude flights at modern jet flight altitudes could be well-protected by Earth’s magnetic field.     

Mid 19th century minimum of galactic cosmic ray flux inferred from Ti in Allegan meteorite

Measurements of ⁴⁴Ti activity in meteorites show that the galactic cosmic ray (GCR) intensity has been declining in the interplanetary space during the past three centuries and has a component of cyclic variation, with periodicity of about 87 years [Taricco, C., Bhandari, N., Cane, D., et al. Galactic cosmic ray flux decline and periodicities in the interplanetary space during the last 3 centuries revealed by ⁴⁴Ti in meteorites. J. Geophys. Res. 111, A08102, 2006.]. In order to verify these results, we have measured ⁴⁴Ti activity in Allegan meteorite which fell in 1899 and in some other meteorites with better precision. The measurements confirm low cosmic ray flux and consequently high solar activity near the middle of 19th century.     

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

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

Energy dependence of the rigidity spectrum of Forbush decrease of galactic cosmic ray intensity

We show that rigidity spectrum of Forbush decrease (Fd) of galactic cosmic ray (GCR) intensity in September 9–23, 2005 clearly depends on energy. We calculated rigidity spectrum of the Fd based on the neutron monitors and Nagoya muon telescope channels’ data divided in three groups according to their cut off rigidities. We found that temporal changes of rigidity spectrum exponent γ are approximately similar for all cut off rigidity groups, but γ values are the larger the higher are cut off rigidities. We conclude that rigidity spectrum of Fd is hard for lower energy range and is soft for the higher energy range. We believe that an energy dependence of the power law rigidity spectrum of Fd is observed owing to the preferential convection–diffusion mechanism during Fd in September 9–23, 2005. It is a reflection of an influence of the temporal changes of the structure of the interplanetary magnetic field (IMF) turbulence in different range of frequency f during Fd. Particularly, a decisive role in formation of the character of the rigidity spectrum belongs to the changes of the exponent ν of the power spectral density (PSD) of the IMF turbulence (PSD ∝ f^−ν). The exponent ν is greater for high frequency region of the IMF turbulence (responsible for scattering of low rigidity particles of GCR), than for low frequency region of the IMF turbulence (being responsible for scattering of higher rigidity particles). Also, we challenge to estimate an existence of slab/2D structure of solar wind turbulence during the Fd in September 9–23, 2005 based on the distribution of average turbulence energy among the IMF’s components.     

Mean energy of response to galactic cosmic ray spectrum: IMP 8 and Climax neutron monitor

Solar modulations of galactic cosmic ray (GCR) intensity contain a wealth of information about their transport in the heliosphere. To extract this information from the data one studies the dependence of the observed modulations on the mean energy of response of detectors providing data for the analyses. There is a great deal of confusion about the detector energy response to GCR spectrum in the literature. We present a preliminary report on the computations of the mean energy of response for the Climax neutron monitor (CL/NM) and IMP 8 cosmic ray nuclear composition instrument to GCR protons for 1973–1998, covering the solar cycles 21 and 22. We find that for penetrating proton channel on IMP 8 the mean energy changes by a factor of over two whereas for the neutron monitor the change is only 21%. However, the corresponding change for the computed modulation function is a factor of about 3.5.     

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.