Modeling the Relationship Between Neutron Counting Rates and Sunspot Numbers Using the Hysteresis Effect

Abstract

Several studies show that temporal variations in the Galactic cosmic ray (GCR) intensity display a distinct 11-year periodicity due to solar modulation of the galactic cosmic rays in the heliosphere. The 11-year periodicity of GCRs is inversely proportional to, but out of phase with, the 11-year solar cycle, implying that there is a time lag between actual solar cycle and the GCR intensity, which is known as the hysteresis effect. In this study, we use the hysteresis effect to model the relationship between neutron counting rates (NCRs), an indicator of the GCR intensity, and sunspot numbers (SSNs) over the period that covers the last four solar cycles (20, 21, 22, and 23). Both linear and ellipse models were applied to SSNs during odd and even cycles in order to calculate temporal variations of NCRs. We find that ellipse modeling provides higher correlation coefficients for odd cycles compared to linear models, e.g. 0.97, 0.97, 0.92, and 0.97 compared to 0.69, 0.72, 0.53, and 0.68 for data from McMurdo, Swarthmore, South Pole, and Thule neutron monitors, respectively, during solar cycle 21 with overall improvement of 31 % for odd cycles. When combined to a continuous model, the better correlation observed for the odd cycles increases the overall correlation between observed and modeled NCRs. The new empirical model therefore provides a better representation of the relationship between NCRs and SSNs. A major goal of the ongoing research is to use the new non-linear empirical model to reconstruct SSNs on annual time scales prior to 1610, where we do not have observational records of SSNs, based on changes in NCRs reconstructed from 10Be in ice cores.