Earth Science Applications of Space Based Geodesy DES7355
![Earth Science Applications of Space Based Geodesy DES-7355 Tu-Th 9: 40 -11: 05 Seminar Earth Science Applications of Space Based Geodesy DES-7355 Tu-Th 9: 40 -11: 05 Seminar](https://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-1.jpg)
Earth Science Applications of Space Based Geodesy DES-7355 Tu-Th 9: 40 -11: 05 Seminar Room in 3892 Central Ave. (Long building) Bob Smalley Office: 3892 Central Ave, Room 103 678 -4929 Office Hours – Wed 14: 00 -16: 00 or if I’m in my office. http: //www. ceri. memphis. edu/people/smalley/ESCI 7355/ESCI_7355_Applications_of_Space_Based_Geodesy. html Class 2 1
![Coordinate systems Simple spherical Geodetic – with respect to ellipsoid normal to surface does Coordinate systems Simple spherical Geodetic – with respect to ellipsoid normal to surface does](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-2.jpg)
Coordinate systems Simple spherical Geodetic – with respect to ellipsoid normal to surface does not intersect origin [in general] ECEF XYZ – earth centered, earth fixed xyz. Is what it says. 2
![GEODETIC COORDINATES: LATITUDE (Herring) 3 GEODETIC COORDINATES: LATITUDE (Herring) 3](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-3.jpg)
GEODETIC COORDINATES: LATITUDE (Herring) 3
![LONGITUDE Longitude measured by time difference of astronomical events (Herring) 4 LONGITUDE Longitude measured by time difference of astronomical events (Herring) 4](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-4.jpg)
LONGITUDE Longitude measured by time difference of astronomical events (Herring) 4
![The “problem” arises because we’re defining the “location” (latitude) based on the orientation of The “problem” arises because we’re defining the “location” (latitude) based on the orientation of](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-5.jpg)
The “problem” arises because we’re defining the “location” (latitude) based on the orientation of the “surface” of the earth (not the ellipsoid) at the point where we want to determine the location. (Assuming vertical perpendicular to “level”. ) “shape” of the surface of the earth - with the variations greatly exaggerated. For now we’re not being very specific about what the surface represents/how it is defined. Image from: http: //kartoweb. itc. nl/geometrics/Reference%20 surfaces/refsurf. html 5
![This means that we have to take the “shape” of the surface into account This means that we have to take the “shape” of the surface into account](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-6.jpg)
This means that we have to take the “shape” of the surface into account in defining our reference frame. We are still not even considering the vertical. We’re still only discussing the problem of 2 -D location on the surface of the earth. Image from: http: //kartoweb. itc. nl/geometrics/Reference%20 surfaces/refsurf. html 6
![Traditional approach was to define local/regional datums (flattening, size, origin – typically not earth Traditional approach was to define local/regional datums (flattening, size, origin – typically not earth](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-7.jpg)
Traditional approach was to define local/regional datums (flattening, size, origin – typically not earth centered, orientation). (Assume gravity perpendicular to "surface" – which is not really the case - since measurements made with a level. ) 7
![These datums were “best fits” for the regions that they covered. They could be These datums were “best fits” for the regions that they covered. They could be](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-8.jpg)
These datums were “best fits” for the regions that they covered. They could be quite bad (up to 1 km error) outside those regions however. 8
![These datums are also not “earth centered” (origin not center of mass of earth). These datums are also not “earth centered” (origin not center of mass of earth).](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-9.jpg)
These datums are also not “earth centered” (origin not center of mass of earth). Converting from one to another not trivial in practice. 9
![Real shape can also have uniqueness problem – using level - more than one Real shape can also have uniqueness problem – using level - more than one](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-10.jpg)
Real shape can also have uniqueness problem – using level - more than one spot with same “latitude”! 10
![“Modern” solution is an earth centered global “best fit” ellipsoid to the shape of “Modern” solution is an earth centered global “best fit” ellipsoid to the shape of](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-11.jpg)
“Modern” solution is an earth centered global “best fit” ellipsoid to the shape of the earth – the geoid. Here we introduce the “thing” that defines the “shape” of the earth – the GEOID. The geoid is what defines the local vertical (and where gravity sneaks in). Image from: http: //kartoweb. itc. nl/geometrics/Reference%20 surfaces/refsurf. html 11
![The geoid is a physical thing (like the topographic surface)– an equipotential of the The geoid is a physical thing (like the topographic surface)– an equipotential of the](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-12.jpg)
The geoid is a physical thing (like the topographic surface)– an equipotential of the gravity field. But we may not be able to “locate” it (can’t see it like we can the topographic surface). So we have to make a “model” for the geoid. Image from: http: //kartoweb. itc. nl/geometrics/Reference%20 surfaces/refsurf. html 12
![Here we introduce the concept of “physical” vs “geometric” position. The geoid (since it Here we introduce the concept of “physical” vs “geometric” position. The geoid (since it](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-13.jpg)
Here we introduce the concept of “physical” vs “geometric” position. The geoid (since it depends on the actual “shape “of the earth, and we will see that it directly effects traditional measurements of latitude) gives a physical definition of position. 13
![Here we introduce the concept of “physical” vs “geometric” position. The ellipsoid gives a Here we introduce the concept of “physical” vs “geometric” position. The ellipsoid gives a](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-14.jpg)
Here we introduce the concept of “physical” vs “geometric” position. The ellipsoid gives a geometric definition of position (and we will see that “modern” positioning – GPS for example – works in this system – even though gravity and other physics effects the system). 14
![Here we introduce the concept of “physical” vs “geometric” position. The horizontal “datum” is Here we introduce the concept of “physical” vs “geometric” position. The horizontal “datum” is](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-15.jpg)
Here we introduce the concept of “physical” vs “geometric” position. The horizontal “datum” is a best fit ellipsoid (to a region or the whole earth) to the shape (geoid) used as a coordinate system for specifying horizontal position. 15
![WHAT ABOUT HEIGHT Geocentric coordinates (f, l, h) (this is based on standard spherical WHAT ABOUT HEIGHT Geocentric coordinates (f, l, h) (this is based on standard spherical](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-16.jpg)
WHAT ABOUT HEIGHT Geocentric coordinates (f, l, h) (this is based on standard spherical coordinate system with h=R-Re, height is clearly defined, simple to understand). h From Kelso, Orbital Coordinate Systems, Part I, Satellite Times, Sep/Oct 1995 16
![WHAT ABOUT HEIGHT For the Ellipsoid coordinates (f, l, h) – Ellipsoidal/Geodetic height. Distance WHAT ABOUT HEIGHT For the Ellipsoid coordinates (f, l, h) – Ellipsoidal/Geodetic height. Distance](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-17.jpg)
WHAT ABOUT HEIGHT For the Ellipsoid coordinates (f, l, h) – Ellipsoidal/Geodetic height. Distance of a point from the ellipsoid measured along the perpendicular from the ellipsoid to this point. h From Kelso, Orbital Coordinate Systems, Part III, Satellite Times, Jan/Feb 1996 17
![What about HEIGHT For the Geoid things get a little more interesting. The height What about HEIGHT For the Geoid things get a little more interesting. The height](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-18.jpg)
What about HEIGHT For the Geoid things get a little more interesting. The height is the distance of a point from the geoid measured along the perpendicular from the geoid (direction from gravity) to this point. Notice that – the height above the geoid (red line) may not be/is not the same as the ellipsoid height (blue line) and that height above the geoid may not be unique Image from: http: //kartoweb. itc. nl/geometrics/Reference%20 surfaces/refsurf. html 18
![What about HEIGHT when we use a level to find the vertical (traditional surveying) What about HEIGHT when we use a level to find the vertical (traditional surveying)](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-19.jpg)
What about HEIGHT when we use a level to find the vertical (traditional surveying) we are measuring with respect to the geoid (what is the “geoid”? ). Image from: http: //kartoweb. itc. nl/geometrics/Reference%20 surfaces/refsurf. html 19
![This brings us to a fundamental problem in Geodesy ---‘"Height" is a common, ordinary This brings us to a fundamental problem in Geodesy ---‘"Height" is a common, ordinary](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-20.jpg)
This brings us to a fundamental problem in Geodesy ---‘"Height" is a common, ordinary everyday word and everyone knows what it means. Or, more likely, everyone has an idea of what it means, but nailing down an exact definition is surprisingly tricky. ‘ 20
![The geoid is the “actual” shape of the earth. Where the word “actual” is The geoid is the “actual” shape of the earth. Where the word “actual” is](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-21.jpg)
The geoid is the “actual” shape of the earth. Where the word “actual” is in quotes for a reason! 21
![The geoid is a representation of the surface the earth would have if the The geoid is a representation of the surface the earth would have if the](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-22.jpg)
The geoid is a representation of the surface the earth would have if the sea covered the earth. This is not the surface one would get if one poured more water on the earth until there is no more dry land! It is the shape a fluid Earth (of the correct volume) would have if that fluid Earth had exactly the same gravity field as the actual Earth. Where did this reference to the gravity field sneak in? 22
![Since water is a fluid, it cannot support shear stresses. This means that the Since water is a fluid, it cannot support shear stresses. This means that the](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-23.jpg)
Since water is a fluid, it cannot support shear stresses. This means that the surface of the sea (or of a lake, or of water in a bucket, etc. ) will be -- perpendicular to the force of gravity -- an equipotential surface (or else it will flow until the surface of the body of water is everywhere in this state). So the definition of the “shape” of the earth, the geoid, is intimately and inseparably tied to the 23 earth’s gravity field.
![This is good gravity is one of the most well understood branches of Physics. This is good gravity is one of the most well understood branches of Physics.](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-24.jpg)
This is good gravity is one of the most well understood branches of Physics. This is bad the gravity field of the earth depends on the details of the mass distribution within the earth (which do not depend on the first principles of physics – the mass distribution of the earth is as we find it!). 24
![The geoid is a representation of the surface the earth would have if a The geoid is a representation of the surface the earth would have if a](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-25.jpg)
The geoid is a representation of the surface the earth would have if a sea covered an earth with the same gravity field. It is the shape a fluid Earth would have if it had exactly the same gravity field as the actual Earth. 25
![The definition is clear concise, and well defined physically. Problems arise when trying to The definition is clear concise, and well defined physically. Problems arise when trying to](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-26.jpg)
The definition is clear concise, and well defined physically. Problems arise when trying to find where this surface actually physically resides due to things like -- currents, winds, tides effecting “sea level” -- where is this imaginary surface located on land? (generally below the land surface – except where the land surface is below sea level, e. g. Death Valley, Dead Sea - it is the level of fluid in channels cut through the land [approximately]. ) 26
![So – what does this surface – the geoid – actually look like? (greatly So – what does this surface – the geoid – actually look like? (greatly](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-27.jpg)
So – what does this surface – the geoid – actually look like? (greatly exaggerated in the vertical) 27
![Shaded, color coded “topographic” representation of the geoid Valleys Hills 28 Shaded, color coded “topographic” representation of the geoid Valleys Hills 28](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-28.jpg)
Shaded, color coded “topographic” representation of the geoid Valleys Hills 28
![Bad joke for the day "What's up? " "Perpendicular to the geoid. " 29 Bad joke for the day "What's up? " "Perpendicular to the geoid. " 29](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-29.jpg)
Bad joke for the day "What's up? " "Perpendicular to the geoid. " 29
![2. Geodesy Shape of the earth / gravity, geoid (physical) reference frames, ellipsoids (geometric) 2. Geodesy Shape of the earth / gravity, geoid (physical) reference frames, ellipsoids (geometric)](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-30.jpg)
2. Geodesy Shape of the earth / gravity, geoid (physical) reference frames, ellipsoids (geometric) From Mulcare or http: //www. ordnancesurvey. co. uk/oswebsite/gps/information/coordinatesystemsinfo/guidecontents/guide 2. html 30
![2. Geodesy How gravity makes it “interesting” Which way is “up”? (how does water 2. Geodesy How gravity makes it “interesting” Which way is “up”? (how does water](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-31.jpg)
2. Geodesy How gravity makes it “interesting” Which way is “up”? (how does water flow? ) What about measurements with light? From Mulcare 31
![What is the Geoid? Since the geoid is a complicated physical entity that is What is the Geoid? Since the geoid is a complicated physical entity that is](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-32.jpg)
What is the Geoid? Since the geoid is a complicated physical entity that is practically indescrible – Find a “best fit” ellipsoid (and look at variations with respect to this ellipsoid). Current NGS definition The equipotential surface of the Earth’s gravity field which best fits, in a least squares sense, global mean sea level. From Mulcare 32
![And now following the axiom that “one person’s noise is another person’s signal” -- And now following the axiom that “one person’s noise is another person’s signal” --](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-33.jpg)
And now following the axiom that “one person’s noise is another person’s signal” -- Geodesy uses gravity to define the geoid (which we will later see is the reference for traditional forms of measuring height). -- Geophysics uses gravity variations, known as anomalies, to learn about density variations in the interior of the earth to interpret figure 33
![One can (some people do) make a career of modeling the “actual “geoid by One can (some people do) make a career of modeling the “actual “geoid by](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-34.jpg)
One can (some people do) make a career of modeling the “actual “geoid by using spherical harmonic expansions of the geoid with respect to the ellipsoidal best fit geoid. There are ~40, 000 terms in the “best” expansions. Famous “pear” shape of earth. 34
![Geodetic Reference Surfaces A beachball globe Mathematical best fit to Earth’s surface: used for Geodetic Reference Surfaces A beachball globe Mathematical best fit to Earth’s surface: used for](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-35.jpg)
Geodetic Reference Surfaces A beachball globe Mathematical best fit to Earth’s surface: used for defining Latitude and Longitude Modeled best fit to “sea surface“ equipotential gravity field used for defining Elevation The real deal Fig from NGS: file: ///C: /Documents%20 and%20 Settings/Bob/My%20 Documents/geodesy/noaa/geo 03_figure. html 35
![Heights and Vertical Datums Define location by triplet - (latitude, longitude, height) hp 36 Heights and Vertical Datums Define location by triplet - (latitude, longitude, height) hp 36](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-36.jpg)
Heights and Vertical Datums Define location by triplet - (latitude, longitude, height) hp 36
![Heights and Vertical Datums More precisely - Geodetic latitude and longitude – referred to Heights and Vertical Datums More precisely - Geodetic latitude and longitude – referred to](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-37.jpg)
Heights and Vertical Datums More precisely - Geodetic latitude and longitude – referred to oblate ellipsoid. Height referred to perpendicular to oblate ellipsoid. (geometrical, is “accessible” by GPS for example). hp This is called “ellipsoidal” height, hp 37
![In traditional surveying – height is measured with respect to mean sea level, which In traditional surveying – height is measured with respect to mean sea level, which](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-38.jpg)
In traditional surveying – height is measured with respect to mean sea level, which serves as the vertical datum (and is accessible at the origin point). Height is measured as distance along the “plumb” line (which is not actually straight) and is called orthometric height (Hp) 38 Jekeli, 2002: http: //www. fgg. uni-lj. si/~/mkuhar/Zalozba/Heights_Jekeli. pdf
![Line follows gradient of level surfaces. www. evergladesplan. org/pm/recover/ recover_docs/mrt/ft_lauderdale. ppt (are not parallel) Line follows gradient of level surfaces. www. evergladesplan. org/pm/recover/ recover_docs/mrt/ft_lauderdale. ppt (are not parallel)](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-39.jpg)
Line follows gradient of level surfaces. www. evergladesplan. org/pm/recover/ recover_docs/mrt/ft_lauderdale. ppt (are not parallel) Little problem – geoid defined by equipotential surface, can’t measure where this is on continents (sometimes even have problems in oceans), can only 39 measure direction of perpendicular to this surface and
![Ellipsoid, Geoid, and Orthometric Heights “h = H + N” P s th’ce r Ellipsoid, Geoid, and Orthometric Heights “h = H + N” P s th’ce r](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-40.jpg)
Ellipsoid, Geoid, and Orthometric Heights “h = H + N” P s th’ce r a E rfa Su Plumb Line Ellipsoid h Q N Mean Sea Level Ocean “Geoid” PO h (Ellipsoid Height) = Distance along ellipsoid normal (Q to P) N (Geoid Height) = Distance along ellipsoid normal (Q to PO) H (Orthometric Height) = Distance along plumb line (PO to P) David B. Zilkoski 138. 23. 217. 17/jwilbur/student_files/ Spatial%20 Reference%20 Seminar/dzilkoski. ppt 40
![Two questions – 1 Given density distribution, can we calculate the gravitational field? Yes Two questions – 1 Given density distribution, can we calculate the gravitational field? Yes](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-41.jpg)
Two questions – 1 Given density distribution, can we calculate the gravitational field? Yes – Newton’s law of universal gravitation 2 Given volume V, bounded by a surface S, and some information about gravity on S, can you find gravity inside V (where V may or may not contain mass)? Qualified yes (need g or normal gradient to potential everywhere on surface) 41
![Potential Fields As was mentioned earlier, the geoid/mean sea level is defined with respect Potential Fields As was mentioned earlier, the geoid/mean sea level is defined with respect](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-42.jpg)
Potential Fields As was mentioned earlier, the geoid/mean sea level is defined with respect to an equipotential surface. So how do we connect what we need (the equipotential surface) with what we have/can measure (direction and magnitude of the force of gravity) Use potential field theory So, first what are Fields? A field is a function of space and/or time. 42
![Examples of scalar fields temperature topography Contours F(x, y)=const Surface plot (“drawing”) Grey (color) Examples of scalar fields temperature topography Contours F(x, y)=const Surface plot (“drawing”) Grey (color)](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-43.jpg)
Examples of scalar fields temperature topography Contours F(x, y)=const Surface plot (“drawing”) Grey (color) scale J vogt -- http: //www. faculty. iu-bremen. de/jvogt/edu/spring 03/Nat. Sci. Lab 2 -Geo. Astro/nslga 2 -lecture 2. pdf 43
![Examples of vector fields streamlines slopes Vector map J vogt -- http: //www. faculty. Examples of vector fields streamlines slopes Vector map J vogt -- http: //www. faculty.](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-44.jpg)
Examples of vector fields streamlines slopes Vector map J vogt -- http: //www. faculty. iu-bremen. de/jvogt/edu/spring 03/Nat. Sci. Lab 2 -Geo. Astro/nslga 2 -lecture 2. pdf 44
![Examples of vector fields streamlines slopes Plot streamlines J vogt -- http: //www. faculty. Examples of vector fields streamlines slopes Plot streamlines J vogt -- http: //www. faculty.](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-45.jpg)
Examples of vector fields streamlines slopes Plot streamlines J vogt -- http: //www. faculty. iu-bremen. de/jvogt/edu/spring 03/Nat. Sci. Lab 2 -Geo. Astro/nslga 2 -lecture 2. pdf 45
![We are interested in Force fields describe forces acting at each point of space We are interested in Force fields describe forces acting at each point of space](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-46.jpg)
We are interested in Force fields describe forces acting at each point of space at a given time Examples: gravity field magnetic field electrostatic field Fields can be scalar, vector or tensor 46
![We know that work is the product of a force applied through a distance. We know that work is the product of a force applied through a distance.](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-47.jpg)
We know that work is the product of a force applied through a distance. If the work done is independent of the path taken from x 0 to x 1, the work done depends only on the starting and ending positions. WBlue=W 0 WRed=W 0 A force with this type of special property is said to be a “conservative” force. WBlack, 2 step=W 0 47
![If we move around in a conservative force field and return to the starting If we move around in a conservative force field and return to the starting](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-48.jpg)
If we move around in a conservative force field and return to the starting point – by using the blue path to go from A to B and then return to A using the red path for example – the work is zero. WBlue=W 0 We can write this as WRed=W 0 WBlack, 2 step=W 0 48
![Important implication of conservative force field A conservative force field is the derivative (gradient Important implication of conservative force field A conservative force field is the derivative (gradient](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-49.jpg)
Important implication of conservative force field A conservative force field is the derivative (gradient in 3 -D) of a scalar field (function)! 49
![This means our work integral is the solution to the differential equation Where we This means our work integral is the solution to the differential equation Where we](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-50.jpg)
This means our work integral is the solution to the differential equation Where we can define a scalar “potential” function U(x) that is a function of position only and Where we have now included an arbitrary constant of integration. The potential function, U(x), is only defined within a constant – this means we can put the position where U(x)=0 where we want. It also makes it hard to determine its “absolute”, as 50 opposed to “relative” value.
![So now we have the pair of equations If you know U(x), you can So now we have the pair of equations If you know U(x), you can](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-51.jpg)
So now we have the pair of equations If you know U(x), you can compute g(x), where I have changed the letter "F" force to “g” for gravity. If you know the force g(x) and that it is conservative, then you can computer U(x) - to within a constant. 51
![U(x) is potential, the negative of the work done to get to that point. U(x) is potential, the negative of the work done to get to that point.](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-52.jpg)
U(x) is potential, the negative of the work done to get to that point. 52
![So to put this to use we now have to --1) Show that gravity So to put this to use we now have to --1) Show that gravity](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-53.jpg)
So to put this to use we now have to --1) Show that gravity is a conservative force and therefore has an associated potential energy function. 2) Determine the gravity potential and gravity force fields for the earth (first approximation – spherical next approximation – ellipsoidal shape due to rotation and then adjust for rotation) 3) Compare with real earth 53
![Newton’s Universal Law of Gravitation 54 Newton’s Universal Law of Gravitation 54](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-54.jpg)
Newton’s Universal Law of Gravitation 54
![In geophysics one of the masses is usually the earth so Figure from Ahern, In geophysics one of the masses is usually the earth so Figure from Ahern,](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-55.jpg)
In geophysics one of the masses is usually the earth so Figure from Ahern, http: //geophysics. ou. edu/gravmag/potential/gravity_potential. html#newton 55
![Figure from Ahern, http: //geophysics. ou. edu/gravmag/potential/gravity_potential. html#newton 56 Figure from Ahern, http: //geophysics. ou. edu/gravmag/potential/gravity_potential. html#newton 56](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-56.jpg)
Figure from Ahern, http: //geophysics. ou. edu/gravmag/potential/gravity_potential. html#newton 56
![Now we can define the potential as the work done to bring a unit Now we can define the potential as the work done to bring a unit](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-57.jpg)
Now we can define the potential as the work done to bring a unit mass from infinity to a distance r (set the work at infinity to zero) 57
![So we can write the force field as the derivative of a scalar potential So we can write the force field as the derivative of a scalar potential](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-58.jpg)
So we can write the force field as the derivative of a scalar potential field in 1 -D going to 3 -D, it becomes a vector equation and we have Which in spherical coordinates is Figure from Ahern, http: //geophysics. ou. edu/gravmag/potential/gravity_potential. html#newton 58
![Apply to our expression for the gravity potential Which agrees with what we know Apply to our expression for the gravity potential Which agrees with what we know](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-59.jpg)
Apply to our expression for the gravity potential Which agrees with what we know Figure from Ahern, http: //geophysics. ou. edu/gravmag/potential/gravity_potential. html#newton 59
![To find the total potential of gravity we have to integrate over all the To find the total potential of gravity we have to integrate over all the](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-60.jpg)
To find the total potential of gravity we have to integrate over all the point masses in a volume. Figure from Ahern, http: //geophysics. ou. edu/gravmag/potential/gravity_potential. html#newton 60
![To find the total potential of gravity we have to integrate over all the To find the total potential of gravity we have to integrate over all the](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-61.jpg)
To find the total potential of gravity we have to integrate over all the point masses in a volume. Figure from Ahern, http: //geophysics. ou. edu/gravmag/potential/gravity_potential. html#newton 61
![If things are spherically symmetric it is easier to work in spherical coordinates Ex: If things are spherically symmetric it is easier to work in spherical coordinates Ex:](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-62.jpg)
If things are spherically symmetric it is easier to work in spherical coordinates Ex: uniform density sphere Figures from: right - Ahern, http: //geophysics. ou. edu/gravmag/potential/gravity_potential. html#newton, left - http: //www. siu. edu/~cafs/surface/file 13. html 62
![Grinding thorugh 63 Grinding thorugh 63](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-63.jpg)
Grinding thorugh 63
![Grinding thorugh 64 Grinding thorugh 64](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-64.jpg)
Grinding thorugh 64
![So for a uniform density sphere The potential and force of gravity at a So for a uniform density sphere The potential and force of gravity at a](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-65.jpg)
So for a uniform density sphere The potential and force of gravity at a point P, a distance s≥R from the center of the sphere, are Figure after Ahern, http: //geophysics. ou. edu/gravmag/potential/gravity_potential. html#newton 65
![Note that in seismology the vector displacement field solution for P waves is also Note that in seismology the vector displacement field solution for P waves is also](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-66.jpg)
Note that in seismology the vector displacement field solution for P waves is also curl free. This means it is the gradient of a scalar field – call it the P wave potential. So one can work with a scalar wave equation for P waves, which is easier than a vector wave equation, and take the gradient at the end to get the physical P wave displacement vector field. (This is how it is presented in many introductory Seismology books such as Stein and Wysession. ) Unfortunately, unlike with gravity, there is no physical interpretation of the P wave potential 66
![Next ex: Force of gravity from spherical shell 67 After Halliday and Resnick, Fundamentals Next ex: Force of gravity from spherical shell 67 After Halliday and Resnick, Fundamentals](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-67.jpg)
Next ex: Force of gravity from spherical shell 67 After Halliday and Resnick, Fundamentals of Physics
![Force of gravity from spherical shell Uniformly dense spherical shell attracts external mass as Force of gravity from spherical shell Uniformly dense spherical shell attracts external mass as](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-68.jpg)
Force of gravity from spherical shell Uniformly dense spherical shell attracts external mass as if all its mass were concentrated at its 68 After Halliday and Resnick, Fundamentals of Physics
![From inside a shell, the lower limit of integration changes to r-R and we From inside a shell, the lower limit of integration changes to r-R and we](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-69.jpg)
From inside a shell, the lower limit of integration changes to r-R and we get zero. R R+r r r-R 69 After Halliday and Resnick, Fundamentals of Physics
![For a solid sphere – we can make it up of concentric shells. Each For a solid sphere – we can make it up of concentric shells. Each](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-70.jpg)
For a solid sphere – we can make it up of concentric shells. Each shell has to have a uniform density, although different shells can have different densities (density a function of radius only – think “earth”). From outside – we can consider all the mass to be concentrated at the center. 70
![Now we need to find the potential and force for our ellispsoid of revolution Now we need to find the potential and force for our ellispsoid of revolution](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-71.jpg)
Now we need to find the potential and force for our ellispsoid of revolution (a nearly spherical body). (note that we are not starting from scratch with a spinning, self gravitating fluid body and figuring out its equilibrium shape – we’re going to find the gravitational potential and force for an almost, but not quite spherical body. ) Discussion after Turcotte, Ahern and Nerem 71
![EARTH’S GRAVITY FIELD 72 EARTH’S GRAVITY FIELD 72](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-72.jpg)
EARTH’S GRAVITY FIELD 72
![Calculate the potential at a point P (outside) due to a nearly spherical body Calculate the potential at a point P (outside) due to a nearly spherical body](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-73.jpg)
Calculate the potential at a point P (outside) due to a nearly spherical body (the earth). Set up the geometry for the problem: For simplicity - put the origin at the center of mass of the body and let P be on an axis. Discussion after Nerem , Turcotte, and Ahern 73
![Calculate the potential at a point P due to a nearly spherical body. 74 Calculate the potential at a point P due to a nearly spherical body. 74](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-74.jpg)
Calculate the potential at a point P due to a nearly spherical body. 74
![Calculate the potential at a point P due to a nearly spherical body. 75 Calculate the potential at a point P due to a nearly spherical body. 75](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-75.jpg)
Calculate the potential at a point P due to a nearly spherical body. 75
![Calculate the potential at a point P due to a nearly spherical body. 76 Calculate the potential at a point P due to a nearly spherical body. 76](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-76.jpg)
Calculate the potential at a point P due to a nearly spherical body. 76
![Calculate the potential at a point P due to a nearly spherical body. 77 Calculate the potential at a point P due to a nearly spherical body. 77](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-77.jpg)
Calculate the potential at a point P due to a nearly spherical body. 77
![Calculate the potential at a point P due to a nearly spherical body. 78 Calculate the potential at a point P due to a nearly spherical body. 78](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-78.jpg)
Calculate the potential at a point P due to a nearly spherical body. 78
![Calculate the potential at a point P due to a nearly spherical body. 79 Calculate the potential at a point P due to a nearly spherical body. 79](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-79.jpg)
Calculate the potential at a point P due to a nearly spherical body. 79
![Calculate the potential at a point P due to a nearly spherical body. 80 Calculate the potential at a point P due to a nearly spherical body. 80](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-80.jpg)
Calculate the potential at a point P due to a nearly spherical body. 80
![Calculate the potential at a point P due to a nearly spherical body. 81 Calculate the potential at a point P due to a nearly spherical body. 81](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-81.jpg)
Calculate the potential at a point P due to a nearly spherical body. 81
![Calculate the potential at a point P due to a nearly spherical body. 82 Calculate the potential at a point P due to a nearly spherical body. 82](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-82.jpg)
Calculate the potential at a point P due to a nearly spherical body. 82
![Calculate the potential at a point P due to a nearly spherical body. 83 Calculate the potential at a point P due to a nearly spherical body. 83](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-83.jpg)
Calculate the potential at a point P due to a nearly spherical body. 83
![Calculate the potential at a point P due to a nearly spherical body. 84 Calculate the potential at a point P due to a nearly spherical body. 84](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-84.jpg)
Calculate the potential at a point P due to a nearly spherical body. 84
![Calculate the potential at a point P due to a nearly spherical body. Potential Calculate the potential at a point P due to a nearly spherical body. Potential](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-85.jpg)
Calculate the potential at a point P due to a nearly spherical body. Potential for sphere plus adjustments for principal moments of inertia and moment of inertia along axis from origin to point of interest, P. This is Mac. Cullagh’s formula for the potential of a nearly spherical body 85
![Calculate the potential at a point P due to a nearly spherical body. For Calculate the potential at a point P due to a nearly spherical body. For](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-86.jpg)
Calculate the potential at a point P due to a nearly spherical body. For a sphere I 1=I 2=I 3=Iop and (which we knew already) 86
![Calculate the potential at a point P due to a nearly spherical body. So Calculate the potential at a point P due to a nearly spherical body. So](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-87.jpg)
Calculate the potential at a point P due to a nearly spherical body. So here’s our semi-final result for the potential of an approximately spherical body Now let’s look at a particular approximately spherical body – the ellipsoid 87
![Calculate the potential at a point P due to a nearly spherical body. 88 Calculate the potential at a point P due to a nearly spherical body. 88](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-88.jpg)
Calculate the potential at a point P due to a nearly spherical body. 88
![Calculate the potential at a point P due to a nearly spherical body. This Calculate the potential at a point P due to a nearly spherical body. This](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-89.jpg)
Calculate the potential at a point P due to a nearly spherical body. This is Mac. Cullagh’s formula for the potential of an an ellipsoid 89
![Calculate the potential at a point P due to a nearly spherical body. So Calculate the potential at a point P due to a nearly spherical body. So](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-90.jpg)
Calculate the potential at a point P due to a nearly spherical body. So the final result for the potential has two parts – the result for the uniform sphere plus a correction for the ellipse 90
![Now we can find the force of gravity This is Mac. Cullagh’s formula for Now we can find the force of gravity This is Mac. Cullagh’s formula for](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-91.jpg)
Now we can find the force of gravity This is Mac. Cullagh’s formula for the gravity of an ellipsoid. 91
![Differential form of Newton’s law So far we’ve looked at the “integral” form for Differential form of Newton’s law So far we’ve looked at the “integral” form for](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-92.jpg)
Differential form of Newton’s law So far we’ve looked at the “integral” form for Newton’s gravitational force law. But we also have Which is a differential equation for the potential U. Can we relate U to the density without the integral? 92
![Poisson’s and Laplace’s equations Start with Gauss’s/Divergence theorem for vector fields Which says the Poisson’s and Laplace’s equations Start with Gauss’s/Divergence theorem for vector fields Which says the](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-93.jpg)
Poisson’s and Laplace’s equations Start with Gauss’s/Divergence theorem for vector fields Which says the flux out of a volume equals the divergence throughout the volume. 93
![Examine field at point M. Point M inside volume Point M outside volume r Examine field at point M. Point M inside volume Point M outside volume r](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-94.jpg)
Examine field at point M. Point M inside volume Point M outside volume r 94 Ahern: http: //geophysics. ou. edu/solid_earth/notes/laplace. html
![Examine field at point M. Point M inside volume Point M outside volume 95 Examine field at point M. Point M inside volume Point M outside volume 95](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-95.jpg)
Examine field at point M. Point M inside volume Point M outside volume 95 Ahern: http: //geophysics. ou. edu/solid_earth/notes/laplace. html
![Examine field at point M. Point M inside volume Point M outside volume 96 Examine field at point M. Point M inside volume Point M outside volume 96](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-96.jpg)
Examine field at point M. Point M inside volume Point M outside volume 96 Ahern: http: //geophysics. ou. edu/solid_earth/notes/laplace. htm
![Examine field at point M. Point M inside volume Point M outside volume So Examine field at point M. Point M inside volume Point M outside volume So](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-97.jpg)
Examine field at point M. Point M inside volume Point M outside volume So the equation for the potential, a scalar field (easier to work with than a vector field) satisfies Poisson’s equation (Lapalce’s equation is a special case of Poisson’s equation). Poisson’s equation is linear, so we can superimpose sol’ns – 97 Ahern: http: //geophysics. ou. edu/solid_earth/notes/laplace. htm
![In the spherical shell example we used the fact that gravity is “linear” i. In the spherical shell example we used the fact that gravity is “linear” i.](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-98.jpg)
In the spherical shell example we used the fact that gravity is “linear” i. e. we get final result by adding up partial results (this is what integration does!) So ellipsoidal earth can be represented as a solid sphere plus a hollow elliposid. Result for the gravity potential and force for an elliposid had two parts – that for a sphere plus an additional term which is due to the mass in the 98
![GRAVITY POTENTIAL All gravity fields satisfy Laplace’s equation in free space or material of GRAVITY POTENTIAL All gravity fields satisfy Laplace’s equation in free space or material of](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-99.jpg)
GRAVITY POTENTIAL All gravity fields satisfy Laplace’s equation in free space or material of density r. If V is the gravitational potential then (Herring) 99
![LINEAR NON-LINEAR • Superposition: break big problems into pieces • No superposition: solve whole LINEAR NON-LINEAR • Superposition: break big problems into pieces • No superposition: solve whole](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-100.jpg)
LINEAR NON-LINEAR • Superposition: break big problems into pieces • No superposition: solve whole problem at once • Smooth, predictable motions • Erratic, aperiodic motion • Response proportional to stimulus • Response need not be proportional to stimulus • Find detailed trajectories of individual particles • Find global, qualitative description of all possible trajectories 100
![Linearity and Superposition Says order you do the “combination” does not matter. Very important Linearity and Superposition Says order you do the “combination” does not matter. Very important](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-101.jpg)
Linearity and Superposition Says order you do the “combination” does not matter. Very important concept. If system is linear you can break it down into little parts, solve separately and combine solutions of parts into solution for whole. 101
![Net force of Gravity on line between Earth and Moon Solve force from Earth Net force of Gravity on line between Earth and Moon Solve force from Earth](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-102.jpg)
Net force of Gravity on line between Earth and Moon Solve force from Earth and force from Moon and add them. Probably did this procedure without even thinking about it. (earth and moon are spherical shells, so g linear inside, 0 in 102
![Net force of Gravity for Earth composed of two spherical shells – the surface Net force of Gravity for Earth composed of two spherical shells – the surface](http://slidetodoc.com/presentation_image/1d960a1c6e01c4a9231ed0e29de6d6fb/image-103.jpg)
Net force of Gravity for Earth composed of two spherical shells – the surface and a concentric "Core" Solve force from Earth and force from Core and add them. Same procedure as before (and same justification) – but probably had to think about it here. (Earth and core again spherical shells so 103
- Slides: 103