Discovery of Intense Gamma-Ray Flashes of Atmospheric Origin

Employing optical lightning data to identify lightning flashes associated to terrestrial gamma-ray flashes

Terrestrial gamma-ray flashes (TGFs) are bursts of energetic X- and gamma-rays emitted from thunderstorms. The Atmosphere-Space Interactions Monitor (ASIM) mounted onto the International Space Station (ISS) is dedicated to measure TGFs and optical signatures of lightning; ISS LIS (Lightning Imaging Sensor) detects lightning flashes allowing for simultaneous measurements with ASIM. Whilst ASIM measures ∼ 300-400 TGFs per year, ISS LIS detects ∼ 10 6 annual lightning flashes giving a disproportion of four orders of magnitude. Based on the temporal evolution of lightning flashes and the spatial pattern of the constituing events, we present an algorithm to reduce the number of space-detected flashes potentially associated with TGFs by ∼ 60% and of associated LIS groups by ∼ 95%. We use ASIM measurements to confirm that the resulting flashes are indeed those associated with TGFs detected at approx. 400 km altitude and thus benchmark our algorithm preserving 70-80% of TGFs from concurrent ASIM-LIS measurements. We compare how the radiance, footprint size and the global distribution of lightning flashes of the reduced set relates to the average of all measured lightning flashes and do not observe any significant difference. Finally, we present a parameter study of our algorithm and discuss which parameters can be tweaked to maximize the reduction efficiency whilst keeping flashes associated to TGFs. In the future, this algorithm will hence be capable of facilitating the search for TGFs in a reduced set of lightning flashes measured from space.

Imaging of 3 bright terrestrial gamma-ray flashes by the atmosphere-space interactions monitor and their parent thunderstorms

The Atmosphere-Space Interactions Monitor (ASIM) on the International Space Station (ISS) includes an instrument designed to geolocate Terrestrial Gamma-ray Flashes (TGF) produced by thunderstorms. It does so with a coded aperture system shadowing the pixelated Low Energy Detector of the Modular X- and Gamma-ray Sensor (MXGS). Additionally, it locates associated lightning flashes with the Modular Multispectral Imaging Array (MMIA). Here we present 3 bright TGFs with very similar duration, fluency and observation distance. The innovative imaging capabilities allow us to determine the TGF position and correlate the TGF-lightning parent event in time and position simultaneously. The accurate position determination and distance to the observer allow us to perform a comparative study of TGF characteristics. These TGFs were produced in association with lightning flashes below the highest cloud tops of developing to mature convective cells. In one event, GLM (Geostationary Lightning Mapper) cloud flash rates were slowing down after the TGF while negative cloud-to-ground flashes suddenly ceased from 10 to 5 min before the TGF. An 8-stroke (strongest: -93 kA) cloud-to-ground flash occurred 31 s before the TGF. Vertical profiles from the ERA5 reanalysis data suggest TGFs may be produced in variety of tropical environments.

Optical observations of thunderstorms from the International Space Station: recent results and perspectives

The International Space Station (ISS) is in the lowest available orbit at ~400 km altitude, bringing instruments as close to the atmosphere as possible from the vantage point of space. The orbit inclination is 51.6°, which brings the ISS over all the low- and mid-latitude regions of the Earth and at all local times. It is an ideal platform to observe deep convection and electrification of thunderstorms, taken advantage of by the Lightning Imaging Sensor (LIS) and the Atmosphere Space Interaction Monitor (ASIM) experiments. In the coming years, meteorological satellites in geostationary orbit (~36,000 km altitude) will provide sophisticated cloud and lightning observations with almost complete coverage of the Earth's thunderstorm regions. In addition, Earth-observing satellite instruments in geostationary- and low-Earth orbit (LEO) will measure more atmospheric parameters at a higher resolution than we know today. The new infrastructure in space offers an opportunity to advance our understanding of the role of thunderstorms in atmospheric dynamics and climate change. Here, we discuss how observations from the ISS or other LEO platforms with instruments that view the atmosphere at slanted angles can complement the measurements from primarily nadir-oriented instruments of present and planned missions. We suggest that the slanted viewing geometry from LEO may resolve the altitude of electrical activity and the cloud structure where they occur, with implications for modelling thunderstorms' effects on the atmosphere's radiative properties and climate balance.

The Temporal Relationship Between Terrestrial Gamma-Ray Flashes and Associated Optical Pulses From Lightning

We present 221 Terrestrial Gamma-ray Flashes (TGFs) and associated optical pulses observed by the Atmosphere-Space Interactions Monitor (ASIM) on board the International Space Station. The events were detected between the end of March 2019 and November 2020 and consist of X- and gamma-ray energy detections, as well as photometer data (180-230, 337, and 777 nm) and optical camera data (337 and 777 nm). Using the available ASIM data and applying a consistency check based on TGF characteristics and lightning detections from lightning radio atmospherics close in time, we determine the most likely position of the TGFs in relation to the photometer field of view (FoV), and the association to the observed optical pulses. Out of the 221 events we find 72 events where the TGF and optical data are determined to be associated and inside the photometer FoV. Using the measured TGF durations and the time between the onsets of the TGFs and optical pulses we find: (a) That the TGF onsets are always before or at the same time as the optical pulse onsets (taking into account cloud scattering). (b) A tendency for longer duration TGFs to have longer delays between onsets. (c) Two groups of events: (a) where there is a possible overlap between the TGFs and the optical emissions, as the TGFs last longer than the delay between onsets and (b) where the TGFs and optical emissions do not overlap, as there are long delays between the onsets, which cannot be explained by cloud scattering.

Energetic Intracloud Lightning in the RELAMPAGO Field Campaign

A particular strength of lightning remote sensing is the variety of lightning types observed, each with a unique occurrence context and characteristically different emission. Distinct energetic intracloud (EIC) lightning discharges-compact intracloud lightning discharges (CIDs) and energetic intracloud pulses (EIPs)-produce intense RF radiation, suggesting large currents inside the cloud, and they also have different production mechanisms and occurrence contexts. A Low-Frequency (LF) lightning remote sensing instrument array was deployed during the RELAMPAGO field campaign in west central Argentina, designed to investigate convective storms that produce high-impact weather. LF data from the campaign can provide a valuable data set for researching the lightning context of EICs in a variety of subtropical convective storms. This paper describes the production of an LF-CID data set in RELAMPAGO and includes a preliminary analysis of CID prevalence. Geolocated lightning events and their corresponding observed waveforms from the RELAMPAGO LF data set are used in the classification of EICs. Height estimates based on skywave reflections are computed, where prefit residual data editing is used to improve robustness against outliers. Even if EIPs occurred within the network, given the low number of very high-peak current events and receiver saturation, automatic classification of EIPs may not be feasible using this data. The classification of CIDs, on the other hand, is straightforward and their properties, for both positive and negative polarity, are investigated. A few RELAMPAGO case studies are also presented, where high variability of CID prevalence in ordinary storms and high-altitude positive CIDs, possibly in overshooting tops, are observed.

Synchrotron mechanism of X-ray and gamma-ray emissions in lightning and spark discharges

X-ray and γ-ray emissions observed in lightning and long sparks are usually connected with the bremsstrahlung of high-energy runaway electrons. Here, an alternative physical mechanism for producing X-ray and gamma-ray emissions caused by the polarization current and associated electromagnetic field moving with relativistic velocity along a curved discharge channel has been proposed. The existence of fast electromagnetic surface waves propagating along the lightning discharge channel at a speed close to the speed of light in vacuum is shown. The possibility of the production of microwave, X-ray and gamma-ray emissions by a polarization current pulse moving along a curved path via synchrotron radiation mechanism is pointed out. The existence of long tails in the power spectrum is shown, which explains observations of photon energies in the range of 10–100 MeV in the terrestrial gamma-ray flashes, as well as measured power spectrum of laboratory spark discharge.

Ionization Waves Enhance the Production of X-rays during Streamer Collisions

Experimental data show that in laboratory sparks, X-rays are produced in time synchronization with the meeting of streamers of opposite polarity just before the final breakdown of the discharge gap. It has been suggested that the electric field enhancement created during the collision of streamers could provide the necessary conditions for electron acceleration, even though some of the theoretical studies show that the duration of the electric field is not long enough to do so. The experimental data on laboratory discharges show that. when streamers of opposite polarity meet each other, a potential or ionization wave that renders the streamer channels conducting is initiated. This paper shows that these ionization waves that convert the discharge channels from weakly conducting to highly conducting are associated with electric fields large enough to accelerate electrons to relativistic energies.

Thunderstorms Producing Sferic-Geolocated Gamma-Ray Flashes Detected by TETRA-II

The terrestrial gamma-ray flash (TGF) and Energetic Thunderstorm Rooftop Array (TETRA-II) detected 22 X-ray/gamma-ray flash events associated with lightning between October 2015 and March 2019 across three ground-based detector locations in subtropical and tropical climates in Louisiana, Puerto Rico, and Panama. Each detector array consists of a set of bismuth germanate scintillators that record X-ray and gamma-ray bursts over the energy range 50 keV-6 MeV (million electron volts). TETRA-II events have characteristics similar to both X-ray bursts associated with lightning leaders and TGFs: sub-millisecond duration, photons up to MeV energies, and association with nearby lightning (typically within 3 km). About 20 of the 22 events are geolocated to individual lightning strokes via spatiotemporally coincident sferics. An examination of radar reflectivity and derived products related to events located within the Next Generation Weather Radar (NEXRAD) monitoring region indicates that events occur within mature cells of severe and non-severe multicellular and squall line thunderstorms, with core echo tops which are at or nearing peak altitude. Events occur in both high lightning frequency thunderstorm cells and low lightning frequency cells. Events associated with high frequency cells occur within 5 min of significant lightning jumps. Among NEXRAD-monitored events, hail is present within 8 km and 5 min of all except a single low-altitude cold weather thunderstorm. An association is seen with maximum thunderstorm development, lightning jumps, and hail cells, indicating that the TETRA-II X-ray/gamma-ray events are associated with the peak storm electrification and development of electric fields necessary for the acceleration of electrons to high energies.

Hazards to Aircraft Crews, Passengers, and Equipment from Thunderstorm-Generated X-rays and Gamma-Rays

Both observational and theoretical research in the area of atmospheric high-energy physics since about 1980 has revealed that thunderstorms produce X-rays and gamma-rays into the MeV region by a number of mechanisms. While the nature of these mechanisms is still an area of active research, enough observational and theoretical data exists to permit an evaluation of hazards presented by ionizing radiation from thunderstorms to aircraft crew, passengers, and equipment. In this paper, we use data from existing studies to evaluate these hazards in a quantitative way. We find that hazards to humans are generally low, although with the possibility of an isolated rare incident giving rise to enough radiation dose to produce noticeable symptoms. On the other hand, unshielded computer memory chips in avionics systems stand a small but non-zero chance of severe damage from thunderstorm-generated radiation and would not leave easily detectable traces of the occurrence. Should a rare phenomenon called ball lightning occur near or within an aircraft, the possibility exists of substantial damage to both equipment and personnel. Overall, radiation hazards from thunderstorms appear to be low, but should be considered and investigated with radiation monitoring equipment on sample flights.

Implications of GNSS-Inferred Tropopause Altitude Associated with Terrestrial Gamma-ray Flashes

The thermal structure of the environmental atmosphere associated with Terrestrial Gamma-ray Flashes (TGFs) is investigated with the combined observations from several detectors (FERMI, RHESSI, and Insight-HXMT) and GNSS-RO (SAC-C, COSMIC, GRACE, TerraSAR-X, and MetOp-A). The geographic distributions of TGF-related tropopause altitude and climatology are similar. The regional TGF-related tropopause altitude in Africa and the Caribbean Sea is 0.1–0.4 km lower than the climatology, whereas that in Asia is 0.1–0.2 km higher. Most of the TGF-related tropopause altitudes are slightly higher than the climatology, while some of them have a slightly negative bias. The subtropical TGF-producing thunderstorms are warmer in the troposphere and have a colder and higher tropopause over land than the ocean. There is no significant land–ocean difference in the thermal structure for the tropical TGF-producing thunderstorms. The TGF-producing thunderstorms have a cold anomaly in the middle and upper troposphere and have stronger anomalies than the deep convection found in previous studies.

Spectral Observations of Optical Emissions Associated With Terrestrial Gamma-Ray Flashes

The Atmosphere-Space Interactions Monitor measures Terrestrial Gamma-Ray Flashes (TGFs) simultaneously with optical emissions from associated lightning activity. We analyzed optical measurements at 180-230, 337, and 777.4 nm related to 69 TGFs observed between June 2018 and October 2019. All TGFs are associated with optical emissions and 90% of them are at the onset of a large optical pulse, suggesting that they are connected with the initiation of current surges. A model of photon delay induced by cloud scattering suggests that the sources of the optical pulses are from 0.7 ms before to 4.4 ms after the TGFs, with a median of -10 ± 80 µs, and 1-5 km below the cloud top. The pulses have rise times comparable to lightning but longer durations. Pulse amplitudes at 337 nm are ∼3 times larger than at 777.4 nm. The results support the leader-streamer mechanism for TGF generation.

Time Evolution of Storms Producing Terrestrial Gamma-Ray Flashes Using ERA5 Reanalysis Data, GPS, Lightning and Geostationary Satellite Observations

In this article, we report the first investigation over time of the atmospheric conditions around terrestrial gamma-ray flash (TGF) occurrences, using GPS sensors in combination with geostationary satellite observations and ERA5 reanalysis data. The goal is to understand which characteristics are favorable to the development of these events and to investigate if any precursor signals can be expected. A total of 9 TGFs, occurring at a distance lower than 45 km from a GPS sensor, were analyzed and two of them are shown here as an example analysis. Moreover, the lightning activity, collected by the World Wide Lightning Location Network (WWLLN), was used in order to identify any links and correlations with TGF occurrence and precipitable water vapor (PWV) trends. The combined use of GPS and the stroke rate trends identified, for all cases, a recurring pattern in which an increase in PWV is observed on a timescale of about two hours before the TGF occurrence that can be placed within the lightning peak. The temporal relation between the PWV trend and TGF occurrence is strictly related to the position of GPS sensors in relation to TGF coordinates. The life cycle of these storms observed by geostationary sensors described TGF-producing clouds as intense with a wide range of extensions and, in all cases, the TGF is located at the edge of the convective cell. Furthermore, the satellite data provide an added value in associating the GPS water vapor trend to the convective cell generating the TGF. The investigation with ERA5 reanalysis data showed that TGFs mainly occur in convective environments with unexceptional values with respect to the monthly average value of parameters measured at the same location. Moreover, the analysis showed the strong potential of the use of GPS data for the troposphere characterization in areas with complex territorial morphologies. This study provides indications on the dynamics of con-vective systems linked to TGFs and will certainly help refine our understanding of their production, as well as highlighting a potential approach through the use of GPS data to explore the lightning activity trend and TGF occurrences.

Glossary on atmospheric electricity and its effects on biology

There is an increasing interest to study the interactions between atmospheric electrical parameters and living organisms at multiple scales. So far, relatively few studies have been published that focus on possible biological effects of atmospheric electric and magnetic fields. To foster future work in this area of multidisciplinary research, here we present a glossary of relevant terms. Its main purpose is to facilitate the process of learning and communication among the different scientific disciplines working on this topic. While some definitions come from existing sources, other concepts have been re-defined to better reflect the existing and emerging scientific needs of this multidisciplinary and transdisciplinary area of research.

Radiation in the Atmosphere—A Hazard to Aviation Safety?

Exposure of aircrew to cosmic radiation has been recognized as an occupational health risk for several decades. Based on the recommendations by the International Commission on Radiological Protection (ICRP), many countries and their aviation authorities, respectively have either stipulated legal radiation protection regulations, e.g., in the European Union or issued corresponding advisory circulars, e.g., in the United States of America. Additional sources of ionizing and non-ionizing radiation, e.g., due to weather phenomena have been identified and discussed in the scientific literature in recent years. This article gives an overview of the different generally recognized sources due to weather as well as space weather phenomena that contribute to radiation exposure in the atmosphere and the associated radiation effects that might pose a risk to aviation safety at large, including effects on human health and avionics. Furthermore, potential mitigation measures for several radiation sources and the prerequisites for their use are discussed.

Electromagnetic Radiation Spectrum of a Composite System

In many physics fields, the radio emission of a composite system composed of a large number of randomly occurring but similar emission sources is encountered. In general, the composite system lasts longer than each individual component and individual source currents vary much more rapidly. This Letter presents a theory to understand the electromagnetic radiation spectrum of such a system. If the temporal distribution of the random occurrence of the component and the distribution to describe the relevant emission source properties are known, the spectrum of the composite system can be readily found from this theory. There are two main terms that define the spectrum: one term results from the coherent summation of the contributions from individual sources and is proportional to the square of the total number of the components in the system; the other term results from an incoherent summation and is proportional to the first power of that number. This can lead to drastically different spectral magnitudes in different spectral regions, typically with the spectral magnitude in the lower frequency region many orders of magnitude stronger than that in the higher frequency region.

Influence of Elastic Scattering on Electron Swarm Distribution in Electrified Gases

The propagation of energetic electrons through air is one key component in the generation of high-energy atmospheric phenomena such as lightning-generated X-ray bursts, terrestrial gamma ray flashes (TGFs), and gamma ray glows. We show here that models for this propagation can be considerably affected by the parameterization of the differential cross section of elastic scattering of electrons on the molecular components of air. We assess existing parameterizations and propose a more accurate one that builds upon the most up-to-date measurements. Then we conclude that by overweighting the forward scattering probability, previous works may have overestimated the production of runaway electrons under high electric fields close to the thermal runaway threshold.

Gamma-Ray and Radio-Frequency Radiation from Thunderstorms Observed from Space and Ground

Terrestrial gamma ray flashes (TGFs) are a class of enigmatic electrical discharges in the Earth’s atmosphere. In this study, we analyze an unprecedentedly large dataset comprised of 2188 TGFs whose signatures were simultaneously measured using space- and ground-based detectors over a five-year period. The Gamma-ray Burst Monitor (GBM) on board the Fermi spacecraft provided the energetic radiation measurements. Radio frequency (RF) measurements were obtained from the Global Lightning Dataset (GLD360). Here we show the existence of two categories of TGFs − those that were accompanied by quasi-simultaneous electromagnetic pulses (EMPs) detected by the GLD360 and those without such simultaneous EMPs. We examined, for the first time, the dependence of the TGF-associated EMP-peak-amplitude on the horizontal offset distance between the Fermi spacecraft and the TGF source. TGFs detected by the GBM with sources at farther horizontal distances are expected to be intrinsically brighter and were found to be associated with EMPs having larger median peak-amplitudes. This provides independent evidence that the EMPs and TGFs are produced by the same phenomenon, rather than the EMPs being from “regular” lightning in TGF-producing thunderstorms.

A terrestrial gamma-ray flash and ionospheric ultraviolet emissions powered by lightning

Terrestrial gamma-ray flashes (TGFs) are transient gamma-ray emissions from thunderstorms, generated by electrons accelerated to relativistic energies in electric fields. Elves are ultraviolet and optical emissions excited in the lower ionosphere by electromagnetic waves radiated from lightning current pulses. We observed a TGF and an associated elve using the Atmosphere-Space Interactions Monitor on the International Space Station. The TGF occurred at the onset of a lightning current pulse that generated an elve, in the early stage of a lightning flash. Our measurements suggest that the current onset is fast and has a high amplitude-a prerequisite for elves-and that the TGF is generated in the electric fields associated with the lightning leader.

Downward Terrestrial Gamma-Ray Flash Observed in a Winter Thunderstorm

During a winter thunderstorm on 24 November 2017, a strong burst of gamma rays with energies up to ∼10  MeV was detected coincident with a lightning discharge, by scintillation detectors installed at the Kashiwazaki-Kariwa Nuclear Power Station at sea level in Japan. The burst had a subsecond duration, which is suggestive of photoneutron production. The leading part of the burst was resolved into four intense gamma-ray bunches, each coincident with a low-frequency radio pulse. These bunches were separated by 0.7-1.5 ms, with a duration of ≪1  ms each. Thus, the present burst may be considered as a "downward" terrestrial gamma-ray flash (TGF), which is analogous to upgoing TGFs observed from space. Although the scintillation detectors were heavily saturated by these bunches, the total dose associated with them was successfully measured by ionization chambers, employed by nine monitoring posts surrounding the power plant. From this information and Monte Carlo simulations, the present downward TGF is suggested to have taken place at an altitude of 2500±500  m, involving 8_{-4}^{+8}×10^{18} avalanche electrons with energies above 1 MeV. This number is comparable to those in upgoing TGFs.

Gamma Ray Glow Observations at 20-km Altitude

In the spring of 2017 an ER-2 aircraft campaign was undertaken over continental United States to observe energetic radiation from thunderstorms and lightning. The payload consisted of a suite of instruments designed to detect optical signals, electric fields, and gamma rays from lightning. Starting from Georgia, USA, 16 flights were performed, for a total of about 70 flight hours at a cruise altitude of 20 km. Of these, 45 flight hours were over thunderstorm regions. An analysis of two gamma ray glow events that were observed over Colorado at 21:47 UT on 8 May 2017 is presented. We explore the charge structure of the cloud system, as well as possible mechanisms that can produce the gamma ray glows. The thundercloud system we passed during the gamma ray glow observation had strong convection in the core of the cloud system. Electric field measurements combined with radar and radio measurements suggest an inverted charge structure, with an upper negative charge layer and a lower positive charge layer. Based on modeling results, we were not able to unambiguously determine the production mechanism. Possible mechanisms are either an enhancement of cosmic background locally (above or below 20 km) by an electric field below the local threshold or an enhancement of the cosmic background inside the cloud but then with normal polarity and an electric field well above the Relativistic Runaway Electron Avalanche threshold.

Catalog of 2017 Thunderstorm Ground Enhancement (TGE) events observed on Aragats

The natural electron accelerator in the clouds above Aragats high-altitude research station in Armenia operates continuously in 2017 providing more than 100 Thunderstorm Ground enhancements (TGEs). Most important discovery based on analysis of 2017 data is observation and detailed description of the long-lasting TGEs. We present TGE catalog for 2 broad classes according to presence or absence of the high-energy particles. In the catalog was summarized several key parameters of the TGEs and related meteorological and atmospheric discharge observations. The statistical analysis of the data collected in tables reveals the months when TGEs are more frequent, the daytime when TGEs mostly occurred, the mean distance to lightning flash that terminates TGE and many other interesting relations. Separately was discussed the sharp count rate decline and following removal of high-energy particles from the TGE flux after a lightning flash. ADEI multivariate visualization and statistical analysis platform make analytical work on sophisticated problems rather easy; one can try and test many hypotheses very fast and come to a definite conclusion allowing crosscheck and validation

Spontaneous Emergence of Space Stems Ahead of Negative Leaders in Lightning and Long Sparks

We investigate the emergence of space stems ahead of negative leaders. These are luminous spots that appear ahead of an advancing leader mediating the leader's stepped propagation. We show that space stems start as regions of locally depleted conductivity that form in the streamers of the corona around the leader. An attachment instability enhances the electric field leading to strongly inhomogeneous, bright, and locally warmer regions ahead of the leader that explain the existing observations. Since the attachment instability is only triggered by fields above 10 kV/cm and internal electric fields are lower in positive than in negative streamers, our results explain why, although common in negative leaders, space stems, and stepping are hardly observed if not absent in positive leaders. Further work is required to fully explain the streamer to leader transition, which requires an electric current persisting for timescales longer than the typical attachment time of electrons, around 100 ns.

Measurement of the Electrical Properties of a Thundercloud Through Muon Imaging by the GRAPES-3 Experiment

The GRAPES-3 muon telescope located in Ooty, India records rapid (∼10  min) variations in the muon intensity during major thunderstorms. Out of a total of 184 thunderstorms recorded during the interval of April 2011-December 2014, the one on December 1, 2014 produced a massive potential of 1.3 GV. The electric field measured by four well-separated (up to 6 km) monitors on the ground was used to help estimate some of the properties of this thundercloud, including its altitude and area that were found to be 11.4 km above mean sea level and ≥380  km^{2}, respectively. A charging time of 6 min to reach 1.3 GV implied the delivery of a power of ≥2  GW by this thundercloud that was moving at a speed of ∼60  km h^{-1}. This work possibly provides the first direct evidence for the generation of gigavolt potentials in thunderclouds that could also possibly explain the production of highest-energy (100 MeV) gamma rays in the terrestrial gamma-ray flashes.

Observationally Weak TGFs in the RHESSI Data

Terrestrial gamma ray flashes (TGFs) are sub-millisecond bursts of high energetic gamma radiation associated with intracloud flashes in thunderstorms. In this paper we use the simultaneity of lightning detections by World Wide Lightning Location Network to find TGFs in the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) data that are too faint to be identified by standard search algorithms. A similar approach has been used in an earlier paper, but here we expand the data set to include all years of RHESSI + World Wide Lightning Location Network data and show that there is a population of observationally weak TGFs all the way down to 0.22 of the RHESSI detection threshold (three counts in the detector). One should note that the majority of these are "normal" TGFs that are produced further away from the subsatellite point (and experience a 1/r 2 effect) or produced at higher latitudes with a lower tropoause and thus experience increased atmospheric attenuation. This supports the idea that the TGF production rate is higher than currently reported. We also show that compared to lightning flashes, TGFs are more partial to ocean and coastal regions than over land.

Ground-Based Observations of Terrestrial Gamma Ray Flashes Associated with Downward-Directed Lightning Leaders

Terrestrial gamma-ray flashes (TGFs) are bursts of gamma-rays initiated in the Earth’s atmosphere. TGFs were serendipitously first observed over twenty years ago by the BATSE gamma ray satellite experiment. Since then, several satellite experiments have shown that TGFs are produced in the upward negative breakdown stage at the start of intracloud lightning discharges. In this proceeding, we present ground-based observation of TGFs produced by downward negative breakdown occurring at the beginning of negative cloud-to-ground flashes.

The Terrestrial gamma-ray flashes discussed in this work were detected between 2014-2017 at ground level by the Telescope Array surface detector (TASD) together with Lightning Mapping Array (LMA) and the slow electric field antenna (SA). The TASD detector is a 700 km2 ultra high energy cosmic ray detector in the southwestern desert of Utah. It is comprised of 507 (3 m2) plastic scintillator detectors on a 1.2 km square grid. The LMA detector, a three-dimensional total lightning location system, is comprised of nine stations located within and around the array. The slow electric field antenna records the electric field change in lightning discharges.

The observed Gamma ray showers were detected in the first 1-2 ms of downward negative breakdown prior to cloud-to-ground lightning strikes. The shower sources were observed by the LMA detector at altitudes of a few kilometers above ground level. The detected energetic burst showers have a footprint on the ground typically ~ 3-5 km in diameter. The bursts comprise of several (2-5) individual pulses, each of which have a span of a few to tens of microseconds and an overall duration of several hundred microseconds. Using a forward-beamed cone of half-angle of 16 degrees, GEANT simulation studies indicate that the showers are consistent with gamma rays of 1012 - 1014 primary photons. We hypothesize that the observed terrestrial gamma-ray flashes are similar to those detected by satellites, but that the ground-based observations are closer to the source and therefore are able to observe weaker sources and report on the structure of the temporal distribution at the source. This result and future studies will enable us to better identify and constrain the mechanisms of downward TGF production.

Broadband longwave radio remote sensing instrumentation

We present the performance characteristics of a high-sensitivity radio receiver for the frequency band 0.5-470 kHz, known as the Low Frequency Atmospheric Weather Electromagnetic System for Observation, Modeling, and Education, or LF AWESOME. The receiver is an upgraded version of the VLF AWESOME, completed in 2004, which provided high sensitivity broadband radio measurements of natural lightning emissions, transmitting beacons, and radio emissions from the near-Earth space environment. It has been deployed at many locations worldwide and used as the basis for dozens of scientific studies. We present here a significant upgrade to the AWESOME, in which the frequency range has been extended to include the LF and part of the medium frequency (MF) bands, the sensitivity improved by 10-25 dB to be as low as 0.03 fT/ Hz , depending on the frequency, and timing error reduced to 15-20 ns range. The expanded capabilities allow detection of radio atmospherics from lightning strokes at global distances and multiple traverses around the world. It also allows monitoring of transmitting beacons in the LF/MF band at thousands of km distance. We detail the specification of the LF AWESOME and demonstrate a number of scientific applications. We also describe and characterize a new algorithm for minimum shift keying demodulation for VLF/LF transmitters for ionospheric remote sensing applications.

High-Energy Emissions Induced by Air Density Fluctuations of Discharges

Bursts of X-rays and γ-rays are observed from lightning and laboratory sparks. They are bremsstrahlung from energetic electrons interacting with neutral air molecules, but it is still unclear how the electrons achieve the required energies. It has been proposed that the enhanced electric field of streamers, found in the corona of leader tips, may account for the acceleration; however, their efficiency is questioned because of the relatively low production rate found in simulations. Here we emphasize that streamers usually are simulated with the assumption of homogeneous gas, which may not be the case on the small temporal and spatial scales of discharges. Since the streamer properties strongly depend on the reduced electric field E/n, where n is the neutral number density, fluctuations may potentially have a significant effect. To explore what might be expected if the assumption of homogeneity is relaxed, we conducted simple numerical experiments based on simulations of streamers in a neutral gas with a radial gradient in the neutral density, assumed to be created, for instance, by a previous spark. We also studied the effects of background electron density from previous discharges. We find that X-radiation and γ-radiation are enhanced when the on-axis air density is reduced by more than ∼25%. Pre-ionization tends to reduce the streamer field and thereby the production rate of high-energy electrons; however, the reduction is modest. The simulations suggest that fluctuations in the neutral densities, on the temporal and spacial scales of streamers, may be important for electron acceleration and bremsstrahlung radiation.

Short-pulse lasers for weather control

Filamentation of ultra-short TW-class lasers recently opened new perspectives in atmospheric research. Laser filaments are self-sustained light structures of 0.1-1 mm in diameter, spanning over hundreds of meters in length, and producing a low density plasma (10¹⁵-10¹⁷ cm^-3) along their path. They stem from the dynamic balance between Kerr self-focusing and defocusing by the self-generated plasma and/or non-linear polarization saturation. While non-linearly propagating in air, these filamentary structures produce a coherent supercontinuum (from 230 nm to 4 µm, for a 800 nm laser wavelength) by self-phase modulation (SPM), which can be used for remote 3D-monitoring of atmospheric components by Lidar (Light Detection and Ranging). However, due to their high intensity (10¹³-10¹⁴ W cm^-2), they also modify the chemical composition of the air via photo-ionization and photo-dissociation of the molecules and aerosols present in the laser path. These unique properties were recently exploited for investigating the capability of modulating some key atmospheric processes, like lightning from thunderclouds, water vapor condensation, fog formation and dissipation, and light scattering (albedo) from high altitude clouds for radiative forcing management. Here we review recent spectacular advances in this context, achieved both in the laboratory and in the field, reveal their underlying mechanisms, and discuss the applicability of using these new non-linear photonic catalysts for real scale weather control.

Spectral Characteristics of VLF Sferics Associated With RHESSI TGFs

We compared the modeled energy spectral density of very low frequency (VLF) radio emissions from terrestrial gamma ray flashes (TGFs) with the energy spectral density of VLF radio sferics recorded by Duke VLF receiver simultaneously with those TGFs. In total, six events with world wide lightning location network (WWLLN) defined locations were analyzed to exhibit a good fit between the modeled and observed energy spectral densities. In VLF range the energy spectral density of the TGF source current moment is found to be dominated by the contribution of secondary low-energy electrons and independent of the relativistic electrons which play their role in low-frequency (LF) range. Additional spectral modulation by the multiplicity of TGF peaks was found and demonstrated a good fit for two TGFs whose VLF sferics consist of two overlapping pulses each. The number of seeding pulses in TGF defines the spectral shape in VLF range, which allows to retrieve this number from VLF sferics, assuming they were radiated by TGFs. For two events it was found that the number of seeding pulses is small, of the order of 10. For the rest of the events the lower boundary of the number of seeding pulses was found to be between 10 to 103.

Electron acceleration during streamer collisions in air

High-voltage laboratory experiments show that discharges in air, generated over a gap of one meter with maximal voltage of 1 MV, may produce X-rays with photon energies up to 1 MeV. It has been suggested that the photons are bremsstrahlung from electrons accelerated by the impulsive, enhanced field during collisions of negative and a positive streamers. To explore this process, we have conducted the first self-consistent particle simulations of streamer encounters. Our simulation model is a 2-D, cylindrically symmetric, particle-in-cell code tracing the electron dynamics and solving the space charge fields, with a Monte Carlo scheme accounting for collisions and ionization. We present the electron density, the electric field, and the velocity distribution as functions of space and time. Assuming a background electric field 1.5 times the breakdown field, we find that the electron density reaches 2·10²¹ m^-3, the size of the encounter region is ∼3·10^-12 m³ and that the field enhances to ∼9 times the breakdown field during ∼10^-11 s. We further find that the radial component becomes comparable to the parallel component, which together with angular scattering leads to an almost isotropic distribution of electrons. This is consistent with laboratory observations that X-rays are emitted isotropically. However, the maximum energy of electrons reached in the simulation is ∼600 eV, which is well below the energies required to explain observations. The reason is that the encounter region is small in size and duration. For the photon energies observed, the field must be enhanced in a larger region and/or for a longer time.

Production mechanisms of leptons, photons, and hadrons and their possible feedback close to lightning leaders

It has been discussed that lightning flashes emit high-energy electrons, positrons, photons, and neutrons with single energies of several tens of MeV. In the first part of this paper we study the absorption of neutron beams in the atmosphere. We initiate neutron beams of initial energies of 350 keV, 10 MeV, and 20 MeV at source altitudes of 4 km and 16 km upward and downward and see that in all these cases neutrons reach ground altitudes and that the cross-section areas extend to several km². We estimate that for terrestrial gamma-ray flashes approximately between 10 and 2000 neutrons per ms and m² are possibly detectable at ground, at 6 km, or at 500 km altitude. In the second part of the paper we discuss a feedback model involving the generation and motion of electrons, positrons, neutrons, protons, and photons close to the vicinity of lightning leaders. In contrast to other feedback models, we do not consider large-scale thundercloud fields but enhanced fields of lightning leaders. We launch different photon and electron beams upward at 4 km altitude. We present the spatial and energy distribution of leptons, hadrons, and photons after different times and see that leptons, hadrons, and photons with energies of at least 40 MeV are produced. Because of their high rest mass hadrons are measurable on a longer time scale than leptons and photons. The feedback mechanism together with the field enhancement by lightning leaders yields particle energies even above 40 MeV measurable at satellite altitudes.

Radio emissions from double RHESSI TGFs

A detailed analysis of Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) terrestrial gamma ray flashes (TGFs) is performed in association with World Wide Lightning Location Network (WWLLN) sources and very low frequency (VLF) sferics recorded at Duke University. RHESSI clock offset is evaluated and found to experience changes on the 5 August 2005 and 21 October 2013, based on the analysis of TGF-WWLLN matches. The clock offsets were found for all three periods of observations with standard deviations less than 100 μs. This result opens the possibility for the precise comparative analyses of RHESSI TGFs with the other types of data (WWLLN, radio measurements, etc.) In case of multiple-peak TGFs, WWLLN detections are observed to be simultaneous with the last TGF peak for all 16 cases of multipeak RHESSI TGFs simultaneous with WWLLN sources. VLF magnetic field sferics were recorded for two of these 16 events at Duke University. These radio measurements also attribute VLF sferics to the second peak of the double TGFs, exhibiting no detectable radio emission during the first TGF peak. Possible scenarios explaining these observations are proposed. Double (multipeak) TGFs could help to distinguish between the VLF radio emission radiated by the recoil currents in the +IC leader channel and the VLF emission from the TGF producing electrons.

Relativistic electrons from sparks in the laboratory

Discharge experiments were carried out at the Eindhoven University of Technology in 2013. The experimental setup was designed to search for electrons produced in meter-scale sparks using a 1 MV Marx generator. Negative voltage was applied to the high voltage (HV) electrode. Five thin (1 mm) plastic detectors (5 cm² each) were distributed in various configurations close to the spark gap. Earlier studies have shown (for HV negative) that X-rays are produced when a cloud of streamers is developed 30-60 cm from the negative electrode. This indicates that the electrons producing the X-rays are also accelerated at this location, that could be in the strong electric field from counterstreamers of opposite polarity. Comparing our measurements with modeling results, we find that ∼300 keV electrons produced about 30-60 cm from the negative electrode are the most likely source of our measurements. A statistical analysis of expected detection of photon bursts by these fiber detectors indicates that only 20%-45% of the detected bursts could be from soft (∼10 keV) photons, which further supports that the majority of detected bursts are produced by relativistic electrons.

Observation of intrinsically bright terrestrial gamma ray flashes from the Mediterranean basin

We present three terrestrial gamma ray flashes (TGFs) observed over the Mediterranean basin by the Reuven Ramaty High Energy Solar Spectroscope Imager (RHESSI) satellite. Since the occurrence of these events in the Mediterranean region is quite rare, the characterization of the events was optimized by combining different approaches in order to better define the cloud of origin. The TGFs on 7 November 2004 and 16 October 2006 came from clouds with cloud top higher than 10-12 km where often a strong penetration into the stratosphere is found. This kind of cloud is usually associated with heavy precipitation and intense lightning activity. Nevertheless, the analysis of the cloud type based on satellite retrievals shows that the TGF on 27 May 2004 was produced by an unusual shallow convection. This result appears to be supported by the model simulation of the particle distribution and phase in the upper troposphere. The TGF on 7 November 2004 is among the brightest ever measured by RHESSI. The analysis of the energy spectrum of this event is consistent with a production altitude ≤12 km, which is in the upper part of the cloud, as found by the meteorological analysis of the TGF-producing thunderstorm. This event must be unusually bright at the source in order to produce such a strong signal in RHESSI. We estimate that this TGF must contain ∼3 × 10¹⁸ initial photons with energy >1 MeV. This is 1 order of magnitude brighter than earlier estimations of an average RHESSI TGF.

Enhanced detection of terrestrial gamma-ray flashes by AGILE

At the end of March 2015 the onboard software configuration of the Astrorivelatore Gamma a Immagini Leggero (AGILE) satellite was modified in order to disable the veto signal of the anticoincidence shield for the minicalorimeter instrument. The motivation for such a change was the understanding that the dead time induced by the anticoincidence prevented the detection of a large fraction of Terrestrial Gamma-Ray Flashes (TGFs). The configuration change was highly successful resulting in an increase of one order of magnitude in TGF detection rate. As expected, the largest fraction of the new events has short duration (<100 μs), and part of them has simultaneous association with lightning sferics detected by the World Wide Lightning Location Network. The new configuration provides the largest TGF detection rate surface density (TGFs/km²/yr) to date, opening prospects for improved correlation studies with lightning and atmospheric parameters on short spatial and temporal scales along the equatorial region.

Modeling the relativistic runaway electron avalanche and the feedback mechanism with GEANT4

This paper presents the first study that uses the GEometry ANd Tracking 4 (GEANT4) toolkit to do quantitative comparisons with other modeling results related to the production of terrestrial gamma ray flashes and high-energy particle emission from thunderstorms. We will study the relativistic runaway electron avalanche (RREA) and the relativistic feedback process, as well as the production of bremsstrahlung photons from runaway electrons. The Monte Carlo simulations take into account the effects of electron ionization, electron by electron (Møller), and electron by positron (Bhabha) scattering as well as the bremsstrahlung process and pair production, in the 250 eV to 100 GeV energy range. Our results indicate that the multiplication of electrons during the development of RREAs and under the influence of feedback are consistent with previous estimates. This is important to validate GEANT4 as a tool to model RREAs and feedback in homogeneous electric fields. We also determine the ratio of bremsstrahlung photons to energetic electrons N_γ /N_e . We then show that the ratio has a dependence on the electric field, which can be expressed by the avalanche time τ(E) and the bremsstrahlung coefficient α(ε). In addition, we present comparisons of GEANT4 simulations performed with a "standard" and a "low-energy" physics list both validated in the 1 keV to 100 GeV energy range. This comparison shows that the choice of physics list used in GEANT4 simulations has a significant effect on the results.

KEY POINTS

Testing the feedback mechanism with GEANT4Validating the GEANT4 programming toolkitStudy the ratio of bremsstrahlung photons to electrons at TGF source altitude.

An altitude and distance correction to the source fluence distribution of TGFs

The source fluence distribution of terrestrial gamma ray flashes (TGFs) has been extensively discussed in recent years, but few have considered how the TGF fluence distribution at the source, as estimated from satellite measurements, depends on the distance from satellite foot point and assumed production altitude. As the absorption of the TGF photons increases significantly with lower source altitude and larger distance between the source and the observing satellite, these might be important factors. We have addressed the issue by using the tropopause pressure distribution as an approximation of the TGF production altitude distribution and World Wide Lightning Location Network spheric measurements to determine the distance. The study is made possible by the increased number of Ramaty High Energy Solar Spectroscopic Imager (RHESSI) TGFs found in the second catalog of the RHESSI data. One find is that the TGF/lightning ratio for the tropics probably has an annual variability due to an annual variability in the Dobson-Brewer circulation. The main result is an indication that the altitude distribution and distance should be considered when investigating the source fluence distribution of TGFs, as this leads to a softening of the inferred distribution of source brightness.

Relativistic runaway ionization fronts

We investigate the first example of self-consistent impact ionization fronts propagating at relativistic speeds and involving interacting, high-energy electrons. These fronts, which we name relativistic runaway ionization fronts, show remarkable features such as a bulk speed within less than one percent of the speed of light and the stochastic selection of high-energy electrons for further acceleration, which leads to a power-law distribution of particle energies. A simplified model explains this selection in terms of the overrun of Coulomb-scattered electrons. Appearing as the electromagnetic interaction between electrons saturates the exponential growth of a relativistic runaway electron avalanche, relativistic runaway ionization fronts may occur in conjunction with terrestrial gamma-ray flashes and thus explain recent observations of long, power-law tails in the terrestrial gamma-ray flash energy spectrum.

Terrestrial gamma-ray flashes as powerful particle accelerators

Strong electric discharges associated with thunderstorms can produce terrestrial gamma-ray flashes (TGFs), i.e., intense bursts of x rays and γ rays lasting a few milliseconds or less. We present in this Letter new TGF timing and spectral data based on the observations of the Italian Space Agency AGILE satellite. We determine that the TGF emission above 10 MeV has a significant power-law spectral component reaching energies up to 100 MeV. These results challenge TGF theoretical models based on runaway electron acceleration. The TGF discharge electric field accelerates particles over the large distances for which maximal voltages of hundreds of megavolts can be established. The combination of huge potentials and large electric fields in TGFs can efficiently accelerate particles in large numbers, and we reconsider here the photon spectrum and the neutron production by photonuclear reactions in the atmosphere.