The cosmic rays differential intensity inside the heliosphere, for energy below 30 GeV/nuc, depends on solar activity and interplanetary magnetic field polarity. This variation, termed solar modulation, is described using a 2D (radius and colatitude) Monte Carlo approach for solving the Parker transport equation that includes diffusion, convection, magnetic drift, and adiabatic energy loss. Since the whole transport is strongly related to the interplanetary magnetic field (IMF) structure, a better understanding of his description is needed in order to reproduce the cosmic rays intensity at the Earth, as well as outside the ecliptic plane. In this work an interplanetary magnetic field model including the standard description on ecliptic region and a polar correction is presented. This treatment of the IMF, implemented in the HelMod Monte Carlo code (version 2.0), was used to determine the effects on the differential intensity of Proton at 1 AU and allowed one to investigate how latitudinal gradients of proton intensities, observed in the inner heliosphere with the Ulysses spacecraft during 1995, can be affected by the modification of the IMF in the polar regions.
1. Introduction
The Solar Modulation, due to the solar activity, affects the Local Interstellar Spectrum (LIS) of Galactic Cosmic Rays (GCR) typically at energies lower than 30 GeV/nucl. This process, described by means of the Parker equation (e.g., see [1, 2] and Chapter 4 of [3]), is originated from the interaction of GCRs with the interplanetary magnetic field (IMF) and its irregularities. The IMF is the magnetic field that is carried outwards during the solar wind expansion. The interplanetary conditions vary as a function of the solar cycle which approximately lasts eleven years. In a solar cycle, when the maximum activity occurs, the IMF reverse his polarity. Thus, similar solar polarity conditions are found almost every 22 years [4]. In the HelMod Monte Carlo code version 1.5 (e.g., see [2]), the “classical” description of IMF, as proposed by Parker [5], was implemented together with the polar corrections of the solar magnetic field suggested subsequently in [6, 7]. This IMF was used inside the HelMod [2] code to investigate the solar modulation observed at Earth and to partially account for GCR latitudinal gradients, that is, those observed with the Ulysses spacecraft [8, 9]. In order to fully account for both the latitudinal gradients and latitudinal position of the proton-intensity minimum observed during the Ulysses fast scan in 1995, the HelMod Code was updated to the version 2.0 to include a new treatment of the parallel and diffusion coefficients following that one described in [10]. In the present formulation, the parallel component of the diffusion tensor depends only on the radial distance from the Sun, while it is independent of solar latitude.
2. The Interplanetary Magnetic Field
Nowadays, we know that there is a Solar Wind plasma (SW) that permeates the interplanetary space and constitutes the interplanetary medium. In IMF models the magnetic-field lines are supposed to be embedded in the nonrelativistic streaming particles of the SW, which carries the field with them into interplanetary space, producing the large scale structure of the IMF and the heliosphere. The “classical” description of the IMF was proposed originally by Parker (e.g., see [2, 5, 11–14] and Chapter 4 of [3]). He assumed (i) a constant solar rotation with angular velocity (ω), (ii) a simple spherically symmetric emission of the SW, and (iii) a constant (or approaching an almost constant) SW speed (Vsw) at larger radial distances (r), for example, for r>rb≈10R⊙ (where R⊙ is the Solar radius), since beyond rb the wind speed varies slowly with the distance. The “classical” IMF can be analytically expressed as [15]
(1)BPar=Ar2(er-Γeφ)[1-2H(θ-θ′)],
where A is a coefficient that determines the IMF polarity and allows |BPar| to be equal to B⊕, that is, the value of the IMF at Earth's orbit as extracted from NASA/GSFC's OMNI data set through OMNIWeb [16, 17]; er and eφ are unit vector components in the radial and azimuthal directions, respectively, θ is the colatitude (polar angle); θ′ is the polar angle determining the position of the Heliospheric Current Sheet (HCS) [18]; H is the Heaviside function: thus, [1-2H(θ-θ′)] allows BPar to change sign in the two regions above and below the HCS [18] of the heliosphere; finally,
(2)Γ=tanΨ=ω(r-rb)sinθVsw,
with Ψ the spiral angle. In the present model ω is assumed to be independent of the heliographic latitude and equal to the sidereal rotation at the Sun's equator. The magnitude of Parker field is thus
(3)BPar=Ar21+Γ2.
In 1989 [6], Jokipii and Kota have argued that the solar surface, where the feet of the field lines lie, is not a smooth surface, but a granular turbulent surface that keeps changing with time, especially in the polar regions. This turbulence may cause the footpoints of the polar field lines to wander randomly, creating transverse components in the field, thus causing temporal deviations from the smooth Parker geometry. The net effect of this is a highly irregular and compressed field line. In other words, the magnitude of the mean magnetic field at the poles is greater than in the case of the smooth magnetic field of a pure Parker spiral. Jokipii and Kota [6] have therefore suggested that the Parker spiral field may be generalized by the introduction of a perturbation parameter [δ(θ)] which amplifies the field strength at large radial distances. With this modification the magnitude of IMF, (3), becomes [6]
(4)BPol=Ar21+Γ2+(rrb)2δ(θ)2.
The difference of the IMF obtained from (3) and (4) is less than ~1% for colatitudes 20°≤θ≤160° (e.g, see Figure 1(a)) and increases for colatitudes approaching the polar regions (e.g., see Figure 2 of [2]).
(a) Maximum percentage difference between BPol and BPar as a function of the solar distance inside the colatitude regions 20°≤θ≤160° (θ20°), 25°≤θ≤155° (θ25°), and 30°≤θ≤150° (θ30°). (b) Maximum value allowed of δm as a function of colatitude: values in the “allowed” region guarantee that the stream line originated from the polar magnetic field do not cross the equatorial plane. Since (5) is symmetric with respect to the solar equatorial plane, values greater than 90° of colatitude lead to same results of those presented.
In the present treatment, the heliosphere is divided into polar regions and a equatorial region where different description of IMF are applied. In the equatorial region the Parker's IMF, Equation (3) is used, while in the polar regions we used a modified IMF that allows a magnitude as in (4)
(5)BPol=Ar2[er+rrbδ(θ)eθ-ω(r-rb)sinθVsweφ]×[1-2H(θ-θ′)]Polarregions,BPar=Ar2[er-ω(r-rb)sinθVsweφ]×[1-2H(θ-θ′)]Equatorialregion,
where equatorial regions are those with colatitude X°≤θ≤(180°-X°). The symbol θX° indicates the corresponding polar regions.
In order to have a divergence-free magnetic-field we require that the perturbation factor [δ(θ)] has to be
(6)δ(θ)=δm[1-2H(θ-θ′)]sinθ,
where δm is the minimum perturbation factor of the field. The perturbation parameter is left to grow with decreasing of the colatitude. However, in their original work, Jokipii and Kota [6] estimated the value of the parameter δ between 10-3 and 3×10-3.
Since the polar field is only a perturbation of the Parker field, it is a reasonable assumption that stream lines of the magnetic field do not cross the equatorial plane, thus, remaining completely contained in the solar hemisphere of injection. This allows one to estimate an upper limit on the possible values of δm (see, Figure 1(b)). Currently, we use δm=1×10-5 by comparing simulations with observations at Earth orbit during Solar Cycle 23 (see, Section 5).
3. The Propagation Model
Parker in 1965 [1] (see also, [2], Chapter 4 of [3] and references therein) treated the propagation of GCRs trough the interplanetary space. He accounted for the so-called adiabatic energy losses, outward convection due to the SW and drift effects. In the heliocentric system the Parker equation is then expressed (e.g., see [2, 19]):
(7)∂U∂t=∂∂xi(KijS∂U∂xj)-∂∂xi[(Vsw,i+vd,i)U]+13∂Vsw,i∂xi∂∂T(αrelTU),
where U is the number density of particles per unit of particle kinetic energy T, at the time t. Vsw,i is the solar wind velocity along the axis xi [20], KijS is the symmetric part of diffusion tensor [1], vd is the drift velocity that takes into account the drift of the particles due to the large scale structure of the magnetic field [21–23], and, finally,
(8)αrel=T+2mrc2T+mrc2,
where mr is the rest mass of the GCR particle. The last term of (7) accounts for adiabatic energy losses [1, 24]. The number density U is related to the differential intensity J as ([2, 25], Chapter 4 of [3] and references therein):
(9)J=vU4π,
where v is the speed of the GCR particle.
Equation (7) was solved using the HelMod code (see the discussion in [2]). This treatment (i) follows that introduced in [10, 26–30] and (ii) determines the differential intensity of GCRs using a set of approximated stochastic differential equations (SDEs) which provides a solution equivalent to that from (7). The equivalence between the Parker equation, that is a Fokker-Planck type equation, and the SDEs is demonstrated in [28, 31]. In the present work, we use a 2D (radius and colatitude) approximation for the particle transport. The model includes the effects of solar activity during the propagation from the effective boundary of the heliosphere down to Earth's position.
The set of SDEs for the 2D approximation of (7) in heliocentric spherical coordinates is(10a)Δr=1r2∂∂r(r2KrrS)Δt-∂∂μ(θ)[KrμS1-μ2(θ)r]Δt+(Vsw+vd,r)Δt+(2KrrS)1/2ωrΔt,(10b)Δμ(θ)=-1r2∂∂r[rKμrS1-μ2(θ)]Δt+∂∂μ(θ)[KμμS1-μ2(θ)r2]Δt-1rvd,μ1-μ2(θ)Δt-2KrμSr[1-μ2(θ)2KrrS]1/2ωrΔt+1r{[1-μ2(θ)]KμμSKrrS-(KrμS)20.5KrrS}1/2ωμΔt,(10c)ΔT=-αrelT3r2∂Vswr2∂rΔt,where μ(θ)=cosθ and ωi is a random number following a Gaussian distribution with a mean of zero and a standard deviation of one. The procedure for determining the SDEs can be found in [2].
In a coordinate system with one axis parallel to the average magnetic field and the other two perpendicular to this the symmetric part of the diffusion tensor KijS is (see e.g., [32]):
(11)KijS=[K||000K⊥,r000K⊥,θ],
with K|| the diffusion coefficient describing the diffusion parallel to the average magnetic field and K⊥,r and K⊥,θ are the diffusion coefficients describing the diffusion perpendicular to the average magnetic field in the radial and polar directions, respectively. In this work K|| is that one proposed by Strauss and collaborators in [10] (see also [33–35]):
(12)K||=β3K0P1GV(1+r1AU),
where K0 is the diffusion parameter—described in Section 2.1 of [2]—which depends on solar activity and polarity, β is the particle speed in unit of speed of light, P=pc/|Z|e is the particle rigidity expressed in GV, and, finally, r is the heliocentric distance from the Sun in AU.
In the current treatment, K|| has a radial dependence proportional to r, but no latitudinal dependence. Mc Donald and collaborators (see, [36]) remarked that (i) a spatial dependence of K||—like the one proposed here—can affect the latitudinal gradients at high latitude and (ii) it is consistent with that originally suggested in [6]. Furthermore, the perpendicular diffusion coefficient is taken to be proportional to K|| with a ratio K⊥,i/K||=0.13 for both r and θi-coordinates. The latter value is discussed in Section 5.
In addition, the practical relationship between K0 and monthly Smoothed Sunspot Numbers (SSN) [37] values—discussed in Section 2.1 of [2]—is currently updated using the most recent data from [38]. As in [2], the K0 data are subdivided into four sets, that is, ascending and descending phases for both negative and positive solar magnetic-field polarities (Table 1). It has to be remarked that after each maximum the sign of the magnetic field (i.e., the A parameter in (5)) is reversed. The updated practical relationships between K0 in AU2GV-1s-1 and SSN values for 1.4≤SSN≤165 for the four periods (from I up to IV, listed in Table 1) are (13a)(I)K0=0.000297-2.9·10-6SSN+8.1·10-9SSN2+1.46·10-10SSN3-8.4·10-13SSN4,(13b)(II)K0={0.000304846-5.8·10-6SSN,ifSSN≤20,0.00195SSN-2.3·10-10SSN2+9.1·10-5,ifSSN>20,(13c)(III)K0=0.0002391-8.453·10-7SSN(13d)(IV)K0=0.000247-1.175·10-6SSN.The rms (root mean square) values of the percentage difference between values obtained with (13a), (13b), (13c), and (13d) from those determined using the procedure discussed in Section 2.1 of [2] applied to the data from [38] were found to be 6.0%, 10.1%, 7.0%, and 13.2% for the period I (ascending phase with A<0), II (descending phase with A<0), III (ascending phase with A>0), and IV (descending phase with A>0), respectively.
Definition of ascending and descending periods.
Period
Years
(I)
A<0 ascending
1964.79–1968.87, 1986.70–1989.54, 2008.95–2009.95
(II)
A<0 descending
1964.53–1964.79, 1979.95–1986.70, 2000.28–2008.95
(III)
A>0 ascending
1976.20–1979.95, 1996.37–2000.28
(IV)
A>0 descending
1968.87–1976.20, 1989.54–1996.37
4. The Magnetic Field in the Polar Regions
Section 2 describes an IMF following the Parker Field with a small region around the poles in which such a field is modified. As already mentioned (see e.g., [6, 39–43]), the correction is needed to better reproduce the complexity of the magnetic field in those regions.
Moreover, Ulysses spacecraft (see e.g., [44–46]) explored the heliosphere outside the ecliptic plane up to ±80° of solar latitude at a solar distance from ~1 up to ~5 AU. Using these observations, the presence of latitudinal gradient in the proton intensity could be determined (e.g., see Figure 2 of [9] and Figure 5 of [8]). The data collected during the latitudinal fast scan (from September 1994 up to August 1995) show (a) a nearly symmetric latitudinal gradient with the minimum near ecliptic plane, (b) a southward shift of the minimum, and (c) an intensity in the North polar region at 80° exceeding the South polar intensity. In [9] a latitudinal gradient of ~0.3%/degree for proton with kinetic energy >0.1 GeV was estimated. While in [8] the analysis to higher energy was extended estimating a gradient of ~0.22%/degree for proton with kinetic energy >2 GeV. The minimum in the charged particle intensity separating the two hemispheres of the heliosphere occurs ~10° South of the heliographic equator [9]. In addition, an independent analysis that takes into account the latitudinal motion of the Earth and IMP8 confirms a significant (~8°±2°) southward offset of the intensity minimum [9] for T>100 MeV proton. Furthermore in [8], a southward offset of about ≈7° is evaluated; this offset of the intensity minimum results to be independent of the particle energy up to 2 GeV. Finally, in [9], the intensity in the North polar region at 80 is observed to exceed the South polar intensity of ~6% for protons with T>100 MeV.
Using the present HelMod code (version 2.0), we could investigate (i) the latitudinal gradient of GCR intensities resulting from solar modulation and (ii) how the magnetic-field structure of the polar regions, as defined in Section 2, is able to influence the GCR spectra on the ecliptic plane. As previously defined, we denote with θX° a polar region of amplitude X° from polar axis, that is, θ<X° and θ>180°-X°. Three regions with X°=20°,25°, and 30° were investigated. Outside any of these regions, the ratio between BPol and BPar in (5) is less than ~1% (see Figure 1) and, thus, it ensures a smooth transition between polar and equatorial regions. For the purpose or this study we consider an energy binning closer to those presented in [8]. The KET instrument [47] collects proton data in three “channels” one with energies ranging from 0.038 GeV up to 2.0 GeV and two for particle with kinetic energy T>0.1 GeV and T>2 GeV, respectively. A successive reanalysis of the collected data allowed the authors to subdivide the 0.25–2 GeV “channel” in three “subchannels” of intermediate energies. Since the Present Model is optimized—as discussed in [2]—for particles with rigidity greater than 1 GV (i.e., ≈0.444 GeV), the present results are compared only with the corresponding “channel” or “subchannels” suited for the corresponding energy range (see Table 2).
Kinetic energy bins (in GeV) selected with HelMod Code and the corresponding proton energy channel for KET instruments on board the Ulysses spacecraft [8].
HelMod Energy bin
KET channel
(0.35–0.98) GeV
(0.4–1.0) GeV
(0.76–2.09) GeV
(0.8–2.0) GeV
(2.09–200) GeV
>2 GeV
At 1 AU and as a function of the solar colatitude, the GCR intensities for protons are shown in Figure 2. For a comparison with Ulysses observations, the modulated intensities of protons—resulting from HelMod code—were investigated from 80° (North) and down to -80° (South). They were obtained using the HelMod code and selected using the energy bins reported in Table 2. In Figure 2, the latitudinal intensity distribution is normalized to the corresponding South Pole intensity. The quoted errors include statistical and systematic errors. The distributions were interpolated using a parabolic function expressed as
(14)I(θlat)=a+c(θlat+d)2,
where I(θlat) is the normalized intensity, θlat is the latitudinal angle (the latitudinal angle is θlat=90°-θ), and a, b, and c are parameters determined from the fitting procedure. The so obtained fitted curves are shown as continuous lines in Figure 2. Furthermore, the latitudinal positions of minimum intensity (θlat,min), percentages of North-South asymmetry of intensities (ΔN-S), and differences in percentage between the maximum and minimum intensities (Δmax) were also determined from the fitting procedure and are listed in Tables 3, 4, and 5, respectively. The quoted errors—following a procedure discussed in [2]—are derived varying the fitted parameters in order to obtain a value of the parameter ηrms,lat two times larger than the one resulting from the best fit (Ibest). ηrms,lat is defined as
(15)ηrms.lat=∑i(ηi,lat/ση,i,lat)2∑i1/ση,i,lat2
with
(16)ηi,lat=I(θi)-Ibest(θi)Ibest(θi),
where θi is the central value of the ith latitudinal bin of the differential intensity distribution and ση,i,lat are the errors including the experimental and Monte Carlo uncertainties. For θ30°, that is, assuming a modified polar magnetic-field for θ<30° and θ>150°, we found a general agreement with Ulysses observations. The position of θlat,min is compatible within the errors with one observed in [8], as well as the values of ΔN-S and Δmax.
Latitudinal positions of minimum proton intensity (θlat,min) (in degrees) as a function of the kinetic energy, using three values for the extension of the polar regions.
Energy range
θlat,min (in degrees)
θ20∘
θ25∘
θ30∘
0.35–0.98 GeV
-5-5+4
-8-4+4
-9-3+3
0.76–2.09 GeV
-5-7+6
-6-6+5
-10-4+4
2.09–200 GeV
-3-7+7
-7-7+7
-9-6+5
Percentages of North-South asymmetry of proton intensities (ΔN-S) as a function of the kinetic energy, using three values for the extension of the polar regions.
Energy range
ΔN-S (%)
θ20∘
θ25∘
θ30∘
0.35–0.98 GeV
-8-6+6
-11-4+5
-13-4+4
0.76–2.09 GeV
-5-6+5
-6-5+5
-9-3+3
2.09–200 GeV
-1-1+3
-3-3+3
-4-2+2
Differences in percentage between the maximum and minimum proton intensities (Δmax) as a function of the kinetic energy, using three values for the extension of the polar regions.
Energy range
Δmax (%)
θ20∘
θ25∘
θ30∘
0.35–0.98 GeV
-34-5+5
-35-4+4
-36-3+4
0.76–2.09 GeV
-22-7+6
-23-4+5
-25-3+3
2.09–200 GeV
-8-3+4
-9-3+3
-10-2+2
Latitudinal relative intensity at r=1AU, obtained at different solar colatitudes for protons in the energy range defined in the Table 2 and using three definitions of polar regions: θ20°, θ25°, and θ30°.
The HelMod code allows one to investigate the relevance of the treatment of the polar region magnetic-field with respect to the resulting modulation of GCR. Thus, the latitudinal normalized intensities were obtained excluding a few (small) regions nearby the poles. This was determined from reducing the latitudinal spatial phase-space admissible for pseudo-particles (see [2]), that is, the latitudinal extension of GCR particles taken into account—to 1°<θ<179° (θR,1), 2°<θ<178° (θR,2) and 10°<θ<170° (θR,10). The so obtained latitudinal gradients are compared with the full latitudinal extension, θR,0 (0°<θ<180°), in Figure 3. By an inspection of Figure 3, one may lead to the conclusion that the GCR diffusion nearby the polar axis has a large impact on the latitudinal gradients in the inner heliosphere. As a consequence, the IMF description in the polar regions is relevant in order to reproduce the observed modulated GCR spectra.
For the polar regions θ<30° and θ>150°, latitudinal relative intensity at r=1AU, accounting for GCR particles with 0°<θ<180° (θR,0), 1°<θ<179° (θR,1), and 2°<θ<178° (θR,2), respectively, are shown as a function of the proton kinetic energy and solar colatitude. Results with 10°<θ<170° (θR,10) are comparable with those obtained with θR,2.
5. Comparison with Observations during Solar Cycle 23
The agreement of HelMod simulated spectra with observations during solar cycle 23 is investigated via quantitative comparisons using (17) and (18). However, since the structure of the heliosphere is different in high and low solar activity the two periods are separately analyzed.
The HelMod Code [2] (version 2.0) allowed us to investigate how the modulated (simulated) differential intensities are affected by the (1) particle drift effect, (2) polar enhancement of the diffusion tensor along the polar direction (K⊥,θ, e.g., see [48]), and, finally, (3) values of the tilt angle (αt) calculated following the approach of the “R” and “L” models [49]. This analysis also allows us to estimate the values of IMF parameters that better describe the modulation along the entire solar cycle. The effects related to particle drift were investigated via the suppression of the drift velocity (No Drift approximation), this accounts for the hypothesis that magnetic drift convection is almost completely suppressed during solar maxima. The differential intensities were calculated for K⊥,μ=ρEK⊥,r with values of ρE of 1, 8, and 10, that is, no enhancement, that suggested in [50] and that suggested in [2, 48] (and reference there in), respectively. Furthermore, the modulated proton spectra were derived from a LIS whose normalization constant depends on the experimental set of data and were already discussed in [2]). In addition, the differential intensities were calculated accounting for particles inside heliospheric regions where solar latitudes are lower than |5.7°|.
During the period of high solar activity for the solar cycle 23, the BESS collaboration took data in the years 1999, 2000, and 2002 (see sets of data in [51]). For period not dominated by high solar activity in solar cycle 23, BESS, AMS, and PAMELA collaborations took data, that is, BESS-1997 [51], AMS-1998 [52], and PAMELA-2006/08 [53].
Following the procedure described in [2], the observation data were compared with those obtained from HelMod code using the error-weighted root mean square (ηrms) of the relative difference (η) between experimental data (fexp) and those resulting from simulated differential intensities (fsim). For each set of experimental data and with the approximations and/or models described above, we determined the quantity
(17)ηrms=∑i(ηi/ση,i)2∑i1/ση,i2,
with
(18)ηi=fsim(Ti)-fexp(Ti)fref(Ti),
where Ti is the average energy of the ith energy bin of the differential intensity distribution and ση,i are the errors including the experimental and Monte Carlo uncertainties; the latter account for the Poisson error of each energy bin. The simulated differential intensities are interpolated with a cubic spline function. The modulation results are studied varying the parameters δm—from 0, that is, the nonmodified Parker IMF, up to 3×10-5, see Figure 1(b)—, K⊥,r/K||=ρk (from 0.10 up to 0.14) and K⊥,θ/K⊥,r=ρE (from 1 up to 10) seeking a set of parameters set that minimize ηrms. In Table 6, the average values of ηrms (in percentage, %) for low solar activity periods are listed. They were obtained in the energy range (Above 30 GeV, the differential intensity is marginally (if at all) affected by modulation.) from 444 MeV up to 30 GeV using the “L” and “R” models for the tilt angle αt and for the No Drift approximation and without any enhancement of the diffusion tensor along the polar direction (K⊥,μ). The results derived with the enhancement of the diffusion tensor along the polar direction indicate that for ρE=8 and 10 one obtains a value of ηrms that is from 1.5 up to 3 times larger with respect the case without enhancement (For a comparison, the scalar approximation presented in [2], that is, assuming that the diffusion propagation is independent of magnetic structure, leads to and average ηrms of ~15%). From inspection of Table 6, one can remark that the drift mechanism leads to a better agreement with experimental data. Furthermore the “R” and “L” models for tilt angles are comparable within the precision of the method (discussed in [2]). The minimum difference with respect to the experimental data occurs when ρk=0.11–0.13 and δm=1.0×10-5 for both “R” and “L” models, with the “L” model slightly preferred to “R”.
Average values (last three columns) of ηrms (in percentage, %) as a function of δm (×10-5) and ρk, for BESS-1997, AMS-1998, PAMELA-2006/08, obtained from (17) without enhancement of the diffusion tensor along the polar direction (ρE=1), using “R” and “L” models for the tilt angle and No Drift approximation. The differential intensities were calculated accounting for particles inside the heliospheric regions for which solar latitudes are lower than |5.7∘|.
δm(×10-5)
ρk
“R” model
“L” model
No drift
0.0
0.10
14.1
11.0
33.4
1.0
0.10
11.7
8.7
33.2
2.0
0.10
11.6
8.3
33.7
3.0
0.10
11.6
8.3
33.7
1.0
0.11
6.4
9.0
27.7
2.0
0.11
7.8
9.0
28.3
3.0
0.11
7.5
8.8
29.2
1.0
0.12
6.3
7.1
23.5
2.0
0.12
6.3
7.3
24.7
3.0
0.12
7.1
6.9
24.4
1.0
0.13
6.3
6.4
20.1
2.0
0.13
6.6
7.6
20.4
3.0
0.13
6.7
7.7
20.5
1.0
0.14
7.3
7.0
15.9
2.0
0.14
7.3
7.2
16.4
3.0
0.14
7.2
6.5
16.8
In Figures 4, 5, and 6, the differential intensities determined with the HelMod code are shown and compared to the experimental data of BESS-1997, AMS-1998, and PAMELA-2006/08, respectively; in these figures, the dashed lines are the LIS as discussed in [2]. These modulated intensities are the ones calculated for a heliospheric region where solar latitudes are lower than |5.7°|, using ρk=0.13, δm=1.0×10-5 and ρE=1 with the “L” model. Finally, one can remark that the present code, combining diffusion and drift mechanisms, is also suited to describe the modulation effect in periods when the solar activity is no longer at the maximum.
Proton differential intensity determined with the HelMod code (continuous line) compared to the experimental data of BESS-1997; the dashed line is the LIS (see the text).
Proton differential intensity determined with the HelMod code (continuous line) compared to the experimental data of AMS-1998; the dashed line is the LIS (see the text).
Proton differential intensity determined with the HelMod code (continuous line) compared to the experimental data of PAMELA 2006/08; the dashed line is the LIS (see the text).
In Table 7 we present the averages ηrms (in percentage, %) during the periods dominated by high solar activity. The simulated differential intensities were obtained for a heliospheric region where solar latitudes are lower than |5.7°| without any enhancement of the diffusion tensor along the polar direction (K⊥,μ). The simulations with ρE= 8 and 10 lead to ηrms comparable with those presented in Table 7. However, using ρE=1 provides a better agreement to experimental data at lower energy. From inspection of Table 7, one can note that “R” and “L” models for tilt angles yield comparable results within the precision of the method. Furthermore the minimum difference with the experimental data occurs when ρk=0.10 and δm=(2.0-3.0)×10-5 with the “L” model slightly preferred to “R”. The HelMod parameter configuration, which minimizes the difference to the experimental data, are reported in Table 8: one may remark that the No Drift approximation is (almost) comparable to a drift treatment for both BESS-2000 and BESS-2002 with data collected during and after the maximum of the solar activity. Apparently, the drift treatment is needed in order to describe BESS-1999 with data taken during a period approaching the solar maximum.
Average values (last three columns) of ηrms (in percentage, %) as a function of δm(×10-5) and ρk, for BESS-1999, BESS-2000, BESS-2002 obtained from (17) without enhancement of the diffusion tensor along the polar direction (ρE=1), using “R” and “L” models for the tilt angle and No Drift approximation. The differential intensities were calculated accounting for particles inside the heliospheric regions for which solar latitudes are lower than |5.7∘|.
δm(×10-5)
ρk
“R” model
“L” model
No drift
0.0
0.10
11.2
10.8
15.4
1.0
0.10
11.0
10.1
15.8
2.0
0.10
9.6
10.0
16.7
3.0
0.10
9.6
10.0
16.7
1.0
0.11
13.4
13.1
16.0
2.0
0.11
12.7
12.9
15.4
3.0
0.11
12.7
12.5
16.2
1.0
0.12
18.7
17.7
13.4
2.0
0.12
18.3
16.9
12.8
3.0
0.12
18.1
17.3
12.8
1.0
0.13
23.3
23.5
14.3
2.0
0.13
25.0
24.7
13.3
3.0
0.13
24.3
24.2
13.1
1.0
0.14
32.3
30.7
18.0
2.0
0.14
32.8
30.8
17.1
3.0
0.14
31.5
30.7
17.9
Average ηrms (in percentage, %), for BESS-1999, BESS-2000, and BESS-2002 obtained from (17) without enhancement of the diffusion tensor along the polar direction (ρE=1), δm=2.0×10-5, ρk=0.10, and using “R” and “L” models for the tilt angle and No Drift approximation. The differential intensities were calculated accounting for particles inside the heliospheric regions for which solar latitudes are lower than |5.7∘|.
Observations
“R” model
“L” model
No drift
BESS-1999
9.3
10.6
25.7
BESS-2000
12.5
12.6
16.7
BESS-2002
6.9
6.7
7.7
In Figures 7, 8, and 9, the differential intensities determined with the HelMod code are shown and compared with the experimental data of BESS-1999, BESS-2000, and BESS-2002, respectively; in these figures, the dashed lines are the LIS as discussed in [2]. These modulated intensities are the ones calculated for a heliospheric region where solar latitudes are lower than |5.7°|, using K⊥,μ=K⊥,r independently of the latitude and including particle drift effects with the values of the tilt angle from the “L” model. Finally, it is concluded that the present code combining diffusion and drift mechanisms is suited to describe the modulation effect in periods with high solar activity [2, 50, 54].
Proton differential intensity determined with the HelMod code (continuous line) compared to the experimental data of BESS-1999; the dashed line is the LIS (see the text).
Proton differential intensity determined with the HelMod code (continuous line) compared to the experimental data of BESS-2000; the dashed line is the LIS (see the text).
Proton differential intensity determined with the HelMod code (continuous line) compared to the experimental data of BESS-2002; the dashed line is the LIS (see the text).
6. Conclusion
In this work an IMF, which combines the Parker Field and its polar modification, is presented. In the polar regions, the Parker IMF was modified with an additional latitudinal components according to those proposed by Jokipii and Kota in [6]. We found the maximum perturbed value with this component yielding, as a physical result, streaming lines completely confined in the solar hemisphere of injection.
The proposed IMF is, then, used within the HelMod Monte Carlo code to determine the effects on the differential intensity of protons at 1 AU as a function of the extension of polar region, in which the modified magnetic-field is employed. We found that a polar region contained within 30° of colatitude is that one ensuring a very smooth transition to the equatorial region and allows to reproduce qualitatively and quantitatively the latitudinal profile of the GCR intensity, and the latitudinal dip shift with respect to the ecliptic plane. Finally we determined how the polar region diffusion is mostly responsible of the proton intensity latitudinal gradient observed in the inner heliosphere with the Ulysses spacecraft during 1995.
Acknowledgments
K. Kudela wishes to acknowledge VEGA Grant agency project 2/0081/10 for support. Finally, the authors acknowledge the use of NASA/GSFC's Space Physics Data Facility's OMNIWeb service and OMNI data.
ParkerE. N.The passage of energetic charged particles through interplanetary spaceBobikP.BoellaG.BoschiniM. J.ConsolandiC.Della TorreS.GervasiM.GrandiD.KudelaK.PensottiS.RancoitaP. G.TM. Tacconi. .Systematic investigation of solar modulation of galactic protons for solar cycle 23 using a monte carlo approach with particle drift effects and latitudinal dependenceheLeroyC.RancoitaP. G.StraussR. D.PotgieterM. S.FerreiraS. E. S.Modeling ground and space based cosmic ray observationsParkerE. N.Dynamics of the interplanetary gas and magnetic fieldsJokipiiJ. R.KotaJ.The polar heliospheric magnetic fieldLangnerU. W.HeberB.DrögeW.FerrandoP.HaasbroekL. J.KunowH.Müller-MellinR.PaizisC.PotgieterM. S.RaviartA.WibberenzG.Spatial variation of >40 MeV/n nuclei fluxes observed during the Ulysses rapid latitude scanSimpsonJ. A.Ulysses cosmic-ray investigations extending from the south to the north polar regions of the Sun and heliosphereStraussR. D.PotgieterM. S.BüschingI.KoppA.Modeling the modulation of galactic and jovian electrons by stochastic processesParkerE. N.Newtonian development of the dynamical properties of ionized gases of low densityParkerE. N.The hydrodynamic theory of solar corpuscular radiation and stellar windsParkerE. N.Sudden expansion of the corona following a large solar flare and the attendant magnetic field and cosmic-ray effectsParkerE. N.HattinghM.BurgerR. A.A new simulated wavy neutral sheet drift modelKingJ. H.PapitashviliN. E.Solar wind spatial scales in and comparisons of hourly Wind and ACE plasma and magnetic field dataNASA-OMNIwebonline database, 2012, http://omniweb.gsfc.nasa.gov/form/dx1.htmlJokipiiJ. R.ThomasB.Effects of drift on the transport of cosmic rays. IV—modulation by a wavy interplanetary current sheetJokipiiJ. R.LevyE. H.HubbardW. B.Effect of particle drift on cosmic-ray transport. I. general properties, application to solar modulationMarschE.AxfordW. I.McKenzieJ. F.DwivediB. N.Solar windJokipiiJ. R.LevyE. H.Effects of particle drifts on the solar modulation of galactic cosmic raysPotgieterM. S.MoraalH.A drift model for the modulation of galactic cosmic raysBurgerR. A.HattinghM.Steady-state drift-dominated modulation models for galactic cosmic raysJokipiiJ. R.E. N. Parker.On the convection, diffusion, and adiabatic deceleration of cosmic rays in the solar windGleesonL. J.AxfordW. I.Solar modulation of galactic cosmic raysYamadaY.YanagitaS.YoshidaT.A stochastic view of the solar modulation phenomena of cosmic raysGervasiM.RancoitaP. G.UsoskinI. G.KovaltsovG. A.Monte-Carlo approach to Galactic cosmic ray propagation in the heliosphereZhangM.A Markov stochastic process theory of cosmic-ray modulationAlanko-HuotariK.UsoskinI. G.MursulaK.KovaltsovG. A.Stochastic simulation of cosmic ray modulation including a wavy heliospheric current sheetPeiC.BieberJ. W.BurgerR. A.ClemJ.A general time-dependent stochastic method for solving Parker's transport equation in spherical coordinatesGardinerC. W.JokipiiJ. R.Propagation of cosmic rays in the solar windPalmerI. D.Transport coefficients of low-energy cosmic rays in interplanetaryPotgieterM. S.FerreiraS. E. S.Effects of the solar wind termination shock on the modulation of Jovian and galactic electrons in the heliosphereDrögeW.Probing heliospheric diffusion coefficients with solar energetic particlesMcDonaldF. B.FerrandoP.HeberB.KunowH.McGuireR.Müller-MellinR.PaizisC.RaviartA.WibberenzG.A comparative study of cosmic ray radial and latitudinal gradients in the inner and outer heliosphereSIDC-teamThe International Sunspot NumberMonthly Report on the International Sunspot Number, online catalogue, 1964–2010UsoskinI. G.BazilevskayaG. A.KovaltsovG. A.Solar modulation parameter for cosmic rays since 1936 reconstructed from ground-based neutron monitors and ionization chambersBaloghA.SmithE. J.TsurutaniB. T.SouthwoodD. J.ForsythR. J.HorburyT. S.The heliospheric magnetic field over the south polar region of the sunMoraalH.Proton modulation near solar minimim periods in consecutive solar cycles6Proceedings of the International Cosmic Ray Conference1990140SmithC. W.BieberJ. W.Solar cycle variation of the interplanetary magnetic field spiralFiskL. A.Motion of the footpoints of heliospheric magnetic field lines at the sun: implications for recurrent energetic particle events at high heliographic latitudesHitgeM.BurgerR. A.Cosmic ray modulation with a Fisk-type heliospheric magnetic field and a latitude-dependent solar wind speedSandersonT. R.MarsdenR. G.WenzelK. P.BaloghA.ForsythR. J.GoldsteinB. E.High-latitude observations of energetic ions during the first Ulysses polar passMarsdenR. G.The 3-D heliosphere at solar maximumBaloghA.MarsdenR. G.SmithE. J.HeberB.PotgieterM. S.FerrandoP.Solar modulation of galactic cosmic rays: the 3D heliospherePotgieterM. S.Heliospheric modulation of cosmic ray protons: role of enhanced perpendicular diffusion during periods of minimum solar modulationHoeksemaJ. T.The large-scale structure of the heliospheric current sheet during the Ulysses epochFerreiraS. E. S.PotgieterM. S.Long-term cosmic-ray modulation in the heliosphereShikazeY.OritoS.HainoS.AbeK.FukeH.HamsT.KimK. C.MakidaY.MatsudaS.MitchellJ. W.MoiseevA. A.NishimuraJ.NozakiM.OrmesJ. F.SanukiT.SasakiM.SeoE. S.StreitmatterR. E.SuzukiJ.TanakaK.YamagamiT.YamamotoA.YoshidaT.YoshimuraK.Measurements of 0.2–20 GeV/n cosmic-ray proton and helium spectra from 1997 through 2002 with the BESS spectrometerAguilarM.AlcarazJ.AllabyJ.The Alpha Magnetic Spectrometer (AMS) on the International Space Station: part I—results from the test flight on the space shuttleAdrianiO.BarbarinoG. C.BazilevskayaG. A.BellottiR.BoezioM.BogomolovE. A.BonechiL.BongiM.BonviciniV.BorisovS.BottaiS.BrunoA.CafagnaF.CampanaD.CarboneR.CarlsonP.CasolinoM.CastelliniG.ConsiglioL.De PascaleM. P.De SantisC.De SimoneN.Di FeliceV.GalperA. M.GillardW.GrishantsevaL.JerseG.KarelinA. V.KoldashovS. V.KrutkovS. Y.KvashninA. N.LeonovA.MalakhovV.MalvezziV.MarcelliL.MayorovA. G.MennW.MikhailovV. V.MocchiuttiE.MonacoA.MoriN.NikonovN.OsteriaG.PalmaF.PapiniP.PearceM.PicozzaP.PizzolottoC.RicciM.RicciariniS. B.RossettoL.SarkarR.SimonM.SparvoliR.SpillantiniP.StozhkovY. I.VacchiA.VannucciniE.VasilyevG.VoronovS. A.YurkinY. T.WuJ.ZampaG.ZampaN.ZverevV. G.PAMELA measurements of cosmic-ray proton and helium spectraNdiitwaniD. C.FerreiraS. E. S.PotgieterM. S.HeberB.Modelling cosmic ray intensities along the Ulysses trajectory