Parametric decay instability (PDI) generated in milliseconds is an important physical phenomenon in ionospheric heating. Usually, numerical simulations are used to study PDI mechanisms. They can intuitively investigate the generation and development process of PDI, which is necessary in experimental studies. When simulating the PDI phenomenon through the explicit finite-difference time-domain (FDTD), the spatial scale spans from kilometers to centimeters, and the time scale needs to meet the Courant–Friedrichs–Lewy condition. Simulating the PDI phenomenon is time-consuming and difficult due to the high spatial resolution and strict restriction on the discrete time step. Although a nested mesh technique can boost the computational efficiency, the application of a parallel strategy is imperative to further improve it. In this study, we present a hybrid Message Passing Interface (MPI)/OpenMP parallelization scheme to solve the above-mentioned problems. This scheme can achieve an adaptive calculation and automatic allocation of MPI tasks and OpenMP threads, proving its flexibility and portability. Under the EISCAT background parameters, the PDI phenomenon was simulated. The results of the wave mode conversion and intense localized turbulence were identical to those of the serial program. Furthermore, a new simulation example and the effect of the cavity depth on electrostatic waves and negative ion density cavity were investigated. By utilizing the proposed parallelization scheme, the simulation time can be reduced from 70 h for the serial program to 3.6 h.

Ground-based high-frequency modulated waves can periodically heat the ionosphere and create “virtual antennas”, which can radiate extremely low frequency (ELF, 0.3–3 kHz) or very low frequency (VLF, 3–30 kHz) waves for long-distance communication. Ionospheric X-mode and O-mode heating experiments using amplitude and beat-wave (BW) modulations were conducted on 21 November 2019. Experimental results were analyzed from multiple perspectives based on data from Dynasonde, a magnetometer, stimulated electromagnetic emissions, an ELF/VLF signal receiver, and ultra-high-frequency radar. The strongest excited ELF/VLF signals in previous BW modulation heating experiments were around 8–12 kHz; however, in this experiment, no signal excited in this frequency range was observed, and the signal with the highest signal/noise ratio was at the frequency of 3517 Hz, which will aid in understanding the best communication frequency under different ionospheric backgrounds. It is well-accepted that the electron temperature changes periodically with the modulation frequency. However, we noted that the electron temperature had insufficient cooling during the O-mode modulated heating process and then increased again, resulting in a continuous electron temperature increase. We found that this was related to the change in ion composition after analyzing ion-line spectra, which will be helpful in studying the effect of modulation heating on the ionosphere background.

To investigate the impact of the ionosphere’s inherent temporal variability and the height-dependent sensitivity of field-aligned scattering on the dynamic characteristics of artificial ionospheric scatter channel, this study incorporates lower-ionospheric absorption effects into a coupled model combining ionospheric heating and field-aligned scattering. Numerical simulations of scattered signal intensity along a typical VHF propagation link at various times reveal a pronounced diurnal variation: signal intensity is significantly higher at night than during the day, exhibiting a maximum attenuation of 26.7 dB. Further analysis indicates that under the combined influence of lower-ionospheric absorption and background electron density gradients, the electron density perturbation within the artificial ionospheric scatterer emerges as the dominant factor governing the channel’s time-variant behavior. The temporal evolution of this perturbation aligns closely with the observed variations in scattered signal intensity. These findings enhance the understanding of the time-variant characteristics of artificial ionospheric scatter channels and provide a theoretical foundation for optimizing channel control strategies and improving resource utilization efficiency.

Abnormal electromagnetic radiation associated with seismic activity has been reported across a wide range of frequencies, but its primary energy is concentrated in the super-low-frequency (SLF) and extremely low-frequency (ELF) bands. To estimate the effect of the seismic radiation source, a radiator in a soil foundation was modeled as a horizontal electric dipole (HED), and the propagation characteristics of the electromagnetic fields were studied in the Earth–ionosphere cavity. The expressions of the electromagnetic fields could be obtained according to the reciprocity theorem. Therefore, a new algorithm named the numerical integral algorithm was proposed, which is suitable for both the SLF and ELF bands. The new algorithm was compared with the asymptotic approximation algorithm when the receiving point was not close to the field source and the antipode. The two algorithms were found to be in excellent agreement, confirming the validity of the new algorithm for SLF and ELF bands.