Articles | Volume 22, issue 5
Nonlin. Processes Geophys., 22, 613–624, 2015
https://doi.org/10.5194/npg-22-613-2015
Nonlin. Processes Geophys., 22, 613–624, 2015
https://doi.org/10.5194/npg-22-613-2015

Research article 19 Oct 2015

Research article | 19 Oct 2015

Using sparse regularization for multi-resolution tomography of the ionosphere

T. Panicciari, N. D. Smith, C. N. Mitchell, F. Da Dalt, and P. S. J. Spencer T. Panicciari et al.
  • Department of Electronic and Electrical Engineering, University of Bath, Bath, BA1 7AY, UK

Abstract. 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 ℓ2 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 ℓ1 minimization technique and wavelet basis functions due to their properties of compact representation. The ℓ1 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 ℓ1 minimization to estimate the coefficients over the ℓ2 minimization. This is in particular true for an uneven observation geometry and especially for multi-resolution CIT.