The study of ionospheric disturbances associated with the two large strike-slip earthquakes in Indonesia was investigated, which are West Sumatra on 2 March 2016 (Mw = 7.8), and Palu on 28 September 2018 (Mw = 7.5). The anomalies were observed by measuring co-seismic ionospheric disturbances (CIDs) using the Global Navigation Satellite System (GNSS). The results show positive and negative CIDs polarization changes for the 2016 West Sumatra earthquake, depending on the position of the satellite line-of-sight, while the 2018 Palu earthquake shows negative changes only due to differences in co-seismic vertical crustal displacement. The 2016 West Sumatra earthquake caused uplift and subsidence, while the 2018 Palu earthquake was dominated by subsidence. TEC anomalies occurred about 10 to 15 min after the two earthquakes with amplitude of 2.9 TECU and 0.4 TECU, respectively. The TEC anomaly amplitude was also affected by the magnitude of the earthquake moment. The disturbance signal propagated with a velocity of ~1–1.72 km s⁻¹ for the 2016 West Sumatra earthquake and ~0.97–1.08 km s⁻¹ for the 2018 Palu mainshock earthquake, which are consistent with acoustic waves. The wave also caused an oscillation signal of ∼4 mHz, and their azimuthal asymmetry of propagation confirmed the phenomena in the Southern Hemisphere. The CID signal could be identified at a distance of around 400–1500 km from the epicenter in the southwestern direction.

TEC prediction models for 24 hours ahead have been developed from JOG2 GPS TEC data during 2016. Eleven month of TEC data were used as a training model of the radial basis function neural network (RBFNN) and 1 month of last data (December 2016) is used for the RBFNN model testing. The RBFNN inputs are the previous 24 hour TEC data and the minimum of Dst index during the previous 24 hours. Outputs of the model are 24 ahead TEC prediction. Comparison of model prediction show that the RBFNN model is able to predict the next 24 hours TEC is more accurate than the TEC GIM model.