The ionosphere is part of the upper atmosphere of Earth where EUV, X-ray radiations, and energetic particles from the sun ionize the atoms and molecules. It is a region where free electrons occur in a high enough density to influence the propagation of radio waves with measurable effects. It is one of the highly variable regions of the Earth-upper atmosphere. Nowadays, society has become reliant on systems that work with communication of the satellites whose performance is dependent on the state of the ionosphere. The ionosphere causes a significant error in single-frequency GPS navigation (Davies, 1990;Tiwari et al., 2013). The temporal and spatial variation of electron density (total electron content, TEC) of the ionosphere affects the GPS signal propagation by introducing delay (slowing down and bending propagation) and scintillation (Jakowski et al., 2012;Opperman et al., 2007;Tulunay et al., 2004Tulunay et al., , 2006. Understanding the ionosphere dynamics and developing a model for nowcasting and forecasting its TEC might mitigate the ionosphere effect. In this regard, the interest in specifying the ionosphere TEC is not only for climatic-like situations but also for weather-like behavior.Global empirical models such as IRI (Bilitza, 2018), NeQuick2 (Nava et al., 2008), and NTCM (Jakowski et al., 2011) are good in capturing the median behavior of the ionosphere. These models have been designed for scientific study and radio communication applications. For example, NeQUick2 is adopted for TEC prediction to improve communication quality for International Telecommunication Union Radiocommunication Sector (ITU-R). It is also approved to be used for single frequency ionospheric error correction in the framework of GALILEO (Nigussie et al., 2012 and references therein). These models have shown different performance in their standard application in describing the observed TEC; for example, Nigussie et al. ( 2013) have compared the performance of NeQuick2 and IRI-2007 models in describing the East-African ionospheric TEC and found, for both models, the modeled and experimental VTEC has shown significant discrepancy. The performance of NTCM in its standard applications has shown also a significant discrepancy to the experimental VTEC (Getahun & Nigussie, 2021). The weak performances of these models in their standard applications indicate limitations in reproducing and forecasting the local and instantaneous conditions, especially in regions of unique irregularities are common (Scherliess et al., 2004;Sojka et al., 2007).