Empirical orthogonal function analysis and modeling of global ionospheric spherical harmonic coefficients

摘要

We developed a global empirical model for computing spherical harmonic (SH) coefficients based on the empirical orthogonal function (EOF) method and periodic functions by utilizing global ionospheric SH coefficients data provided by Center for Orbit Determination in Europe (CODE) during the years 2005-2016. Results show that the first four-order base functions and corresponding associated coefficients can represent 97.15% of the basic characteristics of the original data. The first four-order associated coefficients have noticeable periodic variations, and the correlation between these coefficients and the solar activity intensity is high. By fitting the associated coefficients with the periodic function that takes into account the influence of solar activity, it is possible to establish an EOF model to characterize the variations of the ionospheric SH coefficients with few model parameters and further calculate the global vertical total electron content (VTEC). Relative to the existing global VTEC EOF models, this EOF method can achieve high-accuracy modeling of global VTEC with a smaller number of model coefficients. By comparing the SH coefficients and the global VTEC between the EOF model and CODE, results demonstrate that the EOF model can achieve high accuracy and reliability under different solar activities. The difference between the EOF model and CODE is basically at the same level as that between the individual ionospheric analysis centers and CODE. In the extreme event of magnetic storms, the EOF model can also present higher accuracy than the international reference ionosphere (IRI) model.