Accurate specification of ionization production by energetic electron precipitation is critical for atmospheric chemistry models to assess the resultant atmospheric effects. Recent model‐observation comparison studies have increasingly highlighted the importance of considering precipitation fluxes in the full range of electron energy and pitch angle. However, previous parameterization methods were mostly proposed for isotropically precipitation electrons with energies up to 1 MeV, and the pitch angle dependence has not yet been parameterized. In this paper, we first characterize and tabulate the atmospheric ionization response to monoenergetic electrons with different pitch angles and energies between ∼3 keV and ∼33 MeV. A generalized method that fully accounts for the dependence of ionization production on background atmospheric conditions, electron energy, and pitch angle has been developed based on the parameterization method of Fang et al. (2010, https://doi.org/10.1029/2010GL045406). Moreover, we validate this method using 100 random atmospheric profiles and precipitation fluxes with monoenergetic and exponential energy distributions, and isotropic and sine pitch angle distributions. In a suite of 6,100 validation tests, the error in peak ionization altitude is found to be within 1 km in 91% of all the tests with a mean error of 2.7% in peak ionization rate and 1.9% in total ionization. This method therefore provides a reliable means to convert space‐measured precipitation energy and pitch angle distributions into ionization inputs for atmospheric chemistry models.

Terrestrial Very‐Low‐Frequency (VLF) energy from both lightning discharges and radio transmitters has a role in affecting the energetic electrons in the Van Allen radiation belts, but quantification of these effects is particularly difficult, largely due to the collisional damping experienced in the highly variable electron density in the D‐ and E‐region ionosphere. The Faraday International Reference Ionosphere (FIRI) model was specifically developed by combining lower‐ionosphere chemistry modeling with in situ rocket measurements, and represents to date the most reliable source of electron density profiles for the lower ionosphere. As a full‐resolution empirical model, FIRI is not well suited to D‐ and E‐region ionosphere inversion, and its applicability in transionospheric VLF simulation and in remote sensing of the lower ionosphere is limited. Motivated by how subionospheric VLF remote sensing has been aided by the Wait and Spies (WS) profile (Wait & Spies, 1964), in this study, we parameterize the FIRI profiles and extend the WS profile to the E‐region ionosphere by introducing two new parameters: the knee altitude h_k and the sharpness parameter for the E‐region ionosphere β_E. Using this modified WS profile, we calculate the expected signals at different receiver locations from the NAA, NPM, and NWC transmitters under the full range of possible ionospheric conditions. We also describe and validate a method about how these results can be readily used to translate VLF measurements into estimates of the lower ionosphere electron density. Moreover, we use this method to evaluate the sensitivity of different ground receiver locations in lower‐ionosphere remote sensing.

The Gamma‐ray Burst Monitor (GBM) onboard the Fermi spacecraft has observed many tens of sufficiently bright events, which are suitable for individual analysis. In our previous study, we fit individual, bright terrestrial gamma‐ray flashes (TGFs) with Relativistic Runaway Electron Avalanche (RREA) models for the first time. For relativistic‐feedback‐based models, the TGF‐producing electrons, which are seeded internally by a positive feedback effect, are usually accelerated in a large‐scale field with fully developed RREAs. Alternatively, lightning leader models may apply to either a large‐scale thunderstorm fields with fully developed RREAs or to inhomogeneous fields in front of lightning leaders where RREAs only develop partially. The predictions of the latter, inhomogeneous models for the TGF‐beaming geometry show some differences from estimations of the relativistic feedback models in homogeneous fields. In this work, we analyze a large sample of 66 bright Fermi GBM TGFs in the framework of lightning leader models, making comparisons with previous results from the homogeneous‐field RREA models. In most cases, the spectral analysis does not strongly favor one mechanism over the other, with 59% of the TGF events being best fit with the fully developed RREA mechanism, which corresponds to high‐potential leader models. The majority of the GBM‐measured TGFs can be best fit if the source altitude is below 15 km and 70% of events best fit by leader models cannot be satisfactorily modeled unless a tilted photon beam is used. For several spectrally soft TGFs, the tilted beam low‐potential leader model can best fit the data.

Radiation safety in the Earth's atmosphere is of particular importance to our living environment, especially at aviation altitudes. Aviation radiation has been long known to originate primarily from the galactic and solar system: galactic cosmic rays and solar energetic protons. Recent flight measurements by the Automated Radiation Measurements for Aerospace Safety experiment have uncovered another potential source for aviation radiation: Relativistic Electron Precipitation (REP) from the Van Allen radiation belts. REP can induce radiation at aviation altitudes through bremsstrahlung X‐ray production, which carries radiation down to the stratosphere and even the troposphere. In this study, using a suite of physics‐based Monte Carlo models, we characterize the effective radiation dose produced at altitudes between ground and low‐Earth‐orbit by relativistic precipitation electrons with energies between 100 keV and 10 MeV. We produce a lookup table of atmospheric radiation production that calculates the expected radiation dose for a given precipitation flux, spectrum, and pitch angle distribution. This lookup table provides results that are consistent with X‐ray measurements during radiation belt precipitation by balloon‐borne instruments in the stratosphere, and can be directly used to convert space‐borne measurements of precipitation fluxes into aviation radiation. This study represents our first attempt toward better understanding of REP's role in the atmospheric high‐altitude radiation environment.

Solar flares can severely disturb the Earth's ionosphere, leading to changes of atmospheric chemistry and degradation of Global Navigation Satellite System (GNSS) signals. The key to better assess these effects is to more accurately quantify the electron density enhancements induced by solar flares, which is still uncertain. Thanks to the good coverage and quiet electromagnetic environment provided by the Great Wall Station (GWS, 62.22°S, 58.96°W) in Antarctica, our Very‐Low‐Frequency (VLF) receiver has recorded clean VLF signatures produced by three X‐class solar flares from March to May in 2022. These measurements motivate us to revisit the uncertain question related to the ionospheric effects induced by solar flares. Different from previous studies, we model solar flare events from satellite measurements of X‐ray fluxes to ground‐based measurements of VLF transmitter signals by combining a flare irradiance model, an X‐ray propagation model, a D‐region chemistry model, and a VLF propagation model. For eight transmitting frequencies and propagation paths, modeling results can satisfactorily explain our VLF measurements at GWS during the three X‐class flares, especially so for the X1.3 event occurring on 30 March 2022. Moreover, the electron density results that reconcile space‐ and ground‐based measurements suggest that the VLF reflection height can be lowered to ∼66 km by an X1.1 flare, in line with previous studies. The present set of numerical models can be thus utilized in future studies to improve the nowcasting capability of VLF sensing technique for monitoring solar flares, as well as the accuracy of GNSS navigation/positioning during solar flares.