Intrinsic Identification and Mitigation of Multipath for Enhanced GNSS Positioning

The study of climate change is an important field of research. Monitoring of atmospheric variability especially the tropospheric precipitable water vapor (PWV) is a powerful way to investigate climate change. Global navigation satellite systems (GNSS) provide a good tool for studying atmospheric parameters as GNSS signals along its path from the satellites to the ground based receivers suffer a significant delay due to the refractivity of earth’s atmosphere. GNSS signals are not only delayed but also refracted in the neutral atmosphere. The zenith wet delays (ZWD) caused by the troposphere can be estimated during the geodetic processing of GNSS signals. Since the ZWD is tightly correlated to the PWV, GNSS observations can be used to study the changes of PWV. In this study GNSS data from the Egyptian permanent GNSS network (EPGN), International GNSS Service (IGS), EUREF Permanent Network (EPN) and Scripps Orbit and Permanent Array Center (SOPAC), for the years 2013 and 2014, were used. These GNSS stations provided a good geometrical coverage over the respective region. All GNSS data were processed using Bernese V 5.2 software. The needed product from GNSS data processing is the tropospheric zenith total delay (ZTD) over each GNSS site. The ZTD is divided based on physical parameters into zenith hydrostatic delay (ZHD) and zenith wet delay (ZWD). The ZWD is the basic observable used to calculate PWV. The values of the PWV were calculated and its variation over the study area was investigated. The quality of calculated PWV values from GNSS data evaluated against traditional Radiosonde (RS) measurements. The results of GNSS PWV showed good agreement with RS PWV when the PWV values were low, but this agreement became worse at high PWV values as the differences between the two techniques increased. The correlation coefficient between RS PWV and GNSS PWV varied from 0.31 to 0.84. Standard deviation of the differences between RS PWV and GNSS PWV ranged from 2.369 to 5.973 mm. The PWV estimated from GNSS observations had annual cycle, and its cycle in 2013 was different from that in 2014. The PWV differences between the 2 years clarified that the water vapor content over the east Mediterranean in most of 2014 days was higher than that in 2013. Furthermore, PWV variations were noted on both temporal and spatial scales. The highest temporal variation value was 25.41 mm whereas the maximum value of the spatial variation was 19.67 mm. The present study illustrated the importance of using geodetic networks to provide atmospheric information.

Abstract. For the azimuth observation to be made at its magnetic observatories routinely, Japan Meteorological Agency (JMA) has adopted a traditional method based on Polaris sighting. Due to its difficulty to implement under overcast weather conditions and to its demand on observers, for overtime work into the evening, we are motivated to seek for an alternative method based on the Global Navigation Satellite Systems (GNSS) observation that might potentially relieve those two disadvantages. An experiment is made at Kakioka to assess the eligibility and effectiveness of the GNSS method for JMA's unmanned observatories, Memambetsu and Kanoya. The GNSS observations themselves achieve as high a precision as approximately 1 arcsec, as far as they are analyzed with Static mode. Derived from the results of GNSS observation and some supplementary horizontal angle measurements, the azimuth of the azimuth mark for the absolute measurement is determined with a precision of a few arcsecond, which is comparable to the azimuth precision achieved by the Polaris sighting. However, we end up with their significant difference by about 10 arcsec. We discuss this discrepancy to be possibly due to a local geoid gradient. The Polaris observation is made with a theodolite tilted in the gravitational direction, also known as the vertical line deviation, whereas the GNSS observations are based on the azimuth of the compliant ellipsoid plane.

Abstract. Utilizing ground control points (GCPs) to georeference photogrammetry-based point cloud data is a common practice in unmanned aerial system (UAS) mapping. Direct georeferencing or integrated sensor orientation (ISO) can be used to obtain georeferenced point clouds from UAS without relying heavily on GCPs. However, the accuracy of the point cloud may be impacted by the accuracy of the trajectory solution obtained by GNSS. To improve point cloud accuracy, post-processing kinematic (PPK) solutions can be applied to the UAS trajectory, which may provide higher accuracy than low-accuracy trajectory solutions and minimize the reliance on GCPs. This study compares the accuracy and precision of two different point clouds generated using different methods. One point cloud was generated using traditional photogrammetric methods with low accuracy Global Navigation Satellite System (GNSS) observations from the UAS and GCPs that have an average accuracy of one to two centimeters, while the other was generated using PPK trajectory solution for the UAS’s trajectory with two software: open-source Emlid Studio and the widely used Inertial Explorer. The use of PPK techniques in UAS mapping may have several potential benefits over traditional methods. By correcting the errors in the UAS's trajectory, a user may only need to depend on fewer ground control points, which can reduce the time and cost associated with fieldwork. This is particularly useful in areas that are difficult to access or have limited ground control point options, such as in urban or forested areas. To evaluate performance, a GNSS receiver is used to obtain measurements on checkpoints, which are used to assess the accuracies of the point clouds. In our experiments, the accuracy of the point clouds generated using PPK trajectory solution with high accuracy GCPs was found to be higher than those generated with low accuracy GNSS observations while aided with high accuracy ground control points. While the use of PPK with GCPs is generally expected to provide more accurate and reliable data than low-accuracy GNSS observations even after adjusting with GCPs, the number and distribution of GCPs can still significantly impact overall accuracy. Therefore, careful consideration of the number of GCPs and their placement is essential to achieve the desired level of efficiency and effectiveness in UAS mapping.

This study evaluates the suitability of compact, lowcost Global Navigation Satellite System (GNSS) receiver module in estimating the effects of tropospheric delay on accuracy of GNSS position solution at GNSS laboratory, the University of Burdwan, West Bengal, India with geographical coordinate of 23.2545°N, 87.8467°E, −11.402 m, using compact, low-cost uBlox F9P GNSS module and high cost Leica GR 50 geodetic receiver in March 2023. GNSS observation data were logged by the two receivers for 24 hours at 1Hz dual-frequency and converted into Rinex observation files in GPS-only, GLONASS-only, and GPS+GLONASS hybrid modes. The converted Rinex observation data were uploaded to the Natural Resources of Canada (NRCan) online service to obtain the user fixed values of latitude, longitude and attitude, and the wet and dry tropospheric delay files via email. The converted Rinex observation files were post-processed using GNSS positioning analysis tool that estimate the GNSS positioning accuracy parameters namely Distance Root Mean Square (2DRMS), Circle of Error Probable (CEP), Spherical Error Probable (SEP), and Mean Radian Spherical Error (MRSE) from the users defined known fixed values of latitude, longitude and altitude. The total delay obtained from the two GNSS receivers were compared with their positioning accuracy parameters. The results showed better position solution on observation days with the least total zenith delay in all the constellations considered. The Leica GR 50 shows an average total zenith delay of 2.476 m in GPS mode, 2.486 m in GLONASS mode and 2.4680 m in GPS+GLONASS hybrid mode while the uBlox F9P shows 2.477 m in GPS mode, 2.494 m in GLONASS mode and 2.4684 m in GPS+GLONASS hybrid mode of operation. The GPS+GLONASS hybrid mode of operation reduces the total zenith delay and enhanced the position solution. Hence, there is no distinct difference between the two GNSS receivers, suggesting that compact, low-cost, GNSS receiver modules can be used for tropospheric study.

Precise data processing from the Global Navigation Satellite Systems (GNSS) reference station network is mainly based on a combination of double-differenced carrier phase and code observations. This approach allows most of the measurement errors to be removed or reduced and is characterized as the most accurate method. However, creating observation differences between two receivers and two satellites increases the measurement noise of the observations by a factor of 2. As a result, it increases the impact of the incorrect definition of the noise characteristic on the results of the estimation of the unknowns in the positioning model. This is especially important in Multi-GNSS solutions, which integrate measurements from different systems, for which the stochastic parameters of observation may differ significantly. In this paper, the authors prepared a complex analysis of the noise type in double-differenced GNSS (GPS, GLONASS and Galileo) observations, both carrier phase and code ones, with a 1 s sampling interval. The Autocorrelation Function (ACF) method, the Lomb–Scargle (L-S) periodogram method, and the Allan variance (AVAR) method were used. The results that were obtained for the weekly set of measurement data showed that, depending on the system and type of observation, the noise level and its type are significantly different. Among the code measurements, the lowest noise levels were obtained for the GPS C5Q and Galileo C7Q/C8Q observations, with the standard deviations not exceeding ±10 cm, while the noisiest observations were for the GLONASS C1C and C2C signals, which had standard deviations of about ±90 cm and ±45 cm, respectively. For the carrier phase observations, each signal type was characterized by very similar noise levels of ±1.5–3.5 mm. The ACF analysis showed that 1 Hz double-differenced GNSS data can only be treated as being not correlated to time for carrier phase observations; for code observations, an irrelevant autocorrelation may be considered for measurement intervals greater than 20 s. Depending on the GNSS signals, the spectral index k varies in a range from −1.3 to −0.2 for code data and k = 0.0 in the case of phase data. Using the modified Allan deviation (MDEV) allows for specific noise types for each signal and GNSS system to be determined. All of the code observations were characterized by either flicker PM or white PM. In the case of the phase observations, they were all uniquely characterized by white PM (GPS and Galileo or by white PM and flicker PM (GLONASS).

A wide-range utilisation of Global Navigation Satellite Systems (GNSS) across technology and socioeconomic domains renders the satellite navigation one of the pillars of the modern civilisation. Tackling and overcoming The inherent shortcomings and emerging threats to provision of robust and resilient GNSS Positioning, Navigation and Timing (PNT) services have become a research subject of the utmost importance. An open access to the GNSS positioning estimation process is fundamental for development of advanced methods for robust and resilient GNSS position estimation. Recent introduction of Google Android 7.0 revision allows for direct access to raw GNSS pseudoranges observed by smartphones. Here we address the opportunities for GNSS positioning estimation improvements, given through exploitation of Google Android Location API, in creation of bespoke GNSS position estimation process in navigation application domain of a GNSS Software-Defined Radio (SDR) receivers in smartphones. We present and validate opportunities for improvement of GNSS position estimation process through utilisation of distributed computing architectures, trusted sources of GNSS augmentation data and utilisation of GNSS positioning methods suitable for targeted classes of GNSS applications Authors are members of European GNSS Agency's GNSS Raw Measurements Task Force, and co-authored the white-paper on GNSS-related Google Android Location API applications.

Tropospheric delay is an important error source in Global Navigation Satellite System (GNSS) positioning and can also be used in water vapor monitoring. Many models have been built to correct tropospheric delays or to convert zenith wet delays to precipitable water vapor. However, these models suffer from limited resolutions (spatial resolution lower than 1° and temporal resolution lower than 6 hr), which affects their performance. The release of European Centre for Medium‐Range Weather Forecasts ReAnalysis 5 (ERA5) provides the opportunity to lift this limit. In this study, we use the ERA5 hourly 0.5° × 0.5° data to build a new model over China, which integrates tropospheric delay correction for GNSS positioning and weighted mean temperature calculation for GNSS meteorology. By modeling the diurnal variations of zenith hydrostatic delay, zenith wet delay, and weighted mean temperature and the seasonal variations in their lapse rates, this model has the state‐of‐the‐art spatial resolution of 0.5° × 0.5° and temporal resolution of 1 hr. We validate this new model by the ERA5 data, the radiosonde data, and the GNSS data in comparison with the canonical GPT2w model. The results show that the new model has better accuracies in terms of root‐mean‐square than the GPT2w model in all parameters. Especially, the new model well captures the diurnal variations in tropospheric delay and weighted mean temperature. This new model provides accurate tropospheric delays and weighted mean temperature simultaneously, which enables GNSS receivers to measure precipitable water vapor directly and also benefits GNSS positioning.

The delays of radio signals transmitted by global navigation satellite system (GNSS) satellites and induced by neutral atmosphere, which are usually represented by zenith tropospheric delay (ZTD), are required as critical information both for GNSS positioning and navigation and GNSS meteorology. Establishing a stable and reliable ZTD model is one of the interests in GNSS research. In this study, we proposed a regional ZTD model that makes full use of the ZTD calculated from regional GNSS data and the corresponding ZTD estimated by global pressure and temperature 3 (GPT3) model, adopting the artificial neutral network (ANN) to construct the correlation between ZTD derived from GPT3 and GNSS observations. The experiments in Hong Kong using Satellite Positioning Reference Station Network (SatRet) were conducted and three statistical values, i.e., bias, root mean square error (RMSE), and compound relative error (CRE) were adopted for our comparisons. Numerical results showed that the proposed model outperformed the parameter ZTD model (Saastamoinen model) and the empirical ZTD model (GPT3 model), with an approximately 56%/52% and 52%/37% RMSE improvement in the internal and external accuracy verification, respectively. Moreover, the proposed method effectively improved the systematic deviation of GPT3 model and achieved better ZTD estimation in both rainy and rainless conditions.

It is anticipated that the performance of real‐time (RT) GNSS meteorology can be further improved by incorporating observations from multiple Global Navigation Satellite System (GNSS), including GPS, GLONASS, Galileo, and BeiDou. In this paper, an operational RT system for extracting zenith troposphere delay (ZTD) using a modified version of the Precise Point Positioning With Integer and Zero‐difference Ambiguity Resolution Demonstrator (PPP‐WIZARD) was established. GNSS, including GPS, GLONASS, and Galileo, observation streams were processed using RT Precise Point Positioning (PPP) strategy based on RT satellite orbit/clock products from the Centre National d'Etudes Spatiales. An experiment covering 30 days was conducted, in which the observation streams of 20 globally distributed stations were processed. The initialization time and accuracy of the RT troposphere results using single‐system and multisystem observations were evaluated. The effect of PPP ambiguity resolution was also evaluated. Results reveal that RT troposphere estimates based on single‐system observations can both be applied in weather nowcasting, in which the GPS‐only solution is better than the GLONASS‐only solution. The performance can also be improved by PPP ambiguity resolution and utilizing GNSS observations. Specifically, we notice that ambiguity resolution is more effective in improving the accuracy of ZTD, whereas the initialization process can be better accelerated by GNSS observations. Combining all techniques, the RT troposphere results with an average accuracy of about 8 mm in ZTD can be achieved after an initialization process of approximately 8.5 min, which demonstrates superior results for applying GNSS observations and ambiguity resolution for RT meteorological applications.

Due to the low signal power, the Global Navigation Satellite System (GNSS) signal is vulnerable to interference and even cannot be captured or tracked in harsh environments. As an alternative, the Low Earth Orbit (LEO) satellite has been widely used in the navigation field due to the advantages of low cost and strong signals. It is becoming a significant component of the new combined navigation system with GNSS. The combination of an LEO Doppler signal and GNSS observables can improve the positioning accuracy and high-precision positioning convergence time of the GNSS receiver. However, the GNSS signal receiving capability cannot be improved from this data fusion level. We propose a novel assisted structure where GNSS signal acquisition and Doppler tracking are assisted by LEO Doppler positioning. The receiver uses the LEO signal to achieve Doppler positioning firstly. Then, the coarse position with the GNSS navigation messages received from LEO, as well as the estimated clock information, is used to assist in the acquisition and tracking of GNSS. In this way, the GNSS receiver’s sensitivity can get the benefit from this integrated system. The paper presents the structure of the assisted receiver and analyzes the assisted GNSS signal acquisition and carrier tracking performance in detail. Simulation experiments of this assisted structure are carried out to verify its superiority of acquisition and tracking sensitivity in comparison with standalone GNSS receivers. Theoretical analysis and experimental results show that the proposed acquisition method can achieve 90% detection probability at a carrier-to-noise ratio (C/N0) of 15 dB-Hz, which is about 8 dB higher than the conventional acquisition method without assistance; the proposed tracking method can track weak signals of 5 dB-Hz, which is about 4 dB higher than the conventional method. Therefore, this novel LEO-assisted receiver has significantly improved weak signal acquisition and tracking sensitivity.

Landslide displacement prediction has garnered significant recognition as a pivotal component in realizing successful early warnings and implementing effective control measures. This task remains challenging as landslide deformation involves not only temporal dependency within time series data but also spatial dependence across various regions within landslides. The present study proposes a landslide spatiotemporal displacement forecasting model by introducing attention-based deep learning algorithms based on spatiotemporal analysis. The Maximal Information Coefficient (MIC) approach is employed to quantify the spatial and temporal correlations within the daily data of Global Navigation Satellite System (GNSS) observations. Based on the quantitative spatiotemporal analysis, the proposed prediction model combines a convolutional neural network (CNN) and long short-term memory (LSTM) network to capture spatial and temporal dependencies individually. Spatial–temporal attention mechanisms are implemented to optimize the model. Additionally, we develop a single-point prediction model using LSTM and a multiple-point prediction model using the CNN-LSTM without an attention mechanism to compare the forecasting capabilities of the attention-based CNN-LSTM model. The Outang landslide in the Three Gorges Reservoir Area (TGRA), characterized by a large and active landslide equipped with an advanced monitoring system, is taken as a studied case. The temporal MIC results shed light on the response times of monitored daily displacement to external factors, showing a lagging duration of between 10 and 50 days. The spatial MIC results indicate mutual influence among different locations within the landslide, particularly in the case of nearby sites experiencing significant deformation. The attention-based CNN-LSTM model demonstrates an impressive predictive performance across six monitoring stations within the Outang landslide area. Notably, it achieves a remarkable maximum coefficient of determination (R2) value of 0.9989, accompanied by minimum values for root mean squared error (RMSE), absolute mean error (MAE), and mean absolute percentage error (MAPE), specifically, 1.18 mm, 0.99 mm, and 0.33%, respectively. The proposed model excels in predicting displacements at all six monitoring points, whereas other models demonstrate strong performance at specific individual stations but lack consistent performance across all stations. This study, involving quantitative deformation characteristics analysis and spatiotemporal displacement prediction, holds promising potential for a more profound understanding of landslide evolution and a significant contribution to reducing landslide risk.

The central hypothesis of the Research Unit (RU) New Refined Observations of Climate Change from Spaceborne Gravity Missions (NEROGRAV), funded for the second three years phase by the German Research Foundation DFG, reads: only by concurrently improving and better understanding of sensor data, background models, and processing strategies of satellite gravimetry, the resolution, accuracy, and long-term consistency of mass transport series from satellite gravimetry can be significantly increased; and only in that case the potential of future technological sensor developments can be fully exploited.In continuation of the first RU phase, the individual project Improved Stochastic Modeling in GRACE/GRACE-FO Real Data Processing (ISTORE) 2 works towards the completion of the optimized stochastic modeling for GRACE and GRACE-FO gravity field determination. This includes three main tasks: (1) the extension of the stochastic instrument error models of the first phase including Global Navigation Satellite System (GNSS) observations; (2) the optimization of the combination of the different observations; (3) the inclusion of tidal and temporally changing non-tidal background model error variance-covariance matrices in the adjustment process. Based on recent developments within NEROGRAV we also work on a future GFZ GRACE/GRACE-FO Level-2 data release.This presentation provides an overview of the main outcomes of the advanced processing strategies. We will discuss details on stochastic modeling of accelerometer, inter-satellite ranging, and GNSS observations, and of ocean tide and non-tidal atmospheric-oceanic background models. We present Level-2 results based on the monthly solutions within the three test years 2007, 2014 and 2019 in the spectral and spatial domain in comparison with the standard GFZ GRACE/GRACE-FO RL06 time series.

To achieve an accurate and robust localization result under Global Navigation Satellite System (GNSS) denied environment, coupling GNSS with other sensors tends to be a promising solution, also known as the sensor fusion system. One of the most challenging scenarios in sensor fusion is to couple vision navigation system with GNSS, especially the tightly-coupled fusion. Compared with those loosely coupled ones, it performs the information fusion on a more intrinsic level, which can benefit the accuracy and robustness of the system. However, the major barrier lying in front of tightly coupling method is the way to unify the heterogeneous and asynchronous GNSS observations with vision observations. To solve this problem, we start from the Vison Simultaneous Localization and Mapping (V-SLAM) method, which has already become the mainstream in vision navigation domain. By utilizing existing attempts on graph-optimization-based GNSS data processing, we propose a joint graph optimization combining vision-SLAM with GNSS pseudorange. It is a tightly coupled fusion of GNSS and stereo camera, implemented by considering GNSS satellites as vision feature points. Based on the framework of mature state-of-art pose-graph SLAM, a graph optimization is applied simultaneously on GNSS and vision positioning, to achieve a robust, accurate, yet versatile fusion system.

On-board cooperative spacecraft relative navigation fusing GNSS with vision

<p><strong>Abstract</strong> Global Navigation Satellite Systems (GNSS), which are fundamental to Positioning, Navigation, and Timing (PNT), play also a critical role for the determination of ITRF and the monitoring of global geodetic parameters, such as geocenter and Earth Rotation Parameters (ERP), thanks to the globally distributed network. The Low Earth Orbit (LEO) satellites with onboard GNSS receivers as moving stations can significantly improve the observation geometry of a ground tracking network, and thus enhance the solutions. Previous studies have demonstrated that integrated processing of ground GNSS stations and LEO satellites can improve the GNSS satellite orbits and the determination of the geocenter, especially for spare networks. In this study, we further investigate the contribution of having five LEO satellites (two GRACE-FO satellites and three SWARM satellites) to a rather robust ground GNSS network, i.e., 120 stations. We focus not only on the GNSS satellite orbits, and the geocenter and ERP, but also on some other parameters like ambiguities. We further investigate the impact of the LEO orbit modelling and the weighting strategies of the LEO observations in order to properly consider the fact that onboard GNSS observations can be modelled more precisely than that from ground stations due to the differences in station environment and atmosphere.</p>