Next Article in Journal / Special Issue
Adjustment of Transceiver Lever Arm Offset and Sound Speed Bias for GNSS-Acoustic Positioning
Previous Article in Journal
An End-to-End Conditional Random Fields and Skip-Connected Generative Adversarial Segmentation Network for Remote Sensing Images
Previous Article in Special Issue
On the Study of Influences of Different Factors on the Rapid Tropospheric Tomography
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Priori Solar Radiation Pressure Model for BeiDou-3 MEO Satellites

College of Geology Engineering and Geomatics, Chang’an University, Xi’an 710054, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2019, 11(13), 1605; https://doi.org/10.3390/rs11131605
Submission received: 23 May 2019 / Revised: 3 July 2019 / Accepted: 4 July 2019 / Published: 5 July 2019
(This article belongs to the Special Issue Global Navigation Satellite Systems for Earth Observing System)

Abstract

:
Due to the cuboid satellite body of BeiDou-3 satellites, the accuracy of their orbit showed a trend of systematic variation with the sun-satellite-earth angle (ε) using the Extend CODE Orbit Model (ECOM1). Therefore, an a priori cuboid box-wing model (named the cuboid model) is necessary to compensate ECOM1. Considering that the body-dimensions and optical properties of the BeiDou-3 satellites used to construct the box-wing model have not yet been fully released, the adjustable box-wing model (ABW) was used for precise orbit determination (POD). The a priori cuboid box-wing model was directly estimated by the precision radiation accelerations, obtained from ABW POD. When using ECOM1 model, for 14 < β < 40°, a linear systematic variation of D0 related to the elevation of the sun above the orbital plane (β-angle) with a slope of 0.048 nm/s2/°, was found for C30. After adding the cuboid model to assist ECOM1 (named Cuboid + ECOM1), the slope was reduced to 0.005 nm/s2/°, and for C20 satellite, the standard deviation (STD) of D0 was improved, from 1.28 to 0.85 nm/s2 (34%). For satellite laser ranging (SLR) validation, when using the ECOM1 model, the systematic variation with the ε angle was about 14 cm for C20 and C30. After using the Cuboid + ECOM1 model, the variation was significantly reduced to about 5 cm. For C20 and C21, compared with the ECOM1 model, the root mean square (RMS) of the ECOM2 and Cuboid + ECOM1 model was improved by about 0.54 (10.3%) and 0.43 cm (8.7%). For C29 and C30, the RMS of ECOM2 and Cuboid + ECOM1 model was improved for about 0.7 (10.9%) and 1.6 cm (25.6%). Finally, the RMS of the SLR residuals of 4.37 to 4.88 cm was achieved for BeiDou-3 POD.

Graphical Abstract

1. Introduction

The construction of the BeiDou-3 global navigation system began in November 2017. As of April 2019, 18 Medium Earth Orbit (MEO) satellites, one Geostationary Orbit (GEO) satellite, and one Inclined Geosynchronous Orbit (IGSO) satellite were launched into orbit. By the end of 2020, the full operational capability of the global constellation will be completed, and there will be three GEO satellites, three IGSO satellites and 24 MEO satellites [1,2].
Satellite precise orbit determination (POD) is an indispensable task in the global navigation satellite system (GNSS) service. The main factors affecting the accuracy of satellite orbit determination are the distribution of ground tracking stations, data quality of tracking stations, and accuracy of satellite-end models. Satellite attitude models, the phase center model, and the solar pressure model (SRP) are the main models of the satellite-end, and the SRP is the most important. Extending the CODE Orbit model (ECOM1) with five parameters (D0, Y0, B0, Bc, and Bs) produces an empirical SRP model, which is simple and does not require an a priori model. It was widely adopted for the SRP of GPS and GLONASS for the most of Multi-GNSS Experiment (MGEX) analysis centers [3,4,5]. However, for the Galileo and Quasi Zenith Satellite System (QZSS) satellites, the satellite laser ranging (SLR) validation results of the orbits indicate that there is a systematic error associated with the elevation of the sun above the orbital plane (β-angle), when the ECOM1 model is adopted for POD. The main reason is that the shape of the body of the Galileo and QZSS satellites is not similar to the cubic structure of the early GPS satellites, so there is a significant variation of the perturbation acceleration, generated when the sunlight illuminates different surfaces of the satellite-body, and the significant variation cannot be absorbed and compensated by the parameters of the ECOM1 model. Thus, the systematic errors occurred. The compensation methods for the SRP model of the cuboid structure are: (i) Establishing an a priori cuboid box-wing model to enhance the ECOM1 SRP model [6]; and (ii) proposing a new empirical model, named ECOM2, with nine parameters [7].
Prange et al. verified that the ECOM2 model had a significant improvement for Galileo and QZSS in comparison with the ECOM1 model. However, ECOM2 exhibited a slight degradation for BeiDou-2 IGSO and MEO orbits, unlike the ECOM1 model [8]. For the Galileo, Steigenberger et al. compared the ECOM2 and ECOM1, assisted with an a priori cuboid box-wing model, and found that the latter model had a higher accuracy [9]. For the QZS-1 satellite, a semi-analytical SRP modeling, considering the orbit-normal (ON) and yaw steering (YS) modes, was developed in [10]. Zhao et al. established an a priori box-wing model to improve the orbit accuracy of QZS-1, and a significant improvement was achieved when the satellite was in ON mode [11]. For the BeiDou-2 GEO satellites, Wang et al. derived an empirical a priori model to achieve a 10 cm root mean square (RMS) of SLR validation [12]. For the I2-s satellite of BeiDou-3e, Wang et al. analyzed the orbits, based on the ECOM1 and ECOM2 models, and applied an a priori box-wing model and a Fourier-series model to improve the orbit accuracy [13].
The box-wing model (BW) mainly depends on the attitude model as well as accurate information concerning, for example, the satellite dimensions, optical properties, and mass parameters. However, in the absence of this information, the accuracy of the box-wing SRP model was limited. In this case, Rodriguez-Solano et al. developed an adjustable box-wing model (ABW) for GPS satellites, which directly estimates the optical properties, and obtained a comparable result to that using the ECOM1 model [14,15]. For Galileo, the optical and geometrical parameters, released by the European GNSS Agency (GSA), are validated by Duan et al. [16], and the results of the precise orbit determination were calculated by employing the empirical and analytical parameters of Bury et al. [17]. However, the ABW model is more complicated than ECOM1, and the more estimated parameters caused a higher correlation among them. Thus, the optimum results are only achieved with tight a priori constraints on most estimation parameters. Adopting the ABW model, the optical properties of satellites can be estimated from the GNSS data for a long time [16]. Due to the fact that ABW is a semi-analytical model and has physical characteristics, the ABW model is widely used to establish an a priori model for satellites [11,12].
There are two satellite manufacturers for the Beidou-3 satellites [18], the China Academy of Space Technology (CAST) and the Shanghai Engineering Center for Microsatellites (SECM). The main dimensions of SECM satellites are more cuboid than the BeiDou-2 satellites [19,20], and from the POD from the ECOM1 model, we found that the solar radiation acceleration of the BeiDou-3 MEO satellites incurs systematic errors with respect to the β-angle, which was not found for the BeiDou-2 MEO satellite. Therefore, an a priori cuboid box-wing model is necessary to improve the accuracy of POD.
Based on this background, through precise orbit determination using the ABW model, we obtained the precise solar radiation accelerations. Then, an a priori cuboid box-wing was estimated by the precise accelerations. Finally, the performance of the a priori cuboid box-wing model was validated by the orbit accuracy and clock offset.

2. Precise Orbit Determination for the BeiDou-3 Satellites Using the ECOM1 Model

The accuracy of precision orbit determination mainly depends on the accuracy of the SRP model. At the same time, the SRP model needs to be used together with the attitude model. This work mainly analyzes the orbit determination of the satellite for YS and continuous yaw-steering (CYS) modes [1,21].

2.1. ECOM1 SRP Model

Solar radiation pressure, which is the main source of non-gravitational perturbation of the GNSS satellites [22], and has significant influence on precision orbit determination, depends on the optical properties of the reflecting surface and its orientation with respect to the incident radiation.
In the absence of optical properties and dimensions, as an empirical SRP model, the ECOM1 is widely used among most International GNSS Service (IGS) [23] analysis centers for the processing of GPS and GLONASS, as well as the BeiDou-2 IGSO and MEO satellites [24]. This model introduced the right-handed DYB frame [3].
a D = D 0 a Y = Y 0 a B = B 0 + B c cos u + B s sin u ,
where a D , a Y , a B are accelerations at 1AU in the DYB frame. This model consists of three constant accelerations ( D 0 , Y 0 , B 0 ) and two items B c and B s , and u is the argument angle of the latitude of the satellite.

2.2. ECOM2 SRP Model

Arnold et al. implemented the ECOM2 model [7], to which additional estimation parameters are added in the model of Beutler et al. [3]. Second- and fourth-order harmonic terms are added for the D-component.
a D = D 0 + D c 2 cos 2 Δ u + D s 2 sin 2 Δ u + D c 4 cos 4 Δ u + D s 4 sin 4 Δ u a Y = Y 0 a B = B 0 + B c cos Δ u + B s sin Δ u .
Instead of the argument of latitude u , the angular argument Δ u = u u s is used, where u s stands for the argument of latitude of the sun projected on the satellite’s orbital plane.

2.3. Precise Orbit Determination Strategies

Other models involved in precision orbit determination, such as observation models, error models and parameter estimation models are described in Table 1.
The distribution map of the stations of POD is shown in Figure 1. Among them, 56 MGEX/IGS stations [5] and 17 international GNSS Monitoring and Assessment System (iGMAS) stations provide BeiDou-2, BeiDou-3, and GPS observations [18].
We selected one day (4 April 2019) to analyze the number of observations of BeiDou satellites. From Figure 2, it can be seen that: For the BeiDou-3 MEO, the average number of observations of C20, C21, and C22 was 10487, which had a factor that was of 1.5 times larger than that of other BeiDou-3 MEO satellites. For the BeiDou-2 MEO, the average number of observations was about 11758. In general, the number observations of BeiDou-3 MEO were smaller than that of the BeiDou-2 MEO satellite, especially for the C35 satellite.
When the POD is carried out, the phase center offsets (PCO) of the BeiDou-3 MEO satellites are of crucial importance. However, the accurate PCOs of the BeiDou-3 satellites are still un-released. Therefore, the coarse parameters, recommended by iGMAS in 2018, were used and are shown in Table 2.

2.4. D-Acceleration Analysis

In yaw-steering mode, the solar panel of the satellite is always perpendicular to the sun. The solar radiation acceleration generated by the solar panel is constant regardless of the multiplication distance-factor of the satellite to the sun. Thus, the D-component variation can reflect the solar radiation acceleration variation generated by the satellite body with the β (beta) angle. Therefore, the difference of SRP between the BeiDou-2 and BeiDou-3 satellites can be reflected on the D0 (estimation of D-component of ECOM) variations, where C09 and C12 were selected as representatives for the BeiDou-2 satellite, and C20 and C30 were selected for the BeiDou-3 satellites from the CAST and SECM manufacturers, respectively. The D0 variations are shown in Figure 3.
From Figure 3, it can be seen that the BeiDou-2 satellites had a smaller D0 variation than the BeiDou-3 satellites, except for the larger scatters near the maximum absolute β angle and near the zero, and the range of the D0 variations was smaller than 0.5 nm/s2 for C09 and C12. However, the systematic variation of D0 was about 5 nm/s2 when |β| > 40° for C20.
When |β| < 20°, the variation of D0 was about 1.5 nm/s2 for C30; for 20 < β < 40°, there was a linear systematic of D0, with a range of 1.5 nm/s2. Combined with the function of the D0 variations related to the β-angle, derived by Montenbruck et al. [10] from the Box-wing model, there are four parameters affecting the systematic variations of D0, which are a C α δ , a C ρ , a S α δ , and a S ρ (Section 3.2). Only the partial derivatives of a S α δ and a S ρ were a near linear-slope with respect to the β-angle. Therefore, it can be inferred that the C30 satellite has significant a S α δ and a S ρ parameters, resulting from the cuboid of the satellite bus, like the Galileo satellites.
A similar phenomenon was also found for the other satellites of CAST and SECM. Therefore, an a priori box-wing model is necessary to decrease the D0 variation of the BeiDou-3 satellites.

3. A Priori SRP for the BeiDou-3 Satellites

The box-wing model is an analytical model, with physical properties that reflect the actual solar radiation pressure. Therefore, the box-wing model can be an a priori model to compensate for the inadequacy of the empirical model, like ECOM.

3.1. SRP Acceleration Obtained by the Semi-Analytical Model

The box-wing model needs a priori mass, dimensions, and optical properties. For an individual surface, the acceleration vector of the inertial system can be expressed as [22]:
a = A M S c cos θ [ ( α + δ ) e D + 2 ( δ 3 + ρ cos θ ) e N ] ,
where α , ρ , and δ are the physical optical properties of a satellite surface, named the fraction of the absorbed, specular, and reflected photons; A is the surface area; M is the total mass of the satellite; S = 1367 W/m2 denotes the solar irradiance at the Sun–Earth distance; and c is the vacuum velocity of light. e D is the vector from the satellite to the sun, and e N is the normal vector of the illuminated surface. θ is the angle between e D and e N . The satellite body surfaces are covered by multi-layer insulation (MLI), the energy absorbed by the satellite surface is instantaneously re-radiated back into the space, according to Lambert’s law, and the acceleration vector can be expressed as:
a = A M S c cos θ [ ( α + δ ) ( e D + 2 3 e N ) + 2 ρ cos θ e N ] .
For both Equations (3) and (4), a condition of α + δ + ρ = 1 is implicit.
The performance of the box-wing model mainly depends on the accuracy of the surface areas and their optical properties, which is presently only available for GPS, Galileo and, partly, GLONASS. These problems limit the accuracy of the box-wing model. Rodriguez-Solano et al. used the box-wing model to improve the performance of GPS POD by adjusting the optical properties using the ground tracking data. The estimated parameters are absorption plus diffuse reflection (AD) as well as the specular reflection (R) for the illuminated surfaces (+X, +Z, −Z), the optical properties of the solar panel (SP), the constant acceleration (Y0), and the solar panel rotation lag (SB). Thus, this model was named the adjustable box-wing (ABW) [14]. Since the ABW model directly estimates the optical properties, the accuracy of the a priori optical properties is not very critical, unlike for the box-wing model. The metadata of the BeiDou satellites have so far not yet been published by an official organizer. Therefore, the ABW model is an alternative to replace the box-wing model for the BeiDou-3 satellite POD. To achieve the precise solar radiation accelerations of the box-wing, firstly, the precise accelerations of satellite were integrated from the illuminated panels and body-surfaces with the adjustable optical properties of ABW model. The precision of the accelerations was indirectly assessed by the accuracy of POD.
The strategies of data processing are similar to those presented in Table 1, and only the ABW model, with the nine parameters of +XAD, +XR, +ZAD, +ZR, −ZAD, −ZR, SP, Y0, and SB, was used. Considering the correlations between the nine parameters, the constraint of 0.1 was added to +XAD, +XR, +ZAD, +ZR, −ZAD, and −ZR, and a constraint of 1000 was added to SP, Y0, and SB. For the SECM’s satellites, the areas of the satellites were published in the ION GNSS+ conference in 2018 [19] (Table 3), and only the optical properties were used for BeiDou-2 values as an approximation, which are shown in Table 3.
For the CAST’s satellites, the initial optical properties and areas were used for the BeiDou-2 values instead (Table 4, [20]).

3.2. Estimation of the A Priori Cuboid SRP Model

The box-wing model can be reformulated as a function of the sun–satellite–earth angle ( ε ), where the attitude obeys the nominal yaw-steering law [25]. The contributions of the satellite bus to the acceleration can be isolated as the cube, cuboidness, and z-asymmetry, which was described as the parameters of a C α δ , a S α δ , a A α δ , a C ρ , a S ρ , a A ρ , and superscript C, S, and A represent the cubic, stretched body shape, and asymmetric [6,10]. The detailed expression can be written as:
a C α δ = 1 2 ( a Z α δ + a + X α δ ) , a S α δ = 1 2 ( a Z α δ a + X α δ ) a C ρ = 1 2 ( a Z ρ + a + X ρ ) , a S ρ = 1 2 ( a Z ρ a + X ρ ) a A α δ = 1 2 ( a + Z α δ + a Z α δ ) , a A ρ = 1 2 ( a + Z ρ a Z ρ ) a Z α δ = 1 2 ( a + Z α δ + a Z α δ ) , a Z ρ = 1 2 ( a + Z ρ + a Z ρ ) ,
where a i ρ = S A i M c ρ i , a i α δ = S A i M c ( α i + δ i ) , and i = + X , + Z , Z . The acceleration of box-wing model is converted from the inertial system (Equation (4)) to the DYB frame (ECOM1 model). The box-wing acceleration can be expressed as:
a b o x , D = a C α δ ( | cos ε | + sin ε + 2 3 ) a S α δ ( | cos ε | sin ε 4 3 sin 2 ε + 2 3 ) a A α δ ( cos ε + 2 3 | cos ε | cos ε ) 2 a C ρ ( | cos ε | cos 2 ε + sin 3 ε ) 2 a S ρ ( | cos ε | cos 2 ε sin 3 ε ) 2 a A ρ ( cos 3 ε ) a s p ,
a b o x , B = 4 3 a S α δ ( cos ε sin ε ) 2 3 a A α δ ( | cos ε | sin ε ) 2 a C ρ ( | cos ε | sin ε ) cos ε sin ε ) 2 a S ρ ( | cos ε | + sin ε ) cos ε sin ε ) 2 a A ρ ( cos 2 ε sin ε ) ,
where a s p is the main part of the D-acceleration, contributed and engendered by the solar panels, which was nearly a constant for an arc of the POD in the yaw-steering period. So far, the function expressions for the box-wing accelerations, relative to the ε (eps) angle, are all completed.
There are high correlations of SP with +XAD, +ZAD, −ZAD, +XR, +ZR, and −ZR for the ABW model, and the correlation-coefficients were about −0.56, −0.35, −0.41, −0.55, −0.21, and −0.28 for C20 with the β ≈ 7.2°, and were about −0.59, −0.38, −0.37, −0.56, −0.23, and −0.23 for C30 with the β ≈ 31.8°. Due to the higher correlations between most of ABW parameters, it is difficult to achieve the accurate optical parameters of the BeiDou-3 satellites and establish an a priori box-wing model according to Equations (5). However, the precise solar radiation accelerations of the box-wing can be achieved by a high-accuracy orbit.
The accelerations of C20 (day of year (DOY) 263, 2018; β ≈ −15.8°), C30 (DOY 288, 2018; β ≈ −17.9°), and C12 (DOY 288, 2018; β ≈ −18.8°), which were converted into the DYB frame, are plotted in Figure 4. In the yaw-steering period, the ±Y-surface could not be illuminated. In comparison with the D- and B-accelerations, among the three satellites in Figure 4, when the β angle was about -17°, the ranges of D-variations were about 10, 8, and 4 nm/s2 for C30, C20, and C12, in general. BeiDou-3 was larger by two times, compared with the BeiDou-2 satellite. For the B-variations, the ranges were 10, 10, and 6 nm/s2, and BeiDou-3 was larger by nearly two times, compared with the BeiDou-2 satellite.
Based on the function model of Equations (6) and (7), the parameters of a C α δ , a S α δ , a A α δ , a C ρ , a S ρ , a A ρ , a s p were estimated by a separated least-squares adjustment, where the D- and B-accelerations derived from the ABW model were used as the observations. Therefore, the total parameters are a C α δ , a S α δ , a A α δ , a C ρ , a S ρ , a A ρ (named cuboid parameters) and a s p (named constant parameter). To cover a whole range of the β-angle, long range data processing is necessary. Considering the efficiency of the data processing, we used a step of 10 days to cover the β-angle range. The period is selected from the DOY of 218 in 2018 to the DOY of 058 in 2019. The estimation of the cuboid parameters of C20 and C30 are shown in Figure 5.
For the BeiDou-3 satellites, a clear separation between the scatters of a C α δ and a C ρ , which are the main parameters for the cuboid model, can be seen from Figure 5. For the a S α δ and a S ρ , it can also be seen that the coefficients were significant for the C30 satellites. Nevertheless, the a S α δ and a S ρ were close to zero for C20, and the scatters were still stable. For the a A ρ , the estimations of both C20 and C30 were near to zero. However, for the a A α δ , the scatters were larger than the other cuboid parameters, for both C20 and C30. Considering the poorer accuracy and small magnitude, we ignored a A α δ in the a priori cuboid box-wing model (named the cuboid model). The standard deviation (STD) of the corresponding time series was used to evaluate the precision of the cuboid parameter. For C20, the STD values of the corresponding time series of a C α δ , a S α δ , a A α δ , a C ρ , a S ρ , and a A ρ were 0.08, 0.14, 0.60, 0.11, 0.17, and 0.32 nm/s2. For C30, the STD values of the corresponding time series of a C α δ , a S α δ , a A α δ , a C ρ , a S ρ , and a A ρ were 0.08, 0.14, 0.60, 0.11, 0.17, and 0.32 nm/s2, respectively.
The cuboid estimations of the other BeiDou-3 satellites showed the similar results with C20 and C30, and the final cuboid models of CAST and SECM were obtained by averaging each type of satellite. The final estimations, as well as the STD of the cuboid models, are shown in Table 5.

4. Validations

The validations of the cuboid model were based on the D0 variation, orbit accuracy, and clock offset accuracy of the BeiDou-3 satellites. The D0 variation was used to verify the systematic variation of C30. The orbit accuracy can reflect the performance of the SRP model, and the precision of 24 h orbit overlapping and SLR residuals were selected as the two indicators. Meanwhile, considering that the ECOM2 model can also improve the orbit accuracy of Galileo and QZSS, we also compared the Cuboid + ECOM1 with the ECOM2 model. For the clock offset accuracy, the fitting curve of the clock offset, RMS of the fitting residuals, and STD of the clock offset overlapping were selected as the indicators.

4.1. D0 Variation

The smaller the variation of the estimated D0 parameter with the β-angle, the better the performance of the SRP model, when it is close to the actual of the solar radiation pressure. Therefore, the variation of the D0 parameter with the β-angle can be one of the indicators for evaluating the SRP model. The results are shown in Figure 6.
For the ECOM2 model, the scatters of D0 were larger than those for ECOM1 and Cuboid + ECOM1. The main reason was that the function of the D-accelerations was not a constant of ECOM2 (Equation (2)) for an arc-length. Therefore, only the D0 variations of the ECOM1 and Cuboid + ECOM1 were compared. For C20, the STD of D0 was 1.28 and 0.85 nm/s2 for the ECOM1 and Cuboid + ECOM1 solutions, respectively, and we found that the Cuboid + ECOM1 model can reduce the variation by about 34%. However, for 0 < |β| < 40°, a linear systematic variation, with a slope of 0.028 nm/s2/°, was found. This might be improved by accurate body-dimensions, rather than the “guess” values of the CAST’s satellites. For 14 < β < 40°, a linear systematic variation, with a slope of 0.048 nm/s2/°, was found for C30, when the pure ECOM1 model was adopted for POD. After adding the cuboid model, the linear slope was reduced to 0.005 nm/s2/°. This also proved the previous inference on the D0 variation of C30 (Figure 3), i.e., that the systematic variation of D0 was greatly reduced by the a S α δ and a S ρ parameters of the a priori cuboid model.

4.2. Satellite Orbits

The RMS of the 24 h orbit overlapping, which was the orbit of the middle of the day in the first 3d arc solution compared with those of the first day in the second 3d arc solution, and can also reflect the internal coincidence accuracy of the orbit [29]. A comparison of the RMS of the 24 h orbit overlapping using the three models is shown in Table 6, where pseudo-random noise (PRN) and average (AVG) were used.
The RMS for each of satellite was given for the along-track (A), cross-track (C), and radial (R) components. From Table 6, for the ECOM2 model, the averages of RMS were 8.56, 5.31, and 2.71 cm for the A-, C-, and R-components, respectively. Compared with the ECOM1 model, the degradation was about 1 cm for the R-component. However, the improvement of the C-component was about 1 cm, and the 3-dimension (3D) RMS had a slight improvement of about 0.5 cm. For the Cuboid + ECOM1 model, the average improvements of the ACR-components were 0.99 (11.6%), 1.30 (20.8%), and 0.17 cm (9.5%), and the C-component had a higher improvement than the A- and R-components. In general, an average precision of 7.54, 4.96, and 1.62 cm of the 24 h orbit overlapping can be achieved for BeiDou-3 POD, with a 3D precision of 9.2 cm.
SLR is an independent technique operating in the optical frequency range, which is virtually free from the systematic effects related to ionosphere and troposphere delays, ambiguities, and clocks [7]. Therefore, SLR observations are contaminated by only a few error sources and are widely used to evaluate the accuracy of the GNSS orbit [30]. While the high accuracy of the orbits can be determined by the GNSS technique, some orbit modeling deficiencies remain in the SLR residuals, with the characteristic pattern first noted by Urschl et al. [31]. SLR residuals are the differences between the observed (using SLR) and computed distances (derived from the GNSS technique orbits) at the measurement epochs [8]. The station coordinates of the SLR stations are fixed to the a priori reference frame SLRF2008 [32]. The SLR observations are corrected for the relativistic effects, troposphere delays, and for the offset of the laser retroreflector arrays (LRA), related to the satellites’ centers of mass [7]. For the BeiDou-3 MEO satellites, the LRA offsets are shown in Table 7. Residuals larger than ±0.4 m were removed, as gross error, when data editing.
The one-way SLR residuals are shown in Figure 7. The blue line was the cubic polynomial fitting curve of the SLR residuals. From the fitting curve, it can be seen that: For the ECOM1 model result, a systematic variation between the lower and higher ε -angle was about 14 cm for the C20 and C30 satellites. When the ECOM2 and Cuboid + ECOM1 models were used, the variation was significantly reduced to about 5 cm. Thus, the ECOM2 and Cuboid + ECOM1 models can reduce the systematic variation of the SLR residuals.
Compared with the RMS of the SLR residuals of the ECOM1 model, shown in Table 8, for the CAST’s satellites, the STD was improved by about 0.78 cm (15.4%), and the RMS was slightly improved by about 0.54 cm (10.3%) for the ECOM2 model. For the Cuboid + ECOM1 model, the STD was improved by about 1.0 cm (19.3%), and the RMS was slightly improved by about 0.43 (8.7%). However, the offsets of C20 and C21 were increased for about 0.55 and 1.17 cm for the ECOM2 and Cuboid + ECOM1 models, respectively. The possible reasons were the unmodeled acceleration or unmodeled errors, which might also have been caused by the deviation of the satellite-dimensions of CAST, further leading a bias of the cuboid model. Thus, further work will need to be conducted in this direction.
For the SECM’s satellites, compared with the ECOM1 model, the RMS of the ECOM2 and Cuboid + ECOM1 models was improved by about 0.7 cm (10.9%) and 1.6 cm (25.6%), respectively. The Cuboid + ECOM1 model had a higher accuracy than the ECOM2 model.

4.3. Satellite Clock Offsets

The accuracy of the satellite clock offset is consistent with that of the satellite orbit, and the systematic error of the SRP model can also be reflected in the satellite clock offset. Therefore, the quadratic polynomial fitting was performed on each of the 3d-arc length of the satellite clock offset, and then the residuals were used to analyze the accuracy of the ECOM1, ECOM2, and Cuboid + ECOM1 models. The fitted residuals are shown in Figure 8, and the fitted RMS is shown in Table 9.
Comparing Figure 7 and Figure 8, it can be seen that the clock residuals were consistent with those of the satellite orbit and also had a systematic variation with respect to the ε angle. The ECOM2 and Cuboid + ECOM1 model can reduce the systematic variation of the residuals of C20 and C30. For C20, the fitting curve of the ECOM2 had a smaller intercept and slope than Cuboid + ECOM1. However, the fitting RMS of Cuboid + ECOM1 was slightly smaller than that of ECOM2. For C30, the fitting curve of the Cuboid + ECOM1 had a smaller intercept, slope, and RMS than ECOM2. For the Cuboid + ECOM1 model, the improvements of RMS were about 0.06 (24%) and 0.04 ns (14%) for C30 and C20, respectively. In general, the average RMS of ECOM2 was larger than that of the Cuboid + ECOM1 model. The main reason was that the greater number of parameters of ECOM2 reduces their observability and result in less stable results than those of the ECOM1 model, especially for the small observations of the BeiDou-3 satellites.
Considering that the cuboid model can improve the precision of orbit overlapping, the clock offset overlapping was also analyzed. Considering that there is a smaller number of observations (Figure 2), the ratio of the effective clock offset epochs to the theoretical epochs (864 for three days, with an interval of 300 s) is low, for instance, C24 had a ratio of 71%. Therefore, to increase the number of valid epochs to calculate the STD and mean values of the clock offset overlapping, the overlapping arc was set to 48 h, rather than 24 h. The precision of the 48 h clock overlapping is shown in Table 10. The STD and offset were obtained by a quadratic difference between the satellite clock offset series and reference clock offset series. Compared with the ECOM1 model, the average improvements of the STD and offset were about 0.01 ns (4%) and 0.06 ns (11%), for the ECOM2 model, and were 0.02 ns (13%) and 0.06 ns (16%) for the Cuboid + ECOM1 model, respectively.

5. Conclusions

Considering that the metadata of the BeiDou-3 satellite is not fully released at present, we adopted the adjustable box-wing model (ABW) model, for which the optical properties are not very critical, to obtain the precise solar radiation accelerations. The cuboid parameters were directly estimated by the precise solar radiation accelerations. Compared with the method that uses estimated optical properties, together with mass and body-dimensions, to calculate the cuboid parameters, the directly estimated method can avoid the effect of the propagation of satellite mass and dimension errors.
However, this method for estimating the cuboid parameters is not applicable to the orbit-normal mode (ON), because Equations (6) and (7) of box-wing in the DYB-coordinate system are based on the yaw-steering (YS) mode. Therefore, when estimating the cuboid parameters, the results of |β| < 3° were removed. However, when validating the cuboid model, it was found that the orbital accuracy during the period of |β| < 3° is about 4.8 cm of the RMS of the satellite laser ranging (SLR) residuals, which is consistent with the accuracy of the YS mode. Thus, the a priori cuboid model is also applicable for the continuous yaw-steering (CYS) mode. For details concerning the ON mode, please refer to [10,11,12].
The cuboid model can improve the accuracy of precise orbit determination (POD), and the systematic variation of the SLR residuals between the lower and higher sun-satellite-earth angle ( ε ), was reduced from 14 cm to less than 5 cm. The cuboid model can also improve the clock offset accuracy, and the improvements of the RMS of the quadratic polynomial fitting residuals was about 0.06 (24%) and 0.04 ns (14%) for C30 and C20, respectively.
The phase center offsets (PCO) had a high correlation with the ECOM1 model, especially for the D0 parameter, although the BeiDou-2 IGSO and MEO had a smaller D0 variation than BeiDou-3 MEO, and the horizontal PCO scatters still had a systematic variation with respect to the β-angle [33]. A tight D0 constraint can reduce the high correlation between the horizontal PCO and D0 parameters and obtain more stable horizontal PCO estimations, where D0 can be provided by an a priori box-wing model when estimating Galileo PCOs [34]. Thus, the cuboid model might have a better performance in the estimations of PCO for the BeiDou-3 satellites. Currently, the accuracy of POD of the BeiDou-3 satellite is still limited by the number of tracking stations. Therefore, the accuracy of BeiDou-3 POD will be improved by more stations, and more observations, as well as accurate PCOs in the future.

Author Contributions

X.Y. and G.H. provided the initial idea for this study. X.Y. and L.W. designed the experiments. Z.Q. and S.X. performed the experiments. X.Y. and C.L. analyzed the data and wrote the paper. G.H. and Q.Z. gave valuable advice on the paper writing.

Funding

This work was supported by the National Key R and D Program of China (2018YFC1505102), the Programs of the National Natural Science Foundation of China (41774025, 41731066), the Special Fund for Technological Innovation Guidance of Shaanxi Province (2018XNCGG05), the Special Fund for Basic Scientific Research of Central Colleges (grant no. CHD300102269305, CHD300102268305, Chang’an University), and the Grand Projects of the Beidou-2 System (GFZX0301040308).

Acknowledgments

The IGS MGEX, iGMAS, and ILRS are greatly acknowledged for providing the Multi-GNSS and SLR tracking data. Finally, the authors are also grateful for the comments and remarks of the reviewers, who helped to improve the manuscript significantly.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhao, Q.; Wang, C.; Guo, J.; Wang, B.; Liu, J. Precise orbit and clock determination for BeiDou-3 experimental satellites with yaw attitude analysis. GPS Solut. 2018, 22, 4. [Google Scholar] [CrossRef]
  2. Yang, Y.; Xu, Y.; Li, J.; Yang, C. Progress and performance evaluation of BeiDou global navigation satellite system: Data analysis based on BDS-3 demonstration system. Sci. China Earth Sci. 2018, 61, 614–624. [Google Scholar] [CrossRef]
  3. Beutler, G.; Brockmann, E.; Gurtner, W.; Hugentobler, U.; Mervart, L.; Rothacher, M.; Verdun, A. Extended Orbit Modeling Techniques at the CODE Processing Center of the International GPS Service for Geodynamics (IGS): Theory and Initial Results. Manuscripta Geodaetica 1994, 19, 367–386. [Google Scholar]
  4. Deng, Z.; Fritsche, M.; Uhlemann, M.; Wickert, J.; Schuh, H. Reprocessing of GFZ Multi-GNSS Product GBM. In Proceedings of the IGS Workshop 2016, Sydney, Australia, 8–12 February 2016. [Google Scholar]
  5. Montenbruck, O.; Steigenberger, P.; Prange, L.; Deng, Z.; Zhao, Q.; Perosanz, F.; Romero, I.; Noll, C.; Stürze, A.; Weber, G.; et al. The Multi-GNSS experiment (MGEX) of the International GNSS Service (IGS)–achievements, prospects and challenges. Adv. Space Res. 2017, 59, 1671–1697. [Google Scholar] [CrossRef]
  6. Montenbruck, O.; Steigenberger, P.; Hugentobler, U. Enhanced solar radiation pressure modeling for Galileo satellites. J. Geod. 2015, 89, 283–297. [Google Scholar] [CrossRef]
  7. Arnold, D.; Meindl, M.; Beutler, G.; Dach, R.; Schaer, S.; Lutz, S.; Prange, L.; Sośnica, K.; Mervart, L.; Jäggi, A. CODE’s new solar radiation pressure model for GNSS orbit determination. J. Geod. 2015, 89, 775–791. [Google Scholar] [CrossRef]
  8. Prange, L.; Orliac, E.; Dach, R.; Arnold, D.; Beutler, G.; Schaer, S.; Jäggi, A. CODE’s five-system orbit and clock solution the challenges of multi-GNSS data analysis. J. Geod. 2017, 91, 345–360. [Google Scholar] [CrossRef]
  9. Steigenberger, P.; Montenbruck, O. Galileo status: Orbits, clocks, and positioning. GPS Solut. 2016, 21, 319–331. [Google Scholar] [CrossRef]
  10. Montenbruck, O.; Steigenberger, P.; Darugna, F. Semi-analytical solar radiation pressure modeling for QZS-1 orbit-normal and yaw-steering attitude. Adv. Space Res. 2017, 59, 2088–2100. [Google Scholar] [CrossRef]
  11. Zhao, Q.; Chen, G.; Guo, J.; Liu, J.; Liu, X. An a priori solar radiation pressure model for the QZSS Michibiki satellite. J. Geod. 2018, 92, 109–121. [Google Scholar] [CrossRef]
  12. Wang, C.; Guo, J.; Zhao, Q.; Liu, J. Empirically derived model of solar radiation pressure for BeiDou GEO satellites. J. Geod. 2019, 93, 791. [Google Scholar] [CrossRef]
  13. Wang, C.; Guo, J.; Zhao, Q.; Liu, J. Solar Radiation Pressure Models for BeiDou-3 I2-S Satellite: Comparison and Augmentation. Remote Sens. 2018, 10, 118. [Google Scholar] [CrossRef]
  14. Rodriguez-Solano, C.J.; Hugentobler, U.; Steigenberger, P. Adjustable box-wing model for solar radiation pressure impacting GPS satellites. Adv. Space Res. 2012, 49, 1113–1128. [Google Scholar] [CrossRef]
  15. Rodriguez-Solano, C.J.; Hugentobler, U.; Steigenberger, P.; Allende-Alba, G. Improving the orbits of GPS block IIA satellites during eclipse seasons. Adv. Space Res. 2013, 52, 1511–1529. [Google Scholar] [CrossRef]
  16. Duan, B.; Hugentobler, U.; Selmke, I. The Adjusted Optical Properties for Galileo/BeiDou-2/QZS-1 Satellites and Initial Results on BeiDou-3e and QZS-2 satellites. Adv. Space Res. 2019, 63, 1803–1812. [Google Scholar] [CrossRef]
  17. Bury, G.; Zajdel, R.; Sośnica, K. Accounting for perturbing forces acting on Galileo using a box-wing model. GPS Solut. 2019, 23, 74. [Google Scholar] [CrossRef] [Green Version]
  18. Yan, X.; Huang, G.; Zhang, Q.; Liu, C.; Wang, L.; Qin, Z. Early analysis of precise orbit and clock offset determination for the satellites of the global BeiDou-3 system. Adv. Space Res. 2019, 63, 1270–1279. [Google Scholar] [CrossRef]
  19. Xia, L.; Baojun, L.; Yingchun, L.; Sujie, X.; Tao, B. Satellite Geometry and Attitude Mode of MEO Satellites of BDS-3 Developed by SECM; ION GNSS+: Miami, FL, USA, 2018. [Google Scholar]
  20. Guo, J.; Chen, G.; Zhao, Q.; Liu, J.; Liu, X. Comparison of solar radiation pressure models for BDS IGSO and MEO satellites with emphasis on improving orbit quality. GPS Solut. 2017, 21, 511–522. [Google Scholar] [CrossRef]
  21. Dilssner, F. A Note on the Yaw Attitude Modeling of BeiDou IGSO-6, a Report Dated November 20, 2017. Available online: http://navigation-office.esa.int/attachments_24576369_1_BeiDou_IGSO-6_Yaw_Modeling.pdf (accessed on 20 February 2019).
  22. Milani, A.; Nobili, A.M.; Farinella, P. Non-Gravitational Perturbations and Satellite Geodesy; Adam Hilger: Bristol, UK, 1987. [Google Scholar]
  23. Dow, J.M.; Neilan, R.E.; Rizos, C. The international GNSS service in a changing landscape of global navigation satellite systems. J. Geod. 2009, 83, 689. [Google Scholar] [CrossRef]
  24. Guo, J.; Xu, X.; Zhao, Q.; Liu, J. Precise orbit determination for Quad-constellation satellites at Wuhan University: Strategy, result validation, and comparison. J. Geod. 2016, 90, 143–159. [Google Scholar] [CrossRef]
  25. Montenbruck, O.; Schmid, R.; Mercier, F.; Steigenberger, P.; Noll, C.; Fatkulin, R.; Kogure, S.; Ganeshan, A.S. GNSS satellite geometry and attitude models. Adv. Space Res. 2015, 56, 1015–1029. [Google Scholar] [CrossRef] [Green Version]
  26. Dai, X.; Ge, M.; Lou, Y.; Shi, C.; Wickert, J.; Schuh, H. Estimating the yaw-attitude of BDS IGSO and MEO satellites. J. Geod. 2015, 89, 1005–1018. [Google Scholar] [CrossRef]
  27. Saastamoinen, J. Contributions to the theory of atmospheric refraction, part II: Refraction corrections in satellite Geodesy. Bull. Geod. 1973, 105, 279–298. [Google Scholar] [CrossRef]
  28. Boehm, J.; Niell, A.; Tregoning, P.; Schuh, H. Global Mapping Function (GMF): A New Empirical Mapping Function Based on Numerical Weather Model Data. Geophy. Res. Lett. 2006, 33, L07304. [Google Scholar] [CrossRef]
  29. Zhang, R.; Tu, R.; Liu, J.; Hong, J.; Fan, L.; Zhang, P.; Lu, X. Impact of BDS-3 experimental satellites to BDS-2: Service area, precise products, precise positioning. Adv. Space Res. 2018, 62, 829–844. [Google Scholar] [CrossRef]
  30. Sośnica, K.; Thaller, D.; Dach, R.; Steigenberger, P.; Beutler, G.; Arnold, D.; Jäggi, A. Satellite laser ranging to GPS and GLONASS. J. Geod. 2015, 89, 725–743. [Google Scholar] [CrossRef] [Green Version]
  31. Urschl, C.; Beutler, G.; Gurtner, W.; Hugentobler, U.; Ploner, M. Orbit determination for GIOVE-A using SLR tracking data. In Extending the Range, Proceedings of the 15th International Workshop on Laser Ranging, Canberra, Australia, 15–20 October 2006; Available online: http://cddis.gsfc.nasa.gov/lw15/docs/papers/Orbit%20Determination%20for%20GIOVE-A%20using%20SLR%20Tracking%20Data.pdf (accessed on 2 July 2019).
  32. Pavlis, E. SLRF2008: The ILRS reference frame for SLR POD contributed to ITRF2008. In Proceedings of the 2009 Ocean Surface Topography Science Team Meeting, Seattle, WA, USA, 22–24 June 2009. [Google Scholar]
  33. Huang, G.; Yan, X.; Zhang, Q.; Liu, C.; Wang, L.; Qin, Z. Estimation of antenna phase center offset for BDS IGSO and MEO satellites. GPS Solut. 2018, 22, 49. [Google Scholar] [CrossRef]
  34. Steigenberger, P.; Fritsche, M.; Dach, R.; Schmid, R.; Montenbruck, O.; Uhlemann, M.; Prange, L. Estimation of satellite antenna phase center offsets for galileo. J. Geod. 2016, 90, 773–785. [Google Scholar] [CrossRef]
Figure 1. Distributions of precise orbit determination (POD). Magenta points are international GNSS Monitoring and Assessment System (iGMAS) stations, which can track all BeiDou-2/3 satellites. The others are Multi-GNSS Experiment (MGEX) stations, and the orange point can track BeiDou-3 satellites.
Figure 1. Distributions of precise orbit determination (POD). Magenta points are international GNSS Monitoring and Assessment System (iGMAS) stations, which can track all BeiDou-2/3 satellites. The others are Multi-GNSS Experiment (MGEX) stations, and the orange point can track BeiDou-3 satellites.
Remotesensing 11 01605 g001
Figure 2. The number of observations of BeiDou-2 and BeiDou-3 satellites (4 April 2019). Brown bars are numbers for BeiDou-2 Inclined Geosynchronous Orbit (IGSO) and Medium Earth Orbit (MEO) satellites. Green bars are for BeiDou-3 MEO satellites.
Figure 2. The number of observations of BeiDou-2 and BeiDou-3 satellites (4 April 2019). Brown bars are numbers for BeiDou-2 Inclined Geosynchronous Orbit (IGSO) and Medium Earth Orbit (MEO) satellites. Green bars are for BeiDou-3 MEO satellites.
Remotesensing 11 01605 g002
Figure 3. Comparison of the D0 variation between the BeiDou-2 and BeiDou-3 satellites. (ad) are the D0 variation of C09, C12, C20, and C30, respectively.
Figure 3. Comparison of the D0 variation between the BeiDou-2 and BeiDou-3 satellites. (ad) are the D0 variation of C09, C12, C20, and C30, respectively.
Remotesensing 11 01605 g003
Figure 4. D—(green scatters) and B—(red scatters) solar radiation accelerations derived from the adjustable box-wing (ABW) model. (ac) are the D—accelerations for C30, C20, and C12, respectively; and (df) are the B—accelerations for C30, C20, and C12, respectively.
Figure 4. D—(green scatters) and B—(red scatters) solar radiation accelerations derived from the adjustable box-wing (ABW) model. (ac) are the D—accelerations for C30, C20, and C12, respectively; and (df) are the B—accelerations for C30, C20, and C12, respectively.
Remotesensing 11 01605 g004
Figure 5. Estimations of the cuboid parameters. The black and green squares are for a C α δ and a C ρ , the red and purple dots are for a S α δ and a S ρ , and the blue and brown stars are a A α δ and a A ρ . The left figures are for the estimations for C20 (a), and the right ones are for C30 (b).
Figure 5. Estimations of the cuboid parameters. The black and green squares are for a C α δ and a C ρ , the red and purple dots are for a S α δ and a S ρ , and the blue and brown stars are a A α δ and a A ρ . The left figures are for the estimations for C20 (a), and the right ones are for C30 (b).
Remotesensing 11 01605 g005
Figure 6. Estimations of the D0 parameters. The black, red, and blue dots are for the empirical ECOM2, ECOM1, and Cuboid + ECOM1 models, separately. The top figure (a) is for C20, and the bottom figure (b) is for C30.
Figure 6. Estimations of the D0 parameters. The black, red, and blue dots are for the empirical ECOM2, ECOM1, and Cuboid + ECOM1 models, separately. The top figure (a) is for C20, and the bottom figure (b) is for C30.
Remotesensing 11 01605 g006
Figure 7. One-way satellite laser ranging (SLR) residual of the BeiDou-3 satellites. (a,d) are the ECOM1 results of C20 and C30, (b,e) are the ECOM2 results of C20 and C30, and (c,f) are the Cuboid + ECOM1 results of C20 and C30. The blue line is the cubic polynomial fitting curve of the SLR residuals.
Figure 7. One-way satellite laser ranging (SLR) residual of the BeiDou-3 satellites. (a,d) are the ECOM1 results of C20 and C30, (b,e) are the ECOM2 results of C20 and C30, and (c,f) are the Cuboid + ECOM1 results of C20 and C30. The blue line is the cubic polynomial fitting curve of the SLR residuals.
Remotesensing 11 01605 g007
Figure 8. Residuals of the satellite clock offset, fitted by the quadratic polynomial. (a,d) are the ECOM1 results of C30 and C20, (b,e) are the ECOM2 results of C30 and C20, and (c,f) are the Cuboid + ECOM1 results of C30 and C20, respectively. The blue line is the cubic polynomial fitting curve of the residuals.
Figure 8. Residuals of the satellite clock offset, fitted by the quadratic polynomial. (a,d) are the ECOM1 results of C30 and C20, (b,e) are the ECOM2 results of C30 and C20, and (c,f) are the Cuboid + ECOM1 results of C30 and C20, respectively. The blue line is the cubic polynomial fitting curve of the residuals.
Remotesensing 11 01605 g008
Table 1. Strategies for the precise orbit determination (POD).
Table 1. Strategies for the precise orbit determination (POD).
ItemsDescription
StationsAbout 17 iGMAS stations and 56 IGS/MGEX stations [5]
Time periodFrom 181 in 2018 to 125 in 2019
ObservationZero-difference phase and code observations
Elevation-dependent weight; Elevation cutoff angle is 15°
Data arc72-h orbital arcs
Attitude modelBeiDou-3: YS mode [25] and CYS mode [21]
BeiDou-2: YS mode and ON mode [26]
SRP modelECOM1 ( D 0 , Y 0 , B 0 , B c , B s )
Inter-system biases (ISBs)A constant parameter per arc for each station, and zero-mean constraints were added for all ISB parameters
Ionosphere delayIonosphere-free linear combination
GPS: L1/L2
BeiDou: B1I/B3I
Troposphere modelZenith total delay: 2-h intervals;
SAAS+GMF [27,28];
Horizontal gradients: 6-h intervals
Table 2. The coarse parameters of BeiDou-3 satellite phase center offsets (PCOs). The unit is mm.
Table 2. The coarse parameters of BeiDou-3 satellite phase center offsets (PCOs). The unit is mm.
Manufacturerx-offsety-offsetz-offset
CAST−20001500
SECM−500800
Table 3. Initial optical properties and areas of the BeiDou-3 satellites, made by Shanghai Engineering Center for Microsatellites (SECM).
Table 3. Initial optical properties and areas of the BeiDou-3 satellites, made by Shanghai Engineering Center for Microsatellites (SECM).
SurfaceArea (m2) α δ ρ Mass (kg)
Solar panels14.1000.720.000.281010
+ X 1.2420.350.000.65
X 1.2420.350.000.65
Z 2.5763.7480.350.00
Table 4. Initial optical properties and areas of the BeiDou-2 satellites.
Table 4. Initial optical properties and areas of the BeiDou-2 satellites.
SurfaceArea (m2) α δ ρ Mass (kg)
Solar panels22.7040.720.000.28994
+ X 3.7480.350.000.65
X 3.7480.350.000.65
Z 3.4403.7480.350.00
Table 5. Estimations of the cuboid and constant parameters for the BeiDou-3 satellites. The unit is nm/s2.
Table 5. Estimations of the cuboid and constant parameters for the BeiDou-3 satellites. The unit is nm/s2.
TypesCAST (20,21,22,23,24)SECM (25,26,27,29,30)
ParametersValueSTDValueSTD
a C α δ 5.990.113.020.01
a S α δ −0.320.021.050.01
a A α δ −1.180.110.840.08
a C ρ 11.100.215.610.03
a S ρ −0.530.011.950.01
a A ρ 0.210.080.050.02
a sp 110.620.6959.010.81
Table 6. The precision of the root mean square (RMS) of the 24 h orbit overlapping. The unit is cm.
Table 6. The precision of the root mean square (RMS) of the 24 h orbit overlapping. The unit is cm.
PRNECOM1ECOM2Cuboid + ECOM1
ACRACRACR
C206.985.141.305.923.772.566.183.681.56
C217.265.421.426.213.772.156.163.801.54
C226.905.781.385.813.772.445.983.701.50
C238.986.582.1011.315.752.967.325.861.78
C247.846.442.049.835.332.797.125.381.70
C2510.047.682.209.136.561.987.967.002.02
C2612.348.362.4813.528.252.7912.007.842.64
C277.825.941.548.405.332.717.664.761.38
C287.566.041.667.435.573.436.064.561.38
C298.025.801.588.335.232.718.404.481.52
C307.966.061.808.155.563.066.844.441.34
C329.506.001.946.774.962.887.064.461.50
C3310.486.722.049.274.943.319.764.761.84
C347.845.681.608.605.833.237.104.721.04
Avg8.546.261.798.565.312.747.544.961.62
Table 7. The laser retroreflector arrays (LRA) offsets of the BeiDou-3 satellites. The unit is mm.
Table 7. The laser retroreflector arrays (LRA) offsets of the BeiDou-3 satellites. The unit is mm.
Manufacturersx-offsety-offsetz-offset
SECM664.6424.9642.7
CAST612.0−72.01223.0
Table 8. The one-way SLR residuals of the BeiDou-3 satellites. The unit is cm.
Table 8. The one-way SLR residuals of the BeiDou-3 satellites. The unit is cm.
PRNECOM1ECOM2Cuboid + ECOM1
MeanSTDRMSMeanSTDRMSMeanSTDRMS
C20−1.384.865.05−2.134.124.64−2.774.024.88
C21−1.425.285.47−1.774.464.80−2.364.164.78
C29−1.605.495.72−1.345.235.39−1.054.244.37
C30−3.065.876.62−2.205.155.60−1.644.524.81
Table 9. RMS of the clock offset residuals, fitted by a quadratic polynomial for a 3-day arc. The unit is ns.
Table 9. RMS of the clock offset residuals, fitted by a quadratic polynomial for a 3-day arc. The unit is ns.
PRNECOM1ECOM2Cuboid+ECOM1
C200.280.27 0.24
C210.310.29 0.28
C220.340.350.32
C230.260.270.27
C240.210.220.21
C250.390.390.39
C260.330.310.31
C270.220.220.18
C290.240.220.19
C300.250.220.19
C320.210.210.19
C330.270.290.24
C340.270.260.24
Avg0.270.270.25
Table 10. The precision of the 48 h clock overlapping. The unit is ns.
Table 10. The precision of the 48 h clock overlapping. The unit is ns.
PRNECOM1ECOM2Cuboid+ECOM1
STDMeanSTDMeanSTDMean
C200.130.340.130.330.110.28
C210.120.360.120.310.110.33
C220.120.400.130.320.110.30
C230.210.350.200.340.190.37
C240.170.410.170.340.140.34
C250.250.340.190.380.210.32
C260.250.320.210.340.230.26
C270.160.400.160.320.140.31
C290.160.370.160.380.140.33
C300.180.390.160.350.140.33
C320.180.360.160.280.150.32
C330.180.420.190.320.190.32
C340.160.340.180.250.130.24
Avg0.170.370.170.330.150.31

Share and Cite

MDPI and ACS Style

Yan, X.; Liu, C.; Huang, G.; Zhang, Q.; Wang, L.; Qin, Z.; Xie, S. A Priori Solar Radiation Pressure Model for BeiDou-3 MEO Satellites. Remote Sens. 2019, 11, 1605. https://doi.org/10.3390/rs11131605

AMA Style

Yan X, Liu C, Huang G, Zhang Q, Wang L, Qin Z, Xie S. A Priori Solar Radiation Pressure Model for BeiDou-3 MEO Satellites. Remote Sensing. 2019; 11(13):1605. https://doi.org/10.3390/rs11131605

Chicago/Turabian Style

Yan, Xingyuan, Chenchen Liu, Guanwen Huang, Qin Zhang, Le Wang, Zhiwei Qin, and Shichao Xie. 2019. "A Priori Solar Radiation Pressure Model for BeiDou-3 MEO Satellites" Remote Sensing 11, no. 13: 1605. https://doi.org/10.3390/rs11131605

APA Style

Yan, X., Liu, C., Huang, G., Zhang, Q., Wang, L., Qin, Z., & Xie, S. (2019). A Priori Solar Radiation Pressure Model for BeiDou-3 MEO Satellites. Remote Sensing, 11(13), 1605. https://doi.org/10.3390/rs11131605

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