GPS-derived Total Electron Content (TEC) is an integrated quantity; hence it is difficult to relate the detection of ionospheric perturbations in TEC to a precise altitude. As TEC is weighted by the maximum ionospheric density, the corresponding altitude (hmF2) is, generally, assumed as the perturbation detection altitude. To investigate the validity of this assumption in detail, we conduct an accurate analysis of the GPS-TEC measured early ionospheric signatures related to the vertical surface displacement of the Mw 7.4 Sanriku-Oki earthquake (Sanriku-Oki Tohoku foreshock). Using 3D acoustic ray tracing model to describe the evolution of the propagating seismo-acoustic wave in space and time, we demonstrate how to infer the detection altitude of these early signatures in TEC. We determine that the signatures can be detected at altitudes up to ~130 km below the hmF2. This peculiar behaviour is attributed to the satellite line of sight (LOS) geometry and station location with respect to the source, which allows one to sound the co-seismic ionospheric signatures directly above the rupture area. We show that the early onset times correspond to crossing of the LOS with the acoustic wavefront at lower ionospheric altitudes. To support the proposed approach, we further reconstruct the seismo-acoustic induced ionospheric signatures for a moving satellite in the presence of a geomagnetic field. Both the 3D acoustic ray tracing model and the synthetic waveforms from the 3D coupled model substantiate the observed onset time of the ionospheric signatures. Moreover, our simple 3D acoustic ray tracing approach allows one to extend this analysis to azimuths different than that of the station-source line.

No abstract

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 (~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.

The amplitude asymmetry and initial polarity of seismic induced ionospheric perturbations around the epicenter are considered to be important in providing information about the rupture propagation and related vertical surface deformation. To comprehend this, we study ionospheric perturbations related to the 12 May 2015, Mw 7.3 Nepal earthquake. We model the coseismic slip associated with the event using the interferometric synthetic aperture radar derived surface deformation data. The ionospheric perturbations associated with the obtained surface deformation are explained in terms of rupture propagation, favorable geomagnetic field‐wave coupling, and satellite geometry effects. We discuss the effects of phase cancelation on the perturbation evolution for various receiver satellite line‐of‐sight configurations invoking an elementary version of satellite geometry factor. The present study thus elucidates further the role of nontectonic forcing mechanisms while identifying ground source pattern using the associated ionospheric perturbations.

The oblique‐thrust Kaikoura earthquake of Mw 7.8 that struck New Zealand on 13 November 2016 at 11:02:56 UTC (local time at 00:02:56 a.m. on 14 November 2016) was one of the most geometrically and tectonically complex earthquakes recorded onshore in modern seismology. The event ruptured in the region of multisegmented faults and propagated unilaterally northeastward for more than 170 km from the epicenter. The GPS derived coseismic surface displacements reveal a larger widespread horizontal and vertical coseismic surface offsets of ~6 m and ~2 m, respectively, with two distinct tectonic thrust zones. We study the characteristics of coseismic ionospheric perturbations based on tectonic and nontectonic forcing mechanisms and demonstrate that these perturbations are linked to two distinct surface thrust zones with rotating horizontal reinforcement trending the rupture, rather than merely to the displacements oriented along the rupture propagation direction.