Characteristics of BDS Signal-in-Space User Ranging Errors and Their Effect on Advanced Receiver Autonomous Integrity Monitoring Performance

Signal-In-Space User Range Errors (SIS UREs) are assumed to be overbounded by a normal distribution with a standard deviation represented by the User Range Accuracy (URA). The BeiDou Navigation Satellite System (BDS) broadcast URA is not compatible with the historical SIS URE performance that affects the Advanced Receiver Autonomous Integrity Monitoring (ARAIM) False Alert Probability (Pfa) and availability evaluation. This study compares the BDS broadcast and precise ephemeris from 1 March 2013 to 1 March 2017 to obtain SIS UREs. Through analyzing the statistical characteristics of the SIS UREs, we obtain the standard deviation σURE for the accuracy and continuity and σURA used for the integrity of the SIS UREs. The results show that the broadcast σURA of 2 m cannot completely overbound SIS UREs for all BDS satellites, but the σURA of 2.4 m can. Then, we use the σURA of 2.4 m to evaluate the ARAIM Pfa and availability. The results show that the Pfa may increase to 2 × 10−5 and exceed its limit by an order of magnitude. We also consider the differences between the SIS UREs of Geostationary Earth Orbit (GEO), Inclined Geo-Synchronous Orbit (IGSO), and Medium Earth Orbit (MEO). The results indicate that all Pfa values calculated by the computed σURE are less than the Pfa in the Integrity Support Message (ISM) for the worst-performing GEO satellite. The approximately 55% Pfa calculated by the computed σURE is less than the Pfa in ISM for the worst-performing IGSO satellite. Most Pfa values calculated by the computed σURE is less than the Pfa in the ISM for the worst-performing MEO satellite. For BDS satellites, the Pfa is mainly affected by σURE. When the σURA of 2.4 m is used to evaluate the availability, the computed availability is lower than the availability calculated by the broadcast σURA/σURE and the greatest degradation can reach 25%.


Introduction
Advanced Receiver Autonomous Integrity Monitoring (ARAIM) can provide global users with integrity monitoring service for vertical guidance [1]. The ARAIM algorithm is affected by several error sources. These errors can be divided into errors that arise in navigation satellites, affect signals during propagation, and emerge from the receiver or its surrounding environment [2]. This paper focus on analyzing the errors that arise in navigation satellites, such errors are called Signal-In-Space Errors (SISEs).
SISEs are derived by comparing the broadcast ephemeris with the precise ephemeris that includes satellite orbit and clock errors, signal generation errors and antenna characteristics. The satellite orbit and clock errors are the largest sources of SISEs because the broadcast satellite location and clock states

Statistics and Overbound of BDS SIS URE
The section calculates the BDS SIS UREs and analyzes the overbound of the SIS UREs to obtain the value of the URA.

Data Preprocessing
When comparing the BDS broadcast ephemeris with the precise ephemeris to obtain the SIS UREs, it is necessary to use the useful data and remove the offsets between systems. This section filters out the outliers to obtain effective data and removes the time offset to correct the clock errors.

Data Sources
The broadcast ephemeris file used in this study is in the RINEX V3 format and is generated by the center of the Multi-GNSS Experiment (MGEX). The precise ephemeris contains accurate satellite orbital elements and clock error information. We consider that the precise ephemeris can replace the real satellite orbital elements to calculate the SIS UREs because the precise ephemeris and the real satellite orbital elements are of the same order of magnitude, with errors of only a few centimeters; thus, the precise ephemeris is much more accurate than the broadcast ephemeris [11]. The precise ephemeris is generated by the International GNSS Service (IGS) and MGEX Analysis Center and can be obtained from a public database.

BDS Outlier Filter
Not all data can be statistically analyzed when calculating the SIS UREs. For example, the GPS SPS Performance Standard [12] has defined that SIS UREs must be below 4.42 times URA, and that users cannot use a broadcast ephemeris that is older than 4 h or corrupted. Therefore, the SIS UREs that satisfy any of the following conditions in this analysis are rejected [5]: • The integrity status word is not zero, or the URA is greater than 48 m. • The broadcast ephemeris is not within the 4 h period.

•
The precise ephemeris/clock is lost.

•
The SIS UREs exceed the URA by a factor of more than 4.42.
Long-term stable data are obtained by filtering out outliers according to the abovementioned criteria. Figure 1 shows the status of the 4-year BDS ephemeris. The y-axis represents the Pseudo Random Noise (PRN) and Space Vehicle Number (SVN). Green indicates that both the precise ephemeris and broadcast ephemeris are available. Red indicates that the precise ephemeris is lost but that the broadcast ephemeris is obtained. Purple indicates that the broadcast ephemeris is lost but that the precise ephemeris is available. Gray indicates that both the broadcast ephemeris and precise ephemeris are unavailable. Blue indicates that the broadcast ephemeris is corrupted. Black indicates abnormal clock errors. Satellite C13 is special among all the satellites because there are no ephemeris data nearly half the time. On 20 March 2014, the precise ephemeris of satellite C13 was missing. On 22 October 2014, the broadcast ephemeris of satellite C13 was missing. On 11 October 2016, both the broadcast ephemeris and precise ephemeris can be obtained simultaneously. The satellite MEO05 was replaced by satellite IGSO06 in January 2016, so C13 can represent both satellites. Because of the unique characteristics of C13, this satellite is excluded when calculating the anomalous probability of ephemeris data. There are 455,520 ephemeris data points for the other 13 satellites. There are 23,868 anomalous data points; thus, 5.24% of the ephemeris data are filtered out, and the anomalous data in which the precise ephemeris is lost but there is a broadcast ephemeris appears more than the other types of anomalous data. thus, the precise ephemeris is much more accurate than the broadcast ephemeris [11]. The precise ephemeris is generated by the International GNSS Service (IGS) and MGEX Analysis Center and can be obtained from a public database.

BDS Outlier Filter
Not all data can be statistically analyzed when calculating the SIS UREs. For example, the GPS SPS Performance Standard [12] has defined that SIS UREs must be below 4.42 times URA, and that users cannot use a broadcast ephemeris that is older than 4 h or corrupted. Therefore, the SIS UREs that satisfy any of the following conditions in this analysis are rejected [5]: • The integrity status word is not zero, or the URA is greater than 48 m.

•
The broadcast ephemeris is not within the 4 h period.

•
The precise ephemeris/clock is lost.

•
The SIS UREs exceed the URA by a factor of more than 4.42.
Long-term stable data are obtained by filtering out outliers according to the abovementioned criteria. Figure 1 shows the status of the 4-year BDS ephemeris. The y-axis represents the Pseudo Random Noise (PRN) and Space Vehicle Number (SVN). Green indicates that both the precise ephemeris and broadcast ephemeris are available. Red indicates that the precise ephemeris is lost but that the broadcast ephemeris is obtained. Purple indicates that the broadcast ephemeris is lost but that the precise ephemeris is available. Gray indicates that both the broadcast ephemeris and precise ephemeris are unavailable. Blue indicates that the broadcast ephemeris is corrupted. Black indicates abnormal clock errors. Satellite C13 is special among all the satellites because there are no ephemeris data nearly half the time. On 20 March 2014, the precise ephemeris of satellite C13 was missing. On 22 October 2014, the broadcast ephemeris of satellite C13 was missing. On October 11, 2016, both the broadcast ephemeris and precise ephemeris can be obtained simultaneously. The satellite MEO05 was replaced by satellite IGSO06 in January 2016, so C13 can represent both satellites. Because of the unique characteristics of C13, this satellite is excluded when calculating the anomalous probability of ephemeris data. There are 455,520 ephemeris data points for the other 13 satellites. There are 23,868 anomalous data points; thus, 5.24% of the ephemeris data are filtered out, and the anomalous data in which the precise ephemeris is lost but there is a broadcast ephemeris appears more than the other types of anomalous data

Statistics of BDS SIS URE
where R represents the radial errors, A represents the along-track errors, and C represents the cross-track errors. R, A and C are the three directions of the orbit errors. T represents clock errors. • BDS orbit-only rms SIS UREs. Orbit-only rms SIS UREs are obtained by removing the clock errors according to Equations (1) and (2): • BDS worst-case SIS UREs. The worst-case SIS UREs are the maximum instantaneous SIS UREs calculated by the broadcast ephemeris. The worst-case SIS UREs play a crucial role in the SIS performance evaluation. The worst-case SIS UREs can be calculated either from the instantaneous SIS UREs or from radial errors, along-track errors, cross-track errors and clock errors. The present study uses the latter approach [14]: where γ is the latitude of the edge of satellite coverage.

Analysis of BDS Clock Errors, Orbit Errors and SIS UREs
The time system of the BDS broadcast ephemeris is BeiDou Time (BDT), and the time system of the precise ephemeris is GPS Time (GPST); thus, the time systems must be unified when comparing the broadcast ephemeris with the precise ephemeris. Therefore, we remove the BDT-GPST bias of 14 s [10]. The broadcast ephemeris is the reference single-frequency signal B3I, maintained by the BDS military time-frequency. However, the BDS precision clock is produced by Wuhan University, calculated using dual frequency (B1I and B2I) ionosphere-free observations that use a monitoring station receiver clock as a reference [15]. When calculating the positions of satellites according to the broadcast ephemeris, there are Time Group Delays (TGDs) in signal B1I and signal B2I. The TGDsT GD1 and T GD2 , which are less than one nanosecond, should be removed. In addition, the different timescales of the broadcast ephemeris and precise ephemeris cause a time offset that should be considered [16].
When calculating the orbit errors, the BDS broadcast ephemeris uses the China Geodetic Coordinate System 2000 (CGCS 2000), and the precise ephemeris uses the International Terrestrial Reference Frame (ITRF) when calculating orbit and clock errors. The differences between the CGCS 2000 and ITRF coordinate systems are ignored [17]. The precise ephemeris provides the center-of-mass coordinates of the satellites, and the broadcast ephemeris refers to the antenna phase center. However, the BDS broadcast ephemeris refers to a position close to the center-of-mass [18]. Therefore, the antenna phase center offsets do not have a significant effect and can be ignored.
The daily average clock errors, orbit errors and SIS UREs of the GEO, IGSO and MEO satellites are respectively analyzed using  These figures show that the clock errors of the GEO and IGSO satellites change significantly, and the maximum value of the clock errors can exceed 5 m. However, the trend of the clock errors of the MEO satellites are relatively stable. Therefore, the clock stability of the MEO satellites is better than that of the IGSO and GEO satellites. Comparing the orbit errors of GEO, IGSO and MEO satellites, it can be found that the orbit errors of the GEO satellites are the worst, especially the C04 satellite.  When calculating the orbit errors, the BDS broadcast ephemeris uses the China Geodetic Coordinate System 2000 (CGCS 2000), and the precise ephemeris uses the International Terrestrial Reference Frame (ITRF) when calculating orbit and clock errors. The differences between the CGCS 2000 and ITRF coordinate systems are ignored [17]. The precise ephemeris provides the center-of-mass coordinates of the satellites, and the broadcast ephemeris refers to the antenna phase center. However, the BDS broadcast ephemeris refers to a position close to the center-of-mass [18]. Therefore, the antenna phase center offsets do not have a significant effect and can be ignored.
The daily average clock errors, orbit errors and SIS UREs of the GEO, IGSO and MEO satellites are respectively analyzed using Figures 2, 3 and 4.
These figures show that the clock errors of the GEO and IGSO satellites change significantly, and the maximum value of the clock errors can exceed 5 m. However, the trend of the clock errors of the MEO satellites are relatively stable. Therefore, the clock stability of the MEO satellites is better than that of the IGSO and GEO satellites. Comparing the orbit errors of GEO, IGSO and MEO satellites, it can be found that the orbit errors of the GEO satellites are the worst, especially the C04 satellite.    When calculating the orbit errors, the BDS broadcast ephemeris uses the China Geodetic Coordinate System 2000 (CGCS 2000), and the precise ephemeris uses the International Terrestrial Reference Frame (ITRF) when calculating orbit and clock errors. The differences between the CGCS 2000 and ITRF coordinate systems are ignored [17]. The precise ephemeris provides the center-of-mass coordinates of the satellites, and the broadcast ephemeris refers to the antenna phase center. However, the BDS broadcast ephemeris refers to a position close to the center-of-mass [18]. Therefore, the antenna phase center offsets do not have a significant effect and can be ignored.
The daily average clock errors, orbit errors and SIS UREs of the GEO, IGSO and MEO satellites are respectively analyzed using Figures 2, 3 and 4.
These figures show that the clock errors of the GEO and IGSO satellites change significantly, and the maximum value of the clock errors can exceed 5 m. However, the trend of the clock errors of the MEO satellites are relatively stable. Therefore, the clock stability of the MEO satellites is better than that of the IGSO and GEO satellites. Comparing the orbit errors of GEO, IGSO and MEO satellites, it can be found that the orbit errors of the GEO satellites are the worst, especially the C04 satellite.

Overbound BDS SIS UREs
SIS UREs are usually assumed to be overbounded by a normal distribution with a standard deviation represented by the URA. The URA is an estimate of the UREs caused by the satellite section, which requires overbound SIS UREs with a probability of 68% [19]. The URA is used to describe the accuracy of satellite space signals and is characterized by the User Range Accuracy Index (URAI). The broadcast URA is a relatively conservative estimate of the standard deviation of SIS UREs and is mainly used for integrity monitoring.
The BDS broadcast ephemeris contains a URAI with a value of 2. When the URAI is converted to a URA, the URA is a value range instead of a fixed value. According to the CSNO (2016) BeiDou navigation satellite system signal in space interface control document [20], when 0 < N ≤ 6, URA = 2 N/2+1 . Therefore, when N = 0, 0 < URA ≤ 2.4. Because the broadcast URA is 2 m, the upper URA in the interface control document is 2.4 m.
The cumulative distribution functions of the worst-case UREs are shown in Figure 5. The blue dashed line shows that the URA is 2 m, and the red dashed line shows that the URA is 2.4 m. The black dashed line denotes the 68% overbounds. The required overbound of satellite C01 is the largest and that of satellite C07 is the smallest. When using a broadcast URA of 2 m, the overbound of satellite C01 is 61.47%, which is below 68% and that of satellite C07 is 88.02%. When using the upper URA of 2.4 m, the overbound of satellite C01 is 78.68%, and that of satellite C07 is 92.62%.

Overbound BDS SIS UREs
SIS UREs are usually assumed to be overbounded by a normal distribution with a standard deviation represented by the URA. The URA is an estimate of the UREs caused by the satellite section, which requires overbound SIS UREs with a probability of 68% [19]. The URA is used to describe the accuracy of satellite space signals and is characterized by the User Range Accuracy Index (URAI). The broadcast URA is a relatively conservative estimate of the standard deviation of SIS UREs and is mainly used for integrity monitoring.
The BDS broadcast ephemeris contains a URAI with a value of 2. When the URAI is converted to a URA, the URA is a value range instead of a fixed value. According to the CSNO (2016) BeiDou navigation satellite system signal in space interface control document [20], when 0 < N ≤ 6, URA = 2 N/2+1 . Therefore, when N = 0, 0 < URA ≤ 2.4. Because the broadcast URA is 2 m, the upper URA in the interface control document is 2.4 m.
The cumulative distribution functions of the worst-case UREs are shown in Figure 5. The blue dashed line shows that the URA is 2 m, and the red dashed line shows that the URA is 2.4 m. The black dashed line denotes the 68% overbounds. The required overbound of satellite C01 is the largest and that of satellite C07 is the smallest. When using a broadcast URA of 2 m, the overbound of satellite C01 is 61.47%, which is below 68% and that of satellite C07 is 88.02%. When using the upper URA of 2.4 m, the overbound of satellite C01 is 78.68%, and that of satellite C07 is 92.62%.

Overbound BDS SIS UREs
SIS UREs are usually assumed to be overbounded by a normal distribution with a standard deviation represented by the URA. The URA is an estimate of the UREs caused by the satellite section, which requires overbound SIS UREs with a probability of 68% [19]. The URA is used to describe the accuracy of satellite space signals and is characterized by the User Range Accuracy Index (URAI). The broadcast URA is a relatively conservative estimate of the standard deviation of SIS UREs and is mainly used for integrity monitoring.
The BDS broadcast ephemeris contains a URAI with a value of 2. When the URAI is converted to a URA, the URA is a value range instead of a fixed value. According to the CSNO (2016) BeiDou navigation satellite system signal in space interface control document [20], when 0 < N ≤ 6, URA = 2 N/2+1 . Therefore, when N = 0, 0 < URA ≤ 2.4. Because the broadcast URA is 2 m, the upper URA in the interface control document is 2.4 m.
The cumulative distribution functions of the worst-case UREs are shown in Figure 5. The blue dashed line shows that the URA is 2 m, and the red dashed line shows that the URA is 2.4 m. The black dashed line denotes the 68% overbounds. The required overbound of satellite C01 is the largest and that of satellite C07 is the smallest. When using a broadcast URA of 2 m, the overbound of satellite C01 is 61.47%, which is below 68% and that of satellite C07 is 88.02%. When using the upper URA of 2.4 m, the overbound of satellite C01 is 78.68%, and that of satellite C07 is 92.62%.   Table 1 lists the 68% overbound worst-case UREs. Satellite C01 has the largest 68% overbound, and satellite C07 has the smallest. Therefore, the upper URA of 2.4 m can overbound the worst-case UREs of the 14 BDS satellites. However, the broadcast URA of 2 m cannot overbound the SIS UREs of all BDS satellites such as satellites C01, C02, C03, and C14. Therefore, the broadcast URA is not sufficiently conservative, which is inconsistent with the performance of the worst-case UREs.

Effect of BDS URA on ARAIM
According to the previous section, the BDS broadcast URA is not consistent with the historical performance of SIS UREs, but the URA of 2.4 m can satisfy the requirement. As a result, this section analyzes the Pfa and ARAIM availability when using URAs of 2.4 m and 2 m, respectively.

Effect of BDS URA on Pfa
The ARAIM algorithm relies on the ground monitoring network to observe the state of the satellite and provide this information to aircraft in the form of an ISM to meet the requirements of system integrity and continuity [21]. ISM parameters (either online or offline) can be divided into two types: fault probability and nominal errors. The first group of parameters includes the probability of a satellite fault P sat and the probability of a constellation fault P const . The second group of parameters includes the nominal SISE, which provides the user with three parameters used to overbound nominal pseudorange errors: σ URE , σ URA and the maximum nominal bias b nom [22]. The values of σ URE and σ URA have an effect on the calculation of Pfa.
According to the ARAIM solution separation [23], the coefficient of Pfa allocated to the vertical and horizontal directions are derived as follows: where σ ss,q (k) is the standard deviation of the difference between the all-in-view and the fault tolerant position solutions. T k,q is the threshold of the solution separation test. P FA_HOR is the horizontal Pfa, and P FA_VERT is the vertical Pfa. k represents the k-th failure mode. N f ault_modes is the number of fault modes. The indices q = 1, 2 and 3 designate the east, north and up components, respectively. The standard deviation σ ss,q (k) is calculated by: where G is the geometry matrix in terms of the east, north, and up components with a clock component for each constellation. W is the weighting matrix. C acc is the pseudorange error diagonal covariance matrix for accuracy and continuity: where j is the index of the constellation, and idx k is the indices of the satellites included in fault subset k. C int (j, j) is the pseudorange error diagonal covariance matrix for integrity: where σ tropo,i and σ user,i denote the tropospheric delays and user errors, respectively. They are related to the elevation angle: the calculation model for σ tropo,i is [1]: The calculation model for σ user,i is [1]: where θ is the elevation angle in degrees, f B1I and f B2I are the frequencies of signal B1 and signal B2, σ MP is the multipath error, and σ Noise is the receiver noise.
According to the above formula, the Pfa can be derived as follows: where Q represents the cumulative distribution function of a zero mean unit Gaussian distribution.
To verify the relationship between the σ URA /σ URE and Pfa, we choose 14 BDS satellites for the ARAIM simulation. If only single satellite faults and single constellation faults are considered, then the total number of N fault_modes = 15. When the solution separation test is conducted, the east, north and up components must be considered; thus, the total fault detection quantity is 45. Because the number of tests is enough, it is common to find at least one test for which T k,q is relatively small. The T k,q is the threshold of the solution separation test. If the solution exceeds the threshold, the result of the detection test can be considered a false alert. In this case, the P fa,q (k) changes as σ URE varies. As shown in Figure 6, the change in σ URE causes the probability density function to change, and the Pfa changes if the threshold T k,q is constant. The black and green lines respectively represent the computed UREs that cause the computed σ URE to be less than the broadcast σ URE and greater than the broadcast σ URE , the blue line represents the broadcast URE, and the red region represents the Pfa. When the threshold T k,q is constant, the red area corresponding to the black line is smaller than that corresponding to the  According to the previous section, the performances of the GEO, IGSO and MEO satellites vary. SIS UREs of different satellites can be considered to be uncorrelated in ARAIM [6]. Therefore, the standard deviations σURE of satellites C01, C06 and C14 are selected to consider the change in the Pfa in one day. The SIS UREs are selected to calculate the Pfa on March 1, 2017. The time interval is 1 s, and the total number of samples is 86,400. We take the sample value per minute to obtain σURE.
Because a single satellite is considered, each satellite undergoes a period during each day when it is not within the field of view of the receiver at Wuhan University. Pfa changes only within the field   According to the previous section, the performances of the GEO, IGSO and MEO satellites vary. SIS UREs of different satellites can be considered to be uncorrelated in ARAIM [6]. Therefore, the standard deviations σURE of satellites C01, C06 and C14 are selected to consider the change in the Pfa in one day. The SIS UREs are selected to calculate the Pfa on March 1, 2017. The time interval is 1 s, and the total number of samples is 86,400. We take the sample value per minute to obtain σURE.
Because a single satellite is considered, each satellite undergoes a period during each day when it is not within the field of view of the receiver at Wuhan University. Pfa changes only within the field According to the previous section, the performances of the GEO, IGSO and MEO satellites vary. SIS UREs of different satellites can be considered to be uncorrelated in ARAIM [6]. Therefore, the standard deviations σ URE of satellites C01, C06 and C14 are selected to consider the change in the Pfa in one day. The SIS UREs are selected to calculate the Pfa on 1 March 2017. The time interval is 1 s, and the total number of samples is 86,400. We take the sample value per minute to obtain σ URE .
Because a single satellite is considered, each satellite undergoes a period during each day when it is not within the field of view of the receiver at Wuhan University. Pfa changes only within the field of view of the satellite are considered, which means that changes in the satellite elevation angle must be considered. The mask angle is set to 5 • , and the results are shown in Figure 8. The blue line is the elevation angle of satellite C01, which is 43.4 • -46.1 • because the relative position between the GEO satellites and ground remains unchanged. Therefore, the GEO satellites are readily visible. The black line is the elevation angle of satellite C06, which is 5 • -79.7 • . The elevation angle is below the defilade angle for the 4 h between 9 and 13 h and is thus invisible. The elevation transformation of satellite C06 is approximately symmetrical because the orbit of the IGSO satellite is an "8" shaped orbit, and the 24 h trajectory under the satellite oscillates between north and south within the service area. The green line is the elevation angle of satellite C14, which is 5 • -84.9 • . The satellite is visible only between 0 and 6 h and between 11:50 and 12:50 because the 24 h trajectory under the satellite moves globally.  As the satellite elevation angle changes, the tropospheric delays σtropo,i and user errors σuser,i also change. Figures 9, 10 and 11 show these parameters for satellites C01, C06 and C14, respectively. We use the Beijing Time in one day as the x-axis.  As the satellite elevation angle changes, the tropospheric delays σ tropo,i and user errors σ user,i also change. Figures 9-11 show these parameters for satellites C01, C06 and C14, respectively. We use the Beijing Time in one day as the x-axis.  As the satellite elevation angle changes, the tropospheric delays σtropo,i and user errors σuser,i also change. Figures 9, 10 and 11 show these parameters for satellites C01, C06 and C14, respectively. We use the Beijing Time in one day as the x-axis.     Figure 10 shows that satellite C06 has a σtropo,i of 0.13-1.21 m. At 0-5 h and 17-24 h, the value of σtropo,i changes slowly and approximates a straight line because the elevation angle in these two periods is large. For the same reason, the user errors σuser,i approximate a straight line at 0-6 h and 16-24 h. Figure 11 shows that satellite C14 has a σtropo,i of 0.12-1.226 m and a σuser,i of 0.64-1.85 m. The values of σtropo,i and σuser,i of satellites C06 and C14 are similar but greater than those of C01 because the elevation angle of the GEO satellites is larger and substantially unchanged relative to the receiver. In addition, the value of σuser,i for all three satellites is greater than σtropo,i, which indicates that the effect of σuser,i on the Pfa is greater than the effect of σtropo,i.
The changes in σtropo,i, σuser,i and σURE on March 1, 2017, are used to derive the Pfa. The results are shown in Figure 12 (X-axis is the Beijing Time in one day): satellite C01 has a Pfa of 1.09 × 10 −6 -2.76 × 10 −6 , satellite C06 has a Pfa of 9.38 × 10 −7 -2.70 × 10 −5 , and satellite C14 has a Pfa of 8.45 × 10 −7 -8.73 × 10 −6 . We call the Pfa calculated by the computed σURE the computed Pfa. The computed Pfa of satellite    Figure 10 shows that satellite C06 has a σtropo,i of 0.13-1.21 m. At 0-5 h and 17-24 h, the value of σtropo,i changes slowly and approximates a straight line because the elevation angle in these two periods is large. For the same reason, the user errors σuser,i approximate a straight line at 0-6 h and 16-24 h. Figure 11 shows that satellite C14 has a σtropo,i of 0.12-1.226 m and a σuser,i of 0.64-1.85 m. The values of σtropo,i and σuser,i of satellites C06 and C14 are similar but greater than those of C01 because the elevation angle of the GEO satellites is larger and substantially unchanged relative to the receiver. In addition, the value of σuser,i for all three satellites is greater than σtropo,i, which indicates that the effect of σuser,i on the Pfa is greater than the effect of σtropo,i.
The changes in σtropo,i, σuser,i and σURE on March 1, 2017, are used to derive the Pfa. The results are shown in Figure 12 (X-axis is the Beijing Time in one day): satellite C01 has a Pfa of 1.09 × 10 −6 -2.76 × 10 −6 , satellite C06 has a Pfa of 9.38 × 10 −7 -2.70 × 10 −5 , and satellite C14 has a Pfa of 8.45 × 10 −7 -8.73 × 10 −6 . We call the Pfa calculated by the computed σURE the computed Pfa. The computed Pfa of satellite Figure 11. The tropospheric delays σ tropo,i , user errors σ user,i and σ URE of satellite C14 over one day as seen from the receiver at Wuhan University. Figure 9 shows that satellite C01 has a σ tropo,i of 0.165-0.175 and a σ user,i of 0.65-0.66 m. The value of σ URE is greater than σ tropo,i and σ user,i at all times for satellite C01. Figure 10 shows that satellite C06 has a σ tropo,i of 0.13-1.21 m. At 0-5 h and 17-24 h, the value of σ tropo,i changes slowly and approximates a straight line because the elevation angle in these two periods is large. For the same reason, the user errors σ user,i approximate a straight line at 0-6 h and 16-24 h. Figure 11 shows that satellite C14 has a σ tropo,i of 0.12-1.226 m and a σ user,i of 0.64-1.85 m. The values of σ tropo,i and σ user,i of satellites C06 and C14 are similar but greater than those of C01 because the elevation angle of the GEO satellites is larger and substantially unchanged relative to the receiver. In addition, the value of σ user,i for all three satellites is greater than σ tropo,i , which indicates that the effect of σ user,i on the Pfa is greater than the effect of σ tropo,i .
The changes in σ tropo,i , σ user,i and σ URE on 1 March 2017, are used to derive the Pfa. The results are shown in Figure 12 (X-axis is the Beijing Time in one day): satellite C01 has a Pfa of 1.09 × 10 −6 -2.76 × 10 −6 , satellite C06 has a Pfa of 9.38 × 10 −7 -2.70 × 10 −5 , and satellite C14 has a Pfa of 8.45 × 10 −7 -8.73 × 10 −6 . We call the Pfa calculated by the computed σ URE the computed Pfa. The computed Pfa of satellite C01 is less than the Pfa in the ISM because the computed σ URE is less than the broadcast σ URE . The degradation of Pfa may reach 62.5%. Approximately 55% of the computed Pfa values of satellite C06 are less than the Pfa in the ISM and the largest drop reaches 76.6% because the computed σ URE is close to 0 m. The other computed Pfa values are greater than the Pfa in the ISM because the computed σ URE is close to 4 m. Most of the computed Pfa values of satellite C14 are smaller than the Pfa in the ISM and the greatest drop reaches 78.9%. The computed Pfa drops sharply at approximately 3 h because the computed σ URE decreases to 0 m. The computed Pfa rapidly increases and ultimately exceeds the Pfa in the ISM at 3-6 h, which is related to σ tropo,i and σ user,i . Therefore, the Pfa may be mainly affected by σ URE but also partly affected by σ tropo,i and σ user,i . In the case of a small σ URE , it is possible that σ tropo,i and σ user,i strongly affect the Pfa.

Effect of BDS URA on ARAIM Availability
In the section, the computed σURA of 2.4 m and computed σURE of 1.6 m are used to evaluate the availability. GPS satellites and 14 BDS satellites are used to perform an ARAIM simulation. The parameters in the ISM are set as shown in Table 2.  Figure 13 shows the Vertical Protection Levels (VPLs) and Horizontal Protection Levels (HPLs) with simulation points for the computed σURA/σURE and the broadcast σURA/σURE (960 data points in each figure). The VPLs/HPLs calculated by the computed σURA/σURE are referred to as the actual VPLs/HPLs and the VPLs/HPLs calculated by the broadcast σURA/σURE are referred to as the broadcast VPLs/HPLs. Blue points represent the actual VPLs/HPLs, green points represent the broadcast VPLs/HPLs, and red lines represent the alert values for the HPL and VPL. The actual VPLs and HPLs are greater than the broadcast values. There are 340 actual VPLs and 209 broadcast VPLs exceeding the Vertical Alert limit (VAL); only 13 actual HPLs exceed the Horizontal Alert limit (HAL), and all broadcast HPLs are below the HAL. Therefore, in some cases, the Protection Levels (PLs) calculated by the computed σURA/σURE exceed the Alert Level (AL), but those calculated by the broadcast σURA/σURE are below the AL, which may lead to missing alert. The increase in the number of missing alerts may lead to integrity degradation.

Effect of BDS URA on ARAIM Availability
In the section, the computed σ URA of 2.4 m and computed σ URE of 1.6 m are used to evaluate the availability. GPS satellites and 14 BDS satellites are used to perform an ARAIM simulation. The parameters in the ISM are set as shown in Table 2.  Figure 13 shows the Vertical Protection Levels (VPLs) and Horizontal Protection Levels (HPLs) with simulation points for the computed σ URA /σ URE and the broadcast σ URA /σ URE (960 data points in each figure). The VPLs/HPLs calculated by the computed σ URA /σ URE are referred to as the actual VPLs/HPLs and the VPLs/HPLs calculated by the broadcast σ URA /σ URE are referred to as the broadcast VPLs/HPLs. Blue points represent the actual VPLs/HPLs, green points represent the broadcast VPLs/HPLs, and red lines represent the alert values for the HPL and VPL. The actual VPLs and HPLs are greater than the broadcast values. There are 340 actual VPLs and 209 broadcast VPLs exceeding the Vertical Alert limit (VAL); only 13 actual HPLs exceed the Horizontal Alert limit (HAL), and all broadcast HPLs are below the HAL. Therefore, in some cases, the Protection Levels (PLs) calculated by the computed σ URA /σ URE exceed the Alert Level (AL), but those calculated by the broadcast σ URA /σ URE are below the AL, which may lead to missing alert. The increase in the number of missing alerts may lead to integrity degradation.  Figure 14 shows the degradation of the availability for the computed σURA/σURE relative to the broadcast σURA/σURE. The greatest degradation of the availability is approximately 25%. The degradation of the availability for the southern region in China is within 10% and is smaller than that for the northern region because the trajectories of the five GEO satellites are fixed points on the equator. For these five satellites, there is a relatively high elevation angle at low latitudes, while the elevation angle at high latitudes is relatively low. The Geometric Dilution of Precision (GDOP) of the GEO satellite changes with the variation in the latitude, the value of the GDOP at low latitudes is smaller than that at high latitudes [25].

Conclusion
To study the effect of the computed σURA/σURE on the Pfa and availability, we obtained the corresponding SIS UREs by comparing the broadcast ephemeris with the precise ephemeris from March 1, 2013 to March 1, 2017. For the broadcast ephemeris and precise ephemeris, 23,868 outliers were statistically removed to obtain useful data. Then, the overbound of the SIS UREs is analyzed. The result shows that the broadcast URA of 2 m does not satisfy the performance requirements of  Figure 14 shows the degradation of the availability for the computed σ URA /σ URE relative to the broadcast σ URA /σ URE . The greatest degradation of the availability is approximately 25%. The degradation of the availability for the southern region in China is within 10% and is smaller than that for the northern region because the trajectories of the five GEO satellites are fixed points on the equator. For these five satellites, there is a relatively high elevation angle at low latitudes, while the elevation angle at high latitudes is relatively low. The Geometric Dilution of Precision (GDOP) of the GEO satellite changes with the variation in the latitude, the value of the GDOP at low latitudes is smaller than that at high latitudes [25].  Figure 14 shows the degradation of the availability for the computed σURA/σURE relative to the broadcast σURA/σURE. The greatest degradation of the availability is approximately 25%. The degradation of the availability for the southern region in China is within 10% and is smaller than that for the northern region because the trajectories of the five GEO satellites are fixed points on the equator. For these five satellites, there is a relatively high elevation angle at low latitudes, while the elevation angle at high latitudes is relatively low. The Geometric Dilution of Precision (GDOP) of the GEO satellite changes with the variation in the latitude, the value of the GDOP at low latitudes is smaller than that at high latitudes [25].

Conclusion
To study the effect of the computed σURA/σURE on the Pfa and availability, we obtained the corresponding SIS UREs by comparing the broadcast ephemeris with the precise ephemeris from March 1, 2013 to March 1, 2017. For the broadcast ephemeris and precise ephemeris, 23,868 outliers were statistically removed to obtain useful data. Then, the overbound of the SIS UREs is analyzed. The result shows that the broadcast URA of 2 m does not satisfy the performance requirements of

Conclusions
To study the effect of the computed σ URA /σ URE on the Pfa and availability, we obtained the corresponding SIS UREs by comparing the broadcast ephemeris with the precise ephemeris from 1 March 2013 to 1 March 2017. For the broadcast ephemeris and precise ephemeris, 23,868 outliers were statistically removed to obtain useful data. Then, the overbound of the SIS UREs is analyzed. The result shows that the broadcast URA of 2 m does not satisfy the performance requirements of observed historical SIS UREs and is slightly optimistic, but the upper URA of 2.4 m is more suitable for users. When using the computed σ URA of 2.4 m and σ URE of 1.6 m to calculate the Pfa, the computed Pfa may exceed the Pfa in the ISM by an order of magnitude. An increase in the Pfa will lead to a degradation of the continuity. Considering the effect of σ URE on the Pfa for GEO, IGSO and MEO satellites, the three types of satellites have different effects on the Pfa. The result shows that the computed Pfa values of satellite C01 are less than the Pfa in ISM, approximately 55% of the computed Pfa values of satellite C06 are less than the Pfa in ISM and most of the computed Pfa values of C14 are less than the Pfa in the ISM. The effects of σ tropo,i and σ user,i are less than that of σ URE , but they may greatly affect the Pfa in the case of a small σ URE . The PLs calculated by the computed σ URA /σ URE are larger than those calculated by the broadcast σ URA /σ URE , sometimes even beyond the AL which can lead to missing alerts. The integrity will degrade because of increases in the number of missing alerts. The availability of the computed σ URA /σ URE is less than that of the broadcast σ URA /σ URE . The greatest degradation of the availability is approximately 25%. The degradation of the availability for the southern region in China is less than that for the northern region.
Author Contributions: Z.W. and W.S. conceived of and designed the research. Z.W. and W.S. performed the experiments and analyzed the results. Z.W. and W.S. wrote the paper. R.L. and T.L. were responsible for data acquisitions, data processing. D.S. was responsible for project managing.