GEOPHYSICAL RESEARCH LETTERS, VOL. ???, XXXX, DOI:10.1029/,
1
2
The spectral dependence of terrestrial gamma-ray flashes on source distance 1
1
1
2
B. J. Hazelton, B. W. Grefenstette, D. M. Smith, J. R. Dwyer, X.-M. 3
4
5
6
6
Shao, S. A. Cummer, T. Chronis, E. H. Lay, and R. H. Holzworth
B. J. Hazelton, Department of Physics and Santa Cruz Institute for Particle Physics, University of California, Santa Cruz, 1156 High St., Santa Cruz, CA 95064, USA. (
[email protected]) 1
Department of Physics and Santa Cruz
D R A F T
August 19, 2008, 11:21am
D R A F T
X -2 3
HAZELTON ET AL: TGF SOURCE DISTANCE
Lightning flash locations from the World Wide Lightning Location Net-
4
work (WWLLN) were used to identify storm locations near 362 Terrestrial
5
Gamma-ray Flashes (TGFs) detected by RHESSI from October 2003 to De-
6
cember 2005. The combined spectrum of TGFs with lightning flashes within
Institute for Particle Physics, University of California, Santa Cruz, California, USA. 2
Department of Physics and Space
Sciences, Florida Institute of Technology, Florida, USA. 3
Space and Remote Sensing Sciences, Los
Alamos National Laboratory, Los Alamos, New Mexico, USA. 4
Electrical and Computer engineering
Department, Duke University, Durham, North Carolina, USA. 5
Global Hydrology Climate Center,
NASA Marshall Space Flight Center, Huntsville, Alabama, USA. 6
Department of Earth and Space Science,
University of Washington, Washington, USA.
DRAFT
August 19, 2008, 11:21am
DRAFT
HAZELTON ET AL: TGF SOURCE DISTANCE 7
300 km of the sub-satellite point is found to be much harder than the spec-
8
trum of TGFs with more distant lightning flashes. When these data are com-
9
pared with simulations of relativistic runaway breakdown, it is found that
10
the most likely model has a source altitude of 15 km and a wide beam ge-
11
ometry. Four associations of TGFs with individual sferics localized to po-
12
sitions more than 300 km from the sub-satellite point are reported and it is
13
shown that a narrow beam source at 21 km altitude is unlikely to produce
14
the number of observed high energy photons in these TGFs.
D R A F T
August 19, 2008, 11:21am
X-3
D R A F T
X-4
HAZELTON ET AL: TGF SOURCE DISTANCE
1. Introduction 15
TGFs were first discovered with the BATSE instrument on the Compton Gamma Ray
16
Observatory by Fishman et al. [1994] and over 800 have since been observed by the
17
RHESSI satellite [Smith et al., 2005]. From their first detection, TGFs have been as-
18
sociated with thunderstorms, and some flashes from each satellite have been linked to
19
individual lightning flashes [e.g. Inan et al., 1996, 2006; Cummer et al., 2005; Stanley
20
et al., 2006]. Localizations of lightning flashes associated with RHESSI TGFs indicate
21
that most TGFs occur within ∼300 km of the sub-satellite point [Cummer et al., 2005].
22
Relativistic runaway models have been successful in modeling the combined spectrum
23
of all RHESSI TGFs [Dwyer and Smith, 2005], but details of source beaming and alti-
24
tude cannot be constrained by the combined spectrum because it averages over important
25
differences among individual TGFs. One such difference is the horizontal distance from
26
the TGF source to the sub-satellite point. As this distance increases, the spectrum of
27
photons reaching the satellite is expected to soften (decrease in average photon energy)
28
because the Bremsstrahlung radiation is inherently harder along the beam and Compton
29
scattering of photons out of the beam also removes energy [Østgaard et al., 2008]. The
30
change in the detected TGF spectrum with source distance is model dependent, so the
31
spectral differences between TGFs that are detected close and far from the source can
32
constrain source geometry and altitude. Unfortunately, the spectra of individual TGFs
33
associated with localized sferics cannot be directly fit with model spectra because RHESSI
34
TGFs have fewer than 30 counts on average.
35
To obtain separate spectra of TGFs with close and distant sources, we used the World
D R A F T
August 19, 2008, 11:21am
D R A F T
HAZELTON ET AL: TGF SOURCE DISTANCE
X-5
36
Wide Lightning Location Network (WWLLN) to identify storm locations around the world
37
at the times when TGFs were detected. These storm locations are used to indicate the
38
possible locations of TGF sources and to categorize TGFs based on the distance from the
39
sub-satellite point to the nearest potential source location. WWLLN is a multi-station
40
VLF network that records times and locations of lightning flashes all over the world with
41
a mean accuracy of 15 km [Rodger et al., 2008]. Jacobson et al. [2006] showed that despite
42
its low detection efficiency for a given flash, WWLLN provides a spatially accurate and
43
complete census of storms when averaged over many flashes.
2. Data and Analysis 44
We consider 362 TGFs detected by RHESSI between October 2003 and December 2005
45
for which concurrent WWLLN data are available. This time range is chosen to start
46
after the time of group arrival (TOGA) algorithm was implemented in WWLLN [Rodger
47
et al., 2005] and to end before radiation damage to the RHESSI detectors affected the
48
spectral quality too much (Grefenstette et al., 2008b, in preparation). For each TGF, all
49
WWLLN lightning flashes within ± 20 minutes of the TGF time were accumulated to
50
locate lightning producing storms. The TGFs were divided into two categories: those that
51
had lightning flashes within 300 km of the sub-satellite point (316 TGFs), and those that
52
did not (46 TGFs). To evaluate the spectral differences between TGFs in each category,
53
the combined spectrum for each category was divided by the combined spectrum of all
54
362 TGFs in the sample. The spectrum of TGFs with only distant lightning flashes is
55
much softer than the spectrum of TGFs with closer ones (Figure 1).
D R A F T
August 19, 2008, 11:21am
D R A F T
X-6
HAZELTON ET AL: TGF SOURCE DISTANCE
3. Runaway Breakdown Models 56
A Monte Carlo simulation of runaway breakdown in air similar to the one used in Dwyer
57
and Smith [2005] was used to calculate the spectra of TGFs originating more than and
58
less than 300 km from the sub-satellite point. The simulation includes, in an accurate
59
form, all the important interactions of energetic electrons, positrons and photons in air.
60
A detailed description of the included physics is in Dwyer [2007].
61
In the first phase of the simulation, the runaway avalanche was allowed to develop in a
62
region with a uniform vertical electric field. Outside this region the field was set to zero.
63
The runaway electrons were propagated until they left the avalanche region and came to a
64
stop. The solid angle distribution of the photons produced in this phase of the simulation
65
(Figure 2a) shows that the photons are beamed upward, with the higher energy photons
66
in a narrower cone than the lower energy photons. The photons were then propagated
67
through the atmosphere using GEANT, a standard high-energy particle transport code
68
used in particle physics and astrophysics [GEANT Team, 1993].
69
We also considered a wider beam that might result from divergence in the electric field
70
in the avalanche region. The wider distribution is generated by convolving the original
71
distribution with a gaussian in solid angle. This convolution widens the beam while
72
preserving the energy structure of the original beam (Figure 2b).
73
To compare the models with data, it is important to consider the intensity threshold
74
for detecting the events. As the photons propagate through the atmosphere they are
75
sometimes Compton scattered to very large angles. These Compton scattered photons
76
have a much softer spectrum than the photons in the main part of the beam, so they
D R A F T
August 19, 2008, 11:21am
D R A F T
HAZELTON ET AL: TGF SOURCE DISTANCE
X-7
77
can have a significant effect on the model spectrum. Because they are scattered into a
78
large solid angle, the satellite is more likely to fly into the path of the Compton scattered
79
photons than the main beam. The flux of Comptonized photons is much lower than the
80
flux in the beam, however, so the search algorithm is less likely to identify them as a
81
TGF. In Dwyer and Smith [2005] the model spectra only included emission into angles
82
at which the flash would appear at least half as bright to the satellite as a flash directly
83
below it. This threshold was based on the range of observed fluences in RHESSI TGFs
84
and the assumption that all TGFs have the same intrinsic brightness, so the observed
85
fluence depended only on source distance. A more sophisticated threshold was developed
86
for this paper using the knowledge of source distance provided by WWLLN.
87
First, the number of counts detected by RHESSI from the total of all 362 TGFs was
88
calculated as a function of distance from the sub-satellite point using the distance from the
89
sub-satellite point to the nearest lightning flash for each TGF. Next, the model photons
90
were binned by horizontal distance from the source and the spectrum from each distance
91
bin was convolved with the instrument response to generate the count spectrum that
92
would be detected by RHESSI. Finally, the number of model counts in each distance bin
93
was re-normalized to the number of counts actually detected by RHESSI in that distance
94
bin. One major advantage of this method is that any deadtime effect on the number
95
of detected counts versus distance is duplicated in the models (see section 5). There is
96
some systematic bias in this re-normalization method for placing the source too near the
97
RHESSI sub-satellite point because the horizontal distance to the nearest lightning flash
98
was used for the source distance. The result of this bias is that the calculated model
D R A F T
August 19, 2008, 11:21am
D R A F T
X-8 99
HAZELTON ET AL: TGF SOURCE DISTANCE
spectra are harder than they would be if the actual source distances were known.
100
Three source altitudes (13 km, 15 km, and 21 km) and two types of beaming (narrow and
101
wide) were considered. The models that best fit the combined spectra of all 362 TGFs are
102
the 21 km narrow model and 15 km wide model, as in Dwyer and Smith [2005], except that
103
the wide beam used in that paper was isotropic in a 45◦ half-angle cone. To compare the
104
models with the TGFs separated by source distance, the distance-binned model spectra
105
were binned further to make two spectra for each model: one for TGFs within 300 km of
106
the source and one for TGFs further from the source (Figure 3). The TGFs with close
107
and distant sources are best fit by the 21 km narrow and 21 km wide models respectively;
108
the 15 km wide model is the second-best fit to both categories. Unfortunately, no single
109
model fits all the data perfectly, but the 15 km wide model comes the closest, particularly
110
if all the models are a little too hard as a result of the distance bias mentioned above.
111
To quantify the effect of this bias, the same re-normalizing procedure was carried out,
112
except that the source distance was defined for each TGF using a randomly selected
113
lightning flash that occurred within 20 minutes and 600 km of the TGF, rather than
114
the nearest one. This procedure was performed many times and the results provide an
115
upper limit on how much the models might change if the actual source locations were
116
known, since randomly selecting a lightning flash as the source location will generally
117
overestimate the source distance. This algorithm produced models that were all softer
118
than when the nearest lightning flashes were used, with the largest changes occurring in
119
the narrow models in Figure 3a. The 15 km wide beam model was the best fit for the TGFs
120
with close and distant sources, and it was the second best fit for the combined spectrum
D R A F T
August 19, 2008, 11:21am
D R A F T
HAZELTON ET AL: TGF SOURCE DISTANCE
X-9
121
of all the TGFs, after the 13 km wide model. We conclude that the bias introduced by
122
using the nearest lightning flash affects the narrow models the most and that the effect
123
of removing the bias would be to make the 15 km wide beam model an even better fit to
124
all the data.
4. TGFs with Distant Sferic Localizations 125
We report here on four sferics associated with TGFs and localized to positions more
126
than 300 km from the sub-satellite point (listed in Table 1). Three of these sferics were
127
detected and localized by the Los Alamos Sferic Array (LASA) following the analysis
128
procedures outlined in Stanley et al. [2006] and the fourth was detected and localized
129
both by the Zeus network [Chronis and Anagnostou, 2003] and by magnetic field sensors
130
at Duke University [after Cummer et al., 2005]. The highest energy gamma-rays in these
131
distant localized TGFs provide a useful check on the viability of narrow beam models.
132
A bootstrap Monte Carlo method was used to compare the number of high energy counts
133
observed in each TGF with the number that would be expected from a narrow beam source
134
at 21 km. First, model spectra were generated in 100 km wide distance bins and convolved
135
with the instrument response. Then for each localized TGF, the model spectrum for the
136
appropriate distance bin was combined with the background spectrum for that TGF taken
137
from the RHESSI TGF Catalog (Grefenstette et al., 2008b, in preparation) to generate
138
an expected observed spectrum for the bootstrap Monte Carlo. An energy threshold that
139
contained only the top 5% of model counts was selected and the number of counts in the
140
real TGF with energies above this threshold was compared with the number of above
141
threshold counts in ten million simulated TGFs. The probability of detecting a TGF with
D R A F T
August 19, 2008, 11:21am
D R A F T
X - 10
HAZELTON ET AL: TGF SOURCE DISTANCE
142
at least as many high energy counts as were observed if the source was a narrow beam
143
at 21 km is shown in Table 1. Also shown in Table 1 is the joint probability of detecting
144
at least the number of observed high energy counts in all four localized TGFs. From the
145
probabilities in the table, the likelihood that the 21 km narrow model represents the true
146
source for all TGFs is very low.
147
The radiation damage to the RHESSI detectors has degraded the spectral resolution
148
and decreased the average energy of detected counts in RHESSI TGFs over time. The
149
instrument response used in this analysis is therefore not the true RHESSI response for
150
the last three TGFs in Table 1. The decrease in average TGF count energy with time
151
means that the bootstrap Monte Carlo actually overestimates the probability of detecting
152
high energy counts from the 21 km narrow model, so the true observation probabilities
153
for these TGFs are even lower than reported in Table 1.
5. Discussion 154
Estimating TGF source locations using sferic data has proved a powerful tool to ex-
155
amine models of TGF production. The distance from the estimated source location to
156
the sub-satellite point effectively sorts TGFs into categories with very different spectra,
157
something no other parameter has ever been observed to do (Grefenstette et al., 2008b, in
158
preparation). Information about the source distance for TGFs increases the constraints
159
on the models and provides an improved model threshold algorithm. The results of the
160
modeling are broadly consistent with Dwyer and Smith [2005] and suggest that the 15 km
161
wide beam model is the most likely source geometry, but none of the models considered
162
here are completely successful. The failure of any of the models to fully explain the spec-
D R A F T
August 19, 2008, 11:21am
D R A F T
HAZELTON ET AL: TGF SOURCE DISTANCE
X - 11
163
tra of both TGFs with close sources and TGFs with distant sources suggests that more
164
complicated models, perhaps with non-vertical beams, will need to be considered in the
165
future.
166
The viability of models with vertical narrow beams can be constrained by the number of
167
high energy photons in individual TGFs with associated localized sferics that have source
168
distances greater than 300 km. We find that a 21 km altitude source with the beam profile
169
shown in Figure 2a is unlikely to produce the number of high energy counts observed in
170
the four such TGFs examined here, so models considered in the literature [e.g Carlson
171
et al., 2007; Østgaard et al., 2008] with even narrower beams and higher altitudes are
172
unlikely to produce distant TGFs with the number of observed high energy gamma-rays.
173
It was recently discovered (Grefenstette et al., 2008b, in preparation) that the brightest
174
RHESSI TGFs are likely to be suffering from saturation effects. There are two saturation
175
effects that will change the detected spectrum of TGFs with close sources: deadtime and
176
pileup. Deadtime occurs when a photon enters a RHESSI detector less than 9 µs after
177
a previous one and is not counted because the electronics have not finished processing
178
the previous count [Smith et al., 2002]. Deadtime does not directly change the observed
179
spectrum, but because the brightest part of the TGF is also the hardest [see Grefenstette
180
et al., 2008], the detected spectrum for TGFs affected by deadtime will be softer than
181
the true spectrum. Pileup occurs when two photons enter a RHESSI detector less than
182
1 µs apart and the detector sees only one event with an energy equal to the sum of the
183
energies of the two photons, making the observed spectrum harder. To quantify the spec-
184
tral effects of deadtime and pileup, we ran extensive Monte Carlo simulations similar to
D R A F T
August 19, 2008, 11:21am
D R A F T
X - 12
HAZELTON ET AL: TGF SOURCE DISTANCE
185
those discussed in section 3 with typical TGF time profiles and 50% peak deadtime and
186
compared the resulting spectra to the spectra without saturation effects. We find that
187
above 5 MeV, deadtime results in a 5% decrease in the number of observed counts, while
188
pileup causes a 20% increase, for a net increase in high energy counts of about 15% due
189
to saturation effects. By contrast, the difference between the two best models in Figure
190
3a is more than 30%. The analysis in section 4 and the data in Figure 3b are unlikely to
191
be affected by deadtime or pileup because TGFs with distant sources are expected to be
192
the least saturated.
193
Acknowledgments. We wish to thank Craig Rodger for very useful conversations.
References 194
Carlson, B. E., N. G. Lehtinen, and U. S. Inan (2007), Constraints on terrestrial gamma
195
ray flash production from satellite observation, Geophys. Res. Lett., 34, L08809, doi:
196
10.1029/2006GL029229.
197
Chronis, T. G., and E. N. Anagnostou (2003), Error analysis of a long-range lightning
198
monioting network of ground-based receivers in europe, J. Geophys. Res., 108 (D24),
199
4779, doi:0.1029/2003JD003776.
200
Cummer, S. A., Y. Shai, W. Hu, D. M. Smith, L. I. Lopez, and M. A. Stanley (2005),
201
Measurements and implications of the relationship between lightning and terrestrial
202
gamma ray flashes, Geophys. Res. Lett., 32, L08811, doi:10.1029/2005GL022778.
203
204
Dwyer, J. R. (2007), Relativistic breakdown in planetary atmospheres, Phys. Plasmas, 14, 042901, doi:10.1063/1.2709652.
D R A F T
August 19, 2008, 11:21am
D R A F T
HAZELTON ET AL: TGF SOURCE DISTANCE
X - 13
205
Dwyer, J. R., and D. M. Smith (2005), A comparison between Monte Carlo simulations of
206
runaway breakdown and terrestrial gamma-ray flash observations, Geophys. Res. Lett.,
207
32, L22804, doi:10.1029/2005GL023848.
208
209
210
211
Fishman, G. J., et al. (1994), Discovery of intense gamma-ray flashes of atmospheric origin, Science, 264 (5163), 1313–1316. GEANT Team (1993), GEANT - detector description and simulation tool, CERN program library long write-up W5013, Tech. rep., CERN, Geneva.
212
Grefenstette, B. W., D. M. Smith, J. R. Dwyer, and G. J. Fishman (2008), Time
213
evolution of terrestrial gamma ray flashes, Geophys. Res. Lett., 35, L06802, doi:
214
10.1029/2007GL032922.
215
Inan, U. S., S. C. Reising, G. J. Fishman, and J. M. Horack (1996), On the association of
216
terrestrial gamma-ray bursts with lightning and implications for sprites, Geophys. Res.
217
Lett., 23 (9), 1017–1020.
218
Inan, U. S., M. B. Cohen, R. K. Said, D. M. Smith, and L. I. Lopez (2006), Terres-
219
trial gamma ray flashes and lightning discharges, Geophys. Res. Lett., 33, L18802, doi:
220
10.1029/2006GL027085.
221
Jacobson, A. R., R. H. Holzworth, J. Harlin, R. L. Dowden, and E. H. Lay (2006),
222
Performance assessment of the World Wide Lightning Location Network (WWLLN)
223
using the Los Alamos Sferic Array (LASA) as ground truth, J. Atmos. and Oceanic
224
Tech., 23, 1082–1092.
225
Østgaard, N., T. Gjesteland, J. Stadsnes, P. H. Connell, and B. E. Carlson (2008),
226
Production altitude and time delays of the terrestrial gamma flashes: Revisiting the
D R A F T
August 19, 2008, 11:21am
D R A F T
X - 14
HAZELTON ET AL: TGF SOURCE DISTANCE
227
Burst and Transient Source Experiment spectra, J. Geophys. Res., 113, A02307, doi:
228
10.1029/2007JA012618.
229
Rodger, C. J., J. B. Brundell, and R. L. Dowden (2005), Location accuracy of VLF World-
230
Wide Lightning Location (WWLL) network: post-algorithm upgrade, Ann. Geophys.,
231
23, 277–290.
232
Rodger, C. J., J. B. Brundell, R. H. Holzworth, E. H. Lay, and R. L. Dowden (2008), Im-
233
provements in the WWLLN network: Growing detection efficiencies for ”big lightning”
234
events, paper presented at the Workshop on Coupling of Thunderstorms and Lightning
235
Discharges to Near-Earth Space, Corte, France, 23–27 June.
236
Smith, D. M., L. I. Lopez, R. P. Lin, and C. P. Barrington-Leigh (2005), Ter-
237
restrial gamma-ray flashes observed up to 20 MeV, Science, 307, 1085, doi:
238
10.1126/science.1107466.
239
240
Smith, D. M., et al. (2002), The RHESSI spectrometer, Sol. Phys., 210 (1), 33–60, doi: 10.1023/A:1022400716414.
241
Stanley, M. A., X.-M. Shao, D. M. Smith, L. I. Lopez, M. B. Pongratz, J. Harlin, M. Stock,
242
and A. Regan (2006), A link between terrestrial gamma-ray flashes and intracloud
243
lightning discharges, Geophys. Res. Lett., 33, L06803, doi:10.1029/2005GL025537.
D R A F T
August 19, 2008, 11:21am
D R A F T
X - 15
HAZELTON ET AL: TGF SOURCE DISTANCE
Table 1.
TGFs with sferic localizations more than 300 km from the sub-satellite point.
The last column is the probability of observing a TGF with at least the observed number of high energy counts if the source was a narrow beam at 21 km. The last row shows the joint probability for all four TGFs (see text for details).
Date August 4, 2004
Detection Network
Source Total High Energy Observation Distance (km) Counts Counts Probability
Duke & Zeus
535
18
2
0.22
September 11, 2006
LASA
371
22
5
0.0034
June 11, 2007
LASA
319
22
2
0.28
LASA
374
19
2
0.23
June 16, 2007
4.8 × 10−5
All Events
D R A F T
August 19, 2008, 11:21am
D R A F T
X - 16
Figure 1.
HAZELTON ET AL: TGF SOURCE DISTANCE
Ratios of the spectra of TGFs with close lightning flashes (red) and TGFs with no
lightning flashes within 300 km (blue) to the combined spectrum of all 362 TGFs in the sample. The black line indicates the combined spectrum of all TGFs in the sample. TGFs with only distant lightning flashes have a much softer spectrum than TGFs with close lightning flashes.
D R A F T
August 19, 2008, 11:21am
D R A F T
HAZELTON ET AL: TGF SOURCE DISTANCE
Figure 2.
X - 17
(a) Intrinsic beam shape from relativistic runaway simulations. Note the energy
structure of the beam, with the high energy photons concentrated in a narrower emission cone than the low energy photons. (b) Wide beam shape from the convolution of the intrinsic beam shape with a gaussian in solid angle. The energy structure of the intrinsic beam is preserved.
D R A F T
August 19, 2008, 11:21am
D R A F T
X - 18
Figure 3.
HAZELTON ET AL: TGF SOURCE DISTANCE
(a) The ratio of the spectrum of TGFs with close lightning flashes (316 TGFs) to
the combined spectrum of all 362 TGFs plotted with the runaway breakdown models for close sources. The 21 km narrow model is the best fit, followed by the 15 km wide model. (b) The ratio of the spectrum of TGFs with only distant lightning flashes (46 TGFs) to the combined spectrum of all 362 TGFs plotted with the runaway breakdown models for distant sources. The 21 km wide model is the best fit to the data, followed by the 15 km wide model.
D R A F T
August 19, 2008, 11:21am
D R A F T