Next Article in Journal
Improving Spaceborne GNSS-R Algal Bloom Detection with Meteorological Data
Next Article in Special Issue
Design and Testing of an Adaptive In-phase Magnetometer (AIMAG), the Equatorial-Electrojet-Detecting Fluxgate Magnetometer, for the CAS500-3 Satellite
Previous Article in Journal
Unlocking the Potential of Data Augmentation in Contrastive Learning for Hyperspectral Image Classification
Previous Article in Special Issue
Ionospheric–Thermospheric Responses to Geomagnetic Storms from Multi-Instrument Space Weather Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Method of Ionospheric Inversion Based on Horizontal Constraint and Empirical Orthogonal Function

School of Geography and Remote Sensing, Guangzhou University, Guangzhou 510006, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(12), 3124; https://doi.org/10.3390/rs15123124
Submission received: 28 May 2023 / Revised: 13 June 2023 / Accepted: 13 June 2023 / Published: 15 June 2023
(This article belongs to the Special Issue Satellite Observations of the Global Ionosphere and Plasma Dynamics)

Abstract

:
Tomographic inversion of the ionosphere is a rank-deficient problem. To overcome the above problem, an algebraic reconstruction technique (ART) based on adaptive horizontal constraint and empirical orthogonal function (ARTHCEOF) is proposed. The new algorithm avoids the difficulty of vertically constrained matrix construction and resolves the description of the ionospheric vertical structure by using EOF. To confirm the feasibility and validate the ascendancy of the ARTHCEOF, three algorithms are first tested by using the emulated scheme. The test results show that the ARTHCEOF surpasses the ART and the ART based on the horizontal constraint (ARTHC) in both the inversion accuracy and computational efficiency. Finally, the ARTHCEOF method is applied to invert electron density values using the GNSS measurements during different geomagnetic activities. The tomographic images validate that the ARTHCEOF can reflect ionospheric daily changes in the European region. The altitudinal profiles in a fixed location are illustrated according to the inversion results of ARTHCEOF. Compared with the profiles recorded by the ionosonde station, the altitudinal profiles of ARTHCEOF have a good consistency. In the meantime, the VTEC values are inverted using the CIT results. The differential VTEC values are calculated by means of the inverted VTEC values and ionospheric products of CODE. The differential results further identify the dependability of ARTHCEOF.

1. Introduction

The ionosphere is an important component of the atmosphere above the Earth’s surface. The aggregation of free electrons impacts the propagation speed of satellite signals [1,2]. This leads to an extension of the transmission time of radio waves, which is referred to as ionospheric delay. For a high-accuracy user, the existence of ionospheric delay error affects the accuracy of satellite positioning [3,4]. To meet the demands of high-accuracy users, it is particularly necessary to detect three-dimensional ionospheric fine structures. The constructions of the global navigation satellite system (GNSS) provide a new ionospheric sounding method. Using the GNSS measurements, ionospheric vertical total electron content (VTEC) values are computed. Although the VTEC images can reflect ionospheric horizontal variations, it neglects ionospheric vertical structures since VTEC values are computed using a thin-layer hypothetical ionospheric model [5,6,7].
In 1988, the computerized ionospheric tomography (CIT) technique was introduced to probe ionospheric structure by using the measurements of the navy navigation satellite system (NNSS) [8]. Although NNSS-based CIT successfully captured ionospheric vertical structures, it neglected the ionospheric variations in longitudinal direction [9,10,11]. Using the emulated GNSS data, three-dimensional CIT results were successfully inverted [12]. Then, the GNSS-based CIT technique attracted the attention of some scholars and became the research highlight in the ionospheric field [13,14,15,16,17,18,19]. In general, the GNSS-based CIT technique is a rank-deficient problem due to the inadequacy of GNSS measurements. To surmount the rank deficiency, some algorithms have been proposed [20,21,22,23,24,25,26]. The algebraic reconstruction technique (ART) is very popular due to its computational simplicity. However, some pixels do not have any GNSS observation information, which results in the final solutions of these pixels being the same as their priori information after iteration convergence. Some scholars introduced different horizontal constraints to overcome the deficiency of ART to some degree [27,28,29]. The type of the algorithms is ART based on the different horizontal constraints (ARTHC). However, the studies omit the vertical constraint. This decreases the accuracy of ionospheric inversion. In practice, the vertical constraint matrix is difficult to create due to the complex variation characteristic of altitude. Empirical orthogonal function (EOF) has the advantage of describing the ionospheric vertical structure. Instead of imposing vertical constraint, EOF is used to establish the correlation of adjacent pixels in the vertical direction. The ARTHC based on EOF (ARTHCEOF) is proposed. To examine the ARTHCEOF method, the emulated experiment is made. Finally, the ARTHCEOF method is applied to reconstruct ionospheric electron density (IED) distributions by using the GNSS measurements during different geomagnetic activities. The reconstructed images show the daily variation rules of the IED. The comparisons of ionospheric vertical profiles and the differential VTEC values validate the dependability of the ARTHCEOF.

2. Materials and Methods

2.1. Basic Principle of CIT

For IED tomographic inversion, the input observation data are the ionospheric slant TEC (STEC) value, which is the IED integral along the propagation paths of the GNSS signal. The formula can be given as:
S T E C = p N ( l ) d l
where p is GNSS ray path; N ( l ) represents the IED; d l is an extremely short segment on the ray propagation path. To simplify IED inversion, the reconstructed area is first divided into some small pixels. The IED of each pixel is assumed to be constant in the inverted time period, and then Formula (1) is linearized as:
y i = j = 1 n A i j x j i = 1 , 2 , , m
Taking into account the input STEC noise and discretized error, the matrix expression of Formula (2) is as follows:
Y m × 1 = A m × n X n × 1 + E m × 1
where n stands for the pixel amount; m denotes the amount of input STEC; Y is STEC vector; A is the constructed matrix; X is the column vector of IED; E is the error column vector.

2.2. Theory of ARTHCEOF

ART is one of the classical reconstruction algorithms of IED tomographic inversion. In this work, the iterative initial values are given by using the NeQuick model in order to run ART. It is called a round iteration when all available STEC measurements involve an iteration, and then STEC values are computed by using the iterative result of the last round. The difference can be obtained using the input and computed STEC data. ART uses the difference to correct the IED before iterative convergence. The equation of ART can be given as:
x j ( k + 1 ) = x j ( k ) + γ 0 y i a i x j ( k ) a i 2 a i T
where x j ( k ) is the k th iterative result of the j th pixel; a i represents the i th row of A; γ 0 is the relaxation operator affecting the iterative speed and reconstructed quality of ART; y i indicates the i th row of the input STEC vector
In light of Equation (4), ART can modify the IED priori information of each pixel having observation information. In fact, many pixels do not have any GNSS ray traversing them because of the sparse and uneven distribution of ground-based GNSS stations. Their IED priori values cannot be modified in the iterative process. So, the reconstructed accuracy of ART is affected. To cope with the above problem, some scholars impose various constraints in the horizontal direction [27,28,29]. However, the constraint matrices remain immutably in each round iteration. In addition, vertical constraints are also neglected. To improve the previous work, its elements of the horizontal constraint matrix are anticipated to modify the IED iterative results. The adaptive constrained function is introduced to modify the elements of the horizontal constraint matrix in this work. The horizontal constrained formula can be represented as:
f 1 x 1 + + f d 1 x d 1 x d + f d + 1 f d + 1 + + f n x n = 0
The constraint coefficient f i can be given as:
f j ( k + 1 ) = e ( b c d ( k ) ) 2 / 2 σ 2 c = 1 n e d = 1 n n e ( b c d ( k ) ) 2 / 2 σ 2
b c d ( k ) = D c d × x d ( k ) x c ( k )
where ne is the pixel amount in the longitudinal direction; n n represents the pixel amount in the latitudinal direction; D c d is the space of the center of two pixels; σ is the smoothing operator.
In light of Formula (6), the constraint matrix in the horizontal direction is written as:
F = d i a g ( F 1 , F 2 , , F g )
F g represents the constrained matrix of the g th layer. The constrained formula is as follows:
F X = 0
Combining Formula (3) with Formula (9), a new formula can be obtained.
[ Y 0 ] = [ A F ] X + E
The matrix form of Formula (10) is given as:
B ( m + h ) × n X n × 1 + E ( m + h ) × 1 = L ( m + h ) × 1
where h is the number of constrained equations.
Considering the stratification and gradient features of the ionosphere in altitude, the construction of the vertical constraint matrix is difficult. So, EOF is introduced to describe the vertical variation of the ionosphere. Substituting the orthogonal matrix composed of EOF vertical basis functions into Equation (11), the following equation can be obtained.
B ( m + h ) × n U n × p U n × p T X n × 1 + E ( m + h ) × 1 = L ( m + h ) × 1
p is related to the number of horizontal pixels and EOF. U is a matrix composed of EOF.
In general, the first three orders of EOF can basically describe 95% ionospheric profiles. Figure 1 shows the first three orders of EOF in European regions. From Figure 1, it can be seen that the first three orders of EOF include 79.88 percent, 10.08 percent, and 5.47 percent information on the ionospheric vertical profile, respectively.

3. Results

3.1. Simulated Test of the New CIT Algorithm

To verify the practicability and advantages of ARTHCEOF, numerical examinations of three CIT algorithms are first carried out. The examining results of three algorithms are compared with the emulated IED truth values. In the emulation scheme, the latitudinal range is 36°–60°N and the longitudinal range is 0°–32°E. The scope changes from 100 km to 1000 km in height. Nesterov and Kunitsyn [30] validated that the spatial resolution can be selected as 2°–4° when GNSS stations are dense. To construct the coefficient matrix of CIT, the spatial coordinates of 183 European GNSS stations are used. The geographic locations of the GNSS and ionsonde stations are illustrated in Figure 2, which shows that the selected European GNSS stations are dense. So, the longitudinal step is 4°, and the latitudinal step is 2°. In the altitudinal direction, the interval is 20 km.
The geomagnetic activity is calm at the time period of 13:30–14:00UT on 20 August 2018. The GNSS data of the above time period are selected to examine three algorithms. The IED background values are necessary to run iterative CIT algorithms. In this work, the IED background values are obtained from the NeQuick model. To distinguish the IED background value, the IED true values are given by IRI empirical model. Figure 3 illustrates the images of the IED background information in different transverse sections, and it shows that the imaging distortion is serious.
Figure 4 and Figure 5 compare the inversion results of three tomographic algorithms with those reconstructed by the IRI 2016 model in different cross-sections. The comparisons confirm that the inversion results of the ARTHCEOF method is more consistent with the IED true values than those of the ART and ARTHC, and the image distortions of the ART and ARTHC are more severe than those of the ARTHCEOF. The reason is that the ART and ARTHC cannot accurately describe the vertical variation of the ionospheric structure. The inversion results reveal that the IED values gradually decrease with the increasing latitude and longitude at 14:00UT on 20 August 2018.
To evaluate the reconstructed accuracy, the selection of evaluating indicators is very important. Three evaluating indicators are introduced in this work. One is the IED difference (IEDD), another is the mean absolute error (MAE), and the last is the root mean square error (RMSE). Three indicators are calculated using the following formulas:
IEDD = x j t o m o x j t r u e
MAE = j = 1 n | x j t o m o x j t r u e | / n
RMSE = j = 1 n ( x j t o m o x j t r u e ) 2 / n
where x j t r u e is the emulated IED value of the j th pixel, and x j t o m o is the CIT result of the j th pixel.
The IEDD of three CIT methods is first calculated. Figure 6 counts the number of pixels in each error histogram. Figure 6c reveals that the amount of pixels with small IEDD is significantly larger than those in Figure 6a and Figure 6b. Table 1 shows the error indicators and iterative efficiency of three CIT algorithms after iteration convergence. The census indicates that the inversion accuracy and efficiency of the ARTHCEOF are higher than those of the ART and ARTHC. The emulation results validate that the reconstructed accuracy can be improved by introducing EOF to ARTHC.

3.2. IED Reconstruction during Geomagnetic Quiet Day

To examine the applicability of the ARTHCEOF method, the measured data of 183 GNSS stations shown in Figure 2 are used to compute the input STEC data in the European area. The geomagnetic activity is calm on 24 August 2018, so it is chosen to examine the ARTHCEOF. Carrier phase smoothing pseudorange is introduced to compute the input STEC data. The satellite DCB is corrected through the products released by the Chinese Academy of Science, and the computed formula of the receiver’s DCB is as follows:
D C B r = i = 1 m S T E C i M F ( z ) V T E C C O D E / m
where VTECCODE is the ionospheric VTEC at the puncture point obtained by bilinear interpolation; M F ( z ) is the projection function, which can be represented as:
M F ( z ) = cos ( sin 1 ( R R + H sin α z ) )
where z is the satellite elevation angle; H is the thin-layer height; R is Earth radius; α is the scale factor. To update the accuracy of input STEC measurements, the STEC values of the plasmasphere are removed using the IRI-plas model in this work.
Figure 7 illustrates three-dimensional IED variations in the European area on 24 August 2018. The tomographic results of the ARTHCEOF reveal that European IED values gradually increase as time goes on between 2:00 and 10:00UT, and then the IED begins to decrease. However, an ionospheric anomaly phenomenon occurs at 18:00UT, and the IED values reach the maximum of the day, the phenomenon is consistent with the release’s products of CODE. Figure 7 also reveals that the IED values decrease with an increasing latitude. The IED distributions in the east are generally greater than those in the west.
To confirm the dependability of the inversion results in Figure 7, the vertical profiles of the CIT are compared with those obtained from the PQ052 ionsonde station, which is located at the 50°N and 14.6°E. The compared images are shown in Figure 8. The comparisons show that the vertical profiles of ARTHCEOF coincide with those given by ionsonde data. It demonstrates that the tomographic results are dependable to some extent. However, the above comparisons are made in a fixed geographical location, and they cannot reflect the dependability of all inversion results. To comprehensively evaluate the dependability of ARTHCEOF, VTEC values are computed using the IED inversion results. By means of the computed VTEC values and the VTEC products of CODE, the differences can be obtained. Figure 9 illustrates the differential VTEC images in six time periods. It shows that the maximum difference is 2.5 TECU, and the minimum difference is −1.5 TECU. The statistics prove that the inverse VTEC values of ARTHCEOF are consistent with those released by CODE. Statistical results of differential VTEC thoroughly validate the dependability of ARTHCEOF.

3.3. IED Reconstruction during Geomagnetic Storm Day

Due to the influence of the solar coronal mass ejection, a strong geomagnetic storm erupted on 26 August 2018. Figure 10 shows that the Dst index peak reaches −174nT, and the Kp index reaches 8. So, the GNSS data of this day are selected to further examine the characteristics of the ARTHCEOF.
Figure 11 illustrates the two-dimensional images at 14:00UT on 25–28 August 2018. The fixed longitudinal chain is 14°E. Figure 11 displays that the IED values decrease at 14:00UT on 26 August 2018. This suggests that a negative phase storm occurs, and then the electron densities begin to augment in the recovery phase on 27 August 2018. However, the IED cannot recover to the level before the geomagnetic storm occurs. Then, the IED values return to the level of the initial phase on 28 August 2018.
Figure 12 shows the ionospheric daily variations on 26 August 2018. It reveals that the IED qualities augment between 2:00UT and 10:00UT, and then the IED qualities decrease as time goes on. However, the ionospheric variation exhibits abnormal phenomena at 18:00UT. In the meanwhile, it follows that the peak height is 350 km at 2:00UT, and then the peak heights decrease to 250 km between 6:00UT and 14:00UT. Subsequently, the peak height increases to 280 km at 18:00UT. At 22:00UT, the peak height recovered to 350 km.
Ionospheric VTEC values are inverted using the tomographic results of ARTHCEOF, and then differential VTECs are calculated by means of the computed VTEC values and VTEC products of CODE are calculated. Figure 13 shows the differential VTEC values of the same time periods in Figure 12. From Figure 13, it follows that the maximum absolute IEDD is less than 2 TECU. This validates the accuracy of the ARTHCEOF methods is comparable to the accuracy of CODE. The results further verify that the tomographic results of the ARTHCEOF are reliable.

4. Conclusions

The ARTHCEOF is proposed to reconstruct three-dimensional ionospheric images by using GNSS measurements in the European region. The new algorithm employs EOF to describe the ionospheric vertical structure instead of vertical constraints. It copes with the difficulty of constructing the vertical constraint matrix of the CIT system. The ART, the ARTHC, and the ARTHCEOF are tested using the designed simulation scheme. The numerical scheme demonstrates that the inverted efficiency and accuracy of ARTHCEOF are augmented. Three-dimensional IED images are successfully reconstructed. The CIT results of the ARTHCEOF show the ionospheric daily variation anomaly phenomenon. Finally, an ionospheric storm is used to further test the new algorithm. Compared with the CIT results on 25 August 2018, the ionospheric negative phase storm arises on 26 August 2018. In the meantime, the VTEC are inverted by using the CIT results during the geomagnetic quiet and disturbed days. The differential VTEC results reveal that the maximum values are smaller than 3TECU and 2TECU during geomagnetic quiet and disturbed days, respectively. The differential VTEC results validate the reliability of the ARTHCEOF method.
Although the ARTHCEOF can successfully capture three-dimensional ionospheric structure, the new algorithm has not been tested in other ionospheric events such as seismicity, solar flare and traveling ionospheric disturbance. In future, the ARTHCEOF will be extended to these fields.

Author Contributions

Methodology, D.W.; software, Y.T. and K.X.; validation, Y.T.; formal analysis, D.W.; writing—original draft preparation, D.W.; writing—review and editing, K.X.; funding acquisition, D.W. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the National Natural Science Foundation of China, grant number 42070430, the Natural Science Foundation of Guangdong Province of China, grant number 2022A1515011039. Guangzhou City School Joint Project, grant number 202201020136.

Data Availability Statement

The readers can obtain the data from the corresponding author.

Acknowledgments

Thanks World Data Center and EUREF Permanent GNSS Network for providing ionsonde data and the GNSS data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Su, K.; Jin, S.; Hoque, M.M. Evaluation of Ionospheric Delay Effects on Multi-GNSS Positioning Performance. Remote Sens. 2019, 11, 171. [Google Scholar] [CrossRef] [Green Version]
  2. Wen, D.B.; Xie, K.Y.; Tang, Y.H.; Mei, D.K.; Chen, X.; Chen, H.Q. A New Algorithm for Ill-posed Problem of GNSS-based Ionospheric Tomography. Remote Sens. 2023, 15, 1930. [Google Scholar] [CrossRef]
  3. Chen, C.H.; Saito, A.; Lin, C.H.; Yamamoto, M.; Suzuki, S.; Seemala, G.K. Medium-scale traveling ionospheric disturbances by three-dimensional ionospheric GPS tomography. Earth Planets Space 2016, 68, 32. [Google Scholar] [CrossRef] [Green Version]
  4. Kong, J.; Li, F.; Yao, Y.; Wang, Z.; Peng, W.; Zhang, Q. Reconstruction of 2D/3D ionospheric disturbances in high-latitude and Arctic regions during a geomagnetic storm using GNSS carrier TEC: A case study of the 2015 great storm. J. Geod. 2019, 93, 1529–1541. [Google Scholar] [CrossRef]
  5. Banville, S.; Collins, P.; Zhang, W.; Langley, R.B. Global and Regional Ionospheric Corrections for Faster PPP Convergence. Navigation 2014, 61, 115–124. [Google Scholar] [CrossRef]
  6. Belehaki, A.; Tsagouri, I.; Altadill, D.; Blanch, E.; Borries, C.; Buresova, D.; Chum, J.; Galkin, I.; Juan, J.M.; Segarra, A.; et al. An overview of methodologies for real-time detection, characterisation and tracking of traveling ionospheric disturbances developed in the TechTIDE project. J. Space Weather Spac. 2020, 10, 42. [Google Scholar] [CrossRef]
  7. Ciraolo, L.; Azpilicueta, F.; Brunini, C.; Meza, A.; Radicella, S.M. Calibration errors on experimental slant total electron content (TEC) determined with GPS. J. Geod. 2007, 81, 111–120. [Google Scholar] [CrossRef]
  8. Austen, J.R.; Franke, S.J.; Liu, C.H. Ionospheric imagingusing computerized tomography. Radio Sci. 1988, 23, 299–307. [Google Scholar] [CrossRef]
  9. Raymund, T.D.; Austen, J.R.; Franke, S.J.; Liu, C.H.; Klobuchar, J.A.; Stalker, J. Application of computerized tomography to the investigation of ionospheric structures. Radio Sci. 1990, 25, 771–789. [Google Scholar] [CrossRef]
  10. Fremouw, E.J.; Secan, J.A.; Howe, B.M. Application of stochastic theory to ionospheric tomogrpahy. Radio Sci. 1992, 27, 721–732. [Google Scholar] [CrossRef]
  11. Pryse, S.E.; Kersley, L.; Rice, D.L.; Russell, C.D.; Walker, I.K. Tomographic imaging of the ionospheric mid-latitude trough. Ann. Geophys. 1993, 11, 144–149. [Google Scholar]
  12. Kunitsyn, V.E.; Andreeva, E.S.; Razinkov, O.G. Possibilities of the near-space environment radio tomography. Radio Sci. 1997, 32, 1953–1963. [Google Scholar] [CrossRef]
  13. Giulio, R.; Alejandro, F.; Antonio, R. GPS tomography of the ionospheric electron content with a correlation functional. IEEE Trans. Geosci. Remote Sens. 1998, 36, 143–153. [Google Scholar]
  14. Fridman, S.V.; Nickisch, L.J. Generalization of ionospheric tomography on diverse data sources: Reconstruction of the three-dimensional ionosphere from simultaneous vertical ionograms, backscatter ionograms, and total electron content data. Radio Sci. 2001, 36, 1129–1139. [Google Scholar] [CrossRef]
  15. Ma, X.F. Three-dimensional ionospheric tomography using observation data of GPS ground receivers and ionosonde by neural network. J. Geophys. Res. 2005, 110, A05308. [Google Scholar] [CrossRef]
  16. Vierinen, J.; Norberg, J.; Lehtinen, M.S.; Amm, O.; Roininen, L.; Väänänen, A.; Erickson, P.J.; McKay-Bukowski, D. Beacon satellite receiver for ionospheric tomography. Radio Sci. 2014, 49, 1141–1152. [Google Scholar] [CrossRef]
  17. Wen, D.; Tang, Y.; Chen, X.; Zou, Y. A Double-Adaptive Adjustment Algorithm for Ionospheric Tomography. Remote Sens. 2023, 15, 2307. [Google Scholar] [CrossRef]
  18. Zheng, D.; Zheng, H.; Wang, Y.; Nie, W.; Li, C.; Ao, M.; Hu, W.; Zhou, W. Variable pixel size ionospheric tomography. Adv. Space Res. 2017, 59, 2969–2986. [Google Scholar] [CrossRef]
  19. Saito, S.; Suzuki, S.; Yamamoto, M.; Saito, A.; Chen, C.H. Real-time ionosphere monitoring by three-dimensional tomography over Japan. Navig. J. Inst. Navig. 2017, 64, 504. [Google Scholar] [CrossRef] [Green Version]
  20. Bhuyan, K.; Singh, S.B.; Bhuyan, P.K. Tomographic reconstruction of the ionosphere using generalized singular value decomposition. Curr. Sci. 2002, 83, 1117–1120. [Google Scholar]
  21. Yao, Y.; Tang, J.; Kong, J.; Zhang, L.; Zhang, S. Application of hybrid regularization method for tomographic reconstruction of midlatitude ionospheric electron density. Adv. Space Res. 2013, 52, 2215–2225. [Google Scholar] [CrossRef]
  22. Seemala, G.K.; Yamamoto, M.; Saito, A.; Chen, C.H. Three-dimensional GPS ionospheric tomography over Japan using constrained least squares. J. Geophys. Res. 2014, 119, 3044–3052. [Google Scholar] [CrossRef] [Green Version]
  23. Zhao, H.S.; Yang, L.; Zhou, Y.L.; Ming, D. A AMART Algorithm Applied to Ionospheric Electron Reconstruction. Acta Geod. Cartogr. Sin. 2018, 47, 57–63. [Google Scholar]
  24. Gerzen, T.; Minkwitz, D. Simultaneous multiplicative column-normalized method (SMART) for 3-D ionosphere tomography in comparison to other algebraic methods. Ann. Geophys. 2016, 34, 97–115. [Google Scholar] [CrossRef] [Green Version]
  25. Razin, M.R.G.; Voosoghi, B. Regional application of multi-layer artifificial neural networks in 3-D ionosphere tomography. Adv. Space Res. 2016, 58, 339–348. [Google Scholar] [CrossRef]
  26. Jin, S.; Li, D. 3-D ionospheric tomography from dense GNSS observations based on an improved two-step iterative algorithm. Adv. Space Res. 2018, 62, 809–820. [Google Scholar] [CrossRef]
  27. Hobiger, T.; Kondo, T.; Koyama, Y. Constrained simultaneous algebraic reconstruction technique (C-SART)—A new and simple algorithm applied to ionospheric tomography. Earth Planets Space 2008, 60, 727–735. [Google Scholar] [CrossRef] [Green Version]
  28. Wen, D.B.; Liu, S.J.; Tang, P.Y. Tomographic reconstruction of ionospheric electron density based on constrained algebraic reconstruction technique. GPS Solut. 2010, 14, 375–380. [Google Scholar] [CrossRef]
  29. Li, H.; Yuan, Y.B.; Yan, W.; Li, Z.S. A constrained ionospheric tomography algorithm with smoothing method. Geomat. Inf. Sci. Wuhan Univ. 2013, 38, 412–415. [Google Scholar]
  30. Nesterov, I.A.; Kunitsyn, V.E. GNSS radio tomography of the ionosphere: The problem with essentially incomplete data. Adv. Space Res. 2011, 47, 1789–1803. [Google Scholar] [CrossRef]
Figure 1. Diagrammatic sketch of the first three orders of EOF describing ionospheric information in altitudinal direction.
Figure 1. Diagrammatic sketch of the first three orders of EOF describing ionospheric information in altitudinal direction.
Remotesensing 15 03124 g001
Figure 2. The distribution of the 183 European GNSS stations.
Figure 2. The distribution of the 183 European GNSS stations.
Remotesensing 15 03124 g002
Figure 3. Images of IED background information at 14:00UT at 14°E and 305 km. (a) IED image along 14°E; (b) IED image at 350 km. The IED unit is 1011 el/m3.
Figure 3. Images of IED background information at 14:00UT at 14°E and 305 km. (a) IED image along 14°E; (b) IED image at 350 km. The IED unit is 1011 el/m3.
Remotesensing 15 03124 g003
Figure 4. IED inversion results of three CIT algorithms and IED true values along 14°E. (a) ART; (b) ARTHC; (c) ARTHCEOF; (d) IRI 2016 model. The IED unit is 1011 el/m3.
Figure 4. IED inversion results of three CIT algorithms and IED true values along 14°E. (a) ART; (b) ARTHC; (c) ARTHCEOF; (d) IRI 2016 model. The IED unit is 1011 el/m3.
Remotesensing 15 03124 g004
Figure 5. IED inversion results of three CIT algorithms and IED true values at 350 km. (a) ART; (b) ARTHC; (c) ARTHCEOF; (d) IRI 2016. The IED unit is 1011 el/m3.
Figure 5. IED inversion results of three CIT algorithms and IED true values at 350 km. (a) ART; (b) ARTHC; (c) ARTHCEOF; (d) IRI 2016. The IED unit is 1011 el/m3.
Remotesensing 15 03124 g005
Figure 6. Histogram of the reconstructed error statistics. The unit of differential IED is 1011 el/m3.
Figure 6. Histogram of the reconstructed error statistics. The unit of differential IED is 1011 el/m3.
Remotesensing 15 03124 g006
Figure 7. Three-dimensional IED distribution images during geomagnetic quiet day. The IED unit is 1011 el/m3.
Figure 7. Three-dimensional IED distribution images during geomagnetic quiet day. The IED unit is 1011 el/m3.
Remotesensing 15 03124 g007
Figure 8. Comparisons of ionospheric vertical profiles with those of ionosonde data.
Figure 8. Comparisons of ionospheric vertical profiles with those of ionosonde data.
Remotesensing 15 03124 g008
Figure 9. Differential VTEC images between the inverted VTEC and the products of CODE during geomagnetic quiet day. The unit of TEC is TECU.
Figure 9. Differential VTEC images between the inverted VTEC and the products of CODE during geomagnetic quiet day. The unit of TEC is TECU.
Remotesensing 15 03124 g009
Figure 10. Variations of the Dst and Kp indexes during 24–28 August 2018.
Figure 10. Variations of the Dst and Kp indexes during 24–28 August 2018.
Remotesensing 15 03124 g010
Figure 11. Two-dimensional images along 14°E at 14:00 on 25–28 August 2018.
Figure 11. Two-dimensional images along 14°E at 14:00 on 25–28 August 2018.
Remotesensing 15 03124 g011
Figure 12. IED daily variations along 14°E on 26 August 2018. The IED unit is 1011 el/m3.
Figure 12. IED daily variations along 14°E on 26 August 2018. The IED unit is 1011 el/m3.
Remotesensing 15 03124 g012
Figure 13. Differential VTEC images between the inverted VTEC and the products of CODE during geomagnetic storm day. The unit of TEC is TECU.
Figure 13. Differential VTEC images between the inverted VTEC and the products of CODE during geomagnetic storm day. The unit of TEC is TECU.
Remotesensing 15 03124 g013
Table 1. Census of error indicators and iterative efficiency of three CIT algorithms. The IED unit is 1011 el/m3.
Table 1. Census of error indicators and iterative efficiency of three CIT algorithms. The IED unit is 1011 el/m3.
MethodsARTARTHCARTHCEOF
Maximum absolute IEDD1.481.280.56
MAE0.230.110.06
RMSE0.330.180.14
Iterative round numbers231610
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wen, D.; Tang, Y.; Xie, K. A Novel Method of Ionospheric Inversion Based on Horizontal Constraint and Empirical Orthogonal Function. Remote Sens. 2023, 15, 3124. https://doi.org/10.3390/rs15123124

AMA Style

Wen D, Tang Y, Xie K. A Novel Method of Ionospheric Inversion Based on Horizontal Constraint and Empirical Orthogonal Function. Remote Sensing. 2023; 15(12):3124. https://doi.org/10.3390/rs15123124

Chicago/Turabian Style

Wen, Debao, Yinghao Tang, and Kangyou Xie. 2023. "A Novel Method of Ionospheric Inversion Based on Horizontal Constraint and Empirical Orthogonal Function" Remote Sensing 15, no. 12: 3124. https://doi.org/10.3390/rs15123124

APA Style

Wen, D., Tang, Y., & Xie, K. (2023). A Novel Method of Ionospheric Inversion Based on Horizontal Constraint and Empirical Orthogonal Function. Remote Sensing, 15(12), 3124. https://doi.org/10.3390/rs15123124

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop