Abstract Power line harmonic radiation (PLHR) is a continuous source of electromagnetic energy in the near‐Earth environment. However, it is currently unknown why certain PLHR frequencies are more likely to reach satellites under specific ionospheric conditions. A significant factor likely influencing this behavior is the presence of ions, which can attenuate whistler modes and introduce additional wave modes that can impact PLHR propagation characteristics. Building on previous work that studied very low frequency PLHR propagation through an electron‐only ionosphere, this study employs a three‐dimensional finite‐difference time‐domain model comprised of both electrons and ions to investigate powerline emission and PLHR behavior at 60 Hz and higher harmonics. Two key advancements are introduced relative to the earlier work: (a) two ion profiles are incorporated into the model up to 1,000 km, and their impact is examined; and (b) an enhanced collisional plasma formulation designed for magnetized plasmas containing both electrons and ions is developed for the top of the grid to avoid reflections. The model is validated against analytical solutions and applied to PLHR propagation across polar, mid‐latitude, and equatorial regions. Results demonstrate that the ion composition plays a significant role in wave transmission, particularly when the PLHR source frequency approaches or falls below the ion gyrofrequency cutoff, resulting in notable changes in the electric field amplitudes at certain altitudes. These findings offer insight into the frequency‐selective nature of PLHR detection observed in satellite data. The modeling framework is adaptable, providing a valuable tool for studying electromagnetic wave behavior in different ionospheric environments.
more »
« less
A 3‐D FDTD Methodology for Modeling the Propagation of VLF Whistler Mode PLHR Waves Through the Ionosphere
Abstract The finite‐difference time‐domain (FDTD) method was previously applied to high‐frequency electromagnetic wave propagation through 250 km of theFregion of the ionosphere. That modeling approach was limited to electromagnetic wave propagation above the critical frequency of the ionospheric plasma, and it did not include the lower ionosphere layers or the top of theF‐region. This paper extends the previous modeling methodology to frequencies below the critical frequency of the plasma and to altitudes encompassing the ionosphere. The following changes to the previous work were required to generate this model: (a) theD,Eand top of theFregions of the ionosphere were added; and (b) the perfectly matched layer absorbing boundary on the top side of the grid was replaced with a collisional plasma to prevent reflections. We apply this model to the study of extremely low frequency (ELF) and very low frequency (VLF) electric power line harmonic radiation (PLHR) through the ionosphere. The model is compared against analytical predictions and applied to PLHR propagation in polar, mid‐latitude and equatorial regions. Also, to further demonstrate the advantages of the grid‐based FDTD method, PLHR propagation through a polar cap patch with inhomogeneities is studied. The presented modeling methodology may be applied to additional scenarios in a straightforward manner and can serve as a useful tool for better tracking and studying electromagnetic wave propagation through the ionosphere at any latitude and in the presence of irregularities of any size and shape.
more »
« less
- Award ID(s):
- 2434336
- PAR ID:
- 10630935
- Publisher / Repository:
- American Geophysical Union
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Space Physics
- Volume:
- 130
- Issue:
- 2
- ISSN:
- 2169-9380
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We present FARR (Finite-difference time-domain ARRay), an open source, high-performance, finite-difference time-domain (FDTD) code. FARR is specifically designed for modeling radio wave propagation in collisional, magnetized plasmas like those found in the Earth’s ionosphere. The FDTD method directly solves Maxwell’s equations and captures all features of electromagnetic propagation, including the effects of polarization and finite-bandwidth wave packets. By solving for all vector field quantities, the code can work in regimes where geometric optics is not applicable. FARR is able to model the complex interaction of electromagnetic waves with multi-scale ionospheric irregularities, capturing the effects of scintillation caused by both refractive and diffractive processes. In this paper, we provide a thorough description of the design and features of FARR. We also highlight specific use cases for future work, including coupling to external models for ionospheric densities, quantifying HF/VHF scintillation, and simulating radar backscatter. The code is validated by comparing the simulated wave amplitudes in a slowly changing, magnetized plasma to the predicted amplitudes using the WKB approximation. This test shows good agreement between FARR and the cold plasma dispersion relations for O, X, R, and L modes, while also highlighting key differences from working in the time-domain. Finally, we conclude by comparing the propagation path of an HF pulse reflecting from the bottomside ionosphere. This path compares well to ray tracing simulations, and demonstrates the code’s ability to address realistic ionospheric propagation problems.more » « less
-
Electromagnetic ion cyclotron (EMIC) waves can scatter radiation belt electrons with energies of a few hundred keV and higher. To accurately predict this scattering and the resulting precipitation of these relativistic electrons on short time scales, we need detailed knowledge of the wave field’s spatio-temporal evolution, which cannot be obtained from single spacecraft measurements. Our study presents EMIC wave models obtained from two-dimensional (2D) finite-difference time-domain (FDTD) simulations in the Earth’s dipole magnetic field. We study cases of hydrogen band and helium band wave propagation, rising-tone emissions, packets with amplitude modulations, and ducted waves. We analyze the wave propagation properties in the time domain, enabling comparison within situobservations. We show that cold plasma density gradients can keep the wave vector quasiparallel, guide the wave energy efficiently, and have a profound effect on mode conversion and reflections. The wave normal angle of unducted waves increases rapidly with latitude, resulting in reflection on the ion hybrid frequency, which prohibits propagation to low altitudes. The modeled wave fields can serve as an input for test-particle analysis of scattering and precipitation of relativistic electrons and energetic ions.more » « less
-
Abstract Interplanetary (IP) shock‐driven sudden compression of the Earth's magnetosphere produces electromagnetic disturbances in the polar ionosphere. Several studies have examined the effects of IP shock on magnetosphere‐ionosphere coupling systems using all‐sky cameras and radars. In this study, we examine responses and drivers of the polar ionosphere following an IP shock compression on 16 June 2012. We observe the vertical drift and concurrent horizontal motion of the plasma. Observations from digisonde located at Antarctic Zhongshan station (ZHO) showed an ionospheric thickEregion ionization and associated vertical downward plasma motion atFregion. In addition, horizontal ionospheric convection reversals were observed on the Super Dual Auroral Radar Network ZHO and McMurdo radar observations. Findings suggest that the transient convective reversal breaks the original shear equilibrium, it is expected that the IP shock‐induced electric field triggers an enhanced velocity shear mapping to theEregion. The horizontal motion of the plasma was attributed to only the dusk‐to‐dawn electric field that existed during the preliminary phase of sudden impulse. We also found that ionospheric convection reversals were driven by a downward field‐aligned current. The results of these observations reveal, for the first time, the immediate and direct cusp ionosphere response to the IP shock, which is critical for understanding the global response of the magnetosphere following an abrupt change in Interplanetory Magnetic Field (IMF) and solar wind conditions.more » « less
-
The importance of ElectroMagnetic Ion Cyclotron (EMIC) ultra-low-frequency (ULF) waves (and their Pc1 counterparts) is connected to their critical role in triggering energetic particle precipitation from the magnetosphere to the conjugated ionosphere via pitch angle scattering. In addition, as a prominent element of the ULF zoo, EMIC/Pc1 waves can be considered a perfect tool for the remote diagnosis of the topologies and dynamic properties of near-Earth plasmas. Based on the availability of a comprehensive set of instruments, operating on the ground and in the top-side ionosphere, the present case study provides an interesting example of the evolution of EMIC propagation to both ionospheric hemispheres up to the polar cap. Specifically, we report observations of Pc1 waves detected on 30 March 2021 under low Kp, low Sym-H, and moderate AE conditions. The proposed investigation shows that high-latitude ground magnetometers in both hemispheres and the first China Seismo-Electromagnetic Satellite (CSES-01) at a Low Earth Orbit (LEO) detected in-synch Pc1 waves. In strict correspondence to this, energetic proton precipitation was observed at LEO with a simultaneous appearance of an isolated proton aurora at subauroral latitudes. This supports the idea of EMIC wave-induced proton precipitation contributing to energy transfer from the magnetosphere to the ionosphere.more » « less
An official website of the United States government

