~
Adv. Space Res. Vol. 30, No. 4, pp. 1033-1044, 2002 Published by Elsevier Science Ltd on behalf of COSPAR Printed in Great Britain 0273-1177/02 $22.00 + 0.00
Pergamon
www.elsevier.com/locate/asr
PII: S0273-1177(02)00497-0
A REVIEW OF SOLAR PROTON EVENTS DURING THE 22ND SOLAR CYCLE D. F. Smart and M. A. Shea
Air Force Research Laboratory (VSBS), 29 Randolph Road, Hanscom AFB, MA, 01731, USA. ABSTRACT Solar cycle 22 had significant, large fluence, energetic particle events on a scale reminiscent of the 19th solar cycle. Examination of the characteristics of these large events suggests that some of the old concepts of spectral form, intensity-time envelope and energy extrapolations, used to estimate the dose from large events that occurred during previous solar cycles should be re-evaluated. There has also been a dramatic change in perspective regarding the source of solar protons observed in interplanetary space. Very large fluence events are associated with powerful fast interplanetary shocks. The elemental composition and charge state of these events is suggestive of a dominate source in the solar corona and not from a very hot plasma. Furthermore, there is a strong suggestion that the intensity-time profile observed in space is dominated by the connection of the observer to an interplanetary shock source rather than to a unique location near the surface of the sun. These concepts will be examined from the perspective of energetic panicles contributing to the dose experienced by an astronaut on an interplanetary space mission. Published by Elsevier Science Ltd on behalf of COSPAR.
INTRODUCTION
This paper is oriented toward the viewpoint of radiation effects to human blood-forming organs (vital organ dose). We will discuss the energetic particle radiation capable of penetrating the shell of a space vehicle and the human body to reach the blood forming organs (at least 5 cm of tissue). From this viewpoint, the energetic charged panicle radiation must have a minimum energy of 30 MeV per nucleon. From a solar astrophysics point of view, the increased energetic panicle sensor sensitivity to very low flux levels that has occurred over the past decades provides for hitherto unavailable information on the particle composition and possible origin. However, from a radiation dose point of view, small events of limited duration with maximum particle fluxes of a few per cm 2 per hour are not very significant. Therefore we will not put much emphasis on the very small events routinely measured by the improved instrumentation. OVERVIEW OF SOLAR PROTON EVENTS DURING SOLAR CYCLE 22 Figure 1 illustrates the >30 MeV solar proton events during solar cycle 22. The horizontal lines on this figure indicate one sigma in a log-normal distribution of the solar proton event fluence distribution over the last four solar cycles. This one sigma characterization of the fluence distribution derived by Nymmik (1993) from the data of Feynman et al., (1993) seems to be one reasonable way of categorizing solar proton events. Stassinopoulos et al. (1996) adopted an order of magnitude classification system categorizing solar proton events as a matter of convenience. The lower limit of this figure, an omnidirectional fluence of 105 per cm 2 is comparable to the daily exposure from the integral flux of galactic cosmic ray protons. The solar cycle 22 solar proton events and fluences at energies of>50, >60 and >100 MeV are shown in Figures 2 - 4. In Figure 5 we display some of the very large events of the last solar cycle from the perspective of the size of a cell in the human body. We have used 100 square microns as the unit of area. The scale on the left of this figure indicates the number of proton "hits" per hour during these large events. Also indicated by horizontal lines are the flux levels that correspond to "hits" per day and "hits" per week.
1033
1034
D.F. Smart and M. A. Shea > 80 MeV EVENT I N T E ~ t A T E O PROTON O~(IDIREC'rtONAL R . U B K E > IB0 MeV EVENT INTEGRATED PROTON ~
A
E 0
¢n to
o
101o
•. ' . 1 . ' . ' . ' . = ' ; 2 ! ; : ; ~ . : ' ' . , . ' . ' ' . = . ' . ' . ' . = . ' - . ! ; ' ; = . ' . ' . "
E U
10o
..J
u.
..,...,...,...,...,...,...,..+,...,...
10o
W e-
ii
10a
n
LU O Z u.I ::)
1010
,'+"
]
FLUENCE
107 lo0 lOS
87 88 8g go gl
86
o
LARm=
lo0
0
¢: 107 Z uJ ::) i,. J
g2 93
94
95
YEAR MeV solar proton events during
It~
lo0 10s I . . . . . .
96
86
87
I
It:'. 88 89
i|• . . m• .+
g0
91
92
g3
94
95
96
YEAR
Fig. 1. The > 3 0 solar cycle 22, The horizontal lines indicate one sigma of the log-normal distribution.
Fig. 2. The > 50 MeV solar proton events during solar cycle 2 2 .
• SOMeVEVENTINTEGRATEDPROTONOMNIDBECT1ONAL FLUBCCE ) 100 MIV LrVENT ~ 1 1 D l O P O . .
o
!1.
.
.
,
. .
.
,
.
.
.
,.
.
.
,
,
,
,
,.
,
,
, .
,
,
,
.
.
.
i
.
.
.
1+L
M tO
107
n w
3 u .
87 88
89
90 91
0 z w
lil I i.... g2
93
g4
95
.J U.
96
YEAR
1
0 .
.
.
0 .
.
.
107
I
1. 1~
I
,lll
"" ....... 86 87 88 89
• 90
.L
.I
91
92
I
|=..J 9 -'1
YEAR
.|...
94
95
96
Fig. 4. The > 100 MeV solar proton events during solar cycle 22.
Fig. 3. The > 6 0 MeV solar proton events during solar cycle 2 2 .
o
~
lol
lOa
86
PIIOlrON ~
100 . . . , . . . , . . . , . . . , . . . , . . . , . . . , . . . , . . , . . .
, .
~
> 30 MeV
.
> 50 MeV
1.0 Proton per Day
> 60 MeV > 100 MeV 1 0-4 19
20
21
22
23
24
25
OCTOBER1989
26
27
28
29
30
31
(UT DAY)
Fig. 5. The very large events of October 1989 displayed from the perspective of the size of a cell ( - 1 0 0 square microns).
Reviewof SolarProtonEvents
1035
THE PARTICLE SOURCE Both the classical solar flare source scenario and the shock acceleration scenario have enthusiastic advocates. Many of the "large" solar particle events are associated with "large" solar flares. However, many of the proton events observed at the Earth cannot be unambiguously and confidently time associated with specific solar flares. Furthermore, measurement of the elemental and isotopic composition of solar particles during large flux and fluence events was found to be consistent with the particles having passed through less than 30 mg cm -2 of matter from the acceleration site to their detection location (Mason, 1987). Thus, the solar photosphere is eliminated as the source region of the particles observed in space since there is no evidence of fragmentation due to interaction with the significant mass of the solar atmosphere. The intensity-time profiles of solar energetic particles leave no doubt that interplanetary shocks accelerate ions. The fast coronal mass ejection (CME) generated interplanetary shock is now considered to be the principal source of the ions observed in space (see Reames, 1995, for a review). There is also the concept of a composite solar particle event to which both sources contribute flux (Cliver, 1996). A relatively recent grouping of solar particle events is to classify them according to the associated solar flare X-ray activity (see Table 1). This results in a classification of either "impulsive" (those associated with impulsive soft Xray events) or "gradual" (those associated with long duration soft X-ray events). The "impulsive" solar particle events are the group that appears to be associated with the solar flare process. The long duration (i.e. "gradual") Xray event is believed to be a proxy indicative of a coronal mass ejection, and hence this class of events is assumed to be the result of the shock acceleration process. An extensive study of the "impulsive-flare" associated energetic particle events by Cane et al. (1986) and Reames et al. (1994) shows that these events are generally small flux and fluence events measurable in space over a restricted heliolongitudinal range, about one radian of heliolongitudinal distance centered on the most favorable propagation path from the solar flare location. The charge state of solar energetic ions that are associated with the impulsive soft X-ray flare is consistent with a multi-million degree hot (-2 x 107 OK) plasma source.
Table 1. Properties of Solar Proton Events Associated with Impulsive and Gradual soft X-ray events. (Adapted from Reames, 1995)
Impulsive X-ray event and associated particle event characteristics Particles Electron rich 3He/4He Fe/O H/He Charge state of Iron DurationHours Heliolongitude range
Proton rich ~1 ~1 10 ~20 Days <30° (generally "well connected")
Associated Solar Radio type Associated Solar X-ray emission CME association Associated Interplanetary shock Events per year (At solar maximum)
m, v, (ii) Impulsive
~
~1000
Gradual X-ray event and associated particle event characteristics - 0.0005 ~0.1 ~ 100 ~14 ~ 180° I/, IV Gradual 96 % Yes -I0
1036
D.F. Smart and M, A. Shea
Large solar particle events seem to be associated with the occurrence of a fast coronal mass ejection or, in the absence of actual CM]~ observations, C M E proxies. This association is complicated by the fact that large fast coronal mass ejoctions also seem to have an association with "big" solar flares. The long duration solar X-ray event is considered to be an excellent C M E proxy. Figure 6 illustratesthe fastC M E associated with the 24 October 1989 solar cosmic ray event observed at the Earth. Large solar particle events tend to have a "normal" composition that is relatively consistent with abundances found in the solar corona or in the solar wind, although there are event to event variations, particularly in the helium/hy&ogen ratios. The class of solar particle events associated with a CME/shock source seems to be detectable at very large heliocentric angles from the centroid of the CME/sourcc heliolongitudc position. The largest flux values are observed at locations in space well connected via the interplanetary magnetic fieldtopology to the presumed centroid of the CME/shock source heliolongitude. GOES CORRECTED INTEGRAL FLUX ";"
1 0
$
.....
24 5 Fig. 6. The fast CME associated with the 24 October 1989 solar cosmic ray event.
, .....
, .....
25 5
, .....
, .....
26 5
, .....
, .....
27 5
, .....
, .....
28 5
i .....
, ....
29 5
OCTOBER 1989 (UT Day) Fig. 7. The "classic" solar proton event as illustrated by the 24 October 1989 event.
The charge state of energetic ions in "large" particle events appears to be consistent with coronal temperatures (-106 °K plasma). The initialcharge statemeasurements were made at low energies, of the order ofa M c V (Luhn et al.,1985). The charge stateshave now been measured at energies of 10s o f M c V s (Octliker ctal., 1995). These higher energy charge states are generally consistent with the low energy charge states with an indication that the charge state may be a function of energy. Tylka ct al. (1995) found an ionic charge state of-14 for iron nuclei at -200 M e V per nucleon. The intensity-time profile observed at a specific location in space depends on the "connection" of the observer to the particle source (Cane, 1986). W c will use the terms "near-sun-injection" and "extended interplanetary shock source" to distinguish between types of solar proton events observed in interplanetary space. The term "ncar-suniniection" will refer to the class of solar proton events in which the intensity-time profile at the observer's position is dominated by the particle flux injected near the sun onto the interplanetary field lines passing by the observer. This scenario represents the "classic" solar proton event illustratedby the 24 October 1989 proton event shown in Figure 7. This class of solar particle event often has a relatively "hard" (hard meaning more than an average number of particles at high energies) spectra. The term "extended intcmlanetarv shock source" will refer to the class of solar particle events that does not conform to the "classic" solar proton intensity-time profile, but instead has a characteristic that the flux continues to increase until a maximum flux is observed as a powerful and fast interplanetary shock overtakes the observer, aRcr which the flux decreases. An example of this type of event is the March 1991 solar proton event as observed at the Earth and illustratedin Figure 8. This example is probably typical of this type of event which usually has a relatively "soft" (softmeaning fewer than an average number of particles at high energies) spectra.
Reviewof SolarProtonEvents
1037
FLUX ENHANCEMENTS ASSOCIATED WITH INTERPLANETARY SHOCKS Particle flux enhancements associated with the passage of powerful and fast interplanetary shocks is a concept still being defined. A common misconception is that there will always be a flux enhancement with the passage of an interplanetary shock. This is true only at low energies (~1 MeV and below). At higher energies (>30 MeV) it is the exception when a flux enhancement is observed with the passage of an interplanetary shock. Kallenrode (1993) and Kallenrode et al. (1993) examined the entire HELIOS particle data base and found at 20 MeV, the correlation between the shock jump plasma parameters and the flux increase was approximately zero. Because there are no definitive analytical rules for predicting the flux increase associated with the passage of an interplanetary shock we will use an illustrative example. This example is the last significant particle event of solar cycle 22, which began on 20 February 1994, and shown in Figure 9. This is an event having both an initial near-sun injection and a flux profile dominated by an extended interplanetary shock source. The time associated M4/3B solar flare occurred at 0141 UT on 20 February at a heliographic position (as viewed from the Earth) of North 07 °, West 02 °, near the solar central meridian. This solar event had all of the CME proxies, and it is presumed that a fast CME was launched toward the Earth even though there was no equipment in space to observe the actual CME. The initial particle flux was detected at the Earth shortly after the solar flare electromagnetic emissions were observed. The particle onset times at the Earth were in a time sequence consistent with a near-sun energetic particle injection, the components of which propagated along the interplanetary magnetic field lines from the sun to the Earth at speeds consistent with the particle energy. Inspection of Figure 9 shows that only the highest energy particles recorded at the Earth have the intensity-time profile that is dominated by a nearsun injection. At energies below 30 MeV, the observed intensity-time profile indicates that an extended interplanetary shock source dominated the flux profile observed at the Earth. The maximum low energy particle flux was observed as the fast interplanetary shock passed over the spacecraft and the Earth, initiating a SSC (storm sudden commencement) magnetic storm which began at 0901 UT on 22 February. The hypothesis is that a fast CME propagating toward the Earth was generating a shock wave which was accelerating particles. We can hypothesize that when this shock wave was within 1/2 AU of the Earth, the low energy particle flux observed at the Earth was dominated by the particles accelerated at the shock. At higher energies (>20 MeV), the increased flux due to the extended interplanetary shock source is not apparent until the shock is within 0.2 AU of the observer. COMPOSITE S O L A R P R O T O N EVENT ....
i .....
i .....
i .....
i .....
GOES
i ....
7
data (February
1994)
1 0s 105 ~ . . . . , ! i
1 04
•~
>
E
O
,o L
.o
Source
.'~,
1 01
P=.50
~
,
-J
LL ,.J ; t . . t ~
1%30 2&5
e
m
^ / ] \;. ~, /~/" "~',
~
~v
P4 =uu,v ~ =o v,~.
(D
£
X J 14.
~=~,.source
1 02
O4
, ..... , ....
Inteq~netaryShock
P~=~
:S
u)
v
..... , ..... , .....
Near Sun
104
,oor_
'
/
,o2
P>100
245 250 255
MARCH 1991 (UT DAY) Fig. 8. The e x t e n d e d i n t e r p l a n e t a r y shock t y p e o f e v e n t as illustrated b y t h e M a r c h 1991 solar p r o t o n e v e n t as o b s e r v e d at t h e Earth.
10-4L.--.J
51.0
.....
51.5
, .....
52.0
I .....
52.5
J .....
53.0
J .....
53.5
"
54.0
Day of Year (1994) Fig. 9. The F e b r u a r y 1 9 9 4 e v e n t as o b s e r v e d at the Earth, This " c o m p o s i t e " solar particle e v e n t had b o t h a " n e a r sun i n j e c t i o n " and a f l u x e n h a n c e m e n t a s s o c i a t e d w i t h an i n t e r p l a n e t a r y shock.
1038
D.F. Smart and M. A. Shea
SOLAR ENERGETIC PARTICLE FLUX ANISOTROPY The particleflux anisotropy is generally defined as the ratioof the maximum particleflux divided by the average of the particle flux in all directions. In the isotropic case, there is equal particle flux from all directions. In the anisotropic case, there is an excess of particle flux in one specific direction. If it is assumed that the energetic charged particlesare constrained to travel along an interplanetarymagnetic field line with littlescatteringacross fieldlines,then the particleflux anisotropy also reflectsthe rate at which particlesare injected onto the observer's interplanetaryfieldline with respect to the rate at which the particleflux is transported away from the observer. In the case where there is a "near sun" injectionof energetic charged particles,then the initialparticleflux observed "downstream" will be initiallyanisotropic,with the degree of anisotropy decreasing as particlesthat have passed the observer scatterand are reflectedback by subsequent scatteringin the interplanetarymagnetic field. Particle flux anisotropy is a transient phenomena, and at the energies important to dose calculations, at radial distances from the sun between the Earth and Mars, a duration of a few hours is typical. In the example shown in Figm~ 10, there was no significant anisotropy remaining after 0.1 day. This figure illustrates the proton flux anisotropy at -30 MeV observed by the HELIOS spacecraR at 0.95 AU during a large anisotropic solar particle event in June 1980 (day of year 173). The anisotropy amplitude (denoted by "A" in Figure 10) ranges from extreme (~1) to low anisotropy (-0.2). An artifact of fitting a cosine function to coarsely sectored data results in large anisotropies having values >1 as shown in the first panel of Figure 10.
(1980)
173.077
Fig. 10. observed from the The line direction, indicates
173.126
A=0.2
173.175
The proton flux anisotropy at - 3 0 MeV by the HELLOS spacecraft. The distance center indicates the relative particle flux. labeled "B" indicates the magnetic field and the line to the right of the center the anti-sun direction.
COMPOSITION OF SOLAR ENERGETIC PARTICLE EVENTS The "impulsive-flare-associated" solar energetic particle events seem most closely related to the solar flare particle acceleration process. Initially, these events were considered as "composition anomalies" since there were distinct differences from "large" solar particle events (see Table 1). They were sometimes called "3He rich" or "iron rich" events. In general, the fluxes in the impulsive events are orders of magnitude smaller than "significant" solar particle events such as those illustrated in Figures 1-4. The analysis of Reames et al. (1994) concluded that in the -I - 20 M e V energy range, large solar particleevents tended to have a "normal" composition that was relativelyconsistent(perhaps within a factor of 3), although there are event to event variations, particularly in the helium/hydrogen ratios. The ratio of the abundances of the elemental composition observed in large solar energetic particle events to the abundances of the elemental composition observed in the photosphere appears to be a function of the first ionization potential of the individual elements. Some researchers contend that large solar energetic particle events are a reasonable sample of the solar corona (Breneman and Stone, 1985, Mewaldt and Stone, 1989). The available data reflect the composition at relatively low energies (a few MeV to perhaps 10 MeV). The analysis of Mazur et al. (1992) suggested that there was relative consistency of the elemental composition during an event (up to energies of ~30 MeV). However, the latest analysis of the heavy ion composition at high energies (>100 MeVANuc) for very large flux events (Tylka et al., 1997) suggests that the iron to oxygen ratio is a function of energy with variations larger than those reported by Mazur et al. (1992). SOLAR ENERGETIC PARTICLES - CHARACTERISTICS OBSERVED AT 1 A U The majority of the solarparticledata that existshas been collectat the Earth by ground-based sensors or by Earthorbiting spacecraft. Prior to 1967, most of the solarproton data were derived from analysis of the response of the Earth's ionosphere to solar particleevents; the spectral information derived lacks the accuracy of the spacecraft
Review of Solar Proton Events
1039
data. Most of the spacecraft data collected have been at energies <100 MeV. There is relatively good data coverage for energies in the range of 10, 30, and 60 MeV extending from 1967 to the present. These data are sufficient to derive meaningful statistics for modeling purposes. There is a serious deficiency in the flux and spectral characteristics of solar particle events at energies >100 MeV. There are a number of formulations to model the spectra of small or medium size flux or fluence events, but, in general, these forms do not adequately describe the flux profile encountered in large and very large solar particle events. Particle Flux Evolution and Radiation Dose The possibility that a large very energetic solar cosmic ray event could generate a significant radiation dose to astronauts in a short time interval is a matter of concern. Some of the original radiation dose estimates from the large energetic proton events of the 1950's and 1960's were many rads per hour (Foelsche, 1965). However, more recent dose calculations for these same events (Wilson, et al., 1991) are considerably lower. One of the reasons for these differences is the use of new standards for converting radiation types to dose (ICRU Report 71, 1998). Another reason is a better definition of the solar particle spectra for the more recent events and a subsequent re-evaluation of the spectral characteristics of some of the earlier events (Wilson et al., 1991). Inherent in the early dose estimates was the assumption of "explosive" particle release from the sun and a subsequent "explosive" increase in dose from a very fast increase in high energy particle flux. The spacecraft observations of the large energetic solar proton events of solar cycle 22 do not support the "explosive" particle release and flux increase scenario. Typically, in the 30 to several hundred MeV energy range, the time scale of hours from particle onset to flux maximum was observed. From the analysis of the large events that occurred in the 22nd solar cycle we have learned it is a serious mistake to use one simple value to characterize the spectral characteristics of a large high energy solar cosmic ray event. This is perhaps one of the major reasons that such high dose values were calculated for the events in the 19th solar cycle. Lacking accurate measurements at a number of energies, the early dose calculations were made using one spectral slope determined from the known energies (usually either very high in the GeV range or in the low MeV range), and extrapolating this spectral slope to all other energies. This method typically results in a considerable overestimate of the flux at the other end of the energy spectra thus leading to higher dose calculations. When the 29 September 1989 event occurred, the NOAA Space Environment Center radiation dose prediction was significantly larger than the actual measurement on the Concorde or by subsequent calculation by a number of researchers (O'Brien et al., 1996, 1998; Wilson, et al., 1991; Wilson and Nealy, 1992). The most likely reasons for these differences were the utilization of an outdated radiation transport model and the use of an unrealistic exponential model of the solar particle spectrum to extrapolate from low to high energies. The Flux and Spectral Evolution of the 29 September 1989 event. The 29 September 1989 solar cosmic ray ground-level event was large and significant by every measurement; Smart and Shea (1991) evaluated this as the third largest high energy solar proton event since the first ground-level solar cosmic ray event was identified in 1942. Particles with energies >20 GeV were present (Swinson and Shea, 1990); the event was observed at every operating cosmic ray detector on the Earth. The 29 September 1989 event occurred in the "modem" spacecraft era, and there are direct spacecraft measurements of the particle flux up to energies of about 800 MeV. The particle spectrum of this event is complex and not described by simple parameters. The initial reports of particle energy spectra for this event appeared to be inconsistent. Mathews and Venkatesan (1990) determined the spectral slope at the mean response of the high energy detectors (5 to 7 GeV) whereas Humble et al. (1991) and Smart et al. (1991) determined the spectral slope at 1 GV rigidity. Subsequent analyses of this event (e.g., Cramp et al., 1993) show that the spectral slope is a function of energy and time. The slope is extremely hard at low energies; the magnitude of the slope increases (i.e. becomes "softer") with increasing energy and time. An intensity-time plot of the GOES spacecraft data obtained for the 29 September 1989 event is presented in Figures 11 and 12. (The GOES proton detectors has a serious side penetration problem whereby the high energy protons register in the low energy proton channels. Error flags and corrections for this problem removed the early part of the event for protons below 200 MeV.) The time-evolving spectra observed by the GOES spacecraft has been published by Smart and Shea (1997).
1040
D.F. Smartand M. A. Shea
2 ..... , ..... i ..... , ..... , . " , : . - " , - . - , . " ~
,,,.
102
.~,,,,v 2
lOl
41:z MoV
.~p-..,"-- ~ - ~
s0.s ~ v
~
,-7 ~ ~ Z
142 MeV 0
4 S
X
-/
10"I : ~'
2 10-=
S92 MeV ~ 620 MoV
t~
10-2
X
s
~ -1
IJ.
2
•.':.:
-",,.,,,%.
2
"~
•
1 0-s
S
1 o --~
!
• t4
4 $ 2
2
10-4
440 MoV" - / ~ . , ~
670 MeV
U_
" ~
"l .j
$
10-1
:..,.~,:"~" ....."~'--%.r""
2
~
1 oo
l,,,,,U,,,,.l.,,,,
10o !
(~ E 0
l.,.,,l,,,.,lUU,.,l,,,,,
....
-29.00
I .....
I .....
29.50
I ....
,I
.....
30.00
SEPTEMBER
1989
I .....
I .....
1
I .....
"4
029.00
30.50
.....
Alpha particle
I .....
29.50
SEPTEMBER
(UT day)
Fig. 1 I . The 29 September 1989 solar proton flux from 11 to 800 MeV as observed by the GOES 6 and 7 Spacecraft.
. ,,=,1
I .....
I .....
30,00 1989
I .....
flux
I .....
I .....
30.50 (UT
Day)
Fig. 12. The 29 September 1989 solar alpha particle flux from 41 to 440 MeV/Nucleon as observed by the GOES 6 and 7 Spacecraft.
Fluence Distributions
Analysis of solar proton events measured at the Earth has led Feynman et al. (1993, 1997) to conclude that solar proton fluence during the "active years" of the solar cycle is well approximated by a log-normal function. Other analysis of the same data using exlreme value statistics (Xapsos et al., 1996) suggests that there may be two distributions; however, the application of either method gives similar results for the probability of a very large flux value being encountered during a long duration interplanetary mission. These models can be used to estimate the probable free space exposure to solar protons for various space missions beyond the Earth's geomagnetic shield as shown in Figure 13. Similar models for the heavy ion component of the large solar particle events observed during solar cycle 22, have been published by Tylka et al., (1997). Application of these models can be used to estimate possible solar particle exposure to interplanetary spacecraft from which risk assessments can be derived.
w
> w w o z w o U. z 0
2
[-. 0.1
4
1 N'M Z I " ~
"~
1',2"'%
YR~_ I .4 ;.,,J/ ] J)'% XI ~J "N 2YH--YI I/'P~I,Z" ",Z T',, r'..r' " ~ " ~ ~, 13YR,F4 I ~
i '1
~ 0.01 ? 0 a. 0.001 10o9
f
i
1010 FLUENCE ( ern -2 sr-1)
Fig. 13. The probable exposure to > 3 0 MeV protons during interplanetary missions. From Feynman et al. (1993). Reproduced by permissions of the American Geophysical Union.
1011
Review of Solar Proton Events
Estimations of the fluence distribution of solar proton events over a very long period of time involves measuring the induced radioactivity in moon rocks and meteorites. Analysis of these data, combined with the Earth measurements over the past four solar cycles indicates that the very long term solar proton event fluence distribution may be represented by a broken power law (Reedy, 1996) as illustrated in Figure 14. The fluence of solar proton events observed at the Earth from 1954 through 1991 is in the left part of the figure. The right side of the figure shows the solar proton fiuences estimated from the analysis of induced radioactivity in moon rocks which encompasses a time span of perhaps one million years. The increasing slope of the fiuence distribution for very large events suggests that there may be limits to the acceleration process, and extraordinarily large events arc very rare.
A
Flux Distributions Using the dam from the IMP 4 and IMP 5 spacecraft, van Hollebeke et al. (1975) determined that the proton event peak flux distribution could be fitted to a power law having a slope of -1.15. Cliver et al. (1991), analyzing the solar particle events detected by the IMP 8 spacecraft from November 1977 to October 1983, also fitted the event frequency distribution with a similar power law function. A more recent analysis of solar proton event data (Smart and Shea, 1997) concluded that there has been an unconscious exclusion of the very large events in the previous work due to instrumental limitations or restricted data samples. This more recent analysis of all of the events observed during solar cycles 20, 21 and 22 found two distinct particle flux populations• There was a distribution of the most often occurring events, which for >10 MeV proton events can be described by a functional form having a slope o f 0.47. This population extends from very low fluxes to flux values of approximately 1000 (cm2-s-sr) -l. Then there is a population of rare, very large >10 MeV proton events that can be described by a functional form having a slope of-1.4 as shown in Figure 15. P R O T O N EVENT PEAK FLUX DISTRIBUTION Data from Cycles 20, 21, 22
Cumulative Probabilities. Large Solar Particle Events 10' 10' 10" 10~ 10" 10r'
LL
I@
0 r-
g: .~ 10"
~
,..~ ,,..~ ,,,.,,,~ ,,.,,.I ,..,.I ,,,...,I ,,.
i 1954"1991 ~
' F'°"
1041
101
w ,.J
°
1 o=
~'l
"'''''"
'
''.',!
!
''''",
'
'''"!
' ' ' ' "
'
'''"1
OTONS >10 MeV
10" ~
2
O
UJ "o. L,_
10
10.3
N el)
r~O IZ LIJ > U.I
10 i •
4 3 2
~.
1 O0 =UJ
i
Z
2
m
ftllllll~
Z
10'
IHIII~
lllitfll~
Jrfl~lld
llllll~
I I I I I I I~
Jltltlll~
illillB~
IIt/llld
[flJl~
ld 10" 10~3 10~' 10" F, Event Fluence > 10 MeV (protons/cm2)
Fig. 14. The fluence distribution of Earth-sensed proton events (heavy line) and proton fluence estimates from the analysis of induced radioactivity in moon rocks• (Adapted from Reedy, 1996)
10.7
,A
4
.I
.......
01
,
, ,
2 a45'
...I
........
102
,
2 a=5
,,.,,I
..............
t0=
= a45
I
104
PEAK >1 0 MeV FLUX (Crn2 S Sr)-I Fig. 15. The expected number of proton events during a solar cycle with the > 10 MeV proton flux exceeding a specified flux threshold•
1042
D. E Smart and M. A. Shea
SUMMARY
From the viewpoint of solar particle events observed at the Earth, solar cycle 22 was quite active, having a total proton fluence slightly below the 19th solar cycle. Table 2 gives a summary of the solar proton events observed during the last four solar cycles. Table 3 lists the fluence of the major proton events during the last solar cycle. Table 2. Summary of Solar Proton Fluence for SolarCycles 19 - 22
Cycle No.
Start
No. of Months in Cycle
End
No. of Discrete Proton*
Events
No. of Discrete Proton* Producing Regions
Solar Cycle Integrated Omnidirectional Solar Proton Flucnce >10 M e V > 3 0 M e V (era -2)
19 20 21 22
May Nov Jul O~
54 64 76 86
O~ Jun Sep Sep
64 76 86 96
126 140 123 120
65 72 81 74
47 56 57 45
7.2x10 lo 2,2x10 l0 1,8xlO 10 5.8xi010
(cm -2) 1.Sx10 E0 6.9x109 2.8x109 1.0xl01°
*Proton events with >10 MeV peak flux >I0 (cm2-s-sr)"1
Table 3. Summary of the Largest Solar Proton Fluence Occurrences for Solar Cycle 22
Event Integrated Omnidirectional Solar Proton Fluence
Cycle Begin No. Yr M o Dy
22
89 89 89 89 89 89 90 91 91 91 92 92 94
03 08 09 10 11 12 05 03 06 07 05 10 02
09 12 29 19 27 30 21 22 04 07 09 30 20
End M o Dy
-
03 08 10 10 12 01 05 03 06 07 05 11 02
24 18 02 30 03 02 (90) 31 26 21 12 11 05 22
Totals:
>10 MeV (cm 2) 0.12 x 0.76 x 0.38 x 1.9 x 0.51 x 0.21 x 0.035x 0.96 x 0.32 x 0.11 x 0.066x 0.35 x 0.099x
1010 1010 1010 10 I° 10l° 10l0 101° 1010 1010 1010 101° 10l° 101°
5.82 xl01°
>30 MeV (cm 2) 0.03x 109 1.4 x 109 1.4 x 109 4.2 x 109 0.16x 109 0.13x 109 0.14x 109 1.8 x 109 0.79x 109 0.01x 109 O.02x 109 0.43x 109 O.02x 109 1.Oxl010
Reviewof SolarProtonEvents
1043
The measurements made by the new instrumentation developed during the 22nd solar cycle has resulted in a dramarie change in perspeerive regarding the source of solar proton events observed in interplanetary space. Solar proton events associated with impulsive solar flares seem to be small flux and fluenee events that are relatively unimportant for the dosimelric point of view. The very large fluence events are associated with powerful interplanetary shocks that are related to solar activity generating fast coronal mass ejections. The elemental composition and charge state of the large fluence events is suggestive of a dominate source in the solar corona and not from a very hot plasma. Furthermore, there is a strong suggestion that the intensity-time profile observed in space is dominated by the connection of the observer to an interplanetary shock source rather than to a unique location near the surface of the sun.
REFERENCES Breneman, H.H., and E.C. Stone, Solar Coronal and Photospheric Abundances from Solar Energetic Particle Measurements, Astrophysical Journal, 299, 959-963 (1985). Cane, H.V., R.E. McGuire, and T.T. +ConRosenvinge, Two Classes Of Solar Energetic Particle Events Associated With Impulsive and Long Duration Soft X-ray Flares, Astrophysical Journal, 301,448-459 (1986). Cliver, E., Solar Flare Gamma-ray Emission and Energetic Particles in Space, In High Energy Solar Physics, AIP Conference Proceedings 374, R. Ramaty, N. Mandzhavidze, and C-M. Hau (Eds), American Institute of Physics, New York, 45-60 (1996). Cramp, J.L, M.L. Duldig, and J.E. Humble, The GLE of 29 September 1989, Proceedings of the 23rd International Cosmic Ray Conference, 3, 47-50 (1993). Feynman, J., G. Spitale, J. Wang, and S. Gabriel, Interplanetary Proton Fluence Model, J. Geophys. Res., 98, 13281-13294 (1993). Feynman, J., Proton Fluenee Prediction Models, In Solar-Terrestrial Predictions-V, G. Heckman, K. Maribashi, M.A. Shea, D.F. Smart, and R. Thompson (Eds). 457-469 (1997). Foelsche, T., The Ionizing Radiation in Supersonic Transport Flights, in Second Symposium on Protection Against Radiation in Space, Edited by A. Reetz, Jr., NASA SP-71 (1965). Humble, J.E., M.L. Duldig, D.F. Smart, and M.A. Shea, Detection of 0.5-15 GeV Solar Protons on 29 September 1989 at Australian Stations, Geophys. Res. Letters, 18, 737-740 (199 I). ICRU Report 71, Conversion Coefficients for Use in Radiological Protection Against External Radiation - Joint report of an ICRU/ICRP task group, International Commission of Radiation Units and Measurements, Bethesda, Maryland, USA (1998). Kallenrode, M.-B., Shocks As Mechanisms For The Acceleration and Propagation of Energetic Particles, Adv. Space Res., 13(9), 342-350 (1993). Kallenrode, M.-B., G. Wibberenz, H. Kunow, R. Muller-Merlin, and N. Kontor, Multi-spacecraft Observations of Particle Events and Interplanetary Shocks During November-December 1982, Solar Phys., 147, 377-410 (1993). Lingenfelter, R.E.H.S. Hudson, In Proc. Conf. Ancient Sun, R.O. Pepin, J.A. Eddy, and R.B. Merrill (Eds.), Pergamon Press, New York, 69-79, (1980). Lulm, A., D. Hovestadt, B. Kleckler, M. Scholer, G. Gloeckler, F.M. Impavich, A.B. Galvin, C.Y. Fan, and L.A. Fisk, 19th International Cosmic Ray Conference, 4, 241-244, NASA Conference Publication 2376, Greenbelt, MD, (1985). Mason, G.M., The Composition of Galactic Cosmic Rays and Solar Energetic Particles, Rev. Geophys., 25, 685696, (1987). Mathews, T., and D. Venkatesan, Unique Series of Increases in Cosmic-ray Intensity due to Solar Flares, Nature, 345, 600-602 (1990). Mazur, J.E., G.M. Mason, B. Klecker, and R.E. MeGuire, The Energy Spectra of Solar Flare Hydrogen, Helium, Oxygen and Iron: Evidence For Stochastic Acceleration, Astrophysical Journal, 401, 398-410 (1992). Mewaldt, R A., and E.C. Stone, Isotope Abundances of Solar Coronal Material Derived from Solar Energetic Particle Measurements, Astrophysical Journal, 377, 959-963 (1989). Nymmik, R.A., Average Energy Spectra at Peak Flux and Fluence Values in Solar Cosmic Ray Events, 23rd International Cosmic Ray Conference, 3, Univ. Calgary, Calgary, Alberta, 29-32 (I 993). O'Brien, K., D.F. Smart, and H.H. Sauer, Atmospheric Cosmic Rays and Solar Energetic Particles at Aircraft Altitudes, Proceedings of the 1996 Topical Meeting on Radiation Protection and Shielding, 59-67 (1996).
1044
D.F. Smart and M. A. Shea
O'Brien, K., W. Friedberg, D.F. Smart, and H.H. Sauer, The Atmospheric Cosmic and Solar Energetic Particle Radiation Environment at Aircraft Altitudes, Advances in Space Research, 21(12), 1739-1748 (1998). Oetliker, M., B. Kleckler, D. Hovestadt, M. Scholer, J.B. Blake, M. Looper, and R.A. Mewaldt, Charge States of Heavy Solar Energetic Pm'ticles: Observations with the HILT Sensor on Sampex, 24th International Cosmic Ray Conference, 4, Universita Di Roma "La Sapienza", Rome, 470-474 (1995). Reames, D.V., Solar Energetic Particles: A Paradigm Sl~ft, Rev. Geophys. Supplement, 33, 585-589, 0995). Reames, D.V., J.P. Meyer, and T.T. Von Rosenvinge, Energetic Particle Abundances in Impulsive Solar Flare Events, Astrophys. J. Suppl., 90, 649-667 (1994). Reedy, R.C., Constraints on Solar Particle Events from Comparisons of Recent Events and Million-year Averages, In Solar Drivers Of Interplanetary and Terrestrial Disturbances, K.S. Balasubramanian, S.L. Keil, and R.N. Smartt (Eds.), Astronomical Society of the Pacific, Conference Series, Vol. 95. 429-436 (1996). Shea, M.A., and D.F. Smart, A Summary of Major Solar Proton Events, Solar Physics, 127, 297-320 (1990). Smart, D.F., and M.A. Shea, Comment on the Use of Solar Proton Spectra in Solar Proton Dose Calculations, 27th International Conference on Environmental Systems, SAE Technical Paper Series, 972568 (1997). Smart, D~F., and M.A. Shea, A Comparison of the Magnitude of the 29 September 1989 High Energy Event with Solar Cycles 17, 18, and 19 Events, 22nd International Cosmic Ray Conference, 3, Dublin Institute for Advanced Studies, Dublin, 101-104 (1991). Smart, D.F., M.A. Shea, Margaret D. Wilson, and L.C. Gentile, Solar Cosmic Rays on 29 September 1989; An Analysis Using the World-Wide Network of Cosmic Ray Stations, 22nd International Cosmic Ray Conference, Contributed Papers, 3, Dublin Institute for Advanced Studies, Dublin, 9%100 (1991). Smart, D.F., and M.A. Shea, The >10 Mev Peak Flux Distribution, In Solar-Terrestrial Predictions-V, G. Heckman, K. Maribashi, M.A. Shea, D.F. Smart, and R. Thompson (Eds). 449-452 (1997). Stassinopoulos, E.G., G.J. Brucker, D.W. Nakamura, C.A. Stauffer, G.B. Gee, and J.L. Barth, Solar Flare Proton Evaluation at Geostationary Orbits for Engineering Applications, IEEE Transactions on Nuclear Science, 43, 369-382 (1996). Swinson, D.B., and M.A. Shea, The September 29, 1989 Ground-Level Event Observed at High Rigidity, Ueophys. Res. Letters, 17, 1073-1076 (1990). Tylka, A.J., P.R. Boberg, J.H.Adams Jr., L.P. Beahm, W.F. Dietrich, and T. Kleis, The Mean Ionic Charge State of Solar Energetic Fe Ions above 200 MeV Per Nucleon, Astrophys. J., ~A, L109-L113 (1995). Tylka, AJ., W.F. Dietrich, and P.R. Boberg, Probability Distributions of High-Energy Solar-Heavy-Ion Fluxes from IMP-8: 1973-1996, IEEE Transactions on Nuclear Science, 44, 2140-2149 (1997). Van HoUebeke, M.A., L.S. Ma Sung, and F.B. McDonald, The Variation of Solar Proton Energy Spectra and Size Distribution with Heliolongitude, Solar Phys., 4 I, 189-323 (I 975). Wilson, J.W., L.W. Townsend, W.S. Schimmerling, G.S. Khandelwal, F. Khan, J. Nealy, F.A. Cucinotta, J.L. Shinn, L.C. Simonsen, and J.W. Norbury, Transport Methods and Interactions for Space Radiations, NASA RP-1257 (1991). Wilson, J.W., and J.E. Nealy, Model and Data Base for the Background Radiation Exposure of High Altitude Aircraft, in Proceedings of the Topical Meeting on New Horizons in Radiation Protection and Shielding, American Nuclear Society, Inc., La Grange Park, Illinois, 425-441 (1992). Xapsos, M.A., G.P. Summers, P. Shapiro, and E.A. Burke, New Techniques for Predicting Solar Proton Events, IEEE Trans. Nuc. Sci., 43, 2772-2777 (1996).