Retroreflectors for space traffic management.

Since the beginning of satellite altimetry missions, the ocean surface topography community requires precise and accurate satellite orbits. With start of the satellite Sentinel-3A on February 16, 2016, two radio- and one altimeter onboard the satellite accomplish the Copernicus program with a ocean- and land monitoring mission, planned for a nominal mission lifetime of 7 years. The satellite is orbiting the Earth on a polar, Sun-synchronous trajectory at an altitude of 815 km. For the purpose of precise orbit determination, the satellite is equipped with a geodetic-grade dual-frequency Global Positioning System (GPS) receiver. The GPS measurements are employed together with a set of gravitational and non-gravitational models in a Reduced-Dynamic Orbit Determination (RDOD) approach, which combines the advantages of a dynamic and a kinematic positioning for deriving precise satellite orbits. However, especially the non-gravitational force models require sophisticated modeling techniques. Therefore, a satellite macro model is introduced, which allows a proper modeling of accelerations due to Solar Radiation Pressure (SRP), Earth Radiation Pressure (ERP), and atmospheric drag. Especially the Sun-synchronous orbit, and the huge solar array, which is loosely coupled to the satellite body, makes the precise orbit determination challenging. As members of the Copernicus Quality Working Group, the Astronomical Institute of the University of Bern (AIUB), and the German Aerospace Center (DLR) are, among others, responsible for the orbit validation of Sentinel-3A. Both groups make use of a satellite macro model within a reduced-dynamic approach but differ in the employed software solutions and the pseudo-stochastic modeling. Basically, the pseudo-stochastic parameters allow to compensate potential deficits in the employed force models. Within this poster, the modeling aspects are briefly introduced, followed by comparing the results of both groups. The results include the estimated empirical accelerations, and the estimated satellite orbits of Sentinel-3A. Furthermore, the orbit quality is assessed by Satellite Laser Ranging (SLR), an external and independent tool for orbit validation.

<p>Numerous active low Earth orbiters (LEOs) and Global Navigation Satellite System (GNSS) satellites, including the Galileo constellation, are equipped with laser retroreflectors used for Satellite Laser Ranging (SLR). Moreover, most of LEOs are equipped with GNSS receivers for precise orbit determination. SLR measurements to LEOs, GNSS, and geodetic satellites vary in terms of the number of registered normal points (NPs) or registered satellite passes. In 2016-2018, SLR measurements to LEOs constituted 81% of all NPs and 59% of all registered satellite passes, whereas 10% of NPs and 30% of satellite passes, respectively, were assigned to GNSS. The remaining SLR measurements were completed by geodetic satellites, including LAGEOS-1/2, and LARES-1.</p><p>In this study, we show that the SLR observations to Galileo, passive geodetic and active LEO satellites together with precise GNSS-based orbits of LEOs and Galileo, can be used for the determination of global geodetic parameters, such as geocenter coordinates (GCC) and Earth rotation parameters (ERPs), i.e. pole coordinates, and length-of-day parameter.</p><p>GCC are typically determined using SLR observations to passive geodetic satellites, such as LAGEOS-1/2. Also, the SLR observations to LAGEOS-1/2 together with GNSS and Very Long Baseline Interferometry data are used for the determination of ERPs. Here, we use SLR observations to Galileo, LAGEOS-1/2, LARES-1, Sentinel-3A, SWARM-A/B/C, TerraSAR-X, Jason-2, GRACE-A/B satellites to investigate whether they can be applied for the reference frame realization and for deriving high-quality global geodetic parameters.</p><p>We present various types of solutions to investigate the best solution set-up. The studied solutions differ in terms of solution lengths, the combination of different sets of satellites and the relative weights for the variance scaling factors of technique and satellite-specific normal equations. We compare our results with the standard LAGEOS-based solutions, the combined EOP-14-C04 products and show the consistency of the results.</p>

Due to the continued development of the GLONASS satellites, precise orbit determination (POD) still poses a series of challenges. This study examines the impact of introducing the analytical tube-wing model for GLONASS-M and the box-wing model for GLONASS-K in a series of hybrid POD strategies that consider both the analytical model and a series of empirical parameters. We assess the perturbing accelerations acting on GLONASS spacecraft based on the analytical model. All GLONASS satellites are equipped with laser retroreflectors for satellite laser ranging (SLR). We apply the SLR observations for the GLONASS POD in a series of GNSS + SLR combined solutions. The application of the box-wing model significantly improves GLONASS orbits, especially for GLONASS-K, reducing the STD of SLR residuals from 92.6 to 27.6 mm. Although the metadata for all GLONASS-M satellites reveal similar construction characteristics, we found differences in empirical accelerations and SLR offsets not only between GLONASS-M and GLONASS-M+ but also within the GLONASS-M+ series. Moreover, we identify satellites with inferior orbit solutions and check if we can improve them using the analytical model and SLR observations. For GLONASS-M SVN730, the STD of the SLR residuals for orbits determined using the empirical solution is 48.7 mm. The STD diminishes to 41.2 and 37.8 mm when introducing the tube-wing model and SLR observations, respectively. As a result, both the application of the SLR observations and the analytical model significantly improve the orbit solution as well as reduce systematic errors affecting orbits of GLONASS satellites.

The challenge of precise orbit determination for STSAT-2C using extremely sparse SLR data

Precise Orbit Determination (POD) is an integral part for analyzing measurements from space geodetic techniques such as Satellite Laser Ranging (SLR) and Global Navigation Satellite Systems (GNSS) such as the Global Positioning System (GPS). In the last two decades, POD based on GPS data has furthermore been established as one of the standard techniques to derive trajectories of satellites in the low Earth orbit (LEO) with highest accuracy. Since the launch of dedicated gravity missions, GPS sensors are not only used as a key tracking system for LEO POD, but also for extracting the long wavelength part of the Earth’s gravity field (together with SLR to spherical satellites). This chapter introduces SLR and GNSS measurements collected by the terrestrial networks of the International Laser Ranging Service (ILRS) and the International GNSS Service (IGS) as the observational basis for the realization of a terrestrial reference frame from satellite data. On this foundation, the basic equations and mathematical methods of orbit determination are introduced and extensively discussed. Pseudo-stochastic orbit modeling techniques are eventually presented as a general and efficient concept to determine satellite trajectories of highest quality even in presence of deficient force models, covering the full range between dynamic and purely kinematic solutions. Selected results from the application of the discussed orbit determination techniques are highlighted for GPS LEO data. Special emphasis is also put to present orbit determination in the context of more general orbit determination problems, where satellite trajectories are simultaneously determined with other parameters encompassing (at maximum) all pillars of geodesy, i.e., the shape, rotation, and gravity field of the Earth.

Geostationary satellites (GEOs) play a significant role in the regional satellite navigation system. Simulation experiments show that the clock corrections could be mitigated through a single strategy or double differencing strategies for a navigation constellation, but for the mode of individual GEO orbit determination, high precision orbit and clock correction could not be obtained in the orbit determination based on the pseudorange data. A new GEO combined precise orbit determination (POD) strategy is studied in this paper, which combines pseudorange data and C-band transfer ranging data. This strategy overcomes the deficiency of C-band transfer ranging caused by limited stations and tracking time available. With the combination of transfer ranging and pseudorange data, clock corrections between the GEO and the stations can be estimated simultaneously along with orbital parameters, maintaining self-consistency between the satellite ephemeris and clock correction parameters. The error covariance analysis is conducted to illuminate the contributions from the transfer ranging data and the psudoranging data. Using data collected for a Chinese GEO satellite with 3 C-band transfer ranging stations and 4 L-band pseudorange tracking stations, POD experiments indicate that a meter-level accuracy is achievable. The root-mean-square (RMS) of the post-fit C-band ranging data is about 0.203 m, and the RMS of the post-fit pseudorange is 0.408 m. Radial component errors of the POD experiments are independently evaluated with the satellite laser ranging (SLR) data from a station in Beijing, with the residual RMS of 0.076 m. The SLR evaluation also suggests that for 2-h orbital predication, the predicted radial error is about 0.404 m, and the clock correction error is about 1.38 ns. Even for the combination of one C-band transfer ranging station and 4 pseudorange stations, POD is able to achieve a reasonable accuracy with the radial error of 0.280 m and the 2-h predicted radial error of 0.888 m. Clock synchronization between the GEO and tracking stations is achieved with an estimated accuracy of about 1.55 ns, meeting the navigation service requirements.

Solar radiation pressure (SRP) is an extremely critical perturbative force that affects the GNSS satellites’ precise orbit determination (POD). Its imperfect modelling is one of the main error sources of POD, whose magnitude is even to10−9 m/s2. The shadow factor (i.e., eclipse factor) is one crucial parameter of SRP, generally estimated by the cylindrical model, the conical model, or shadow models considering the Earth’s oblateness and the atmospheric effect, such as the Perspective Projection Method atmosphere (PPMatm) model and Solar radiation pressure with Oblateness and Lower Atmospheric Absorption, Refraction, and Scattering Curve Fit (SOLAARS-CF) model. This paper applies the former four shadow models to determine the corresponding precise orbit using BeiDou satellites’ ground-based observation, and then compared and assessed the orbit accuracy through Satellite Laser Ranging (SLR) validation and Inter-Satellite Link (ISL) check. The results show that the PPMatm model’s accuracy is equivalent to the SOLAARS-CF model. Compared with the conical shadow model, SLR validations show the orbit accuracy from the PPMatm and SOLAARS-CF model can be generally improved by 2–10 mm; ISL range check shows that the Root Mean Square (RMS) can be decreased by 2–7 mm. These results show that the shadow model in GNSS POD should fully consider the Earth’s oblateness and the atmospheric effect, especially for the perturbative acceleration higher than 10–10 m/s2.Graphical

Precise orbit determination for geostationary satellites with multiple tracking techniques

TerraSAR-X, to be launched at the end of May, 2007, carries the tracking, occultation and ranging (TOR) Category A payload instrument package. The TOR consists of a high-precision dual-frequency GPS receiver, called integrated GPS occultation receiver (IGOR), for precise orbit determination and atmospheric sounding and a laser retro-reflector (LRR) serving as target for the global satellite laser ranging (SLR) ground station network. The TOR is supplied by the GeoForschungsZentrum Potsdam (GFZ) Germany, and the Center for Space Research (CSR), Austin, Texas. The objective of the German/US collaboration is twofold: provision of atmospheric profiles for use in numerical weather predictions and climate studies from the occultation data and precision SAR data processing based on precise orbits and atmospheric products. GFZ also supplies the TOR instrumentation for the follow-on mission TanDEM-X. For the scientific objectives of the TanDEM- X mission, i.e. bi-static SAR together with TerraSAR-X, the dual- frequency GPS receiver is of vital importance for the millimeter level determination of the baseline between the two spacecrafts. In the following, the TOR instrumentation is characterized by the features of and accuracies achievable with the GPS receiver and the LRR. The data flow is shown, as GFZ operates a high latitude data receiving station for fast data access, and runs a data system for preprocessing, archiving and retrieval. With this data system, higher level products are generated for TerraSAR- X, including near real-time orbits by GFZ and post-processed precise science orbits by GFZ and CSR. In the case of TanDEM- X, GFZ operationally determines the inter-satellite baseline on the millimeter level. CSR produces calibration orbits to support the GFZ operational activity. Orbit results in terms of latency and accuracy from current missions are given and achievable results for TerraSAR-X and TanDEM-X are discussed.

Satellite laser ranging (SLR) to the satellites of the global navigation satellite systems (GNSS) provides substantial and valuable information about the accuracy and quality of GNSS orbits and allows for the SLR-GNSS colocation in space. In the framework of the NAVSTAR-SLR experiment two GPS satellites of Block-IIA were equipped with laser retroreflector arrays (LRAs), whereas all satellites of the GLONASS system are equipped with LRAs in an operational mode. We summarize the outcome of the NAVSTAR-SLR experiment by processing 20 years of SLR observations to GPS and 12 years of SLR observations to GLONASS satellites using the reprocessedmicrowave orbits provided by the center for orbit determination in Europe (CODE). The dependency of the SLR residuals on the size, shape, and number of corner cubes in LRAs is studied. We show that the mean SLR residuals and the RMS of residuals
\ndepend on the coating of the LRAs and the block or type of GNSS satellites. The SLR mean residuals are also a function of the equipment used at SLR stations including the single-photon and multi-photon detection modes.We also showthat the SLR observations to GNSS satellites are important to validate GNSS orbits and to assess deficiencies in the solar radiation pressure models. We found that the satellite signature effect, which is defined as a spread of optical pulse signals due to reflection from multiple reflectors, causes the variations of mean SLR residuals of up to 15mmbetween the observations at nadir angles of 0◦ and 14◦. in case of multiphoton SLR stations. For single-photon SLR stations this effect does not exceed 1 mm. When using the new empirical
\nCODE orbit model (ECOM), the SLR mean residual falls into the range 0.1–1.8 mm for high-performing singlephoton SLR stations observing GLONASS-M satellites with uncoated corner cubes. For best-performing multi-photon stations the mean SLR residuals are between −12.2 and
\n−25.6 mm due to the satellite signature effect.

ADEOS-II precise orbit determinations with GPS and SLR

Precise orbit determination and rapid orbit recovery supported by time synchronization

In the last decade, we have been witnessing a rapid development of the constellations of Global and Regional Navigation Satellite Systems (GNSS/RNSS). Besides the well-known GPS and GLONASS, newly developed systems such as Galileo, BeiDou, QZSS and NAVIC have become increasingly important. All satellites of new GNSS are equipped with laser retroreflector arrays (LRA) dedicated to Satellite Laser Ranging (SLR). SLR allows, e.g., an independent validation of microwave-based orbit products. Therefore, a fully operational online service called the multi-GNSS Orbit Validation Visualizer Using SLR (GOVUS) has been developed allowing for near real-time analysis of the quality of multi-GNSS orbits. The mean offsets of SLR residuals for Center for Orbit Determination in Europe (CODE) orbits in 2016 are at the level of −8, −38, −14, and −107 mm, for BeiDou, Galileo, GLONASS, and QZSS, respectively, with the standard deviations of 66, 36, 29, and 100 mm. Moreover, GOVUS can be used as a database containing information on equipment used at SLR stations and multi-GNSS satellite parameters. This paper includes a comprehensive description of the functionality and the structure of the developed service with exemplary analyses. The paper points out the most critical issues, limitations and challenges of multi-GNSS and SLR tracking network in the context of the SLR orbit validation. The goal of the paper and GOVUS itself is to determine: (1) what is the current quality of multi-GNSS orbits validated using SLR results; (2) what kinds of systematic errors can affect GNSS orbits and SLR observations; and (3) how to provide the online analysis tools to the broadest possible multi-GNSS community. The service has been officially operating since March 2017 as the Associate Analysis Center of the International Laser Ranging Service (ILRS ACC).

The inter-satellite link (ISL) equipment of the third-generation BeiDou Navigation Satellite System (BDS-3) can enhance the accuracy of precise orbit determination (POD) for BDS-3 satellites. However, the hardware delays within ISL observations can impact the measurement accuracy of ISL observations. This study evaluates the characteristics of hardware delays in ISL observations for BDS-3 medium Earth orbit (MEO) satellites and investigates their effects on POD. Analysis shows that most links exhibit slow variations in hardware delays, with standard deviations (STD) of 6.4 cm and 3.6 cm over 45 days for the C19-C25 and C20-C21 links, respectively. Conversely, some links like C30-C23 and C30-C33 show more pronounced variations, with STD values reaching 17.7 cm and 14.1 cm within three days. Links involving satellites C23 and C30 display higher variability and significant jumps in hardware delays compared with other links. The larger variations in hardware delays of these links can be absorbed by other parameters, thereby affecting the accuracy of parameter estimation. To mitigate their impact on POD, a segmented estimation strategy for hardware delays is proposed. This strategy divides a single link into multiple segments based on a set interval length and estimates the hardware delay parameters for each segment separately, thereby absorbing the variations in hardware delays or unmodeled errors. The effectiveness of this strategy is demonstrated through ISL POD residual analysis, comparison with L-band orbit determination results, and validation using Satellite Laser Ranging (SLR) residuals. The root mean square (RMS) of ISL observation residuals decreased by approximately 44.3%, from 6.1 cm to 3.4 cm. Comparisons with L-band orbits confirmed an approximate 28.2% reduction in the 3D RMS of orbit discrepancies, from 26.5 cm to 19.0 cm. Additionally, the RMS of the overall SLR residuals showed a slight decrease. Overall, the improved segmented estimation strategy effectively reduces the influence of hardware delay variations on BDS-3 orbit determination, notably enhancing the orbit determination accuracy of BDS-3 satellites, especially for C23 and C30, which experience larger hardware delay fluctuations.

Satellite laser ranging (SLR) has proven a very efficient method for contributing to the tracking of altimetric satellites and determining accurately their orbit although hampered by the non-worldwide coverage and the meteorological conditions. Indeed, in some cases it is the only method available to determine the satellite orbit (e.g., the orbits of the ERS-1 and Geosat-Follow-On missions). Moreover, any operational and non-weather dependent techniques, like GPS, DORIS, PRARE, can exhibit systematic errors in positioning and orbitography. A comparison with SLR results allows to evidence such errors and vice versa. For doing that, two different approaches for determining precise orbits can be considered: one based on global orbit determination, the other on a short-arc technique used to locally improve a global orbit determined by another tracking techniques, such as DORIS or GPS. We can thus validate a global orbit and achieve orbit quality control to a level of 2 to 3 centimeters at present and expect to achieve a level of 1 to 2 centimeters in the near future. Errors induced by station coordinates or by the gravity field (geographically correlated errors, for example) can be estimated from SLR tracking data. Colocation experiments with different techniques in the same geodetic site play also a key role to ensure precise relationships between the geodetic reference frames linked to each technique. In particular, the role of the SLR technique is to strengthen the vertical component (including velocity) of the positioning, which is crucial for altimetry missions. The role of SLR data in the modelling of the first terms of the gravity field has finally to be emphasized, which is of primary importance in orbitography, whatever the tracking technique used. Another application of SLR technology is the satellite altimeter calibration. Examples of past calibration and future experiments are given, including the accuracy we can expect from the Jason-1 and EnviSat space oceanography missions.