Seismic waves excited before and after earthquakes can make vertical motion of the Earth’s surface. The perturbations can propagate to the middle and upper atmosphere and cause variations and oscillations of atmospheric characteristics at different altitudes. Among the mechanisms of such ionospheric perturbations could be propagation of acoustic-gravity waves (AGWs) in the atmosphere caused by seismic excitations at the Earth’s surface. In the present study, we use the recently developed nonlinear numerical Whole-altitude Acoustic-Gravity Wave Model (WAGWM) to simulate propagation of supersonic wave modes from the ground to the upper atmosphere.