GNSS-Based Non-Negative Absolute Ionosphere Total Electron Content, its Spatial Gradients, Time Derivatives and Differential Code Biases: Bounded-Variable Least-Squares and Taylor Series

Global navigation satellite systems (GNSS) allow estimating total electron content (TEC). However, it is still a problem to calculate absolute ionosphere parameters from GNSS data: negative TEC values could appear, and most of existing algorithms does not enable to estimate TEC spatial gradients and TEC time derivatives. We developed an algorithm to recover the absolute non-negative vertical and slant TEC, its derivatives and its gradients, as well as the GNSS equipment differential code biases (DCBs) by using the Taylor series expansion and bounded-variable least-squares. We termed this algorithm TuRBOTEC. Bounded-variable least-squares fitting ensures non-negative values of both slant TEC and vertical TEC. The second order Taylor series expansion could provide a relevant TEC spatial gradients and TEC time derivatives. The technique validation was performed by using independent experimental data over 2014 and the IRI-2012 and IRI-plas models. As a TEC source we used Madrigal maps, CODE (the Center for Orbit Determination in Europe) global ionosphere maps (GIM), the IONOLAB software, and the SEEMALA-TEC software developed by Dr. Seemala. For the Asian mid-latitudes TuRBOTEC results agree with the GIM and IONOLAB data (root-mean-square was < 3 TECU), but they disagree with the SEEMALA-TEC and Madrigal data (root-mean-square was >10 TECU). About 9% of vertical TECs from the TuRBOTEC estimates exceed (by more than 1 TECU) those from the same algorithm but without constraints. The analysis of TEC spatial gradients showed that as far as 10–15° on latitude, TEC estimation error exceeds 10 TECU. Longitudinal gradients produce smaller error for the same distance. Experimental GLObal Navigation Satellite System (GLONASS) DCB from TuRBOTEC and CODE peaked 15 TECU difference, while GPS DCB agrees. Slant TEC series indicate that the TuRBOTEC data for GLONASS are physically more plausible.

[1]  J. Gonzalez-Esparza,et al.  Features of TEC behaviour over the low-latitude North-American region during the period of medium solar activity , 2017 .

[2]  Jingnan Liu,et al.  The Barcelona ionospheric mapping function (BIMF) and its application to northern mid-latitudes , 2018, GPS Solutions.

[3]  Anthony J. Mannucci,et al.  Automated daily processing of more than 1000 ground‐based GPS receivers for studying intense ionospheric storms , 2005 .

[4]  Bodo W. Reinisch,et al.  Global Ionospheric Radio Observatory (GIRO) , 2011 .

[5]  Yu. V. Yasyukevich,et al.  Variability of GPS/GLONASS differential code biases , 2015 .

[6]  P. Stark.  Bounded-Variable Least-Squares: an Algorithm and Applications , 2008 .

[7]  Markus Rothacher,et al.  Mapping and predicting the ionosphere , 1998 .

[8]  I. V. Zhivetiev,et al.  Global electron content: a new conception to track solar activity , 2008 .

[9]  E. Afraimovich,et al.  Using GPS-GLONASS-GALILEO data and IRI modeling for ionospheric calibration of radio telescopes and radio interferometers , 2008 .

[10]  Artem M. Padokhin,et al.  Influence of GPS/GLONASS differential code biases on the determination accuracy of the absolute total electron content in the ionosphere , 2015, Geomagnetism and Aeronomy.

[11]  Guanyi Ma,et al.  Derivation of TEC and estimation of instrumental biases from GEONET in Japan , 2002 .

[12]  Murat Durmaz,et al.  Regional vertical total electron content (VTEC) modeling together with satellite and receiver differential code biases (DCBs) using semi-parametric multivariate adaptive regression B-splines (SP-BMARS) , 2015, Journal of Geodesy.

[13]  M. Klimenko,et al.  Investigating range error compensation in UHF radar through IRI-2007 real-time updating: Preliminary results , 2015 .

[14]  W. E. Gordon,et al.  Incoherent Scattering of Radio Waves by Free Electrons with Applications to Space Exploration by Radar , 1958, Proceedings of the IRE.

[15]  Irina Zakharenkova,et al.  Ionospheric response to the 2015 St. Patrick's Day storm: A global multi‐instrumental overview , 2015 .

[16]  Shuanggen Jin,et al.  GPS ionospheric tomography: A comparison with the IRI-2001 model over South Korea , 2007 .

[17]  M. Aquino,et al.  On the estimate and assessment of the ionospheric effects affecting low frequency radio astronomy measurements , 2011, 2011 XXXth URSI General Assembly and Scientific Symposium.

[18]  Thomas L. Gaussiran,et al.  Ionospheric Data Assimilation Three‐Dimensional (IDA3D): A global, multisensor, electron density specification algorithm , 2004 .

[19]  Qiang Zhang,et al.  BDS Satellites and Receivers DCB Resolution , 2014 .

[20]  Jin Wang,et al.  Mitigation of Short-Term Temporal Variations of Receiver Code Bias to Achieve Increased Success Rate of Ambiguity Resolution in PPP , 2020, Remote. Sens..

[21]  J. Bernard Minster,et al.  GPS detection of ionospheric perturbations following the January 17, 1994, Northridge Earthquake , 1995 .

[22]  Zishen Li,et al.  Two-step method for the determination of the differential code biases of COMPASS satellites , 2012, Journal of Geodesy.

[23]  Chris Rizos,et al.  The International GNSS Service in a changing landscape of Global Navigation Satellite Systems , 2009 .

[24]  Jaume Sanz,et al.  New approaches in global ionospheric determination using ground GPS data , 1999 .

[25]  E. Afraimovich,et al.  Global electron content during solar cycle 23 , 2008 .

[26]  Gabor E. Lanyi,et al.  A comparison of mapped and measured total ionospheric electron content using global positioning system and beacon satellite observations , 1988 .

[27]  Wu Chen,et al.  Revisit the calibration errors on experimental slant total electron content (TEC) determined with GPS , 2018, GPS Solutions.

[28]  D. Venkata Ratnam,et al.  Determination of NavIC differential code biases using GPS and NavIC observations , 2020 .

[29]  Artem M. Padokhin,et al.  Ionospheric radio tomography based on the GPS/GLONASS navigation systems , 2011 .

[30]  Shuanggen Jin,et al.  M_DCB: Matlab code for estimating GNSS satellite and receiver differential code biases , 2012, GPS Solutions.

[31]  V. Lebedev,et al.  Recording and control digital systems of the Irkutsk Incoherent Scatter Radar , 2009 .

[32]  Anthea J. Coster,et al.  Automated GPS processing for global total electron content data , 2006 .

[33]  Weidong Zhu,et al.  Investigation and Validation of the Time-Varying Characteristic for the GPS Differential Code Bias , 2019, Remote. Sens..

[34]  Anthony J. Mannucci,et al.  A global mapping technique for GPS‐derived ionospheric total electron content measurements , 1998 .

[35]  Yu. V. Yasyukevich,et al.  Ionospheric TEC estimation with the signals of various geostationary navigational satellites , 2015, GPS Solutions.

[36]  C. B. Erol,et al.  Regularized estimation of TEC from GPS data for certain midlatitude stations and comparison with the IRI model , 2007 .

[37]  J. Klobuchar.  Ionospheric Time-Delay Algorithm for Single-Frequency GPS Users , 1987, IEEE Transactions on Aerospace and Electronic Systems.

[38]  Dieter Bilitza,et al.  Towards ISO Standard Earth Ionosphere and Plasmasphere Model , 2012 .

[39]  Dorota A. Grejner-Brzezinska,et al.  Efficient GPS receiver DCB estimation for ionosphere modeling using satellite-receiver geometry changes , 2008 .

[40]  L. Gonz'alez,et al.  First joint observations of space weather events over Mexico , 2018, Annales Geophysicae.

[41]  Richard B. Langley,et al.  The nature of GPS differential receiver bias variability: An examination in the polar cap region , 2015 .

[42]  C. Mitchell,et al.  A three-dimensional time-dependent algorithm for ionospheric imaging using GPS , 2003 .

[43]  L. Rolland,et al.  Parameters of seismic source as deduced from 1 Hz ionospheric GPS data: Case study of the 2011 Tohoku‐oki event , 2013 .

[44]  Yulia S. Tumanova,et al.  Efficiency of updating the ionospheric models using total electron content at mid- and sub-auroral latitudes , 2019, GPS Solutions.

[45]  Michael S. Waterman,et al.  A Restricted Least Squares Problem , 1974 .

[46]  S. V. Voeykov,et al.  A review of GPS/GLONASS studies of the ionospheric response to natural and anthropogenic processes and phenomena , 2013 .

[47]  Torben Schüler,et al.  Single-frequency single-site VTEC retrieval using the NeQuick2 ray tracer for obliquity factor determination , 2013, GPS Solutions.

[48]  Cathryn N. Mitchell,et al.  Measurement of Ionospheric Total Electron Content Using Single‐Frequency Geostationary Satellite Observations , 2019, Radio Science.

[49]  V. E. Kunitsyn,et al.  GNSS radio tomography of the ionosphere: The problem with essentially incomplete data , 2011 .

[50]  S. Schaer.  Mapping and predicting the Earth's ionosphere using the Global Positioning System. , 1999 .

[51]  Vladimir Truhlik,et al.  The International Reference Ionosphere 2012 – a model of international collaboration , 2014 .

[52]  Bodo W. Reinisch,et al.  New Digisonde for research and monitoring applications , 2009 .

[53]  E. Sardón,et al.  Estimation of total electron content using GPS data: How stable are the differential satellite and receiver instrumental biases? , 1997 .

[54]  Maorong Ge,et al.  Eliminating negative VTEC in global ionosphere maps using inequality-constrained least squares , 2013 .

[55]  A. Garcia-Rigo,et al.  The IGS VTEC maps: a reliable source of ionospheric information since 1998 , 2009 .

[56]  Sang Jeong Lee,et al.  The influence of grounding on GPS receiver differential code biases , 2018, Advances in Space Research.

[57]  Geoffrey Blewitt,et al.  An Automatic Editing Algorithm for GPS data , 1990 .