Combined Global Navigation Satellite System Research Articles

ARCHITECTURAL AND CIRCUIT MODEL OF NAVIGATION DEVICE FOR ULTRA-LIGHT LAUNCH VEHICLE

The paper considered the architectural and circuitry model of a navigation device for an ultra-lightweight launch vehicle designed to improve the accuracy and reliability of positioning and orientation in flight conditions. The study utilized modern approaches to the integration of global navigation satellite system (GNSS) and inertial navigation system (INS) to achieve high accuracy and stability of the system. The key components of the device including UM980 module, LSM6DSV32X inertial measurement module, STM32 F7 series microcontroller, E220-400T30S radio module, LMR50410XDBVR converter based power supply and RS422 interface were described in detail. The UM980 module, including SAW and LNA filter blocks, demonstrated the ability to efficiently process incoming RF signals. The NebulasIV main function block combined GNSS BB, GNSS RF, interfaces and power management unit (PMU) modules, supporting multiple communication interfaces such as UART, SPI and I²C, as well as RTK and PPS signals. The study analyzed the integration and accuracy aspects of the module in navigation systems. The signal processing scheme of the inertial navigation system module was also presented, including digital filters, free-fall and activity/inactivity detection functions, and interfaces for data communication. User bias functions and filter modes are considered to ensure the accuracy and stability of the navigation device. The work lays the groundwork for further improvements in navigation technology, especially in the context of increasing availability and accuracy for small scale space launches.

Combined GNSS reflectometry–refractometry for automated and continuous in situ surface mass balance estimation on an Antarctic ice shelf

Abstract. Reliable in situ surface mass balance (SMB) estimates in polar regions are scarce due to limited spatial and temporal data availability. This study aims at deriving automated and continuous specific SMB time series for fast-moving parts of ice sheets and shelves (flow velocity > 10 m a−1) by developing a combined global navigation satellite system (GNSS) reflectometry and refractometry (GNSS-RR) method. In situ snow density, snow water equivalent (SWE), and snow deposition or erosion are estimated simultaneously as an average over an area of several square meters and independently on weather conditions. The combined GNSS-RR method is validated and investigated regarding its applicability to a moving, high-latitude ice shelf. A combined GNSS-RR system was therefore installed in November 2021 on the Ekström ice shelf (flow velocity ≈ 150 m a−1) in Dronning Maud Land, Antarctica. The reflected and refracted GNSS observations from the site are post-processed to obtain snow accumulation (deposition and erosion), SWE, and snow density estimates with a 15 min temporal resolution. The results of the first 16 months of data show a high level of agreement with manual and automated reference observations from the same site. Snow accumulation, SWE, and density are derived with uncertainties of around 9 cm, 40 kg m−2 a−1, and 72 kg m−3, respectively. This pilot study forms the basis for extending observational networks with GNSS-RR capabilities, particularly in polar regions. Regional climate models, local snow modeling, and extensive remote sensing data products will profit from calibration and validation based on such in situ time series, especially if many such sensors will be deployed over larger regional scales.

Investigation of determining the accuracy of spatial vectors by the satellite method in a real time mode

Abstract The study of determining the accuracy of spatial vectors by the global navigation satellite system (GNSS) in real time (RTK) was conducted. The possibility of construction of precision geodetic networks by the combined method of static and RTK GNSS measurements which correspond to the set accuracy and reach the maximum economic efficiency is investigated. A technique providing the densification of GNSS network and the use of two simultaneously operating GNSS receivers (Rover) is proposed. The research was carried out at the points of the GNSS network of Dnister Pumped Storage Power Plant (PSPP) (Ukraine). As a result of comparison of reference and measured elements of vectors, it was found that the average absolute error in determining the spatial distance of 14 vectors was 5.3 mm. Rejection of vectors with a closed horizon reduced the error to 2.1 mm. The vectors are most accurately determined from two satellite systems (GPS, GLONASS) and from a single base station. The recommended distance to the base station is within one kilometer. Increasing the accuracy by 75 % in determining the vector by the proposed method in RTK mode is also shown. As a result of a posteriori optimization of combined GNSS networks, high accuracy on the reliability of vectors determined by the method in RTK mode was confirmed. The technique can be used to construct precision networks, to carry out repeated measurements for the monitoring of large engineering structures with an open horizon.

Ground Deformation Modelling of the 2020 Mw6.9 Samos Earthquake (Greece) Based on InSAR and GNSS Data

Modelling of combined Global Navigation Satellite System (GNSS) and Interferometric Synthetic Aperture Radar (InSAR) data was performed to characterize the source of the Mw6.9 earthquake that occurred to the north of Samos Island (Aegean Sea) on 30 October 2020. Pre-seismic analysis revealed an NNE–SSW extensional regime with normal faults along an E–W direction. Co-seismic analysis showed opening of the epicentral region with horizontal and vertical displacements of ~350 mm and ~90 mm, respectively. Line-of-sight (LOS) interferometric vectors were geodetically corrected using the GNSS data and decomposed into E–W and vertical displacement components. Compiled interferometric maps reveal that relatively large ground displacements had occurred in the western part of Samos but had attenuated towards the eastern and southern parts. Alternating motions occurred along and across the main geotectonic units of the island. The best-fit fault model has a two-segment listric fault plane (average slip 1.76 m) of normal type that lies adjacent to the northern coastline of Samos. This fault plane is 35 km long, extends to 15 km depth, and dips to the north at 60° and 40° angles for the upper and lower parts, respectively. A predominant dip-slip component and a substantial lateral one were modelled.

Intrinsic Identification and Mitigation of Multipath for Enhanced GNSS Positioning.

In global navigation satellite system (GNSS)-based positioning and applications, multipath is by far the most obstinate impact. To overcome paradoxical issues faced by current processing approaches for multipath, this paper employs an intrinsic method to identify and mitigate multipath based on empirical mode decomposition (EMD) and Hilbert–Huang transform (HHT). Frequency spectrum and power spectrum are comprehensively employed to identify and extract multipath from complex data series composed by combined GNSS observations. To systematically inspect the multipath from both code range and carrier phase, typical kinds of combinations of the GNSS observations including the code minus phase (CMP), differential correction (DC), and double differential (DD) carrier phase are selected for the suggested intrinsic approach to recognize and mitigate multipath under typical positioning modes. Compared with other current processing algorithms, the proposed methodology can deal with multipath under normal positioning modes without recourse to the conditions that satellite orbits are accurately repeated and surrounding environments of observing sites remain intact. The method can adaptively extract and eliminate multipath from solely the GNSS observations using intrinsic decomposition mechanism. From theoretical discussions and validating tests, it is found that both code and carrier phase multipath can be identified and distinguished from ionospheric delay and other impacts using the EMD based techniques. The resultant positioning accuracy is therefore improved to an obvious extent after the removal of the multipath. Overall, the proposed method can form an extensive and sound technical frame to enhance localization accuracy under typical GNSS positioning modes and harsh multipath environments.

Assessing the Performance of Multi-GNSS PPP-RTK in the Local Area

The benefits of an increased number of global navigation satellite systems (GNSS) in space have been confirmed for the robustness and convergence time of standard precise point positioning (PPP) solutions, as well as improved accuracy when (most of) the ambiguities are fixed. Yet, it is still worthwhile to investigate fast and high-precision GNSS parameter estimation to meet user needs. This contribution focuses on integer ambiguity resolution-enabled Precise Point Positioning (PPP-RTK) in the use of the observations from four global navigation systems, i.e., GPS (Global Positioning System), Galileo (European Global Navigation Satellite System), BDS (Chinese BeiDou Navigation Satellite System), and GLONASS (Global’naya Navigatsionnaya Sputnikova Sistema). An undifferenced and uncombined PPP-RTK model is implemented for which the satellite clock and phase bias corrections are computed from the data processing of a group of stations in a network and then provided to users to help them achieve integer ambiguity resolution on a single receiver by calibrating the satellite phase biases. The dataset is recorded in a local area of the GNSS network of the Netherlands, in which 12 stations are regarded as the reference to generate the corresponding corrections and 21 as the users to assess the performance of the multi-GNSS PPP-RTK in both kinematic and static positioning mode. The results show that the root-mean-square (RMS) errors of the ambiguity float solutions can achieve the same accuracy level of the ambiguity fixed solutions after convergence. The combined GNSS cases, on the contrary, reduce the horizontal RMS of GPS alone with 2 cm level to GPS + Galileo/GPS + Galileo + BDS/GPS + Galileo + BDS + GLONASS with 1 cm level. The convergence time benefits from both multi-GNSS and fixing ambiguities, and the performances of the ambiguity fixed solution are comparable to those of the multi-GNSS ambiguity float solutions. For instance, the convergence time of GPS alone ambiguity fixed solutions to achieve 10 cm three-dimensional (3D) positioning accuracy is 39.5 min, while it is 37 min for GPS + Galileo ambiguity float solutions; moreover, with the same criterion, the convergence time of GE ambiguity fixed solutions is 19 min, which is better than GPS + Galileo + BDS + GLONASS ambiguity float solutions with 28.5 min. The experiments indicate that GPS alone occasionally suffers from a wrong fixing problem; however, this problem does not exist in the combined systems. Finally, integer ambiguity resolution is still necessary for multi-GNSS in the case of fast achieving very-high-accuracy positioning, e.g., sub-centimeter level.

A Sensitivity Study of POD Using Dual-Frequency GPS for CubeSats Data Limitation and Resources

Making use of dual-frequency (DF) global navigation satellite system (GNSS) observations and good dynamic models, the precise orbit determination (POD) for the satellites on low earth orbits has been intensively investigated in the last decades and has achieved an accuracy of centimeters. With the rapidly increasing number of the CubeSat missions in recent years, the POD of CubeSats were also attempted with combined dynamic models and GNSS DF observations. While comprehensive dynamic models are allowed to be used in the postprocessing mode, strong constraints on the data completeness, continuity, and restricted resources due to the power and size limits of CubeSats still hamper the high-accuracy POD. An analysis of these constraints and their impact on the achievable orbital accuracy thus needs to be considered in the planning phase. In this study, with the focus put on the use of DF GNSS data in postprocessing CubeSat POD, a detailed sensitivity analysis of the orbital accuracy was performed w.r.t. the data continuity, completeness, observation sampling interval, latency requirements, availability of the attitude information, and arc length. It is found that the overlapping of several constraints often causes a relatively large degradation in the orbital accuracy, especially when one of the constraints is related to a low duty-cycle of, e.g., below 40% of time. Assuming that the GNSS data is properly tracked except for the assumed constraints, and using the International GNSS Service (IGS) final products or products from the IGS real-time service, the 3D orbital accuracy for arcs of 6 h to 24 h should generally be within or around 1 dm, provided that the limitation on data is not too severe, i.e., with a duty-cycle not lower than 40% and an observation sampling interval not larger than 60 s.

Influence of station density and multi-constellation GNSS observations on troposphere tomography

Abstract. Troposphere tomography, using multi-constellation observations from global navigation satellite systems (GNSSs), has become a novel approach for the three-dimensional (3-D) reconstruction of water vapour fields. An analysis of the integration of four GNSSs (BeiDou, GPS, GLONASS, and Galileo) observations is presented to investigate the impact of station density and single- and multi-constellation GNSS observations on troposphere tomography. Additionally, the optimal horizontal resolution of the research area is determined in Hong Kong considering both the number of voxels divided, and the coverage rate of discretized voxels penetrated by satellite signals. The results show that densification of the GNSS network plays a more important role than using multi-constellation GNSS observations in improving the retrieval of 3-D atmospheric water vapour profiles. The root mean square of slant wet delay (SWD) residuals derived from the single-GNSS observations decreased by 16 % when the data from the other four stations are added. Furthermore, additional experiments have been carried out to analyse the contributions of different combined GNSS data to the reconstructed results, and the comparisons show some interesting results: (1) the number of iterations used in determining the weighting matrices of different equations in tomography modelling can be decreased when considering multi-constellation GNSS observations and (2) the reconstructed quality of 3-D atmospheric water vapour using multi-constellation GNSS data can be improved by about 11 % when compared to the SWD estimated with precise point positioning, but this was not as high as expected.

Positioning stability improvement with inter-system biases on multi-GNSS PPP

Abstract The availability of multiple signals from different Global Navigation Satellite System (GNSS) constellations provides opportunities for improving positioning accuracy and initial convergence time. With dual-frequency observations from the four constellations (GPS, GLONASS, Galileo, and BeiDou), it is possible to investigate combined GNSS precise point positioning (PPP) accuracy and stability. The differences between GNSS systems result in inter-system biases (ISBs). We consider several ISB values such as GPS-GLONASS, GPS-Galileo, and GPS-BeiDou. These biases are compliant with key parameters defined in the multi-GNSS PPP processing. In this study, we present a unified PPP method that sets ISB values as fixed or constant. A comprehensive analysis that includes satellite visibility, position dilution of precision, position accuracy is performed to evaluate a unified PPP method with constrained cut-off elevation angles. Compared to the conventional PPP solutions, our approach shows more stable positioning at a constrained cut-off elevation angle of 50 degrees.

GNSS related periodic signals in coordinate time-series from Precise Point Positioning

In Global Navigation Satellite System (GNSS) coordinate time-series unrecognized errors and unmodelled (periodic) effects may bias nonlinear motions induced by geophysical signals. Hence, understanding and mitigating these errors is vital to reducing biases and on revealing subtle geophysical signals. To assess the nature of periodic signals in coordinate time-series Precise Point Positioning (PPP) solutions for the period 2008–2015 are generated. The solutions consider Global Positioning System (GPS), GLObalnaya NAvigatsionnaya Sputnikovaya Sistema (GLONASS) or combined GPS+GLONASS (GNSS) observations. We assess the periodic signals of station coordinates computed using the combined International GNSS Service (IGS) and four of its Analysis Centers (ACs) products. Furthermore, we make use of different filtering methods to investigate the sources of the periodic signals. A faint fortnightly signal in our PPP solution based on Jet Propulsion Laboratory (JPL) products and the existence of an 8 d period for those ACs generating combined GPS+GLONASS products are the main features in the GPS-only solutions. The existence of the 8 d period in the GPS-only solution indicates that GPS orbits computed in a combined GNSS solution contain GLONASS-specific signals. The GLONASS-only solution shows highly elevated powers at the third draconitic harmonic (∼120 d period), at the 8 d period and its harmonics (4 d, 2.67 d) besides the well-known annual, semi-annual and other draconitic harmonics. We show that the GLONASS constellation gaps before December 2011 contribute to the power at some of the frequencies. However, the well-known fortnightly signal in GPS-only solutions is not discernible in the GLONASS-only solution. The combined GNSS solution contains periodic signals from both systems, with most of the powers being reduced when compared to the single-GNSS solutions. A 52 per cent reduction for the horizontal components and a 36 per cent reduction for the vertical component are achieved for the fortnightly signal from the GNSS solution compared to the GPS-only solution. Comparing the results of the employed filtering methods reveals that the source of most of the powers of draconitic and fortnightly signals are satellite-induced with a non-zero contribution of site-specific errors.

Precise Point Positioning using Multi-Constellation GNSS Observations for Kinematic Applications

AbstractTraditional precise point positioning (PPP) is commonly based on un-differenced ionosphere-free linear combination of Global Positioning System (GPS) observations. Unfortunately, for kinematic applications, GPS often experiences poor satellite visibility or weak satellite geometry in urban areas. To overcome this limitation, we developed a PPP model, which combines the observations of three global navigation satellite systems (GNSS), namely GPS, GLONASS and Galileo. Both un-differenced and between-satellite single-difference (BSSD) ionosphere-free linear combinations of pseudorange and carrier phase GNSS measurements are processed. The performance of the combined GNSS PPP solution is compared with the GPS-only PPP solution using a real test scenario in downtown Kingston, Ontario. Inter-system biases between GPS and the other two systems are also studied and obtained as a byproduct of the PPP solution. It is shown that the addition of GLONASS observations improves the kinematic PPP solution accuracy in comparison with that of GPS-only solution. However, the contribution of adding Galileo observations is not significant due to the limited number of Galileo satellites launched up to date. In addition, BSSD solution is found to be superior to that of traditional un-differenced model.

Galileo orbit determination using combined GNSS and SLR observations

The first two Galileo In-Orbit Validation satellites were launched in October 2011 and started continuous signal transmission on all frequencies in early 2012. Both satellites are equipped with two different types of clocks, namely rubidium clocks and hydrogen masers. Based on two test periods, the quality of the Galileo orbit determination based on Global Navigation Satellite System (GNSS) and Satellite Laser Ranging (SLR) observations is assessed. The estimated satellite clock parameters are used as quality indicator for the orbits: A bump at orbital periods in the Allan deviation indicates systematic errors in the GNSS-only orbit determination. These errors almost vanish if SLR observations are considered in addition. As the internal consistency is degraded by the combination, the offset of the SLR reflector is shifted by +5 cm, resulting in an improved orbit consistency as well as accuracy. Another approach to reduce the systematic errors of the GNSS-only orbit determination employs constraints for the clock estimates with respect to a linear model. In general, one decimeter orbit accuracy could be achieved.

Instantaneous GPS–Galileo Attitude Determination: Single-Frequency Performance in Satellite-Deprived Environments

New and modernized global navigation satellite systems (GNSSs) are paving the way for an increasing number of applications in positioning, navigation, and timing (PNT). A combined GNSS constellation will significantly increase the number of visible satellites and, thus, will improve the geometry of observed satellites, enabling improvements in navigation solution availability, reliability, and accuracy. In this paper, a global positioning system (GPS) +Galileo robustness analysis is carried out for instantaneous single-frequency GNSS attitude determination. Precise attitude determination using multiple GNSS antennas mounted on a platform relies on successful resolution of the integer carrier-phase ambiguities. The multivariate-constrained least squares ambiguity decorrelation adjustment (MC-LAMBDA) method has been developed to resolve the integer ambiguities of the nonlinearly constrained GNSS attitude model that incorporates the known antenna geometry. In this paper, the method is used to analyze the attitude determination performance of a combined GPS +Galileo system. Special attention is thereby given to the GPS and Galileo intersystem biases (ISBs). The attitude determination performance is evaluated using GPS/Galileo data sets from a hardware-in-the-loop experiment and two real-data campaigns. In the hardware-in-the-loop experiment, a full GPS/Galileo constellation is simulated, and performance analyses are carried out under various satellite-deprived environments, such as urban canyons, open pits, and other satellite outages. In the first real-data experiment, single-frequency GPS data, combined with the data of Galileo in-orbit validation element (GIOVE) satellites GIOVE-A/GIOVE-B (the two experimental Galileo satellites), are used to analyze the two constellation attitude solutions. In the second real-data experiment, we present the results based on single-frequency data from one of the Galileo IOV satellites, combined with the data of GIOVE-A and GPS. We demonstrate and quantify the improved availability, reliability, and accuracy of attitude determination using the combined constellation.

Correction of static pressure on a research aircraft in accelerated flight using differential pressure measurements

Abstract. A method is described that estimates the error in the static pressure measurement on an aircraft from differential pressure measurements on the hemispherical surface of a Rosemount model 858AJ air velocity probe mounted on a boom ahead of the aircraft. The theoretical predictions for how the pressure should vary over the surface of the hemisphere, involving an unknown sensitivity parameter, leads to a set of equations that can be solved for the unknowns – angle of attack, angle of sideslip, dynamic pressure and the error in static pressure – if the sensitivity factor can be determined. The sensitivity factor was determined on the University of Wyoming King Air research aircraft by comparisons with the error measured with a carefully designed sonde towed on connecting tubing behind the aircraft – a trailing cone – and the result was shown to have a precision of about ±10 Pa over a wide range of conditions, including various altitudes, power settings, and gear and flap extensions. Under accelerated flight conditions, geometric altitude data from a combined Global Navigation Satellite System (GNSS) and inertial measurement unit (IMU) system are used to estimate acceleration effects on the error, and the algorithm is shown to predict corrections to a precision of better than ±20 Pa under those conditions. Some limiting factors affecting the precision of static pressure measurement on a research aircraft are discussed.

Enhancing mobile applications by combination of GNSS and W-CDMA

The location requirements for emergency callers outside urban areas can hardly be fulfilled without global navigation satellite systems (GNSS). Consequently, interest in positioning techniques based on use of a GNSS such as GPS or on the cellular network infrastructure itself is growing rapidly in the mobile-telephone community. Moreover, the increasing demand for commercial location-based services (LBS) has driven cellular-phone and network manufacturers to focus on positioning solutions which are even more accurate than the regulatory mandates for positioning of emergency callers. One example of these upcoming LBS is our PARAMOUNT project, which aims at improving user-friendly info-mobility services for hikers and mountaineers by combining wireless communications (GMTS), satellite navigation (GNSS) and geographic information systems (GIS), based on a mobile client/server architecture. The availability of mobile phones or PDAs with combined GNSS and cellular network-based wireless communication on a high integration level is one primary demand of such LBS applications. Based on this, we will give some initial answers to the question of whether mobile handset architecture synergies exist for the combination of GNSS with wireless location in CDMA cellular wireless networks. In order to identify synergies, we will outline similarities and differences between wireless communication and satellite navigation. In this respect, we pay particular attention to the so-called RAKE receiver architecture employed in mobile CDMA cellular handsets. Our initial investigations will show that the RAKE receiver architecture, on which mobile CDMA cellular handsets are based, will most likely be the one most suitable for achieving synergies between the two positioning techniques within the same mobile handset architecture. Consequently, several receiver components could be used to handle both types of signals (navigation and communications), resulting in a reduction of manufacturing costs and in a decrease in energy consumption.