Pseudo-Stochastic Orbit Modeling Techniques for Low-Earth Orbiters

Total Page:16

File Type:pdf, Size:1020Kb

Pseudo-Stochastic Orbit Modeling Techniques for Low-Earth Orbiters View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by RERO DOC Digital Library J Geod (2006) 80: 47–60 DOI 10.1007/s00190-006-0029-9 ORIGINAL ARTICLE A. Jäggi · U. Hugentobler · G. Beutler Pseudo-stochastic orbit modeling techniques for low-Earth orbiters Received: 1 March 2005 / Accepted: 24 January 2006 / Published online: 30 March 2006 © Springer-Verlag 2006 Abstract The Earth’s non-spherical mass distribution and Keywords Low-Earth orbiter (LEO) · Precise orbit atmospheric drag cause the strongest perturbations on very determination (POD) · Pseudo-stochastic orbit modeling · low-Earth orbiting satellites (LEOs). Models of gravitational GPS and non-gravitational accelerations are utilized in dynamic precise orbit determination (POD) with GPS data, but it is 1 Introduction also possible to derive LEO positions based on GPS pre- cise point positioning without dynamical information. We Gravitational forces related to the Earth’s oblateness and the use the reduced-dynamic technique for LEO POD, which non-gravitational forces due to residual atmospheric densities combines the geometric strength of the GPS observations are responsible for the strongest perturbations acting on sat- with the force models, and investigate the performance of ellites in very low-Earth orbits. These forces produce large different pseudo-stochastic orbit parametrizations, such as periodic variations in the orbit and a decay in the orbital instantaneous velocity changes (pulses), piecewise constant height. Radial orbit differences caused by atmospheric drag accelerations, and continuous piecewise linear accelerations. after one day may reach 150 m for an initial height of 500 km The estimation of such empirical orbit parameters in a stan- (Beutler 2004) during normal solar activity. In addition, ma- dard least-squares adjustment process of GPS observations, jor perturbations induced by higher-order terms of the Earth’s together with other relevant parameters, strives for the high- gravity field and minor perturbations induced by solar radia- est precision in the computation of LEO trajectories. We used tion pressure pose a challenge for precise orbit determination the procedures for the CHAMP satellite and found that the (POD) for low-Earth orbiters (LEOs). orbits may be validated by means of independent SLR mea- The determination of precise orbits using space geo- surements at the level of 3.2 cm RMS. Validations with inde- detic techniques, such as Satellite Laser Ranging (SLR), has pendent accelerometer data revealed correlations at the level opened a new era for LEO POD and related applications. It of 95% in the along-track direction. As expected, the empir- has been shown that geophysically relevant information such ical parameters compensate to a certain extent for deficien- as geocenter variations or gravity field coefficients can be ex- cies in the dynamic models. We analyzed the capability of tracted by combining many years of observations to geodetic pseudo-stochastic parameters for deriving information about satellites (Reigber 1989). LEO POD is also a necessity for the mismodeled part of the force field and found evidence that missions of Earth observing satellites, e.g., the determina- the resulting orbits may be used to recover force field param- tion of ocean topography from the TOPEX/Poseidon mis- eters, if the number of pseudo-stochastic parameters is large sion (Fu et al. 1994). This was the first extensive use of the enough. Results based on simulations showed a significantly Global Positioning System (GPS) with a spaceborne receiver better performance of acceleration-based orbits for gravity for LEO POD (Bertiger et al. 1994). The geometric strength field recovery than for pulse-based orbits, with a quality com- of continuously collected GPS observations allowed for orbit parable to a direct estimation if unconstrained accelerations determination to closer than 3 cm in the radial direction. This are set up every 30 s. mission stimulated a number of follow-up satellite missions to carry on-board GPS receivers for POD, but also for other purposes such as atmospheric sounding (e.g., Kursinski et al. B · · A. Jäggi ( ) U. Hugentobler G. Beutler 1997). Astronomical Institute, University of Bern, Sidlerstrasse 5, 3012 Bern, Switzerland The launch of the CHAllenging Minisatellite Payload E-mail: [email protected] (CHAMP) on July 15, 2000, opened a new era in process- Tel.: +41-31-6318592 ing GPS data from spaceborne receivers. The combined Fax: +41-31-6313869 analysis of non-conservative accelerations, which were 48 A. Jäggi et al. (k) (k) directly measured with unprecedented accuracy by an on- with initial conditions r (t0) = r (a, e, i,,ω,T0; t0), board accelerometer (Touboul et al. 1999), and the GPS k = 0, 1 (level of time differentiation). The parameters tracking data proved the feasibility to separate gravitational a, e, i,,ω,T0 are the six Keplerian elements pertaining to from non-gravitational perturbations for gravity field esti- epoch t0. q1,...,qd denote additional dynamical parame- mation (Reigber et al. 2002). The large number of high qual- ters considered as unknowns, e.g., scaling factors of analyti- ity global gravity field models underscores the success of cally known accelerations, which describe deterministically CHAMP (e.g., Reigber et al. 2003, 2004). Current (GRACE) the perturbing acceleration acting on the satellite. and future (GOCE) gravity field recovery missions, which Let us assume that an a priori orbit r0(t) is available, e.g., aim to recover the gravity field with unprecedented accuracy, from a GPS code solution. Dynamic orbit determination may therefore pose the highest requirements on POD, e.g., 2 cm then be considered as an orbit improvement process, i.e., the RMS in each direction for GOCE (ESA 1999). actual orbit r(t) is expressed as a truncated Taylor series with There are many different studies for CHAMP POD based respect to the unknown orbit parameters pi about the a priori on GPS tracking data (e.g., Švehla and Rothacher 2002; van orbit, which is represented by the parameter values pi0: den IJssel et al. 2003; König et al. 2004). Most of them are n ∂ ( ) based on a technique called “reduced-dynamic” (Wu et al. r(t) = r (t) + r0 t · (p − p ), (2) 0 ∂pi i i0 1991), which exploits the geometric strength of the GPS i=1 observations to reduce the dependency on dynamic models where n = 6 + d denotes the total number of unknown orbit (see Sects. 5 and 6). Due to CHAMP’s very low altitude parameters, i.e., the six initial conditions for position and (below 450 km), it is common practice to solve for a large velocity and d dynamical parameters. Equation (2) allows us number of empirical parameters to compensate for defi- to improve the a priori orbit provided that the unknown orbit ciencies in the dynamic models. An alternative kinematic parameters and the partial derivatives of the a priori orbit with approach was successfully demonstrated for CHAMP by respect to those parameters are known. Švehla and Rothacher (2004), which yields precise ephem- eris entirely by geometric means. Such positions are attractive mainly for gravity field estimation because they are indepen- 2.1 Variational equations dent of any dynamic models. Gerlach et al. (2003) used kine- matic positions and accelerometer data and showed that grav- Let us assume that p is one of the parameters defining the ity field models can be estimated with a quality comparable to initial values or the dynamics in Eq. (1) and that the par- the official CHAMP models by means of the energy integral tial derivative of the a priori orbit r0(t) with respect to the method (O’Keefe 1957). The use of precise ephemeris for parameter p is designated by the function gravity field recovery had a considerably stimulating impact . ∂r0(t) on several groups (e.g., Mayer-Gürr et al. 2004) due to less zp(t) = . (3) demanding computational resources as in the case of classical ∂p numerical integration techniques (e.g., Visser et al. 2003b). The initial value problem associated with the partial deriv- In the following we consider pseudo-stochastic orbit ative in Eq. (3) is referred to as the system of variational modeling techniques for reduced-dynamic LEO POD as a equations in this article, and is obtained by taking the partial powerful and efficient method to derive most precise satel- derivative of Eq. (1). The result may be written as lite trajectories in the very low-Earth orbit (e.g., Visser and ∂ IJssel 2003a). Furthermore, we discuss different orbit repre- f1 z¨ p = A0 · zp + A1 · z˙ p + , (4) sentations and question whether such methodologies can be ∂p used to generate reduced-dynamic orbits, which could serve where the 3×3 matrices A and A are defined by as input for a subsequent gravity field estimation process. 0 1 Section 2 reviews the methods of dynamic LEO POD. ∂ fi ∂ fi A0[i,k] = ; A1[i,k] = ; (5) Section 3 presents two widely-used pseudo-stochastic orbit ∂r0,k ∂r˙0,k parametrizations and develops the mathematical background where fi denotes the i-th component of the total accelera- for a third, more refined orbit parametrization. CHAMP orbit ∈{ , , ,,ω, } results are shown in Sect. 4. Section 5 analyzes the presented tion f in Eq. (1). For p a e i T0 Eq. (4) is a lin- methods in a more detailed way in a simulation environ- ear, homogeneous, second-order differential equation system with initial values zp(t0) = 0 and z˙ p(t0) = 0, whereas for ment and focuses, together with Sect. 6, on a possible use of ∈{ ,..., } reduced-dynamic orbits for gravity field estimation. p q1 qd Eq. (4) is inhomogeneous with zero initial values. Note that in the latter case, the homogeneous part of Eq.
Recommended publications
  • Orbit Determination Using Modern Filters/Smoothers and Continuous Thrust Modeling
    Orbit Determination Using Modern Filters/Smoothers and Continuous Thrust Modeling by Zachary James Folcik B.S. Computer Science Michigan Technological University, 2000 SUBMITTED TO THE DEPARTMENT OF AERONAUTICS AND ASTRONAUTICS IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF MASTER OF SCIENCE IN AERONAUTICS AND ASTRONAUTICS AT THE MASSACHUSETTS INSTITUTE OF TECHNOLOGY JUNE 2008 © 2008 Massachusetts Institute of Technology. All rights reserved. Signature of Author:_______________________________________________________ Department of Aeronautics and Astronautics May 23, 2008 Certified by:_____________________________________________________________ Dr. Paul J. Cefola Lecturer, Department of Aeronautics and Astronautics Thesis Supervisor Certified by:_____________________________________________________________ Professor Jonathan P. How Professor, Department of Aeronautics and Astronautics Thesis Advisor Accepted by:_____________________________________________________________ Professor David L. Darmofal Associate Department Head Chair, Committee on Graduate Students 1 [This page intentionally left blank.] 2 Orbit Determination Using Modern Filters/Smoothers and Continuous Thrust Modeling by Zachary James Folcik Submitted to the Department of Aeronautics and Astronautics on May 23, 2008 in Partial Fulfillment of the Requirements for the Degree of Master of Science in Aeronautics and Astronautics ABSTRACT The development of electric propulsion technology for spacecraft has led to reduced costs and longer lifespans for certain
    [Show full text]
  • Assessment of Earth Gravity Field Models in the Medium to High Frequency Spectrum Based on GRACE and GOCE Dynamic Orbit Analysis
    Article Assessment of Earth Gravity Field Models in the Medium to High Frequency Spectrum Based on GRACE and GOCE Dynamic Orbit Analysis Thomas D. Papanikolaou 1,2,* and Dimitrios Tsoulis 3 1 Geoscience Australia, Cnr Jerrabomberra Ave and Hindmarsh Drive, Symonston, ACT 2609, Australia 2 Cooperative Research Centre for Spatial Information, Docklands, VIC 3008, Australia 3 Department of Geodesy and Surveying, Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece; [email protected] * Correspondence: [email protected] Received: 22 October 2018; Accepted: 23 November 2018; Published: 27 November 2018 Abstract: An analysis of current static and time-variable gravity field models is presented focusing on the medium to high frequencies of the geopotential as expressed by the spherical harmonic coefficients. A validation scheme of the gravity field models is implemented based on dynamic orbit determination that is applied in a degree-wise cumulative sense of the individual spherical harmonics. The approach is applied to real data of the Gravity Field and Steady-State Ocean Circulation (GOCE) and Gravity Recovery and Climate Experiment (GRACE) satellite missions, as well as to GRACE inter-satellite K-band ranging (KBR) data. Since the proposed scheme aims at capturing gravitational discrepancies, we consider a few deterministic empirical parameters in order to avoid absorbing part of the gravity signal that may be included in the monitored orbit residuals. The present contribution aims at a band-limited analysis for identifying characteristic degree ranges and thresholds of the various GRACE- and GOCE-based gravity field models. The degree range 100–180 is investigated based on the degree-wise cumulative approach.
    [Show full text]
  • The Evolution of Earth Gravitational Models Used in Astrodynamics
    JEROME R. VETTER THE EVOLUTION OF EARTH GRAVITATIONAL MODELS USED IN ASTRODYNAMICS Earth gravitational models derived from the earliest ground-based tracking systems used for Sputnik and the Transit Navy Navigation Satellite System have evolved to models that use data from the Joint United States-French Ocean Topography Experiment Satellite (Topex/Poseidon) and the Global Positioning System of satellites. This article summarizes the history of the tracking and instrumentation systems used, discusses the limitations and constraints of these systems, and reviews past and current techniques for estimating gravity and processing large batches of diverse data types. Current models continue to be improved; the latest model improvements and plans for future systems are discussed. Contemporary gravitational models used within the astrodynamics community are described, and their performance is compared numerically. The use of these models for solid Earth geophysics, space geophysics, oceanography, geology, and related Earth science disciplines becomes particularly attractive as the statistical confidence of the models improves and as the models are validated over certain spatial resolutions of the geodetic spectrum. INTRODUCTION Before the development of satellite technology, the Earth orbit. Of these, five were still orbiting the Earth techniques used to observe the Earth's gravitational field when the satellites of the Transit Navy Navigational Sat­ were restricted to terrestrial gravimetry. Measurements of ellite System (NNSS) were launched starting in 1960. The gravity were adequate only over sparse areas of the Sputniks were all launched into near-critical orbit incli­ world. Moreover, because gravity profiles over the nations of about 65°. (The critical inclination is defined oceans were inadequate, the gravity field could not be as that inclination, 1= 63 °26', where gravitational pertur­ meaningfully estimated.
    [Show full text]
  • Draft American National Standard Astrodynamics
    BSR/AIAA S-131-200X Draft American National Standard Astrodynamics – Propagation Specifications, Test Cases, and Recommended Practices Warning This document is not an approved AIAA Standard. It is distributed for review and comment. It is subject to change without notice. Recipients of this draft are invited to submit, with their comments, notification of any relevant patent rights of which they are aware and to provide supporting documentation. Sponsored by American Institute of Aeronautics and Astronautics Approved XX Month 200X American National Standards Institute Abstract This document provides the broad astrodynamics and space operations community with technical standards and lays out recommended approaches to ensure compatibility between organizations. Applicable existing standards and accepted documents are leveraged to make a complete—yet coherent—document. These standards are intended to be used as guidance and recommended practices for astrodynamics applications in Earth orbit where interoperability and consistency of results is a priority. For those users who are purely engaged in research activities, these standards can provide an accepted baseline for innovation. BSR/AIAA S-131-200X LIBRARY OF CONGRESS CATALOGING DATA WILL BE ADDED HERE BY AIAA STAFF Published by American Institute of Aeronautics and Astronautics 1801 Alexander Bell Drive, Reston, VA 20191 Copyright © 200X American Institute of Aeronautics and Astronautics All rights reserved No part of this publication may be reproduced in any form, in an electronic retrieval
    [Show full text]
  • The Joint Gravity Model 3
    Journal of Geophysical Research Accepted for publication, 1996 The Joint Gravity Model 3 B. D. Tapley, M. M. Watkins,1 J. C. Ries, G. W. Davis,2 R. J. Eanes, S. R. Poole, H. J. Rim, B. E. Schutz, and C. K. Shum Center for Space Research, University of Texas at Austin R. S. Nerem,3 F. J. Lerch, and J. A. Marshall Space Geodesy Branch, NASA Goddard Space Flight Center, Greenbelt, Maryland S. M. Klosko, N. K. Pavlis, and R. G. Williamson Hughes STX Corporation, Lanham, Maryland Abstract. An improved Earth geopotential model, complete to spherical harmonic degree and order 70, has been determined by combining the Joint Gravity Model 1 (JGM 1) geopotential coef®cients, and their associated error covariance, with new information from SLR, DORIS, and GPS tracking of TOPEX/Poseidon, laser tracking of LAGEOS 1, LAGEOS 2, and Stella, and additional DORIS tracking of SPOT 2. The resulting ®eld, JGM 3, which has been adopted for the TOPEX/Poseidon altimeter data rerelease, yields improved orbit accuracies as demonstrated by better ®ts to withheld tracking data and substantially reduced geographically correlated orbit error. Methods for analyzing the performance of the gravity ®eld using high-precision tracking station positioning were applied. Geodetic results, including station coordinates and Earth orientation parameters, are signi®cantly improved with the JGM 3 model. Sea surface topography solutions from TOPEX/Poseidon altimetry indicate that the ocean geoid has been improved. Subset solutions performed by withholding either the GPS data or the SLR/DORIS data were computed to demonstrate the effect of these particular data sets on the gravity model used for TOPEX/Poseidon orbit determination.
    [Show full text]
  • Statistical Orbit Determination
    Preface The modem field of orbit determination (OD) originated with Kepler's inter­ pretations of the observations made by Tycho Brahe of the planetary motions. Based on the work of Kepler, Newton was able to establish the mathematical foundation of celestial mechanics. During the ensuing centuries, the efforts to im­ prove the understanding of the motion of celestial bodies and artificial satellites in the modem era have been a major stimulus in areas of mathematics, astronomy, computational methodology and physics. Based on Newton's foundations, early efforts to determine the orbit were focused on a deterministic approach in which a few observations, distributed across the sky during a single arc, were used to find the position and velocity vector of a celestial body at some epoch. This uniquely categorized the orbit. Such problems are deterministic in the sense that they use the same number of independent observations as there are unknowns. With the advent of daily observing programs and the realization that the or­ bits evolve continuously, the foundation of modem precision orbit determination evolved from the attempts to use a large number of observations to determine the orbit. Such problems are over-determined in that they utilize far more observa­ tions than the number required by the deterministic approach. The development of the digital computer in the decade of the 1960s allowed numerical approaches to supplement the essentially analytical basis for describing the satellite motion and allowed a far more rigorous representation of the force models that affect the motion. This book is based on four decades of classroom instmction and graduate- level research.
    [Show full text]
  • The Celestial Mechanics Approach: Theoretical Foundations
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by RERO DOC Digital Library J Geod (2010) 84:605–624 DOI 10.1007/s00190-010-0401-7 ORIGINAL ARTICLE The celestial mechanics approach: theoretical foundations Gerhard Beutler · Adrian Jäggi · Leoš Mervart · Ulrich Meyer Received: 31 October 2009 / Accepted: 29 July 2010 / Published online: 24 August 2010 © Springer-Verlag 2010 Abstract Gravity field determination using the measure- of the CMA, in particular to the GRACE mission, may be ments of Global Positioning receivers onboard low Earth found in Beutler et al. (2010) and Jäggi et al. (2010b). orbiters and inter-satellite measurements in a constellation of The determination of the Earth’s global gravity field using satellites is a generalized orbit determination problem involv- the data of space missions is nowadays either based on ing all satellites of the constellation. The celestial mechanics approach (CMA) is comprehensive in the sense that it encom- 1. the observations of spaceborne Global Positioning (GPS) passes many different methods currently in use, in particular receivers onboard low Earth orbiters (LEOs) (see so-called short-arc methods, reduced-dynamic methods, and Reigber et al. 2004), pure dynamic methods. The method is very flexible because 2. or precise inter-satellite distance monitoring of a close the actual solution type may be selected just prior to the com- satellite constellation using microwave links (combined bination of the satellite-, arc- and technique-specific normal with the measurements of the GPS receivers on all space- equation systems. It is thus possible to generate ensembles crafts involved) (see Tapley et al.
    [Show full text]
  • Non-Averaged Regularized Formulations As an Alternative to Semi-Analytical Orbit Propagation Methods
    Celestial Mechanics and Dynamical Astronomy manuscript No. (will be inserted by the editor) Non-averaged regularized formulations as an alternative to semi-analytical orbit propagation methods Davide Amato ¨ Claudio Bombardelli ¨ Giulio Ba`u ¨ Vincent Morand ¨ Aaron J. Rosengren Received: date / Accepted: date Abstract This paper is concerned with the comparison of semi-analytical and non-averaged propagation methods for Earth satellite orbits. We analyse the total integration error for semi-analytical methods and propose a novel decomposition into dynamical, model truncation, short-periodic, and numerical error com- ponents. The first three are attributable to distinct approximations required by the method of averaging, which fundamentally limit the attainable accuracy. In contrast, numerical error, the only component present in non-averaged methods, can be significantly mitigated by employing adaptive numerical algorithms and regularized formulations of the equations of motion. We present a collection of non-averaged methods based on the integration of existing regularized formulations of the equations of motion through an adaptive solver. We implemented the collection in the orbit propagation code THALASSA, which we make publicly available, and we compared the non-averaged methods to the semi-analytical method implemented in the orbit prop- agation tool STELA through numerical tests involving long-term propagations (on the order of decades) of LEO, GTO, and high-altitude HEO orbits. For the test cases considered, regularized non-averaged methods were found to be up to two times slower than semi-analytical for the LEO orbit, to have comparable speed for the GTO, and to be ten times as fast for the HEO (for the same accuracy).
    [Show full text]
  • Why Atens Enjoy Enhanced Accessibility for Human Space Flight
    (Preprint) AAS 11-449 WHY ATENS ENJOY ENHANCED ACCESSIBILITY FOR HUMAN SPACE FLIGHT Daniel R. Adamo* and Brent W. Barbee† Near-Earth objects can be grouped into multiple orbit classifications, among them being the Aten group, whose members have orbits crossing Earth's with semi-major axes less than 1 astronomical unit. Atens comprise well under 10% of known near-Earth objects. This is in dramatic contrast to results from recent human space flight near-Earth object accessibility studies, where the most favorable known destinations are typically almost 50% Atens. Geocentric dynamics explain this enhanced Aten accessibility and lead to an understanding of where the most accessible near-Earth objects reside. Without a com- prehensive space-based survey, however, highly accessible Atens will remain largely un- known. INTRODUCTION In the context of human space flight (HSF), the concept of near-Earth object (NEO) accessibility is highly subjective (Reference 1). Whether or not a particular NEO is accessible critically depends on mass, performance, and reliability of interplanetary HSF systems yet to be designed. Such systems would cer- tainly include propulsion and crew life support with adequate shielding from both solar flares and galactic cosmic radiation. Equally critical architecture options are relevant to NEO accessibility. These options are also far from being determined and include the number of launches supporting an HSF mission, together with whether consumables are to be pre-emplaced at the destination. Until the unknowns of HSF to NEOs come into clearer focus, the notion of relative accessibility is of great utility. Imagine a group of NEOs, each with nearly equal HSF merit determined from their individual characteristics relating to crew safety, scientific return, resource utilization, and planetary defense.
    [Show full text]
  • Highly-Reduced Dynamic Orbits and Their Use for Global Gravity Field
    HIGHLY-REDUCED DYNAMIC ORBITS AND THEIR USE View metadata, citation and similarFOR papers GLOBAL at core.ac.uk GRAVITY FIELD RECOVERY: brought to you by CORE A SIMULATION STUDY FOR GOCE provided by RERO DOC Digital Library A. JÄGGI1, H. BOCK1, R. PAIL2 AND H. GOIGINGER2 1 University of Bern, Astronomical Institute, Sidlerstrasse 5, 3012 Bern, Switzerland ([email protected], [email protected]) 2 TU Graz, Institute of Navigation and Satellite Geodesy, Steyrergasse 30, 8010 Graz, Austria ([email protected], [email protected]) Received: November 2, 2007; Revised: February 4, 2008; Accepted: March 17, 2008 ABSTRACT The so-called highly reduced-dynamic (HRD) orbit determination strategy and its use for the determination of the Earth’s gravitational field are analyzed. We discuss the functional model for the generation of HRD orbits, which are a compromise of the two extreme cases of dynamic and purely geometrically determined kinematic orbits. For gravity field recovery the energy integral approach is applied, which is based on the law of energy conservation in a closed system. The potential of HRD orbits for gravity field determination is studied in the frame of a simulated test environment based on a realistic GOCE orbit configuration. The results are analyzed, assessed, and compared with the respective reference solutions based on a kinematic orbit scenario. The main advantage of HRD orbits is the fact that they contain orbit velocity information, thus avoiding numerical differentiation on the orbit positions. The error characteristics are usually much smoother, and the computation of gravity field solutions is more efficient, because less densely sampled orbit information is sufficient.
    [Show full text]
  • Improvements in Autonomous GPS Navigation of Low Earth Orbit Satellites
    Ph.D. Thesis Doctoral Program in Aerospace Science & Technology Improvements in autonomous GPS navigation of Low Earth Orbit satellites Pere Ramos-Bosch Advisors: Dr. Manuel Hern´andezPajares [1] Dr. J.Miguel Juan Zornoza [1] Dr. Oliver Montenbruck [2] [1] Research group of Astronomy and Geomatics (gAGE) Depts. of Applied Mathematics IV and Applied Physics Universitat Polit`ecnicade Catalunya (UPC), Spain [2] German Space Operations Center (DLR), Germany September 10, 2008 ii Acknowledgements I would like to thank my advisors, Dr. Manuel Hern´andezPajares, Dr. Miguel Juan Zornoza and Dr. Oliver Montenbruck for their enthusiastic support and useful discussions in the whole development of this thesis. Without their unconditional aid that work would had not been possible. I would like to thank the funding support for both the thesis and a stage at the German Space Operations Center (Deutsches Zentrum f¨urLuft- und Raumfahrt - DLR) from the “Departament d’Universitats, Recerca i Societat de la Informaci´ode la Generalitat de Catalunya” and from the European Social Fund. This work started in the frame of the Neural Networks for Radion- avigation project funded by ESA (grant number ITT4584). I would like to thank the project team for their support, and to the Institute of Navigation for their Student Sponsorship to present the initial version of this work at ION GNSS 2006, held in Fort Worth, USA, in September 2006. Special thanks to Bruce Haines of JPL for kindly providing attitude data for JASON satellite. I would also like to thank all the people in the group of Astronomy and GEomatics (gAGE), it has been a real pleasure to share with you all kind of discussions.
    [Show full text]
  • Kinematic Orbit Determination of Low Earth Orbit Satellites Using Gps and Galileo Observations
    KINEMATIC ORBIT DETERMINATION OF LOW EARTH ORBIT SATELLITES USING GPS AND GALILEO OBSERVATIONS A THESIS SUBMITTED TO THE GRADUATE SCHOOL OF NATURAL AND APPLIED SCIENCES OF MIDDLE EAST TECHNICAL UNIVERSITY BY OZAN KILIÇ IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF MASTER OF SCIENCE IN GEODETIC AND GEOGRAPHIC INFORMATION TECHNOLOGIES JULY 2019 Approval of the thesis: KINEMATIC ORBIT DETERMINATION OF LOW EARTH ORBIT SATELLITES USING GPS AND GALILEO OBSERVATIONS submitted by OZAN KILIÇ in partial fulfillment of the requirements for the degree of Master of Science in Geodetic and Geographic Information Technologies Department, Middle East Technical University by, Prof. Dr. Halil Kalıpçılar Dean, Graduate School of Natural and Applied Sciences Prof. Dr. Zuhal Akyürek Head of Department, Geodetic and Geographic Inf. Tech. Prof. Dr. Mahmut Onur Karslıoglu˘ Supervisor, Geodetic and Geographic Inf. Tech., METU Examining Committee Members: Prof. Dr. Zuhal Akyürek Civil Engineering, METU Prof. Dr. Mahmut Onur Karslıoglu˘ Geodetic and Geographic Information Technologies, METU Prof. Dr. Haluk Özener KOERI, Bogazici University Assoc. Prof. Dr. Ugur˘ Murat Leloglu˘ Geodetic and Geographic Information Technologies, METU Assist. Prof. Dr. Kamil Teke Geomatics Engineering, Hacettepe University Date: I hereby declare that all information in this document has been obtained and presented in accordance with academic rules and ethical conduct. I also declare that, as required by these rules and conduct, I have fully cited and referenced all material and results that are not original to this work. Name, Surname: Ozan Kılıç Signature : iv ABSTRACT KINEMATIC ORBIT DETERMINATION OF LOW EARTH ORBIT SATELLITES USING GPS AND GALILEO OBSERVATIONS Kılıç, Ozan M.S., Department of Geodetic and Geographic Information Technologies Supervisor: Prof.
    [Show full text]