An evaluation of potential solar radio emission power threat on GPS and GLONASS performance

The L-band solar radio emission has recently been regarded as a potential threat to stable GPS and GLONASS performance. However, the threat has not been completely investigated or assessed so far. We evaluate in detail the occurrence of GPS/GLONASS signal tracking failures under the direct exposure of wideband solar radio emission. By means of theoretical analysis, we found that the solar radio emission power level of 1,000?sfu (solar flux units) or higher can cause GPS/GLONASS signal tracking failures especially at L2 frequency. In order to prove this evaluation, we investigated GPS/GLONASS signal tracking failures at L1 and L2 frequencies during power solar flares X6.5 (December 6, 2006) and X3.4 (December 13, 2006). Comparing these events with weaker solar flare X17.2 on October 28, 2003, we found that L2 signal tracking failures appeared when the solar radio emission power exceeds 1,000?sfu. Therefore, our theoretical and experimental results confirm the earlier results by other authors.     

Climatological study of ionospheric irregularities over the European mid-latitude sector with GPS

High-frequency variability of the ionosphere, or irregularities, constitutes the main threat for real-time precise positioning techniques based on Global Navigation Satellite Systems (GNSS) measurements. Indeed, during periods of enhanced ionospheric variability, GNSS users in the field—who cannot verify the integrity of their measurements—will experience positioning errors that can reach several decimeters, while the nominal accuracy of the technique is cm-level. In the frame of this paper, a climatological analysis of irregularities over the European mid-latitude region is presented. Based on a 10 years GPS dataset over Belgium, the work analyzes the occurrence rate (as a function of the solar cycle, season and local time) as well as the amplitude of ionospheric irregularities observed at a single GPS station. The study covers irregularities either due to space weather events (solar origin) or of terrestrial origin. If space weather irregularities are responsible for the largest effects in terms of ionospheric error, their occurrence rate highly depends on solar activity. Indeed, the occurrence rate of ionospheric irregularities is about 9 % during solar maximum, whereas it drops to about 0 % during medium or low solar activity periods. Medium-scale ionospheric disturbances (MSTIDs) occurring during daytime in autumn/winter are the most recurrent pattern of the time series, with yearly proportions slightly varying with the solar cycle and an amplitude of about 10 % of the TEC background. Another recurrent irregularity type, though less frequent than MSTIDs, is the noise-like variability in TEC observed during summer nighttime, under quiet geomagnetic conditions. These summer nighttime irregularities exhibit amplitudes ranging between 8 and 15 % of the TEC background.     

地基GNSS大气水汽探测遥感研究进展和展望

大气水汽是表征极端天气事件和气候变化的重要参数,准确监测与分析水汽含量对于精准预测各类灾害性天气事件与研究气候变化具有显著意义。作为新兴的大气水汽探测方法,GNSS大气水汽探测技术得到了广泛的关注与应用研究,随着多频多模GNSS系统的发展,全球服务能力的逐步完备和地面基础设施的不断加强,地基GNSS大气水汽探测遥感技术水平得到显著提升,为基于空间大数据揭示气候变化、极端天气过程提供了强有力的数据支撑和发展契机。本文首先系统阐述了GNSS大气水汽探测遥感技术及其应用的发展过程;然后介绍了近年来包括对流层延迟、大气可降水量等多类型GNSS大气参数高精度反演的研究进展,特别是对GNSS大气反演在极端天气短临预报及气候变化现象解释两个方向的研究工作进行了科学探析;最后,阐明了GNSS大气水汽探测遥感技术面临的主要挑战及未来研究展望。     

Using DORIS measurements for modeling the vertical total electron content of the Earth’s ionosphere

The Doppler orbitography and radiopositioning integrated by satellite (DORIS) system was originally developed for precise orbit determination of low Earth orbiting (LEO) satellites. Beyond that, it is highly qualified for modeling the distribution of electrons within the Earth’s ionosphere. It measures with two frequencies in L-band with a relative frequency ratio close to 5. Since the terrestrial ground beacons are distributed quite homogeneously and several LEOs are equipped with modern receivers, a good applicability for global vertical total electron content (VTEC) modeling can be expected. This paper investigates the capability of DORIS dual-frequency phase observations for deriving VTEC and the contribution of these data to global VTEC modeling. The DORIS preprocessing is performed similar to commonly used global navigation satellite systems (GNSS) preprocessing. However, the absolute DORIS VTEC level is taken from global ionospheric maps (GIM) provided by the International GNSS Service (IGS) as the DORIS data contain no absolute information. DORIS-derived VTEC values show good consistency with IGS GIMs with a RMS between 2 and 3 total electron content units (TECU) depending on solar activity which can be reduced to less than 2 TECU when using only observations with elevation angles higher than \(50^\circ \) . The combination of DORIS VTEC with data from other space-geodetic measurement techniques improves the accuracy of global VTEC models significantly. If DORIS VTEC data is used to update IGS GIMs, an improvement of up to 12  % can be achieved. The accuracy directly beneath the DORIS satellites’ ground-tracks ranges between 1.5 and 3.5 TECU assuming a precision of 2.5 TECU for altimeter-derived VTEC values which have been used for validation purposes.     

The ionosphere: effects,GPS modeling and the benefits for space geodetic techniques

The main goal of this paper is to provide a summary of our current knowledge of the ionosphere as it relates to space geodetic techniques, especially the most informative technology, global navigation satellite systems (GNSS), specifically the fully deployed and operational global positioning system (GPS). As such, the main relevant modeling points are discussed, and the corresponding results of ionospheric monitoring are related, which were mostly computed using GPS data and based on the direct experience of the authors. We address various phenomena such as horizontal and vertical ionospheric morphology in quiet conditions, traveling ionospheric disturbances, solar flares, ionospheric storms and scintillation. Finally, we also tackle the question of how improved knowledge of ionospheric conditions, especially in terms of an accurate understanding of the distribution of free electrons, can improve space geodetic techniques at different levels, such as higher-order ionospheric effects, precise GNSS navigation, single-antenna GNSS orientation and real-time GNSS meteorology.     

Variational data assimilation for deriving global climate analyses from GNSS radio occultation data

A comprehensive global navigation satellite system (GNSS) based radio occultation (RO) data set is available for meteorology and climate applications since the start of GNSS RO measurements aboard the CHAllenging Mini-satellite Payload (CHAMP) satellite in February 2001. Global coverage, all-weather capability, long-term stability and accuracy not only makes this innovative use of GNSS signals a valuable supplement to the data set assimilated into numerical weather prediction (NWP) systems but also an excellent candidate for global climate monitoring. We present a 3D variational data assimilation (3D-Var) scheme developed to derive consistent global analysis fields of temperature, specific humidity, and surface pressure from GNSS RO data. The system is based on the assimilation of RO data within 6 h time windows into European Centre for Medium-Range Weather Forecasts (ECMWF) short-term (24 h, 30 h) forecasts, to derive climatologic monthly mean fields. July 2003 was used as a test-bed for assessing the system’s performance. The results show good agreement with climatologies derived from RO data only and recent NWP impact studies. These findings are encouraging for future developments to apply the approach for longer term climatologic analyses, validation of other data sets, and atmospheric variability studies.     

Detecting discontinuities in GNSS coordinate time series with STARS: case study,the Bologna and Medicina GPS sites

Global navigation satellite systems (GNSS) data are a fundamental source of information for achieving a better understanding of geophysical and climate-related phenomena. However, discontinuities in the coordinate time series might be a severe limiting factor for the reliable estimate of long-term trends. A methodological approach has been adapted from Rodionov (Geophys Res Lett 31:L09204, 2004; Geophys Res Lett 31:L12707, 2006) and from Rodionov and Overland (J Marine Sci 62:328–332, 2005) to identify both the epoch of occurrence and the magnitude of jumps corrupting GNSS data sets without any a priori information on these quantities. The procedure is based on the Sequential t test Analysis of Regime Shifts (STARS) (Rodionov in Geophys Res Lett 31:L09204, 2004). The method has been tested against a synthetic data set characterized by typical features exhibited by real GNSS time series, such as linear trend, seasonal cycle, jumps, missing epochs and a combination of white and flicker noise. The results show that the offsets identified by the algorithm are split into 48 % of true-positive, 28 % of false-positive and 24 % of false-negative events. The procedure has then been applied to GPS coordinate time series of stations located in the southeastern Po Plain, in Italy. The series span more than 15 years and are affected by offsets of different nature. The methodology proves to be effective, as confirmed by the comparison between the corrected GPS time series and those obtained by other observation techniques.     

Global morphology of ionospheric F-layer scintillations using FS3/COSMIC GPS radio occultation data

We report on the FormoSat-3/Constellation Observing System for Meteorology, Ionosphere and Climate (FS3/COSMIC) limb-viewing observations of GPS L-band scintillations since mid-2006 and propose to study global F-layer irregularity morphology. The FS3/COSMIC has generally performed more than 1000 ionospheric radio occultation (RO) observations per day. We reprocess 1-Hz amplitude data and obtain complete limb-viewing profiles of the undersampling (sampling frequency lower than Fresnel frequency) S4 scintillation index from about 80% of the RO observations. There are a few percent of FS3/COSMIC RO observations having greater than 0.09 undersampling S4max values on average. However, seven identified areas, Central Pacific Area (?20° to 20° dip latitude, 160°E–130°W), South American Area (?20° to 20° dip latitude, 100°W–30°W), African Area (?20° to 20° dip latitude, 30°W–50°E), European Area (30°–55°N, 0°–55°E), Japan Sea Area (35°–55°N, 120°–150°E), Arctic Area (>65° dip latitude), and Antarctic Area (<?65° dip latitude), have been designated to have a much higher percentage of strong limb-viewing L-band scintillations. During the years in most of the last sunspot cycle from mid-2006 to the end 2014, the scintillation climatology, namely, its variations with each identified area, season, local time, magnetic activity, and solar activity, have been documented.     

Estimation and evaluation of real-time precipitable water vapor from GLONASS and GPS

The revitalized Russian GLONASS system provides new potential for real-time retrieval of zenith tropospheric delays (ZTD) and precipitable water vapor (PWV) in order to support time-critical meteorological applications such as nowcasting or severe weather event monitoring. In this study, we develop a method of real-time ZTD/PWV retrieval based on GLONASS and/or GPS observations. The performance of ZTD and PWV derived from GLONASS data using real-time precise point positioning (PPP) technique is carefully investigated and evaluated. The potential of combining GLONASS and GPS data for ZTD/PWV retrieving is assessed as well. The GLONASS and GPS observations of about half a year for 80 globally distributed stations from the IGS (International GNSS Service) network are processed. The results show that the real-time GLONASS ZTD series agree quite well with the GPS ZTD series in general: the RMS of ZTD differences is about 8 mm (about 1.2 mm in PWV). Furthermore, for an inter-technique validation, the real-time ZTD estimated from GLONASS-only, GPS-only, and the GPS/GLONASS combined solutions are compared with those derived from very long baseline interferometry (VLBI) at colocated GNSS/VLBI stations. The comparison shows that GLONASS can contribute to real-time meteorological applications, with almost the same accuracy as GPS. More accurate and reliable water vapor values, about 1.5–2.3 mm in PWV, can be achieved when GLONASS observations are combined with the GPS ones in the real-time PPP data processing. The comparison with radiosonde data further confirms the performance of GLONASS-derived real-time PWV and the benefit of adding GLONASS to stand-alone GPS processing.     

Ionospheric applications of the scintillation and tomography receiver in space (CITRIS) mission when used with the DORIS radio beacon network

The scintillation and tomography receiver in space (CITRIS) instrument will orbit the Earth near 560 km altitude to detect signals from the ground-based array of more than 50 DORIS UHF/S-band radio beacons established at sites around the world by the French Centre National d‘Etudes Spatiales (CNES) and the Institut Géographique National (IGN). The CITRIS receiver is on the US Air Force Space Test Program satellite STPSAT1, which is scheduled for launch in November 2006. CITRIS will record ionospheric total electron content (TEC) and radio scintillations with a unique ground-to-space geometry. The new instrument has been developed to study the ionosphere using data obtained with the UHF and S-band radio transmissions from the DORIS beacons because ionospheric radio scintillations can seriously degrade the performance of many space-geodetic systems, including the DORIS precise satellite orbitography system and GNSS (Global Navigation Satellite Systems). The ionospheric data will be based on radio signals sampled at a rate of 200 Hz by the CITRIS receiver. Numerical models have been used to predict that the DORIS signals measured by CITRIS may have 30 dB fluctuations in amplitude and 30 rad in phase as the satellite flies over kilometer-scale ionospheric structures. The data from the space-based CITRIS receiver will help update and validate theories on the generation and effect of ionospheric irregularities known to influence radio systems. By using simultaneous beacon transmissions from DORIS on the ground and from low-Earth-orbit beacons in space, the concept of reciprocity in a non-bilateral propagation medium like the ionosphere will be tested. Computer simulations are used to predict the magnitude of amplitude and phase scintillations that are expected to be recorded with the CITRIS instrument.     

Remote leveling of tide gauges using GNSS reflectometry: case study at Spring Bay,Australia

We further developed a new approach using GNSS reflectometry to determine the leveling connection between a tide gauge and a GNSS antenna. This approach includes the optimization of the unknown receiver bandwidth and the estimation of frequency changes in the signal-to-noise ratio (SNR) oscillation through an extended Kalman filter/smoother algorithm. We also corrected the geometric bending of the GNSS signals due to tropospheric refraction using local meteorological observations. Using 3 weeks of SNR data in Spring Bay, Australia, from a GNSS antenna placed sideways (i.e., ground plane orientated vertically and directed in azimuth toward the sea surface) to improve the SNR interference near the horizon, we obtained mean leveling differences of approximately 5 mm, with an RMS of approximately 3 cm level with respect to the nominal leveling from classical surveying techniques. SNR data from three different receiver manufacturers, coupled to the same antenna, provided similar leveling results. With a second antenna in the usual upright configuration, we obtained mean leveling differences of 1–2 cm and a RMS of about 10 cm. In the upright configuration, the leveling differences may include errors in the GNSS antenna phase center calibration, which are avoided in our technique but not in the classical surveying techniques. These results demonstrate the usefulness of the reflectometry technique to obtain precisely and remotely the leveling between a GNSS antenna and a tide gauge. In addition, this technique can be applied continuously, providing an independent and economical means to monitor the stability of the tide gauge zero.     

Investigation on the effect of strong out-of-band signals on global navigation satellite systems receivers

The mitigation of radio frequency interference (RFI) has a fundamental role in global navigation satellite system (GNSS) applications, especially when a high level of availability is required. Several electromagnetic sources, in fact, might degrade the performance of the global positioning system (GPS) and Galileo receivers, and their effects can be either in-band (i.e., secondary harmonics generated by transmitters of other communication systems due to non-linearity distortions) or out-of-band (i.e., strong signals that occupy frequency bandwidths very close to GNSS bands). We investigated the effects of real out-of-band signals on GNSS receivers and analyzed the impact on the overall receiver chain in order to evaluate the impact of the interference source. In particular, the analysis focuses on the spectrum at the front-end output, on the automatic gain control (AGC) behavior, as well as on the digital processing stages (signal acquisition and tracking) at the analog digital converter (ADC) output. This study refers to several experiments and data collections performed in interfered areas of downtown Torino (Italy). The obtained results underline how digital/analog TV transmissions represent a potential interference source for GNSS applications and might be critical for the safety of life services.     

A preliminary evaluation of the performance of multiple ionospheric models in low- and mid-latitude regions of China in 2010–2011

Ionospheric delay is a dominant error source in Global Navigation Satellite System (GNSS). Single-frequency GNSS applications require ionospheric correction of signal delay caused by the charged particles in the earth’s ionosphere. The Chinese Beidou system is developing its own ionospheric model for single-frequency users. The number of single-frequency GNSS users and applications is expected to grow fast in the next years in China. Thus, developing an appropriate ionospheric model is crucially important for the Chinese Beidou system and worldwide single-frequency Beidou users. We study the performance of five globally accessible ionospheric models Global Ionospheric Map (GIM), International Reference Ionosphere (IRI), Parameterized Ionospheric Model (PIM), Klobuchar and NeQuick in low- and mid-latitude regions of China under mid-solar activity condition. Generally, all ionospheric models can reproduce the trend of diurnal ionosphere variations. It is found that all the models have better performances in mid-latitude than in low-latitude regions. When all the models are compared to the observed total electron content (TEC) data derived from GIM model, the IRI model (2012 version) has the best agreement with GIM model and the NeQuick has the poorest agreement. The RMS errors of the IRI model using the GIM TEC as reference truth are about 3.0–10.0 TECU in low-latitude regions and 3.0–8.0 TECU in mid-latitude regions, as observed during a period of 1 year with medium level of solar activity. When all the ionospheric models are ingested into single-frequency precise point positioning (PPP) to correct the ionospheric delays in GPS observations, the PIM model performs the best in both low and mid-latitudes in China. In mid-latitude, the daily single-frequency PPP accuracy using PIM model is ~10 cm in horizontal and ~20 cm in up direction. At low-latitude regions, the PPP error using PIM model is 10–20 cm in north, 30–40 cm in east and ~60 cm in up component. The single-frequency PPP solutions indicate that NeQuick model has the lowest accuracy among all the models in both low- and mid-latitude regions of China. This study suggests that the PIM model may be considered for single-frequency GNSS users in China to achieve a good positioning accuracy in both low- and mid-latitude regions.     

GNSS satellite transmit power and its impact on orbit determination

Antenna thrust is a small acceleration acting on Global Navigation Satellite System satellites caused by the transmission of radio navigation signals. Knowledge about the transmit power and the mass of the satellites is required for the computation of this effect. The actual transmit power can be obtained from measurements with a high-gain antenna and knowledge about the properties of the transmit and receive antennas as well as losses along the propagation path. Transmit power measurements for different types of GPS, GLONASS, Galileo, and BeiDou-2 satellites were taken with a 30-m dish antenna of the German Aerospace Center (DLR) located at its ground station in Weilheim. For GPS, total L-band transmit power levels of 50–240 W were obtained, 20–135 W for GLONASS, 95–265 W for Galileo, and 130–185 W for BeiDou-2. The transmit power differs usually only slightly for individual spacecraft within one satellite block. An exception are the GLONASS-M satellites where six subgroups with different transmit power levels could be identified. Considering the antenna thrust in precise orbit determination of GNSS satellites decreases the orbital radius by 1–27 mm depending on the transmit power, the satellite mass, and the orbital period.     

Ionospheric TEC estimation with the signals of various geostationary navigational satellites

GPS Solutions - With the development of receiver equipment and GNSS and SBAS constellations, the coherent dual-frequency L-band transmissions are now available from a number of geostationary...     

Brazilian active GNSS networks as systems for monitoring the ionosphere

This research shows the viability of using Global Navigation Satellite System (GNSS) stations from Brazilian active networks in monitoring the ionosphere. Various indexes of ionospheric irregularities and scintillation of GNSS signals, estimated in real-time and post-processed from GNSS data, are explored for this purpose. This way, an increase in the spatial resolution of ionospheric information is provided, allowing the generation of maps of scintillation and irregularities in observing the spatial and temporal behavior of the layer’s activity cycle, since the number of ionosondes, imagers, and radars is insufficient for monitoring the irregularities in Brazil. Experiments to evaluate the estimates of the indexes are performed for periods of high and low variability of electrons. Three Brazilian networks are used: the Brazilian Network for Continuous Monitoring (RBMC), the GNSS Active Network of Sao Paulo State (GNSS-SP), and CIGALA/CALIBRA. The results are compared with data from ionosondes and PolaRxS-PRO Septentrio receivers, proving compatible with moderate to high correlations. An analysis of the seasonal variation during the peak of solar cycle 24 is carried out. The maps allow identifying the displacement of ionospheric irregularities along the magnetic equator over Brazil, from northeast to southwest, starting at 7:00 pm and ending at 2:00 am local time. Real-time monitoring is carried out for the summer solstice in the southern hemisphere, and results are consistent with those from the post-processed mode. The indexes and maps can be applied to the analysis of GNSS positioning. Real-time ionospheric information can be used in important practical applications because the displacement monitoring of irregularities allows prior knowledge of whether there will be a deterioration of positioning accuracy in a certain region.     

GNSS卫星精密定轨综述:现状、挑战与机遇

GNSS卫星精密轨道是高精度GNSS应用的基础与前提,GNSS卫星精密定轨技术也一直都是卫星导航领域的研究重点与热点。本文首先介绍了GNSS星座与跟踪数据概况,梳理了精密定轨函数模型、动力学模型及随机模型构建过程中的关键问题,归纳了低轨星载观测和星间链路观测等多源数据增强GNSS精密定轨的研究进展;然后,从应用的角度总结了当前GNSS精密轨道产品的基本状态,并进行了精度评估;最后,讨论了GNSS精密定轨在大网快速解算、多层次观测数据融合、太阳光压模型精化及高精度实时定轨等方面所面临的挑战,并展望了低轨星座、光钟、激光链路等新技术给GNSS精密定轨带来的机遇。     

Tropospheric refractivity and zenith path delays from least-squares collocation of meteorological and GNSS data

Precise positioning requires an accurate a priori troposphere model to enhance the solution quality. Several empirical models are available, but they may not properly characterize the state of troposphere, especially in severe weather conditions. Another possible solution is to use regional troposphere models based on real-time or near-real time measurements. In this study, we present the total refractivity and zenith total delay (ZTD) models based on a numerical weather prediction (NWP) model, Global Navigation Satellite System (GNSS) data and ground-based meteorological observations. We reconstruct the total refractivity profiles over the western part of Switzerland and the total refractivity profiles as well as ZTDs over Poland using the least-squares collocation software COMEDIE (Collocation of Meteorological Data for Interpretation and Estimation of Tropospheric Pathdelays) developed at ETH Zürich. In these two case studies, profiles of the total refractivity and ZTDs are calculated from different data sets. For Switzerland, the data set with the best agreement with the reference radiosonde (RS) measurements is the combination of ground-based meteorological observations and GNSS ZTDs. Introducing the horizontal gradients does not improve the vertical interpolation, and results in slightly larger biases and standard deviations. For Poland, the data set based on meteorological parameters from the NWP Weather Research and Forecasting (WRF) model and from a combination of the NWP model and GNSS ZTDs shows the best agreement with the reference RS data. In terms of ZTD, the combined NWP-GNSS observations and GNSS-only data set exhibit the best accuracy with an average bias (from all stations) of 3.7 mm and average standard deviations of 17.0 mm w.r.t. the reference GNSS stations.     

An efficient solution of real-time data processing for multi-GNSS network

Global navigation satellite systems (GNSS) are acting as an indispensable tool for geodetic research and global monitoring of the Earth, and they have been rapidly developed over the past few years with abundant GNSS networks, modern constellations, and significant improvement in mathematic models of data processing. However, due to the increasing number of satellites and stations, the computational efficiency becomes a key issue and it could hamper the further development of GNSS applications. In this contribution, this problem is overcome from the aspects of both dense linear algebra algorithms and GNSS processing strategy. First, in order to fully explore the power of modern microprocessors, the square root information filter solution based on the blocked QR factorization employing as many matrix–matrix operations as possible is introduced. In addition, the algorithm complexity of GNSS data processing is further decreased by centralizing the carrier-phase observations and ambiguity parameters, as well as performing the real-time ambiguity resolution and elimination. Based on the QR factorization of the simulated matrix, we can conclude that compared to unblocked QR factorization, the blocked QR factorization can greatly improve processing efficiency with a magnitude of nearly two orders on a personal computer with four 3.30 GHz cores. Then, with 82 globally distributed stations, the processing efficiency is further validated in multi-GNSS (GPS/BDS/Galileo) satellite clock estimation. The results suggest that it will take about 31.38 s per epoch for the unblocked method. While, without any loss of accuracy, it only takes 0.50 and 0.31 s for our new algorithm per epoch for float and fixed clock solutions, respectively.     

Applying GNSS and DTM Technologies to Monitor the Ice Balance of the Horcones Inferior Glacier, Aconcagua Region, Argentina

This work provides results and analysis of the behavior in a reconstituted debris covered glacier, Horcones Inferior (HIG). This glacier is located at 32° 41′ S and 69° 57′ W in the Aconcagua region, Mendoza, Argentina. The HIG has experienced surges events in 1984 and 2003. GNSS techniques and Digital Terrain Models (DTMs) derived from ASTER optical imagery, were used to detect movements and altimetric changes of the HIG in the 2001–2012 period. Based on a GNSS semi - continuous station, N and E mean velocities of 0.9 cm/d and 1.1 cm/d and in U direction were 1.1 cm/d and 1.5 cm/d were obtained respectively in 2009–2010 period. Additionally, GNSS profiles showed a velocity range between 0.5 cm/d and 2.8 cm/d. Finally, through a DTM comparison the altimetric change before and after the last surge in 2003 was obtained. The applied techniques allowed accurate and reliable change detection of the HIG.