Inverting Coseismic TEC Disturbances for Neutral Atmosphere Pressure Wave

Research from the past 20 years has shown that we can detect coseismic disturbances in the total electron content (TEC) using global navigation space systems (GNSS). In the near field, TEC disturbances are created by the direct wave from rupture on the surface. This pressure wave travels through the neutral atmosphere to the ionosphere within about 10 minutes. This provides the opportunity to almost immediately characterize the source of the acoustic disturbance on the surface using methods from seismology. In populated areas, this could provide valuable information to first responders. To retrieve the surface motion amplitude information we must account for changes in the waveform caused by the geomagnetic field, motion of the satellites and the geometry of the satellites and receivers. One method is to use a transfer function to invert for the neutral atmosphere pressure wave. Gómez et al (2015) first employed an analytical model to invert for acoustic waves produced by Rayleigh waves propagating along the Earth's surface. Here, we examine the same model in the near field using the TEC disturbances from the direct wave produced by rupture at the surface. We compare results from the forward model against a numerical model that has been shown to be in good agreement with observations from the 2011 Van (Turkey) earthquake. We show the forward model predictions using both methods for the Van earthquake. We then analyze results for hypothetical events at different latitudes and discuss the reliability of the analytical model in each scenario. Gómez, D., R. Jr. Smalley, C. A. Langston, T. J. Wilson, M. Bevis, I. W. D. Dalziel, E. C. Kendrick, S. A. Konfal, M. J. Willis, D. A. Piñón, et al. (2015), Virtual array beamforming of GPS TEC observations of coseismic ionospheric disturbances near the Geomagnetic South Pole triggered by teleseismic megathrusts, J. Geophys. Res. Space Physics, 120, 9087-9101, doi:10.1002/2015JA021725.