A MARS-based method for estimating regional 2-D ionospheric VTEC and receiver differential code bias

The geometry-free linear combination of dual-frequency GNSS reference station ground observations are currently used to build the Vertical Total Electron Content (VTEC) model of the ionosphere. As it is known, besides ionospheric delays, there are differential code bias (DCB) of satellite (SDCB) and receiver (RDCB) in the geometry-free observation equation. The SDCB can be obtained using the International GNSS Service (IGS) analysis centers, but the RDCB for regional and local network receivers are not provided. Therefore, estimating the RDCB and VTEC model accurately and simultaneously is a critical factor investigated by researchers. This study uses Multivariate Adaptive Regression Splines (MARS) to estimate the VTEC approximate model and then substitutes this model in the observation equation to form the normal equation. The least squares method is used to solve the RDCB and VTEC model together. The research findings show that this method has good modeling effectiveness and the estimated RDCB has good reliability. The estimated VTEC model applied to GPS single-frequency precise point positioning has better positioning accuracy in comparison to the IGS global ionosphere map (GIM).     

The spatio-spectral localization approach to modeling VTEC over the western part of the USA using GPS observations

The precise ionosphere modeling is crucial and remains a challenge for GPS positioning and navigation, as well as many other Earth observation systems. In this research, two approaches have been proposed to model the vertical total electron content (VTEC) of the ionosphere using spherical slepian function. For the two-dimensional case, VTEC has been modeled in a sun-fixed reference frame and the three-dimensional approach based on the system of the three-dimensional base functions has been defined as the tensor product of the spherical slepian function for the longitude and latitude in an Earth-fixed reference frame, and the polynomial B-spline function for time. Rather than the spherical harmonics, the spherical slepian functions can be employed to produce the locally and globally orthogonal bases to optimally represent the data in any arbitrary region up to a given degree. The spherical slepian functions have been applied to the real data obtained from the ground-based GPS observation across the western part of the USA.     

Using GPS-TEC data to calibrate VTEC computed with the IRI model over Nigeria

This paper presents the development of a Total Electron Content (TEC) map for the Nigerian ionosphere. In this work, TEC measurements obtained from the AFRL-SCINDA GPS (Air Force Research Laboratory-Scintillation Network Decision Aid, Global Positioning System) equipment installed at Nsukka (6.87°N, 7.38°E) are used to adapt the International Reference Ionosphere (IRI) model for the Nigerian Ionosphere. The map is being developed as a computer program (implemented in the MATLAB programming language) that shows spatial and temporal representations of TEC for the Nigerian ionosphere. The method is aimed at showing how the IRI model can be used to estimate VTEC over wide areas by incorporating GPS measurements. This method is validated by using GPS VTEC data collected from a station in Ilorin (8.50°N, 4.55°E).     

Solar activity variations of ionosonde measurements and modeling results

The time series of hourly electron density profiles N(h) obtained at several mid-latitude stations in Europe have been used to obtain N(h) profiles on a monthly basis and to extract both the expected bottomside parameters and a proxy of the ionospheric variability as functions of time and height. With these data we present advances on a “Local Model” technique for the parameters B0 and B1, its applicability to other ionospheric stations, to other bottomside ionospheric parameters, and to modeling the time/height variability of the profile. The Local Model (LM) is an empirical model based on the experimental results of the solar activity dependence of the daily and seasonal behavior of the above parameters. The LM improves the IRI-2001 prediction of the B0 and B1 by factor of two at mid-latitudes. Moreover, the LM can be used to simulate other ionospheric parameters and to build mean N(h) profiles and the deviations from them. The modeling of both the average N(h) profiles and their deviations is an useful tool for ionospheric model users who want to know both the expected patterns and their deviations.     

Regional 4-D modeling of the ionospheric electron density from satellite data and IRI

Accurate knowledge of the electron density is the key point in correcting ionospheric delays of electromagnetic measurements and in studying ionosphere physics. During the last decade Global Navigation Satellite Systems (GNSS) have become a promising tool for monitoring ionospheric parameters such as the total electron content (TEC). In this contribution we present a four-dimensional (4-D) model of the electron density consisting of a given reference part, i.e., the International Reference Ionosphere (IRI), and an unknown correction term expanded in terms of multi-dimensional base functions. The corresponding series coefficients are calculable from the satellite measurements by applying parameter estimation procedures. Since satellite data are usually sampled between GPS satellites and ground stations, finer structures of the electron density are modelable just in regions with a sufficient number of ground stations. The proposed method is applied to simulated geometry-free GPS phase measurements. The procedure can be used, for example, to study the equatorial anomaly.     

GPS single-frequency orbit determination for low Earth orbiting satellites

The determination of high-precision orbits for Low Earth Orbiting (LEO) satellites (e.g., CHAMP, GRACE, MetOp-A) is based on dual-frequency tracking data from on-board GPS (Global Positioning System) receivers. The two frequencies allow it to eliminate the first order ionosphere effects. Data screening and precise orbit determination (POD) procedures are optimized under the assumption of the availability of two frequencies.     

Analysis of electron content variations over Japan during solar minimum: Observations and modeling

Comparative analysis of GPS TEC data and FORMOSAT-3/COSMIC radio occultation measurements was carried out for Japan region during period of the extremely prolonged solar minimum of cycle 23/24. COSMIC data for different seasons corresponded to equinox and solstices of the years 2007–2009 were analyzed. All selected electron density profiles were integrated up to the height of 700 km (altitude of COSMIC satellites), the monthly median estimates of Ionospheric Electron Content (IEC) were retrieved with use of spherical harmonics expansion. Monthly medians of TEC values were calculated from diurnal variations of GPS TEC estimates during considered month. Joint analysis of GPS TEC and COSMIC data allows us to extract and estimate electron content corresponded to the ionosphere (its bottom and topside parts) and the plasmasphere (h > 700 km) for different seasons of 2007–2009. Percentage contribution of ECpl to GPS TEC indicates the clear dependence from the time and varies from a minimum of about 25–50% during day-time to the value of 50–75% at night-time. Contribution of both bottom-side and topside IEC has minimal values during winter season in compare with summer season (for both day- and night-time). On average bottom-side IEC contributes about 5–10% of GPS TEC during night and about 20–27% during day-time. Topside IEC contributes about 15–20% of GPS TEC during night and about 35–40% during day-time. The obtained results were compared with TEC, IEC and ECpl estimates retrieved by Standard Plasmasphere–Ionosphere Model that has the plasmasphere extension up to 20,000 km (GPS orbit).     

South American regional ionospheric maps computed by GESA: A pilot service in the framework of SIRGAS

SIRGAS (Geocentric Reference Frame for the Americas) is an international enterprise of the geodetic community that aims to realize the Terrestrial Reference Frame in the America’s countries. In order to fulfill this commitment, SIRGAS manages a network of continuously operational GNSS receivers totalling around one hundred sites in the Caribbean, Central, and South American region. Although the network was not planed for ionospheric studies, its potential to be used for such a purpose was recently recognized and SIRGAS started a pilot experiment devoted to establish a regular service for computing and releasing regional vertical TEC (vTEC) maps based on GNSS data. Since July, 2005, the GESA (Geodesia Espacial y Aeronomía) laboratory belonging to the Facultad de Ciencias Astronómicas y Geofísicas of the Universidad Nacional de La Plata computes hourly maps of vertical Total Electron Content (vTEC) in the framework of the SIRGAS pilot experiment. These maps exploit all the GNSS data available in the South American region and are computed with the LPIM (La Plata Ionospheric Model). LPIM implements a de-biasing procedure that improves data calibration in relation to other procedures commonly used for such purposes. After calibration, slant TEC measurements are converted to vertical and mapped using local-time and modip latitude. The use of modip latitude smoothed the spatial variability of vTEC, especially in the South American low latitude region and hence allows for a better vTEC interpolation. This contribution summarizes the results obtained by GESA in the framework of the SIRGAS pilot experiment.     

Negative VTEC with spherical harmonic expansion model under different solar activity conditions

The Earth’s ionosphere can be described by a spherical harmonic (SH) expansion up to a specific degree. However, there exist negative vertical total electron content (VTEC) values in the global ionosphere map (GIM) with the SH expansion model. In this contribution, we specifically investigated the negative VTEC values that are induced by the SH expansion model and validated the performance of the inequality-constrained least squares (ICLS) method in eliminating the negative VTEC values. The GPS data from 2004 to 2017 was selected to cover one solar cycle and the experiments under different solar activity conditions were analyzed. The results in our work show that the occurrence of the negative VTEC values is attributed to the deficiency of the SH expansion model when the VTEC itself is small instead of the unevenly distribution of the GNSS stations. The negative VTEC values appear periodically in the temporal domain, showing apparently one year and half year periods. During one year, two peaks in June and December can be observed in the time series of the negative VTEC values. The number of negative VTEC values in June is obvious larger than that in December. During one solar cycle, the number of negative VTEC values under quiet solar activity condition is obvious larger than that under strong solar activity condition. In the spatial domain, the appearance of the negative VTEC values is strongly related with the movement of the subsolar point. In the latitude of the subsolar point has the largest magnitude, the negative values will appear on the opposite hemisphere and the further from the subsolar point the more negative values. The maximum number of the negative VTEC values in the southern hemisphere appears in June, while the peak value in the northern hemisphere appears in December. The maximum number of negative VTEC values in the southern hemisphere is generally larger than that in the northern hemisphere. In addition, the negative VTEC values are distributed both at middle latitude and high latitude in the southern hemisphere, while they are mainly distributed at high latitude in the northern hemisphere. When the ICLS method is used, the negative VTEC values can be eliminated efficiently and it has nearly no influence on the positive VTEC values. The ICLS method can also improve the receiver’s differential code bias (DCB) and significantly decrease the unreasonable negative slant TEC (STEC) values along the lines of sight. Using the final GIM product of the Jet Propulsion Laboratory (JPLG) as a reference, the root mean square (RMS) of the ICLS solution shows maximum 25%, 20% and 45% improvement relative to the least squares (LS) solution at northern high latitude, southern middle latitude and southern high latitude, respectively.     

TEC derived from some GPS stations in Nigeria and comparison with the IRI and NeQuick models

Total electron content (TEC) measured simultaneously using Global Positioning System (GPS) ionospheric monitors installed at some locations in Nigeria during the year 2011 (Rz = 55.7) was used to study the diurnal, seasonal, and annual TEC variations. The TEC exhibits daytime maximum, seasonal variation and semiannual variations. Measured TEC were compared with those predicted by the improved versions of the International Reference Ionosphere (IRI) and NeQuick models. The models followed the diurnal and seasonal variation patterns of the observed values of TEC. However, IRI model produced better estimates of TEC than NeQuick at all locations.     

基于BDS/GPS和经验模型的区域电离层建模

基于北斗卫星导航系统（BDS）和全球定位系统（GPS）实测电离层穿刺点（IPP）数据,结合国际参考电离层（IRI）经验模型历史数据,提出一种对区域二维电离层总电子含量（TEC）进行高精度建模的方法.针对缺乏穿刺点的区域内短时间电离层建模时精度较低且各时段穿刺点空间分布不同的问题,该方法使用IRI模型在建模区域内均匀添加虚拟穿刺点数据,并根据与实测穿刺点的距离,使用构造的权重计算公式赋予其动态权重值,通过加权最小二乘法进行球谐模型参数解算.与欧洲定轨中心（CODE）发布的全球电离层图（GIM）进行数据比对发现,相对于只使用BDS/GPS实测穿刺点数据的建模方法,利用本文建模方法计算获得的垂直总电子含量（VTEC）值对缺乏实测穿刺点的区域精度有明显的提升.     

Ionospheric correction using GPS Klobuchar coefficients with an empirical night-time delay model

We proposed an ionospheric correction approach called NKlob to mitigate the ionospheric delay errors. NKlob is a modification of the original GPS Ionospheric Correction Algorithm (ICA), which uses an empirical night-time model depending on the time, geomagnetic location and periodicities of the ionospheric behavior to replace the night-time constant delay in GPS ICA. Performance of NKlob was evaluated by the independent total electron contents (TECs) derived from Global Ionospheric Maps (GIMs) of the International GNSS Services (IGS) and Jason-2 altimetry satellite during 2013–2017. Compared to GIM TECs, NKlob corrects 51.5% of the ionospheric delay errors, which outperforms GPS ICA by 6.3%. Compared to Jason-2 TECs, NKlob mitigates the ionospheric errors by 58.1%, which is approximately 3.7% better than that of GPS ICA. NKlob shows significant improvement in low-latitude and equatorial regions with respect to GPS ICA, meanwhile exhibiting comparable performance at middle and high latitudes. Since NKlob only requires slight technical changes at the processing level of GPS receivers, we suppose that it can be easily implemented for better ionospheric delay corrections of real-time GPS single-frequency applications.     

Comparison of peak characteristics of F2 ionospheric layer over Tehran region at a low solar activity period with IRI-2001 and IRI-2007 models predictions

In the present work values of peak electron density (NmF2) and height of F2 ionospheric layer (hmF2) over Tehran region at a low solar activity period are compared with the predictions of the International Reference Ionosphere models (IRI-2001 and IRI-2007). Data measured by a digital ionosonde at the ionospheric station of the Institute of Geophysics, University of Tehran from July 2006 to June 2007 are used to perform the calculations. Formulations proposed by  and  are utilized to calculate the hmF2. The International Union of Radio Science (URSI) and International Radio Consultative Committee (CCIR) options are employed to run the IRI-2001 and IRI-2007 models. Results show that both IRI-2007 and IRI-2001 can successfully predict the NmF2 and hmF2 over Tehran region. In addition, the study shows that predictions of IRI-2007 model with CCIR coefficient has closer values to the observations. Furthermore, it is found that the monthly average of the percentage deviation between the IRI models predictions and the values of hmF2 and NmF2 parameters are less than 10% and 21%, respectively.     

Dual-frequency carrier smoothed code filtering with dynamical ionospheric delay modeling

In GNSS applications, carrier-smoothed-code is a widely used technique to combine code pseudo-range and carrier phase measurements. A dynamical ionospheric delay modeling method is proposed based on Kalman filter and least-squares theory. The level of the process noise is adaptively tuned along with the real-time KF state estimation, based on the online variance component estimation method. Meanwhile, the correlations of the time differenced carrier phase measurements are considered. This approach avoids overly optimistically evaluating the estimate and improves the transient accuracy of the estimates. A real GPS dataset is employed to check the performance of the proposed method under different conditions. The results show that the new algorithm can model the ionospheric delay variation well with different sampling intervals or even in ionospheric abnormal environment. The positioning accuracy can be confirmed, about 21%, 35% and 16% better are obtained in the N, E, and U direction than raw dataset.     

SHAKING: Adjusted spherical harmonics adding KrigING method for near real-time ionospheric modeling with multi-GNSS observations

Precise positioning based on Global Navigation Satellite System (GNSS) technique requires high accuracy ionospheric total electron content (TEC) correction models to account for the ionospheric path delay errors. We present an adjusted Spherical Harmonics Adding KrigING method (SHAKING) approach for regional ionospheric vertical TEC (VTEC) modeling in real time. In the proposed SHAKING method, the VTEC information over the sparse observation data area is extrapolated by the Adjusted Spherical Harmonic (ASH) function, and the boundary distortion in regional VTEC modeling is corrected by the stochastic VTEC estimated using Kriging interpolation. Using real-time GPS, GLONASS and BDS-2/3 data streams of the Crust Movement Observation Network of China (CMONOC), the SHAKING-based regional ionospheric VTEC maps are re-constructed over China and its boundary regions. Compared to GNSS VTECs derived from the independent stations, the quality of SHAKING solution improves by 13–31% and 6–33% with respect to the ASH-only solution during high and low geomagnetic periods, respectively. Compared to the inverse distance weighting (IDW) generated result, significant quality improved of SHAKING-based VTEC maps is also observed, especially over the edge areas with an improvement of 60–80%. Overall, the proposed SHAKING method exhibits notable advantage over the existing regional VTEC modeling techniques, which can be used for regional TEC modeling and associated high-precision positioning applications.     

Comparing IRI and a regional model with ionosonde measurements in Pakistan

We have used the technique of expansion in Empirical Orthogonal Functions (EOFs) to develop regional models of the critical frequencies of E and F₂ layers (f_oE, f_oF₂), peak height (h_mF₂), and semi-thickness of F₂ layer (Y_mF₂) over Pakistan. In the present study levels of solar activity specified by Smoothed Sunspot Number (R) from 10 to 200 are taken into account. The magnetic dip angle for the model ranges from 30° to 60°. We have compared the regional model and the International Reference Ionosphere (IRI) with measurements of three ionosondes in Pakistan. The model parameters f_oE and f_oF₂ are found overall comparable to the observed hourly median values during daytime at Karachi (geographic latitude = 24.95°N, longitude = 67.13°E, magnetic inclination = 37°), Multan (30.18°N, 71.48°E, 45°) and Islamabad (33.75°N, 73.13°E, 51.5°) during the years 1988, 1996 and 2000. For h_mF₂ the computed values by regional and IRI model for the year 1995 are found close to each other. However, for Y_mF₂the results are better during daytime as compared to nighttime.