Adapting empirical solar radiation pressure model for BDS-3 medium Earth orbit satellites

Abstract

For the precise orbit determination (POD) of global navigation satellite systems (GNSS) constellation, it is very difficult to precisely model the solar radiation pressure (SRP) force acting on GNSS satellites. For GPS satellites, the ECOM model developed by the Center for Orbit Determination in Europe has been utilized by most of International GNSS Service (IGS) analysis centers. However, it should be adapted and optimized to the characteristics of satellites of each GNSS system or even individual satellites. It was extended to the ECOM2 model for GLONASS satellites and then for Galileo satellites by employing a box–wing model. Since November 2020, the third generation of the BeiDou satellite system (BDS-3) has been in its full operation and there are about 200 globally distributed IGS ground stations tracking BDS-3 signals, which creates a great potential to evaluate and optimize its SRP modeling. From the POD processing carried out in this study, we found significant fluctuations of up to 20 cm in overlapping orbit differences for satellites over eclipses in the radial direction and of about 20 and 50 cm in the cross and along directions for ECOM2 and ECOM models. Then, based on numerical analyses we demonstrate that the fourth- and sixth-order sine terms in the Sun direction can significantly reduce the overlapping orbit differences of ECOM. Therefore, an adapted SRP model by adding the fourth- and sixth-order sine periodical terms in the Sun direction to the ECOM model is presented. The adapted model is then validated for BDS-3 POD and orbit prediction. Results show that fluctuations in the amplitude of overlapping estimated orbits using ECOM models are reduced from 20 to < 10 cm in the radial-track component and satellite laser ranging residuals are reduced to half by the adapted SRP model. For the predicted BDS-3 satellite orbits, the RMS values over deep eclipses can be improved from about 7, 14 and 26 cm to about 3, 5 and 12 cm, in the radial, cross and along directions, respectively, compared to the ECOM model.