Dynamic and non-dynamic ERS-1 radial orbit improvement from ERS-1/TOPEX dual-satellite altimetry

Dual-satellite altimeter crossover differences between ERS-1 and TOPEX/Poseidon have been included as supplementary tracking data in ERS-1 orbit computations from SLR and single-satellite crossover differences. It was found that including the dual-satellite crossover differences slightly improves the ERS-1 radial orbit accuracy of about 12 cm for orbits computed with the JGM-2 gravity field and also leads to a better ‘centering’ of the ERS-1 orbit in the terrestrial reference frame defined for TOPEX/Poseidon. In addition to this dynamic orbit improvement technique, a non-dynamic technique has been investigated that removes the larger part of the ERS-1 radial orbit error from the dual-satellite crossover difference residuals. For ERS-1 orbits computed with the GEM-T2 gravity field, it was found that the non-dynamic technique could improve the radial orbit accuracy from 140 cm to the same level of accuracy as the ERS-1 JGM-2 orbits.     

Precision orbit determination for TOPEX/Poseidon using TDRSS doppler tracking data

Precision orbit determination on the TOPEX/Poseidon (T/P) altimeter satellite is now being routinely achieved with sub-5cm radial and sub-15 cm total positioning accuracy using state-of-the-art modeling with precision tracking provided by a combination of: (a) global Satellite Laser Ranging (SLR) and Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS), or (b) the Global Positioning System (GPS) Constellation which provides pseudo-range and carrier phase observations. The geostationary Tracking and Data Relay Satellite System (TDRSS) satellites are providing the operational tracking and communication support for this mission. The TDRSS Doppler data are of high precision (0.3 mm/s nominal noise levels). Unlike other satellite missions supported operationally by TDRSS, T/P has high quality independent tracking which enables absolute orbit accuracy assessments. In addition, the T/P satellite provides extensive geometry for positioning a satellite at geostationary altitude, and thus the TDRSS-T/P data provides an excellent means for determining the TDRS orbits. Arc lengths of 7 and 10 days with varying degrees of T/P spacecraft attitude complexity are studied. Sub-meter T/P total positioning error is achieved when using the TDRSS range-rate data, with radial orbit errors of 10.6 cm and 15.5 cm RMS for the two arcs studied. Current limitations in the TDRSS precision orbit determination capability include mismodeling of numerous TDRSS satellite-specific dynamic and electronic effects, and in the inadequate treatment of the propagation delay and bending arising from the wet troposphere and ionosphere.     

Gravity field model adjustment from ERS-1 and TOPEX altimetry and SLR tracking using an analytical orbit perturbation theory

The JGM-2 gravity field model has been adjusted using 70 days of ERS-1 and 19 10-day repeat cycles of TOPEX/Poseidon SLR and single satellite altimeter crossover differences. In addition, dual satellite altimeter crossover differences between ERS-1 and TOPEX for the selected 70-day period of ERS-1 and TOPEX repeat cycle 18 have been used in the JGM-2 model adjustment. In the computation of the normal equations, use was made of the analytical Lagrange linear perturbation theory. The single satellite altimeter crossover difference rms was brought down from 15.2 to 14.1 cm for ERS-1 and from 10.6 to 10.5 cm for TOPEX. The dual satellite altimeter crossover difference rms was reduced from 17.7 to 16.9 cm. Furthermore, the weighted rms of fit of SLR measurements was brought down from 16.0 to 14.4 cm for ERS-1, and from 5.4 to 5.0 cm for TOPEX/Poseidon.     

Towards development of a consistent orbit series for TOPEX,Jason-1, and Jason-2

The TOPEX/Poseidon, Jason-1 and Jason-2 set of altimeter data now provide a time series of synoptic observations of the ocean that span nearly 17 years from the launch of TOPEX in 1992. The analysis of the altimeter data including the use of altimetry to monitor the global change in mean sea level requires a stable, accurate, and consistent orbit reference over the entire time span. In this paper, we describe the recomputation of a time series of orbits that rely on a consistent set of reference frames and geophysical models. The recomputed orbits adhere to the IERS 2003 standards for ocean and earth tides, use updates to the ITRF2005 reference frame for both the SLR and DORIS stations, apply GRACE-derived models for modeling of the static and time-variable gravity, implement the University College London (UCL) radiation pressure model for Jason-1, use improved troposphere modeling for the DORIS data, and apply the GOT4.7 ocean tide model for both dynamical ocean tide modeling and for ocean loading. The new TOPEX orbits have a mean SLR fit of 1.79 cm compared to 2.21 cm for the MGDR-B orbits. These new TOPEX orbits agree radially with independent SLR/crossover orbits at 0.70 cm RMS, and the orbit accuracy is estimated at 1.5–2.0 cm RMS over the entire TOPEX time series. The recomputed Jason-1 orbits agree radially with the Jason-1 GDR-C orbits at 1.08 cm RMS. The GSFC SLR/DORIS dynamic and reduced-dynamic orbits for Jason-2 agree radially with independent orbits from the CNES and JPL at 0.70–1.06 cm RMS. Applying these new orbits, and using the latest altimeter corrections for TOPEX, Jason-1, and Jason-2 from September 1992 to May 2009, we find a global rate in mean sea level of 3.0 ± 0.4 mm/yr.     

GSFC DORIS contribution to ITRF2008

The NASA GSFC DORIS analysis center has provided weekly DORIS solutions from November 1992 to January 2009 (839 SINEX files) of station positions and Earth Orientation Parameters for inclusion in the DORIS contribution to ITRF2008. The NASA GSFC GEODYN orbit determination software was used to process the orbits and produce the normal equations. The weekly SINEX gscwd10 submissions included DORIS data from Envisat, TOPEX/Poseidon, SPOT-2, SPOT-3, SPOT-4, SPOT-5. The orbits were mostly seven days in length (except for weeks with data gaps or maneuvers). The processing used the GRACE-derived EIGEN-GL04S1 gravity model, updated modeling for time-variable gravity, the GOT4.7 ocean tide model and tuned satellite-specific macromodels for SPOT-2, SPOT-3, SPOT-4, SPOT-5 and TOPEX/Poseidon. The University College London (UCL) radiation pressure model for Envisat improves nonconservative force modeling for this satellite, reducing the median residual empirical daily along-track accelerations from 3.75 × 10⁻⁹ m/s² with the a priori macromodel to 0.99 × 10⁻⁹ m/s² with the UCL model. For the SPOT and Envisat DORIS satellite orbits from 2003 to 2008, we obtain average RMS overlaps of 0.8–0.9 cm in the radial direction, 2.1–3.4 cm cross-track, and 1.7–2.3 cm along-track. The RMS orbit differences between Envisat DORIS-only and SLR & DORIS orbits are 1.1 cm radially, 6.4 cm along-track and 3.7 cm cross-track and are characterized by systematic along-track mean offsets due to the Envisat DORIS system time bias of ±5–10 μs. We obtain a good agreement between the geometrically-determined geocenter parameters and geocenter parameters determined dynamically from analysis of the degree one terms of the geopotential. The intrinsic RMS weekly position repeatability with respect to the IDS-3 combination ranges from 2.5 to 3.0 cm in 1993–1994 to 1.5 cm in 2007–2008.     

Impact of ITRS 2014 realizations on altimeter satellite precise orbit determination

This paper evaluates orbit accuracy and systematic error for altimeter satellite precise orbit determination on TOPEX, Jason-1, Jason-2 and Jason-3 by comparing the use of four SLR/DORIS station complements from the International Terrestrial Reference System (ITRS) 2014 realizations with those based on ITRF2008. The new Terrestrial Reference Frame 2014 (TRF2014) station complements include ITRS realizations from the Institut National de l’Information Géographique et Forestière (IGN) ITRF2014, the Jet Propulsion Laboratory (JPL) JTRF2014, the Deutsche Geodätisches Forschungsinstitut (DGFI) DTRF2014, and the DORIS extension to ITRF2014 for Precise Orbit Determination, DPOD2014. The largest source of error stems from ITRF2008 station position extrapolation past the 2009 solution end time. The TRF2014 SLR/DORIS complement impact on the ITRF2008 orbit is only 1–2 mm RMS radial difference between 1992–2009, and increases after 2009, up to 5 mm RMS radial difference in 2016. Residual analysis shows that station position extrapolation error past the solution span becomes evident even after two years, and will contribute to about 3–4 mm radial orbit error after seven years. Crossover data show the DTRF2014 orbits are the most accurate for the TOPEX and Jason-2 test periods, and the JTRF2014 orbits for the Jason-1 period. However for the 2016 Jason-3 test period only the DPOD2014-based orbits show a strong and statistically significant margin of improvement. The positive results with DTRF2014 suggest the new approach to correct station positions or normal equations for non-tidal loading before combination is beneficial. We did not find any compelling POD advantage in using non-linear over linear station velocity models in our SLR & DORIS orbit tests on the Jason satellites. The JTRF2014 proof-of-concept ITRS realization demonstrates the need for improved SLR+DORIS orbit centering when compared to the Ries (2013) CM annual model. Orbit centering error is seen as an annual radial signal of 0.4 mm amplitude with the CM model. The unmodeled CM signals show roughly a 1.8 mm peak-to-peak annual variation in the orbit radial component. We find the TRF network stability pertinent to POD can be defined only by examination of the orbit-specific tracking network time series. Drift stability between the ITRF2008 and the other TRF2014-based orbits is very high, the relative mean radial drift error over water is no larger than 0.04 mm/year over 1993–2015. Analyses also show TRF induced orbit error meets current altimeter rate accuracy goals for global and regional sea level estimation.     

Derniers resultats de calcul d'orbite SPOT et TOPEX/POSEIDON obtenus avec DORIS

DORIS is a tracking system developed by CNES to support precise orbit computation of the US-French TOPEX/POSEIDON project. Moreover DORIS data are currently processed to compute SPOT2 and SPOT3 orbits. Although the SPOT satellites are at 800 km altitudes, their orbits reach the decimetric level, to be compared to the 3–4 cm RMS on the T/P radial component.For each type of orbit, there is an adaptive period of a few months which is used for improving the precision. The paper describes what has been done to reduce some items in the total error budget. The latest results will be presented as well as the criteria which are settled to characterize the improvements.Finally, the future developments of DORIS and the potential projects flowing DORIS will be outlined such as ENVISAT and TOPEX/POSEIDON Follow On.     

TOPEX/Poseidon precise orbit determination using GPS double difference phase observables

Data from the experimental onboard GPS receiver were used to accurately compute the orbit of TOPEX/Poseidon. This represents a unique opportunity to intercompare with two other classical tracking techniques (SLR and DORIS). A review of the methodology used is given together with current results.     

Generating precise and homogeneous orbits for Jason-1 and Jason-2

Driven by the GMES (Global Monitoring for Environment and Security) and GGOS (Global Geodetic Observing System) initiatives the user community has a strong demand for high-quality altimetry products. In order to derive such high-quality altimetry products, precise orbits for the altimetry satellites are a necessity. With the launch of the TOPEX/Poseidon mission in 1992 a still on-going time series of high-accuracy altimetry measurements of ocean topography started, continued by the altimetry missions Jason-1 in 2001 and Jason-2/OSTM in 2008. This paper contributes to the on-going orbit reprocessing carried out by several groups and presents the efforts of the Navigation Support Office at ESA/ESOC using its NAPEOS software for the generation of precise and homogeneous orbits referring to the same reference frame for the altimetry satellites Jason-1 and Jason-2. Data of all three tracking instruments on-board the satellites (beside the altimeter), i.e. GPS, DORIS, and SLR measurements, were used in a combined data analysis. About 7 years of Jason-1 data and more than 1 year of Jason-2 data were processed. Our processing strategy is close to the GDR-C standards. However, we estimated slightly different scaling factors for the solar radiation pressure model of 0.96 and 0.98 for Jason-1 and Jason-2, respectively. We used 30 s sampled GPS data and introduced 30 s satellite clocks stemming from ESOC’s reprocessing of the combined GPS/GLONASS IGS solution. We present the orbit determination results, focusing on the benefits of adding GPS data to the solution. The fully combined solution was found to give the best orbit results. We reach a post-fit RMS of the GPS phase observation residuals of 6 mm for Jason-1 and 7 mm for Jason-2. The DORIS post-fit residuals clearly benefit from using GPS data in addition, as the DORIS data editing improves. The DORIS observation RMS for the fully combined solution is with 3.5 mm and 3.4 mm, respectively, 0.3 mm better than for the DORIS-SLR solution. Our orbit solution agrees well with external solutions from other analysis centers, as CNES, LCA, and JPL. The orbit differences between our fully combined orbits and the CNES GDR-C orbits are of about 0.8 cm for Jason-1 and at 0.9 cm for Jason-2 in the radial direction. In the cross-track component we observe a clear improvement when adding GPS data to the POD process. The 3D-RMS of the orbit differences reveals a good orbit consistency at 2.7 cm and 2.9 cm for Jason-1 and Jason-2. Our resulting orbit series for both Jason satellites refer to the ITRF2005 reference frame and are provided in sp3 file format on our ftp server.     

Precision orbit determination performance for CryoSat-2

In this paper we discuss our efforts to perform precision orbit determination (POD) of CryoSat-2 which depends on Doppler and satellite laser ranging tracking data. A dynamic orbit model is set-up and the residuals between the model and the tracking data is evaluated. The average r.m.s. of the 10?s averaged Doppler tracking pass residuals is approximately 0.39?mm/s; and the average of the laser tracking pass residuals becomes 1.42?cm. There are a number of other tests to verify the quality of the orbit solution, we compare our computed orbits against three independent external trajectories provided by the CNES. The CNES products are part of the CryoSat-2 products distributed by ESA. The radial differences of our solution relative to the CNES precision orbits shows an average r.m.s. of 1.25?cm between Jun-2010 and Apr-2017. The SIRAL altimeter crossover difference statistics demonstrate that the quality of our orbit solution is comparable to that of the POE solution computed by the CNES. In this paper we will discuss three important changes in our POD activities that have brought the orbit performance to this level. The improvements concern the way we implement temporal gravity accelerations observed by GRACE; the implementation of ITRF2014 coordinates and velocities for the DORIS beacons and the SLR tracking sites. We also discuss an adjustment of the SLR retroreflector position within the satellite reference frame. An unexpected result is that we find a systematic difference between the median of the 10 s Doppler tracking residuals which displays a statistically significant pattern in the South Atlantic Anomaly (SSA) area where the median of the velocity residuals varies in the range of ?0.15 to +0.15?mm/s.     

Evaluation of JGM 2 geopotential errors from Geosat, TOPEX/Poseidon and ERS-1 crossover altimetry

World-ocean distribution of the crossover altimetry data from Geosat, TOPEX/Poseidon (T/P) and the ERS 1 missions have provided strong independent evidence that NASA's/CSR's JGM 2 geopotential model (70 × 70 in spherical harmonics) yields accurate radial ephemerides for these satellites. In testing the sea height crossover differences found from altimetry and JGM 2 orbits for these satellites, we have used the sea height differences themselves (of ascending minus descending passes averaged at each location over many exact repeat cycles) and the Lumped Latitude Coefficients (LLC) derived from them. For Geosat we find the geopotential-induced LLC errors (exclusive of non-gravitational and initial state discrepancies) mostly below 6 cm, for TOPEX the corresponding errors are usually below 2 cm, and for ERS 1 (35-day cycle) they are generally below 5 cm. In addition, we have found that these observations agree well overall with predictions of accuracy derived from the JGM 2 variance-covariance matrix; the corresponding projected LLC errors for Geosat, T/P, and ERS 1 are usually between 1 and 4 cm, 1 – 2 cm, and 1 – 4 cm, respectively (they depend on the filtering of long-periodic perturbations and on the order of the LLC). This agreement is especially impressive for ERS 1 since no data of any kind from this mission was used in forming JGM 2.

The observed crossover differences for Geosat, T/P and ERS 1 are 8, 3, and 11 cm (rms), respectively. These observations also agree well with prediction of accuracy derived from the JGM 2 variance-covariance matrix; the corresponding projected crossover errors for Geosat and T/P are 8 cm and 2.3 cm, respectively. The precision of our mean difference observations is about 3 cm for Geosat (approx. 24,000 observations), 1.5 cm for T/P (approx. 6,000 observations) and 5 cm for ERS 1 (approx. 44,000 observations). Thus, these “global” independent data should provide a valuable new source for improving geopotential models. Our results show the need for further correction of the low order JGM 2 geopotential as well as certain resonant orders for all 3 satellites.     

Vertical TEC in the equatorial anomaly region from satellite altimetry data

The high-precision demands imposed by the ocean altimetry community of the late 1980 resulted in the TOPEX/Poseidon mission. This mission was the first to carry as its main instrument a dual-frequency sea-altimeter on board a satellite. This instrument together with other state-of-the-art technologies involved in the mission, led to sea-height determinations with precision better than 2 cm. As a by-product, the TOPEX/Poseidon mission provided vertical TEC determinations that since they became available, have demonstrated to be a powerful tool for ionospheric studies.     

Improved kHz-SLR tracking techniques and orbit quality analysis for LEO missions

Satellite gravity field missions such as CHAMP, GRACE and GOCE are designed as low Earth orbiting spacecraft (LEO) with orbit heights of about 250–500 km. The challenging mission objectives require a very precise knowledge of the satellite orbit position in space. For these missions precise orbit information is typically provided by GPS satellite-to-satellite tracking (SST) observations supported by satellite laser ranging (SLR).     

Characteristics of GOCE orbits based on Satellite Laser Ranging

The Gravity field and steady-state Ocean Circulation Explorer (GOCE) was the first European Space Agency’s (ESA) Earth Explorer core mission. Through its extremely low, about 260?km above the Earth, circular, sun-synchronous orbit, the satellite gained high spatial resolution and accuracy gravity gradient, and ocean circulation data. Global Positioning System (GPS) receivers, mounted on the spacecraft, allowed the determination of reduced-dynamic and kinematic GOCE orbits, whereas Laser Retroreflector Array (LRA) dedicated to Satellite Laser Ranging (SLR) allowed an independent validation of GPS-derived orbits. In this paper, residuals between different GPS-based orbit types and SLR observations are used to investigate the sensitivity and the influence of solar, geomagnetic, and ionospheric activities on the quality of kinematic and reduced-dynamic GOCE orbits. We also analyze the quality of data provided by individual SLR sites, by detecting time biases using ascending and descending sun-synchronous GOCE orbit passes, and the residual analysis of the measurement characteristics, i.e., the dependency of SLR residuals as a function of nadir and horizontal angles. Results show a substantial vulnerability of kinematic orbit solutions to the solar F10.7 index and the ionospheric activity measured by the variations of the Total Electron Content (TEC) values. The sensitivity of kinematic orbits to the three-hour-range KP index is rather minor. The reduced-dynamic orbits are almost insensitive to indices describing ionospheric, solar, and geomagnetic activities. The investigation of individual SLR sites shows that some of them are affected by time bias errors, whereas other demonstrate systematics, such as a dependency between observation residuals and the satellite nadir angle or the horizontal azimuth angle from the SLR station to the direction of the satellite.     

ERS-1 radial position refinement by dual crossover analysis with TOPEX/Poseidon

Achieving orbital accuracies in the radial direction for ERS-1 commensurate with those for TOPEX/Poseidon is of utmost importance for the integration of the two altimeter data sets. This paper outlines a procedure whereby the radial orbit error for ERS-1 is recovered as a time series expansion in the form of a finite Fourier series with additional terms for atmospheric drag, solar radiation pressure, and initial state vector mismodelling. Using a least squares collocation method with constraints derived from the JGM2 gravity field co-variance matrix, the radial error is recovered using both dual crossovers and ERS-1 single satellite crossovers. Aggregate arcs are then used to derive the ERS-1 orbit error over the repeat period of 35 days. The results are presented in the improvement of fit in the dual crossover, ERS-1 crossover and altimetry data sets as well as the recovery of an altimeter bias for the two satellites.     

DPOD2005: An extension of ITRF2005 for Precise Orbit Determination

For Precise Orbit Determination of altimetry missions, we have computed a data set of DORIS station coordinates defined for specific time intervals called DPOD2005. This terrestrial reference set is an extension of ITRF2005. However, it includes all new DORIS stations and is more reliable, as we disregard stations with large velocity formal errors as they could contaminate POD computations in the near future. About 1/4 of the station coordinates need to be defined as they do not appear in the original ITRF2005 realization. These results were verified with available DORIS and GPS results, as the integrity of DPOD2005 is almost as critical as its accuracy. Besides station coordinates and velocities, we also provide additional information such as periods for which DORIS data should be disregarded for specific DORIS stations, and epochs of coordinate and velocity discontinuities (related to either geophysical events, equipment problem or human intervention). The DPOD model was tested for orbit determination for TOPEX/Poseidon (T/P), Jason-1 and Jason-2. Test results show DPOD2005 offers improvement over the original ITRF2005, improvement that rapidly and significantly increases after 2005. Improvement is also significant for the early T/P cycles indicating improved station velocities in the DPOD2005 model and a more complete station set. Following 2005 the radial accuracy and centering of the ITRF2005-original orbits rapidly degrades due to station loss.     

The effect of geocenter motion on Jason-2 orbits and the mean sea level

We compute a series of Jason-2 GPS and SLR/DORIS-based orbits using ITRF2005 and the std0905 standards ( Lemoine et al., 2010). Our GPS and SLR/DORIS orbit data sets span a period of 2 years from cycle 3 (July 2008) to cycle 74 (July 2010). We extract the Jason-2 orbit frame translational parameters per cycle by the means of a Helmert transformation between a set of reference orbits and a set of test orbits. We compare the annual terms of these time-series to the annual terms of two different geocenter motion models where biases and trends have been removed. Subsequently, we include the annual terms of the modeled geocenter motion as a degree-1 loading displacement correction to the GPS and SLR/DORIS tracking network of the POD process. Although the annual geocenter motion correction would reflect a stationary signal in time, under ideal conditions, the whole geocenter motion is a non-stationary process that includes secular trends. Our results suggest that our GSFC Jason-2 GPS-based orbits are closely tied to the center of mass (CM) of the Earth consistent with our current force modeling, whereas GSFC’s SLR/DORIS-based orbits are tied to the origin of ITRF2005, which is the center of figure (CF) for sub-secular scales. We quantify the GPS and SLR/DORIS orbit centering and how this impacts the orbit radial error over the globe, which is assimilated into mean sea level (MSL) error, from the omission of the annual term of the geocenter correction. We find that for the SLR/DORIS std0905 orbits, currently used by the oceanographic community, only the negligence of the annual term of the geocenter motion correction results in a – 4.67 ± 3.40 mm error in the Z-component of the orbit frame which creates 1.06 ± 2.66 mm of systematic error in the MSL estimates, mainly due to the uneven distribution of the oceans between the North and South hemisphere.     

Precise orbit determination and baseline consistency assessment for Swarm constellation

Precise orbit determination (POD) and precise baseline determination (PBD) of Swarm satellites with 4 years of data are investigated. Ambiguity resolution (AR) plays a crucial role in achieving the best orbit accuracy. Swarm POD and PBD based on single difference (SD) AR and traditional double difference (DD) AR methods are explored separately. Swarm antenna phase center variation (PCV) corrections are developed to further improve the orbit determination accuracy. The code multipath of C1C, C1W and C2W observations is first evaluated and clear variations in code noise related to different receiver settings are observed. Carrier phase residuals of different time periods and different loop tracking settings of receiver are studied to explore the effect of ionospheric scintillation on POD. The reduction of residuals in the polar and geomagnetic equator regions confirms the positive impact of the updated carrier tracking loops (TLs) on POD performance. The SD AR orbits and orbits with float ambiguity (FA) are compared with the Swarm precise science orbits (PSOs). An average improvement of 27 %, 4 % and 16 % is gained in along-track, cross-track and radial directions by fixing the ambiguity to integer. For Swarm-A/B and Swarm-B/C formations, specific days are selected to perform the DD AR-based POD during which the average distance of the formation satellites is less than 5000 km. Satellite laser ranging (SLR) observations are employed to validate the performance of FA, SD AR and DD AR orbits. The consistency between the SD AR orbits and SLR data is at a level of 10 mm which shows an improvement of 25 % when comparing with the FA results. An SLR residuals reduction of 15 % is also achieved by the DD AR solution for the selected days. Precise relative navigation is also an essential aspect for spacecraft formation flying missions. The closure error method is proposed to evaluate the baseline precision in three dimensions. A baseline precision of 1–3 mm for Swarm-A/C formation and 3–5 mm for Swarm-A/B and Swarm-B/C satellite pairs is verified by both the consistency check and closure error method.     

Achievable debris orbit prediction accuracy using laser ranging data from a single station

Earlier studies have shown that an orbit prediction accuracy of 20 arc sec ground station pointing error for 1–2 day predictions was achievable for low Earth orbit (LEO) debris using two passes of debris laser ranging (DLR) data from a single station, separated by about 24 h. The accuracy was determined by comparing the predicted orbits with subsequent tracking data from the same station. This accuracy statement might be over-optimistic for other parts of orbit far away from the station. This paper presents the achievable orbit prediction accuracy using satellite laser ranging (SLR) data of Starlette and Larets under a similar data scenario as that of DLR. The SLR data is corrupted with random errors of 1 m standard deviation so that its accuracy is similar to that of DLR data. The accurate ILRS Consolidated Prediction Format orbits are used as reference to compute the orbit prediction errors. The study demonstrates that accuracy of 20 arc sec for 1–2 day predictions is achievable.