초록
The SRP (Solar Radiation Pressure) model has always been an issue in the dynamic GPS (Global Positioning System) orbit determination. The widely used CODE (Center for Orbit Determination in Europe) model and its variants have nine parameters to estimate the solar radiation pressure from the Sun and to absorb the remaining forces. However, these parameters show a very high correlation with each other and, therefore, only several of them are estimated at most of the IGS (International GNSS Service) analysis centers. In this study, we attempted to numerically verify the correlation between the parameters. For this purpose, a bi-directional, multi-step numerical integrator was developed. The correlation between the SRP parameters was analyzed in terms of post-fit residuals of the orbit. The integrated orbit was fitted to the IGS final orbit as external observations. On top of the parametric analysis of the SRP parameters, we also verified the capabilities of orbit prediction at later time epochs. As a secondary criterion for orbit quality, the positional discontinuity of the daily arcs was also analyzed. The resulting post-fit RMSE (Root-Mean-Squared Error) shows a level of 4.8 mm on average and there is no significant difference between block types. Since the once-per-revolution parameters in the Y-axis are highly correlated with those in the B-axis, the periodic terms in the D- and Y-axis are constrained to zero in order to resolve the correlations. The 6-hr predicted orbit based on the previous day yields about 3 cm or less compared to the IGS final orbit for a week, and reaches up to 6 cm for 24 hours (except for one day). The mean positional discontinuity at the boundary of two 1-day arcs is on the level of 1.4 cm for all non-eclipsing satellites. The developed orbit integrator shows a high performance in statistics of RMSE and positional discontinuity, as well as the separations of the dynamic parameters. In further research, additional verification of the reference frame for the estimated orbit using SLR is necessary to confirm the consistency of the orbit frames.