Ionospheric assimilation of radio occultation and ground-based GPS data using non-stationary background model error covariance

Ionospheric data assimilation is a powerful approach to reconstruct the 3-D distribution of the ionospheric electron density from various types of observations. We present a data assimilation model for the ionosphere, based on the Gauss–Markov Kalman filter with the International Reference Ionosphere (IRI) as the background model, to assimilate two different types of total electron content (TEC) observations from ground-based GPS and space-based FORMOSAT-3/COSMIC (F3/C) radio occultation. Covariance models for the background model error and observational error play important roles in data assimilation. The objective of this study is to investigate impacts of stationary (location-independent) and non-stationary (location-dependent) classes of the background model error covariance on the quality of assimilation analyses. Location-dependent correlations are modeled using empirical orthogonal functions computed from an ensemble of the IRI outputs, while location-independent correlations are modeled using a Gaussian function. Observing System Simulation Experiments suggest that assimilation of TEC data facilitated by the location-dependent background model error covariance yields considerably higher quality assimilation analyses. Results from assimilation of real ground-based GPS and F3/C radio occultation observations over the continental United States are presented as TEC and electron density profiles. Validation with the Millstone Hill incoherent scatter radar data and comparison with the Abel inversion results are also presented. Our new ionospheric data assimilation model that employs the location-dependent background model error covariance outperforms the earlier assimilation model with the location-independent background model error covariance, and can reconstruct the 3-D ionospheric electron density distribution satisfactorily from both ground- and space-based GPS observations.