Early Alert of Solar Radiation Hazard Paul Evenson (firstname.lastname@example.org), T. Kuwabara, J. W. Bieber, J. Clem, R. Pyle Bartol Research Institute and Department of Physics and Astronomy, University of Delaware.
Paul Evenson(email@example.com), T. Kuwabara, J. W. Bieber, J. Clem, R. Pyle
Bartol Research Institute and Department of Physics and Astronomy,University of Delaware
Abstract We describe our automated Ground Level Enhancement (GLE) Alarm that is now operational in beta testing. Using our neutron monitor stations at Fort Smith, Inuvik, McMurdo, Nain, Newark, Peawanuck, South Pole, and Thule, automated e-mail alarms can be sent to subscribers within minutes of large GLE onsets. In brief, the system generates a WATCH when one station observes an increase of 4% or more in a three minute moving average, a WARNING when two stations observe such an increase, and an ALERT when three or more stations record such an increase. (A higher threshold for percent increase may be imposed for some stations with lower count rates.) A scientific discussion of the alarm system can be found in the article "Development of a Ground Level Enhancement Alarm System Based upon Neutron Monitors", T. Kuwabara, J. W. Bieber, J. Clem, P. Evenson, and R. Pyle, Space Weather, 4, S10001, doi: 10.1029/2006SW000223, 2006. A web page specifically devoted to the GLE alarm, including graphical displays is available at: http://www.bartol.udel.edu/~takao/neutronm/glealarm/index.html
Algorithm for issuing an alarm
We define three levels of alarm (Watch, Warning, and Alert) based on the number of stations that record a significant intensity increase.
A threshold level Ithis set for the cosmic ray intensity increase, and an alarm is generated when the number of stations that exceed the threshold level is
1 - Watch, 2 - Warning, 3 –Alert.
Baseline to calculate the intensity increase is defined by two parameters tband t0. Then, by using a trailing moving average value for the current count rate, the intensity at time t = tis calculated each minute from the observed count rate N(t) averaged over the preceding tcminutes, shown as left equation.
To reduce fluctuations and accurately detect the GLE, all parameters (threshold and baseline) were optimized as by backtesting against past neutron monitor data as
tc=3 min,tb=75min, t0=10min, Ith = 4%.
At this condition, the false alarm rate for Watch, Warning, and Alertwas ~40/yr, less than 1/yr, and 0/yr respectively, during the 4.4-year period of our backtesting study,
Condition for issuing three levels of alarm. Intensity increase recorded at three stations during a schematic GLE are illustrated.
tc : average time
tb : duration of the
t0 : time interval between the
baseline and current time
Ith : threshold level
GLE event on January 20, 2005
Alarm issued for January 20, 2005 event
On Jan 20, 2005 event, alarm times for Warning and Alert are 6:49 and 6:50 respectively; a separate Watch alarm was not generated. The SEC alert was issued at 7:02 from >10 MeV data, and 7:04 from >100 MeV data.
Proton flux observed in the >10 MeV channel exceeds 100 pfu at 6:55. This flux level is defined as a “Moderate Storm” (S2 on NOAA Space Weather Scale for Solar Radiation Storms) that has the possibility to cause single-event upsets aboard satellites, and our system can produce alarm before this.
Space Weather Message Code: ALTPX2
Serial Number: 27
Issue Time: 2005 Jan 20 0702 UTC
ALERT: Proton Event 10MeV Integral Flux exceeded 100pfu
Begin Time: 2005 Jan 20 0701 UTC
NOAA Scale: S2 - Moderate
GOES proton flux and cosmic ray intensity increase. In upper panel, black arrows indicate the onset time of SEP event(>100 MeV), and the start time of the Moderate storm. Vertical lines drawn in lower panels also show the onset time of GLE. Colored arrows indicate the time when each alarm is issued (proton monitor) or generated (neutron monitor).
Space Weather Message Code: ALTPC0
Serial Number: 17
Issue Time: 2005 Jan 20 0704 UTC
ALERT: Proton Event 100MeV Integral Flux exceeded 1pfu
Begin Time: 2005 Jan 20 0701 UTC
Comparison between the alert times from our system
and alarm issue times from proton data at SEC/NOAA for nine GLE events.
GLE event on January 20, 2005. Upper panel shows the low energy proton integral flux recorded by the GOES-11 satellite (solid line >10MeV, dotted line >100MeV). Lower panel shows neutron rates detected in several neutron monitors. These data are normalized to average counting rate at Inuvik (05:30-06:30).
The GLE (Ground Level Enhancement) of January 20, 2005 was the largest in half a century. The onset of the intensity increase in neutron monitors is earlier than that of the low energy proton flux, and time to reach maximum intensity is shorter for neutron monitors.
The GLE particles have large mean free paths and travel almost at the speed of light, and can be detected with high accuracy by detectors with large volume such as ground-based neutron monitors. A GLE alarm can provide a very useful early warning of an impending solar radiation storm. We developed a system that watches (in real time) for count rate increases in our data, and gives an alarm when a GLE is detected. We compare the alarm time produced in our system with that of the system operated by SEC/NOAA by using the GOES proton data.
Number of Minutes by which each alert precedes Moderate Storm (>100pfu in >10 MeV channel)
Number of Minutes by which GLE Alert precedes earliest SEC Proton Alert
Eight neutron monitors are used in this work. The station name, detector type, average cosmic ray background count rate (in 2005),
geographical latitude, longitude, and altitude are listed.
*Bare : Neutron counters without the usual lead shielding that respond to a slightly lower energy primary cosmic ray than the standard monitor.