Appendix to Section 3: a Brief Overview of Geodetics

Total Page:16

File Type:pdf, Size:1020Kb

Appendix to Section 3: a Brief Overview of Geodetics Appendix to Section 3: A brief overview of Geodetics MAE 5540 - Propulsion Systems Geodesy • Navigation Geeks do Calculations in Geocentric (spherical) Coordinates • Map Makers Give Surface Data in Terms !earth of Geodetic (elliptical) Coordinates • Need to have some idea how to relate one to another -- science of geodesy MAE 5540 - Propulsion Systems How Does the Earth Radius Vary with Latitude? k z ? 2 ? 2 REarth = x + y2 + z = r + z r R " Ellipse: y j ! r 2 + z 2 = 1 j 2 Req Req 1 - eEarth y x i "a" "b" r Prime Meridian ! x MAE 5540 - Propulsion Systems How Does the Earth Radius vary with Latitude? r 2 + z 2 = 1 2 Req Req 1 - eEarth ! 2 2 2 2 2 1 - eEarth r + z = Req 1 - eEarth ! z 2 2 2 2 z 2 r Req = r + 2 = r 1 + 2 1 - eEarth 1 - eEarth Rearth z " r 2 2 2 z 2 2 r = Rearth cos " r = tan " MAE 5540 - Propulsion Systems How Does the Earth Radius vary with Latitude? 2 2 Req 2 tan ! 2 = cos ! 1 + 2 = Rearth 1 - eEarth 2 2 2 1 - eEarth cos ! + sin ! 2 = 1 - eEarth Inverting .... 2 2 2 2 2 2 cos ! + sin ! - eEarth cos ! 1 - eEarth cos ! 2 = 2 1 - eEarth 1 - eEarth R earth ! = 2 Req 1 - eEarth 2 2 1 - eEarth cos ! MAE 5540 - Propulsion Systems Earth Radius vs Geocentric Latitude R earth ! = Polar Radius: 6356.75170 km 2 Req 1 - eEarth Equatorial Radius: 6378.13649 km 2 2 1 - eEarth cos ! 2 2 2 e = 1 - b = a - b = Earth a a2 6378.13649 2 6378.136492 - = 0.08181939 6378.13649 MAE 5540 - Propulsion Systems Earth Radius vs Geocentric Latitude (concluded) 6380. 6375. 6370. Radius, Km 6365. 6360. 6355. -100 -50 0 50 100 Geocentric Latitude, deg. MAE 5540 - Propulsion Systems Earth Radius … alternate formula • Earth radius as Function of Latitude MAE 5540 - Propulsion Systems What is the mean radius of the earth? RE = Req 2 ! 2 dA = " [ RE cos ( !) ] 1 - eEarth ! 2 2 1 - eEarth cos ! RE cos(!) • IAU Convention: ! Based on Earth's Volume RE d! cos (!) ! RE ! Sphere Volume: ! 3 4 " R = V 3 Emean E 2 RE cos(!) x RE d! cos (!) dV = " [ ] ! MAE 5540 - Propulsion Systems What is the Earth's Mean Radius? (continued) • Earth's (ellipsoid) Volume ! 2 V = R cos 3d = E ! E" " " -! 2 ! 2 1 - e2 3/2 R3 Earth cos3 d = eq ! 2 2 " " 1 - eEarth cos " -! 2 3 4 ! 2 R 1- e eq 3 earth MAE 5540 - Propulsion Systems What is the Earth's Mean Radius? (continued) • Based on Volume 3 Ellipsoid Volume: 4 ! 1- e2 R 3 eq 3 Sphere Volume: 4 ! R 3 sphere " R R = 1- 0.081819392 1/6 6378.13649 = 6371.0002 km sphere # mean • Mean Radius we have been using is for a Sphere with same volume as the Earth ME = $E VE "gravitational radius" MAE 5540 - Propulsion Systems Earth Radius vs Geocentric Latitude (concluded) “Gaussian Surface” 6380. 6375. 6370. Radius, Km 6365. 6360. 6355. -100 -50 0 50 100 Geocentric Latitude, deg. MAE 5540 - Propulsion Systems Geocentric vs Geodetic Coordinates Mean North (Celestial) Pole • Map makers Geocentric vs Geodetic C2 ooordinates define a new Req 1 - e(earth) latitude which is Mean Greenwich the angle that Meridian RE normal to the ! Earth's surface ! makes with the Geocenter R respect to the !' eq equatorial plane " Req • Geodetic latitude Equator MAE 5540 - Propulsion Systems Geocentric vs Geodetic Coordinates (contined) Mean •Since the Earth is North (Celestial) Elliptical only along the Pole z-axis ... geodetic and geocentric 2 h longitude are Req 1 - e(earth) identical Mean Greenwich Meridian RE • Altitude is an ! extension of the Geocenter ! line of latitude Req geodetic) !' " Req Equator MAE 5540 - Propulsion Systems Geocentric vs Geodetic Coordinates (contined) •Complex nonlinear Mean equations describe North (Celestial) relationship between Pole geocentric and geodetic latitude 2 h Req 1 - e(earth) Mean • Derivation requires Greenwich Meridian RE Extensive ! Knowledge of Spherical ! Trigonometry Geocenter R !' eq " Req Equator MAE 5540 - Propulsion Systems Geocentric vs Geodetic Coordinates (contined) Geocentric Cartesian Coordinates * rtarget xtarget = R' +h cos !' cos " !' h y = R' +h cos ' sin z target !' ! " 2 !' ztarget = R' 1 - e +h sin !' ! !' earth R' ' Req ! R' = !' 2 1 - e sin2 !' !' earth rtargetr * We would be here all week if I try to derive this "Radius of Curvature" Req MWAhEit m55o4re0, S- tPeprohepnu Als.i, oand S Hyasetreinmgs, Edward A., Jr., FORTRAN Program for the Analysis of Ground Based Range Tracking Data--Usage and Derivations, NASA TM 104201, December, 1992 Geocentric vs Geodetic Coordinates (contined) Geocentric Polar Coordinates r 2 2 2 target R = x + y + z target target target target h z -1 ytarget target = tan " xtarget !' ! z -1 target R' !target = tan 2 2 !' xtarget + ytarget !' rtargetr "Radius of Curvature" Req MAE 5540 - Propulsion Systems Geocentric vs Geodetic Coordinates (concluded) Inverse Relationships, non-linear no direct solution 2 2 xtarget + ytarget h = - R' cos !' !' y = tan-1 target "target xtarget ztarget 1 ' = -1 # ! target tan 2 2 R' xtarget + ytarget 1 - e2 !' earth R' + h !' target MWAhEit m55o4re0, S- tPeprohepnu Als.i, oand S Hyasetreinmgs, Edward A., Jr., FORTRAN Program for the Analysis of Ground Based Range Tracking Data--Usage and Derivations, NASA TM 104201, December, 1992 Pulling it all together • Given geodetic coordinates -- compute geocentric i) Compute geocedtric cartesian coordinates Range, runway threshholds, radar antennae, beacon .... x = R' +h cos !' cos " target !' y = R' +h cos !' sin " target !' 2 ztarget = R' 1 - e +h sin !' !' earth Req R' = !' 2 1 - e sin2 !' MAE 5540 - Propulsion Systems earth Pulling it all together (continued) • Given geodetic coordinates -- compute geocentric ii) Compute Geocentric polar coordinates next 2 2 2 Rtarget = xtarget + ytarget + ztarget -1 ytarget target = tan ! xtarget -1 ztarget "target = tan 2 + 2 xtarget ytarget MAE 5540 - Propulsion Systems Pulling it all together (concluded) • Given geocentric (usually x,y,z) coordinates -- compute geodetic Req R' = ' 2 ! 1 - e sin2 !' GPS, INS, TLE's .... earth No explicit solution: 2 2 xtarget + ytarget requires h = - R' 1) series expansion solution, cos !' !' 2) numerical iteration, -1 ytarget 3) or a special solution "target = tan called "Ferrari's method"** xtarget z -1 target 1 !'target = tan # x 2 + y 2 R' target target 1 - e2 !' earth R' + h !' target **NASA Technical Paper 3430, Whitmore and Haering, FORTRAN Program forAnalyzing Ground-Based Tracking Data: Usage and Derivations, Version 6.2, 1995 MAE 5540 - Propulsion Systems Numerical Example • Edwards Air Force Base, Radar Site #34 Numerical example !' = 34.96081° " = –117.91150° h = 2563.200 ft • Find corresponding geocentric cartesian and polar coordinates MAE 5540 - Propulsion Systems Numerical Example (cont'd) • Compute Local Radius of Curvature Req R' = = !' 2 1 - e sin2 !' earth 6378.13649 km = 1 - 0.08181939 sin 34.96081 " # 2 180 6392.187109 km MAE 5540 - Propulsion Systems Numerical Example (cont'd) • Compute X and Y (geocentric) r = R' +h cos ' = target !' ! 6392.1871 + 2536.2 " 3. 048 " 10-4 cos 34.96081 " # = 180 5239.3131 km xtarget = rtarget cos $ = ytarget = rtarget sin $ = 5239.3131 km " -117.91150 " # = 5216.0074 km " sin -117.91150 " # = cos 180 180 -2452.5602 km -4629.83218 km MAE 5540 - Propulsion Systems Numerical Example (cont'd) • Compute z (geocentric) 2 ztarget = R' 1 - e +h sin !' = !' earth 6392.1871 km 0.081819392 2536.2 " 3. 048 " 10-4 sin 34.96081 " # = 1 - + 180 3638.7480 km MAE 5540 - Propulsion Systems Numerical Example (cont'd) • Compute Geocentric Polar Coordinates 2 2 2 Rtarget = 2452.5602 + 4629.83218 + 3638.7480 = 6378.94104km -1 ytarget = tan = !target xtarget 180 # -1 -4629.83218 = -117.9115° " tan -2452.5602 z -1 target = $target = tan 2 2 xtarget + ytarget 180 # -1 3638.7480 = 34.7803° " tan 5239.3131 MAE 5540 - Propulsion Systems Numerical Example (cont'd) • Compute Local Earth Radius and Geocentric Distance Above Geoid RE = Req = ! 2 1 - eEarth 2 2 1 - eEarth cos ! 6378.13649 km = 1 - 0.081819392 0.081819392 2 34.7083 " # 1 - cos 180 6364.23 km D = R - R = geoid target E! 6378.94104 - 6364.23 km = 14.7 km = 48228.25 ft MAE 5540 - Propulsion Systems Numerical Example (concluded) • Comparison Geodetic Geocentric !' = 34.96081° ! = 34.7803° " = –117.91150° " = –117.91150° h = 2563.200 ft Dgeoid = 48228.25 ft • Earth Oblateness is NOT trivial, and in the REAL World -- it must be accounted for MAE 5540 - Propulsion Systems Appendix II: Rigorous Derivation of Realizable Launch Inclination • Launch Initial Conditions North Pole Greenwich Meridian • Position: ", Latitude z #, Longitude (Inertial) Launch meridian h, Altitude Equator y Re " x Vernal Equinox !g 95 MAE 5540 - Propulsion Systems • Launch Initial Conditions North Pole Greenwich Meridian z Launch meridian • Launch Initial Conditions (Inertial Coordinates) Equator y Re " x Vernal Equinox !g Sidereal hour angle 96 MAE 5540 - Propulsion Systems Sidereal Hour Angle Ignore for now! 97 MAE 5540 - Propulsion Systems Computing the Hour Angle 98 MAE 5540 - Propulsion Systems.
Recommended publications
  • Navigation and the Global Positioning System (GPS): the Global Positioning System
    Navigation and the Global Positioning System (GPS): The Global Positioning System: Few changes of great importance to economics and safety have had more immediate impact and less fanfare than GPS. • GPS has quietly changed everything about how we locate objects and people on the Earth. • GPS is (almost) the final step toward solving one of the great conundrums of human history. Where the heck are we, anyway? In the Beginning……. • To understand what GPS has meant to navigation it is necessary to go back to the beginning. • A quick look at a ‘precision’ map of the world in the 18th century tells one a lot about how accurate our navigation was. Tools of the Trade: • Navigations early tools could only crudely estimate location. • A Sextant (or equivalent) can measure the elevation of something above the horizon. This gives your Latitude. • The Compass could provide you with a measurement of your direction, which combined with distance could tell you Longitude. • For distance…well counting steps (or wheel rotations) was the thing. An Early Triumph: • Using just his feet and a shadow, Eratosthenes determined the diameter of the Earth. • In doing so, he used the last and most elusive of our navigational tools. Time Plenty of Weaknesses: The early tools had many levels of uncertainty that were cumulative in producing poor maps of the world. • Step or wheel rotation counting has an obvious built-in uncertainty. • The Sextant gives latitude, but also requires knowledge of the Earth’s radius to determine the distance between locations. • The compass relies on the assumption that the North Magnetic pole is coincident with North Rotational pole (it’s not!) and that it is a perfect dipole (nope…).
    [Show full text]
  • Chapter 11. Three Dimensional Analytic Geometry and Vectors
    Chapter 11. Three dimensional analytic geometry and vectors. Section 11.5 Quadric surfaces. Curves in R2 : x2 y2 ellipse + =1 a2 b2 x2 y2 hyperbola − =1 a2 b2 parabola y = ax2 or x = by2 A quadric surface is the graph of a second degree equation in three variables. The most general such equation is Ax2 + By2 + Cz2 + Dxy + Exz + F yz + Gx + Hy + Iz + J =0, where A, B, C, ..., J are constants. By translation and rotation the equation can be brought into one of two standard forms Ax2 + By2 + Cz2 + J =0 or Ax2 + By2 + Iz =0 In order to sketch the graph of a quadric surface, it is useful to determine the curves of intersection of the surface with planes parallel to the coordinate planes. These curves are called traces of the surface. Ellipsoids The quadric surface with equation x2 y2 z2 + + =1 a2 b2 c2 is called an ellipsoid because all of its traces are ellipses. 2 1 x y 3 2 1 z ±1 ±2 ±3 ±1 ±2 The six intercepts of the ellipsoid are (±a, 0, 0), (0, ±b, 0), and (0, 0, ±c) and the ellipsoid lies in the box |x| ≤ a, |y| ≤ b, |z| ≤ c Since the ellipsoid involves only even powers of x, y, and z, the ellipsoid is symmetric with respect to each coordinate plane. Example 1. Find the traces of the surface 4x2 +9y2 + 36z2 = 36 1 in the planes x = k, y = k, and z = k. Identify the surface and sketch it. Hyperboloids Hyperboloid of one sheet. The quadric surface with equations x2 y2 z2 1.
    [Show full text]
  • Radar Back-Scattering from Non-Spherical Scatterers
    REPORT OF INVESTIGATION NO. 28 STATE OF ILLINOIS WILLIAM G. STRATION, Governor DEPARTMENT OF REGISTRATION AND EDUCATION VERA M. BINKS, Director RADAR BACK-SCATTERING FROM NON-SPHERICAL SCATTERERS PART 1 CROSS-SECTIONS OF CONDUCTING PROLATES AND SPHEROIDAL FUNCTIONS PART 11 CROSS-SECTIONS FROM NON-SPHERICAL RAINDROPS BY Prem N. Mathur and Eugam A. Mueller STATE WATER SURVEY DIVISION A. M. BUSWELL, Chief URBANA. ILLINOIS (Printed by authority of State of Illinois} REPORT OF INVESTIGATION NO. 28 1955 STATE OF ILLINOIS WILLIAM G. STRATTON, Governor DEPARTMENT OF REGISTRATION AND EDUCATION VERA M. BINKS, Director RADAR BACK-SCATTERING FROM NON-SPHERICAL SCATTERERS PART 1 CROSS-SECTIONS OF CONDUCTING PROLATES AND SPHEROIDAL FUNCTIONS PART 11 CROSS-SECTIONS FROM NON-SPHERICAL RAINDROPS BY Prem N. Mathur and Eugene A. Mueller STATE WATER SURVEY DIVISION A. M. BUSWELL, Chief URBANA, ILLINOIS (Printed by authority of State of Illinois) Definitions of Terms Part I semi minor axis of spheroid semi major axis of spheroid wavelength of incident field = a measure of size of particle prolate spheroidal coordinates = eccentricity of ellipse = angular spheroidal functions = Legendse polynomials = expansion coefficients = radial spheroidal functions = spherical Bessel functions = electric field vector = magnetic field vector = back scattering cross section = geometric back scattering cross section Part II = semi minor axis of spheroid = semi major axis of spheroid = Poynting vector = measure of size of spheroid = wavelength of the radiation = back scattering
    [Show full text]
  • John Ellipsoid 5.1 John Ellipsoid
    CSE 599: Interplay between Convex Optimization and Geometry Winter 2018 Lecture 5: John Ellipsoid Lecturer: Yin Tat Lee Disclaimer: Please tell me any mistake you noticed. The algorithm at the end is unpublished. Feel free to contact me for collaboration. 5.1 John Ellipsoid In the last lecture, we discussed that any convex set is very close to anp ellipsoid in probabilistic sense. More precisely, after renormalization by covariance matrix, we have kxk2 = n ± Θ(1) with high probability. In this lecture, we will talk about how convex set is close to an ellipsoid in a strict sense. If the convex set is isotropic, it is close to a sphere as follows: Theorem 5.1.1. Let K be a convex body in Rn in isotropic position. Then, rn + 1 B ⊆ K ⊆ pn(n + 1)B : n n n Roughly speaking, this says that any convex set can be approximated by an ellipsoid by a n factor. This result has a lot of applications. Although the bound is tight, making a body isotropic is pretty time- consuming. In fact, making a body isotropic is the current bottleneck for obtaining faster algorithm for sampling in convex sets. Currently, it can only be done in O∗(n4) membership oracle plus O∗(n5) total time. Problem 5.1.2. Find a faster algorithm to approximate the covariance matrix of a convex set. In this lecture, we consider another popular position of a convex set called John position and its correspond- ing ellipsoid is called John ellipsoid. Definition 5.1.3. Given a convex set K.
    [Show full text]
  • The Global Positioning System the Global Positioning System
    The Global Positioning System The Global Positioning System 1. System Overview 2. Biases and Errors 3. Signal Structure and Observables 4. Absolute v. Relative Positioning 5. GPS Field Procedures 6. Ellipsoids, Datums and Coordinate Systems 7. Mission Planning I. System Overview ! GPS is a passive navigation and positioning system available worldwide 24 hours a day in all weather conditions developed and maintained by the Department of Defense ! The Global Positioning System consists of three segments: ! Space Segment ! Control Segment ! User Segment Space Segment Space Segment ! The current GPS constellation consists of 29 Block II/IIA/IIR/IIR-M satellites. The first Block II satellite was launched in February 1989. Control Segment User Segment How it Works II. Biases and Errors Biases GPS Error Sources • Satellite Dependent ? – Orbit representation ? Satellite Orbit Error Satellite Clock Error including 12 biases ? 9 3 Selective Availability 6 – Satellite clock model biases Ionospheric refraction • Station Dependent L2 L1 – Receiver clock biases – Station Coordinates Tropospheric Delay • Observation Multi- pathing Dependent – Ionospheric delay 12 9 3 – Tropospheric delay Receiver Clock Error 6 1000 – Carrier phase ambiguity Satellite Biases ! The satellite is not where the GPS broadcast message says it is. ! The satellite clocks are not perfectly synchronized with GPS time. Station Biases ! Receiver clock time differs from satellite clock time. ! Uncertainties in the coordinates of the station. ! Time transfer and orbital tracking. Observation Dependent Biases ! Those associated with signal propagation Errors ! Residual Biases ! Cycle Slips ! Multipath ! Antenna Phase Center Movement ! Random Observation Error Errors ! In addition to biases factors effecting position and/or time determined by GPS is dependant upon: ! The geometric strength of the satellite configuration being observed (DOP).
    [Show full text]
  • Download/Pdf/Ps-Is-Qzss/Ps-Qzss-001.Pdf (Accessed on 15 June 2021)
    remote sensing Article Design and Performance Analysis of BDS-3 Integrity Concept Cheng Liu 1, Yueling Cao 2, Gong Zhang 3, Weiguang Gao 1,*, Ying Chen 1, Jun Lu 1, Chonghua Liu 4, Haitao Zhao 4 and Fang Li 5 1 Beijing Institute of Tracking and Telecommunication Technology, Beijing 100094, China; [email protected] (C.L.); [email protected] (Y.C.); [email protected] (J.L.) 2 Shanghai Astronomical Observatory, Chinese Academy of Sciences, Shanghai 200030, China; [email protected] 3 Institute of Telecommunication and Navigation, CAST, Beijing 100094, China; [email protected] 4 Beijing Institute of Spacecraft System Engineering, Beijing 100094, China; [email protected] (C.L.); [email protected] (H.Z.) 5 National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100094, China; [email protected] * Correspondence: [email protected] Abstract: Compared to the BeiDou regional navigation satellite system (BDS-2), the BeiDou global navigation satellite system (BDS-3) carried out a brand new integrity concept design and construction work, which defines and achieves the integrity functions for major civil open services (OS) signals such as B1C, B2a, and B1I. The integrity definition and calculation method of BDS-3 are introduced. The fault tree model for satellite signal-in-space (SIS) is used, to decompose and obtain the integrity risk bottom events. In response to the weakness in the space and ground segments of the system, a variety of integrity monitoring measures have been taken. On this basis, the design values for the new B1C/B2a signal and the original B1I signal are proposed, which are 0.9 × 10−5 and 0.8 × 10−5, respectively.
    [Show full text]
  • Introduction to Aberrations OPTI 518 Lecture 13
    Introduction to aberrations OPTI 518 Lecture 13 Prof. Jose Sasian OPTI 518 Topics • Aspheric surfaces • Stop shifting • Field curve concept Prof. Jose Sasian OPTI 518 Aspheric Surfaces • Meaning not spherical • Conic surfaces: Sphere, prolate ellipsoid, hyperboloid, paraboloid, oblate ellipsoid or spheroid • Cartesian Ovals • Polynomial surfaces • Infinite possibilities for an aspheric surface • Ray tracing for quadric surfaces uses closed formulas; for other surfaces iterative algorithms are used Prof. Jose Sasian OPTI 518 Aspheric surfaces The concept of the sag of a surface 2 cS 4 6 8 10 ZS ASASASAS4 6 8 10 ... 1 1 K ( c2 1 ) 2 S Sxy222 K 2 K is the conic constant K=0, sphere K=-1, parabola C is 1/r where r is the radius of curvature; K is the K<-1, hyperola conic constant (the eccentricity squared); -1<K<0, prolate ellipsoid A’s are aspheric coefficients K>0, oblate ellipsoid Prof. Jose Sasian OPTI 518 Conic surfaces focal properties • Focal points for Ellipsoid case mirrors Hyperboloid case Oblate ellipsoid • Focal points of lenses K n2 Hecht-Zajac Optics Prof. Jose Sasian OPTI 518 Refraction at a spherical surface Aspheric surface description 2 cS 4 6 8 10 ZS ASASASAS4 6 8 10 ... 1 1 K ( c2 1 ) 2 S 1122 2222 22 y x Aicrasphe y 1 x 4 K y Zx 28rr Prof. Jose Sasian OTI 518P Cartesian Ovals ln l'n' Cte . Prof. Jose Sasian OPTI 518 Aspheric cap Aspheric surface The aspheric surface can be thought of as comprising a base sphere and an aspheric cap Cap Spherical base surface Prof.
    [Show full text]
  • Geodynamics and Rate of Volcanism on Massive Earth-Like Planets
    The Astrophysical Journal, 700:1732–1749, 2009 August 1 doi:10.1088/0004-637X/700/2/1732 C 2009. The American Astronomical Society. All rights reserved. Printed in the U.S.A. GEODYNAMICS AND RATE OF VOLCANISM ON MASSIVE EARTH-LIKE PLANETS E. S. Kite1,3, M. Manga1,3, and E. Gaidos2 1 Department of Earth and Planetary Science, University of California at Berkeley, Berkeley, CA 94720, USA; [email protected] 2 Department of Geology and Geophysics, University of Hawaii at Manoa, Honolulu, HI 96822, USA Received 2008 September 12; accepted 2009 May 29; published 2009 July 16 ABSTRACT We provide estimates of volcanism versus time for planets with Earth-like composition and masses 0.25–25 M⊕, as a step toward predicting atmospheric mass on extrasolar rocky planets. Volcanism requires melting of the silicate mantle. We use a thermal evolution model, calibrated against Earth, in combination with standard melting models, to explore the dependence of convection-driven decompression mantle melting on planet mass. We show that (1) volcanism is likely to proceed on massive planets with plate tectonics over the main-sequence lifetime of the parent star; (2) crustal thickness (and melting rate normalized to planet mass) is weakly dependent on planet mass; (3) stagnant lid planets live fast (they have higher rates of melting than their plate tectonic counterparts early in their thermal evolution), but die young (melting shuts down after a few Gyr); (4) plate tectonics may not operate on high-mass planets because of the production of buoyant crust which is difficult to subduct; and (5) melting is necessary but insufficient for efficient volcanic degassing—volatiles partition into the earliest, deepest melts, which may be denser than the residue and sink to the base of the mantle on young, massive planets.
    [Show full text]
  • Reference Systems for Surveying and Mapping Lecture Notes
    Delft University of Technology Reference Systems for Surveying and Mapping Lecture notes Hans van der Marel ii The front cover shows the NAP (Amsterdam Ordnance Datum) ”datum point” at the Stopera, Amsterdam (picture M.M.Minderhoud, Wikipedia/Michiel1972). H. van der Marel Lecture notes on Reference Systems for Surveying and Mapping: CTB3310 Surveying and Mapping CTB3425 Monitoring and Stability of Dikes and Embankments CIE4606 Geodesy and Remote Sensing CIE4614 Land Surveying and Civil Infrastructure February 2020 Publisher: Faculty of Civil Engineering and Geosciences Delft University of Technology P.O. Box 5048 Stevinweg 1 2628 CN Delft The Netherlands Copyright ©2014­2020 by H. van der Marel The content in these lecture notes, except for material credited to third parties, is licensed under a Creative Commons Attributions­NonCommercial­SharedAlike 4.0 International License (CC BY­NC­SA). Third party material is shared under its own license and attribution. The text has been type set using the MikTex 2.9 implementation of LATEX. Graphs and diagrams were produced, if not mentioned otherwise, with Matlab and Inkscape. Preface This reader on reference systems for surveying and mapping has been initially compiled for the course Surveying and Mapping (CTB3310) in the 3rd year of the BSc­program for Civil Engineering. The reader is aimed at students at the end of their BSc program or at the start of their MSc program, and is used in several courses at Delft University of Technology. With the advent of the Global Positioning System (GPS) technology in mobile (smart) phones and other navigational devices almost anyone, anywhere on Earth, and at any time, can determine a three–dimensional position accurate to a few meters.
    [Show full text]
  • Geodetic Position Computations
    GEODETIC POSITION COMPUTATIONS E. J. KRAKIWSKY D. B. THOMSON February 1974 TECHNICALLECTURE NOTES REPORT NO.NO. 21739 PREFACE In order to make our extensive series of lecture notes more readily available, we have scanned the old master copies and produced electronic versions in Portable Document Format. The quality of the images varies depending on the quality of the originals. The images have not been converted to searchable text. GEODETIC POSITION COMPUTATIONS E.J. Krakiwsky D.B. Thomson Department of Geodesy and Geomatics Engineering University of New Brunswick P.O. Box 4400 Fredericton. N .B. Canada E3B5A3 February 197 4 Latest Reprinting December 1995 PREFACE The purpose of these notes is to give the theory and use of some methods of computing the geodetic positions of points on a reference ellipsoid and on the terrain. Justification for the first three sections o{ these lecture notes, which are concerned with the classical problem of "cCDputation of geodetic positions on the surface of an ellipsoid" is not easy to come by. It can onl.y be stated that the attempt has been to produce a self contained package , cont8.i.ning the complete development of same representative methods that exist in the literature. The last section is an introduction to three dimensional computation methods , and is offered as an alternative to the classical approach. Several problems, and their respective solutions, are presented. The approach t~en herein is to perform complete derivations, thus stqing awrq f'rcm the practice of giving a list of for11111lae to use in the solution of' a problem.
    [Show full text]
  • The Global Positioning System II Field Experiments
    The Global Positioning System II Field Experiments Geo327G/386G: GIS & GPS Applications in Earth Sciences Jackson School of Geosciences, University of Texas at Austin 10/8/2020 5-1 Mexico DGPS Field Campaign Cenotes in Tamaulipas, MX, near Aldama Geo327G/386G: GIS & GPS Applications in Earth Sciences 10/8/2020 Jackson School of Geosciences, University of Texas at Austin 5-2 Are Cenote Water Levels Related? Geo327G/386G: GIS & GPS Applications in Earth Sciences 10/8/2020 Jackson School of Geosciences, University of Texas at Austin 5-3 DGPS Static Survey of Cenote Water Levels Geo327G/386G: GIS & GPS Applications in Earth Sciences 10/8/2020 Jackson School of Geosciences, University of Texas at Austin 5-4 Determining Orthometric Heights ❑Ortho. Height = H.A.E. – Geoid Height Height above MSL (Orthometric height) H.A.E. Geoid height Earth Surface = H.A.E. Geoid Ellipsoid Geoid Height Geo327G/386G: GIS & GPS Applications in Earth Sciences 10/8/2020 Jackson School of Geosciences, University of Texas at Austin 5-5 Determining Orthometric Heights ❑Ortho. Height = (H.A.E. – Geoid Height) Need: 1) Ellipsoid model – GRS80 – NAVD88 ❑reference stations: HARN (+ 2 cm), CORS (+ ~2 cm) 2) Geoid model – GEOID99 ( + 5 cm for US) Procedure: Base receiver at reference station, rover at point of interest a) measure HAE, apply DGS corrections b) subtract local Geoid Height Geo327G/386G: GIS & GPS Applications in Earth Sciences 10/8/2020 Jackson School of Geosciences, University of Texas at Austin 5-6 Sources of Error ❑Geoid error – model less well constrained in areas of few gravity measurement ❑NAVD88 error – benchmark stability, measurement errors ❑GPS errors – need precise ephemeri, tropospheric delay model, equipment (antennae should be same for base and rover) Geo327G/386G: GIS & GPS Applications in Earth Sciences 10/8/2020 Jackson School of Geosciences, University of Texas at Austin 5-7 Static Carrier-phase solutions obtained by: ❑Commercial post-processing software ❑e.g.
    [Show full text]
  • Gravity, Geoid and Earth Observation IAG Commission 2: Gravity Field, Chania, Crete, Greece, 23-27 June 2008
    S.P. Mertikas (Ed.) Gravity, Geoid and Earth Observation IAG Commission 2: Gravity Field, Chania, Crete, Greece, 23-27 June 2008 Series: International Association of Geodesy Symposia, Vol. 135 ▶ State of the art scientific achievements of gravity field research prospects These Proceedings include the written version of papers presented at the IAG International Symposium on "Gravity, Geoid and Earth Observation 2008". The Symposium was held in Chania, Crete, Greece, 23-27 June 2008 and organized by the Laboratory of Geodesy and Geomatics Engineering, Technical University of Crete, Greece. The meeting was arranged by the International Association of Geodesy and in particular by the IAG Commission 2: Gravity Field. The symposium aimed at bringing together geodesists and geophysicists working in the 2010, XXXIV, 702 p. 340 illus. general areas of gravity, geoid, geodynamics and Earth observation. Besides covering the traditional research areas, special attention was paid to the use of geodetic methods for: Earth observation, environmental monitoring, Global Geodetic Observing System (GGOS), Printed book Earth Gravity Models (e.g., EGM08), geodynamics studies, dedicated gravity satellite Hardcover missions (i.e., GOCE), airborne gravity surveys, Geodesy and geodynamics in polar regions, and the integration of geodetic and geophysical information. ▶ 299,99 € | £249.99 | $379.99 ▶ *320,99 € (D) | 329,99 € (A) | CHF 354.00 eBook Available from your bookstore or ▶ springer.com/shop MyCopy Printed eBook for just ▶ € | $ 24.99 ▶ springer.com/mycopy Order online at springer.com ▶ or for the Americas call (toll free) 1-800-SPRINGER ▶ or email us at: [email protected]. ▶ For outside the Americas call +49 (0) 6221-345-4301 ▶ or email us at: [email protected].
    [Show full text]