The simulations presented here are based on the observational data of lightning electric currents associated with the eruption of the Hunga Tonga volcano in January 2022. The response of the lithosphere (Earth)–atmosphere–ionosphere–magnetosphere system to unprecedented lightning currents is theoretically investigated at low frequencies, including ultra low frequency (ULF), extremely low frequency (ELF), and very low frequency (VLF) ranges. The electric current source due to lightning near the location of the Hunga Tonga volcano eruption has a wide-band frequency spectrum determined in this paper based on a data-driven approach. The spectrum is monotonous in the VLF range but has many significant details at the lower frequencies (ULF, ELF). The decreasing amplitude tendency is maintained at frequencies exceeding 0.1 Hz. The density of effective lightning current in the ULF range reaches the value of the order of 10−7 A/m2. A combined dynamic/quasi-stationary method has been developed to simulate ULF penetration through the lithosphere (Earth)–atmosphere–ionosphere–magnetosphere system. This method is suitable for the ULF range down to 10−4 Hz. The electromagnetic field is determined from the dynamics in the ionosphere and from a quasi-stationary approach in the atmosphere, considering not only the electric component but also the magnetic one. An analytical/numerical method has been developed to investigate the excitation of the global Schumann resonator and the eigenmodes of the coupled Schumann and ionospheric Alfvén resonators in the ELF range and the eigenmodes of the Earth–ionosphere waveguide in the VLF range. A complex dispersion equation for the corresponding disturbances is derived. It is shown that oscillations at the first resonance frequency in the Schumann resonator can simultaneously cause noticeable excitation of the local ionospheric Alfvén resonator, whose parameters depend on the angle between the geomagnetic field and the vertical direction. VLF propagation is possible over distances of 3000–10,000 km in the waveguide Earth–ionosphere. The results of simulations are compared with the published experimental data.
Read full abstract