*Yongxin Gao1, Ting Li1, Chieh-Hung Chen2, Yang-Yi Sun2, Xuemin Zhang3, Jann-Yenq Liu4
(1.School of Civil Engineering, Hefei University of Technology, 2.Institute of Geophysics and Geomatics,China University of Geosciences, 3.Institute of Earthquake Forecasting, China Earthquake Administration, 4.Graduate Institute of Space Science, National Central University,)
Keywords:Earthquake, Ionosphere/atmosphere interactions, Acoustic-gravity waves
We proposed an analytically-based method to simulate the acoustic-gravity waves generated by a point force on the Earth’s surface in the horizontally-stratified atmosphere-Earth structure. The method solves the linear momentum, continuity and adiabaticity equations in the atmosphere and elastodynamic equations in the solid Earth in the frequency-wavenumber domain. The time-domain waveforms are obtained by wavenumber integration and fast Fourier transform with respect to the frequency. Numerical simulations are conducted to investigate the properties of the acoustic-gravity waves, including both the high-frequency acoustic-mode waves and low-frequency gravity-mode waves. Simulations of the high-frequency responses show that disturbances in the atmosphere with three apparent horizontal velocities can be identified. They are, namely, the direct acoustic-mode wave generated by source traveling with the sound speed, the head wave generated by the seismic P wave traveling with apparent horizontal speed identical to the P velocity, and the head wave generated by the Rayleigh wave with a horizontal speed same to the Rayleigh wave velocity. Simulations of the low-frequency responses show that the gravity-mode wave and Lamb wave can be identified. The gravity-mode wave travels with a speed lower than the sound seed and does not reach everywhere, especially the area directly above the source. The Lamb wave travels along the Earth surface with a speed of about 310 m/s and its energy decays with the altitude. We also apply our method to explaining the Doppler sounding data observed in Taiwan area during the 2011 Tohoku M 9 earthquake, and find good agreement between the predicted signals and observed data in the arrival time and wave envelope associated with the Rayleigh wave.
This study was supported by the National Natural Science Foundation of China (grant 42174084) and Fundamental Research Funds for the Central Universities of China (grant JZ2021HGPB0058).