The impact of orbital errors on the estimation of satellite clock errors and PPP

Precise satellite orbit and clocks are essential for providing high accuracy real-time PPP (Precise Point Positioning) service. However, by treating the predicted orbits as fixed, the orbital errors may be partially assimilated by the estimated satellite clock and hence impact the positioning solutions. This paper presents the impact analysis of errors in radial and tangential orbital components on the estimation of satellite clocks and PPP through theoretical study and experimental evaluation. The relationship between the compensation of the orbital errors by the satellite clocks and the satellite-station geometry is discussed in details. Based on the satellite clocks estimated with regional station networks of different sizes (∼100, ∼300, ∼500 and ∼700 km in radius), results indicated that the orbital errors compensated by the satellite clock estimates reduce as the size of the network increases. An interesting regional PPP mode based on the broadcast ephemeris and the corresponding estimated satellite clocks is proposed and evaluated through the numerical study. The impact of orbital errors in the broadcast ephemeris has shown to be negligible for PPP users in a regional network of a radius of ∼300 km, with positioning RMS of about 1.4, 1.4 and 3.7 cm for east, north and up component in the post-mission kinematic mode, comparable with 1.3, 1.3 and 3.6 cm using the precise orbits and the corresponding estimated clocks. Compared with the DGPS and RTK positioning, only the estimated satellite clocks are needed to be disseminated to PPP users for this approach. It can significantly alleviate the communication burdens and therefore can be beneficial to the real time applications.

Pertubações orbitais devidas a maré terrestre

Pós-graduação em Física - FEG

Determinação de órbitas com o GPS através de mínimos quadrados recursivo com rotações de Givens

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Enhancement of Trajectory Determination of Orbiter Spacecraft by Using Pairs of Planetary Optical Images

The subject of the present thesis is about the enhancement of orbiter spacecraft navigation capabilities obtained by the standard radiometric link, taking advantage of an imaging payload and making use of a novel definition of optical measurements. An ESA Mission to Mercury called BepiColombo, was selected as a reference case for this study, and in particular its Mercury Planetary Orbiter (MPO), because of the presence of SIMBIO-SYS, an instrument suite part of the MPO payload, capable of acquiring high resolution images of the surface of Mercury. The use of optical measurements for navigation, can provide complementary informations with respect to Doppler, for enhanced performances or a relaxation of the radio tracking requisites in term of ground station schedule. Classical optical techniques based on centroids, limbs or landmarks, were the base to a novel idea for optical navigation, inspired by concepts of stereoscopic vision. In brief, the relation between two overlapped images acquired by a nadir pointed orbiter spacecraft at different times, was defined, and this information was then formulated into an optical measurement, to be processed by a navigation filter. The formulation of this novel optical observable is presented, moreover the analysis of the possible impact on the mission budget and images scheduling is addressed. Simulations are conducted using an orbit determination software already in use for spacecraft navigation in which the proposed optical measurements were implemented and the final results are given.

Design of a GNSS receiver for the ESEO mission

This thesis was carried out inside the ESA's ESEO mission and focus in the design of one of the secondary payloads carried on board the spacecraft: a GNSS receiver for orbit determination. The purpose of this project is to test the technology of the orbit determination in real time applications by using commercial components. The architecture of the receiver includes a custom part, the navigation computer, and a commercial part, the front-end, from Novatel, with COCOM limitation removed, and a GNSS antenna. This choice is motivated by the goal of demonstrating the correct operations in orbit, enabling a widespread use of this technology while lowering the cost and time of the device’s assembly. The commercial front-end performs GNSS signal acquisition, tracking and data demodulation and provides raw GNSS data to the custom computer. This computer processes this raw observables, that will be both transferred to the On-Board Computer and then transmitted to Earth and provided as input to the recursive estimation filter on-board, in order to obtain an accurate positioning of the spacecraft, using the dynamic model. The main purpose of this thesis, is the detailed design and development of the mentioned GNSS receiver up to the ESEO project Critical Design Review, including requirements definition, hardware design and breadboard preliminary test phase design.

Earth gravity field recovery using GPS, GLONASS, and SLR satellites

The time variable Earth’s gravity field provides the information about mass transport within the system Earth, i.e., the relationship of mass transport between atmosphere, oceans, and land hydrology. We recover the low-degree parameters of the time variable gravity field using microwave observations from GPS and GLONASS satellites and from SLR data to five geodetic satellites, namely LAGEOS-1/2, Starlette, Stella, and AJISAI. GPS satellites are particularly sensitive to specific coefficients of the Earth's gravity field, because of the deep 2:1 orbital resonance with Earth rotation (two revolutions of the GPS satellites per sidereal day). The resonant coefficients cause, among other, a “secular” drift (actually periodic variations of very long periods) of the semi-major axes of up to 5.3 m/day of GPS satellites. We processed 10 years of GPS and GLONASS data using the standard orbit models from the Center of Orbit Determination in Europe (CODE) with a simultaneous estimation of the Earth gravity field coefficients and other parameters, e.g., satellite orbit parameters, station coordinates, Earth rotation parameters, troposphere delays, etc. The weekly GNSS gravity solutions up to degree and order 4/4 are compared to the weekly SLR gravity field solutions. The SLR-derived geopotential coefficients are compared to monthly GRACE and CHAMP results.

Homogeneous reprocessing of GPS, GLONASS and SLR observations

The International GNSS Service (IGS) provides operational products for the GPS and GLONASS constellation. Homogeneously processed time series of parameters from the IGS are only available for GPS. Reprocessed GLONASS series are provided only by individual Analysis Centers (i. e. CODE and ESA), making it difficult to fully include the GLONASS system into a rigorous GNSS analysis. In view of the increasing number of active GLONASS satellites and a steadily growing number of GPS+GLONASS-tracking stations available over the past few years, Technische Universität Dresden, Technische Universität München, Universität Bern and Eidgenössische Technische Hochschule Zürich performed a combined reprocessing of GPS and GLONASS observations. Also, SLR observations to GPS and GLONASS are included in this reprocessing effort. Here, we show only SLR results from a GNSS orbit validation. In total, 18 years of data (1994–2011) have been processed from altogether 340 GNSS and 70 SLR stations. The use of GLONASS observations in addition to GPS has no impact on the estimated linear terrestrial reference frame parameters. However, daily station positions show an RMS reduction of 0.3 mm on average for the height component when additional GLONASS observations can be used for the time series determination. Analyzing satellite orbit overlaps, the rigorous combination of GPS and GLONASS neither improves nor degrades the GPS orbit precision. For GLONASS, however, the quality of the microwave-derived GLONASS orbits improves due to the combination. These findings are confirmed using independent SLR observations for a GNSS orbit validation. In comparison to previous studies, mean SLR biases for satellites GPS-35 and GPS-36 could be reduced in magnitude from −35 and −38 mm to −12 and −13 mm, respectively. Our results show that remaining SLR biases depend on the satellite type and the use of coated or uncoated retro-reflectors. For Earth rotation parameters, the increasing number of GLONASS satellites and tracking stations over the past few years leads to differences between GPS-only and GPS+GLONASS combined solutions which are most pronounced in the pole rate estimates with maximum 0.2 mas/day in magnitude. At the same time, the difference between GLONASS-only and combined solutions decreases. Derived GNSS orbits are used to estimate combined GPS+GLONASS satellite clocks, with first results presented in this paper. Phase observation residuals from a precise point positioning are at the level of 2 mm and particularly reveal poorly modeled yaw maneuver periods.

Processing 20 years of SLR observations to GNSS satellites

We process 20 years of SLR observations to GPS and GLONASS satellites using the reprocessed 3-day and 1-day microwave orbits provided by the Center for Orbit Determination in Europe (CODE) for the period 1994-2013. We study the dependency of the SLR residuals on the type, size, and a number of corner cubes in satellite laser reflector arrays (LRA). We show that the mean SLR residuals and the RMS of residuals depend on the coating of LRA and the block or type of GNSS satellites. The SLR mean residuals are also a function of the equipment used at SLR stations including detector types and detecting modes.

CODE’s new ultra-rapid orbit and ERP products for the IGS

The International GNSS Service (IGS) issues four sets of so-called ultra-rapid products per day, which are based on the contributions of the IGS Analysis Centers. The traditional (“old”) ultra-rapid orbit and earth rotation parameters (ERP) solution of the Center for Orbit Determination in Europe (CODE) was based on the output of three consecutive 3-day long-arc rapid solutions. Information from the IERS Bulletin A was required to generate the predicted part of the old CODE ultra-rapid product. The current (“new”) product, activated in November 2013, is based on the output of exactly one multi-day solution. A priori information from the IERS Bulletin A is no longer required for generating and predicting the orbits and ERPs. This article discusses the transition from the old to the new CODE ultra-rapid orbit and ERP products and the associated improvement in reliability and performance. All solutions used in this article were generated with the development version of the Bernese GNSS Software. The package was slightly extended to meet the needs of the new CODE ultra-rapid generation.

Impact of the arc length on GNSS analysis results

Homogeneously reprocessed combined GPS/GLONASS 1- and 3-day solutions from 1994 to 2013, generated by the Center for Orbit Determination in Europe (CODE) in the frame of the second reprocessing campaign REPRO-2 of the International GNSS Service, as well as GPS- and GLONASS-only 1- and 3-day solutions for the years 2009 to 2011 are analyzed to assess the impact of the arc length on the estimated Earth Orientation Parameters (EOP, namely polar motion and length of day), on the geocenter, and on the orbits. The conventional CODE 3-day solutions assume continuity of orbits, polar motion components, and of other parameters at the day boundaries. An experimental 3-day solution, which assumes continuity of the orbits, but independence from day to day for all other parameters, as well as a non-overlapping 3-day solution, is included into our analysis. The time series of EOPs, geocenter coordinates, and orbit misclosures, are analyzed. The long-arc solutions were found to be superior to the 1-day solutions: the RMS values of EOP and geocenter series are typically reduced between 10 and 40 %, except for the polar motion rates, where RMS reductions by factors of 2–3 with respect to the 1-day solutions are achieved for the overlapping and the non-overlapping 3-day solutions. In the low-frequency part of the spectrum, the reduction is even more important. The better performance of the orbits of 3-day solutions with respect to 1-day solutions is also confirmed by the validation with satellite laser ranging.

Atmospheric, oceanic and hydrological polar motion excitation mechanisms based on a combination of geometric and gravimetric space observations

The goal of our study is to determine accurate time series of geophysical Earth rotation excitations to learn more about global dynamic processes in the Earth system. For this purpose, we developed an adjustment model which allows to combine precise observations from space geodetic observation systems, such as Satellite Laser Ranging (SLR), Global Navigation Satellite Systems (GNSS), Very Long Baseline Interferometry (VLBI), Doppler Orbit determination and Radiopositioning Integrated on Satellite (DORIS), satellite altimetry and satellite gravimetry in order to separate geophysical excitation mechanisms of Earth rotation. Three polar motion time series are applied to derive the polar motion excitation functions (integral effect). Furthermore we use five time variable gravity field solutions from Gravity Recovery and Climate Experiment (GRACE) to determine not only the integral mass effect but also the oceanic and hydrological mass effects by applying suitable filter techniques and a land-ocean mask. For comparison the integral mass effect is also derived from degree 2 potential coefficients that are estimated from SLR observations. The oceanic mass effect is also determined from sea level anomalies observed by satellite altimetry by reducing the steric sea level anomalies derived from temperature and salinity fields of the oceans. Due to the combination of all geodetic estimated excitations the weaknesses of the individual processing strategies can be reduced and the technique-specific strengths can be accounted for. The formal errors of the adjusted geodetic solutions are smaller than the RMS differences of the geophysical model solutions. The improved excitation time series can be used to improve the geophysical modeling.

Resonance studies of artificial earth satellites

Orbit determination from artificial satellite observations is a key process in obtaining information about the Earth and its environment. A study of the perturbations experienced by these satellites enables knowledge to be gained of the upper atmosphere, the gravity field, ocean tides, solid-Earth tides and solar radiation. The gravity field is expressed as a double infinite series of associated Legendre functions (tesseral harmonics). In contemporary global gravity field models the overall geoid is well determined. An independent check on these gravity field harmonics of a particular order may be made by analysis of satellites that pass through resonance of that order. For such satellites the perturbations of the orbital elements close to resonance are analysed to derive lumped harmonic coefficients. The orbital parameters of 1984-106A have been determined at 43 epochs, during which time the satellite was close to 14th order resonance. Analysis of the inclination and eccentricity yielded 6 lumped harmonic coefficients of order 14 whilst analysis of the mean motion yielded additional pairs of lumped harmonics of orders 14, 28 and 42, with the 14th order harmonics superseding those obtained from analysis of the inclination. This thesis concentrates in detail on the theoretical changes of a near-circular satellite orbit perturbed by the Earth's gravity field under the influence of minimal air-drag whilst in resonance with the Earth. The satellite 1984-106A experienced the interesting property of being temporarily trapped with respect to a secondary resonance parameter due to the low air-drag in 1987. This prompted the theoretical investigation of such a phenomenon. Expressions obtained for the resonance parameter led to the determination of 8 lumped harmonic coefficients, coincidental to those already obtained. All the derived lumped harmonic values arc used to test the accuracy of contemporary gravity field models and the underlying theory in this thesis.

Modeling, Simulation, and Characterization of Space Debris in low-Earth Orbit

Every space launch increases the overall amount of space debris. Satellites have limited awareness of nearby objects that might pose a collision hazard. Astrometric, radiometric, and thermal models for the study of space debris in low-Earth orbit have been developed. This modeled approach proposes analysis methods that provide increased Local Area Awareness for satellites in low-Earth and geostationary orbit. Local Area Awareness is defined as the ability to detect, characterize, and extract useful information regarding resident space objects as they move through the space environment surrounding a spacecraft. The study of space debris is of critical importance to all space-faring nations. Characterization efforts are proposed using long-wave infrared sensors for space-based observations of debris objects in low-Earth orbit. Long-wave infrared sensors are commercially available and do not require solar illumination to be observed, as their received signal is temperature dependent. The characterization of debris objects through means of passive imaging techniques allows for further studies into the origination, specifications, and future trajectory of debris objects. Conclusions are made regarding the aforementioned thermal analysis as a function of debris orbit, geometry, orientation with respect to time, and material properties. Development of a thermal model permits the characterization of debris objects based upon their received long-wave infrared signals. Information regarding the material type, size, and tumble-rate of the observed debris objects are extracted. This investigation proposes the utilization of long-wave infrared radiometric models of typical debris to develop techniques for the detection and characterization of debris objects via signal analysis of unresolved imagery. Knowledge regarding the orbital type and semi-major axis of the observed debris object are extracted via astrometric analysis. This knowledge may aid in the constraint of the admissible region for the initial orbit determination process. The resultant orbital information is then fused with the radiometric characterization analysis enabling further characterization efforts of the observed debris object. This fused analysis, yielding orbital, material, and thermal properties, significantly increases a satellite’s Local Area Awareness via an intimate understanding of the debris environment surrounding the spacecraft.

Lightcurve inversion for asteroid spins and shapes

Knowledge of the physical properties of asteroids is crucial in many branches of solar-system research. Knowledge of the spin states and shapes is needed, e.g., for accurate orbit determination and to study the history and evolution of the asteroids. In my thesis, I present new methods for using photometric lightcurves of asteroids in the determination of their spin states and shapes. The convex inversion method makes use of a general polyhedron shape model and provides us at best with an unambiguous spin solution and a convex shape solution that reproduces the main features of the original shape. Deriving information about the non-convex shape features is, in principle, also possible, but usually requires a priori information about the object. Alternatively, a distribution of non-convex solutions, describing the scale of the non-convexities, is also possible to be obtained. Due to insufficient number of absolute observations and inaccurately defined asteroid phase curves, the $c/b$-ratio, i.e., the flatness of the shape model is often somewhat ill-defined. However, especially in the case of elongated objects, the flatness seems to be quite well constrained, even in the case when only relative lightcurves are available. The results prove that it is, contrary to the earlier misbelief, possible to derive shape information from the lightcurve data if a sufficiently wide range of observing geometries is covered by the observations. Along with the more accurate shape models, also the rotational states, i.e., spin vectors and rotation periods, are defined with improved accuracy. The shape solutions obtained so far reveal a population of irregular objects whose most descriptive shape characteristics, however, can be expressed with only a few parameters. Preliminary statistical analyses for the shapes suggests that there are correlations between shape and other physical properties, such as the size, rotation period and taxonomic type of the asteroids. More shape data of, especially, the smallest and largest asteroids, as well as the fast and slow rotators is called for in order to be able to study the statistics more thoroughly.

High-accuracy same-beam VLBI observations using Shanghai and Urumqi telescopes

The same-beam VLBI observations of Rstar and Vstar, which were two small satellites of Japanese lunar mission, SELENE, were successfully performed by using Shanghai and Urumqi 25-m telescopes. When the separation angle between Rstar and Vstar was less than 0.1 deg, the differential phase delay of the X-band signals between Rstar and Vstar on Shanghai-Urumqi baseline was obtained with a very small error of 0.15 mm rms, which was reduced by 1-2 order compared with the former VLBI results. When the separation angle was less than 0.56 deg, the differential phase delay of the S-band signals was also obtained with a very small error of several mm rms. The orbit determination for Rstar and Vstar was performed, and the accuracy was improved to a level of several meters by using VLBI and Doppler data. The high-accuracy same-beam differential VLBI technique is very useful in orbit determination for a spacecraft, and will be used in orbit determination for Mars missions of China Yinghuo-1 and Russia Phobos-grunt.