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/s²/°, was found for C30. After adding the cuboid model to assist ECOM1 (named Cuboid + ECOM1), the slope was reduced to 0.005 nm/s²/°, and for C20 satellite, the standard deviation (STD) of D0 was improved, from 1.28 to 0.85 nm/s² (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.

Satellite antenna phase center offsets (PCOs) and phase variations (PVs) for BeiDou-3 satellites are estimated based on the tracking data of the Multi-GNSS Experiment (MGEX) and the international GNSS Monitoring and Assessment System (iGMAS) network. However, when estimating the (PCOs) of BeiDou-3 medium Earth orbit (MEO) satellites by pure Extending the CODE Orbit Model (ECOM1), the x-offset estimations of the PCOs have a systematic variation of about 0.4 m with the elevation of the Sun above the orbital plane (β-angle). Thus, a priori box-wing solar radiation pressure (SRP) model of BeiDou-3 MEO was assisted with ECOM1. Then, the satellite type-specific PCOs and common PVs were obtained. The estimations of PCOs and PVs were compared with the MGEX PCOs from the precise orbit and clock offset. When the MGEX PCOs were used, the root mean square (RMS) of 24 h overlap was 6.76, 4.36, 1.46 cm, in along-track, cross-track, and radial directions, respectively; the RMS and standard deviations (STD) of the 24 h clock offset overlap were 0.28 and 0.15 ns; the fitting RMS of the 72 h clock offset of the quadratic polynomial was 0.243 ns. After comparing this with the estimated PCOs and PVs, the RMS of the 24 h orbit overlap was decreased by 6.5 mm (10.54%), 1.8 mm (4.4%), and 1.1 mm (8.03%) in the along-track, cross-track, and radial directions, respectively; the RMS and STD of the 24 h clock offset overlap were decreased by 0.024 ns (8.6%) and 0.020 ns (13.1%), respectively; the fitting RMS of the 72 h clock offset of the quadratic polynomial was reduced by about 0.016 ns (6.5%).

In order to provide better service for the Asia-Pacific region, the BeiDou navigation satellite system (BDS) is designed as a constellation containing medium earth orbit (MEO), geostationary earth orbit (GEO), and inclined geosynchronous orbit (IGSO). However, the multi-orbit configuration brings great challenges for orbit determination. When orbit maneuvering, the orbital elements of the maneuvered satellites from broadcast ephemeris are unusable for several hours, which makes it difficult to estimate the initial orbit in the process of precise orbit determination. In addition, the maneuvered force information is unknown, which brings systematic orbit integral errors. In order to avoid these errors, observation data are removed from the iterative adjustment. For the above reasons, the precise orbit products of maneuvered satellites are missing from IGS (international GNSS (Global Navigation Satellite System) service) and iGMAS (international GNSS monitoring and assessment system). This study proposes a method to determine the precise orbits of maneuvered satellites for BeiDou GEO and IGSO. The initial orbits of maneuvered satellites could be backward forecasted according to the precise orbit products. The systematic errors caused by unmodeled maneuvered force are absorbed by estimated pseudo-stochastic pulses. The proposed method for determining the precise orbits of maneuvered satellites is validated by analyzing data of stations from the Multi-GNSS Experiment (MGEX). The results show that the precise orbits of maneuvered satellites can be estimated correctly when orbit maneuvering, which could supplement the precise products from the analysis centers of IGS and iGMAS. It can significantly improve the integrality and continuity of the precise products and subsequently provide better precise products for users.

Because there are different types of BeiDou constellations with participating geostationary orbit (GEO) and inclined geosynchronous orbit (IGSO) satellites, the maneuvering frequency of BeiDou satellites is higher than that of other navigation systems. The satellite orbital maneuvers lead to orbital parameter failure for several hours from broadcast ephemeris. Due to the missing initial orbit, the maneuvering thrust, and the period of orbital maneuvering, the orbit products of maneuvering satellites cannot be provided by the International Global Navigation Satellite System (GNSS) Service (IGS) and International GNSS Monitoring and Assessment System (iGMAS). In addition, the period of unhealthy status and the orbital parameters of maneuvering satellites in broadcast ephemeris are unreliable, making the detection of orbital maneuver periods more difficult. Here, we develop a method to detect orbital maneuver periods involving two key steps. The first step is orbit prediction of maneuvering satellites based on precise orbit products. The second step is time period detection of orbit maneuvering. The start time detection factor is calculated by backward prediction orbit and pseudo-range observations, and the end time detection factor is calculated by forward prediction orbit and pseudo-range observations. Data of stations from the Multi-GNSS Experiment (MGEX) and iGMAS were analyzed. The results show that the period of orbit maneuvering could be detected accurately for BeiDou GEO and IGSO satellites. In addition, the orbital maneuver period of other GNSS medium Earth orbit (MEO) satellites could also be determined by this method. The results of period detection for orbit maneuvering provide important reference information for precision orbit and clock offset determination during satellite maneuvers.