Lightning superbolts are the most powerful and rare lightning events with intense optical emission, first identified from space. Superbolt events occurred in 2010-2018 could be localized by extracting the high energy tail of the lightning stroke signals measured by the very low frequency ground stations of the World-Wide Lightning Location Network. Here, we report electromagnetic observations of superbolts from space using Van Allen Probes satellite measurements, and ground measurements, and with two events measured both from ground and space. From burst-triggered measurements, we compute electric and magnetic power spectral density for very low frequency waves driven by superbolts, both on Earth and transmitted into space, demonstrating that superbolts transmit 10-1000 times more powerful very low frequency waves into space than typical strokes and revealing that their extreme nature is observed in space. We find several properties of superbolts that notably differ from most lightning flashes; a more symmetric first ground-wave peak due to a longer rise time, larger peak current, weaker decay of electromagnetic power density in space with distance, and a power mostly confined in the very low frequency range. Their signal is absent in space during day times and is received with a long-time delay on the Van Allen Probes. These results have implications for our understanding of lightning and superbolts, for ionosphere-magnetosphere wave transmission, wave propagation in space, and remote sensing of extreme events.
Lightning superbolts are rare and extreme events that were first identified from optical stroke data measured by a photometer on board the Vela satellites1,2, yielding between 1011 and 1013 W per stroke. With such high radiated power, the temperature in the core channel of the stroke must exceed the commonly accepted maximum temperature of the lightning return stroke3,4 (i.e., ~3 × 105 K)5, creating a debate concerning our current understanding of the energy balance in a lightning discharge3,6. Recently, radio frequency (RF) superbolts were geographically localized7 using the tail of the occurrence distribution in lightning energy8, defined as above 1 MJ (1000 times greater than the 1 kJ mean), measured by the very low frequency (VLF) ground stations of the World-Wide Lightning Location Network (WWLLN)9,10,11,12. Interestingly, the distribution of superbolt locations and occurrence times was not equivalent to that of ordinary lightning: instead, superbolts were found to occur over oceans and seas at a much higher rate, and more often in winter7. The north Atlantic (west of Europe) and the Mediterranean Sea have some of the highest wintertime occurrence rates of superbolts7 (see Fig. 1 and methods subsection Ground-based measurements of superbolts). As with any cloud-to-ground (CG) lightning flash, superbolts emit electromagnetic radiation in the very low frequency band that propagates within the Earth-ionosphere waveguide and escapes to the magnetosphere as whistler-mode waves along Earth’s magnetic field lines13.
Here, we show superbolt VLF electromagnetic (EM) power density in space. We combine space and ground-based measurements in a unique manner to follow electromagnetic superbolt signals from Earth to space over thousands of kilometers, to widely characterize their VLF electric and magnetic wave power density in space and on Earth, to compute ground-space transmitted power ratio, and to extract statistical electromagnetic properties of lightning superbolts never before reported. We conclude that, in addition to a location bias, superbolts exhibit other properties that differ from ordinary lightning, deepening the mystery associated with these extreme events.
Superbolt electromagnetic power
Superbolts are first located world-wide by making use of the WWLLN database, with an identification of superbolts from their energy greater than 1 MJ from synchronized measurements by WWLLN ground receivers7. In Europe, we make use of electric field measurements from a ground measurement campaign (ECLAIR)14,15 conducted by CEA as well as data recorded by the stations of the French lightning location network, Météorage (MTRG)16 (see methods subsection Ground-based measurements of superbolts). In space, we use two instruments17,18 on board the NASA Van Allen Probes mission19 to identify superbolt electric and magnetic wave fields (see methods subsection Space measurements of superbolts) from high-definition recordings. Superbolts signals are then gathered and retained after a selective screening process (see methods subsection Superbolt detection and selection). High-resolution electric power spectral density (PSD) from a selection of 6 superbolts is shown in Fig. 2, from space (a1–d1) and ground-based (e1–f1) measurements. The black line in Fig. 2 (a2–f2) is the wave electric (magnetic) intensity, E2 in V2/m2 (B2 in pT2), more commonly called, the wave power, which computation from the PSD is explained in the method section. In space, burst-mode spectrograms (Fig. 2 a1–d1) show superbolt VLF waves have a clear whistler-mode descending tone over a time period of seconds, as with all lightning-generated whistlers20, due to frequency-dependent dispersion through the ionosphere and plasmasphere. Their PSD frequency range spans ~0.1 to ~10 kHz, with most power above 2 kHz21 and within the first second. The sharp rising tone just prior to the descending tone shape (e.g., Fig. 2 d1 for frequencies within ~102 and ~104 Hz at t ~ 0.4) is an anti-aliasing filter that should be disregarded. Superbolts are observed in space with a time delay of 0.1 to 0.4 s from the superbolt WWLLN recorded time that marks t = 0. The delay corresponds to the time the wave takes to propagate away from the stroke location, along the Earth-ionosphere waveguide, and then through the ionosphere into the magnetosphere along field lines and ultimately to the satellite. In some cases, the wave propagates to the magnetically conjugate footprint and reflects back before being seen at the spacecraft (see discussion below related to the superbolt of Fig. 3). Event times of all WWLLN events in the satellite burst window are marked with dashed vertical lines in Fig. 2 (bottom line plots, with the superbolt time indicated in red and non-superbolt lightning signals in green). An empirical estimate of the median squared electric field at the satellite location is calculated and reported for each of these strokes22,23 based on both the estimated WWLLN total energy on the ground and the distance of the event (circles in the bottom line plots). The non-superbolt median squared electric field estimates (green circles) are lower than the superbolt estimates (red circle), by >3 orders of magnitude, confirming that the detected space-based signals only correspond to the WWLLN superbolt, and are not caused by other lightning.
Simultaneous ground-based and space measurement of a superbolt
The closest of the two superbolts (1.2 MJ) that we observe synchronously from Earth and space is represented in Fig. 3 (more detail in Supplementary Fig. 6 and Supplementary Fig. 7). The waveforms (3c and 3f) are decomposed by Fourier transform to produce electric field (3a) and magnetic field (Supplementary Fig. 6e) PSDs. Other lightning activity (Fig. 3c) is composed of nine lightning events identified with WWLLN, listed in Supplementary Table 4, and plotted on a map in Supplementary Fig. 7f. Among these nine normal lightning events, event #3 is a strong 95.8 kJ stroke occurring ~2 s later. This lightning activity is also visible in space with whistler-mode waves (3a), showing good temporal alignment with the WWLLN signal and representing significant wave power (line plot in 3b and WWLLN estimated power) for times between 2 and 3 s. Nevertheless, this activity does not mask or perturb the superbolt burst-mode signal which persists for ~1 s of the first ~1.4 s of data. The superbolt signal is found to be composed of one main wave (with a 0.4 s delay from WWLLN detection time) and a second similarly powerful wave 0.2 s later (3a, 3c). The time it takes for the wave to travel from the satellite to the magnetically conjugate point along the field line at L ~ 2.43, reflect, and return at the equator is 0.2 s. Computing the Poynting flux (see Supplementary Method 1), we find a change of sign between these two waves, which indicates the second wave reaches the spacecraft from the opposite side of the field line. In the absence of any other event attested by WWLLN and since the time delay is consistent with one reflection, this suggests the two superbolt’s waves have entered the magnetosphere at different latitudes surrounding the source, followed then different paths, with one wave reaching first the spacecraft and the second reaching the spacecraft 0.2 s later after bouncing at the conjugate point of the field line. This is one of the plausible and classic scenarios obtained from ray-tracing simulations24. The peak squared electric field measured by an ECLAIR ground station for this event is 20 (V/m)2 1387 km (3e) from the source. The average squared electric field over its duration (1.5 ms) is 4 (V/m)2, about 400 times greater than cloud-to-ground lightning measured on the ground (cf. discussion below). Quantifying the power spectral density on the ground is critical for validating ray-tracing simulations that predict the wave propagation from the source to space. The ground electric field peak (Fig. 3f and Supplementary Fig. 7e) is well resolved (at 0.08 µs from the 12.5 MHz acquisition) and presents a longer-than-average rise time (the time for the signal to go from 50% to 90% of peak magnitude), which is confirmed to be specific to superbolts from rise time and fall time statistics below. This leads to a somewhat more symmetric first peak. In space, the northern magnetic footprint of the Van Allen Probes North (3c) is 3193 km away from the superbolt location and the measured peak squared electric field is 6 × 10−7 (V/m)2, for a mean value of 8.5 × 10−8 (V/m)2 over 1 s, about 200 times higher than for normal lightning21. The root mean-squared (rms) amplitude of the magnetic field (3c), essential to compute pitch angle scattering of radiation belt electrons caused by lightning-generated whistlers in space (see review ref. 25), is 19 pT, 19 times larger than for normal lightning21. However, this is not as large as some of the superbolts we discuss next.
Propagation and attenuation of superbolt electromagnetic signal
The method displayed in Fig. 3 (Supplementary Figs. 6 and 7) was applied to each of the 66 superbolts observed in burst mode in 2012–2018 (cf. list and characteristics in Supplementary Tables 5 and 6), to analyze electric and magnetic field magnitudes in space. Ground squared electric field and PSD is similarly extracted from the 368 ECLAIR ground measurements of superbolts26. Mean-squared electric and magnetic fields are computed over a time period starting at the WWLLN superbolt time and lasting 1.5 ms on the ground and 1 s in space. These means are then plotted in Fig. 4, scaled by the WWLLN energy, and plotted against distance from the source in order to present a unique global view of superbolts’ electromagnetic mean power density on Earth (Fig. 4a) and in Space (Fig. 4b, c). We find that ground EM power decay follows an inverse power law with a power ~2 over distance, which confirms a well-known far-field decay for EM waves in free space. In space, at the equator, superbolt EM power decays with an inverse power law with a power ~1.6 over distance. In comparison, normal lightning VLF power decays with a power law varying between ~1.727 and ~2.323 when measured from space at low altitude (~700 km) (see also ref. 28). VLF power loss between ground and space is expected due to wave power spreading out globally in the Earth-ionosphere waveguide and attenuation of the VLF signal by the atmosphere/ionosphere during propagation. The variability of power in space among individual superbolts with similar energy and distance is also found to be large (1–2 orders of magnitude). This variability is consistent with the great variability in L-shell, longitude, local time, and season of all lightning-generated waves21. However, such a variability limits the utility of the power density laws of Fig. 4 and put forward the need of reliable full modeling29.
Ground-to-space transmission factor
Using these empirical fall-off relations with distance, we can rescale backward the space electric PSD and forward propagate the ground electric PSD observed in Fig. 3 to the same altitude (here 300 km where the law of Fig. 4b starts) and take their ratio to obtain a transmission factor from ground to space per frequency for the VLF range reported in Supplementary Table 7. The transmission factor is ~10−8 in the VLF range, which represents the lower bound of recent computations29, noticing the latter are widely spread over orders of magnitude according to simulation parameters. The transmission factor of the second event observed synchronously is also reported in Supplementary Table 7, smaller due to the far distance of the superbolt. These two unique and rare synchronous cases we report can serve for modeling the VLF wave emission by ray-tracing, starting from the source in the atmosphere and reaching the magnetic equator in space, which has important applications for radiation belt modeling30.
Ground statistics of superbolts
Statistics of all superbolts measured with the ground network are plotted in Fig. 5a–e, along with statistics obtained from all 3349 lightning strokes measured from ECLAIR ground stations in the same period (09/2012-06/2013) (in green) in order to highlight the unique characteristics of superbolts (in red and blue). Superbolt median peak current estimate31 is 363 kA, ~10 times higher than normal lightning (panel 5a). Fifteen percent of superbolts are positive cloud-to-ground (+CG) flashes, compared with 28% of +CG lightning triggered during the whole campaign. While a typical rise time (Fig. 5b and Supplementary Fig. 8c) of normal lightning is 1–2 µs (up to 5 µs), superbolts last longer, 5–6 µs (up to 15 µs). The decay time of the ground wave (from maximum to zero) is between 50 and 90 µs for normal negative cloud-to-ground (-CG) lightning, giving a well-known asymmetric shape to the ground wave5 for normal strokes. On the contrary, superbolts have a fast decay time of 10–20 µs (panel 5c), comparable to the rise time (Fig. 5b and Supplementary Fig. 8c), leading to a more symmetric ground wave (as seen in Fig. 3f and Supplementary Fig. 7e). This remarkable symmetry of the discharge, existing, for instance, for narrow bipolar pulses5, could help reveal physical differences between superbolts and other lightning. The median of the ground wave electric field (reported at ground distance of 100 km) is +240 V/m to be compared with + /−20 V/m for normal lightning (panel 5d). Median squared electric field, computed over 1.5 ms is 2.6 (V/m)2 on the ground, about 100 times higher than for normal strokes (panel 5e). We find on average 38% of the total power (2 kHz–5 MHz) is in the VLF range (2–12 kHz) for normal lightning while this percentage reaches 68% for superbolts (cf. distribution in Supplementary Fig. 9a). This percentage increases as the lightning power increases (Supplementary Fig. 9b). This suggests that is the longer duration of powerful lightning that populates the lowest frequencies of the VLF range.
Space statistics of superbolts
In space, another remarkable characteristic in these statistics is that superbolts are rare during day time, with 5/66 events occurring from ~9am to ~5 pm local time (Fig. 6a, pink histogram). Yet there is no sensitivity to day/night measured on the ground (6a, black/blue lines). This local time dependence is similar to the one found for non-superbolt lightning-generated waves (cf. second panel of Fig. 4 in ref. 21) and supports previous findings that the day time ionosphere makes the propagation of any lightning-generated waves into the magnetosphere more difficult. Superbolts detected in space from local burst-mode measurements are measured at a median satellite L-shell of 1.6 (~4000 km altitude at the equator). The mean [median] distance between WWLLN superbolt location and the closest magnetic footprint of the Van Allen Probes is 5924 km [5023 km]. They have a long-time delay (0.1–0.8 s) that corresponds to propagation time in the Earth-ionosphere waveguide and along field lines in the magnetosphere to the satellite, with some possible reflection between the conjugate points. The time delay, δt, follows therefore a linear scaling with respect to L-shell, δt = 0.277L – 0.225 in second, that is established in (panel 6b). Further ray-tracing studies should be able to answer whether this delay is due to ducted wave propagation along field lines in density depletions/enhancements (ducts) or to unducted wave propagation32, which is a general open problem of wave propagation25,33. Of the 10,724 superbolt events identified by WWLLN in 2012–2018, 431 WWLLN-detected superbolts were also captured by EMFISIS lower-resolution survey-mode data between 2012 and 2018. However, Van Allen Probes survey-mode measurements collect 0.5-s long records every 6 s. As superbolts signals in space typically exceed 1-s duration, the survey measurements usually do not contain the entire signal, thus truncating the total measured power. We show in Fig. 6(c, d) that survey-mode superbolt average power (black lines) is statistically larger than normal lightning, but not as powerful as burst-mode determined average power (pink histogram). This artificial bias is revealed by the burst-mode squared electric/magnetic field (pink, Fig. 6c/6d), which is found to be 10–1000 times more powerful than normal lightning (with global statistics21 scaled in green for comparison). The three strongest events have a ~250 pT (6 mV/m) wave magnetic (electric) rms amplitude in Fig. 6c, d. The mean magnetic (electric) rms amplitude of the 66 superbolts is 83 pT (873 µV/m) in Fig. 6c, d. In comparison, the mean magnetic (electric) rms amplitudes of lightning-generated waves computed from survey measurements of all types of lightning strokes (possibly including superbolts) is ~83 (~43) times lower, with 1 ± 1.6 pT (19 ± 59 μV/m)21.
The maximum of the burst-mode squared electric/magnetic field within the superbolt time window is provided in Supplementary Fig. 10. The comparison of Supplementary Fig. 10 with Fig. 6(c, d) shows the peak power exceeds by a factor ~10 (13) the mean magnetic (electric) power. Mean peak superbolt magnetic (electric) amplitude is 0.2 nT (1.9 mV2/m2) (Supplementary Fig. 10). Peak values are useful for nonlinear computations briefly mentioned below. Intercomparison of superbolt’s electromagnetic power with other natural waves such that whistler-mode hiss waves20,34,35,36,37,38,39,40,41,42, whistler-mode chorus waves20,43, and ElectroMagnetic Ion Cyclotron (EMIC) waves44 is left for future work.
Given spatial and temporal coverage of the Van Allen Probes data set, as well as the burst trigger algorithms, these events likely represent (to order-of-magnitude) the maximum power attainable for lightning-driven VLF waves near the Van Allen Probes. Such values may be useful for defining the dynamic range of future wave instruments or as an upper bound on lightning VLF wave power in simulations or numeric calculations. Extremely powerful lightning events including superbolts (>100 pT, representing 0.02% of lightning-generated waves) have been identified as strongly contributing to the global rms magnetic amplitude of lightning-generated waves for L < 2 (e.g., 44% at L = 1.1)21. These results suggest thus that the tail of the power density distribution of any type of waves (such as whistler-mode or EMIC waves) needs to be derived from, or at least confirmed by, from burst-mode measurements in order not to under-estimate the true wave power that is required for accurate radiation belt modeling (while today most statistical models are built from survey-mode data21,34,35,37,39,41,42). Finally, superbolt’s extreme power is a good candidate of highly powerful waves for which the quasilinear approximation (commonly used to compute wave-particle interactions in the radiation belts45,46,47,48,49,50) should break51,52, and to study further various nonlinear wave-particle interaction models51,52,53,54,55,56,57,58.
This study shows simultaneous measurements of a superbolt on the ground and the intense VLF signals in space driven by lightning. Ground and space observations demonstrate how different and extreme superbolts can be compared with normal lightning, for reasons not yet established. This study is a significant characterization of the previously unknown superbolts’ electromagnetic power, that should guide modeling and understanding of lightning electrodynamics, atmospheric discharges, and wave transmission from Earth to space, with applications in remote sensing, and wave modeling in space for radiation belt physics. Simultaneous optical and electromagnetic observations59,60 should be critical to help reveal more mysteries of superbolts.
Ground-based measurements of superbolts
The north Atlantic in western of Europe and the Mediterranean Sea in southern-western of Europe have some of the highest wintertime occurrence rates of superbolts7. This is a region covered by Commissariat à l’Energie Atomique et aux Energies Alternatives (CEA) ground stations14,15 in France. The ground stations measure electric fields from a few hundred Hertz to 5 MHz14 (VLF to High Frequency; HF) with vertical dipole whip antennas, at a sampling rate of 12.5 MHz.
The WWLLN network identifies 10,724 superbolts world-wide (gray dots, Supplementary Fig. 1) with a radiated energy greater than 1 MJ from synchronized measurements by >7 ground receivers7 (with a WWLLN residual time of group arrival for all stations of less than 30 μs) for the years 2012–2018, out of 1.5 × 109 total strokes detected (see Supplementary Method 2 for the influence of the WWLLN minimum station number detecting a stroke and the residual time on the total number of superbolts). Thus, the occurrence frequency of superbolts in this WWLLN data set is seven events per million CG lightning events, yielding an average of 3.5 superbolt events per day world-wide (with a four times higher probability of occurrence in northern winter), and implying a very low probability to observe them simultaneously from space and on the ground. Among these 10,724 superbolts, 4034 were identified in Europe (Fig. 1), of which 384 occurred between 09/2012 and 06/2013, corresponding to a ground measurement campaign (ECLAIR) conducted by CEA. Of these, 368 high-resolution superbolt ground waveforms were recorded during ECLAIR26, with 86 of them also detected by the stations of the French lightning location network, Météorage (MTRG)16, within 1500 km of the ECLAIR network center (Fig. 1). These ground measurements (including WWLLN) provide superbolt occurrence time (with time accuracy <1 ms), geographic location, and WWLLN energy. They also provide ECLAIR high-resolution waveforms with rise and decay times, and ECLAIR and MTRG peak current estimates31 (see Supplementary Method 1 for more information).
Space measurements of superbolts
In space, we use the Electric and Magnetic Field Instrument Suite and Integrated Science (EMFISIS) instrument suite17 and the Electric Fields and Waves (EFW) instrument18 on board the two probes (Probe A and Probe B) of the NASA Van Allen Probes spacecraft19 to measure electric and magnetic wave fields near the magnetic equatorial plane (+/−20° magnetic latitude). The instruments provide similar burst data, but capture records at different times. These ground and space-based measurements are used to study the symmetry of the first ground-wave peak, determine the time it takes for the signal to propagate to the satellite, quantify the loss of signal strength with distance, compute the frequency-dependent transmission of the wave power, and generate statistics (see Supplementary Method 1 for more information).
Superbolt detection and selection
Of the 10,724 WWLLN superbolts, 1143 occurred during a time when Van Allen Probes EFW or EMFISIS were capturing burst data that records electromagnetic waves at high temporal resolution, called the waveform, although most of these waveforms did not measure a superbolt due to the satellite’s location. Only 212 of these burst captures occurred when the Van Allen Probes were at L-shell, L, less than three where most of lightning-generated waves are observed21. For facilitating further selection, composite figures (such as Figs. 2, 3 and Supplementary Figs. 6 and 7 discussed in the article) were generated for each of the 212 superbolts and used to screen them. Among them, 81 were retained because the on board burst system was triggered by identified whistler-mode waves and no unidentified (or other) radio wave power corrupted (partially or entirely) the burst signal. Additional selection criteria were then applied for retaining only clearly identified superbolts: (1) the lightning VLF signal was fully captured by the burst data (uninterrupted); (2) no other lightning flashes mask the superbolt signal (i.e., no lightning of higher or comparable estimated signal power22,23 at the satellite occurred within the burst window); (3) and the distance (direct line), d, from the superbolt location to the magnetic footpoint of the field line on which the satellite sits is <8000 km. Note that although we are requiring that the satellite L is less than three, there is no restriction on the superbolt location. For superbolts with d > 8000 km, we identified and retained some additional weak signals by hand that satisfied criteria (1) to (2). This process ensures the least possible contamination of the superbolt wave signal by other lightning. At the end of the selection process, only 66 burst captures of superbolt remained (38 from EFW and 28 from EMFISIS, with 1 recorded by both) out of the 1143 candidate events during which burst data was also collected. Among the 66 selected burst-mode signals, 33 were located by WWLLN to be over Europe (cf. Fig. 1), and only two of these were simultaneously recorded by ECLAIR, with one too far (9435 km away) from the closest Van Allen Probes magnetic footprint to be clearly identified. For these two events there is high-resolution data available from both the ground and space instrumentation that permits a detailed study linking the superbolt properties on the ground and in space, discussed below. Datasets are summarized in Supplementary Table 1.
Superbolt electromagnetic power
High-resolution power spectral density (PSD) from a selection of six superbolts is shown in Fig. 2, computed from either the burst-mode electric field (a1–d1) or the ground-based (e1–f1) electric field waveforms. These waveforms (displayed in Supplementary Figs. 2, 3, 4) were decomposed by Fourier transform to produce electric field PSDs. The respective Van Allen Probes magnetic field PSDs are shown in Supplementary Fig. 5. Characteristics of all superbolts discussed specifically in the article are gathered in Supplementary Table 2 and 3. The black line in Fig. 2 (a2–f2) (see also the respective magnetic field intensity in Supplementary Fig. 5) represents the evolution of the square of the electric (magnetic) field wave, E2 in V2/m2 (B2 in pT2), a.k.a. the wave electric (magnetic) intensity or squared amplitude or, more commonly called, the wave power, that is the frequency integral of the electric (magnetic) field PSD over the VLF range in space (~2–10 kHz) and (~2 kHz–5 MHz) on the ground.
ECLAIR ground waveforms26 are available at https://zenodo.org/record/3952425. Météorage data may be ordered from https://www.météorage.com. Van Allen Probes field data are available from the EFW and EMFISIS team websites, which one can link to here: http://rbspgway.jhuapl.edu. The WWLLN data are available from any WWLLN host or may be ordered from links (http://wwlln.net). The satellite Situation Center Locator operated online by NASA provides Van Allen Probes trajectories (https://sscweb.gsfc.nasa.gov/).
Turman, B. N. Detection of lightning superbolts. J. Geophys. Res. 82, 2566–2568 (1977).
Turman, B. N. Analysis of lightning data from the DMSP satellite. J. Geophys. Res. 83, 5019–5024 (1978).
Ripoll, J.-F., Zinn, J., Jeffery, C. A. & Colestock, P. L. On the dynamics of hot air plasmas related to lightning discharges: 1. Gas dynamics. J. Geophys. Res. Atmos. 119, https://doi.org/10.1002/2013JD020067 (2014a).
Ripoll, J.-F., Zinn, J., Colestock, P. L. & Jeffery, C. A. On the dynamics of hot air plasmas related to lightning discharges: 2. Electrodynamics. J. Geophys. Res.: Atmos. 119, 9218–9235 (2014).
Rakov, V. A. & Uman, M. A. Lightning: Physics and Effects (Cambridge Univ. Press, Cambridge, 2007).
Borovsky, J. E. Lightning energetics: Estimates of energy dissipation in channels, channel radii, and channel-heating risetimes. J. Geophys. Res. 103, 11,537–11,553 (1998).
Holzworth, R. H., McCarthy, M. P., Brundell, J. B., Jacobson, A. R. & Rodger, C. J. Global distribution of superbolts. J. Geophys. Res.: Atmos. 124, 9996–10,005 (2019).
Hutchins, M. L., Holzworth, R. H., Rodger, C. J. & Brundell, J. B. Far-field power of lightning strokes as measured by the World Wide Lightning Location Network. J. Atmos. Ocean. Technol. 29, 1102–1110 (2012b).
Dowden, R. L., Brundell, J. B. & Rodger, C. J. VLF lightning location by time of group arrival (TOGA) at multiple sites. J. Atmos. Sol.-Terrestrial Phys. 64, 817–830 (2002).
Lay, E. H. et al. WWLL global lightning detection system: regional validation study in Brazil. Geophys. Res. Lett. 31, L03102 (2004).
Rodger C. J., Brundell, J. B., Holzworth, H. & Lay, E. H. Growing Detection Efficiency of the World Wide Lightning Location Network, CP1118, Coupling of Thunderstorms and Lightning Discharges to Near-Earth (eds Crosby, N. B., Huang, T. Y. & Rycroft, M. J.) (American Institute of Physics, 2009).
Holzworth, R. H. et al. Lightning-generated whistler waves observed by probes on the Communication/Navigation Outage Forecast System satellite at low latitudes. J. Geophys. Res. 116, A06306 (2011).
Helliwell, R. A. Whistlers and Related Ionospheric Phenomena (Stanford University Press, 1965).
Farges, T. & Blanc, E. Lightning and TLE electric fields and their impact on the ionosphere. C. R. Phys. 12, 171–179 (2011).
Kolmašová, I. et al. Subionospheric propagation and peak currents of preliminary breakdown pulses before negative cloud-to-ground lightning discharges. Geophys. Res. Lett. 43, 1382–1391 (2016).
Pédeboy, S., Defer, E., & Schulz, W. Performance of the EUCLID network in cloud lightning detection in the South-East France. 8th HyMeX Workshop, Valletta, Malta (2014).
Kletzing, C. A. et al. The electric and magnetic field instrument suite and integrated science (EMFISIS) on van allen probes. Space Sci. Rev. 179, 127–181 (2013).
Wygant, J. R. et al. The electric field and waves instruments on the radiation belt storm probes mission. Space Sci. Rev. 179, 183–220 (2013).
Mauk, B. H. et al. Science objectives and rationale for the radiation belt storm probes mission. Space Sci. Rev. 179, 3–27 (2013).
Hospodarsky, G. B. et al. (2016). in Magnetosphere-ionosphere coupling in the solar system (eds Chappell, C. R. et al.) 127–143 (John Wiley, 2016).
Ripoll, J.-F. et al. Analysis of electric and magnetic lightning-generated wave amplitudes measured by the Van Allen Probes. Geophys. Res. Lett. 47, e2020GL087503 (2020).
Burkholder, B. S., Hutchins, M. L., McCarthy, M. P., Pfaff, R. F. & Holzworth, R. H. Attenuation of lightning-produced sferics in the Earth-ionosphere waveguide and low-latitude ionosphere. J. Geophys. Res. Space Phys. 118, 3692–3699 (2013).
Ripoll, J.-F., Farges, T., Lay, E. H., & Cunningham, G. S. (2019). Local and statistical maps of lightning-generated wave power density estimated at the Van Allen Probes footprints from the World-Wide Lightning Location Network database. Geophys. Res. Lett. 46, https://doi.org/10.1029/2018GL081146.
Bortnik, J. Precipitation of Radiation Belt Electrons by Lightning-generated Magnetospherically Reflecting Whistlers. Ph.D. thesis, Stanford University (2004). https://vlf.stanford.edu/pubs/theses
Ripoll, J.-F. et al. Particle dynamics in the Earth’s radiation belts: review of current research and open questions. J. Geophys. Res. 125, e2019JA026735 (2020).
Farges, T. ECLAIR Superbolt Waveforms, https://doi.org/10.5281/zenodo.3952425 (2021).
Fiser, J. C., Diendorfer, G., Parrot, M. & Santolik, O. Whistler intensities above thunderstorms. Ann. Geophys. 28, 37–46 (2010).
Jacobson, A. R., Holzworth, R. H., Pfaff, R. F. & McCarthy, M. P. Study of oblique whistlers in the low-latitude ionosphere, jointly with the C/NOFS satellite and the World-Wide Lightning Location Network. Ann. Geophys. 29, 851–863 (2011).
Jacobson, A. R., Holzworth, R. H., Pfaff, R. & Heelis, R. Coordinated satellite observations of the very low frequency transmission through the ionospheric D layer at low latitudes, using broadband radio emissions from lightning. J. Geophys. Res. 123. https://doi.org/10.1002/2017JA024942 (2018)
Starks, M. J., Albert, J. M., Ling, A., O’Malley, S. & Quinn, R. A. VLF transmitters and lightning-generated whistlers: 1. Modeling waves from source to space. J. Geophys. Res. 125, e2019JA027029 (2020).
Météorage. Calculation principle of peak intensity of stroke current. Technical Note, https://www.meteorage.com/sites/default/files/2018-12/TN%20-%20Calculation%20principle%20of%20peak%20intensity%20of%20stroke%20current_EN.pdf. (2015)
Sonwalkar, V. S. In Geospace Electromagnetic Waves and Radiation (eds LaBelle, J. W. & Treumann, R. A.) Ch. 6, p. 141 (Springer, 2006) https://doi.org/10.1007/3-540-33203-0_6
Clilverd, M. A. et al. Ground-based transmitter signals observed from space: ducted or nonducted? J. Geophys. Res. 113, A04211 (2008).
Meredith, N. P. et al. Global model of plasmaspheric hiss from multiple satellite observations. J. Geophys. Res. 123, 4526–4541 (2018).
Malaspina, D. M. et al. The distribution of plasmaspheric hiss wave power with respect to plasmapause location. Geophys. Res. Lett. 43, 7878–7886 (2016).
Tsurutani, B. T., Falkowski, B. J., Pickett, J. S., Santolik, O. & Lakhina, G. S. Plasmaspheric hiss properties: observations from Polar. J. Geophys. Res. 120, 414–431 (2015).
Hartley, D. P., Kletzing, C. A., Santolík, O., Chen, L. & Horne, R. B. Statistical properties of plasmaspheric hiss from Van Allen Probes observations. J. Geophys. Res. 123, 2605–2619 (2018).
Gurnett, D. A. et al. The Polar plasma wave instrument. Space Sci. Rev. 71, 597–622 (1995).
Falkowski, B. J., Tsurutani, B. T., Lakhina, G. S. & Pickett, J. S. Two sources of dayside intense, quasi-coherent plasmaspheric hiss: a new mechanism for the slot region? J. Geophys. Res. 122, 1643–1657 (2017).
Tsurutani, B. T. et al. Plasmaspheric hiss: coherent and intense. J. Geophys. Res. 123, 10,009–10,029 (2018).
Malaspina, D. M., Ripoll, J.-F., Chu, X., Hospodarsky, G. & Wygant, J. Variation in plasmaspheric hiss wave power with plasma density. Geophys. Res. Lett. 45, https://doi.org/10.1029/2018GL078564 (2018)
Shi, R. et al. Properties of whistler mode waves in Earth’s plasmasphere and plumes. J. Geophys. Res. 124, https://doi.org/10.1029/2018JA026041 (2019)
Cattell, C. et al. Discovery of very large amplitude whistler- mode waves in Earth’s radiation belts. Geophys. Res. Lett. 35, L01105 (2008).
Engebretson, M. J. et al. Van Allen probes, NOAA, GOES, and ground observations of an intense EMIC wave event extending over 12 h in magnetic local time. J. Geophys. Res. 120, 5465–5488 (2015).
Abel, B. & Thorne, R. M. Electron scattering and loss in Earth’s inner magnetosphere, 1: Dominant physical processes. J. Geophys. Res. 103, 2385–2396 (1998).
Kim, K.-C., Shprits, Y., Subbotin, D. & Ni, B. Understanding the dynamic evolution of the relativistic electron slot region including radial and pitch angle diffusion. J. Geophys. Res. 116, A10214 (2011).
Ripoll, J.-F., Chen, Y., Fennell, J. F. & Friedel, R. H. W. On long decays of electrons in the vicinity of the slot region observed by HEO3. J. Geophys. Res. 120, 460–478 (2014).
Ripoll, J.-F. et al. Effects of whistler mode hiss waves in March 2013. J. Geophys. Res. 122, 7433–7462 (2017).
Ripoll, J.‐F. et al. Observations and Fokker-Planck simulations of the L-shell, energy, and pitch angle structure of Earth’s electron radiation belts during quiet times. J. Geophys. Res. 124, 1125–1142 (2019).
Zhao, H. et al. Plasmaspheric hiss waves generate a reversed energy spectrum of radiation belt electrons. Nat. Phys. https://doi.org/10.1038/s41567-018-0391-6 (2019).
Albert, J. M. & Bortnik, J. Nonlinear interaction of radiation belt electrons with electromagnetic ion cyclotron waves. Geophys. Res. Lett. 36, L12110 (2009).
Tao, X., Bortnik, J., Thorne, R. M., Albert, J. M. & Li, W. Effects of amplitude modulation on nonlinear interactions between elec- trons and chorus waves. Geophys. Res. Lett. 39, L06102 (2012).
Nunn, D. & Omura, Y. A computational and theoretical investigation of nonlinear wave-particle interactions in oblique whistlers. J. Geophys. Res. 120, 2890–2911 (2015).
Hsieh, Y. K., Kubota, Y. & Omura, Y. Nonlinear evolution of radiation belt electron fluxes interacting with oblique whistler mode chorus emissions. J. Geophys. Res. 125, e2019JA027465 (2020).
da Silva, C. L. et al. Test-particle simulations of linear and nonlinear interactions between a 2-D whistler-mode wave packet and radiation belt electrons. Geophys. Res. Lett. 45, 5234–5245 (2018).
Denton, R. E. et al. Pitch angle scattering of sub-MeV relativistic electrons by electromagnetic ion cyclotron waves. J. Geophys. Res. 124, 5610–5626 (2019).
Tsurutani, B. T. et al. Low frequency (f<200 Hz) polar plasmaspherichiss: coherent and intense. J. Geophys. Res. 124,10063–10084 (2019).
Delzanno, G. L. & Roytershteyn, V. High-frequency plasma waves and pitch angle scattering induced by pulsed electron beams. J. Geophys. Res. 124, https://doi.org/10.1029/2019JA027046 (2019).
Lefeuvre, F. et al. TARANIS—a satellite project dedicated tothe physics of TLEs and TGFs. Space Sci. Rev. 137, 301–315 (2008).
Farges, T., Hébert, P., Le Mer-Dachard, F., Ravel, K. & Gaillac, S. MicroCameras and photometers (MCP) on board the TARANIS satellite, XVI International Conference on Atmospheric Electricity, 17–22 June, Nara, Japan (2018).
The authors thank the World Wide Lightning Location Network, a collaboration among over 50 universities and institutions, for providing the lightning location and energy data used in this paper. The work of J.-F.R., T.F., G.S.C., and E.H.L. was performed under the auspices of an agreement between CEA/DAM (Commissariat à l’Energie Atomique, Direction des Applications Militaires) and NNSA/DP (National Nuclear Security Administration, Defense Program) on cooperation on fundamental science. The authors acknowledge the International Space Sciences Institute (ISSI). The work of GSC was supported in part by the Defense Threat Reduction Agency (DTRA). Thanks also to DTRA for supporting a collaborative visit by JFR at Los Alamos National Laboratory in 2018. The work of DM was supported by the National Science Foundation under award 1841011. The authors thank the entire Van Allen Probes team, and especially the EFW and EMFISIS teams for their support.
The authors declare no competing interests.
Peer review information Nature Communications thanks Robert Holzworth and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Ripoll, JF., Farges, T., Malaspina, D.M. et al. Electromagnetic power of lightning superbolts from Earth to space. Nat Commun 12, 3553 (2021). https://doi.org/10.1038/s41467-021-23740-6
This article is cited by
Advanced Composites and Hybrid Materials (2022)
The Van Allen Probes Electric Field and Waves Instrument: Science Results, Measurements, and Access to Data
Space Science Reviews (2022)