SH 4.6.4

NEUTRON MONITOR SENSITIVITY TO SOLAR MODULATION CHANGES: ALTITUDE VS. CUTOFF RIGIDITY

K. R. Pyle
Laboratory for Astrophysics and Space Research, Enrico Fermi Institute,
University of Chicago, 933 E. 56th St., Chicago, IL 60637-1460, USA

ABSTRACT
This paper examines the effects of altitude and vertical cutoff rigidity upon the sensitivity of neutron monitors to changes in the solar modulation of cosmic rays. The counting rate response of a given monitor is determined by a) the geomagnetic latitude effect, which controls the access of primaries to the top of the atmosphere, and b) the amount of atmosphere above the monitor, which determines the properties of the cascade at the monitor (the specific yield function, first defined by Simpson et al., 1953). The magnitude of the effect of solar modulation on the level of the Galactic flux at Earth above ~1 GV is inversely related to particle rigidity. It is often assumed that monitors with the lowest geomagnetic vertical cutoff rigidity (Vcr) should measure the greatest variation in counting rate. However, for a given Vcr, when the nucleonic yield functions are taken into account, a monitor which is at high altitude exhibits significantly greater sensitivity to changes in the galactic flux level than one near sea level. In this paper we compare the Climax neutron monitor (Vcr = 3.0 GV, Alt. = 3400 M) with Deep River (Vcr = 1.0 GV, Alt. = 145 M), and show that the counting-rate variation at Climax is 16 percent larger than at Deep River (~30% vs. ~26% over an 11-year solar cycle), although Deep River has a much lower Vcr than Climax. Similar differences in response sensitivity are found for shorter time intervals.

INTRODUCTION
The study of the long- and short-term modulation of galactic cosmic rays above ~3 GV has relied for many decades upon the world-wide network of neutron monitors, which measures the nucleonic cascades in the Earth's atmosphere and provides the only long-term measurements of the high-energy cosmic ray flux. These monitors, first developed by John Simpson around 1948-49 (Simpson et al., 1953), consist of an array of BF3 neutron counters surrounded by a structure composed of lead (to generate local neutron showers) and paraffin or polyethylene (to moderate the neutron velocity so that they can be measured efficiently by the counters). Simpson's original monitor design (Simpson, 1958), used throughout the world during and after the 1957-58 International Geophysical Year (IGY), became known as the IGY monitor. This design was later enlarged during the 1960s (supermonitors, which use much larger BF3 tubes than the IGY design). The Climax monitor, which has been in continuous operation since 1950, provides the longest baseline of low-rigidity (>3 GV) galactic cosmic ray data in the world. Our other station, on Mt. Haleakala in Hawaii, when combined with the now-closed Huancayo monitor, provides the longest base-line of high-rigidity (>13 GV) data. The University of Chicago monitors have been described recently in Pyle (1993). Our neutron monitor data can be accessed at the following URL: http://astro.uchicago.edu/ home/web/pyle/neutron.html..

The counting-rate response of a given monitor (IGY or supermonitor) is determined by a) the geomagnetic latitude effect, which controls the access of primaries to the top of the atmosphere, where the cascade begins, b) the design (type and number of tubes), and c) the pressure-altitude of the monitor, which determines the properties of the cascade at the monitorŐs elevation. The development of the nucleonic cascade begins in the upper atmosphere; it grows and is then attenuated with increasing atmospheric depth.

At neutron monitor energies ( 500 MV) the magnitude of the effect of solar modulation on the level of the Galactic flux at Earth is inversely related to particle energy or rigidity. For instance, a monitor which has a geomagnetic vertical cutoff rigidity (VCR) of 2 GV would be expected to exhibit a larger percentage variation for a given 'modulation eventŐ than a monitor with a VCR of 3 GV all other factors being equal. It is therefore often assumed that monitors with the lowest vertical cutoff rigidity (Vcr) (i.e. those at high geomagnetic latitude) should measure the greatest variation in counting rate. However, when the simultaneous development and attenuation of the cascade (the Specific Yield Function, Simpson et al. 1953) is taken into account, the altitude of the monitor can be seen to play a very important role, especially at low to medium Vcrs (i.e. <5-6 GV). In fact, as will be demonstrated by several examples below, a mid-latitude high-altitude monitor such as ours at Climax, Colorado (Elevation 3400 M, Pcv ~3 GV) is one of the most sensitive in the world to changes in modulation.

COSMIC RAY YIELD FUNCTIONS
Historically, primary proton yield functions have been calculated using either latitude surveys (transporting a monitor through a range of magnetic latitudes and altitudes - e.g. Simpson, 1948; Simpson et al., 1953; Rose et al., 1956; Lockwood and Webber, 1967; Raubenheimer and Stoker, 1974; Bieber and Evenson, 1997) or Monte-Carlo simulations of the cascade produced by the interaction of a primary proton in the upper atmosphere (e.g. Webber and Quenby, 1958; Debrunner and Flückiger, 1971a, 1971b).

It is not the purpose of this paper to examine in detail the various measurements and calculations of the primary proton yield function, but rather to call attention to the fact that the shape of the yield function varies appreciably with altitude, and that this is often the determining factor for sensitivity. All studies of the yield function exhibit this shape change with altitude, and we have chosen just one typical example for illustration in Figure 1, adapted from DeBrunner and Flückiger (1971b). In it we show the yield functions for an IGY monitor at two altitudes: sea level and 650 g-cm-2. The latter pressure-altitude is close to that for our Climax neutron monitor (685 g-cm-2). The major point to be made here is that the higher pressure-altitude curve has a maximum at a significantly lower rigidity (~4.5 GV) than the sea level curve (~7 GV). The annotations that we have added to Figure 1 show the vertical cutoff rigidities for 4 selected neutron monitors (Climax, Haleakala, Deep River, McMurdo, and South Pole). Simply put, the counting rates registered by these monitors are the integrals to the right of the vertical lines; it is apparent from these curves that a high altitude mid-latitude monitor such as Climax responds to a significantly lower mean primary rigidity than any sea-level monitor, including polar monitors. However, equatorial high-altitude monitors (such as our IGY and supermonitor on Mt. Haleakala in Hawaii, would not exhibit significantly larger sensitivity than a sea-level equatorial monitor. In the following section, we will provide examples of the "sensitivity" of the Climax monitor to changes in the modulation level by comparisons with low-altitude monitors over several time-scales, from several days (Forbush decreases) to variations over several 11-year solar cycles.

DATA COMPARISONS
Data used in this paper are taken from the following neutron monitors: Climax, Colorado, USA, operated by the Univ. of Chicago; Deep River, Canada, operated by the National Research Council, Canada; McMurdo and South Pole, Antarctica, operated by the Univ. of Delaware. Non-Chicago data were obtained from World Data Center A, Boulder, Colorado, USA. Due to lack of space in this paper, other high-altitude, sensitive monitors can not be included, notably Mt. Washington, Jungfraujoch, and others.

In Figure 2 we have plotted a typical Forbush Decrease (2/26/92) as seen at Climax and Deep River, along with their ratio. For this event, the Climax "sensitivity" is ~16% greater than Deep River; the ratio of the two counting rates drops.

Figure 3 is a comparison of these same two stations over the nearly 4 decades of simultaneous operation (unfortunately, the Deep River monitor was shut down in 1995). All data has been normalized to the 1965 level, defined as 100%.

In this Figure we have used Climax data which has been adjusted after 1980/81 to correct normalization problems which occurred at that time. (Our thanks to H. Ahluwahlia for calling our attention to this problem). Relative to the previously published counting rates for Climax, this renormalization results in an adopted multiplicative factor of 0.9891 from 16 Sept. 1980 through 31 Aug. 1981, and 1.0121 thereafter. These corrections will be published in the near future.

The top panel displays monthly averages of the counting rates, while the bottom panel displays the ratio between the two monitors. Over the period 1958-1995 the ratio varies only ~4%, while the individual counting rates vary by more than 30%. In Figure 4 we plot the two sets of monthly averages against each other and show the best-fit straight line. Based on this, the Climax monitor is again 16% more "sensitive" to changes in solar modulation. This sensitivity appears to be independent of the phase of the solar cycle, in that the fit appears to equally good at both solar minimum and maximum.

As a final illustration of the importance of atmospheric mass to the sensitivity to modulation of neutron monitors, we compare two monitors with essentially the same Pcv but very different elevations (McMurdo and South Pole), which demonstrates the altitude effect very clearly. In Figure 5 we show counting rates from these two monitors during the same Forbush decrease used in Figure 2. Again, the increased sensitivity with elevation is very apparent.

SUMMARY
We have demonstrated the importance of elevation or pressure-altitude in determining the sensitivity of a neutron monitor to changes in modulation, especially for monitors with geomagnetic cutoffs below the peak of the yield curve, about ~5-7 GV. For higher values of the cutoff rigidity, the sensitivity of the monitor is much less dependent on altitude.

ACKNOWLEDGMENTS
To J. Simpson, for many helpful comments. Research supported by NSF Grant ATM-96-13963.

REFERENCES
Bieber, J. W., and Evenson, P., Proc. XXV ICRC (Durban), Paper SH 4.2.3 (1997).
Debrunner, H. and Flückiger, E., Helv. Phys. Acta 44, 241 (1971).
Debrunner, H. and Flückiger, E., Proc. 12th ICRC (Hobart), 3, 911 (1971).
Lockwood, J. A. and Webber, W. R., J. Geophys. Res. 72, 3395 (1967).
Pyle, K. R., Proc. 23rd ICRC (Calgary), 3, 609 (1993).
Raubenheimer, B. C. and Stoker, P. H., J. Geophys. Res. 79, 5069 (1974).
Rose, D. C., Fenton, K. B., Katzman, J. and Simpson, J. A., Can. J. Phys., 34, 968 (1956).
Simpson, J. A., Phys. Rev. 73, 1389 (1948).
Simpson, J. A., Fonger, W. and Treiman, S. B., Phys. Rev. 90, 934 (1953).
Simpson, J. A., Annals of the Intl. Geophy. Yr, IV, Parts I-III, 351, Pergamon Press, London, 1958.
Webber, W. R. and Quenby, J. J., Phil. Mag. 4, 654 (1959).