Using sparse regularization for multi-resolution tomography of the ionosphere

• Main Authors: T. Panicciari, N. D. Smith, C. N. Mitchell, F. Da Dalt, P. S. J. Spencer
• Format: Article
• Language: English
• Published: Copernicus Publications 2015-10-01
• Series: Nonlinear Processes in Geophysics
• Online Access: http://www.nonlin-processes-geophys.net/22/613/2015/npg-22-613-2015.pdf
• ISSN: 1023-5809; 1607-7946

Computerized ionospheric tomography (CIT) is a technique that allows reconstructing the state of the ionosphere in terms of electron content from a set of slant total electron content (STEC) measurements. It is usually denoted as an inverse problem. In this experiment, the measurements are considered coming from the phase of the GPS signal and, therefore, affected by bias. For this reason the STEC cannot be considered in absolute terms but rather in relative terms. Measurements are collected from receivers not evenly distributed in space and together with limitations such as angle and density of the observations, they are the cause of instability in the operation of inversion. Furthermore, the ionosphere is a dynamic medium whose processes are continuously changing in time and space. This can affect CIT by limiting the accuracy in resolving structures and the processes that describe the ionosphere. Some inversion techniques are based on ℓ₂ minimization algorithms (i.e. Tikhonov regularization) and a standard approach is implemented here using spherical harmonics as a reference to compare the new method. A new approach is proposed for CIT that aims to permit sparsity in the reconstruction coefficients by using wavelet basis functions. It is based on the ℓ₁ minimization technique and wavelet basis functions due to their properties of compact representation. The ℓ₁ minimization is selected because it can optimize the result with an uneven distribution of observations by exploiting the localization property of wavelets. Also illustrated is how the inter-frequency biases on the STEC are calibrated within the operation of inversion, and this is used as a way for evaluating the accuracy of the method. The technique is demonstrated using a simulation, showing the advantage of ℓ₁ minimization to estimate the coefficients over the ℓ₂ minimization. This is in particular true for an uneven observation geometry and especially for multi-resolution CIT.