Issue 
J. Space Weather Space Clim.
Volume 6, 2016



Article Number  A15  
Number of page(s)  7  
DOI  https://doi.org/10.1051/swsc/2016011  
Published online  04 March 2016 
Research Article
An empirical model of heliospheric cosmic ray modulation on longterm time scale
^{1}
ReSoLVE Centre of Excellence, University of Oulu, 90014
Oulu, Finland
^{2}
Sodankylä Geophysical Observatory, University of Oulu, Oulu, FIN90014, Finland
^{*} Corresponding author: Eleanna.Asvestari@oulu.fi
Received:
5
November
2015
Accepted:
5
February
2016
Galactic Cosmic Rays (GCRs) entering the heliosphere are subject to modulation processes due to variable solar magnetic activity. Finding a relationship between cosmicray variations and the heliospheric parameters is important for reconstruction of solar activity in the past. Here, we develop a semiempirical model describing the heliospheric modulation of GCRs in terms of heliospheric parameters such as the open solar magnetic flux, the tilt angle of the heliospheric current sheet and the polarity of the large scale solar magnetic field. Our model is fitted using annual data obtained for the period 1976–2013, which includes the very weak solar minimum during 2008–2010. The model shows a good agreement with the data, and therefore, can be used for reconstructions of the modulation potential at different levels of solar activity. The model’s validity is also tested using the cosmogenic radionuclides ^{14}C and ^{10}Be stored in terrestrial archives. The tilt angle used to fit the parameters in our semiempirical modulation model is reconstructed by a mathematical model described here.
Key words: Heliosphere / Cosmic ray / Current sheet / Solar activity / Cosmogenic isotopes
© E. Asvestari and I.G. Usoskin, Published by EDP Sciences 2016
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
A theory describing Galactic Cosmic Rays (GCRs) transport in the heliosphere (Parker 1965; Krymskij 1969) includes four main processes: diffusion of particles along and perpendicular to the heliospheric magnetic field (HMF), outward convection and adiabatic deceleration in the expanding solar wind (SW), gradient and curvature drifts in the large scale HMF and drift due to the tilted heliospheric current sheet (HCS). This interplay results in spatial and temporal variations of GCRs intensity and energy.
Despite essential progress in theoretical development of the heliosperic transport models (see a review by Potgieter 2013 and references therein), a semiempirical approach (e.g. Cliver 1993; Belov 2000; AlankoHuotari et al. 2006) is still useful in many practical studies. Fully developed numerical 3D timedependent models of cosmic rays heliospheric transport do exist (e.g. Jokipii & Thomas 1981; Kota & Jokipii 1983; Hattingh & Burger 1995; Potgieter et al. 2001; Ferreira & Potgieter 2004; Potgieter 2013). However, such models include parameters that cannot be directly measured, such as, e.g., the diffusion tensor or the spectrum of turbulence of the solar wind magnetic inhomogeneities. Therefore, it is difficult to apply realistic computations using the sophisticated 3D models which are often based on ad hoc parameterisations (e.g. Ferreira & Potgieter 2004). As an alternative way, a number of (semi)empirical models are based on observations, aiming to associate cosmicray variations with various solar and heliospheric parameters. Some of these empirical models are limited either by using cosmic rays intensity at a fixed energy (Stozhkov et al. 2004), or time series obtained by a single neutron monitor (Sabbah & Rybansk 2006), or neglecting possible longterm solar variability (Belov et al. 2006).
In this work, we use an approach similar to that by AlankoHuotari et al. (2006), improved by using revised relations and the extended database of the heliospheric parameters for the period 1951–2013. Earlier models were developed based on observations corresponding to the period of high levels of solar activity, the socalled Modern Grand maximum of solar activity (Solanki et al. 2004). Accordingly, they may not well represent a possible centennial trend and periods with low solar activity. We fit our model for the period including the latest very weak solar cycle, which is similar to low activity levels in the past, such as the Dalton Minimum (DM) ca. 1800 or the Modern minimum ca. 1900 (McCracken & Beer 2014). Thus, our model aims at reconstructing the GCRs behaviour during periods with different solar activity levels.
We reconstruct, by applying our model, the modulation potential since 1616, a parameter that describes the mean energy loss of GCR particles within the heliosphere due to solar modulation, and using that, we model the global production of radionuclides ^{14}C and ^{10}Be, adopting the method described in Kovaltsov et al. (2012) and Kovaltsov & Usoskin (2010), respectively. The modelled series are next compared with the records from terrestrial archives such as tree rings (Roth & Joos 2013) and ice cores (Berggren et al. 2009).
2. Reconstruction of the HCS tilt angle
The tilt of the HCS is an important factor affecting modulation of galactic cosmic rays in the heliosphere (e.g., Jokipii & Thomas 1981; Kota & Jokipii 1983). Therefore, it is important to know its variability in order to study large scale HMF dynamics and the heliospheric modulation of cosmic rays on long time scale including the centennial trends. Although continuous observations of the tilt angle are made by the Wilcox Solar Observatory (WSO) since 1976, for the purpose of centennial reconstructions of the CR modulation it is important to know the tilt angle variability on longer time scale. Here we develop a model to reconstruct the HCS tilt angle using the phase of the solar cycle. This model describes the cyclic behaviour of the tilt angle (Hoeksema 1991; Suess et al. 1993) and its dependence solely on the phase of the solar cycle.
Variability of the HCS tilt angle follows the 11year solar cycle along with the sunspot number (Fig. 1). However, although the latter varies significantly from cycle to cycle (as reflected in different cycle magnitudes varying by a factor of 2), the tilt angle exhibits roughly the same variations over all cycles, in accord with the idea of a regular cyclic behaviour (Suess et al. 1993; Cliver & Ling 2001; AlankoHuotari et al. 2007), depending only on the solar cycle phase and not on its strength. This makes it possible to reconstruct the tilt angle variability using only the solar cycle phase. Moreover, strong variations in sunspot number amplitude during the solar maxima indicate that the sunspot number cannot be a good direct proxy for the tilt angle.
Fig. 1. Annual variations of the yearly averaged sunspot numbers (Solar Influences Data Analysis Center (SIDC) at the Royal Observatory of Belgium) (version 1.0) (blue solid line) and the HCS tilt angle, α, (green dashed line) for the period 1976–2014. The correlation coefficient between the two is . 
Figure 2 shows a superposed plot of the tilt angle as a function of the cycle phase for solar cycles 21 through 24. By comparing the tilt angle values from cycle to cycle, one can notice that the tilt angle has a similar behaviour during the ascending phase with a fast and smooth increase, for all the cycles. During the descending phase, the tilt angle fluctuates essentially, and these variations differ from cycle to cycle. However, the decrease is gradual in all cycles supporting the idea of a cyclic behaviour. Cliver & Ling (2001) concluded, from a similar analysis, that the descending phase of the tilt angle cycle is more gradual for odd cycles, but that is not observed in Figure 2, and moreover the statistic is too low (only two odd cycles) to substantiate such a conclusion.
Fig. 2. Yearly averaged HCS tilt angle variations based on the model (Eq. (1)) and the WSO observations of the photospheric magnetic field in relation to the cycle phase, for solar cycles 21–24 as indicated in the legend. 
AlankoHuotari et al. (2007) developed an empirical relation describing the cyclic behaviour of the tilt angle for the period 1976–2005. With the WSO tilt angle calculations extending until present, we revisited the model and fitted new parameters. Here we used, as the HCS tilt angle, the maximum latitudinal extent of the HCS, as provided by the WSO (http://wso.stanford.edu/Tilts.html). WSO offers two different approaches to estimate these values, using their photospheric magnetic field observations: One is the socalled “classic” model, and the other is the newer radial model. The “classic” model applies lineofsight boundary conditions to the photosphere and therefore, it requires polar field corrections. Such corrections are not necessary for the radial model, since it applies radial conditions. Also, the radial model computes results with a higher source surface radius, and accordingly, the current sheet appears flatter and has a lower maximum extent. Here we use the values defined by the radial model, similar to AlankoHuotari et al. (2007), in view of its more robust definition (Wang & Sheeley 1992). Averages of the values produced by WSO over a calendar year are derived and considered as annual values.
The bestfit empirical model, describing the annual value of the HCS tilt angle based on the progress of the solar cycle, takes the form:(1)where X_{i} = i/N, i = [1, N] is the year in the cycle and N is the cycle length in years. The value X_{i} represents the phase of a cycle. We note that this relation is applied only to annual values. The length of the solar cycle 24 was defined, assuming the minimum of the cycle 24 in 2021 (Uzal et al. 2012). For this model, the ascending and maximum phases are considered to be shorter than the descending phase, reflecting the asymmetry of the tilt angle cyclic shape (Hathaway 2010). Although observational limitations restrict the maximum angle to be 70° (Suess et al. 1993), this does not affect the use of this model’s reconstructions for the study of GCRs heliospheric modulation, since convectivediffusive processes are the main contributors to cosmic rays modulation for tilt angle values greater than 50° (Lopate & Simpson 1991; Potgieter & Le Roux 1992; Cliver et al. 2013). The cyclic behaviour of the model is consistent with the cyclic shape followed by the WSO estimated values (Fig. 2).
In order to test how well the tilt angle values, derived by this empirical relation, agree with the ones obtained by WSO, the root mean square (RMS) error between the two sets was estimated. The model was applied for the years 1976–2013, and the corresponding R_{rms} is 6.6°. We note that this is a significant improvement with respect to the earlier model (AlankoHuotari et al. 2007) which yields R_{rms} = 11.2° for annual values. The fact that the model with more data gives a better fit to the data confirms its convergence and validity. Figure 3 shows the annual variations of the (radial) tilt angle provided by WSO and the one calculated here (Eq. (1)). The two quantities are well correlated () with the best correlation found during the ascending phase (Fig. 4). Some indirect sparse tilt angle estimates, based on image analysis of solar eclipses from 1870 through 2002, are done by Pishkalo (2006). We also show them in the same figure for comparison. They all appear to be in good agreement (see Fig. 3), confirming the robustness of the model.
Fig. 3. Temporal variations of the HCS tilt angle based on the radial model by WSO (blue line) and reconstructed here using the empirical model (Eq. (1)) (magenta steplike line). Black stars represent reconstructed HCS tilt angles by an analysis of total solar eclipse images (courtesy of M.I. Pishkalo). 
Fig. 4. Scatter plot of the modelled here versus the WSO provided values of the annually averaged HCS tilt angle for the period 1976–2014 (as shown in Fig. 3). The blue solid line marks the diagonal. 
3. Reconstruction of the modulation potential
Observable heliospheric parameters involved in the modulation processes described in Parker’s theory are the solar wind velocity and density, the HCS tilt angle, the polarity, the strength and the level of turbulence of the HMF. These parameters vary with solar activity and the solar cycle phase.
In our model, we consider the heliospheric parameters that satisfy the following two main criteria (AlankoHuotari et al. 2006). Firstly, they should describe the global heliosphere and secondly, they should be recorded or calculated continuously for a long period. Parameters that satisfy these criteria are the open solar magnetic flux (F), the tilt angle of the HCS (α) and the HMF polarity (p). The solar wind velocity is not considered within these parameters since its variations at 1 AU, in the ecliptic plane, do not represent the global behaviour of the solar wind, especially at high heliolatitudes. The open solar magnetic flux is used in the model as a global parameter of the HMF, instead of the HMF strength, B, measured in the ecliptic plane, since the latter is local and poorly represents the global HMF distribution. The time variations of these parameters are shown in Figure 5.
Fig. 5. Annual values of the heliospheric parameters. From top to bottom are the modulation potential ϕ (MV), open solar magnetic flux F (10^{15} Wb), calculated tilt angle α (deg) and polarity p for the period 1951–2013. 
Parker’s transport equation can be reduced, under some simplifying yet realistic assumptions, to the socalled force field approximation (Gleeson & Axford 1968; CaballeroLopez & Moraal 2004), where the modulation is described by a single parameter, the modulation potential, ϕ, which parametrises the shape of the GCR energy spectrum (see formalism in Usoskin et al. 2005) and is expected to be inversely proportional to the diffusion coefficient, κ, of the heliospheric transport of GCRs, to some power n.
Taking into consideration that the diffusion coefficient, κ, is supposed to be roughly inversely related to the strength of the interplanetary magnetic field (Wibberenz et al. 2001b), which in turn is directly related, but not necessarily linearly proportional, to the OSF, one can assume that the modulation potential is roughly proportional to the OSF, ϕ ∝ F ^{ n } (Usoskin et al. 2002). In some studies n is taken to be equal to unity (e.g. Kta & Jokipii 2001), resulting in a linear relationship between the modulation potential and OSF. On the other hand, some studies suggest that the relation is not exactly linear or quasilinear (Mursula et al. 2003; AlankoHuotari et al. 2006). AlankoHuotari et al. (2006) studied the correlation between ϕ and F over periods of flat or highly tilted HCS and showed that the dependence of ϕ on F varies with the tilt angle. This is supported by other studies according to which n can differ from unity (e.g. CaballeroLopez & Moraal 2004), having a dependence on the solar cycle phase (Wibberenz et al. 2001b; Ferreira & Potgieter 2004) and vary with the HCS tilt angle (Wibberenz et al. 2001a; Ferreira & Potgieter 2004).
The role of the HMF polarity in the GCR modulation within the heliosphere varies with the phase of solar cycle and subsequently with the flatness/waviness of the HCS. During periods of minimum solar activity, when the HCS is flat, particles drift along the equatorial plane and away from the Sun for positive polarity periods, while the opposite occurs for negative polarity periods. For the negative polarity periods the modulation becomes stronger with increasing tilt angle (AlankoHuotari et al. 2007; Dorman 2006 and references therein). During positive polarity periods an increasing tilt angle leads to declining modulation. However, the effects of the HCS on the modulation potential are not significant when the solar cycle reaches its maximum.
Accordingly, we use a model relating the heliospheric parameters to the modulation potential in the following form:(2)where the midterm defines the dependence of ϕ on the OSF, that is slightly modulated by the HCS tilt angle, and the last term parametrises the drift effect of the HCS. In this model ϕ_{0}, α_{0}, n and β are free parameters.
In order to find the bestfit parameters that minimise the mean square differences between the observed and the modelled ϕ, we performed a nonlinear model fitting using the annual modulation potential values, calculated from ground based cosmic ray observations (Usoskin et al. 2005, 2011), the OSF reconstructed by Lockwood et al. (2013a, 2013b, 2014a, 2014b) and the tilt angle as described in Section 2. The polarity used for the fitting is defined as p = 1 for positive polarity periods, p = −1 for negative polarity periods and p = 0 for the years when polarity reversal occurs. Based on the Wilcox Solar Observatory (WSO) polar field observations since 1976, the polarity reversal years are taken to be 1980, 1990, 2000, 2013. For the years prior to 1976, when no direct data are available, we assume that the polarity reverses at the year of sunspot maximum. The fitting was done for the period 1951–2013, and the best fitted parameters are ϕ_{0} = 1473.9 MV, α_{0} = 150°, n = 1.03 and β = 0.095.
The annual variation of the computed modulation potential is shown in Figure 6 (magenta dashed curve). In the same figure the annual reconstructed values of the modulation potential by ground based cosmic ray data are plotted (blue curve). The correlation coefficient between the two is R = 0.88 ± 0.03 (Fig. 7). The modelled curve appears to follow well the one based on neutron monitor observations. However, a discrepancy appears during the maximum of the solar cycle 22. This may be related to the high plasma flow pressure during years 1991 and 1992. Though solar wind streams are an important factor of cosmic rays modulation, they could not be evaluated in the past and therefore it is difficult to include them in the model, which is intended for longterm studies.
Fig. 6. Temporal variations of the modulation potential reconstructed by ground based observations of cosmic rays (blue curve) and the modelled (magenta dashed curve) over the period 1951–2013. 
Fig. 7. Scatter plot of the data shown in Figure 6. The solid line represents the diagonal. The correlation coefficient, R, between the two parameters is approximately 0.88 ± 0.03. 
The discrepancy during the maximum of solar cycle 24 can be explained in terms of the polarity dependence of the modulation potential in Eq. (2), and the polar field reversal around that period. According to WSO, the polarity reversed in the end of 2013. Therefore, in our model p was set to zero on that year and negative prior. Both from a theoretical and observational perspective, for large tilt angles, i.e. around solar maxima, GCR particles of all energies experience larger modulation during negative polarity epochs comparing to positive polarity epochs. And as we expect, our model yields high modulation, for p = −1. However, based on neutron monitor measurements the modulation potential was lower. This can be explained considering the unusual polarity reversal during this maximum. More precisely, the maximum of solar cycle 24 was characterised by a rather slow and intermittent reversal of the polar field, with northsouth asymmetry, meaning that the Northern polar field reversed on November 2012 and the Southern on March 2014, with both polar fields being simultaneously positive for more than a year (Sun et al. 2015). Therefore, although the polarity is officially reversed from negative to positive at the end of 2013, it was positive for a rather long time since 2012, leading to the lower level of modulation recorded by neutron monitors.
Here, we explain the discrepancy over the last solar maximum considering the parameters involved in our model, and more specifically the polarity and the model’s response to such a change in the input data. Considering that during periods of high solar activity, other timedependent transient phenomena, such as CMEs and propagating interaction regions (Burlaga et al. 1993; Wang et al. 2006), have an essential contribution to GCRs modulation, it is worth mentioning that they might be a source of differences between the modelled and observed series in Figure 6. More sophisticated models are based on these observable parameters (e.g. Burlaga et al. 1993) and describe well the modulation during solar maxima over the past few solar cycles. Although we acknowledge their importance, such phenomena were not recorded in the past, over time scales for which we want to reconstruct the modulation potential and thus they are not considered by our model.
4. Centennial reconstructions and radionuclide production
Cosmic rays particles interact in the atmosphere with the nuclei of atmospheric gases, leading to the production of cosmogenic radionuclides, such as ^{14}C and ^{10}Be. Solar variability and the geomagnetic field affect the flux of cosmic rays entering the atmosphere and subsequently, the radionuclide production rate (Beer et al. 2012). After their production, the radionuclides, following different distribution and deposition processes, are eventually stored in terrestrial archives such as tree rings or ice cores. The radionuclide signal reflects mainly production (Berggren et al. 2009), making them a useful tool of solar variability reconstructions in the past (Usoskin 2013). However, contribution from local climate to ^{10}Be deposition may be essential at interannualdecadal time scales (Pedro et al. 2006, 2011; Usoskin et al. 2009).
We reconstructed the modulation potential since 1616, using the semiempirical modulation model (Eq. (2)), where for the OSF we used the one published by Lockwood et al. (2014a) and Lockwood & Owens (2014), the polarity p was used as defined in Section 3 and the tilt angle α used, is a series reconstructed using Eq. (1). To define the length of the solar cycles since 1616, and subsequently the cycle phase input in Eq. (1), we considered the sunspot number maxima and minima as listed in SIDC, version 1.0. The obtained modulation potential series was applied first to the global radiocarbon production rate model (Kovaltsov et al. 2012), and then to the ^{10}Be flux model (Kovaltsov & Usoskin 2010). In Figure 8, we compare the calculated global radiocarbon production rate with the reconstruction by Roth & Joos (2013) based on tree ring records. The two curves agree well prior to 1900. Due to fossil fuel burning (the socalled Suess effect), since the late 19th century, and atomic bomb testing after 1950, the ^{14}C records cannot be used for the modern epoch.
Fig. 8. Annual variations of radiocarbon production reconstructed by Roth & Joos (2013) (magenta curve, with 95% of confidence interval shade) and modelled here (blue curve). 
We also compared (Fig. 9) the calculated ^{10}Be flux, with the one reconstructed using the ice core records from Greenland NGRIP (Berggren et al. 2009). Our computations show similar longterm trends to the reconstruction by Berggren et al. (2009). It is worthy to note that the two curves show an excellent agreement with each other over the 20th century. Prior to this period, the two curves appear not to show similar solar cycle variations, which is likely due to regional climatic effects (Usoskin et al. 2009) leading to strong fluctuations of the ^{10}Be deposition (Pedro et al. 2006, 2011).
Fig. 9. Annual variations of radionuclide ^{10}Be flux by Berggren et al. (2009) (magenta curve) and modelled here (blue curve). Although the two curves show a very good agreement after 1900, they vary from each other significantly before that period. This is likely due to atmospheric distribution and deposition processes that affect the stored amount of the radionuclide in ice. 
5. Conclusions
In this paper we introduce an empirical model to reconstruct the HCS tilt angle annual variations and a semiempirical model for centenial reconstructions of the modulation potential. The HCS tilt angle model depends only on the solar cycle phase and describes the cyclic behaviour of the HCS tilt angle, reflecting its asymmetric shape, with short length and fast ascending phase in contrast to the gradual descending phase. The reconstructed series agree well with the observations by WSO as well as with sparse estimates from image analysis of solar eclipses (Pishkalo 2006). Due to observational constraints there is an upper limit in the tilt angle maximum at 70°, but this has no effect on the application of the reconstructed tilt angle in studying the heliospheric modulation of GCRs.
The semiempirical modulation model relates the GCR’s modulation potential with those heliosperic parameters measured at 1 AU and describes the global heliospheric solar activity. These are the open solar magnetic flux, the HCS tilt angle and the solar magnetic field polarity. The model fitting is done for the period 1951–2013 and gives a good correlation between the modelled and the reconstructed by ground based measurements’ modulation potential. However, there appears to be a small deviation during the solar maximum of cycle 22, which is likely related to high pressure of solar wind plasma flow velocity. Such events are hard to predict and subsequently to include in the model, but their effects may be significant. The second discrepancy occurring during the maximum of solar cycle 24 can be explained in terms of the unusually extended polar field reversal, with both northern and southern polar fields being simultaneously positive for over a year, leading to a higher flux of GCR particles at Earth.
To test the validity of our model, we reconstructed the modulation potential since the 1616, which was subsequently used to compute the global radiocarbon production rate and the ^{10}Be flux. The computations show similar longterm variations with the global radionuclides production records from terrestrial archives such as tree rings and ice cores which validate the approach. The period studied includes different levels of solar activity and therefore the semiempirical model can be considered as a good approach for reconstructions in centennial and millennia scales.
Acknowledgments
This work was done in the framework of the ReSoLVE Center of Excellence (Academy of Finland, Project No. 272157). We thank M.I. Pishkalo for providing us with the observed HCS tilts from image analysis of the solar corona during total solar eclipses.
The editor thanks Charles W. Smith and one anonymous referee for their assistance in evaluating this paper.
References
 AlankoHuotari, K., K. Mursula, I.G. Usoskin, and G.A. Kovaltsov. Global heliospheric parameters and cosmicray modulation: an empirical relation for the last decades. Sol. Phys., 238, 391–404, 2006. [NASA ADS] [CrossRef] [Google Scholar]
 AlankoHuotari, K., I.G. Usoskin, K. Mursula, and G.A. Kovaltsov. Cyclic variations of the heliospheric tilt angle and cosmic ray modulation. Adv. Space Res., 40, 1064–1069, 2007. [CrossRef] [Google Scholar]
 Beer, J., K. McCracken, and R. von Steiger. Cosmogenic radionuclides. SpringerVerlag, Berlin Heidelberg, 2012. [Google Scholar]
 Belov, A. Large scale modulation: view from the Earth. Space Sci. Rev., 93, 79–105, 2000. [CrossRef] [Google Scholar]
 Belov, A.V., R.T. Gushchina, V.N. Obridko, B.D. Shelting, and V.G. Yanke. Longterm variations of galactic cosmic rays in the past and future from observations of various solar activity characteristics. J. Atmos. Sol. Terr. Phys., 68, 1161–1166, 2006. [CrossRef] [Google Scholar]
 Berggren, A.M., J. Beer, G. Possnert, A. Aldahan, P. Kubik, M. Christl, S.J. Johnsen, J. Abreu, and B.M. Vinther. A 600year annual ^{10}Be record from the NGRIP ice core, Greenland. Geophys. Res. Lett., 36, L11801, 2009. [Google Scholar]
 Burlaga, L.F., F.B. McDonald, and N.F. Ness. Cosmic ray modulation and the distant heliospheric magnetic field – Voyager 1 and 2 observations from 1986 to 1989. J. Geophys. Res., 98, 1–11, 1993. [NASA ADS] [CrossRef] [Google Scholar]
 CaballeroLopez, R.A., and H. Moraal. Limitations of the force field equation to describe cosmic ray modulation. J. Geophys. Res., 109, A01101, 2004. [Google Scholar]
 Cliver, E.W. The shapes of galactic cosmic ray intensity maxima and the evolution of the heliospheric current sheet. J. Geophys. Res., 98, 17435–17442, 1993. [CrossRef] [Google Scholar]
 Cliver, E.W., and A.G. Ling. 22 year patterns in the relationship of sunspot number and tilt angle to cosmicray intensity. Astrophys. J. Lett., 551, L189–L192, 2001. [Google Scholar]
 Cliver, E.W., I.G. Richardson, and A.G. Ling. Solar drivers of 11yr and longterm cosmic ray modulation. Space Sci. Rev., 176, 3–19, 2013. [Google Scholar]
 Dorman, L. Cosmic ray interactions, propagation, and acceleration in space plasmas, vol. 339, Springer, Dordrecht, The Netherlands, 2006. [Google Scholar]
 Ferreira, S.E.S., and M.S. Potgieter. Longterm cosmicray modulation in the heliosphere. Astrophys. J., 603, 744–752, 2004. [Google Scholar]
 Gleeson, L.J., and W.I. Axford. Solar modulation of galactic cosmic rays. Astrophys. J., 154, 1011, 1968, DOI: 10.1086/149822. [CrossRef] [Google Scholar]
 Hathaway, D.H. The solar cycle. Living Rev. Sol. Phys., 7, 1, 2010. [Google Scholar]
 Hattingh, M., and R.A. Burger. Some properties of a fully threedimensional drift model for the modulation of galactic cosmic rays. International Cosmic Ray Conference, 4, 337, 1995. [Google Scholar]
 Hoeksema, J.T. Largescale solar and heliospheric magnetic fields. Adv. Space Res., 11, 15–24, 1991. [NASA ADS] [CrossRef] [Google Scholar]
 Jokipii, J.R., and B. Thomas. Effects of drift on the transport of cosmic rays. IV – Modulation by a wavy interplanetary current sheet. Astrophys. J., 243, 1115–1122, 1981. [Google Scholar]
 Kota, J., and J.R. Jokipii. Effects of drift on the transport of cosmic rays. VI – A threedimensional model including diffusion. Astrophys. J., 265, 573–581, 1983. [CrossRef] [Google Scholar]
 Kovaltsov, G.A., and I.G. Usoskin. A new 3D numerical model of cosmogenic nuclide ^{10}Be production in the atmosphere. Earth Planet. Sci. Lett., 291, 182–188, 2010. [NASA ADS] [CrossRef] [Google Scholar]
 Kovaltsov, G.A., A. Mishev, and I.G. Usoskin. A new model of cosmogenic production of radiocarbon ^{14}C in the atmosphere. Earth Planet. Sci. Lett., 337, 114–120, 2012. [Google Scholar]
 Krymskij, G.F. Modulation of cosmic rays in interplanetary space, Nauka, Moskva, 152, 1969. [Google Scholar]
 Kta, J., and J.R. Jokipii. Cosmic ray transport in a heliospheric magnetic field with nonpolar coronal holes. Space Sci. Rev., 97, 327–330, 2001. [CrossRef] [Google Scholar]
 Lockwood, M., and M.J. Owens. Centennial variations in sunspot number, open solar flux and streamer belt width: 3. Modeling. J. Geophys. Res., 119, 5193–5209, 2014. [Google Scholar]
 Lockwood, M., L. Barnard, H. Nevanlinna, M.J. Owens, R.G. Harrison, A.P. Rouillard, and C.J. Davis. Reconstruction of geomagnetic activity and nearEarth interplanetary conditions over the past 167 yr – Part 1: a new geomagnetic data composite. Ann. Geophys., 31, 1957–1977, 2013a. [Google Scholar]
 Lockwood, M., L. Barnard, H. Nevanlinna, M.J. Owens, R.G. Harrison, A.P. Rouillard, and C.J. Davis. Reconstruction of geomagnetic activity and nearEarth interplanetary conditions over the past 167 yr – Part 2: a new reconstruction of the interplanetary magnetic field. Ann. Geophys., 31, 1979–1992, 2013b. [Google Scholar]
 Lockwood, M., H. Nevanlinna, L. Barnard, M.J. Owens, R.G. Harrison, A.P. Rouillard, and C.J. Scott. Reconstruction of geomagnetic activity and nearEarth interplanetary conditions over the past 167 yr – Part 4: nearEarth solar wind speed, IMF, and open solar flux. Ann. Geophys., 32, 383–399, 2014a. [Google Scholar]
 Lockwood, M., H. Nevanlinna, M. Vokhmyanin, D. Ponyavin, S. Sokolov, et al. Reconstruction of geomagnetic activity and nearEarth interplanetary conditions over the past 167 yr – Part 3: improved representation of solar cycle 11. Ann. Geophys., 32, 367–381, 2014b. [Google Scholar]
 Lopate, C., and J.A. Simpson. The physics of cosmic ray modulation – heliospheric propagation during the 1987 minimum. J. Geophys. Res., 96, 15877, 1991. [CrossRef] [Google Scholar]
 McCracken, K.G., and J. Beer. Comparison of the extended solar minimum of 2006–2009 with the Spoerer, Maunder, and Dalton Grand Minima in solar activity in the past. J. Geophys. Res., 119, 2379–2387, 2014. [Google Scholar]
 Mursula, K., I.G. Usoskin, and G.A. Kovaltsov. Reconstructing the longterm cosmic ray intensity: linear relations do not work. Ann. Geophys., 21, 863–867, 2003. [CrossRef] [Google Scholar]
 Parker, E.N. The passage of energetic charged particles through interplanetary space. Planet. Space Sci., 13, 9–49, 1965. [Google Scholar]
 Pedro, J., T. van Ommen, M. Curran, V. Morgan, A. Smith, and A. McMorrow. Evidence for climate modulation of the ^{10}Be solar activity proxy. J. Geophys. Res., 111, D21105, 2006. [CrossRef] [Google Scholar]
 Pedro, J.B., U.E. Heikkil, A. Klekociuk, A.M. Smith, T.D. van Ommen, and M.A.J. Curran. Beryllium10 transport to Antarctica: results from seasonally resolved observations and modeling. J. Geophys. Res., 116, D23120, 2011. [CrossRef] [Google Scholar]
 Pishkalo, M.I. Reconstruction of the heliospheric current sheet tilts using sunspot numbers. Sol. Phys., 233, 277–290, 2006. [CrossRef] [Google Scholar]
 Potgieter, M. Solar modulation of cosmic rays. Living Rev. Sol. Phys., 10, 3, 2013, DOI: 10.12942/lrsp20133. [Google Scholar]
 Potgieter, M.S., and J.A. Le Roux. The simulated features of heliospheric cosmicray modulation with a timedependent drift model. I – General effects of the changing neutral sheet over the period 1985–1990. Astrophys. J., 386, 336–346, 1992. [CrossRef] [Google Scholar]
 Potgieter, M.S., R.A. Burger, and S.E.S. Ferreira. Modulation of cosmic rays in the heliosphere from solar minimum to maximum: a theoretical perspective. Space Sci. Rev., 97, 295–307, 2001. [NASA ADS] [CrossRef] [Google Scholar]
 Roth, R., and F. Joos. A reconstruction of radiocarbon production and total solar irradiance from the Holocene ^{14}C and CO_{2} records: implications of data and model uncertainties. Climate of the Past, 9, 1879–1909, 2013. [Google Scholar]
 Sabbah, I., and M. Rybansk. Galactic cosmic ray modulation during the last five solar cycles. J. Geophys. Res., 111, A01105, 2006. [CrossRef] [Google Scholar]
 Solanki, S.K., I.G. Usoskin, B. Kromer, M. Schssler, and J. Beer. Unusual activity of the Sun during recent decades compared to the previous 11,000 years. Nature, 431, 1084–1087, 2004. [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Stozhkov, Y.I., V.P. Okhlopkov, and N.S. Svirzhevsky. Cosmic ray fluxes in present and past times. Sol. Phys., 224, 323–333, 2004. [CrossRef] [Google Scholar]
 Suess, S.T., D.J. McComas, and J.T. Hoeksema. Prediction of the heliospheric current sheet tilt – 1992–1996. Geophys. Res. Lett., 20, 161–164, 1993. [CrossRef] [Google Scholar]
 Sun, X., J.T. Hoeksema, Y. Liu, and J. Zhao. On polar magnetic field reversal and surface flux transport during solar cycle 24. Astrophys. J., 798, 114, 2015. [Google Scholar]
 Usoskin, I.G. A history of solar activity over millennia. Living Rev. Sol. Phys., 10, 1, 2013. [Google Scholar]
 Usoskin, I.G., K. Mursula, S.K. Solanki, M. Schssler, and G.A. Kovaltsov. A physical reconstruction of cosmic ray intensity since 1610. J. Geophys. Res., 107, 1374, 2002. [NASA ADS] [CrossRef] [Google Scholar]
 Usoskin, I.G., K. AlankoHuotari, G.A. Kovaltsov, and K. Mursula. Heliospheric modulation of cosmic rays: monthly reconstruction for 1951–2004. J. Geophys. Res., 110, A12108, 2005. [NASA ADS] [CrossRef] [Google Scholar]
 Usoskin, I.G., K. Horiuchi, S. Solanki, G.A. Kovaltsov, and E. Bard. On the common solar signal in different cosmogenic isotope data sets. J. Geophys. Res., 114, A03112, 2009. [Google Scholar]
 Usoskin, I.G., G.A. Bazilevskaya, and G.A. Kovaltsov. Solar modulation parameter for cosmic rays since 1936 reconstructed from groundbased neutron monitors and ionization chambers. J. Geophys. Res., 116, A02104, 2011. [NASA ADS] [CrossRef] [Google Scholar]
 Uzal, L.C., R.D. Piacentini, and P.F. Verdes. Predictions of the maximum amplitude, time of occurrence, and total length of solar cycle 24. Sol. Phys., 279, 551–560, 2012. [CrossRef] [Google Scholar]
 Wang, Y.M., and N.R. Sheeley Jr. On potential field models of the solar corona. Astrophys. J., 392, 310–319, 1992. [Google Scholar]
 Wang, Y.M., N.R. Sheeley Jr., and A.P. Rouillard. Role of the Sun’s nonaxisymmetric open flux in cosmicray modulation. Astrophys. J., 644, 638–645, 2006. [Google Scholar]
 Wibberenz, G., H.V. Cane, I.G. Richardson, and T.T. von Rosenvinge. The influence of tilt angle and magnetic field variations on cosmic ray modulation. Space Sci. Rev., 97, 343–347, 2001a. [CrossRef] [Google Scholar]
 Wibberenz, G., S.E.S. Ferreira, M.S. Potgieter, and H.V. Cane. Timedependent 2D model compared with observations during the 1974 mini cycle. Space Sci. Rev., 97, 373–376, 2001b. [CrossRef] [Google Scholar]
Cite this article as: Asvestari E & Usoskin IG. An empirical model of heliospheric cosmic ray modulation on longterm time scale. J. Space Weather Space Clim., 6, A15, 2016, DOI: 10.1051/swsc/2016011.
All Figures
Fig. 1. Annual variations of the yearly averaged sunspot numbers (Solar Influences Data Analysis Center (SIDC) at the Royal Observatory of Belgium) (version 1.0) (blue solid line) and the HCS tilt angle, α, (green dashed line) for the period 1976–2014. The correlation coefficient between the two is . 

In the text 
Fig. 2. Yearly averaged HCS tilt angle variations based on the model (Eq. (1)) and the WSO observations of the photospheric magnetic field in relation to the cycle phase, for solar cycles 21–24 as indicated in the legend. 

In the text 
Fig. 3. Temporal variations of the HCS tilt angle based on the radial model by WSO (blue line) and reconstructed here using the empirical model (Eq. (1)) (magenta steplike line). Black stars represent reconstructed HCS tilt angles by an analysis of total solar eclipse images (courtesy of M.I. Pishkalo). 

In the text 
Fig. 4. Scatter plot of the modelled here versus the WSO provided values of the annually averaged HCS tilt angle for the period 1976–2014 (as shown in Fig. 3). The blue solid line marks the diagonal. 

In the text 
Fig. 5. Annual values of the heliospheric parameters. From top to bottom are the modulation potential ϕ (MV), open solar magnetic flux F (10^{15} Wb), calculated tilt angle α (deg) and polarity p for the period 1951–2013. 

In the text 
Fig. 6. Temporal variations of the modulation potential reconstructed by ground based observations of cosmic rays (blue curve) and the modelled (magenta dashed curve) over the period 1951–2013. 

In the text 
Fig. 7. Scatter plot of the data shown in Figure 6. The solid line represents the diagonal. The correlation coefficient, R, between the two parameters is approximately 0.88 ± 0.03. 

In the text 
Fig. 8. Annual variations of radiocarbon production reconstructed by Roth & Joos (2013) (magenta curve, with 95% of confidence interval shade) and modelled here (blue curve). 

In the text 
Fig. 9. Annual variations of radionuclide ^{10}Be flux by Berggren et al. (2009) (magenta curve) and modelled here (blue curve). Although the two curves show a very good agreement after 1900, they vary from each other significantly before that period. This is likely due to atmospheric distribution and deposition processes that affect the stored amount of the radionuclide in ice. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.