Abstract:This work presents statistical models of the variability of plasma in the topside ionosphere based on observations made by the European Space Agency’s (ESA) Swarm satellites. The models were developed in the “Swarm Variability of Ionospheric Plasma” (Swarm-VIP) project within the European Space Agency's Swarm+4D-Ionosphere framework. The configuration of the Swarm satellites, their near-polar orbits and the data products developed, enable studies of the spatial variability of the ionosphere at multiple scale sizes. The statistical modelling technique of Generalised Linear Modelling was used to create models of both the electron density and measures of the variability of the plasma structures at horizontal spatial scales between 20 km and 100 km. Despite being developed using the Swarm data, the models provide predictions that are independent of these data. Separate models were created for low, middle, auroral and polar latitudes. The models make predictions based on heliogeophysical variables, which act as proxies for the solar and geomagnetic processes. The first and most significant term in the majority of the models was a proxy for solar activity. The most common second term varied with the latitudinal region. This was the Solar Zenith Angle (SZA) in the polar region, a measure of latitude in the auroral region, solar time in the mid-latitude region and a measure of latitude in the equatorial region. Other, less significant terms in the models covered a range of proxies for the solar wind, geomagnetic activity and location. In this paper the formulation, optimisation and evaluation of these models is discussed. The models show very little bias, with a mean error of zero to two decimal places in 14 out of 20 cases. The models capture some, but not all, of the trends present in the data, with Pearson correlation coefficients of up to 0.75 between the observations and the model predictions. The models also capture some, but not all, of the variability of the ionospheric plasma, as indicated by the precision, which ranged between 0.20 and 0.83. The addition of the thermospheric density as an explanatory variable in the models improved the precision in the polar and auroral regions. It is suggested that, if the thermosphere could be observed at a higher spatial resolution, then even more of the variability of the plasma structures could be captured by statistical models. The formulation and optimisation of the models are presented in this paper. The capability of the model in reproducing the expected climatological features of the topside ionosphere, in supporting GNSS-based ionospheric observations and the performance of the model against TIE-GCM, are provided in a companion paper (Spogli et al., 2023).

Development and Improvement of the International Reference Ionosphere with special emphasis on the topside and extension to the plasmasphere

The International Reference Ionosphere model: A review and description of an ionospheric benchmark

Validation of the IRI-2020 topside ionosphere options through in-situ electron density observations by low-Earth-orbit satellites

PyIRI: Whole‐Globe Approach to the International Reference Ionosphere Modeling Implemented in Python

On the Development of Global Plasma-Ionosphere FDTD Algorithms for Electromagnetic Calculations in the Earth-ionosphere System

Methods for Correction of the Altitude Electron-Density Profiles of the IRI Model for the Bottomside Ionosphere from Satellite Sounding Data

Ionospheric Nowcasting Over Italy Through Data Assimilation: A Synergy Between IRI UP and IONORING

Multi-source data ingestion for IRI-2020 model: a combination of ground-based and space-borne observations

IRIS 2+ : A Comprehensive Database of Stratified Thermodynamic Models in the Low Solar Atmosphere

Evaluation of GNSS-TEC Data-Driven IRI-2016 Model for Electron Density

A Dual-Layer Ionosphere Model Based on 3-D Ionospheric Constraint

A PHYSICAL MODEL FOR IONOSPHERIC VERTICAL PROFILE AND COMPARISON BETWEENIT AND IRI-90

IRI model in the problem of predicting ionospheric radio-wave propagation under conditions of high solar activity

Assessment of the accuracy of IRI-2016 model at the equatorial-, low-, mid-, and high-latitudes ionosphere along 30-45 o E meridian lines during minimum (2009) and maximum (2013) phases of solar cycle 24

Empirical Models of foF2 and hmF2 Reconstituted by Global Ionosonde and Reanalysis Data and COSMIC Observations

Statistical Models of the Variability of Plasma in the Topside Ionosphere: 1. Development and Optimisation

A new transformative framework for data assimilation and calibration of physical ionosphere-thermosphere models

Empirical Orthogonal Function Analysis and Modeling of the Topside Ionospheric and Plasmaspheric TECs

PyIRTAM: A New Module of PyIRI for IRTAM Coefficients

Understanding the ionosphere thermosphere response to solar and magnetospheric drivers: status, challenges and open issues

Integrating plasmasphere, ionosphere and thermosphere observations and models into a standardised open access research environment: The PITHIA-NRF international project