Next Article in Journal
Dynamics and Embedded Internet of Things Input Shaping Control for Overhead Cranes Transporting Multibody Payloads
Next Article in Special Issue
GNSS Code Multipath Mitigation by Cascading Measurement Monitoring Techniques
Previous Article in Journal
Sparse Method for Direction of Arrival Estimation Using Denoised Fourth-Order Cumulants Vector
Previous Article in Special Issue
Energy Efficient GNSS Signal Acquisition Using Singular Value Decomposition (SVD)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Real-Time Monitoring for BDS Signal-In-Space Anomalies Using Ground Observation Data

1
Chinese Antarctic Center of Surveying and Mapping, Wuhan University, Wuhan 430079, China
2
State Key Laboratory of Geodesy and Earth’s Dynamics, Institute of Geodesy and Geophysics, Chinese Academy of Sciences, Wuhan 430077, China
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(6), 1816; https://doi.org/10.3390/s18061816
Submission received: 9 May 2018 / Revised: 1 June 2018 / Accepted: 1 June 2018 / Published: 4 June 2018
(This article belongs to the Special Issue GNSS and Fusion with Other Sensors)

Abstract

:
Signal-in-space (SIS) User Range Error (URE) is one of the major error sources for BeiDou Navigation Satellite System (BDS) applications and can reach tens of meters or even more. Therefore, real-time monitoring of SIS anomalies has a great realistic significance to guarantee the safety of users. According to an analysis of the BDS navigation messages, it showed that the User Range Accuracy (URA) index could not reflect the change of URE when it was abnormal. The conventional models using the relationship between URA and URE to monitor SIS anomalies are not suitable to the present BDS. Therefore, we use a prior information of SIS URE derived from ground observational data instead of URA to monitor BDS SIS anomalies. In order to realize the corresponding functions, we analysed the distribution of SIS UREs and obtained their prior models. Then, the monitoring threshold is determined using the prior models and a confidence interval instead of URA. The scheme was tested by applying to BDS SIS anomalies monitoring based on 13 ground tracking stations. The performance of this method was assessed by comparison with the satellite-health indicators from broadcast ephemeris. The results confirm that the method developed in this paper can rightly and timely detect abnormal SIS.

1. Introduction

The BeiDou Navigation Satellite System (BDS) was officially put into operation on 27 December 2012 and can overlap the Earth surface in 2020 s [1]. The service modes of BDS include the open service (OS) and the authorized service (AS) [2]. The Positioning, Navigation and Timing (PNT) performance of OS has basically reached or exceeded the design-performance requirements [3,4].
For most BDS OS users, real-time satellite orbits and clocks are derived from predicted ephemeris and clock parameters in navigation messages broadcast by BDS satellites. Signal-in-Space (SIS) errors arise primarily from imperfect estimation of a satellite orbit and clock error and are usually undetectable and uncorrectable for stand-alone OS users [5]. It is mainly described by the parameter of User Range Error (URE) that is the pseudorange inaccuracy attributable to the ground control and space vehicles and the one of major error sources affecting Global Navigation Satellite System (GNSS) applications [6]. The URE are dominated by the BDS Space and Ground Control Segment and do not include the error budget components assigned to the BDS User Segment such as the tropospheric or ionospheric modelling errors, multipath effects, and receiver noise. The constellation of BDS is different from other satellite navigation systems, including three types of orbits: medium earth orbits (MEO), inclined geostationary (IGSO), and geostationary orbits (GEO). SIS UREs of different type of BDS satellites are different [7,8]. According to the statistical analysis and evaluation on the broadcast ephemeris data, the SIS URE accuracy of the BDS is better than 2.5 m compared against the precise ephemeris data [9].
Nominally, OS users can assume that broadcast navigation message is reliable and the URE derived from a healthy SIS is at the meter level. In practice, however, SIS anomalies occurred occasionally and UREs of tens of meters or even more were observed, which could lead to hazardous misleading position solutions for unaugmented receivers [6]. Therefore, real-time monitoring of SIS anomalies has great importance for ensuring the safety of BDS users.
The study of monitoring SIS anomalies has long been conducted for Global Positioning System (GPS) [10,11,12,13,14], however, few studies have focused on BDS. Present detecting methods mainly depend on the relationship between User Range Accuracy (URA) in the broadcast navigation messages and URE to monitor SIS condition [6,13,14]. The URA is a conservative representation of the standard deviation of URE at the worst cast location on the earth [15]. In the case of GPS, the SPS SIS standard assures that the URE should be less than 4.42 times the upper bound on the URA value ( U R A U B ) for any healthy SIS during normal operation [16]. However, according to statistical analysis of BDS navigation messages for a long time, it showed that the BDS URA index values were constant, even in the case of satellite failure [17]. In the case of BDS, the specific definitions of the signal in SIS accuracy index parameters of broadcast navigation messages are not published in the current ICDs [18]. Therefore, the URA of BDS is unsuitable to construct a threshold to monitor the BDS SIS anomaly. Besides, the empirical values are also used as thresholds to detect BDS anomalies [19]. However, some error sources, such as multipath effect and receiver noise, are also closely related to monitoring stations, so the random model based on general hypothesis has defects that may result in false detection and false alarm in the practice use. In those cases, knowledge of characteristics of the BDS SIS URE has a great importance for monitoring the SIS anomalies.
The objective of the work presented here is to propose a scheme to monitor the BDS SIS anomalies using ground observation data. This program constructed a random model that was consistent with the actual URE based on the features of the BDS constellation and the observation conditions of the regional monitoring stations. According to the random models and a given confidence interval to obtain the threshold of URE instead of URA, the method is called BDS Ground SIS Monitoring (BDSGSISM). In this contribution, the BDS SIS UREs are derived from BDS satellites tracked by 13 ground stations of Multi-GNSS Experiment campaign (MGEX) of International GNSS Service (IGS) in the Asia-pacific region via deducting non-SIS errors from the total pseudorange errors in Day of Year (DOY) 36–56, 2017. In order to ensure the originality of broadcast navigation messages, the broadcast navigation messages are cleaned based on majority voting [13]. Then, we analyse SIS URE time series and calculate their distribution parameters. Based on the prior information of SIS URE, BDSGSISM is used to monitor the BDS healthy status for 27 days from 1 November 2015 to 27 November 2015 and presents the potential SIS anomalies during that time.
The remainder of this contribution is organized as follows: Section 2 introduces the methodology of BDSGSISM, the data source, the real-time broadcast ephemeris and the computation of SIS UREs. Section 3 analyses the characteristics of BDS SIS UREs. Section 4 presents two cases of SIS anomaly monitoring, and Section 5 provides concluding remarks.

2. Materials and Methods

2.1. BDS Ground SIS Monitoring Method

Satellite-health status has specific markers in the broadcast ephemeris, however, if there is any sudden and unexpected failure, the Operational Control Segment (OCS) will guarantee a failure warning in the next broadcast ephemeris. This is insufficient for the users of BDS. Traditional SIS anomalies are monitored on the basis of the relationship between SIS UREs and URA as follows: if the SIS UREs are greater than 4.42 × U R A U B , it can be determined that the satellite service is not guaranteed [12]. However, the URA index in BDS broadcast ephemeris is not equivalent to that of GPS [17]. Figure 1 shows the time series of U R A U B and satellite healthy flag (0: health, 1: anomaly) from BDS broadcast ephemeris for the C04, C09 and C14 during the year 2017. The U R A U B remains constant even if the satellites fail. In this case, a prior model of the SIS UREs of each satellite can be obtained through analysing the SIS UREs for a period of time from the selected tracking stations. Then, the thresholds of SIS UREs can be determined using the prior model and confidence level, instead of the URA. The specific processes are as follows:
  • Combining cleanly BDS ephemeris based on real-time raw data from tracking network.
  • The SIS UREs of all the satellites are calculated by the selected tracking stations’ observations.
  • The mean value and the standard deviation of the SIS UREs of each satellite observed by each station are calculated, as the prior information.
  • The threshold of the SIS UREs is calculated using Equation (1) based on a given confidence level and the prior information obtained in step 3.
    P r ( | U R E U R E ¯ | κ ( P r ) σ U R E ) = α
    where U R E ¯ is the mean value of the SIS UREs, σ U R E is the standard deviation of the SIS UREs, and κ ( P r ) is the quantile of the corresponding confidence α . In this step, we assume that the SIS UREs obey a Gaussian distribution.
In this paper, the monitored satellite-health status is represented by one of the following three numbers: −1 (abnormal), 0 (cannot be determined), and 1 (normal). Determining the health status of a satellite requires more than three tracking stations observing the satellites simultaneously; otherwise, the status of the satellite cannot be determined and is flagged 0. If more than three tracking stations detecting anomalies (URE exceeds a given confidence interval) at the same time then it is determined that the satellite cannot be used, and the status is flagged with −1. If a satellite is observed by more than three tracking stations at the same time and the SIS UREs for fewer than three tracking stations exceed the not-to-exceed (NTE) limit, then it is determined that the satellite is working normally, and the status is flagged 1.
The method may experience a missed alarm in the following situations: (1) when the satellite orbit is abnormal, but the projection of the orbital deviation is small in the monitoring service area, therefore, the monitoring stations cannot detect the abnormality. The effect of the anomaly is also not serious enough to cause dangerous misleading information for the users in the service area; (2) when a navigation satellite fails and misleads to users in the service area, but less than three monitoring stations simultaneously detect this anomaly in the service area. The probability of such event can be expressed as:
P = P 0 + P 1 + P 2
where P i ( i = 0 , 1 , 2 ) is the probability of i stations that monitored the abnormal behavior of satellite when it suffered a failure. Then the probability of effectively monitoring the satellite failure is 1 - P . The probability of a missed alarm is related to the tracking station density and the receiver stability. The denser the network of monitoring stations, the lower the probability of a missing alarm will be.

2.2. Data Source

In order to verify the reliability of the BDSGSISM, the BDS observation data from MGEX of International GNSS Service (IGS) for the years 2015 and 2017 have been used. As BDS mainly serves the Asia-pacific region, we selected 13 stations over the region and the nearby area as the experimental data sources. The specific information of these tracking stations including the tracking-station name, the receiver type, the antenna type and the BDS satellite is summarized in Table 1 and the distribution is shown in Figure 2.
The MGEX tracking station has the following advantages: (1) the selected tracking station can provide dual-frequency observation on B1 and B2 bands; (2) the MGEX tracking station can provide 1 Hz observation data, which is conducive to monitoring short-time satellite anomalies; and (3) we can obtain relatively long-time BDS observation data. To verify the reliability of the algorithm, we re-encoded broadcast ephemeris and observation data that have been broadcasted and transmitted in the form of real-time streams.

2.3. Combined Real-Time Broadcast Ephemeris

Due to accidental bad receiver data and various hardware/software bugs, a small proportion of the navigation data from the tracking stations have defects, such as losses, duplications, inconsistencies, discrepancies, and errors. Therefore, to accurately and effectively monitor BDS SIS anomalies, the navigation messages must be guaranteed authenticity. In this paper, we use a systematic methodology to combine BDS ephemeris using real-time data of simulation from multiple tracking stations.
Currently, the combined broadcast ephemeris provided by IGS and other institutions mostly take pseudorandom noise (PRN) and Time of Clock (TOC) (the meaning of these broadcast ephemeris parameters is same as in the reference [1], similarly hereinafter) as keywords to search and integrate ephemeris, but existing research [13] have shown that these two parameters have errors. In this paper, the broadcast ephemeris parameters are divided into robust parameters and fragile parameters [13]. The robust parameters are utilized to identify the equivalence of two navigation messages as follows: two navigation messages are deemed identical if and only if they agree on all the robust parameters, although their fragile parameters could be different. Most orbital and clock parameters in navigation messages are usually reported correctly, and even when errors occur, only a few stations agree on the same incorrect value [13]. These parameters are referred to as robust parameters, including: C r s , Δ n , M 0 , C u c , e , C u s , A , C i c , Ω , C i s , i 0 , C r c , ω , and Ω ˙ . On the contrary, other parameters, such as P R N , T O C , I D O E , and i ˙ are more likely to be erroneous, and when errors occur, several stations may make the same mistakes. These parameters are referred to as fragile parameters.
The broadcast ephemeris can be combined using the algorithm described below.
  • Broadcast ephemeris from a number of tracking stations are simultaneously received and added into the set O. Ensure that each satellite can be observed by a number of tracking stations (generally more than three tracking stations) at the same time.
  • For each navigation message e in O, if there is already a navigation message f in the set P (its database of different robust parameters) having the same robust parameters as e, then add the fragile parameters of e into f ’s database; otherwise, add e into P.
  • For each navigation message f in P, apply majority vote to each fragile parameter (except the Transmission Time of Message (TTOM)) according to f’s database, and record the number of stations that report f.

2.4. Computation of SIS UREs

The observation equation for the pseudorange observables can be modelled as:
P r , i s = ρ r s + c ( δ t r δ t s ) + T r s + μ i I r , 1 s + ε
where superscript s, subscripts r and i denote a specific satellite, receiver, and frequency band, respectively; P r , i s is the raw code observation; ρ r s is the geometric distances; δ t r and δ t s are the clock offsets for receiver and satellite, respectively; T r s is the tropospheric delay; I r , 1 s is the slant ionospheric delay on B1 signal and μ i = ( f 1 / f i ) 2 is a coefficient related to the frequency, in which f i denotes the frequency of Bi signal; ε are the observation noise and unmodeled effects.
For the real-time users computing satellite orbits and clock offsets based on broadcast navigation messages, the geometric distance and satellite clock offsets can be written as:
{ ρ r s = ρ r s ^ + φ s δ t s = δ t s ^ + ϕ s
where ρ r s ^ is the geometric distance derived from broadcast ephemeris, φ s is the broadcast ephemeris error projected onto the line-of-sight (LOS) from a BDS satellite to a receiver, δ t s ^ is the broadcast satellite clock offset and ϕ s is the broadcast clock error.
Based on Equations (3) and (4), we can obtain the SIS URE, that is [6]
S I S   U R E r , i s = c ϕ s φ s = ρ r s ^ + c ( δ t r δ t s ^ ) + T r s + μ i I r , 1 s + ε P r , i s .
In order to obtain the SIS URE, T r s is derived from Saastamoinen model [20], ρ r s ^ is computed from the broadcast ephemeris and the receiver position obtained from the “igs.snx” file (The file contains a summary of the station logs and is maintained at the IGS central bureau.) and δ t s ^ is computed based on broadcast navigation messages. However, the slant ionospheric delay and receiver clock offset are estimated using least squares estimation. The detailed derivation is as follows:
[ c δ t r I r , 1 1 I r , 1 2 I r , 1 n ] T = ( A T W A ) 1 A T W L
where A = [ 1 μ 1 0 0 1 μ 2 0 0 1 0 μ 1 0 1 0 μ 2 0 1 0 0 μ 1 1 0 0 μ 2 ] , L = [ P r , 1 1 ρ r 1 ^ + c δ t 1 ^ T r 1 P r , 2 1 ρ r 1 ^ + c δ t 1 ^ T r 1 P r , 1 2 ρ r 2 ^ + c δ t 2 ^ T r 2 P r , 1 2 ρ r 2 ^ + c δ t 2 ^ T r 2 P r , 1 n ρ r n ^ + c δ t n ^ T r n P r , 1 n ρ r n ^ + c δ t n ^ T r n ] , W = d i a g ( [ sin 2 ( e l r 1 ) sin 2 ( e l r 1 ) sin 2 ( e l r 2 ) sin 2 ( e l r 2 ) sin 2 ( e l r n ) sin 2 ( e l r n ) ] ) , d i a g ( ) is a function of diagonal matrices, and e l is the satellite elevation angle.

3. BDS SIS URE Analysis and Discussion

For the BDS, and SIS URE the pseudorange inaccuracy, attributable to the ground control and the space vehicles, is one of the major error sources affecting BDS application. SIS URE refers to errors caused by satellite segment, including satellite ephemeris and clock errors, satellite antenna variations, and signal imperfections. However, SIS anomalies are mainly caused by satellite ephemeris and clock errors because antenna variations and signal imperfections are at a level of millimeter or centimeter [21]. As SIS anomalies occur occasionally and UREs of tens of metres or even more have been observed, knowledge of the SIS UREs is of great importance for developing GNSS SIS-monitoring systems.
The SIS UREs can be utilized to identify the satellite fault, and their description and assessment is very significant for BDS-SIS monitoring [7]. At present, most of the existing research assumes that the SIS UREs follow a Gaussian distribution with a mean value of 0 [22]. However, some researches have revealed that the code-peseudorange measurements of the MEO and IGSO satellites of the BDS exist nonzero mean values and have a certain relationship with the satellite elevation angle [23,24]. We have calculated the BDS SIS UREs based on the elevation weighting, as shown in Figure 3, Figure 4 and Figure 5.
Observations from 21-day intervals of DOY 036–056, 2017 were selected as the core data sets for studying BDS SIS UREs. The Figure 3, Figure 4 and Figure 5 show three types of URE time series (C04, C09, and C14) observed by three stations (NNOR, XMIS, and GMSD) with different latitudes. According to the time series of SIS UREs, which are very similar across all stations, there are obvious characteristics: the time series of C04 and C09 are dominated by some sinusoidal waves, however, the time series of C14 are not continuous because the satellite cannot be continuously observed by tracking stations. According to a comparison between broadcast ephemeris data and precise ephemeris product, the broadcast orbit errors exit significant periodic fluctuation in radial, along, and cross direction [7]. Therefore, the fluctuations of SIS UREs are mainly caused by the broadcast orbit errors.
Although the SIS UREs are generally assumed to be zeroes-mean, the reality may be different [23,25]. To monitor the SIS anomalies, it is necessary to eliminate the periodic system error in the SIS UREs, so that they satisfy a Gaussian distribution. In the section, we utilize the Fourier series to model the SIS UREs for each satellite observed by selected tracking stations. As to SIS UREs after eliminating the trend term, the traditional method can be used to monitor the health status of the satellite. Based on the long-term observation data, the period and amplitude of the UREs are analyzed, and the finite-term Fourier trigonometric-series model is established for the non-random term of the UREs. In this paper, we use the URE data during the recent 14 days to obtain the Fourier-series model and forecast the next seven days of the UREs. Assuming that the URE time series ( y ) can be expressed in the form of Equation (7) in the time domain, it can be represented by a set of independent forms, such as the Fourier trigonometric-series function, given in Equation (8).
y T = [ y 1 y 2 y m ]
{ E { y ( t ) } = A 0 + k = 1 q A k x k A k = [ cos ω k t 1 sin ω k t 1 cos ω k t 2 sin ω k t 2 cos ω k t m sin ω k t m ] , x k = [ a k b k ] .
where A 0 is the mean value of the coefficient and is a constant term; q is the number of trigonometric functions; ω k is the frequency of trigonometric function k, corresponding to the period of time series y T ; and a k , b k are the model coefficients to be estimated. If a k , b k and ω k are known, then y can be calculated at any time.
The prior models of SIS UREs can be written as:
{ U R E ¯ = 1 n i = 1 n ( U R E i U R E i f ) σ U R E = 1 n i = 1 n ( U R E i U R E i f ) 2
where U R E i f is computed from Fourier trigonometric-series function.
Representative results of such Fourier-series models are presented in Figure 6. The figure shows the results for three different types of satellites (PRN C04, C09, and C14) on B1 band, observed by GMSD station in 2017 (DOY 036–056). The red thin lines in Figure 6 represent the Fourier-series models for the SIS UREs (blue dots). The model parameters are determined according to the UREs during the 14 days (DOY 036–049), and used to forecast the next seven days (DOY 050–056).It can be seen that the Fourier-series models can fit the SIS UREs variation trend very well.
By analyzing the SIS UREs (from which the periodic system error has been eliminated) from October 2016 to March 2017, there is a certain difference among the means and the standard deviations of the different satellites at different frequencies observed by different stations, as shown in Figure 7, Figure 8 and Figure 9. Besides, the statistical results of the SIS UREs (m) after eliminating the trend term are shown in Table 2. The overall level of the difference depends on certain characteristics of the signal itself and the tracking performance of the specific receiver and the antenna quality with respect to its multipath sensitivity [6]. Furthermore, the surroundings of the receiving antenna and its multipath contamination affect the SIS UREs values [26,27]. The computation method of SIS UREs also results in the differences at different frequencies. In addition, the observation accuracy of different BDS satellite is also different, varying from 0.5 m to 1.8 m [7]. Therefore, various models instead of the same model are preferable to monitor SIS UREs for different stations, satellites and frequencies.

4. Case Studies of BDS SIS Anomalies

SIS anomaly monitoring requires real-time data streams of the tracking stations, however, for various reasons, the current MGEX tracking stations have not provided the BDS real-time data streams. In order to verify the reliability of the proposed algorithm, we re-encoded the broadcast ephemeris and the observation data that have been broadcasted and transmitted in the form of real-time streams. Through processing the simulated real-time data streams, some potential SIS anomalies are found and listed in Table 3. The parameter κ ( P r ) in Equation (1) is set to 4.42 based on the relationship between URA and URE as well as integrity monitoring requirements. An SIS anomaly is claimed when both of the following conditions are fulfilled: (1) a satellite is observed by more than three tracking stations and its SIS UREs derived from more than three stations exceed the given confidence interval at the same time; (2) the unhealthy flag of the broadcast ephemeris is not set on.
For the monitoring results in Table 3, in-depth case studies of the C02 anomaly on 12 December 2015 are presented. Figure 10 shows the health statuses of C02, monitored by BDSGSISM (blue asterisks), compared with the health statuses given by the broadcast ephemeris (red dots). As illustrated in Figure 10, the C02 satellite occurred two periods of abnormal performances on that day based on the information provided by the BDSGSISM as well as the broadcast ephemeris. Although the two results showed that each anomaly lasted for one hour, starting times of each anomaly are not consistent for two methods. This example indicates that at least one monitoring result is wrong. Therefore, we analysed the SIS URE time series of C02 and found that C02 could have been simultaneously observed by six to seven tracking stations at the two periods of anomaly, and which all monitored these two anomalies. The monitoring results of BDSGSISM are better than those from the broadcast ephemeris since it timely and accurately detected the anomaly.
Figure 11 shows the SIS UREs time series of C02 observed by a MGEX station, NNOR on DOY 346, 2015. The left panel shows the SIS URE time series of C02 whole day tracked by NNOR, and the right panel shows the SIS URE of C02 by NNOR with the anomalous periods removed. The first anomaly experienced by NNOR started at 11:00 and ended at 12:00. The UREs suddenly increased and exceeded the confidence interval during this period. The anomaly may have resulted from the unplanned operations. Therefore, broadcast ephemeris cannot alarm on time. When anomalies were detected, the operation control centre immediately adjusted the health status of the satellite, but the TTA (Time to Alarm) of the broadcast ephemeris was too long, which may impact some applications requiring the safety. The BDSGSISM can make up for this shortcoming.
According to the analysis of the monitoring results, it can be found that in the BDS broadcast ephemeris there exist missed detection as well as false alarms. Figure 12 shows the health statuses of C10 monitored by BDSGSISM (blue asterisks) compared with the health statuses provided by the broadcast ephemeris (red dots) in 25 December 2015 (DOY 359, 2015). We can find that the health statuses monitored by BDSGSISM do not match well with those provided by the broadcast ephemeris. According to the broadcast ephemeris, the anomaly started at 08:00 and ended at 19:00, indicating that the anomaly lasted for 11 h. However, the monitoring result of BDSGSISM shows that the anomaly started at 08:46 and ended at 17:00, indicating that the anomaly lasted for only 8.24 h. Therefore, we analysed the SIS UREs observed by the tracking-station network. Between 08:00 and 08:46, there were 9 to 11 tracking stations to observe C10 and all show it to be operating normally. At about 08:46, the UREs began to increase with time, and exceeded confidence intervals, as three stations simultaneously monitored the anomaly, and this lasted for only 8.24 h, four to eight stations simultaneously monitored the anomaly during the period. At about 17:00, the satellite restored healthy. As illustrated in Figure 13, we can find that the SIS UREs exceeded the confidence interval at 08:46 and gradually increased until returning to normal at 17:00.

5. Conclusions

BDS SIS anomaly monitoring is essential to guarantee the safety of the user’s life and property. SIS anomaly is one of the major risk sources for user’s PNT services. However, the traditional monitoring method of SIS anomaly cannot be applied to BDS due to invalid URA parameter in BDS broadcast ephemeris. Therefore, it is essential to take an effective approach to monitor the SIS anomaly for the wide application of BDS. In this paper, we propose BDSGSISM method based on the prior information (mean and variance) of SIS UREs to monitor BDS SIS anomalies.
In order to guarantee authenticity of broadcast ephemeris, the cleaned navigation messages are combined based on broadcast ephemeris received from a number of tracking stations. According to analysing BDS SIS UREs, it shows that the SIS UREs do not obey a Gaussian distribution due to systemic errors. Therefore, we utilize Fourier series to model the SIS UREs to eliminate the periodic system error, making them obey a Gaussian distribution and satisfy the precondition of SIS-anomaly monitoring. Besides, we calculate the SIS URE’s means and variances for all satellites observed by selected tracking stations and find that there are differences in these values for different satellites. Finally, the NTE threshold is determined by the prior information of SIS UREs and a given confidence interval, instead of URA.
To verify the reliability of BDSGSISM, the collected BDS data of 57 days (DOY 305–361, 2015) were analysed using BDSGSISM and presented the potential signal-in-space anomalies. According to a comparison between the satellite health states obtained by BDSGSISM monitoring and those given by broadcast ephemeris, it can be found that there exist false alarms and missed detections in the results obtained by broadcast ephemeris, however, BDSGSISM can quickly and accurately find the abnormal satellite and overcome the shortcoming.

Author Contributions

H.J. and H.W. provided the initial idea for this study; H.J., H.W. conceived and designed the experiments; H.J. analyzed the data; H.J. and H.W. wrote the main manuscript; Z.W. and Y.Y. reviewed the paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (41474029, 41776195, 41574015), the Science and Technology Innovation Fund of Chinese Academy of Sciences (CXJJ-14-M18), the National Basic Research Program of China (2012CB825604), and the CAS/SAFEA International Partnership Program for Creative Research Teams (KZZD-EW-TZ-05). We also would like to acknowledge the International GNSS Service (IGS) for providing the GPS data and the products used in our experiment.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. CSNO. Beidou Navigation Satellite System Signal in Space Interface Control Document, Open Service Signal (Version 2.0); China Satellite Navigation Office (CSNO): Beijing, China, 2013. [Google Scholar]
  2. Han, C.; Yang, Y.; Cai, Z. Beidou navigation satellite system and its time scales. Metrologia 2011, 48, S213. [Google Scholar] [CrossRef]
  3. Shi, C.; Zhao, Q.L.; Hu, Z.G.; Liu, J.N. Precise relative positioning using real tracking data from compass geo and igso satellites. GPS Solut. 2013, 17, 103–119. [Google Scholar] [CrossRef]
  4. Yang, Y.; Li, J.; Wang, A.; Xu, J.; He, H.; Guo, H.; Shen, J.; Dai, X. Preliminary assessment of the navigation and positioning performance of beidou regional navigation satellite system. Sci. China-Earth Sci. 2014, 57, 144–152. [Google Scholar] [CrossRef]
  5. Cohenour, C.; Graas, F.V. GPS orbit and clock error distributions. Navigation 2011, 58, 17–28. [Google Scholar] [CrossRef]
  6. Heng, L.; Gao, G.X.; Walter, T.; Enge, P. Automated verification of potential GPS signal-in-space anomalies using ground observation data. In Proceedings of the IEEE/ION PLANS 2012, Myrtle Beach, SC, USA, 24–26 April 2012; pp. 1111–1118. [Google Scholar]
  7. Chen, L.; Jiao, W.; Huang, X.; Geng, C.; Ai, L.; Lu, L.; Hu, Z. Study on signal-in-space errors calculation method and statistical characterization of beidou navigation satellite system. In Lecture Notes in Electrical Engineering, Proceedings of the China Satellite Navigation Conference (CSNC), Wuhan, China, 15–17 May 2013; Springer: Berlin, Germany, 2013; Volume 243, pp. 423–434. [Google Scholar]
  8. Shi, C.; Zheng, F.; Lou, Y.; Gu, S.; Zhang, W.; Dai, X.; Li, X.; Guo, H.; Gong, X.; Shi, C. National bds augmentation service system (nbass) of China: Progress and assessment. Remote Sens. 2017, 9, 837. [Google Scholar]
  9. Liu, S.A.; Jia, X.L.; Ma, Z.X. The analysis of GNSS SIS accuracy. In Proceedings of the China Satellite Navigation Conference, Xi’an, China, 13–15 May 2015; Sun, J., Liu, J., Fan, S., Lu, X., Eds.; Springer: New York, NY, USA, 2015; Volume 340, pp. 375–386. [Google Scholar]
  10. Wu, T.; Peck, S. An analysis of satellite integrity monitoring improvement for was. In Proceedings of the 15th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GPS 2002), Portland, OR, USA, 24–27 September 2002; pp. 756–765. [Google Scholar]
  11. DiLellio, J.A. GPS Navigation System with Integrity and Reliability Monitoring Channels. US Patent 7400292, 15 July 2008. [Google Scholar]
  12. Gao, G.X.; Tang, H.C.; Blanch, J.; Lee, J.; Walter, T.; Enge, P. Methodology and case studies of signal-in-space error calculation top-down meets bottom-up. In Proceedings of the 22nd International Technical Meeting of the Satellite Division of the Institute of Navigation (Ion Gnss 2009), Savannah, GA, USA, 22–25 September 2009; pp. 2824–2831. [Google Scholar]
  13. Heng, L.; Gao, G.X.; Walter, T.; Enge, P.; Inst, N. GPS ephemeris error screening and results for 2006–2009. In Proceedings of the 2010 International Technical Meeting of the Institute of Navigation—ITM 2010, San Diego, CA, USA, 25–27 January 2010; pp. 1014–1022. [Google Scholar]
  14. Lee, Y.C. Ground monitoring schemes for GPS integrity channel. IEEE Trans. Aerosp. Electron. Syst. 1990, 26, 481–489. [Google Scholar] [CrossRef]
  15. Jan, S.S.; Tao, A.L. The open service signal in space navigation data comparison of the global positioning system and the beidou navigation satellite system. Sensors 2014, 14, 15182–15202. [Google Scholar] [CrossRef] [PubMed]
  16. Heng, L.; Gao, G.X.; Walter, T.; Enge, P. GPS signal-in-space integrity performance evolution in the last decade. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 2932–2946. [Google Scholar] [CrossRef]
  17. Dins, A.; Ping, Y.; Schipper, B. Statistical characterization of beidou and GPS sis errors in the Asian region. In Proceedings of the 2015 IEEE/AIAA 34th, Digital Avionics Systems Conference, New York, NY, USA, 13–18 September 2015. [Google Scholar]
  18. CNSO. Beidou Navigation Satellite System Signal in Space Interface Control Document Open Service Signal B2A (Version 1.0); China Satellite Navigation Office: Beijing, China, 2017. [Google Scholar]
  19. Ye, F.; Yuan, Y.; Tan, B.; Ou, J. A robust method to detect beidou navigation satellite system orbit maneuvering/anomalies and its applications to precise orbit determination. Sensors 2017, 17, 1129. [Google Scholar] [CrossRef] [PubMed]
  20. Saastamoinen, J.; Richard, J.R. Simultaneous drying and pyrolysis of solid fuel particles. Combust. Flame 1996, 106, 288–300. [Google Scholar] [CrossRef]
  21. Schmid, R.; Rothacher, M.; Thaller, D.; Steigenberger, P. Absolute phase center corrections of satellite and receiver antennas. GPS Solut. 2005, 9, 283–293. [Google Scholar] [CrossRef]
  22. Shively, C.A. A gaussian mixture model for error distributions used in assessing raim performance. In Proceedings of the 2012 International Technical Meeting of the Institute of Navigation, Newport Beach, CA, USA, 30–31 January 2012; pp. 1590–1602. [Google Scholar]
  23. Montenbruck, O.; Hauschild, A.; Steigenberger, P.; Hugentobler, U.; Teunissen, P.; Nakamura, S. Initial assessment of the compass/beidou-2 regional navigation satellite system. GPS Solut. 2013, 17, 211–222. [Google Scholar] [CrossRef]
  24. Wanninger, L.; Beer, S. Beidou satellite-induced code pseudorange Variations: Diagnosis and therapy. GPS Solut. 2015, 19, 639–648. [Google Scholar] [CrossRef]
  25. Montenbruck, O.; Steigenberger, P.; Hauschild, A. Broadcast versus precise Ephemerides: A multi-gnss perspective. GPS Solut. 2015, 19, 321–333. [Google Scholar] [CrossRef]
  26. Dong, D.; Wang, M.; Chen, W.; Zeng, Z.; Song, L.; Zhang, Q.; Cai, M.; Cheng, Y.; Lv, J. Mitigation of multipath effect in gnss short baseline positioning by the multipath hemispherical map. J. Geod. 2016, 90, 255–262. [Google Scholar] [CrossRef]
  27. Panicciari, T.; Soliman, M.A.; Moura, G. Evaluation of a gnss receiver performance in different multipath environments with a novel real-time multipath simulation system. In Proceedings of the 30th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GNSS + 2017), Portland, OR, USA, 25–29 September 2017; pp. 427–435. [Google Scholar]
Figure 1. Time series of User Range Accuracy (URA) and satellite healthy flag (0: health, 1: anomaly) from BeiDou Navigation Satellite System (BDS) broadcast ephemeris for the C04 (top), C09 (middle), and C14 (bottom) during the year 2017.
Figure 1. Time series of User Range Accuracy (URA) and satellite healthy flag (0: health, 1: anomaly) from BeiDou Navigation Satellite System (BDS) broadcast ephemeris for the C04 (top), C09 (middle), and C14 (bottom) during the year 2017.
Sensors 18 01816 g001
Figure 2. Multi-GNSS Experiment campaign (MGEX) tracking network used in this study.
Figure 2. Multi-GNSS Experiment campaign (MGEX) tracking network used in this study.
Sensors 18 01816 g002
Figure 3. Signal-in-Space-User Range Error (SIS-URE) time series of Station NNOR, C04 (top), C09 (middle) and C14 (bottom), DOY 036–064, 2017.
Figure 3. Signal-in-Space-User Range Error (SIS-URE) time series of Station NNOR, C04 (top), C09 (middle) and C14 (bottom), DOY 036–064, 2017.
Sensors 18 01816 g003
Figure 4. SIS-URE time series of Station XMIS, C04 (top), C09 (middle) and C14 (bottom), DOY 036–056, 2017.
Figure 4. SIS-URE time series of Station XMIS, C04 (top), C09 (middle) and C14 (bottom), DOY 036–056, 2017.
Sensors 18 01816 g004
Figure 5. SIS-URE time series of Station GMSD, C04 (top), C09 (middle) and C14 (bottom), DOY 036–056, 2017.
Figure 5. SIS-URE time series of Station GMSD, C04 (top), C09 (middle) and C14 (bottom), DOY 036–056, 2017.
Sensors 18 01816 g005
Figure 6. Examples of SIS-URE time series of Station GMSD, C04 (top), C09 (middle) and C14 (bottom), DOY 036–056, 2017, SIS UREs (blue) and Fourier models (red) on B1 band.
Figure 6. Examples of SIS-URE time series of Station GMSD, C04 (top), C09 (middle) and C14 (bottom), DOY 036–056, 2017, SIS UREs (blue) and Fourier models (red) on B1 band.
Sensors 18 01816 g006
Figure 7. The histograms of C04 SIS UREs (m) after eliminating the trend term covering fourteen experimental days, overlaid with empirical normal probability-density function (PDF) curves (in red). Six subplots, arranged in two rows and three columns, show the results for three tracking stations (left column, NNOR; middle column, XMIS; right column, GMSD) and two frequencies (upper row, B 1 band; bottom row, B 2 band).
Figure 7. The histograms of C04 SIS UREs (m) after eliminating the trend term covering fourteen experimental days, overlaid with empirical normal probability-density function (PDF) curves (in red). Six subplots, arranged in two rows and three columns, show the results for three tracking stations (left column, NNOR; middle column, XMIS; right column, GMSD) and two frequencies (upper row, B 1 band; bottom row, B 2 band).
Sensors 18 01816 g007
Figure 8. The histograms of C09 SIS UREs (m) after eliminating the trend term covering fourteen experimental days, overlaid with empirical normal probability-density function (PDF) curves (in red). Six subplots, arranged in two rows and three columns, show the results for three tracking stations (left column, NNOR; middle column, XMIS; right column, GMSD) and two frequencies (upper row, B 1 band; bottom row, B 2 band).
Figure 8. The histograms of C09 SIS UREs (m) after eliminating the trend term covering fourteen experimental days, overlaid with empirical normal probability-density function (PDF) curves (in red). Six subplots, arranged in two rows and three columns, show the results for three tracking stations (left column, NNOR; middle column, XMIS; right column, GMSD) and two frequencies (upper row, B 1 band; bottom row, B 2 band).
Sensors 18 01816 g008
Figure 9. The histograms of C14 SIS UREs (m) after eliminating the trend term covering 14 experimental days, overlaid with empirical normal probability-density function (PDF) curves (in red). Six subplots, arranged in two rows and three columns, show the results for three tracking stations (left column, NNOR; middle column, XMIS; right column, GMSD) and two frequencies (upper row, B 1 band; bottom row, B 2 band).
Figure 9. The histograms of C14 SIS UREs (m) after eliminating the trend term covering 14 experimental days, overlaid with empirical normal probability-density function (PDF) curves (in red). Six subplots, arranged in two rows and three columns, show the results for three tracking stations (left column, NNOR; middle column, XMIS; right column, GMSD) and two frequencies (upper row, B 1 band; bottom row, B 2 band).
Sensors 18 01816 g009
Figure 10. Satellite-health flag as a function of BDT, determined using two methods (blue *, BDS Ground SIS Monitoring (BDSGSISM); red dots, broadcast ephemeris) on day 346 of 2015 for BDS satellite C02.
Figure 10. Satellite-health flag as a function of BDT, determined using two methods (blue *, BDS Ground SIS Monitoring (BDSGSISM); red dots, broadcast ephemeris) on day 346 of 2015 for BDS satellite C02.
Sensors 18 01816 g010
Figure 11. SIS-URE after eliminating the trend term time series as a function of BDT on day 346 of 2015 for BDS satellite C02 tracked by NNOR. Subplot (left) corresponds to the tracking data of satellite C02 by NNOR for the whole day, and subplot (right) corresponds to the tracking data of satellite C02 by NNOR with the anomalous periods removed.
Figure 11. SIS-URE after eliminating the trend term time series as a function of BDT on day 346 of 2015 for BDS satellite C02 tracked by NNOR. Subplot (left) corresponds to the tracking data of satellite C02 by NNOR for the whole day, and subplot (right) corresponds to the tracking data of satellite C02 by NNOR with the anomalous periods removed.
Sensors 18 01816 g011
Figure 12. Satellite-health flag as a function of BDT, determined using two methods (blue *, BDSGSISM; red dots, broadcast ephemeris) on day 359 of 2015 for BDS satellite C10.
Figure 12. Satellite-health flag as a function of BDT, determined using two methods (blue *, BDSGSISM; red dots, broadcast ephemeris) on day 359 of 2015 for BDS satellite C10.
Sensors 18 01816 g012
Figure 13. SIS-URE after eliminating the trend term time series as a function of BDT on day 359 of 2015 for BDS satellite C10 tracked by XMIS. Subplot (left) corresponds to the tracking data of satellite C10 by XMIS for the whole day, and subplot (right) corresponds to the tracking data of satellite C10 by XMIS with the anomalous periods removed.
Figure 13. SIS-URE after eliminating the trend term time series as a function of BDT on day 359 of 2015 for BDS satellite C10 tracked by XMIS. Subplot (left) corresponds to the tracking data of satellite C10 by XMIS for the whole day, and subplot (right) corresponds to the tracking data of satellite C10 by XMIS with the anomalous periods removed.
Sensors 18 01816 g013
Table 1. MGEX tracking-station information in this study.
Table 1. MGEX tracking-station information in this study.
Tracking StationsReceiver TypeAntenna TypePRN of Tracking
CAS1Trimble NETR9LEIAR25.R3C01 C03 C06 C07 C08 C09 C10 C11 C12 C14
DJIGTrimble NETR9TRM59800.00C02 C03 C05 C06 C07 C08 C09 C10 C11 C12 C14
GMSDTrimble NETR9TRM59800.00C01 C02 C03 C04 C06 C07 C08 C09 C10 C11 C12 C14
KARRTrimble NETR9TRM59800.00C01 C03 C04 C05 C06 C07 C08 C09 C10 C11 C12 C14
KRGGLeica GR10LEIAR25.R4C02 C03 C05 C06 C07 C08 C09 C10 C11 C12 C14
KZN2Trimble NETR9TRM59800.00C02 C05 C06 C07 C08 C09 C10 C11 C12 C14
JFNGTrimble NETR9TRM59800.00C01 C02 C03 C04 C05 C06 C07 C08 C09 C10 C11 C12 C14
NNORSEPT POLARX4SEPCHOKE_MCC01 C02 C03 C04 C05 C06 C07 C08 C09 C10 C11 C12 C14
POHNTrimble NETR9TRM59800.00C01 C03 C04 C06 C07 C08 C09 C10 C11 C12 C14
PTVLTrimble NETRTRM59800.00C01 C03 C04 C06 C07 C08 C09 C10 C11 C12 C14
REUNTrimble NETR9TRM55971.00C02 C03 C05 C06 C07 C08 C09 C10 C11 C12 C14
UNX3SEPT ASTERX3LEIAR25.R3C01 C03 C04 C06 C07 C08 C09 C10 C11 C12 C14
XMISTrimble NETR9TRM59800.00C01 C03 C04 C05 C06 C07 C08 C09 C10 C11 C12 C14
Table 2. Statistical results of the SIS UREs (m) after eliminating the trend term covering 14 experimental days.
Table 2. Statistical results of the SIS UREs (m) after eliminating the trend term covering 14 experimental days.
StationFrequencyC04C09C14
Mean/mSTD/mMean/mSTD/mMean/mSTD/m
NNORB100.35−0.020.46−0.010.53
B200.210.010.210.010.32
XMISB100.3800.410.010.40
B200.2300.2400.24
GMSDB10.050.3300.370.020.32
B20.060.2000.21−0.010.19
Table 3. List of Potential SIS Anomalies from 1 November 2015 to 27 December 2015.
Table 3. List of Potential SIS Anomalies from 1 November 2015 to 27 December 2015.
SatellitePRNDateStart TimeDuration (Minutes)
G04C046 November 201507:2932
I01C067 November 201503:2341
I02C077 November 201511:1750
M03C118 November 201515:1544
G06C0211 November 201504:0851
G03C0311 November 201507:4021
G04C0414 November 201502:2993
I05C1014 November 201512:0020
I02C0714 November 201515:3047
M03C1114 November 201517:55365
M03C1115 November 201500:01180
G06C0215 November 201505:1842
G06C0219 November 201506:1051
G06C0221 November 201507:0753
G06C0227 November 201501:4429
G06C021 December 201500:2581
I01C062 December 201501:42112
G06C027 December 201520:0852
G03C038 December 201507:5961
G04C048 December 201521:2338
G06C029 December 201519:5763
I01C0611 December 201506:4062
I02C0711 December 201509:0284
G06C0212 December 201511:00120
G06C0213 December 201501:5028
G04C0414 December 201508:2732
G06C0215 December 201507:1248
M04C1215 December 201507:3825
G06C0216 December 201506:3331
G05C0518 December 201504:1543
G06C0220 December 201510:5329
I01C0623 December 201504:4432
M04C1223 December 201507:1023
I02C0725 December 201512:4529
G06C0227 December 201505:2325

Share and Cite

MDPI and ACS Style

Jiang, H.; Wang, H.; Wang, Z.; Yuan, Y. Real-Time Monitoring for BDS Signal-In-Space Anomalies Using Ground Observation Data. Sensors 2018, 18, 1816. https://doi.org/10.3390/s18061816

AMA Style

Jiang H, Wang H, Wang Z, Yuan Y. Real-Time Monitoring for BDS Signal-In-Space Anomalies Using Ground Observation Data. Sensors. 2018; 18(6):1816. https://doi.org/10.3390/s18061816

Chicago/Turabian Style

Jiang, Hu, Haitao Wang, Zemin Wang, and Yunbin Yuan. 2018. "Real-Time Monitoring for BDS Signal-In-Space Anomalies Using Ground Observation Data" Sensors 18, no. 6: 1816. https://doi.org/10.3390/s18061816

APA Style

Jiang, H., Wang, H., Wang, Z., & Yuan, Y. (2018). Real-Time Monitoring for BDS Signal-In-Space Anomalies Using Ground Observation Data. Sensors, 18(6), 1816. https://doi.org/10.3390/s18061816

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