Gravity An Introduction

Similar documents
Fluid Physics 8.292J/12.330J % (1)

CARIBBEAN EXAMINATIONS COUNCIL CARIBBEAN SECONDARY EDUCATION EXAMINATION ADDITIONAL MATHEMATICS. Paper 02 - General Proficiency

Lecture 7: Polar representation of complex numbers

ARISTOTELIAN PHYSICS

Fundamental Concepts: Surfaces and Curves

METHOD OF FUNDAMENTAL SOLUTIONS FOR HELMHOLTZ EIGENVALUE PROBLEMS IN ELLIPTICAL DOMAINS

Problem Cosider the curve give parametrically as x = si t ad y = + cos t for» t» ß: (a) Describe the path this traverses: Where does it start (whe t =

: Transforms and Partial Differential Equations

The Mathematical Model and the Simulation Modelling Algoritm of the Multitiered Mechanical System

(c) Write, but do not evaluate, an integral expression for the volume of the solid generated when R is

lil fit c tai an.ie't 1111 At I 6qei ATA I Atb Y Ex Find the linear regression today power law example forhw i gl.es 6 inner product spaces

AP Calculus BC Review Applications of Derivatives (Chapter 4) and f,

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,

Mathematics Extension 2

Signal Processing in Mechatronics. Lecture 3, Convolution, Fourier Series and Fourier Transform

2C09 Design for seismic and climate changes

The axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c.

mx bx kx F t. dt IR I LI V t, Q LQ RQ V t,

Vibratory Motion. Prof. Zheng-yi Feng NCHU SWC. National CHung Hsing University, Department of Soil and Water Conservation

Signal Processing in Mechatronics

CALCULUS BASIC SUMMER REVIEW

Streamfunction-Vorticity Formulation

Voltage controlled oscillator (VCO)

Section 1.1. Calculus: Areas And Tangents. Difference Equations to Differential Equations

On the Construction of a Synthetic Earth Gravity Model

Finally, we show how to determine the moments of an impulse response based on the example of the dispersion model.

Salmon: Lectures on partial differential equations. 3. First-order linear equations as the limiting case of second-order equations

COMPARISON OF SOME METHODS FOR MODIFYING STOKES' FORMULA IN THE GOCE ERA

Statistical Noise Models and Diagnostics

Notes on the GSW function gsw_geostrophic_velocity (geo_strf,long,lat,p)

Chapter 4. Fourier Series

Mathematics Extension 1

Two or more points can be used to describe a rigid body. This will eliminate the need to define rotational coordinates for the body!

The AMSU Observation Bias Correction and Its Application Retrieval Scheme, and Typhoon Analysis

GOCE. Gravity and steady-state Ocean Circulation Explorer

Castiel, Supernatural, Season 6, Episode 18

Lecture 7: Density Estimation: k-nearest Neighbor and Basis Approach

Quantization and Special Functions

A widely used display of protein shapes is based on the coordinates of the alpha carbons - - C α

GULF MATHEMATICS OLYMPIAD 2014 CLASS : XII

On the Weak Localization Principle of the Eigenfunction Expansions of the Laplace-Beltrami Operator by Riesz Method ABSTRACT 1.

Higher Course Plan. Calculus and Relationships Expressions and Functions

PHYS-3301 Lecture 9. CHAPTER 5 Wave Properties of Matter and Quantum Mechanics I. 5.3: Electron Scattering. Bohr s Quantization Condition

Assignment 1 : Real Numbers, Sequences. for n 1. Show that (x n ) converges. Further, by observing that x n+2 + x n+1

Complex Analysis Spring 2001 Homework I Solution

Problem 1. Problem Engineering Dynamics Problem Set 9--Solution. Find the equation of motion for the system shown with respect to:

φ φ sin φ θ sin sin u = φθu

Synopsis of Euler s paper. E Memoire sur la plus grande equation des planetes. (Memoir on the Maximum value of an Equation of the Planets)

Hilbert Space and Least-squares Collocation

AP Calculus BC 2005 Scoring Guidelines

Chapter 2 Systems and Signals

Lecture 4 Conformal Mapping and Green s Theorem. 1. Let s try to solve the following problem by separation of variables

ES.182A Topic 40 Notes Jeremy Orloff

Statistical Fundamentals and Control Charts

AP Calculus BC 2011 Scoring Guidelines Form B

A Slight Extension of Coherent Integration Loss Due to White Gaussian Phase Noise Mark A. Richards

Free Surface Hydrodynamics

DETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS.

Linear Regression Models

THE KALMAN FILTER RAUL ROJAS

1.3 Convergence Theorems of Fourier Series. k k k k. N N k 1. With this in mind, we state (without proof) the convergence of Fourier series.

8. Applications To Linear Differential Equations

ELEG 4603/5173L Digital Signal Processing Ch. 1 Discrete-Time Signals and Systems

Formula List for College Algebra Sullivan 10 th ed. DO NOT WRITE ON THIS COPY.

*X203/701* X203/701. APPLIED MATHEMATICS ADVANCED HIGHER Numerical Analysis. Read carefully

Mathematics Extension 2

EDEXCEL NATIONAL CERTIFICATE UNIT 4 MATHEMATICS FOR TECHNICIANS OUTCOME 4 - CALCULUS

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

6.003 Homework #3 Solutions

Patterns in Complex Numbers An analytical paper on the roots of a complex numbers and its geometry

AH Checklist (Unit 2) AH Checklist (Unit 2) Proof Theory

6.003 Homework #12 Solutions

567. Research of Dynamics of a Vibration Isolation Platform

Paper-II Chapter- Damped vibration

SLIP TEST 3 Chapter 2,3 and 6. Part A Answer all the questions Each question carries 1 mark 1 x 1 =1.

ROBUST ATTITUDE DETERMINATION AND CONTROL SYSTEM OF MICROSATELLITE. Principle of the control system construction

ECON 3150/4150, Spring term Lecture 3

ECE4270 Fundamentals of DSP. Lecture 2 Discrete-Time Signals and Systems & Difference Equations. Overview of Lecture 2. More Discrete-Time Systems

Solutions to quizzes Math Spring 2007

Math 113 Exam 3 Practice

HE ATOM & APPROXIMATION METHODS MORE GENERAL VARIATIONAL TREATMENT. Examples:

Topic 10: Introduction to Estimation

Phys. 201 Mathematical Physics 1 Dr. Nidal M. Ershaidat Doc. 12

MID-YEAR EXAMINATION 2018 H2 MATHEMATICS 9758/01. Paper 1 JUNE 2018

PARTIAL DIFFERENTIAL EQUATIONS SEPARATION OF VARIABLES

Chapter 10 Partial Differential Equations and Fourier Series

LECTURE 14. Non-linear transverse motion. Non-linear transverse motion

Similarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall

6.003 Homework #12 Solutions

Bayesian Methods: Introduction to Multi-parameter Models

Chapter 10 Advanced Topics in Random Processes

Session 5. (1) Principal component analysis and Karhunen-Loève transformation

[ ] sin ( ) ( ) = 2 2 ( ) ( ) ( ) ˆ Mechanical Spectroscopy II

Limitation of Applicability of Einstein s. Energy-Momentum Relationship

Discrete-Time Signals and Systems. Discrete-Time Signals and Systems. Signal Symmetry. Elementary Discrete-Time Signals.

2. Neutronic calculations at uranium powered cylindrical reactor by using Bessel differential equation

MCT242: Electronic Instrumentation Lecture 2: Instrumentation Definitions

Holographic Renyi Entropy

Physics Supplement to my class. Kinetic Theory

MIDTERM 2 CALCULUS 2. Monday, October 22, 5:15 PM to 6:45 PM. Name PRACTICE EXAM

Transcription:

Gravity A Itroductio Reier Rummel Istitute of Advaced Study (IAS) Techische Uiversität Müche Lecture Oe 5th ESA Earth Observatio Summer School 2-13 August 2010, ESA-ESRIN, Frascati / Italy

gravitatio ad space sciece micro-g eviromet: e.g. boilig water fudametal physics: e.g. equivalece priciple earth scieces e.g. gravitatioal field of earth moo ad plaets

global map of gravity aomalies

the pla Lecture Oe: Theoretical basics of gravitatio as applied to the earth [gravitatioal law, properties, mathematical represetatio] The laguage for the two other lectures Lecture Two: The role of the earth s gravitatioal field i earth scieces [gravity aomalies, geoid as a referece, temporal variatios] Lecture Three: Priciples of satellite gravimetry; i their logic derived from free fall tests i a laboratory o earth [the orbit, priciples of GRACE, ESA s missio GOCE ad satellite gradiometry]

itroductio to gravitatio Newto s law of gravitatio: mm mm F G e G x x A B A B A 2 AB 3 A B AB x x A B Newto s secod law: F A m I A a A F A A gravitatioal acceleratio: m m aa G e m A I A B 2 AB AB AB e AB B

itroductio to gravitatio from sigle mass, to may masses, to a cotiuum a G e d B A 2 AB B AB Fudametal properties of Newto s law of gravitatio: cetral force actio = reactio iverse square distace superpositio of all partial forces istataeous

itroductio to gravitatio a G e d B A 2 AB B AB a A is a vector field i space with the followig properties: a 0 curl free a V A A A i.e. there exists a gravitatioal potetial V A ad i outer space, we get: a source free V 2 0 0

itroductio to gravitatio a G e d B A 2 AB B AB a A is a vector field with the followig properties: a 0 curl free a V A A A i.e. there exists a gravitatioal potetial V A ad i outer space, we get: a source free V 2 0 0 B VA G d AB B V V V x y z 2 2 2 0 2 2 2

example: satellite orbit 2 V 2 x a ice applicatio: a satellite orbit x A aa AV " perturbatios " ad iitial coditios : x ; x 0 0

example: satellite orbit a ice applicatio: a satellite orbit sphere flatteed sphere real earth Kelpleria ellipse precessig ellipse precessig ellipse plus gravitatioal code

example: satellite orbit What about the attractio of su, moo ad plaets? Aswer: They determie the earth s orbit about the su Tides are acceleratio relative to the earth s ceter of mass Marshak, 2005

gravitatio ad gravity o the surface of the rotatig earth oe measures: gravity= gravitatio + cetrifugal acceleratio g a z ad W V Z z a g

size of gravity sigals gravity (i laboratory at TU Müche) 9.807 246 72 m/s 2 statioary variable 10 0 spherical Earth 10-3 flatteig & cetrifugal acceleratio 10-4 moutais, valleys, ocea ridges, subductio 10-5 desity variatios i crust ad matle 10-6 salt domes, sedimet basis, ores 10-7 tides, atmospheric pressure 10-8 temporal variatios: oceas, hydrology 10-9 ocea topography, polar motio 10-10 geeral relativity

geometry of the earth s gravity field gravity vector plumb lie geoid W W cost 0. level surface

gravity related quatities ) <1m to 2m

gravity related quatities map with geoid heights relative to the GRS80 ellipsoid 80.00 60.00 50.00 40.00 20.00 0.00 0.00-20.00-40.00-50.00-60.00-80.00-150.00-100.00-50.00 0.00 50.00 100.00 150.00-100.00

series represetatio of gravitatioal field V V V x y z 2 2 2 0 2 2 2 Laplace equatio (PDE) solutio i Cartesia coordiates (after determiatio):,, 2 2 0 exp k exp VA x y z V k z c i kx y with z 0 k solutio i spherical coordiates (after determiatio): 1 R VA,, r V0 Pm Cm cosm Sm sim 0 r m0 1 1 R R V0 C P cos m S P si m m m m m m0 m r m r almost a Fourier series

series represetatio of gravitatioal field log 500km 333km 250km 200km 20000km [ km] max spectral represetatio of the earth s gravitatioal field: triagular plot of the spherical harmoic coefficiets C ; S m m

series represetatio of gravitatioal field C 00 A 00 C 11 S 11 C 10 S 22 S 21 C 20 C 21 C 22 S 33 S 32 S 31 C 30 C 31 C 32 C 33 S 44 S 43 S 42 S 41 C 40 C 41 C 42 C 43 C 44 S 55 S 54 S 53 S 52 S 51 C 50 C 51 C 52 C 53 C 54 C 55 sectorial tesseral zoal tesseral sectorial surface spherical harmoic fuctios: Ym, Pm cos m sim

series represetatio of gravitatioal field Degree C 00 S 11 C 10 C 11 S 22 S 21 C 21 scale origi orietatio S m Order m C m

series represetatio of gravitatioal field sigal degree variace 2 2 ( m m) m0 c C S ave c 21 i aalogy to sigal processig: characteristics of sigal ad oise Here: degree variaces correspod to power spectral desity

series represetatio of gravitatioal field power law by WM Kaula error c 1.610 / 10 3 2 sigal ad error degree variaces (dimesioless)

series represetatio of temporal variatios of gravitatioal field 10 20 30 40 50 10-3 10-3 Degree Stadard Deviatio i Geoid Height [m] Atmosphere daily ECMWF Ocea daily MIT Hydrology mothly Europe GRACE Error Predictio 10-4 10-4 10-5 10-5 10-6 10 20 30 40 10-6 50 Degree rapid time variable geoid sigals (RMS) [ to be divided by the earth radius i order to arrive at dimesioless uits]

series represetatio of temporal variatios of gravitatioal field 10 20 30 40 50 10-3 10-3 Degree Stadard Deviatio i Geoid Height [m] Ocea semi-aual Ocea aual GRACE Error Predictio Atmosphere aual 10-4 10-4 10-5 10-5 10-6 10-6 10 20 30 40 50 Degree slow geoid time variable sigals (RMS) [ to be divided by the earth radius i order to arrive at dimesioless uits]

three levels of gravity quatities o earth ad i space 1 R VA,, r V0 Pm Cm cosm Sm sim 0 r m0 1 1 R R V C P 0 cos m S P si m m m m m m0 m r m r r 1 r 2 δv δ r V δ rr V disturbace potetial or geoid gravity disturbaces or gravity aomalies gravity gradiets or torsio balace

various gravity quatities o earth ad i space Gravity model EGM08, D/O1000

h = 0km δv δv r δv rr Gravity model EGM08, D/O1000

satellite altitude: r earth s Surface: R δv δ r V δ rr V 1 R 2 R δv δ r V δ rr V 1 r 2 r 1 R r 2 R r 3 R r three levels of gravity quatities o earth ad i space

various gravity quatities o earth ad i space Gravity model EGM08, D/O1000

δv r h = 400km h = 250km h = 0km Gravity model EGM08, D/O1000

T r T rr T 1 R 2 R T r T rr T 1 r 2 r 1 R r 2 R r 3 R r Gravity model EGM08, D/O1000

various gravity quatities o earth ad i space x y z x y V ik [E] 0.5 0 0.5 z

satellite-to-satellite trackig low-low satellite-to-satellite trackig high-low satellite gradiometry satellite altitude r r 1 r 2 δv δ r V δ rr V R r 1 R r 2 R r 3 earth s surface R R 1 R 2 δv δ r V δ rr V disturbace potetial or geoid gravity disturbaces or gravity aomalies gravity gradiets or torsio balace

summary of lecture Oe 1. Newto s law of gravitatio describes all its relevat properties such as iverse square distace, priciple of superpositio, its statioary part beig vorticity free, ad source free outside the earth (Laplace equatio) 2. Gravity is the sum of gravitatio ad the cetrifugal part 3. Satellite orbits are essetially described by gravitatio 4. Tides are a accelertatio (a force) relative to the earth s ceter of mass 5. The global gravitatioal field is represeted as a series of spherical harmoics, beig a solutio of Laplace partial differetial equatio (Dirichlet) 6. Spherical harmoics o a sphere are aalogous to a Fourier series i a plae 7. Therefore there exists a closed theory of sigal ad oise processig 8. With icreasig distace from the earth sphere the series coefficiets are dampeig out per degree like (R/R+h) +1 9. With each radial derivative of the gravitatioal potetial the series coefficiets are amplified per degree like (+1) 10. The strategy of satellite missios GRACE ad GOCE rests o the priciple of compesatig the damperig effect by amplificatio (see 8. ad 9.)