[1] In this paper we present modeling studies of air heating by electrical discharges in a wide range of pressures. The developed model is capable of quantifying the different contributions for heating of air at the particle level and rigorously accounts for the vibration-dissociation-vibration coupling. The model is validated by calculating the breakdown times of short air gaps and comparing to available experimental data. Detailed discussion on the role of electron detachment in the development of the thermal-ionizational instability that triggers the spark development in short air gaps is presented. The dynamics of fast heating by quenching of excited electronic states is discussed and the scaling of its main channels with ambient air density is quantified. The developed model is employed to study the streamer-to-leader transition process and to obtain its scaling with ambient air density. Streamer-to-leader transition is the name given to a sequence of events occurring in a thin plasma channel through which a relatively strong current is forced through, culminating in heating of ambient gas and increase of the electrical conductivity of the channel. This process occurs during the inception of leaders (from sharp metallic structures, from hydrometeors inside the thundercloud, or in virgin air) and during their propagation (at the leader head or during the growth of a space leader). The development of a thermal-ionizational instability that culminates in the leader formation and propagation is characterized by a change in air ionization mechanism from electron impact to associative ionization and by contraction of the plasma channel. The introduced methodology for estimation of leader speeds shows that the propagation of a leader is limited by the air heating of every newly formed leader section. It is demonstrated that the streamer-to-leader transition time has an inverse-squared dependence on the ambient air density at near-ground pressures, in agreement with similarity laws for Joule heating in a streamer channel. Model results indicate that a deviation from this similarity scaling occurs at very low air densities, where the rate of electronic power deposition is balanced by the channel expansion, and air heating from quenching of excited electronic states is very inefficient. These findings place a limit on the maximum altitude at which a hot and highly conducting lightning leader channel can be formed in the Earth's atmosphere, result which is important for understating of the gigantic jet (GJ) discharges between thundercloud tops and the lower ionosphere. Simulations of leader speeds at GJ altitudes demonstrate that initial speeds of GJs are consistent with the leader propagation mechanism. The simulation of a GJ, escaping upward from a thundercloud top, shows that the lengthening of the leader streamer zone, in a medium of exponentially decreasing air density, determines the existence of an altitude at which the streamer zones of GJs become so long that they dynamically extend (jump) all the way to th...

To date the true nature of initial breakdown pulses (IBPs) and narrow bipolar events (NBEs) in lightning discharges remains a mystery. Recent experimental evidence has correlated IBPs to the initial development of lightning leaders inside the thundercloud. NBE wideband waveforms resemble classic IBPs in both amplitude and duration. Most NBEs are quite peculiar in the sense that very frequently they occur in isolation from other lightning processes. The remaining fraction, 16% of positive polarity NBEs, according to Wu et al. (2014), happens as the first event in an otherwise regular intracloud lightning discharge. These authors point out that the initiator type of NBEs has no difference with other NBEs that did not start lightning, except for the fact that they occur deeper inside the thunderstorm (i.e., at lower altitudes). In this paper, we propose a new physical mechanism to explain the source of both IBPs and NBEs. We propose that IBPs and NBEs are the electromagnetic transients associated with the sudden (i.e., stepwise) elongation of the initial negative leader extremity in the thunderstorm electric field. To demonstrate our hypothesis a novel computational/numerical model of the bidirectional lightning leader tree is developed, consisting of a generalization of electrostatic and transmission line approximations found in the literature. Finally, we show how the IBP and NBE waveform characteristics directly reflect the properties of the bidirectional lightning leader (such as step length, for example) and amplitude of the thunderstorm electric field.

Electromagnetic ion cyclotron (EMIC) waves have long been considered to be a significant loss mechanism for relativistic electrons. This has most often been attributed to resonant interactions with the highest amplitude waves. But recent observations have suggested that the dominant energy of electrons precipitated to the atmosphere may often be relatively low, less than 1 MeV, whereas the minimum resonant energy of the highest amplitude waves is often greater than 2 MeV. Here we use relativistic electron test particle simulations in the wavefields of a hybrid code simulation of EMIC waves in dipole geometry in order to show that significant pitch angle scattering can occur due to interaction with low‐amplitude short‐wavelength EMIC waves. In the case we examined, these waves are in the H band (at frequencies above the He+ gyrofrequency), even though the highest amplitude waves were in the He band frequency range (below the He+ gyrofrequency). We also present wave power distributions for 29 EMIC simulations in straight magnetic field line geometry that show that the high wave number portion of the spectrum is in every case mostly due to the H band waves. Though He band waves are often associated with relativistic electron precipitation, it is possible that the He band waves do not directly scatter the sub‐megaelectron volts (sub‐MeV) electrons, but that the presence of He band waves is associated with high plasma density which lowers the minimum resonant energy so that these electrons can more easily resonate with the H band waves.

Simultaneous data from two interferometers separated by 16 km and synchronized within 100 ns were collected for a thunderstorm near Langmuir Lab on October 23, 2018. Analysis via triangulation followed by a least squares fit to time of arrival across all six antennae produced a three‐dimensional interferometer (3DINTF) data set. Simultaneous Lightning Mapping Array data enabled an independent calculation of 3DINTF accuracy, yielding a median location uncertainty of 200 m. This is the most accurate verified result to date for a two‐station interferometer. The 3D data allowed profiling the velocity of multiple dart leaders and K leaders that followed the same channel. 3D velocities calculated from the in‐cloud initiation site to ground ranged from 3 × 106 to 20 × 106 m/s. Average velocity generally increased with subsequent leaders, consistent with increased conditioning of the channel. Also, all leaders showed a factor of 2–3 decrease in velocity as they proceeded over 15 km of channel. We speculate that the velocity decrease is consistent with energy lost in the reionization of the channel at the leader tip. This paper includes an appendix providing details of the triangulation technique used.

We present evidence of ionospheric optical signatures of lightning, known as elves, which sometimes occur in pairs separated in time by ∼80–160 μs. We demonstrate that these “elve doublets” are the ionospheric signature of compact intracloud discharges (CIDs), which are extremely powerful, compact discharges that are thought to occur near the tops of thunderclouds. In this paper, using simple geometric calculations and full electromagnetic simulations, we show that CIDs from altitudes 14–22 km explain the time delay observed in elve doublets, consistent with typical CID altitudes. Furthermore, we show that the relative brightness of the first and second elves in the doublet is likely related to the orientation of the CID, and angles of 5°– 20° with respect to the vertical are consistent with the observed brightness ratios.