Evaluation of the Multipath Environment Using Electromagnetic-Absorbing Materials at Continuous GNSS Stations

o date, no universal modelling technique is available to mitigate the effect of site-specific multipaths in high-precision global navigation satellite system (GNSS) data processing. Multipaths affect both carrier-phase and code/pseudorange measurements, and the errors can propagate and cause position biases. This paper presents the use of an Eccosorb AN-W-79 microwave-absorbing material mounted around a GNSS antenna that reflects less than −17 dB of normal incident energy above a frequency of 600 MHz. To verify the feasibility and effectiveness of the Eccosorb, we installed two close stations by continuously operating multi-GNSS (BeiDou, GLONASS, Galileo and GPS) in a challenging location. One station is equipped with the Eccosorb AN-W-79, covering a square area of 3.35 m2 around the antenna, and the second station operates without it. The standard deviation reductions from single point positioning estimates are significant for all the individual GNSS solutions for the station equipped with microwave-absorbing material. The reductions are as follows: for GPS, between 15% and 23%; for Galileo, between 22% and 45%; for GLONASS, 22%; and for BeiDou, 4%. Furthermore, we assess the influence of multipaths by analysing the linear combinations of code and carrier phase measurements for various GNSS frequencies. The Galileo code multipath shows a reduction of more than 60% for the station with microwave-absorbing material. For GLONASS, particularly for the GLOM3X and GLOM1P code multipath combinations, the reduction reaches 50%, depending on the observation code types. For BeiDou, the reduction is more than 30%, and for GPS, it reaches between 20% and 40%. The Eccosorb AN-W-79 microwave-absorbing material shows convincing results in reducing the code multipath noise level. Again, using microwave-absorbing material leads to an improvement between 15% and 60% in carrier phase cycle slips. The carrier-phase multipath contents on the post-fit residuals from the processed GNSS solutions show a relative RMS reduction of 13% for Galileo and 9% for GLONASS and GPS when using the microwave-absorbing material. This study also presents power spectral contents from residual signal-to-noise ratio time series using Morlet wavelet transformation. The power spectra from the antenna with the Eccosorb AN-W-79 have the smallest magnitude, demonstrating the capacity of microwave-absorbing materials to lessen the multipath influence while not eliminating it

Low-cost GNSS for Meteorology: A receiver comparison

GNSS, particularly the GPS and GLONASS satellite constellations are commonly used for accurate estimation of atmospheric parameters. The high accuracy is accomplished by sophisticated analysis methods and expensive high-end receivers and antennas along with precise products and bias corrections. The recent market introduction of low-cost dual frequency receivers that can produce raw data from multiple constellations offers an insight into the potential use of these receivers for meteorological applications. Here we demonstrate that GPS and GLONASS measurements from a low-cost dual-frequency receiver can be used to estimate the Zenith Total Delay ( ZTD) commensurate to meteorological applications

A New Global Vertical Land Movement Data Set from the TIGA Combined

Globally averaged sea level has been estimated from the network of tide gauges installed around the world since the 19th century. These mean sea level (MSL) records provide sea level relative to a nearby tide gauge benchmark (TGBM), which allows for the continuation of the instrumental record in time. Any changes in the benchmark levels, induced by vertical land movements (VLM) affect the MSL records and hence sea level estimates. Over the last two decades sea level has also been observed using satellite altimeters. While the satellite observations are globally more homogeneous providing a picture of sea level not confined to coastlines, they require the VLM-corrected MSL records for the bias calibration of instrumental drifts. Without this calibration altimeter instruments from different missions cannot be combined. GPS has made it possible to obtain highly accurate estimates of VLM in a geocentric reference frame for stations at or close to tide gauges. Under the umbrella of the International GNSS Service (IGS), the Tide Gauge Benchmark Monitoring (TIGA) Working Group (WG) has been established to apply the expertise of the GNSS community to solving issues related to the accuracy and reliability of the vertical component to provide estimates of VLM in a well-defined global reference frame. To achieve this objective, five TIGA Analysis Centers (TACs) contributed re-processed global GPS network solutions to TIGA, employing the latest bias models and processing strategies in accordance with the second re-processing campaign (repro2) of the IGS. These solutions include those of the British Isles continuous GNSS Facility – University of Luxembourg consortium (BLT), the German Research Centre for Geosciences (GFZ) Potsdam, the German Geodetic Research Institute (DGF) at the Technical University of Munich, Geoscience Australia (AUT) and the University of La Rochelle (ULR). In this study we present to the sea level community an evaluation of the VLM estimates from the first combined solution from the IGS TIGA WG. The TAC solutions include more than 700 stations and span the common period 1995-2014. The combined solution was computed by the TIGA Combination Centre (TCC) at the University of Luxembourg, which used the Combination and Analysis of Terrestrial Reference Frame (CATREF) software package for this purpose. This first solution forms Release 1.0 and further releases will be made available after further reprocessing campaigns. We evaluate the combined solution internally using the TAC solutions and externally using solutions from the IGS and the ITRF2008. The derived VLM estimates have undergone an initial evaluation and should be considered as the primary TIGA product for the sea level community to correct MSL records for land level change

First Results of the GNSS Reprocessing at the University of Luxembourg for an Updated Global Vertical Land Movement Data Set

The University of Luxembourg (UL) is currently contributing to the most recent reprocessing effort of the International GNSS Service (IGS) Tide Gauge Benchmark Monitoring Working Group (TIGA-WG) with multi-constellation GNSS solutions, including GPS, GLONASS and Galileo. As part of this new reprocessing and reanalysis effort of GNSS data including stations at or near tide gauges worldwide, several model enhancements consistent with the IGS’s recent effort have been incorporated. During 1994 to 2022 the network generally contains data from over 700 stations. The IGS has placed high importance on unifying processing standards because homogeneous and consistent reprocessing of all GNSS data over the complete time span is necessary for estimating useful geophysical parameters, such as long-term trends in station positions. In addition to the reprocessing of the data, the time series analysis strategy is crucial for deriving accurate long-term estimates. In order to obtain the best parameter estimates and the most realistic uncertainties, it is anticipated that a number of stochastic and deterministic models will be fitted to the position time series. Additionally, the stochastic properties of the series will be investigated. Existing automated processes will be updated with the most recent developments in geodetic time series analysis due to the vast number of stations. Through the full reprocessing of all GNSS observations and the state-of-the-art analysis of the daily position time series, this study will be able to obtain highly accurate estimates of horizontal and vertical land movements that can be employed for the most challenging applications such as correcting coastal sea level records for a better understanding in their changes and constraining glacial isostatic adjustment models. During this presentation, we will provide details on the current reprocessing, present preliminary results and a first cross-evaluation of the vertical land movement estimates

Tracking hurricanes Harvey and Irma using GPS tropospheric products

The 2017 Hurricanes season was one of the most powerful severe weather events producing catastrophic socio-economic and environmental effects on the east coast of the United States. Therefore, tracking their path accurately is extremely useful. Today Global Navigation Satellite Systems (GNSS) tropospheric products, such as Zenith Wet Delays (ZWD), and Integrated Water Vapor (IWV) are used as complementary data sets in Numerical Weather Prediction (NWP) models. In this study, we employed GPS-derived IWV and horizontal tropospheric gradient information to monitor and investigate the complicated characteristics of hurricane events in their spatial and temporal distribution using a dense ground network of GPS stations. Our results show that a surge in GPS-derived IWV occurred several hours prior to the manifestation of the major hurricanes Harvey and Irma. We used the derived GPS-derived IWV information as input to spaghetti lines weather models, allowing us to predict the paths of Harvey and Irma hurricanes. As such, a parameter directly estimated from GPS can provide an additional resource for improving the monitoring of hurricane path

Signal Obstructions at GNSS Stations: Benefits From Multi-GNSS Observations

The current accuracy of IGS products, few centimeter level, requires amongst other things that the location for GNSS antennas are nearly optimal for GNSS observations. This includes a low multipath environment and little to no signal obstructions. However, this is not guaranteed for every station especially in urban areas and mountainous regions. As some applications such as GNSS for sea level studies or to monitor landslides require GNSS antennas to be installed at a specific site, it is clear that the environment might not be favourable for GNSS observations. In this study, we investigate the effect of signal obstructions on station positions, specifically the height component, based on simulated obstruction scenarios using a modified Bernese GNSS Software version 5.2 (BSW52). The behaviours of different obstruction scenarios and the impact of multi-GNSS (GPS+GLONASS for now) observations for both clear and obstructed stations are discussed

Multi-GNSS Slant Wet Delay Retrieval Using Multipath Mitigation Maps

The conventional Global Navigation Satellite System (GNSS) processing is typically contaminated with errors due to atmospheric variabilities, such as those associated with the mesoscale phenomena. These errors are manifested in the parameter estimates, including station coordinates and atmospheric products. To enhance the accuracy of these GNSS products further, a better understanding of the local-scale atmospheric variability is necessary. As part of multi-GNSS processing, station coordinates, carrier phase ambiguities, orbits, zenith total delay (ZTD) and horizontal gradients are the main parameters of interest. Here, ZTD is estimated as the average zenith delay along the line-of-sight to every observed GNSS satellite mapped to the vertical while the horizontal gradients are estimated in NS and EW directions and provide a means to partly account for the azimuthally inhomogeneous atmosphere. However, a better atmospheric description is possible by evaluating the slant path delay (SPD) or slant wet delay (SWD) along GNSS ray paths, which are not resolved by ordinary ZTD and gradient analysis. SWD is expected to provide better information about the inhomogeneous distribution of water vapour that is disregarded when retrieving ZTD and horizontal gradients. Usually, SWD cannot be estimated directly from GNSS processing as the number of unknown parameters exceeds the number of observations. Thus, SWD is generally calculated from ZTD for each satellite and may be dominated by un-modelled atmospheric delays, clock errors, unresolved carrier-phase ambiguities and near-surface multipath scattering. In this work, we have computed multipath maps by stacking individual post-fit carrier residuals incorporating the signals from four GNSS constellations, i.e. BeiDou, Galileo, Glonass and GPS. We have selected a subset of global International GNSS Service (IGS) stations capable of multi-GNSS observables located in different climatic zones. The multipath effects are reduced by subtracting the stacked multipath maps from the raw post-fit carrier phase residuals. We demonstrate that the multipath stacking technique results in significantly reduced variations in the one-way post-fit carrier phase residuals. This is particularly evident for lower elevation angles, thus, producing a retrieval method for SWD that is less affected by site-specific multipath effects. We show a positive impact on SWD estimation using our multipath maps during increased atmospheric inhomogeneity as induced by severe weather events

Multi-GNSS Precise Point Positioning with Ambiguity Resolution using IGS Products: A Cross-Evaluation of Three Software Packages and Combinations of Constellation Data for Geophysical Investigations

Precise point positioning (PPP) is a GNSS processing strategy for the estimation of station coordinates, receiver clocks, and zenith total delays with high efficiency, precision, and accuracy. However, to achieve the highest quality in PPP solutions, precise satellite orbits and clocks products, and various bias information for code and carrier observations, necessary for PPP integer ambiguity resolution (AR), are essential. These critical products are routinely provided by the International GNSS Service (IGS) and its analysis centers. Due to the undifferenced nature of PPP solutions, they are particularly affected by error sources and in case of insufficient mitigation models, residual errors may propagate into the solutions and contaminate the time series of the estimates. Also, inconsistencies between model implementations in software used during product generation and the ones employed by users during product application, lead to artifacts in the PPP solutions. On the other hand, PPP AR does compensate somewhat and improve the solutions despite the frequently inseparable sources of residual errors. In this study, we present an analysis of single and multi-GNSS PPP AR solutions using observations from the GPS, GLONASS, and Galileo constellations. We employed three independent software packages, the Bernese GNSS Software v5.4, PRIDE v2.2, and GipsyX v2.1, each employing their recommended set of products and processing settings while attempting to keep settings as consistent as possible between the software packages and processing runs. We processed data from 100 globally distributed IGS stations, carefully selected for known quality and network geometry, for 2019.0 – 2023.5. In our analysis, we evaluate the impact of different single- and multi-constellation GNSS PPP-AR solutions in terms of time series characteristics, mean coordinate differences, and solution repeatabilities for the different software packages

Total Impact of Periodic Terms and Coloured Noise on Velocity Estimates

The uncertainties of velocity estimates for position time series of Global Navigation Satellite System (GNSS) stations are mainly affected by a misfit of the deterministic model applied to this data. Insufficiently modelled seasonal signals will propagate into the stochastic model and falsify the results of the noise analysis besides the velocity estimates and their uncertainties. In this presentation we derived the General Dilution of Precision (GDP) of velocity uncertainties. We define this dilution as the ratio between the uncertainties of velocities determined when different deterministic and stochastic models are applied. In this way we discuss, referring to previously published results, how insufficiently modelled seasonal signals influence station velocity uncertainties with white and coloured noise. Using simulated and real data from selected (115) IGS (International GNSS Service) stations we show that the noise character affects GNSS data more than seasonals for time series longer than 9 years

Error analysis of Tide Gauge Benchmark Monitoring (TIGA) Analysis Center stacked solutions

In 2013 the International GNSS Service (IGS) Tide Gauge Benchmark Monitoring (TIGA) Working Group (WG) started their reprocessing campaign, which proposes to re-analyze all relevant Global Positioning System (GPS) observations from 1995 to the end of 2013. This re-processed dataset will provide high quality estimates of land motions, enabling regional and global high-precision geophysical/geodetic studies. Several of the individual TIGA Analysis Centers (TACs) have completed processing the full history of GPS observations recorded by the IGS global network, as well as, many other GPS stations at or close to tide gauges, which are available from the TIGA data center at the University of La Rochelle (www.sonel.org). The TAC solutions contain a total of over 700 stations. This study focuses on the evaluations of any systematic error present in the three TIGA analysis center (TAC) SINEX solutions: the British Isles continuous GNSS Facility – University of Luxembourg consortium (BLT), the GeoForschungsZentrum (GFZ) Potsdam, and of the University of La Rochelle (ULR). We have analyzed the residual position time series of the individual TAC a combination of automatic and manual discontinuity identification, applying a post-seismic deformation model adopted from ITRF2014 for those stations that are affected by earthquakes, followed by the stacking of the daily solution of the individual TAC into a long term linear frame. We have carried out the error analysis using the Combination and Analysis of Terrestrial Reference Frame (CATREF) software package. The TIGA Combination Centre (TCC) at the University of Luxembourg (UL) is responsible for providing a combined solution with a global set of vertical land movement estimates