Impact of Approximate Implementation manner in Klobuchar Model

This paper focuses on the approximations that John A. Klobuchar made in mid 70s in his famous algorithm of ionospheric correction model for single frequency GPS receiver. At that time Klobuchar used a system of fixed geomagnetic north pole coordinates which are not accurate nowadays according to the International Geomagnetic Reference Field and to the World Magnetic Model because the geomagnetic poles move slowly. In addition, Klobuchar had to do other trigonometry simplifications in his implementation to avoid sophisticated computations. In order to evaluate this approximate implementation in a single frequency GPS receiver, ionospheric time and range delay are estimated on the entire day of January 1st 2010, using a different implementation in MATLAB. The required GPS data is obtained from recorded RINEX files at UDMC near DAMASCUS, SYRIA. In this comparative study, we reformulated the standard equations of Klobuchar model and examined the influence of its approximations on the ionospheric range delay and found a non- negligible bias in order of ten centimeters, whereas the influence of the movement of the geomagnetic poles was in order of few centimeters.

Galileo Ionospheric Correction Algorithm Integration into the Open-Source GNSS Laboratory Tool Suite (gLAB)

Users of the global navigation satellite system (GNSS) operating with a single-frequency receiver must use an ionospheric correction algorithm (ICA) to account for the delay introduced on radio waves by the upper atmosphere. Galileo, the European GNSS, uses an ICA named NeQuick-G. In an effort to foster the adoption of NeQuick-G by final users, two implementations in C language have been recently made available to the public by the European Space Agency (ESA) and the Joint Research Centre (JRC) of the European Commission (EC), respectively. The aim of the present contribution is to compare the slant total electron content (STEC) predictions of the two aforementioned implementations of NeQuick-G. For this purpose, we have used actual multi-constellation and multi-frequency data for several hundreds of stations distributed worldwide belonging to the Multi GNSS Experiment (MGEX) network of the International GNSS Service (IGS). For each first day of the month during year 2019, the STECs of the two NeQuick-G versions were compared in terms of accuracy, consistency, availability, and execution time. Our study concludes that both implementations of NeQuick-G perform equivalently. Indeed, in over 99.998% of the 2125 million STECs computed, the output is exactly coincident. In contrast, 0.002% of the whole set of STECs for those rays are tangent to the Earth, the behavior of both implementations differs. We confirmed the discrepancy by processing radio-occultation actual measurements from a COSMIC-2 low Earth orbit satellite. We selected the JRC version of the Galileo ICA to be integrated into the GNSS LABoratory (gLAB) tool suite, because its open license and its processing speed (it is 13.88% faster than the ESA version). NeQuick-G outperforms the GPS ICA in STEC residuals up to 12.15 TECUs (percentile 96.23th) and in the 3D position errors, up to 5.76 m (percentile 99.18th) for code-pseudorange positioning.

Postseismic deformation following the 2016 Kumamoto earthquake detected by ALOS-2/PALSAR-2

I have been conducting a study of postseismic deformation following the 2016 Kumamoto earthquake using ALOS-2/PALSAR-2 acquired till 2018. I apply ionospheric correction to interferograms of ALOS-2/PALSAR-2. L-band SAR gives us high coherence enough to reveal surface deformation even in vegetated or mountainous area for pairs of images acquired more than 2 years. Postseismic deformation following the Kumamoto earthquake exceeds 10 cm during 2 years at some spots in and around Kumamoto city and Aso caldera. Westward motion of ~ 6 cm/year was dominant on the southeast side of the Hinagu fault, while westward shift was detected on both sides of the Futagawa fault. The area of latter deformation seems to have correlation with distribution of pyroclastic flow deposits. Significant uplift was found around the eastern Futagawa fault and on the southwestern frank of Aso caldera, whose rate reaches 4 cm/year. There are sharp changes across several coseismic surface ruptures such as Futagawa, Hinagu, and Idenokuchi faults. Rapid subsidence between Futagawa and Idenokuchi faults also found. It is confirmed that local subsidence continued along the Suizenji fault, which newly appeared during the mainshock in Kumamoto City. Subsidence with westward shift of up to 4 cm/year was also found in Aso caldera. Time constant of postseismic decay ranges from 1 month to 600 days at selected points, but that postseismic deformation during the first epochs or two is dominant at point in the Kumamoto Plain. This result suggests multiple source of deformation. Westward motion around the Hinagu fault may be explained with right lateral afterslip on the shallow part of this fault. Subsidence along the Suizenji fault can be attributed to normal faulting on dipping westward. Deformation around the Hinagu and Idenokuchi faults cannot be explained with right lateral afterslip of Futagawa fault, which requires other sources. Deformation in northern part of Aso caldera might be the result of right lateral afterslip on a possible buried fault.

NeQuick-G and Android Devices: A Compromise between Computational Burden and Accuracy

Ionospheric delay is one of the largest errors affecting Global Navigation Satellite System (GNSS) positioning in open-sky conditions, and different methods are currently available for mitigating ionospheric effects including dual-frequency measurements and corrections from augmentation systems. For single-frequency standalone receivers, the most widely used approach to correct ionospheric delays is to rely on a model. In this respect, Klobuchar and NeQuick-G Ionospheric Correction Algorithms (ICAs) are the approaches adopted by GPS and Galileo, respectively. While the latter outperforms the Klobuchar model, it requires a significantly higher computational load, which can limit its exploitation in some market segments such as smartphones. In order to foster adoption of the NeQuick-G model in this type of device, a smart application of NeQuick-G is proposed. The solution relies on the assumption that ionospheric delays are practically constant over short time intervals. Thus, the update rate of the ionospheric correction computation can be significantly reduced. This solution was implemented, tested, and evaluated using real data collected with a static smartphone in an ad hoc set-up. The impact of reducing the ionospheric correction update rate has been evaluated in terms of processing time, of ionospheric correction deviations and in the Ranging Error (RE) and position domains. The analysis shows that a significant reduction of the processing time can be obtained with negligible degradation of the navigation solution.

Postseismic Deformation Following the 2016 Kumamoto Earthquake Detected by ALOS-2/PALSAR-2

I have been conducting a study of postseismic deformation following the 2016 Kumamoto earthquake using ALOS-2/PALSAR-2 acquired till 2018. I apply ionospheric correction to interferograms of ALOS-2/PALSAR-2. L-band SAR gives us high coherence enough to reveal surface deformation even in vegetated or mountainous area for pairs of images acquired more than 2 years. Postseismic deformation following the Kumamoto earthquake exceeds 10 cm during two years at some spots in and around Kumamoto city and Aso caldera. Westward motion of ~6 cm/yr was dominant on the southeast side of the Hinagu fault, while westward shift was detected on both side of the Futagawa fault. The area of latter deformation seems to have correlation with distribution of pyroclastic flow deposits. Significant uplift was found around the eastern Futagawa fault and on the southwestern frank of Aso caldera, whose rate reaches 4 cm/yr. There are sharp changes across several coseismic surface ruptures such as Futagawa, Hinagu, and Idenokuchi faults. Rapid subsidence between Futagawa and Idenokuchi faults also found. It is confirmed that local subsidence continued along the Suizenji fault, which newly appeared during the mainshock in Kumamoto City. Subsidence with westward shift of up to 4 cm/yr was also found in Aso caldera.Time constant of postseismic decay ranges from 1 month to 600 days at selected points, but that postseismic deformation during the first epochs or two are dominant at point in the Kumamoto Plain. This result suggests multiple source of deformation. Westward motion around the Hinagu fault may be explained with right lateral afterslip on the shallow part of this fault. Subsidence along the Suizenji fault can be attributed to normal faulting on dipping westward. Deformation around the Hinagu and Idenokuchi faults cannot be explained with right-lateral afterslip of Futagawa fault, which requires other sources. Deformation in northern part of Aso caldera might be the result of right lateral afterslip on a possible buried fault.

Postseismic Deformation Following the 2016 Kumamoto Earthquake Detected by ALOS-2/PALSAR-2

I have been conducting a study of postseismic deformation following the 2016 Kumamoto earthquake using ALOS-2/PALSAR-2 acquired till 2018. I apply ionospheric correction to interferograms of ALOS-2/PALSAR-2. L-band SAR gives us high coherence enough to reveal surface deformation even in vegetated or mountainous area for pairs of images acquired more than 2 years. Postseismic deformation following the Kumamoto earthquake exceeds 10 cm during two years at some spots in and around Kumamoto city and Aso caldera. Westward motion of ~6 cm/yr was dominant on the southeast side of the Hinagu fault, while westward shift was detected on both side of the Futagawa fault. The area of latter deformation seems to have correlation with distribution of pyroclastic flow deposits. Significant uplift was found around the eastern Futagawa fault and on the southwestern frank of Aso caldera, whose rate reaches 4 cm/yr. There are sharp changes across several coseismic surface ruptures such as Futagawa, Hinagu, and Idenokuchi faults. Rapid subsidence between Futagawa and Idenokuchi faults also found. It is confirmed that local subsidence continued along the Suizenji fault, which newly appeared during the mainshock in Kumamoto City. Subsidence with westward shift of up to 4 cm/yr was also found in Aso caldera. Time constant of postseismic decay ranges from 1 month to 600 days at selected points, but that postseismic deformation during the first epochs or two are dominant at point in the Kumamoto Plain. This result suggests multiple source of deformation. Westward motion around the Hinagu fault may be explained with right lateral afterslip on the shallow part of this fault. Subsidence along the Suizenji fault can be attributed to normal faulting on dipping westward. Deformation around the Hinagu and Idenokuchi faults cannot be explained with right-lateral afterslip of Futagawa fault, which requires other sources. Deformation in northern part of Aso caldera might be the result of right lateral afterslip on a possible buried fault.