Mapping the Impact of Non-Tectonic Forcing mechanisms on GNSS measured Coseismic Ionospheric Perturbations

Global Navigation Satellite System (GNSS) measured Total Electron Content (TEC) is now widely used to study the near and far-field coseismic ionospheric perturbations (CIP). The generation of near field (similar to 500-600 km surrounding an epicenter) CIP is mainly attributed to the coseismic crustal deformation. The azimuthal distribution of near field CIP may contain information on the seismic/tectonic source characteristics of rupture propagation direction and thrust orientations. However, numerous studies cautioned that before deriving the listed source characteristics based on coseismic TEC signatures, the contribution of non-tectonic forcing mechanisms needs to be examined. These mechanisms which are operative at ionospheric altitudes are classified as the i) orientation between the geomagnetic field and tectonically induced atmospheric wave perturbations ii) orientation between the GNSS satellite line of sight (LOS) geometry and coseismic atmospheric wave perturbations and iii) ambient electron density gradients. So far, the combined effects of these mechanisms have not been quantified. We propose a 3D geometrical model, based on acoustic ray tracing in space and time to estimate the combined effects of non-tectonic forcing mechanisms on the manifestations of GNSS measured near field CIP. Further, this model is tested on earthquakes occurring at different latitudes with a view to quickly quantify the collective effects of these mechanisms. We presume that this simple and direct 3D model would induce and enhance a proper perception among the researchers about the tectonic source characteristics derived based on the corresponding ionospheric manifestations.