The ionospheric delay is of paramount importance to radio communication, satellite navigation and positioning. It is necessary to predict high-accuracy ionospheric peak parameters for single frequency receivers. In this study, the state-of-the-art artificial neural network (ANN) technique optimized by the genetic algorithm is used to develop global ionospheric models for predicting foF2 and hmF2. The models are based on long-term multiple measurements including ionospheric peak frequency model (GIPFM) and global ionospheric peak height model (GIPHM). Predictions of the GIPFM and GIPHM are compared with the International Reference Ionosphere (IRI) model in 2009 and 2013 respectively. This comparison shows that the root-mean-square errors (RMSEs) of GIPFM are 0.82 MHz and 0.71 MHz in 2013 and 2009, respectively. This result is about 20%-35% lower than that of IRI. Additionally, the corresponding hmF2 median errors of GIPHM are 20% to 30% smaller than that of IRI. Furthermore, the ANN models present a good capability to capture the global or regional ionospheric spatial-temporal characteristics, e.g., the equatorial ionization anomaly and Weddell Sea anomaly. The study shows that the ANN-based model has a better agreement to reference value than the IRI model, not only along the Greenwich meridian, but also on a global scale. The approach proposed in this study has the potential to be a new three-dimensional electron density model combined with the inclusion of the upcoming Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC-2) data.Remote Sens. 2020, 12, 866 2 of 17 space-borne and ground-based instruments and data processed techniques, more attention has been paid to the development and improvement of these ionospheric models [3].In the previous studies, the artificial neural network (ANN), nonlinear least squares and AdaBoost techniques have been used to predict ionospheric variations with a satisfactory accuracy [4][5][6][7]. Among them, the ANN technique has proven to be a successful tool in the modeling of ionospheric variations as well as solving the forecast problems in many geophysical applications over a single station, regional area and global scale [3,[8][9][10]. The neural network was first applied to build the foF2 model over the Grahamstown station (33 • S, 26 • E) using the data range set of 1973-1983 [11]. This model can predict the daily noon foF2 value with a root-mean-square error (RMSE) of 0.95 MHz and the monthly averaged RMSE value was 0.48 MHz. Afterward, this approach has been widely applied in the Asia/Pacific sector [5], equatorial region [12], European region [13] and the auroral zone [14]. These studies showed an improvement in the accuracy of the ANN-based models over the International Reference Ionosphere (IRI) model.With the advancement of data acquisition, more ionospheric measurements have become available for the global modelling of foF2. For instance, Oyeyemi et al. [15] used the foF2 data from 59 globally distributed ionospheric stations to establish th...