Next Article in Journal
A Survey of Energy-Efficient Communication Protocols with QoS Guarantees in Wireless Multimedia Sensor Networks
Next Article in Special Issue
Research on the Algorithm Model for Measuring Ocean Waves Based on Satellite GPS Signals in China
Previous Article in Journal
Faster R-CNN and Geometric Transformation-Based Detection of Driver’s Eyes Using Multiple Near-Infrared Camera Sensors
Previous Article in Special Issue
The Marine Boundary Layer Height over the Western North Pacific Based on GPS Radio Occultation, Island Soundings, and Numerical Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

BeiDou Augmented Navigation from Low Earth Orbit Satellites

1
GNSS Research Center, Wuhan University, 129 Luoyu Road, Wuhan 430079, China
2
College of Geomatics, Shandong University of Science and Technology, No. 579 Qianwangang Road, Qingdao 266590, China
3
Beijing Institute of Tracking and Telecommunications Technology, No. 26 Beiqing Road, Beijing 100094, China
4
Collaborative Innovation Center of Geospatial Technology, Wuhan University, 129 Luoyu Road, Wuhan 430079, China
*
Authors to whom correspondence should be addressed.
Sensors 2019, 19(1), 198; https://doi.org/10.3390/s19010198
Submission received: 16 November 2018 / Revised: 20 December 2018 / Accepted: 26 December 2018 / Published: 7 January 2019
(This article belongs to the Special Issue High-Precision GNSS in Remote Sensing Applications)

Abstract

:
Currently, the Global Navigation Satellite System (GNSS) mainly uses the satellites in Medium Earth Orbit (MEO) to provide position, navigation, and timing (PNT) service. The weak navigation signals limit its usage in deep attenuation environments, and make it easy to interference and counterfeit by jammers or spoofers. Moreover, being far away to the Earth results in relatively slow motion of the satellites in the sky and geometric change, making long time needed for achieved centimeter positioning accuracy. By using the satellites in Lower Earth Orbit (LEO) as the navigation satellites, these disadvantages can be addressed. In this contribution, the advantages of navigation from LEO constellation has been investigated and analyzed theoretically. The space segment of global Chinese BeiDou Navigation Satellite System consisting of three GEO, three IGSO, and 24 MEO satellites has been simulated with a LEO constellation with 120 satellites in 10 orbit planes with inclination of 55 degrees in a nearly circular orbit (eccentricity about 0.000001) at an approximate altitude of 975 km. With simulated data, the performance of LEO constellation to augment the global Chinese BeiDou Navigation Satellite System (BeiDou-3) has been assessed, as one of the example to show the promising of using LEO as navigation system. The results demonstrate that the satellite visibility and position dilution of precision have been significantly improved, particularly in mid-latitude region of Asia-Pacific region, once the LEO data were combined with BeiDou-3 for navigation. Most importantly, the convergence time for Precise Point Positioning (PPP) can be shorted from about 30 min to 1 min, which is essential and promising for real-time PPP application. Considering there are a plenty of commercial LEO communication constellation with hundreds or thousands of satellites, navigation from LEO will be an economic and promising way to change the heavily relay on GNSS systems.

1. Introduction

The Global Navigation Satellite System (GNSS) represents a constellation of satellites providing signals from space, transmitting positioning and timing data with global coverage [1]. A GNSS receiver employs trilateration to determine its position on or near the Earth’s surface by timing signals from four or more GNSS satellites. Currently, GNSS is well recognized as the major enabler of ‘precision’, and widely used for the precise timing and positioning. There are two fully operational GNSS systems at present: the United States’ Global Positioning System (GPS), and the Russian Federation’s Global Navigation Satellite System (GLONASS). Two systems under development, the Chinese Beidou Navigation Satellite System (BeiDou), and the European Union’s Galileo system, which are both expected to achieve full global coverage capability by 2020. In addition, the Japanese Quasi-Zenith Satellite System (QZSS) and Indian Navigation with Indian Constellation (NAVIC) are two regional systems to provide the augmentation service. The Multi-GNSS constellations provide the orbit as well as frequency diversity to strengthen the observation geometry and reduce the transmitting errors, resulting in improvement of positioning accuracy and stability [2].
The above-mentioned GNSS systems are consisted of the satellites at Geostationary orbit (GEO), Inclined Geostationary orbit (IGSO) or Medium Earth Orbit (MEO), and all of them are 20×~30× further from the Earth compared to the LEO. Hence, the energy of navigation signals is significantly spread out during their passage to the Earth, resulting in too weak signals to be used reliably in deep attenuation environments such as indoors or in city canyons. Moreover, the GNSS signals are susceptible to interference and counterfeit by jammers or spoofers [3]. From the view of high accuracy positioning with carrier-phase measurements, being far away to the Earth results in relative slow motion of the satellites in the sky and geometric change with respect to the ground stations. Hence, a quite long time (normally up to more than 30 min) is needed for the ambiguities to be reliably converged to the correct integer values [2]. By using the satellites at Lower Earth Orbiter (LEO) as the navigation satellites, these disadvantages can be improved [4,5,6].
LEO satellites (with altitude from 160 km to 2000 km) are 20 times closer to the Earth than GNSS at MEO, leading to stronger radio signals on the ground for navigation because of the lower path loss of the signal. The swift motion of LEO provides the rapid geometric change for rapid integer ambiguity resolution. The other advantages of LEO used as navigation satellites include whitening multipath and effective Doppler positioning, as the rapid motion of LEO makes the reflected signals are no longer effectively static over short averaging times and strengthens the Doppler shift [7]. However, being them too close to the Earth, the ground converge of LEO satellite is less than one-tenth of that of MEO. Hence, the hundreds of LEO satellites are needed to match the coverage of GNSS [3]. It is extremely expensive to deploy such a massive number of LEO satellites to only provide navigation service. Thanks to the development of commercial communication LEO constellation, it is promising and possible to integrate the navigation components to the communication satellites. The summary of some deployed or proposed commercial communication LEO constellations is given in Table 1.
The Iridium satellite constellation was proposed in 1987 and became operational in 1998 to provide voice and data coverage to satellite phones. The constellation consists of 66 active satellites and additional backup satellites. The operational satellites are deployed in six orbital planes spaced 30° apart, with 11 satellites in each plane. The orbit altitude and inclination are approximately 781 km and 86.4°, respectively [8]. The second generation Iridium-NEXT satellites began to be deployed into the existing constellation in January 2017. Somewhat similar to the Iridium satellite constellation, Globalstar is also a LEO satellite constellation for satellite phone and low-speed data communications, the first generation Globalstar constellation have 48 satellites with an orbital height of approximately 1400 km and an inclination of 52°. A Globalstar second-generation constellation consists of 24 LEO satellites, which were fully launched in 2013.
To keep up with the rising demand for broadband, several new LEO constellations are being proposed. Proposed in 2013, LeoSat is designed to consist of 108 satellites in 1400 km orbit and use Ka-band for provide secure broadband links [9]. The Canadian Telesat proposed a LEO communication system consists of a constellation of 117 satellites in 11 orbital planes, with six planes (12 satellites per plane) inclined 99.5 degrees in a circular orbit at an approximate altitude of 1000 km and five planes (nine satellites per plane) inclined 37.4 degrees in a circular orbit at an approximate altitude of 1248 km [10]. In January 2015, OneWeb announced a partnership with Virgin and Qualcomm to produce a satellite internet constellation of approximately 720 satellites expected to provide global Internet broadband service to individual consumers as early as 2019. These satellites evenly distributed in 18 near-polar orbital planes, at an approximate altitude of 1200 km. Moreover, the system will be expended to accommodate additional 1260, which would double the number of orbital planes from 18 to 36, and increase the maximum number of satellites per plane from 40 to 55 [11]. SpaceX, with support from Google, announced a similar ambition for StarLink constellation with 4425 LEO satellites operating in 83 orbital planes at altitudes ranging from 1100 km to 1325 km in March 2018, and this new approval gives the green light for the planned 7518 satellites in the Starlink constellation [12,13]. The prototype test-flight satellites were successfully launched on 22 February 2018. In August of 2015, Samsung expressed interest with a proposal for a LEO constellation of 4600 [14]. Boeing joined the race in June of 2016 announcing plans for a LEO constellation of 2956 satellites to provide broadband internet access globally. Boeing plans to deploy the first portion of the system, 1396 satellites operating at an altitude of 1200 km, and to subsequently increase the constellation to the total of 2956 satellites [15]. Moreover, China Aerospace Science and Industry Corporation (CASIC) have launched a project named ‘Hongyun’ to build a global wideband communications LEO satellite constellation consisting of 156 satellites with altitude 1000 km. The first phase of the project is expected to complete by 2021, and the first experimental satellite has been launched in 2017 [16]. Meanwhile, in late-February 2018, the China Aerospace Science and Technology Corporation (CASC) announced plans to build the ‘Hongyan’ LEO constellation of 324 small satellites for global communications and other services. The Hongyan constellation with 1100 km altitude is targeted to be operational by 2021 [17].
Currently, the Iridium constellation is already broadcasting structured independent navigation signals, and provides independent positioning, navigation and timing (PNT) service based on Satelles’ satellite time and location (STL) signal technology to augment or serve as a back-up to existing GNSS PNT solutions [18]. STL is designed such that a receiver can reliably decode the bursts and perform precise Doppler and range measurements at attenuation of up to 39 dB relative to unobstructed reception. This is sufficient to penetrate buildings and other occlusions, providing coverage in most deep indoor and urban canyon environments. The testing results show that STL always experienced strong signals and can provide sub-microsecond timekeeping even in deep attenuation environments such as indoors. In addition, after 10 min of convergence, the STL solution had converged to an accuracy of better than 35 m with same accuracy in vertical and horizontal direction [19].
Used as the independent navigation satellite system, a PNT service comparable to GPS could be achieved although with the added benefits of LEO, including stronger signals and rapid changes in geometry. Besides, the LEO can also be used to augment the GNSS to provide combined solution with GNSS and LEO navigation signals. The aim of this study is to assess the possibility and performance of LEO constellation used to augment the BeiDou navigation system for precise positioning. In Section 2, the advantage of LEO navigation system in positioning accuracy and integer ambiguity resolution will be analyzed theoretically. In Section 3, the strategy used for simulating the observations from GNSS and LEO constellation as well as Precise Point Positioning (PPP) with GPS-only or combined GNSS and LEO will be presented. The results will be analyzed and discussed in Section 4 with emphasizing on the convergence of ambiguity as well as accuracy after convergence. Finally, the study will be concluded in Section 5.

2. Advantages of Navigation with LEO Constellation

2.1. Positioning Error

Like GNSS systems, the users will also use the trilateration to compute their position with navigation service provided by LEO constellation. Hence, statistically, the point positioning accuracy can be described by the product of the user equivalent range error (UERE), which combines all measurement and modeling errors of an individual pseudo-range, and the dilution of precision (DOP), which maps such errors to the position uncertainty for a given number and geometric distribution of observed navigation satellites [20]. The relation can be expressed as following.
σ 3 D = P D O P σ U E R E
where σ 3 D is the user 3D root-mean-square (RMS) positioning error, PDOP is the position DOP, and σ U E R E is the UERE.
The DOP factors are derived from the terms associated with the user satellite geometry. Therefore, only geometry maps range errors into user space. This means that DOP is the geometry factor of the visual satellites related to their distribution. Whether the satellites are distributed more homogeneously, it will be helpful to the position estimation for the observation information gathered from more directions. Usually, the DOP value is resolved into horizontal (HDOP) and vertical dilution of precision (VDOP), these give us the relationship between UERE and errors in the vertical and horizontal directions. Thanks to the fast motion and numerous satellites, the DOP will be significantly improved by LEO constellation. The simulated results will be presented in Section 4.1.
The UERE is further divided into a signal-in-space [user] range error (SIS URE or SISRE) and a user equipment error (UEE). While the UEE term comprises all receiver-related contributions such as noise, multipath and uncorrected atmospheric errors, the SISRE describes the statistical uncertainty of the modeled pseudo-range due to errors in the broadcast orbit and clock information. It is driven by the space segment characteristics (e.g., clock stability and predictability of orbital motion), as well as the control segment capabilities (orbit and clock determination performance, distribution of monitoring stations and upload capacity). The orientation of the line-of-sight vector impacts the orbit errors Δ r = ( Δ r R , Δ r A , Δ r C ) in the radial ( Δ r R ), along-track ( Δ r A ) and cross-track ( Δ r C ) directions on the line-of-sight range, and hence the user location. The average contribution over all points of the Earth within the visibility cone of the satellite is considered normally for a statistical description. The orbit-only contribution to the SISRE can then be described as a weighted average:
S I S R E ( o r b ) = R M S [ w R 2 Δ r R 2 + w A , C 2 ( Δ r A 2 + Δ r C 2 ) ] = w R 2 R 2 + w A , C 2 ( A 2 + C 2 )
of the RMS errors R = R M S ( Δ r R ) , A = R M S ( Δ r A ) and C = R M S ( Δ r C ) scale like text along the three axes. The altitude of the GNSS satellite determines the weight factors w R and w A , C [21]. For MEO, the w R is closed to 1, and it indicates that the radial orbits largely project onto the user. However, for LEO orbits, the weight factors become almost equally important to the SISRE [3], as their cones should be wider to provide larger converge. The combined orbit and clock SISRE is obtained from the expression:
S I S R E = ( w R Δ r R Δ c d t ) 2 + w A , C ( A 2 + C 2 )
the Δ c d t represents the error of the broadcast clock offset. It takes into consideration that the modeled pseudo-range is derived from the difference of the radial orbit error and the clock offset. The analysis performed by Montenbruck et al. [21] indicates that the SISRE for GPS, GLONASS, Galileo, QZSS, and BeiDou regional system (BeiDou-2) are 0.71 m, 1.97 m, 1.64 m, 0.57 m, and 1.46 m, respectively. For the UEE term, although it varies with different GNSS system and correction model used, the uncertainty is about 1.2 m, 0.2 m, and 1.0 m for ionosphere and troposphere modeling erros as well as Equipment Receiver Noise & Multipath. In total, the UERE is about 1.77 m, 2.51 m, 2.01 m, 1.66 m, and 1.96 m for GPS, GLONASS, Galileo, QZSS, and BeiDou-2. For the LEO navigation constellation, as the same errors are exist in receiver side as well as when the signals pass through the atmosphere, the UEE for LEO constellation may have similar magnitude as that of GNSS. Hence, the clock and orbit accuracy of LEO constellation are essential for SISRE of the LEO navigation constellation. For LEO navigation constellation, three potential data can be used for orbit and clock determination, e.g., ground tracking data, intersat link data, and onboard GNSS tracking data. With ground tracking or intersat link data, as GNSS satellites, the accuracy of orbit and clock are rely on the predicted one. However, as closer to the satellites, the perturbation acting on the LEO are hard to model, resulting in poor prediction ability for longer orbit. Hence, the frequency updating of the broadcast ephemerics of LEO is needed. As to the analysis by Reid [22], the 3D RMS of orbit error can reach to about 1.5 m. Once the onboard GNSS data available, the real-time onboard solution can be obtained with accuracy about 1.5 m in 3D. For the clock, the worst uncertainty is about 6.0 ns (approximate 1.8 m), that because the LEO use the chip-scale atomic clocks (CSAC) for the restriction of power consumption, size, weight, and cost etc.

2.2. Integer Ambiguity Resolution

For high accuracy positioning using the carrier-phase measurements (e.g., PPP), the positioning accuracy is rely on if the ambiguities can be correctly fixed to be integer. Hence, the convergence speed of ambiguity solution has essential impact on the positioning solution. Normally ambiguities are solved at the same time with the other parameters such as station coordinates, receiver clocks, troposphereic delay, and so on. With the change of observation geometry, the multipath and atmospheric delay is estimated quickly, making that the ambiguities can be separated with those parameters. In this case, a reliable ambiguity solution can be obtained. The probability of ambiguity solution is dependent on the satellite types and altitude and the empirical formula can be expressed as [23]:
P ( a m b i g u i t y ) = f ( s ˙ , P D O P ) = κ s ˙ P D O P
where κ coefficient related to frequency, atmosphere and so on; P probability; s ˙ change rate of line-of-sight between tracking station and satellite. The ambiguity solution related to types of satellite orbits using GPS L1 frequency as an example is approximately shown in Table 2.
The results demonstrated that the probability of ambiguity solution is dependent on the PDOP and the change rate of line-of-sight between tracking station and satellite. For the satellite at lower orbit, the swift motion makes that the rate of line-of-sight between tracking station and satellite change rapidly, and improve the geometry and PDOP. In this case, the ambiguity can be converged to the integer in quite short period. However, for higher altitude, longer time is needed for integer ambiguity resolution. In particular, for GEO, the ambiguity cannot be solved if there are no multi-frequency observations.

3. Strategy for Data Simulation and Procession

In order to assess the independent and integrated navigation performance of LEO constellation, the ground observations for BeiDou global (BeiDou-3) constellation and a LEO constellation were simulated. According to the Interface Control Document [24], the BeiDou-3 constellation will consist of three GEO, three IGSO, and 24 MEO satellites. The three GEO satellites will operate in the orbit at an altitude of 35,786 km, and are located at 80°E, 110.5°E, and 140°E, respectively. The IGSO satellites with same altitude as GEO have an inclination of the orbital planes of 55 degrees. The MEO satellites will operate in orbit with altitude of 21,528 km and inclination of 55 degrees [24]. The simulated LEO navigation system consists of a constellation of 120 satellites in 10 orbit planes with inclination of 55 degree in a nearly circular orbit (eccentricity about 0.000001) at an approximate altitude of 975 km. As the main aim of the proposed LEO constellation is to argument the BeiDou-3 system, in particularly providing the high accuracy positioning in China. With the initial orbit elements, the 24 h orbits for these BeiDou-3 and LEO satellites can be generated with numerical integration based on the force models summarized on Table 3. Figure 1 shows the distributions of BeiDou-3 and LEO constellation in the space.
Three ground stations (CENT located in Wuhan, latitude 30.52779, longitude 114.35686; POTS located in Potsdam, latitude 52.37929, longitude 13.06609; NTUS located in Singapore, latitude 1.34580, longitude 103.67995) are selected for simulating the ground tracking data to demonstrate the positioning performance of the BeiDou-3 and LEO constellation in different regions. Figure 2 shows the distribution. The observations are simulated based on the following observation equations of code and phase:
ρ r , i s = R r s + c d t r c d t s + m r s T r + f i 2 f 1 2 I r s + b r , i b i S + δ P C O , i + δ P C V , i + δ r e l + δ m u l , ρ r , i s + ε ρ r , i s φ r , i s = R r s + c d t r c d t s + m r s T r f i 2 f 1 2 I r s + λ i N r s + B r , i B i S + δ P C O , i + δ P C V , i + δ r e l + δ m u l , φ r , i s + δ w i n d u p + ε φ r , i s
where the subscript i denote the number of frequency with frequency value f i . The subscript r and the superscript s denote receiver and satellite, respectively. ρ r , i s and φ r , i s are code and phase range at frequency i which need to be simulated. R r s denotes satellite-to-receiver distance. d t r and d t s denote receiver and satellite clock. T r and I r s denote tropospheric and ionospheric delay. b r , i and b i S indicate the biases for code in the receiver and satellites sides at frequency i , while B r , i and B i S are the corresponding biases for phase. These biases terms are omitted in the simulation, resulting in no biases between different observations and system. The term λ i denotes wavelength and N r s means integer ambiguity. δ P C O , i and δ P C V , i are the phase centre offset (PCO) and phase centre variation (PCV) at frequency i , respectively. δ r e l and δ w i n d u p are the relative effect and phase wind-up. δ m u l , ρ r , i s and δ m u l , φ r , i s are the multipath error for code and phase measurement. The terms ε ρ r , i s and ε φ r , i s denote the code and phase observation errors, respectively.
The simulation is mainly to calculate all the components on the right side of observation equations. The satellite-to-receiver distance is computed with the receiver and satellite positions. Specifically, this distance is associated with the mass centers of the satellite and the receiver. Thus, this distance needed to corrected by the PCO and PCV of receiver and satellite. The receiver clock error of the ground tracking stations are with uncertainty about 10 us (3000 m), while the satellite clock are interpolated from the precise clock error file. The ionospheric delay is ignored as we use the ionosphere-free combination observations, and the second and third ionospheric delays are too small to consider in the data processing. For the tropospheric delay, we use the Saastamoinen model [27] and global mapping functions (GMF) [28] for the ground tracking network. Moreover, relativistic correction, phase windup correction, and tidal displacement are all taken into consideration. The ambiguities were set as zero. For GNSS satellites, the nominal yaw-steering attitude mode is used [21], which the LEO’s attitude is assumed along the along-track, cross-track, and radial direction. The noises for code and phase are set as 1 m and 2 mm. Both of the LEO and GNSS satellites used the same frequency for the simulation. Table 4 summarizes the configuration.
Once the observations are generated, PPP are performed for the station to investigate the accuracy and convergence time achieved by BeiDou-3 only and BeiDou-3/LEO combined solutions. In this case, the epoch-wise station coordinates, epoch-wide receiver clock, hourly troposphere zenith delays, and ambiguity are estimated. As no inter-system bias was simulated, this term were not be estimated when BeiDou-3 is combined with LEO for PPP [29]. The ionosphere-free linear combination was used, with 2 mm observational standard deviation applied in the weighting of both the BeiDou-3 and LEO phase observations, and 1 m for both the BeiDou-3 and LEO code observations. Once the 3D positioning accuracy within ten continuous epochs is less than 10 cm, the positioning is treated as convergence.

4. Results and Discussion

In this section, the performance of LEO constellation used to argument the BeiDou-3 is assessed with several metrics including satellite visibility, PDOP, accuracy and convergence of positioning.

4.1. Satellite Visibility

The satellite visibility of BeiDou-3 only as well as BeiDou-3 and LEO constellation is computed. With dividing the global to the grids with 5° in latitude and 5° longitude, the number of satellites tracked by the location in the center of grid is calculated with 0 mask elevation and 30 s data sampling interval. Afterword, the average number for 1 day is used as the indicator of satellite visibility. Figure 3 demonstrates the global distribution of the satellite visibility for BeiDou-3 only (Figure 3a) and BeiDou-3 as well as LEO constellation (Figure 3b). For BeiDou-3, in general, the satellite visibility is distributed symmetrically in Northern hemisphere and Southern hemispheres, but unevenly in the Eastern hemisphere and Western hemisphere. More than 12 satellites can be tracked in Asia-Pacific region, where the 3 GEO and 3 IGSO satellites are located. In particular, for the low latitude region in the Asia-Pacific region, up to 15 satellites can be tracked. However, the number of tracked satellites are lower in the western hemisphere, particularly for the mid-latitude region of north and South American, there are only 8 satellite can be tracked in averaged. For the region with latitude beyond 55°, about 10 satellites are in view, as the orbital inclination of BeiDou-3 MEO satellites are set as 55°.
Once the 120 LEO satellites have been incorporated, although the global distribution is still uneven in Western and Eastern hemisphere, the satellite visibility has been improved significantly, particularly in the mid-latitude region, as the low-orbit navigation satellites are designed to have an orbital inclination of 55°. In China, at least 18 satellites can be tracked, particularly it can be reach to more than 20 in the north part of China. The great satellite visibility can improve the accuracy obviously, and explain the reason for selection of such kind of LEO constellation configuration. Compared to that of BeiDou-3 only, the satellite visibility in the mid-latitude region of north and South American has been improved, and more then 14 satellites can be tracked. However, the contribution of LEO to the satellite visibility is limited, as the limitation of orbital inclination as well as ground coverage of LEO satellites. Table 5 lists the statistical results for the two configurations. In general, once the LEO satellites have been incorporated, the maximum number of visible satellites has increased from 15.5 to 21.6. Meanwhile, the minimum number of visible satellites has increased from 7.6 to 10.1. Most significant improvement can be observed by the averaged visible satellite, which is increased from 10.7 to 16.3, indicating that the LEO has greatly improved the performance of navigation services worldwide.

4.2. PDOP

The PDOP indicates the positioning accuracy theoretically, and is related to the number of satellite tracked as well as the distribution of satellites. In general, the smaller the PDOP value is, the better accuracy for positioning is. Similar as calculation of satellite visibility, the averaged PDOP with 24 h data for the location in the center of each 5° × 5° grid is obtained with with 0 mask elevation and 30 s data sampling interval. Figure 4 demonstrates the global distribution of the PDOP for BeiDou-3 only and BeiDou-3 (Figure 4a) as well as LEO constellation (Figure 4b). Similar as satellite visibility, the PDOP is distributed symmetrically in Northern hemisphere and Southern hemispheres, but unevenly in the Eastern hemisphere and Western hemisphere. For BeiDou-3, the PDOP value reaches a minimum of 1.2–1.3 in the low latitudes of Asia-Pacific region, where is the great tracking ability for the three GEO and three IGSO satellites. In most areas of China, the PDOP value is about 1.5. As less satellites can be tracked, the mid-latitude regions of the north and south Americas have the largest PDOP up to 2.2. The poor positioning ability (with PDOP up to 2.2) in the north and south pole region is related to the configuration of BeiDou-3 satellites, as above mentioned, the 55° orbit inclination of MEO satellites limits the satellite visibility in these regions, resulting in larger PDOP.
Once the LEO satellites have been used to augment the BeiDou-3, the navigation service performance in mid-latitude regions has been significantly improved. The PDOP values in the Asia-Pacific region were around 1, and low-latitude and mid-latitude show similar performance. For PDOP in the mid-latitude region of north and south Americas, it has been reduced to approximate 1.2 from abound 1.8. Thanks to the 55° orbital inclination used for LEO satellites, the PDOP has been improved to about 1.2 in the region with latitude about 55° and −55°. However, limited improvement is identified in low latitude region of north and south Americas as well as the pole regions. In statistics, as listed in Table 6, the maximum PDOP value was reduced from 2.24 to 1.93, the minimum PDOP value was decreased from 1.19 to 1.01, and the average PDOP value was reduced from 1.63 to 1.22. This again confirms that the LEO constellation can greatly improve the performance of BeiDou-3.

4.3. SPP

The standard point positioning (SPP) with pseudo-range measurements was carried to investigate the performance of incorporating LEO constellation to augment BeiDou-3 for the selected station. The simulated code measurements with 1 m random noise were processed with precise satellites’ orbit and clock for observation simulation. The ionosphere delay has been corrected with ionosphere-free combination, and the wet and dry troposphere delay was corrected with Saastamoinen model and global mapping functions (GMF). In this case, only four epoch-wise parameters including coordinates and receiver clock were estimated. Two solutions with BeiDou-3 only or combined BeiDou-3 and LEO measurements were performed.
In general, the inclusion of LEO data has limited effect on SPP. Taking CENT station for example, the positioning accuracy has been improved from 105.3 to 100.7 cm in earth direction, 99.1 to 93.2 in north direction, as well as 293.7 to 261.3 cm in Up direction, once LEO data have been used. Table 7 shows the SPP precision of these three stations in north, east and up direction.

4.4. PPP

In order to investigate the contribution of LEO to accuracy and convergence of high-accuracy positioning with carrier-phase measurements, two PPP experiments have been performed with BeiDou-3 or combined BeiDou-3 and LEO measurements. The simulated code and phase measurements with 1 m and 2 mm random noise were processed with precise satellites’ orbit and clock. The ionosphere delay has been corrected with ionosphere-free combination, and the dry troposphere delay was corrected with the Saastamoinen model and global mapping functions (GMF), and the wet delay has been estimated per hour. Besides, coordinates and receiver clock, the constant and inter-system bias were estimated.
Figure 5 shows the positioning error in north, east, and up direction for the two solutions of the selected station. In general, the inclusion of LEO can significantly accelerate the positioning convergence. The convergence time can be shorted to 1 min 9 s from 28 min 35 s. The results indicates the most great advantage of LEO constellation, which can greatly accelerating the convergence speed of the ambiguities, and providing promising way for real-time PPP application.
Figure 6 further show the variations of kinematic positioning errors of CENT station in 1 h after the convergence of PPP for BeiDou-3 only and BeiDou-3 as well as LEO combined solutions. It is easy to notice that the positioning accuracy is within 5 cm in horizontal direction, whereas the it is within 8 cm in vertical direction. However, once the LEO data used, the variations of kinematic positioning errors become stable, and it is below 5 cm in all of the three direction. Statically, the RMS of position accuracy has been improved from 1.12, 0.74, and 2.05 cm to 0.65, 0.59, 1.23 cm in east, north, and up direction, once the LEO satellites are used to augment the BeiDou-3. The simulation results of all these three stations are listed in Table 8. This clearly demonstrate the advantage of LEO used as the navigation system.

5. Conclusions

Currently, the Global Navigation Satellite System (GNSS) mainly use the satellites in Medium Earth Orbit (MEO) to provide position, navigation, and timing (PNT) service. The weak navigation signals limit its usage in deep attenuation environments, and make it easy to interference and counterfeit by jammers or spoofers. Moreover, being far away to the Earth results in relative slow motion of the satellites in the sky and geometric change, making long time needed for achieved centimeter positioning accuracy. By using the satellites at lower earth orbit (LEO) as the navigation satellites, these disadvantages can be addressed.
LEO satellites (with altitudes from 160 km to 2000 km) are 20 times closer to the Earth than GNSS at MEO, leading to stronger radio signals on the ground for navigation. The swift motion of LEO provides the rapid geometric change, which resulting quicker PDOP change and rapid convergence for integer ambiguity. We compared and analyze this theoretically. We have simulated the BeiDou-3 constellation consisting three GEO, three IGSO, and 24 MEO satellites and LEO navigation system consisting of a constellation of 120 satellites in 10 orbit planes with inclination of 55 degrees in a nearly circular orbit (eccentricity about 0.000001) at an approximate altitude of 975 km to assess the integrated navigation performance of LEO constellation. The results confirm the theoretical analysis, and demonstrate the satellite visibility and PDOP can be improved. More importantly, the accuracy and convergence of PPP can be improved and significantly reduced to about 1 min.
Although the simulation used in this study is used for augmenting BeiDou-3 with LEO constellation. The LEO constellation can also be used as independent navigation system, as Iridium does. However, as too close to the Earth, the ground converge of LEO satellites is less than one-tenth of that of MEO ones. Hence, hundreds of LEO satellites would be needed to match the coverage of GNSS. Thanks to the development of commercial communication LEO constellation, it is promising and possible to integrate the navigation components to the communication satellites to provide independent navigation from LEO. Currently, many plans have been proposed to achieve this around the world, such as Chinese Hongyun and Hongyan. Once it has been implemented, the real-time navigation with cm level accuracy will be achieved easily, and definitely significantly change our daily life.

Author Contributions

M.S. and X.S. wrote the main manuscript text. X.S. conducted the experiments. Q.Z. and J.L. conceived of the idea and reviewed the experimental methods and results. All authors reviewed the manuscript.

Funding

This work was supported by the Shandong Provincial Natural Science Foundation, China (Grant number ZR2018PD006); Scientific Research Foundation of Shandong University of Science and Technology for Recruited Talents (Grant number 2017RCJJ072).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Teunissen, P.J.G.; Montenbruck, O. Springer Handbook of Global Navigation Satellite Systems; Springer International Publishing Switzerland: Basel, Switzerland, 2017. [Google Scholar] [CrossRef]
  2. Li, X.; Ge, M.; Dai, X.; Ren, X.; Fritsche, M.; Wickert, J.; Schuh, H. Accuracy and reliability of multi-GNSS real-time precise positioning: GPS, GLONASS, BeiDou, and Galileo. J. Geod. 2015, 89, 607–635. [Google Scholar] [CrossRef] [Green Version]
  3. Reid, T.G.; Neish, A.M.; Walter, T.F.; Enge, P.K. Leveraging Commercial broadband LEO constellation for navigation. In Proceedings of the ION GNSS+ 2016, Portland, OR, USA, 12–16 September 2016; pp. 2300–2314. [Google Scholar]
  4. Li, X.; Ma, F.; Li, X.; Lv, H.; Bian, L.; Jiang, Z.; Zhang, X. LEO constellation-augmented multi-GNSS for rapid PPP convergence. J. Geod. 2018, 1–16. [Google Scholar] [CrossRef]
  5. Ge, H.; Li, B.; Ge, M.; Zang, N.; Nie, L.; Shen, Y.; Schuh, H. Initial assessment of Precise Point Positioning with LEO Enhanced Global Navigation Satellite Systems (LeGNSS). Remote Sens. 2018, 10, 984. [Google Scholar] [CrossRef]
  6. Ke, M.; Lv, J.; Chang, J.; Dai, W.; Tong, K.; Zhu, M. Integrating GPS and LEO to accelerate convergence time of precise point positioning. In Proceedings of the 2015 International Conference on Wireless Communications & Signal Processing (WCSP), Nanjing, China, 15–17 October 2015. [Google Scholar]
  7. Lawrence, D.; Cobb, H.S.; Gutt, G.; Connor, M.O.; Reid, T.G.; Walter, T.; Whelan, D. Ionnovation: Navigation from LEO; GPS World: Cleveland, OH, USA, 2017. [Google Scholar]
  8. Fossa, C.E.; Raines, R.A.; Gunsch, G.H.; Temple, M.A. An overview of the IRIDIUM (R) low Earth orbit (LEO) satellite system. In Proceedings of the IEEE 1998 National Aerospace and Electronics Conference (NAECON 1998), Dayton, OH, USA, 17 July 1998; pp. 152–159. [Google Scholar] [CrossRef]
  9. De Selding, P.B. LeoSat Corporate Broadband Constellation Sees GEO Satellite Operators as Partners. Space News. 14 June 2016. Available online: https://spacenews.com/leosat-corporate-broadband-constellation-sees-geo-satellite-operators-as-partners/ (accessed on 17 October 2018).
  10. De Selding, P.B. Telesat: LEO Gives More User Bandwidth than GEO HTS. Space Intel Report. 2017. Available online: https://www.spaceintelreport.com/telesat-leo-gives-more-user-bandwidth-than-geo-hts/ (accessed on 17 October 2018).
  11. De Selding, P.B. Virgin, Qualcomm Invest in OneWeb Satellite Internet Venture. Space News. 15 January 2015. Available online: http://spacenews.com/virgin-qualcomm-invest-in-global-satellite-internet-plan/sthash.XtMF5rpE.dpuf (accessed on 12 February 2017).
  12. De Selding, P.B. SpaceX To Build 4000 Broadband Satellites in Seattle. Space News. 19 January 2015. Available online: http://spacenews.com/spacexopening-seattle-plant-to-build-4000-broadband-satellites/ (accessed on 12 February 2017).
  13. Via Satellite. SpaceX Receives FCC Approval to Launch 7518 Starlink Satellites. 16 November 2018. Available online: https://www.satellitetoday.com/broadband/2018/11/16/spacex-receives-fcc-approval-to-deploy-7518-satellites/ (accessed on 14 December 2018).
  14. Via Satellite. Samsung Exec Envisions LEO Constellation for Satellite Internet Connectivity. 18 August 2015. Available online: https://www.satellitetoday.com/telecom/2015/08/18/samsung-exec-envisions-leo-constellation-for-satellite-internet-connectivity/ (accessed on 10 October 2018).
  15. De Selding, P.B. Boeing Proposes Big Satellite Constellations in V- and Cbands. Space News. 23 June 2016. Available online: http://spacenews.com/boeingproposes-big-satellite-constellations-in-v-and-c-bands/ (accessed on 17 October 2018).
  16. CNAGA (China National Administration of GNSS and Applications). CASIC Plans to Launch 156 Small Satellites for the Hongyun Program. 2017. Available online: http://en.chinabeidou.gov.cn/c/393.html (accessed on 3 March 2018).
  17. Jones, A. China to Build 300-Satellite Hongyan Communications Constellation in Low-Earth Orbit. GBtimes. 2018. Available online: https://gbtimes.com/china-to-build-300-satellite-hongyan-communications-constellation-in-low-earth-orbit (accessed on 17 October 2018).
  18. Lawrence, D.; Cobb, H.S.; Gutt, G.; Tremblay, F.; Laplante, P.; O’Connor, M. Test Results from a LEO-Satellite-Based Assured Time and Location Solution. In Proceedings of the 2016 International Technical Meeting of the Institute of Navigation, Monterey, CA, USA, 25–28 January 2016. [Google Scholar]
  19. Cobb, S.; Lawrence, D.; Gutt, G.; O’Conner, M. Differential and Rubidium Disciplined Test Results from an Iridium-Based Secure Timing Solution. In Proceedings of the 2017 International Technical Meeting of the Institute of Navigation, Monterey, CA, USA, 30 January–2 February 2017. [Google Scholar]
  20. Misra, P.; Enge, P. Global Positioning System: Signals, Measurements, and Performance, 2nd ed.; Ganga-Jamuna Press: Lincoln, MA, USA, 2011. [Google Scholar]
  21. Montenbruck, O.; Steigenber, P.; Hauschild, A. Broadcast versus precise ephemerides: A multi-GNSS perspective. GPS Solut. 2015, 12, 321–333. [Google Scholar] [CrossRef]
  22. Reid, T.G. Orbital Diversity for Global Navigation Satellite Systems. Ph.D. Thesis, Stanford University, Stanford, CA, USA, 2017. [Google Scholar]
  23. Su, H. Precise orbit determination of Global Navigation Satellite System of Second Generation (GNSS-2). Ph.D. Thesis, University FAF Munich, Neubiberg, Germany, 2000. [Google Scholar]
  24. China Satellite Navigation Office (CSNO). BeiDou Navigation Satellite System Signal in Space Interface Control Document Open Service Signal (Version 2.1); China Satellite Navigation Office: Beijing, China, 2016. [Google Scholar]
  25. Petit, G.; Luzum, B. IERS Conventions (2010) (IERS Technical Note; 36); Verlag des Bundesamts für Kartographie und Geodäsie: Frankfurt am Main, Germany, 2010; p. 179. ISBN 3-89888-989-6. [Google Scholar]
  26. Desai, S.D. Observing the pole tide with satellite altimetry. J. Geophys. Res. 2002, 107, 3186. [Google Scholar] [CrossRef]
  27. Saastamoinen, J. Atmospheric Correction for Troposphere and Stratosphere in Radio Ranging of Satellites. Geophysical Monograph; American Geophysical Union: Washington, DC, USA, 1972; pp. 247–252. [Google Scholar]
  28. Boehm, J.; Niell, A.; Tregoning, P.; Schuh, H. Global Mapping Function (GMF): A new empirical mapping function based on numerical weather model data. Geophys. Res. Lett. 2006. [Google Scholar] [CrossRef]
  29. Li, B.; Ge, H.; Ge, M.; Nie, L.; Shen, Y.; Schuh, H. LEO enhanced Global Navigation Satellite System (LeGNSS) for real-time precise positioning services. Life Sci. Space Res. 2019, 63, 73–93. [Google Scholar] [CrossRef]
Figure 1. The satellite distribution of the BeiDou-3 and LEO constellation.
Figure 1. The satellite distribution of the BeiDou-3 and LEO constellation.
Sensors 19 00198 g001
Figure 2. The distribution of selected ground stations (a) and ground track of LEO constellation (b).
Figure 2. The distribution of selected ground stations (a) and ground track of LEO constellation (b).
Sensors 19 00198 g002
Figure 3. Global distribution for satellite visibility of BeiDou-3-only (a) and BeiDou-3/LEO combined constellation (b).
Figure 3. Global distribution for satellite visibility of BeiDou-3-only (a) and BeiDou-3/LEO combined constellation (b).
Sensors 19 00198 g003aSensors 19 00198 g003b
Figure 4. Global distribution for PDOP of BeiDou-3-only (a) and BeiDou-3/LEO combined constellation (b).
Figure 4. Global distribution for PDOP of BeiDou-3-only (a) and BeiDou-3/LEO combined constellation (b).
Sensors 19 00198 g004
Figure 5. Time series of kinematic PPP positioning errors for BeiDou-3 and BeiDou-3 as well as LEO combined solutions. (a) CENT BeiDou-3 only; (b) CENT BeiDou-3 + LEO; (c) POTS BeiDou-3 only; (d) POTS BeiDou-3 + LEO; (e) NTUS BeiDou-3 only; (f) NTUS BeiDou-3 + LEO.
Figure 5. Time series of kinematic PPP positioning errors for BeiDou-3 and BeiDou-3 as well as LEO combined solutions. (a) CENT BeiDou-3 only; (b) CENT BeiDou-3 + LEO; (c) POTS BeiDou-3 only; (d) POTS BeiDou-3 + LEO; (e) NTUS BeiDou-3 only; (f) NTUS BeiDou-3 + LEO.
Sensors 19 00198 g005aSensors 19 00198 g005bSensors 19 00198 g005c
Figure 6. Time series of kinematic PPP positioning errors of CENT station in 1 h after the convergence of PPP for BeiDou-3 (a) and BeiDou-3 as well as LEO combined solutions (b), respectively.
Figure 6. Time series of kinematic PPP positioning errors of CENT station in 1 h after the convergence of PPP for BeiDou-3 (a) and BeiDou-3 as well as LEO combined solutions (b), respectively.
Sensors 19 00198 g006
Table 1. The summary of some deployed or proposed commercial communication LEO constellations.
Table 1. The summary of some deployed or proposed commercial communication LEO constellations.
ConstellationNo. SatsAltitude [km]Inclination [°]YearCountry
Iridium6678186.41998USA
Globalstar481400522000USA
LeoSat1081400Not defined, yet2019/2020USA
Telesat1171000~124599.52021Canada
Hongyun1561000Not defined, yet2022China
Hongyan3241100Not defined, yet2023China
OneWeb6481200882019USA/UK
Boeing2956120045, 55, 88Not defined, yetUSA
SpaceX Starlink75181110~132553, 53.8, 74, 81, 702020USA
Astrome Technology1501400Not defined, yet2020India
Samsung4600<1500Not defined, yetNot definedKorea
Table 2. Ambiguity Solution Related to the Altitude of Satellite.
Table 2. Ambiguity Solution Related to the Altitude of Satellite.
AltitudeConvergence Time for Ambiguity of Float SolutionTime for Integer Ambiguity Resolution
1000 km1 min10 min
10,000 km7 min1 h
20,000 km20 min4 h
IGSO 36,000 km2 h25 h
GEO 36,000 km+∞+∞
Table 3. Dynamic models used for generating the 24 h orbits for BeiDou-3 and LEO constellation.
Table 3. Dynamic models used for generating the 24 h orbits for BeiDou-3 and LEO constellation.
Geophysical ModelsDescription
BeiDou-3LEO
StaticEGM2008 up to degree and order 12Static part of EIGEN-6C up to degree and order 150
TemporalNoneTemporal part of EIGEN-6C up to degree and order 50
Secular rates for low degree coefficientsIERS Conventions 2010 [25]IERS Conventions 2010 [25]
n-bodyMoon, Sun, Mercury, Venus, Mars, Jupiter, Saturn, Uranus, Neptune, Pluto
JPL DE405
Moon, Sun, Mercury, Venus, Mars, Jupiter, Saturn, Uranus, Neptune, Pluto
JPL DE405
Solid Earth TidesIERS Conventions 2010 [25]IERS Conventions 2010 [25]
Ocean TidesNoneFES2004
Ocean pole tidesNoneDesai [26]
Relativistic effectsIERS Conventions 2010 [25]IERS Conventions 2010 [25]
Satellite surface models and Attitude
Atmospheric dragNoneDTM94 with Box-wing model
Solar Radiation Pressure5-parameter ECOM modelBox-wing
AttitudeNominal yaw-steering modelconsistent with local orbital reference frame
Reference frame
Inertial frameJ2000.0J2000.0
Earth tide and Ocean loadingIERS Conventions 2010 [25]IERS Conventions 2010 [25]
Precession/NutationIAU 2000AIAU 2000A
EOPIERS EOP 08 C04 (IAU2000A)IERS EOP 08 C04 (IAU2000A)
Table 4. The configures for observation simulation.
Table 4. The configures for observation simulation.
ConstellationBeiDou-3LEO
Cone angleGEO/IGSO: 1065
MEO: 15
PCO(600, 0, 11,000) mm(0, 0, 0)
PCV00
Satellite clockClock fileClock file
Mask elevation55
PCO00
PCV00
Receiver clock00
Solid/Ocean/Pole tideIERS Conventions 2010IERS Conventions 2010
Troposphere delaySaastamonion for dray and wet delaySaastamonion for dray and wet delay
GMFGMF
Ionosphere delayNoNo
Phase wide-upYesYes
RelativityYesYes
Code noise1 m1 m
Phase noise2 mm2 mm
Ambiguity00
Table 5. Statistical results for satellite visibility.
Table 5. Statistical results for satellite visibility.
ConstellationMax. # of Tracked SatellitesMin. # of Tracked SatellitesAvg. # of Tracked Satellites
3GEO + 3IGSO + 24MEO15.57.610.7
3GEO + 3IGSO + 24MEO + 120LEO21.610.116.3
Table 6. Statistical results for PDOP.
Table 6. Statistical results for PDOP.
ConstellationMax. PDOPMin. PDOPAvg. of PDOP
3GEO + 3IGSO + 24MEO2.241.191.63
3GEO + 3IGSO + 24MEO + 120LEO1.931.011.22
Table 7. SPP precision of both BDS and BDS + LEO (cm).
Table 7. SPP precision of both BDS and BDS + LEO (cm).
StationBDSBDS + LEO
NorthEastUpNorthEastUp
CENT99.1105.3293.793.2100.7261.3
POTS152.7131.8290.7121.387.4229.9
NTUS91.6114.9217.284.5101.5199.2
Table 8. PPP precision of both BDS and BDS + LEO after the ambiguity convergence (cm).
Table 8. PPP precision of both BDS and BDS + LEO after the ambiguity convergence (cm).
StationBDSBDS + LEO
NorthEastUpNorthEastUp
CENT0.741.122.050.590.651.23
POTS0.930.991.960.500.721.44
NTUS0.540.681.330.630.491.23

Share and Cite

MDPI and ACS Style

Su, M.; Su, X.; Zhao, Q.; Liu, J. BeiDou Augmented Navigation from Low Earth Orbit Satellites. Sensors 2019, 19, 198. https://doi.org/10.3390/s19010198

AMA Style

Su M, Su X, Zhao Q, Liu J. BeiDou Augmented Navigation from Low Earth Orbit Satellites. Sensors. 2019; 19(1):198. https://doi.org/10.3390/s19010198

Chicago/Turabian Style

Su, Mudan, Xing Su, Qile Zhao, and Jingnan Liu. 2019. "BeiDou Augmented Navigation from Low Earth Orbit Satellites" Sensors 19, no. 1: 198. https://doi.org/10.3390/s19010198

APA Style

Su, M., Su, X., Zhao, Q., & Liu, J. (2019). BeiDou Augmented Navigation from Low Earth Orbit Satellites. Sensors, 19(1), 198. https://doi.org/10.3390/s19010198

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop