Solar Radiation Pressure Parameters Research Articles

Impacts of inter-satellite links on the ECOM model performance for BDS-3 MEO satellites

Inter-satellite link (ISL) plays an essential role in current and future Global Navigation Satellite System (GNSS). In this study, we investigate the impact of ISL observations on precise orbit determination for BeiDou-3 Navigation Satellite System (BDS-3) Medium Earth Orbit (MEO) satellites based on different Extended CODE Orbit Models (ECOM). Thanks to the better observation geometry of the Ka-band ISL data compared to the L-band data for BDS-3 MEO satellites, the ISL solution substantially reduces Orbit Boundary Discontinuity (OBD) errors, except for C30, which suffers from unstable Ka-band hardware delay. From the external quality analysis, ISL significantly enhances the reliability of the orbit of MEO satellites manufactured by the China Academy of Space Technology (CAST). The standard deviation (STD) of the satellite laser ranging (SLR) residuals is approximately 2.5 cm, and the root mean square (RMS) is reduced by 10–23% compared to L-band solutions. Besides, the Sun-elongation angle dependent systematic error in SLR residuals nearly vanishes based on the reduced 5-parameter ECOM (ECOM1) or extended 7-parameter ECOM (ECOM2) with ISL data. This is because the ISL reduces the correlation between state parameters and solar radiation pressure (SRP) parameters as well as those among SRP parameters, leading to a more accurate estimation of both orbit and SRP perturbations, particularly those along B direction. This confirms that the deficiency of the SRP models for BDS-3 CAST satellites can be compensated by using better observation geometry from ISL data. On the other hand, for the satellite manufactured by Shanghai Engineering Center for Microsatellites (SECM), the ISL allows for a more accurate estimation of the Bc1 parameter in the ECOM1 model. This only reduces linear systematic error, possibly because the impact generated by the satellite bus cannot be entirely absorbed by the B-direction parameters.

A Study of Solar Radiation Pressure Model in BDS-3 Precise Orbit Determination

For the BDS (BeiDou Navigation Satellite System) satellites, the solar radiation pressure is the major non-gravitational perturbation the satellites suffer. Due to the influence of multiple factors, accurate modeling the solar radiation pressure acting on these satellites is difficult to accomplish. The solar radiation pressure perturbation is an essential source of error in the precise orbit determination (POD) and orbit prediction processes of the BDS satellites. As the ECOMC (Empirical CODE Orbit Model 1 and 2 Combined) model adopts the advantageous features of the ECOM1 (Empirical CODE Orbit Model 1) model and the ECOM2 (Empirical CODE Orbit Model 2) model, a greater number of parameters have been introduced into the model, leading to strong correlation between the ECOMC model parameters. In response to the deficiency of the ECOMC model, the authors collected the BDS-3 satellite precise ephemeris data that were provided by Wuhan University (WHU) from January 2019 to April 2022. And by these ephemeris data, the ECOMC model’s 13 parameters were obtained through the dynamic orbit fitting method. By conducting these time sequence diagrams analysis of the ECOMC model parameters, this thesis affords the selection strategies of the solar radiation pressure parameters of the ECOMC model for BDS-3 inclined geosynchronous orbit (IGSO) and medium earth orbit (MEO) satellite. And the rationality of the selection strategies of the solar radiation pressure parameters for the ECOMC model was verified by orbit fitting and orbit prediction experiments. The results show that using the improved ECOMC model gains the best orbit fitting effects for BDS-3 IGSO and MEO satellite, and can also enhance the accuracy of medium and long-term orbit prediction of BDS-3 IGSO and MEO satellite.

Evaluation of strategies for the ultra-rapid orbit prediction of BDS GEO satellites

ABSTRACT The quality of BeiDou Navigation Satellite System (BDS) Geostationary Earth Orbit (GEO) ultra-rapid products is unsatisfactory because GEO satellites are nearly stationary relative to ground stations. To optimize the quality of these ultra-rapid orbit products, we investigated the effects of the fitting arc length, an a priori Solar-Radiation Pressure (SRP) model, and the along-track empirical acceleration on the prediction of BDS GEO satellite orbits. The predicted orbit arcs of 24-h were evaluated through comparisons with the corresponding observed orbit arc and Satellite Laser Ranging (SLR) observations. In both eclipse and non-eclipse seasons, accuracy of the orbit predictions obtained using a 48-h fitting arc length were better than those obtained using 24-h and 72-h fitting arc lengths. Although the overlapping precision of predicted orbits exhibited no obvious improvement when an a priori SRP model was employed, the systematic bias in the SLR residuals was significantly reduced. Specifically, the mean value of SLR residuals decreased from −0.248 m to −0.024 m during non-eclipse seasons and from −0.333 m to −0.041 m during eclipse seasons, respectively. In addition, when an empirical acceleration in the along-track direction was introduced, the three-Dimensional Root-Mean-Square (3D RMS) of overlapping orbits during eclipse seasons decreased from 2.964 to 1.080 m, which is comparable to that during non-eclipse seasons. Furthermore, the Standard Deviation (STD) of SLR residuals decreased from 0.419 to 0.221 m during eclipse seasons. The analysis of SRP estimates shows that the stability of SRP parameters was significantly enhanced after the introduction of along-track empirical acceleration in eclipse seasons. The optimal BDS GEO ultra-rapid orbit prediction products were yielded by using a 48-h fitting arc length, an a priori SRP model and an along-track empirical acceleration.

A Hybrid ECOM Model for Solar Radiation Pressure Effect on GPS Reference Orbit Derived by Orbit Fitting Technique

A hybrid ECOM (Empirical CODE Orbit Model) solar radiation pressure (SRP) model, which is termed ECOMC in this work, is proposed for global navigation satellite system (GNSS) orbit modeling. The ECOMC is mainly parameterized by both ECOM1 and ECOM2 models. The GNSS orbit mainly serves as a reference datum not only for its ranging measurement but also for the so-called precise point positioning (PPP) technique. Compared to a complex procedure of orbit determination with real tracking data, the so-called orbit fitting technique simply uses satellite positions from GNSS ephemeris as pseudo-observations to estimate the initial state vector and SRP parameters. The accuracy of the reference orbit is mainly dominated by the SRP, which is usually handled by either ECOM1 or ECOM2. However, the reference orbit derived by ECOM1 produces periodic variations on orbit differences with respect to International GNSS Service (IGS) final orbit for GPS IIR satellites. Such periodic variations are removed from a reference orbit formed using the ECOM2 model, which, however, yields large cross-track orbit errors for the IIR and IIF satellites. Such large errors are attributed to the fact that the ECOM2 intrinsically lacks 1 cycle per revolution (CPR) terms, which stabilize the estimations of the even-order CPR terms in the satellite-Sun direction when the orbit fitting is used. In comparison, a reference orbit constructed with the ECOMC model is free of both the periodic variations from the ECOM1 and the large cross-track orbit errors from the ECOM2. The above improvements from the ECOMC are associated with (1) the even CPR terms removing the periodic variations and (2) the 1 CPR terms compensating for the force mismodeling at ∆u = 90° and 270°, where the ∆u is the argument of the latitude of the satellite with respect to the Sun. The parameter correlation analysis also presents that the direct SRP estimation is sensitive to the 1 and 2 CPR terms in the ECOMC case. In addition, the root-mean-square (RMS) of orbit difference with respect to IGS orbit is improved by ~40%, ~10%, and ~50% in the radial, along-track, and cross-track directions, respectively, when the SRP model is changed from the ECOM2 to the ECOMC. The orbit accuracy is assessed through orbit overlaps at day boundaries. The accuracy improvements of the ECOMC-derived orbit over the ECOM2-derived orbit in the radial, along-track, and cross-track directions are 13.2%, 14.8%, and 42.6% for the IIF satellites and 7.4%, 7.7%, and 35.0% for the IIR satellites. The impact of the reference orbit using the three models on the PPP is assessed. The positioning accuracy derived from the ECOMC is better than that derived from the ECOM1 and ECOM2 by approximately 13% and 20%, respectively. This work may serve as a reference for forming the GNSS reference orbit using the orbit fitting technique with the ECOMC SRP model.

A long‐term broadcast ephemeris model for extended operation of GNSS satellites

GNSS positioning relies on orbit and clock information, which is predicted on ground and transmitted by the individual satellites as part of their broadcast navigation message. The predictions are typically refreshed at least once per day and partitioned into short-arc ephemeris data sets covering a representative validity period of 0.5-4 h. For an increased autonomy of either the space or user segment, the capability to predict a GNSS satellite orbit over extended periods of up to two weeks is studied. A tailored force model for numerical orbit propagation is proposed that offers high accuracy but can still be used in real-time environments. Aside from Earth-orientation parameters, six state vector components and three empirical solar radiation pressure parameters are employed for each satellite and adjusted to past orbits. Using the Galileo constellation with its high-grade hydrogen maser clocks as an example, global average signal-in-space range errors of less than 25 m RMS and 3D position errors of less than about 50 m are demonstrated after two week predictions in 95% of all test cases over a half-year period. The autonomous orbit prediction model thus enables adequate quality for a rapid first fix or contingency navigation in case of lacking ground segment updates.

An empirical solar radiation pressure model for satellites moving in the orbit-normal mode

We present a family of empirical solar radiation pressure (SRP) models suited for satellites orbiting the Earth in the orbit normal (ON) mode. The proposed ECOM-TB model describes the SRP accelerations in the so-called terminator coordinate system. The choice of the coordinate system and the SRP parametrization is based on theoretical assumptions and on simulation results with a QZS-1-like box-wing model, where the SRP accelerations acting on the solar panels and on the box are assessed separately. The new SRP model takes into account that in ON-mode the incident angle of the solar radiation on the solar panels is not constant like in the yaw-steering (YS) attitude mode. It depends on the elevation angle of the Sun above the satellite’s orbital plane. The resulting SRP vector acts, therefore, not only in the Sun-satellite direction, but has also a component normal to it. Both components are changing as a function of the incident angle. ECOM-TB has been used for precise orbit determination (POD) for QZS-1 and BeiDou2 (BDS2) satellites in medium (MEO) and inclined geosynchronous Earth orbits (IGSO) based on IGS MGEX data from 2014 and 2015. The resulting orbits have been validated with SLR, long-arc orbit fits, orbit misclosures, and by the satellite clock corrections based on the orbits. The validation results confirm that—compared to ECOM2—ECOM-TB significantly (factor 3–4) improves the POD of QZS-1 in ON-mode for orbits with different arc lengths (one, three, and five days). Moderate orbit improvements are achieved for BDS2 MEO satellites—especially if ECOM-TB is supported by pseudo-stochastic pulses (the model is then called ECOM-TBP). For BDS2 IGSOs, ECOM-TB with its 9 SRP parameters appears to be over-parameterized. For use with BDS2 IGSO spacecraft we therefore developed a minimized model version called ECOM-TBMP, which is based on the same axis decomposition as ECOM-TB, but has only 2 SRP parameters and is supported by pseudo-stochastic parameters, as well. This model shows a similar performance as ECOM-TB with short arcs, but an improved performance with (3-day) long-arcs. The new SRP models have been activated in CODE’s IGS MGEX solution in Summer 2018. Like the other ECOM models the ECOM-TB derivatives might be used together with an a priori model.

Prediction versus real-time orbit determination for GNSS satellites

To serve real-time users, the IGS (International GNSS Service) provides GPS and GLONASS Ultra-rapid (IGU) orbits with an update of every 6 h. Following similar procedures, we produce Galileo and BeiDou predicted orbits. Comparison with precise orbits from the German Research Centre for Geosciences (GFZ) and Satellite Laser Ranging (SLR) residuals show that the quality of Galileo and BeiDou 6-h predicted orbits decreases more rapidly than that for GPS satellites. Particularly, the performance of BeiDou IGSO and MEO 6-h predicted orbits is 5–6 times worse than the corresponding estimated orbits when satellites are in the eclipse seasons. An insufficient number and distribution of tracking stations, as well as an imperfect solar radiation pressure (SRP) model, limit the quality of Galileo and BeiDou orbit products. Rather than long time prediction, real-time orbit determination by means of a square root information filter (SRIF) produces precise orbits every epoch. By setting variable processing noise on SRP parameters, the filter has the capability of accommodating satellite maneuvers and attitude switches automatically. An epoch-wise ambiguity resolution procedure is introduced to estimate better real-time orbit products. Results show that the real-time estimated orbits are in general better than the 6-h predicted orbits if sufficient observations are available after real-time data preprocessing. On average, 3D RMS values of the real-time estimated orbits reduce by about 30%, 60% and 40% over the 6 h predicted orbits for GPS, BeiDou IGSO and BeiDou MEO eclipsing satellites, respectively. Galileo satellites did not enter into the eclipse season during the experimental period, the standard derivation (STD) of SLR residuals for the real-time estimated orbits are almost the same as for the post-processed orbits.

Comparison of Ultra-Rapid Orbit Prediction Strategies for GPS, GLONASS, Galileo and BeiDou.

Currently, ultra-rapid orbits play an important role in the high-speed development of global navigation satellite system (GNSS) real-time applications. This contribution focuses on the impact of the fitting arc length of observed orbits and solar radiation pressure (SRP) on the orbit prediction performance for GPS, GLONASS, Galileo and BeiDou. One full year’s precise ephemerides during 2015 were used as fitted observed orbits and then as references to be compared with predicted orbits, together with known earth rotation parameters. The full nine-parameter Empirical Center for Orbit Determination in Europe (CODE) Orbit Model (ECOM) and its reduced version were chosen in our study. The arc lengths of observed fitted orbits that showed the smallest weighted root mean squares (WRMSs) and medians of the orbit differences after a Helmert transformation fell between 40 and 45 h for GPS and GLONASS and between 42 and 48 h for Galileo, while the WRMS values and medians become flat after a 42 h arc length for BeiDou. The stability of the Helmert transformation and SRP parameters also confirmed the similar optimal arc lengths. The range around 42–45 h is suggested to be the optimal arc length interval of the fitted observed orbits for the multi-GNSS joint solution of ultra-rapid orbits.

Determining characteristics of artificial near-Earth objects using observability analysis

Observability analysis is a method for determining whether a chosen state of a system can be determined from the output or measurements. Knowledge of state information availability resulting from observability analysis leads to improved sensor tasking for observation of orbital debris and better control of active spacecraft. This research performs numerical observability analysis of artificial near-Earth objects. Analysis of linearization methods and state transition matrices is performed to determine the viability of applying linear observability methods to the nonlinear orbit problem. Furthermore, pre-whitening is implemented to reformulate classical observability analysis. In addition, the state in observability analysis is typically composed of position and velocity; however, including object characteristics beyond position and velocity can be crucial for precise orbit propagation. For example, solar radiation pressure has a significant impact on the orbit of high area-to-mass ratio objects in geosynchronous orbit. Therefore, determining the time required for solar radiation pressure parameters to become observable is important for understanding debris objects. In order to compare observability analysis results with and without measurement noise and an extended state, quantitative measures of observability are investigated and implemented.

A new analytical solar radiation pressure model for current BeiDou satellites: IGGBSPM

An analytical solar radiation pressure (SRP) model, IGGBSPM (an abbreviation for Institute of Geodesy and Geophysics BeiDou Solar Pressure Model), has been developed for three BeiDou satellite types, namely, geostationary orbit (GEO), inclined geosynchronous orbit (IGSO) and medium earth orbit (MEO), based on a ray-tracing method. The performance of IGGBSPM was assessed based on numerical integration, SLR residuals and analyses of empirical SRP parameters (except overlap computations). The numerical results show that the integrated orbit resulting from IGGBSPM differs from the precise ephemerides by approximately 5 m and 2 m for GEO and non-GEO satellites, respectively. Moreover, when IGGBSPM is used as an a priori model to enhance the ECOM (5-parameter) model with stochastic pulses, named ECOM + APR, for precise orbit determination, the SLR RMS residual improves by approximately 20–25 percent over the ECOM-only solution during the yaw-steering period and by approximately 40 percent during the yaw-fixed period. For the BeiDou GEO01 satellite, improvements of 18 and 32 percent can be achieved during the out-of-eclipse season and during the eclipse season, respectively. An investigation of the estimated ECOM D0 parameters indicated that the β-angle dependence that is evident in the ECOM-only solution is no longer present in the ECOM + APR solution.

An impact analysis of arc length on orbit prediction and clock estimation for PPP ambiguity resolution

Real-time precise orbit and clock products are crucial for real-time high-precision precise point positioning (PPP). The key for precise orbit prediction is to precisely determine the initial conditions (satellite position, velocity at the initial epoch) and the solar radiation pressure (SRP) parameters using an arc of known precise orbits. The optimal arc length will be investigated for precise orbit prediction, with which the high-accuracy satellite clocks can be estimated. The orbit prediction accuracies with different arc lengths and conditions over time are also investigated. The results indicate that the predicted orbit accuracy is governed by the SRP parameters. In addition, the predicted orbits and estimated clocks are evaluated quantitatively with respect to PPP solutions during in and out of eclipses seasons, based on which the optimal arc length is identified. The numerical results show that high-accuracy orbital prediction can be achieved using an arc length between 36 and 48 h and generally an arc length of 42 h is considered to be optimal. With this optimal arc length, the estimated clock accuracy is 0.048 and 0.122 ns for a prediction interval of 6 and 24 h, respectively. The average PPP ambiguity fix-rates with a prediction interval of 6 h are 80.2 and 84.7 % over an observation period of 30 min and 1 h, respectively.

Precise orbit determination of BeiDou constellation based on BETS and MGEX network.

Chinese BeiDou Navigation Satellite System is officially operational as a regional constellation with five Geostationary Earth Orbit (GEO) satellites, five Inclined Geosynchronous Satellite Orbit (IGSO) satellites and four Medium Earth Orbit (MEO) satellites. Observations from the BeiDou Experimental Tracking Stations (BETS) and the IGS Multi-GNSS Experiment (MGEX) network from 1 January to 31 March 2013 are processed for orbit determination of the BeiDou constellation. Various arc lengths and solar radiation pressure parameters are investigated. The reduced set of ECOM five-parameter model produces better performance than the full set of ECOM nine-parameter model for BeiDou IGSO and MEO. The orbit overlap for the middle days of 3-day arc solutions is better than 20 cm and 14 cm for IGSO and MEO in RMS, respectively. Satellite laser ranging residuals are better than 10 cm for both IGSO and MEO. For BeiDou GEO, the orbit overlap of several meters and satellite laser ranging residuals of several decimetres can be achieved.

Impact of orbit modeling on DORIS station position and Earth rotation estimates

The high precision of estimated station coordinates and Earth rotation parameters (ERP) obtained from satellite geodetic techniques is based on the precise determination of the satellite orbit. This paper focuses on the analysis of the impact of different orbit parameterizations on the accuracy of station coordinates and the ERPs derived from DORIS observations. In a series of experiments the DORIS data from the complete year 2011 were processed with different orbit model settings. First, the impact of precise modeling of the non-conservative forces on geodetic parameters was compared with results obtained with an empirical-stochastic modeling approach. Second, the temporal spacing of drag scaling parameters was tested. Third, the impact of estimating once-per-revolution harmonic accelerations in cross-track direction was analyzed. And fourth, two different approaches for solar radiation pressure (SRP) handling were compared, namely adjusting SRP scaling parameter or fixing it on pre-defined values.Our analyses confirm that the empirical-stochastic orbit modeling approach, which does not require satellite attitude information and macro models, results for most of the monitored station parameters in comparable accuracy as the dynamical model that employs precise non-conservative force modeling. However, the dynamical orbit model leads to a reduction of the RMS values for the estimated rotation pole coordinates by 17% for x-pole and 12% for y-pole.The experiments show that adjusting atmospheric drag scaling parameters each 30min is appropriate for DORIS solutions. Moreover, it was shown that the adjustment of cross-track once-per-revolution empirical parameter increases the RMS of the estimated Earth rotation pole coordinates. With recent data it was however not possible to confirm the previously known high annual variation in the estimated geocenter z-translation series as well as its mitigation by fixing the SRP parameters on pre-defined values.

Evaluation of GPS orbit prediction strategies for the IGS Ultra-rapid products

The International GNSS Service (IGS) provides Ultra-rapid GPS & GLONASS orbits every 6 h. Each product is composed of 24 h of observed orbits with predicted orbits for the next 24 h. We have studied how the orbit prediction performance varies as a function of the arc length of the fitted observed orbits and the parameterization strategy used to estimate the empirical solar radiation pressure (SRP) effects. To focus on the dynamical aspects of the problem, nearly ideal conditions have been adopted by using IGS Rapid orbits and known earth rotation parameters (ERPs) as observations. Performance was gauged by comparison with Rapid orbits as truth by examining WRMS and median orbit differences over the first 6-h and the full 24-h prediction intervals, as well as the stability of the Helmert frame alignment parameters. Two versions of the extended SRP orbit model developed by the Centre for Orbit Determination in Europe (CODE) were tested. Adjusting all nine SRPs (offsets plus once-per-revolution sines and cosines in each satellite-centered frame direction) for each satellite shows smaller mean sub-daily, scale, and origin translation differences. On the other hand, eliminating the four once-per-revolution SRP parameters in the sun-ward and the solar panel axis directions yields orbit predictions that are much more rotationally stable. We found that observed arc lengths of 40---45 h produce the most stable and accurate predictions during 2010. A combined strategy of rotationally aligning the 9 SRP results to the 5 SRP frame should give optimal predictions with about 13 mm mean WRMS residuals over the first 6 h and 50 mm over 24 h. Actual Ultra-rapid performance will be degraded due to the unavoidable rotational errors from ERP predictions.

Orbit determination and prediction accuracy analysis for a regional tracking network

China’s COMPASS satellite navigation system relies on a regional tracking network to provide navigation services. Limited by its geographic border, the regional network is able to cover only 30% of the medium-earth-orbits (MEO). Accuracy of determined and predicted orbits is not able to satisfy system requirements if the tracking data processing strategy for global tracking network processing is used for the regional network. Two major error sources for orbital prediction are accuracy of initial orbital elements and dynamical modeling. To achieve better prediction accuracy, we propose a two-step orbit determination and prediction strategy. For step 1, only solar radiation pressure (SRP) parameters are estimated along with the orbital elements and other parameters; for step 2, all parameters are estimated but the SRP parameters are tightly constrained to their step 1 estimates. Experimenting with data from a regional GPS network, we conclude for orbital prediction using the proposed two-step strategy, the average user range error (URE) for 24-h prediction arcs is better than 0.6 m.

DORIS processing at the European Space Operations Centre

This paper gives an overview of the DORIS related activities at the Navigation Support Office of the European Space Operations Centre. The DORIS activities were started in 2002 because of the launch of the Envisat satellite where ESOC is responsible for the validation of the Envisat Precise Orbits and a brief overview of the key Envisat activities at ESOC is given. Typical orbit comparison RMS values between the CNES POE (GDR-C) and the ESOC POD solution is 6.5, 18.8 and 23.1 mm in radial-, along- and cross-track direction. In the framework of the generation of the ITRF2008 ESOC participated in the reprocessing of all three space geodetic techniques; DORIS, SLR, and GPS. Here the main results of our DORIS reprocessing, in the framework of the International DORIS Service (IDS), are given. The WRMS of the weekly ESOC solution (esawd03) for the 2004–2009 period compared to the IDS-1 combined solution is of the order of 12 mm. Based on the long time series of homogeneously processed data a closer look is taken at the estimated solar radiation pressure parameters of the different satellites used in this DORIS analysis. The main aim being the stabilization of the Z-component of the geocentre estimates. We conclude that the ESOC participation to the IDS ITRF2008 contribution has been beneficial for both ESOC and the IDS. ESOC has profited significantly from the very open and direct communications and comparisons that took place within the IDS during the reprocessing campaign.

Sensitivity of GPS and GLONASS orbits with respect to resonant geopotential parameters

The Center for Orbit Determination in Europe (CODE) has been involved in the processing of combined GPS/GLONASS data during the International GLONASS Experiment (IGEX). The resulting precise orbits were analyzed using the program SORBDT. Introducing one satellite’s positions as pseudo-observations, the program is capable of fitting orbital arcs through these positions using an orbit improvement procedure based on the numerical integration of the satellite’s orbit and its partial derivative with respect to the orbit parameters. For this study, the program was enhanced to estimate selected parameters of the Earth’s gravity field. The orbital periods of the GPS satellites are —in contrast to those of the GLONASS satellites – 2:1 commensurable (P Sid:P GPS) with the rotation period of the Earth. Therefore, resonance effects of the satellite motion with terms of the geopotential occur and they influence the estimation of these parameters. A sensitivity study of the GPS and GLONASS orbits with respect to the geopotential coefficients reveals that the correlations between different geopotential coefficients and the correlations of geopotential coefficients with other orbit parameters, in particular with solar radiation pressure parameters, are the crucial issues in this context. The estimation of the resonant geopotential terms is, in the case of GPS, hindered by correlations with the simultaneously estimated radiation pressure parameters. In the GLONASS case, arc lengths of several days allow the decorrelation of the two parameter types. The formal errors of the estimates based on the GLONASS orbits are a factor of 5 to 10 smaller for all resonant terms.

A new solar radiation pressure model for GPS

The largest error source in GPS orbits is due to the effects of the solar radiation pressure. Over the last few years many improvements were made in modeling the orbits of GPS satellites within the IGS. However, most improvements were achieved by increasing the number of estimated orbit and/or solar radiation pressure parameters. This increased number of estimated satellite parameters weakens the solutions of all estimated parameters. Due to correlations, the additional parameters may cause biases in other estimated quantities like, e.g., the length of day. In this paper we present a recently developed solar radiation pressure model for the GPS satellites. This model is based on experiences and results acquired at the Center of Orbit Determination in Europe (CODE) in the context of its IGS activities since June 1992. The new solar radiation pressure model outperforms the existing ROCK models by almost an order of magnitude. It also allows a reduction of the number of orbit parameters that have to be estimated.

A New Solar Radiation Pressure Model for GPS Satellites

The largest error in currently used GPS orbit models is due to the effect of solar radiation pressure. Over the last few years many improvements were made in modeling the orbits of GPS satellites within the International GPS Service (IGS). Howeer, most improvements were achieved by increasing the number of estimated orbit and/or solar radiation pressure parameters. This increase in the number of estimated satellite parameters weakens the solutions of all estimated parameters (not only orbit parameters). Because of correlations the additional orbit parameters may introduce biases in other estimated quantities, for example the length of day. We present a recently developed solar radiation pressure model for the GPS satellites. This model is based on experiences and results gained at the Center for Orbit Determination in Europe (CODE) in the context of its IGS activities since June 1992. The performance of the new model is almost an order of magnitude better than that of the existing ROCK models. It also allows a reduction of the number of orbit parameters that have to be estimated. © 1999 John Wiley & Sons, Inc.