The precise orbit determination (POD) for BeiDou satellites is usually limited by the insufficient quantity and poor distribution of ground tracking stations. To cope with this problem, this study used the GPS and BeiDou joint POD method based on Chinese national continuous operating reference stations (CNCORS) and IGS/MGEX stations. The results show that the 3D RMS of the differences of overlapping arcs is better than 22 cm for geostationary orbit (GEO) satellites and better than 10 cm for inclined geosynchronous orbit (IGSO) and medium earth orbit (MEO) satellites. The radial RMS is better than 2 cm for all three types of BeiDou satellites. The results of satellite laser ranging (SLR) residuals show that the RMS of the IGSO and MEO satellites is better than 5 cm, whereas the GEO satellite has a systematic bias. This study investigates the contributions of CNCORS to the POD of BeiDou satellites. The results show that after the incorporation of CNCORS, the precision of overlapping arcs of the GEO, IGSO, and MEO satellites is improved by 15.5%, 57.5%, and 5.3%, respectively. In accordance with the improvement in the precision of overlapping arcs, the accuracy of the IGSO and MEO satellites assessed by the SLR is improved by 30.1% and 4.8%, respectively. The computation results and analysis demonstrate that the inclusion of CNCORS yields the biggest contribution in the improvement of orbit accuracy for IGSO satellites, when compared to GEO satellites, while the orbit improvement for MEO satellites is the lowest due to their global coverage.

High precision real-time orbit of navigation satellites are usually predicted based on batch estimation solutions, which is highly dependent on the accuracy of the dynamic model. However, for the BDS satellites, the accuracy and reliability of the predicted orbit usually decrease due to the inaccurate dynamic model or orbit maneuvers. To improve this situation, the sequential estimation Square Root Information Filtering (SRIF) was applied to determine the real-time BDS orbits. In the filter algorithm, usually a long period is required for the orbit to converge to the final accuracy, due to lake of accurate initial state. This paper focuses on the impact of the initial state with different a priori Standard Deviation (STD) on the BDS orbit convergence performance in both normal and abnormal periods. For the normal period, the Ultra-Rapid (UR) orbit products and the Broadcast Ephemerides (BRDC) used as initial orbits are discussed respectively. For the abnormal period, orbit maneuver is analyzed. Experimental results show that a proper a priori STD of initial state can significantly accelerate the orbit convergence, while a loose a priori STD takes more than 10 h to converge in the radial direction for the BDS GEO/IGSO/MEO satellites. When the UR orbit product is used as the initial orbit, the orbit of the IGSO/MEO satellites can converge to decimeter-level immediately. When the BRDC product is used, the accuracy of meter-level can be obtained for the IGSO/MEO immediately, and converge to decimeter-level in about 6 h. For the period after the orbit maneuver, the real-time orbit accuracy can reach meter-level in about 6 h after the first group of broadcast ephemerides is received.