Primary Cosmic Ray Spectra in Primary Cosmic Ray Spectra in the Planet Atmospheres the Planet Atmospheres Marusja Buchvarova Marusja Buchvarova 1 , Peter Velinov , Peter Velinov 2 (1) Space Research Institute – Bulgarian Academy of Sciences, BULGARIA (2) Solar-Terrestrial Influences Laboratory (STIL), Bulgarian Academy of Sciences, BULGARIA E-mail: [email protected]
18
Embed
Primary Cosmic Ray Spectra in the Planet Atmospheres Marusja Buchvarova 1, Peter Velinov 2 (1) Space Research Institute – Bulgarian Academy of Sciences,
PARTICLE POPULATIONS ARE DIVERSE: 1) Galactic CR (GCR), with energy range from ~ 10s MeV to 10 6 GeV, which are accelerated in the space of our Galaxy. 2) Metagalactic CR with energies 10 6 GeV GeV, accelerated in the Metagalaxy. 3) Solar CR, with energy range from ~ 10s of MeV to 100s of MeV, accelerated on the Sun 4) Anomalous CR, with energy from 1 MeV ~ 100 MeV, accelerated in the interplanetary space. CR are important for the Sun-Earth connections, because they possess maximal penetration capability in comparison with the other radiations. CR maintain ionization in ionosphere, atmosphere, hydrosphere, cryosphere and lithosphere of the Earth and planets. INTRODUCTION
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Primary Cosmic Ray Spectra in the Primary Cosmic Ray Spectra in the Planet Atmospheres Planet Atmospheres
Marusja Buchvarova Marusja Buchvarova 11, Peter Velinov, Peter Velinov22
(1) Space Research Institute – Bulgarian Academy of Sciences, BULGARIA
(2) Solar-Terrestrial Influences Laboratory (STIL), Bulgarian Academy of Sciences, BULGARIA
INTRODUCTION MODELING COSMIC RAY DIFFERENTIAL SPECTRUM MODELING COSMIC RAY INTEGRAL SPECTRUM IN OUTER PLANETS COMPARISON OF THE MODELING COSMIC RAY SPECTRUM WITH THEORETICAL MODELS: • THE FORCE FIELD APPROXIMATION• 2D SOLUTIONS OF COSMIC RAY TRANSPORT EQUATION
APPLICATION FOR CALCULATION OF THE ELECTRON PRODUCTION RATE q(h) PROFILES DUE TO PARTICLES WITH ENERGY INTERVALS OF GALACTIC CR AND ANOMALOUS CR COMPONENT
OUTLINE OF THE TALK:
PARTICLE POPULATIONS ARE DIVERSE:
1) Galactic CR (GCR), with energy range from ~ 10s MeV to 106 GeV, which are accelerated in the space of our Galaxy.
2) Metagalactic CR with energies 106 GeV - 1012 GeV, accelerated in the Metagalaxy.
3) Solar CR, with energy range from ~ 10s of MeV to 100s of MeV, accelerated on the Sun
4) Anomalous CR, with energy from 1 MeV ~ 100 MeV,accelerated in the interplanetary space.
CR are important for the Sun-Earth connections, because they possess maximal penetration capability in comparison with the other radiations. CR maintain ionization in
ionosphere, atmosphere, hydrosphere, cryosphere and lithosphere of the Earth and planets.
INTRODUCTION
The intensity of the GCR with E< 20 GeV shows an inverse relationship to the 11 – years solar cycle.
Galactic cosmic ray particles - anti-correlated with SA
Solar Activity
Solar particles – correlated with SA
INTRODUCTION
This picture is from F. Nichitiu, CANADA
The expression for the differential spectrum (energy range E from ~ 30 MeV to 100 GeV) of the protons and other groups of cosmic ray nuclei on account of the anomalous cosmic rays (energy range E from 1 MeV to about 30 MeV) is (Velinov, 2000; Buchvarova, 2006):
Here the first term presents the galactic CR, and the last term takes into account anomalous CR (Cummings and Stone, 1987). The members with tanh are smoothing functions (Velinov, 2002). Here we take Kp=25.298 GeV2.75/(s.m2ster.MeV) and p = 2.75 for protons. The used parameters for the alpha particles are K = 1.145 GeV2.68/(s.m2 ster.MeV) and =2.68 . The normalization constants Kp and K are chosen to match the modulated data near to 100 GeV/ /nucl, where the modulation effect is negligible. The parameters α, β, x and y are related with modulation levels in corresponding energy intervals. The dimensionless parameter = 100 is inversely proportional to the length of the smoothing interval between the two addends. The physical meaning of (GeV) is the energy at which the differential spectrum of GCR crosses the differential spectrum of ACR (Buchvarova, 2006). The differential spectrum is given as the number of particles observed per (m2 s.ster.MeV/n) [part/ m2s.ster.MeV/n].The calculation of parameters α, β, x, y and is performed by Levenberg-Marquardt algorithm (Press et al., 1991), applied to the special case of a least squares. The described programme is realized in algorithmic language C++ .
Experimental data (Ei, Di) for protons and helium nuclei was got from Hillas (1972) for 20 solar cycle.
MODELING COSMIC RAY MODELING COSMIC RAY DIFFERENTIAL DIFFERENTIAL SPECTRSPECTRUMUM
MODELING COSMIC RAY DIFFERENTIAL SPECTRUM
10-3 10-2 10-1 100 101 10210-5
10-4
10-3
10-2
10-1
100
101
102
1
3
11975
solar maximum
IMAX92CAPRICE94AMS98
solar minimum
E, GeV
D(E)
, par
t/(m
2.s.st
r.MeV
)
1 2 3 4 5 6 7 8 9 10 11
Fig.3.The modelled spectrum D(E) of CR protons for eleven levels of solar activity and measurements: period near to solar maximum - ■ IMAX92 and periods near to solar minimum – ●CAPRICE94 and ▲ AMS98.
The results from the differential energy spectrum D(E) of primary protons and helium for solar minimum and maximum for the Earth. The modelled spectra are compared with the measurements for the period near to solar maximum - ■ IMAX92 (Menn et al., 2000) and periods near to solar minimum – ● CAPRICE94 (Boezio et al., 1999) and ▲ AMS98 (Alcaraz et al., 2000a, b). This data practically coincides with our results for these periods.
MODELING COSMIC RAY MODELING COSMIC RAY DIFFERENTIAL DIFFERENTIAL SPECTRSPECTRUMUM
10-3 10-2 10-1 100 101 10210-6
10-5
10-4
10-3
10-2
10-1
100
101
102
32
1
1
IMAX92CAPRICE94AMS98
solar minimum
solar maximum
E, GeV/nucl
1 2 3
D(E)
, part
/(m2.s
.str.M
eV/n
ucl)
Fig. 4. The modelled spectrum D(E) of CR helium nuclei for three levels of solar activity and measurements: period near
to solar maximum - ■ IMAX92 and periods near to solar minimum – ● CAPRICE94 and ▲ AMS98. Curve 1 relates to solar maximum, 2 – to comparatively average level of the solar activity and 3 – to solar minimum.
COMPARISON OF THE MODELING COSMIC RAY SPECTRUM WITH THE FORCE FIELD APPROXIMATION AND 2D SOLUTIONS OF COSMIC RAY TRANSPORT EQUATION
Under some simplifying assumption, CR transport equation can be reduced to Force Field approximation. The force field parameterization of cosmic ray nuclei at 1 AU is given as:
D (E, Φ) is differential intensity of cosmic ray nuclei and DLIS(E+Φ) - LIS of cosmic ray nuclei (was taken from Burger
et al., 2000 ). E is the kinetic energy (in MeV per nucleon) of cosmic nuclei with charge number Z and mass number A, Φ=(Ze/A)ϕ - modulation strength (in MeV), and Е0= 938 MeV is the proton’s rest mass energy.
We calculate differential intensity of galactic protons and alpha particles from Eq.(2) at given values of the
modulation potential ϕ. Received spectra are fitted to Eq. (1). On the base on the calculated coefficients from our model for the proton and alpha particles, D(E) is estimated for different values of modulation parameter ϕ. The results from the differential energy spectrum D(E) of primary protons for four values of modulation parameter: = 400, 550, 700 and 1200 MV are shown in Fig.4. These results are given for helium nuclei In Fig. 5. The number ratio α/p = 0.05 in LIS (Usoskin et al., 2005) is used for alpha particles .
COMPARISON WITH THE FORCE FIELD APPROXIMATION
COMPARISON WITH THE FORCE FIELD APPROXIMATION
10-1 100 101 102
10-1
100
101
102
103
IMAX92 CAPRICE94 AMS98
E, GeV
D (E),
partic
le /( m
2.sr.
s.GeV
) φ=1200 GV φ=700 GV φ= 550 GV φ= 400 GV
Fig. 5. The differential spectra D(E) of GCR protons for modulation levels = 400, 550, 700 and 1200 MV. Measurements with IMAX92 are in the vicinity of 700 MV, CAPRCE94 and AMS98 of 550 MV.
Fig. 6. The differential spectra D(E) of GCR helium nuclei for modulation levels = 400, 550, 700 and 1200 MV. Measurements with IMAX92 are in the vicinity of 700 MV and CAPRCE94 and AMS98 of 550 MV.
Figs. 4 and 5 show that our empirical model (1) well agrees with the data from force field approximation (2) for the protons and alpha particles. The measurements with IMAX 92 are in the vicinity of 700 MV, while CAPRICE 94 and AMS 98 of 550 MV. On the fit the standard deviations for protons and alpha particles are in the range of 1.5%.
10-1 100 101 10210-4
10-3
10-2
10-1
100
101
102
103
IMAX92 CAPRICE94 AMS98
E, GeV
φ=1200 GV φ=700 GV φ= 550 GV φ= 400 GV
D (E),
partic
le /( m
2.sr.
s.GeV
)
Comparison of modeling spectrum with two-dimensional model without drift of Potgieter and Moraal (Astroph.J., 1985)
1 0 - 3 1 0 - 2 1 0 - 1 1 0 0 1 0 1 1 0 21 0 - 5
1 0 - 4
1 0 - 3
1 0 - 2
1 0 - 1
1 0 0
1 0 1
1 0 2
2
P o t g i e t e r a n d M o r a a l , ( 1 9 8 5 ) : K I I = 2 . 2 8 x 1 0 2 2 c m 2 s - 1
K I I = 0 . 6 x 1 0 2 2 c m 2 s - 1
1
4
1 086
s o l a r m a x i m u m
s o l a r m i n i m u m
E , G e V
D (
E),
pa
rt/(
m2
.s.s
tr.M
eV
)
1 2 3 4 5 6 7 8 9 1 0
Fig.7. The modelled spectrum D(E) of CR protons for eleven levels of solar activity is compare with solutions (Potgieter, Moraal, Ap.J., 1985) of the cosmic ray transport equation without drift for nearly flat neutral sheet with parameters : (K┴)0 = 2x1020cm2s-1, rb = 50AU,
(KII)0 = 2.28x1022cm2s-1 (near solar minimum)(KII) 0 = 0.6 x1022 cm2 s-1) (near solar maximum )
ra, AU is planetary average distance of the planet from the Sun
PEUV is parameter of decreasing of solar UV radiation, proportional to 1/ ra2
PCR is parameter of CR increasing intensity because of solar modulation.
For 1 AU, the energy flux with wave length below 100 nm is about 3.3 erg/cm2.s Galactic CR energy flux is 210-2 erg/cm2.s
The solar EUV radiation is weaker than the galactic cosmic ray intensity for Jupiter, Saturn, Uranus and Neptune. This comparison shows the importance of the galactic
cosmic rays in the formation of outer planet ionospheres. We assume mean differential gradient of GCR in the interplanetary space as 3%/AU
Values of planetary average distances ra from the Sun, parameter PEUV of decreasing of solar EUV radiation, and parameter PCR of increasing of galactic CR intensity
COSMIC RAY INTEGRAL SPECTRA
Using the values of the parameters PGCR of increasing of galactic CR intensity for outer planets and the received values of coefficients from Eq.(1) for the Earth the integral spectra of the outer planets are obtained. At integration we use equation (3):
D(E) - differential spectrum of galactic and anomalous CR. D(>E) - integral spectrum, expressed by the number of particles per unit solid angle,
square centimetre, and second, with total energies at least E. The integration on E begins from the energy, corresponding to the geomagnetic cut-off rigidity in the point of measurements. (The low’energy part of the spectrum (below some tens of GeV) is dependent of the geographycal position).
In our model integral spectra are computed only in first approximation. We assume mean differential gradient of GCR as 3%/AU for all rigidity, irrespectively of the distance in the heliosphere or solar activity level. The measurements give radial gradient ~1.5%/AU at solar minimum and 3%/AU at solar maximum in inner heliosphere .
(3) )( dEEDEDE
1 0 - 1 1 0 0 1 0 1 1 0 2
1 0 - 2
1 0 - 1
1 0 0
B o b ic e t a l. , ( 2 0 0 6 ) : ( K I I = 2 . 0 x 1 0 2 2 c m 2 s -1 )
E a r th J u p ite r
E a r t h ( s o la r m a x . ) J u p i t e r ( s o la r m a x . ) E a r t h ( s o la r m in . ) J u p i t e r ( s o la r m in . ) S c h o p p e r V o y a g e r 2 ( 1 9 7 9 )
D(>
E),
pa
rtic
le/(
cm
2.s
ec
.ste
r)
R , G V
Fig. 8. The modeled integral spectra D(>E) of CR protons for maximum and minimum levels of solar activity for Earth and Jupiter. The computations are compared with data: ᵡ - Schopper (1967) for Earth, ▼- Voyager 2 (5 July 1979 – near to solar maximum) for Jupiter and theoretical results of 2D stochastic model built from Bobik et al., (ECRS, 2007).
It is seen from Fig. 8 that in Bobik model, the integral spectrum of the Earth almost tally with Jupiter’s integral spectrum at solar minimum. It is due to lower average values of integral radial gradient in this model. Actually at high K radial gradient has lower values in the inner heliosphere. (The radial intensity gradient has lower values in higher energy also). In Bobik model for positive solar period (A > 0) are used the following values (KII)0 = 2x1022 cm2s-1 , tilt angle ϑ = 300 and the ratio between the parallel and the perpendicular diffusion coefficient is (K┴)0 = 0.025. Burgers’s model (Burger et al., J. Geophys. Res., 2000) is used as local spectrum of protons. The model with these parameters reproduces the spectrum measured by AMS – 01 (corresponding to a modulation strength 510 - 550 MV (in agreement with independent estimation)).
COSMIC RAY INTEGRAL SPECTRA
COSMIC RAY INTEGRAL SPECTRA
1 0 - 1 1 0 0 1 0 1
1 0 - 2
1 0 - 1
1 0 0
B o b ic e t a l. , ( 2 0 0 6 ) : ( K I I = 2 . 0 x 1 0 2 2 c m 2 s -1 )
S a tu r n U r a n u sN e p tu n e
D(>
E),
pa
rtic
le/(
cm
2.s
ec
.ste
r) S a t u r n ( s o la r m a x . ) U r a n u s ( s o la r m a x . ) N e p t u n e ( s o la r m a x . ) S a t u r n ( s o la r m in . ) U r a n u s ( s o la r m in . ) N e p t u n e ( s o la r m in . )
R , G V
Fig. 9. The modeled integral spectra D(>E) of CR protons for maximum and minimum levels of solar activity for Saturn, Uranus and Neptune. The results are compared with computations of Bobik et al., (ECRS, 2006) for 2D transport model with drift.
The observed CR spectrum in the energy interval from 30 MeV/n to 3.102 GeV/n has important contribution for the physical processes in the ionospheric CR-layers of the planetary ionospheres. The GCR modulation is essential
in this energy interval. The energy interval E = 1 – 30 MeV/n is related to the polar cap region in Earth’s ionosphere. The particles
from this interval can penetrate directly in the cusp regions, where they cause enhanced ionization, heating and excitation of the upper and middle atmosphere.
Application of the model for calculation of the electron production rate q(h) profiles in energy intervals of galactic CR and anomalous CR component
The differential D(E) spectra (1) of galactic and anomalous CR are used for computation of the electron production rate profiles in the atmospheres and ionospheres, at which the ACR component is also taken (Mateev L., Bul. Geophys.J, 1997; Velinov et al., JASR, 2008). In Fig. 10 electron production rate q(h) profiles are presented. The computational data present the ionization by galactic CR and their anomalous CR component, which has significant influence above 80 km. The profile with maximal values refers to the polar cap region, the one with minimal values presents the electron production rate at the equator. The positive deviation in the upper part of the polar cap profile shows clearly the contribution of the anomalous CR component
Application of the model for calculation of the electron production rate q(h) profiles due to particles of energy intervals of galactic CR and anomalous CR component.
0.8 1.0 1.2 1.4Relative Electron Production Rate q_ E / q_ S
0
100
200
300
400
5000 4 0 9 0
Rela
tive
Altit
ude
abov
e 1
bar l
evel
, km 15
16
17
18
19 Log_
10 N
eutr
al N
umbe
r De
nsity
, cm
-3
0.8 1.0 1.2 1.4Relative E lectron Production Rate q_ E / q_ S
0
100
200
300
400
5000 4 0 9 0
Relat
ive A
ltitu
de a
bove
1 b
ar le
vel ,
km
15
16
17
18
19 Log_
10 N
eutra
l Num
ber D
ensit
y, c
m-3
The relative electron production rate profiles (for elliptical per spherical geometry) in Kronian ionosphere for latitudes: 0o - equator, 40o - middle
latitudes, and 90o - polar regions
The lowest value for electron production rate profiles are obtained in case of spherical geometry. The observed difference is the smallest for the equator and the biggest for the polar regions.
THANK YOU FOR YOUR ATTENSIONTHANK YOU FOR YOUR ATTENSION