A. Aragon-Angel

Comparative testing of four ionospheric models driven with GPS measurements

In the context of the European Space Agency/European Space Operations Centre funded Study A¢Â�Â�GNSS Contribution to Next Generation Global Ionospheric Monitoring,A¢Â�Â� four ionospheric models based on GNSS data (the Electron Density Assimilative Model, EDAM; the Ionosphere Monitoring Facility, IONMON v2; the Tomographic Ionosphere model, TOMION; and the Neustrelitz TEC Models, NTCM) have been run using a controlled set of input data. Each model output has been tested against differential slant TEC (dSTEC) truth data for high (May 2002) and low (December 2006) sunspot periods. Three of the models (EDAM, TOMION, and NTCM) produce dSTEC standard deviation results that are broadly consistent with each other and with standard deviation spreads of ~1 TECu for December 2006 and ~1.5 TECu for May 2002. The lowest reported standard deviation across all models and all stations was 0.99 TECu (EDAM, TLSE station for December 2006 night). However, the model with the best overall dSTEC performance was TOMION which has the lowest standard deviation in 28 out of 52 test cases (13 stations, two test periods, day and night). This is probably related to the interpolation techniques used in TOMION exploiting the spatial stationarity of vertical TEC error decorrelation.

Enhanced Precise Point Positioning for GNSS Users

This paper summarizes the main results obtained during the development of an Enhanced Precise Point Positioning (EPPP) Global Navigation Satellite Systems multifrequency user algorithm. The main innovations include the application of precise ionospheric corrections to facilitate the resolution of undifferenced carrier phase ambiguities, ambiguity validation, and integrity monitoring. The performance of the EPPP algorithm in terms of accuracy, convergence time, and integrity is demonstrated with actual GPS and simulated Galileo data. This can be achieved with very limited bandwidth requirements for EPPP users (less than 300 b/s for dual-frequency GPS data).

Wide Area RTK: A satellite navigation system based on precise real‐time ionospheric modelling

[1] The Wide Area Real Time Kinematic (WARTK) is an augmentation system concept for multi-frequency users based on precise real-time ionospheric modeling. It is able to provide a high accuracy and integrity GNSS positioning service over continental areas using the infrastructure of a network of permanent ground monitor stations, such as the European Geostationary Navigation Overlay Service (EGNOS) network of Ranging and Integrity Monitoring Stations (RIMS) in Europe. In this way, it allows an additional benefit to be obtained from these reference stations, that is, the network has the potential to support two independent systems: a satellite-based augmentation system, such as EGNOS, and a high-precision positioning service, based on WARTK. Indeed, thanks to the accuracy of the ionospheric corrections provided, WARTK users have available in real-time an extra constraint per satellite between the carrier phase ambiguities, which helps solve them quickly. Once such ambiguities have been solved, the GNSS user obtains navigation accurate to within 20 cm at the 95th percentile (about 10 cm RMS). Moreover, this precise positioning is achieved in a few minutes (with two frequency signals) or in a single epoch, after initial convergence of the tropospheric delay (with three frequency signals), even up to hundreds of kilometers away from the nearest reference station. While previous WARTK research has been devoted to implementing the concept and assessing its feasibility, considering in particular the accuracy achievable, the work reported in this paper focused on consolidating the results by analyzing a large and representative data set, and on deeper analysis of the integrity issue. It was carried out in the context of the Multi-constellation Regional System (MRS) project, within the European Space Agency GNSS Evolution Programme, with the aim of designing a high accuracy service for GPS and/or Galileo. Three months of actual data, from more than 25 permanent GPS stations in Europe, have been processed (some of them as a roving user), for high-, mid- and low-solar cycle conditions (in 2002, 2004 and 2006 respectively). In addition, several ionospheric storms occurred during the selected periods, with Dst values reaching up to � 150 nT. Results based on these data show that user domain integrity was maintained for baselines of up to 400 km. At the 95th percentile, the daily horizontal and vertical position errors were 20 and 30 cm, respectively, and the corresponding protection levels were about 1 and 2 m. The convergence time was around 5 minutes with actual GPS constellation data. The benefits of using a multi-constellation system were also studied, with simulated GPS and three-frequency Galileo data, showing that it is possible to reduce the convergence time to a few seconds.

Global prediction of the vertical total electron content of the ionosphere based on GPS data

[1] Although vertical total electron content (VTEC) forecasting is still an open subject of research, the use of predictions of the ionospheric state at a scale of several days is an area of increased interest. A global VTEC forecast product for two days ahead, which is based exclusively on actual Global Positioning System (GPS) data, has been developed in the frame of the International Global Navigation Satellite Systems (GNSS) Service (IGS) Ionospheric Working Group (IGS Iono‐WG). The UPC ionospheric VTEC prediction model is based on the Discrete Cosine Transform (DCT), which is widely used in image compression (for instance, in JPEG format). Additionally, a linear regression module is used to forecast the time evolution of each of the DCT coefficients. The use of the DCT coefficients is justified because they represent global features of the whole two‐ dimensional VTEC map/image. Also, one can therefore introduce prior information affecting the VTEC, for instance, smoothness or the distribution of relevant features in different directions. For this purpose, the use of a long time series of final/rapid UPC VTEC maps is required. Currently, the UPC Predicted product is being automatically generated in test mode and is made available through the main IGS server for public access. This product is also used to generate two days ahead preliminary combined IGS Predicted product. Finally, the results presented in this work suggest that the two days ahead UPC Predicted product could become an official IGS product in the near future.

Distribution and mitigation of higher‐order ionospheric effects on precise GNSS processing

Higher-order ionospheric effects (I2+) are one of the main limiting factors in very precise Global Navigation Satellite Systems (GNSS) processing, for applications where millimeter accuracy is demanded. This paper summarizes a comprehensive study of the I2+ effects in range and in GNSS precise products such as receiver position and clock, tropospheric delay, geocenter offset, and GNSS satellite position and clock. All the relevant higher-order contributions are considered: second and third orders, geometric bending, and slant total electron content (dSTEC) bending (i.e., the difference between the STEC for straight and bent paths). Using a realistic simulation with representative solar maximum conditions on GPS signals, both the effects and mitigation errors are analyzed. The usage of the combination of multifrequency L band observations has to be rejected due to its increased noise level. The results of the study show that the main two effects in range are the second-order ionospheric and dSTEC terms, with peak values up to 2 cm. Their combined impacts on the precise GNSS satellite products affects the satellite Z coordinates (up to +1 cm) and satellite clocks (more than ±20 ps). Other precise products are affected at the millimeter level. After correction the impact on all the precise GNSS products is reduced below 5 mm. We finally show that the I2+ impact on a Precise Point Positioning (PPP) user is lower than the current uncertainties of the PPP solutions, after applying consistently the precise products (satellite orbits and clocks) obtained under I2+ correction

A linear scale height Chapman model supported by GNSS occultation measurements

GNSS radio occultations allow the vertical sounding of the Earths atmosphere, in particular the ionosphere. The physical observables estimated with this technique permit to test theoretical models of the electron density such as, for example, the Chapman and the Vary-Chap models. The former is characterized by a constant scale height, whereas the latter considers a more general function of the scale height with respect to height. We propose to investigate the feasibility of the Vary-Chap model where the scale height varies linearly with respect to height. In order to test this hypothesis, the scale height data provided by radio occultations from a receiver on board a low Earth orbit (LEO) satellite, obtained by iterating with a local Chapman model at every point of the topside F2 layer provided by the GNSS satellite occultation, are fitted to height data by means of a linear least squares fit (LLS). Results, based on FORMOSAT-3/COSMIC GPS occultation data inverted by means of the Improved Abel transform inversion technique (which takes into account the horizontal electron content gradients), show that the scale height presents a more clear linear trend above the F2 layer peak height, hm, which is in good agreement with the expected linear temperature dependence. Moreover, the parameters of the linear fit obtained during four representative days for all seasons, depend significantly on local time and latitude, strongly suggesting that this approach can significantly contribute to build realistic models of the electron density directly derived from GNSS occultation data.

Ionospheric and plasmaspheric electron contents inferred from radio occultations and global ionospheric maps

We introduce a methodology to extract the separate contributions of the ionosphere and the plasmasphere to the vertical total electron content, without relying on a fixed altitude to perform that separation. The method combines two previously developed and tested techniques, namely, the retrieval of electron density profiles from radio occultations using an improved Abel inversion technique and a two-component model for the topside ionosphere plus protonosphere. Taking measurements of the total electron content from global ionospheric maps and radio occultations from the Constellation Observing System for Meteorology, Ionosphere, and Climate/FORMOSAT-3 constellation, the ionospheric and plasmaspheric electron contents are calculated for a sample of observations covering 2007, a period of low solar and geomagnetic activity. The results obtained are shown to be consistent with previous studies for the last solar minimum period and with model calculations, confirming the reversal of the winter anomaly, the hemispheric asymmetry of the semiannual anomaly, and the existence in the plasmasphere of an annual anomaly in the South American sector of longitudes. The analysis of the respective fractional contributions from the ionosphere and the plasmasphere to the total electron content shows quantitatively that during the night the plasmasphere makes the largest contribution, peaking just before sunrise and during winter. On the other hand, the fractional contribution from the ionosphere reaches a maximum value around noon, which is nearly independent of season and geomagnetic latitude.

Electron density extrapolation above F2 peak by the linear Vary‐Chap model supporting new Global Navigation Satellite Systems‐LEO occultation missions

The new radio-occultation (RO) instrument on-board the future EUMETSAT Polar System 2nd Generation (EPS-SG) satellites, flying at a height of 820 km, is primarily focusing on neutral atmospheric profiling. It will also provide an opportunity for RO ionospheric sounding, but only below impact heights of 500 km, in order to guarantee a full data gathering of the neutral part. This will leave a gap of 320 km, which impedes the application of the direct inversion techniques to retrieve the electron density profile. To contribute to overcome this challenge, we have looked for new ways (accurate and simple) of extrapolating the electron density (also applicable to other Low-Earth-Orbiting, LEO, missions like CHAMP): a new Vary-Chap Extrapolation Technique (VCET). VCET is based on the scale height behaviour, linearly dependent on the altitude above hmF2 . This allows the electron density profile extrapolation for impact heights above its peak height (this is the case for EPS-SG), up to the satellite orbital height. VCET has been assessed with more than 3700 complete electron density profiles obtained in 4 representative scenarios of FORMOSAT-3/COSMIC occultations, in solar maximum and minimum conditions, and geomagnetically disturbed conditions, by applying an updated Improved Abel Transform Inversion technique to dual-frequency GPS measurements. It is shown that VCET performs much better than other classical Chapman models, with 60% of occultations showing relative extrapolation errors below 20%, in contrast with conventional Chapman model extrapolation approaches with 10% or less of the profiles with relative error below 20%.

Toward a SIRGAS Service for Mapping the Ionosphere’s Electron Density Distribution

SIRGAS is responsible of the terrestrial reference frame of Latin America and the Caribbean. To fulfil this commitment it manages a continuously operational GNSS network with more than 200 receivers. Although that network was not planed for ionospheric studies, SIRGAS attempted to exploit it by establishing, in early 2008, a regular service for computing regional maps of the vertical Total Electron Content. This paper describes an effort for developing a new SIRGAS product, concretely, a 4-dimensional (space and time) representation of the free electron distribution in the ionosphere. The working methodology is based on the ingestion of dual-frequency GNSS observations into a global electron density model in order to update its parameters. Preliminary results are presented and their quality is assessed by comparing the electron density computed with the methodology here described and the one estimated from totally independent observations. A preliminary analysis reveals that the performance of the electron density model improves by a factor greater than 2 after data ingestion.

A method for scintillation characterization using geodetic receivers operating at 1 Hz

Ionospheric scintillation produces strong disruptive effects on global navigation satellite system (GNSS) signals, ranging from degrading performances to rendering these signals useless for accurate navigation. The current paper presents a novel approach to detect scintillation on the GNSS signals based on its effect on the ionospheric-free combination of carrier phases, i.e. the standard combination of measurements used in precise point positioning (PPP). The method is implemented using actual data, thereby having both its feasibility and its usefulness assessed at the same time. The results identify the main effects of scintillation, which consist of an increased level of noise in the ionospheric-free combination of measurements and the introduction of cycle-slips into the signals. Also discussed is how mis-detected cycle-slips contaminate the rate of change of the total electron content index (ROTI) values, which is especially important for low-latitude receivers. By considering the effect of single jumps in the individual frequencies, the proposed method is able to isolate, over the combined signal, the frequency experiencing the cycle-slip. Moreover, because of the use of the ionospheric-free combination, the method captures the diffractive nature of the scintillation phenomena that, in the end, is the relevant effect on PPP. Finally, a new scintillation index is introduced that is associated with the degradation of the performance in navigation.