Introduction to the Use of Gravity Measurements in Ge111a

Total Page:16

File Type:pdf, Size:1020Kb

Introduction to the Use of Gravity Measurements in Ge111a Introduction to the use of gravity measurements in Ge111A • Background & basic intuition • Data reduction • Use of the instrument See Reynolds for detailed quantitative discussion What and Why Gravity measures the distribution of mass •Variations in density •Variations in interface locations •Very non-unique => use in conjunction with other constraints (geologic, seismic,…) Targets include: •Topographic compensation - the crustal root (long λ) •Fault structures/offsets (best for dip-slip faults or large offset strike-slip) •Aquifers •Tunnels •Buried objects (tunnels, ore deposits,…) •Basin sediment fill depth (together with magnetic field variations) Data Availability Regional: Carry the gravimeter around •Land, boat, plane •Arbitrary resolution Global: Satellite tracking/altimetry (over water)/other •~2 km resolution over oceans •~100 km resolution from harmonic models •None of these useful for our purposes Basic intuition Everything we need can be worked out from basic potential field theory (see Reynolds or Turcotte & Schubert): ∇2 A = 0 where A is the gravitational potential ∂A g = − z ∂z where gz is the vertical component of the gravity field Units: 1 mGal = 10-5 m/s2 1 Gal = 1 cm/s2 (named after Galileo) 1 g.u. = 0.1 mGal = 10-6 m/s2 Attraction from harmonic distribution of mass I Z X ∆ρhk cos(kx)= ∆σ k cos(kx) Why? •Intuition •Decompose any mass distribution into individual harmonic (Fourier) contributions and sum linearly •Example: Topography, density variations on a 2D or 3D grid Attraction from harmonic distribution of mass II We won’t derive solutions today, but you can imagine integrating Z Gδm δg = rˆ X r 2 ∆ρhk cos(kx)= ∆σ k cos(kx) over our given harmonic mass distribution. The vertical component is: −kz −kz gz = 2πG∆ρhk cos(kx)e = 2πG∆σ k cos(kx)e where k = 2π/λ, h<<λ The e-kz term corresponds to the 1/r2 decay with distance (the peaks and troughs tend to cancel when viewed from increasingly large distances) Attraction from harmonic distribution of mass III At long λ, k⇒0 gz = 2πG∆ρh = 2πG∆σ •No z-dependence •Attraction from an infinite sheet of mass •Will be used in making the Bouguer correction For z≠0 •If z increases (or k increases ⇔ λ decreases): e-kz ⇒0 •Upward continuation: e.g., ∆Z=λ, e-2π=0.002 (faint signals) •Downward continuation: e.g., ∆Z=-λ, e2π=500 (amplification of short λ noise) Standard gravity corrections • earth tide correction, which corrects for tidal effects of moon & sun • instrument drift correction to compensate for drift in the instrument's spring; • the latitude correction, which incorporates the variation of the Earth's gravity with latitude; • the free-air correction, which accounts for the variation in gravity due to elevation relative to sea level; • the Bouguer correction, which corrects for the attraction of material between the station and sea level; • the curvature correction, which corrects the Bouguer correction for the effect of the Earth's curvature; • the terrain correction, which removes the effect of topography to a radial distance of 166.7 km; • the isostatic correction, which removes long-wavelength variations in the gravity field inversely related to topography. • The Eötvös correction for a moving platform Free Air Correction •Bring all observations to a common horizontal reference surface (E.g., a datum or geoid) •Effect of distance not mass •Results in “Free Air Gravity” •Commonly used for airborne and marine data Gm Gm g = ; g = 0 R2 h (R+ h)2 δg g − g 2g FA = h 0 ≈ 0 = −0.3086 mGals/m h h R •The L&R gravimeter has a sensitivity of ~3µGal δg R h = FA ≈1cm 2g0 •Careful with surveying (know your target) Bouguer Correction Effect of mass (not distance) - Assume λ>>h, ignores dh/dx +… gB = 2πG∆ρh ∆ρ h h x What density to use? Standard is 2670 kg/m3 You can also use this to estimate density •Nettleton’s Approach - assume all variations are due to topography •Clever differencing of gravity profiles is sometimes useful More sophisticated approaches •Terrain corrections (graphical) •Talwani (2D and 2.5D prisms - Matlab) •Parker’s Method (finite amplitude harmonic summation) •Common Shapes (see textbook) The LaCoste & Romberg Model G Gravimeter Read Reynolds p. 50-51 for how this works. •Needs power to keep isothermal - keep plugged in somewhere •Highly sensitive spring and mass system (read Reynolds) - Do not move it without clamping the mass Getting data •Level gravimeter (you want all of gz) •Unclamp •Get a reading by adjusting the spring s.t. mass is in the “right place” - can be tedious if you have changed location a lot •A manual raw data point is also derived by reading off the “counts” and converting to mGals using a calibrated table (or Matlab) •Typically, take a manual reading at a base station, then lock the beam and use the feedback system to record values relative to the base station reading •Reclamp •Return to base station every hour for tide/drift correction. The next 3 figures relate to the problem of a terrain correction. The circular chart (hammer chart) is used to assess the effect of mass at different distances from the measurement point. See Reynolds p. 61-65 for more details. Isostatic residual gravity map (Langenheim & Jachens, 2002, BSSA v. 92(4) p. 1606-1620). Magnetic anomaly map (Langenheim & Jachens 2002) Remember from last week: magnetic anomaly strength goes as 1/r3; gravity anomaly goes as 1/r2. Magnetic potential goes as 1/r2; gravitational potential goes as 1/r. Therefore the magnetic potential and the gravity anomaly have similar form (if they arise from an object of uniform density & magnetization). See figure to the left showing that the magnetic anomaly looks like the derivative of the gravity anomaly. If you integrate the magnetic anomaly you can make a “pseudogravity anomaly” which can help with your 3-E structural interpretation. Read about this and see Fig. 6 in the Langenheim & Jachens paper. Figures at left from Langenheim & Jachens 2002 Homework Millikan library, Caltech Floor Gravity (mGal) 1. The measurements listed at right were made sequentially in the Caltech Base 0.990 library building during the course of Third -1.891 one hour. “Base” is outside the door Sixth -5.817 at ground level. What is distance Ninth -9.700 between each floor? How high is the Roof -11.263 library? Show your work. Explain what assumptions you are making. Base 0.970 2. Look up the coordinates (lat-long-elevation) of Caltech and Bombay Beach, on the E side of the Salton Sea. How will the gravity change between Caltech and our field area? Again, show your work..
Recommended publications
  • The Evolution of Earth Gravitational Models Used in Astrodynamics
    JEROME R. VETTER THE EVOLUTION OF EARTH GRAVITATIONAL MODELS USED IN ASTRODYNAMICS Earth gravitational models derived from the earliest ground-based tracking systems used for Sputnik and the Transit Navy Navigation Satellite System have evolved to models that use data from the Joint United States-French Ocean Topography Experiment Satellite (Topex/Poseidon) and the Global Positioning System of satellites. This article summarizes the history of the tracking and instrumentation systems used, discusses the limitations and constraints of these systems, and reviews past and current techniques for estimating gravity and processing large batches of diverse data types. Current models continue to be improved; the latest model improvements and plans for future systems are discussed. Contemporary gravitational models used within the astrodynamics community are described, and their performance is compared numerically. The use of these models for solid Earth geophysics, space geophysics, oceanography, geology, and related Earth science disciplines becomes particularly attractive as the statistical confidence of the models improves and as the models are validated over certain spatial resolutions of the geodetic spectrum. INTRODUCTION Before the development of satellite technology, the Earth orbit. Of these, five were still orbiting the Earth techniques used to observe the Earth's gravitational field when the satellites of the Transit Navy Navigational Sat­ were restricted to terrestrial gravimetry. Measurements of ellite System (NNSS) were launched starting in 1960. The gravity were adequate only over sparse areas of the Sputniks were all launched into near-critical orbit incli­ world. Moreover, because gravity profiles over the nations of about 65°. (The critical inclination is defined oceans were inadequate, the gravity field could not be as that inclination, 1= 63 °26', where gravitational pertur­ meaningfully estimated.
    [Show full text]
  • Student Academic Learning Services Pounds Mass and Pounds Force
    Student Academic Learning Services Page 1 of 3 Pounds Mass and Pounds Force One of the greatest sources of confusion in the Imperial (or U.S. Customary) system of measurement is that both mass and force are measured using the same unit, the pound. The differentiate between the two, we call one type of pound the pound-mass (lbm) and the other the pound-force (lbf). Distinguishing between the two, and knowing how to use them in calculations is very important in using and understanding the Imperial system. Definition of Mass The concept of mass is a little difficult to pin down, but basically you can think of the mass of an object as the amount of matter contain within it. In the S.I., mass is measured in kilograms. The kilogram is a fundamental unit of measure that does not come from any other unit of measure.1 Definition of the Pound-mass The pound mass (abbreviated as lbm or just lb) is also a fundamental unit within the Imperial system. It is equal to exactly 0.45359237 kilograms by definition. 1 lbm 0.45359237 kg Definition of Force≡ Force is an action exerted upon an object that causes it to accelerate. In the S.I., force is measured using Newtons. A Newton is defined as the force required to accelerate a 1 kg object at a rate of 1 m/s2. 1 N 1 kg m/s2 Definition of ≡the Pound∙ -force The pound-force (lbf) is defined a bit differently than the Newton. One pound-force is defined as the force required to accelerate an object with a mass of 1 pound-mass at a rate of 32.174 ft/s2.
    [Show full text]
  • What Is Fluid Mechanics?
    What is Fluid Mechanics? A fluid is either a liquid or a gas. A fluid is a substance which deforms continuously under the application of a shear stress. A stress is defined as a force per unit area Next, what is mechanics? Mechanics is essentially the application of the laws of force and motion. Fluid statics or hydrostatics is the study of fluids at rest. The main equation required for this is Newton's second law for non-accelerating bodies, i.e. ∑ F = 0 . Fluid dynamics is the study of fluids in motion. The main equation required for this is Newton's second law for accelerating bodies, i.e. ∑ F = ma . Mass Density, ρ, is defined as the mass of substance per unit volume. ρ = mass/volume Specific Weight γ, is defined as the weight per unit volume γ = ρ g = weight/volume Relative Density or Specific Gravity (S.g or Gs) S.g (Gs) = ρsubstance / ρwater at 4oC Dynamic Viscosity, µ µ = τ / (du/dy) = (Force /Area)/ (Velocity/Distance) Kinematic Viscosity, v v = µ/ ρ Primary Dimensions In fluid mechanics there are only four primary dimensions from which all other dimensions can be derived: mass, length, time, and temperature. All other variables in fluid mechanics can be expressed in terms of {M}, {L}, {T} Quantity SI Unit Dimension Area m2 m2 L2 Volume m3 m3 L3 Velocity m/s ms-1 LT-1 Acceleration m/s2 ms-2 LT-2 -Angular velocity s-1 s-1 T-1 -Rotational speed N Force kg m/s2 kg ms-2 M LT-2 Joule J energy (or work) N m, kg m2/s2 kg m2s-2 ML2T-2 Watt W Power N m/s Nms-1 kg m2/s3 kg m2s-3 ML2T-3 Pascal P, pressure (or stress) N/m2, Nm-2 kg/m/s2 kg m-1s-2 ML-1T-2 Density kg/m3 kg m-3 ML-3 N/m3 Specific weight kg/m2/s2 kg m-2s-2 ML-2T-2 a ratio Relative density .
    [Show full text]
  • Department of the Interior U.S. Geological Survey Two FORTRAN
    Department of the Interior U.S. Geological Survey Two FORTRAN programs for the reduction of atmospheric data collected by aircraft along EDM survey lines by Elliot T. Endol Gene Y. Iwatsubo Lyn J. Topinka Open-File Report 85-279 This report is preliminary and has not been edited or reviewed for conformity with Geological Survey standards and nomenclature. Any use of trade names and trade marks in this publication is for descriptive purposes only and does not constitute endorsement by the U.S. Geological Survey. David A. Johnston Cascades Volcano Observatory 5400 MacArthur Blvd. Vancouver, Washington 98661 1985 Abstract The Cascades Volcano Observatory had a requirement for a computer program to reduce atmospheric data collected for electronic distance measurement surveys. Existing programs could not be easily adapted for this simple need. The program AVEINDEX.FOR described in this report fills this requirement. AVEINDEX.FOR, written in FORTRAN-77 standards, is modular in structure. The modular structure allows easy program modification. Twelve supporting subroutines and functions for AVEINDEX.FOR and one data entry program are also described. Program comments accompany each program and most subroutines. TABLE OF CONTENTS Page INTRODUCTION................................................... 1 PRIMARY ASSUMPTIONS OR APPROXIMATIONS.......................... 2 SPECIFIC EQUATIONS USED FOR CALCULATIONS....................... 4 MAIN PROGRAMS.................................................. 7 LISTING OF MAIN PROGRAMS AND SUBROUTINES......................
    [Show full text]
  • Preparation of Papers for AIAA Technical Conferences
    Gravity Modeling for Variable Fidelity Environments Michael M. Madden* NASA, Hampton, VA, 23681 Aerospace simulations can model worlds, such as the Earth, with differing levels of fidel- ity. The simulation may represent the world as a plane, a sphere, an ellipsoid, or a high- order closed surface. The world may or may not rotate. The user may select lower fidelity models based on computational limits, a need for simplified analysis, or comparison to other data. However, the user will also wish to retain a close semblance of behavior to the real world. The effects of gravity on objects are an important component of modeling real-world behavior. Engineers generally equate the term gravity with the observed free-fall accelera- tion. However, free-fall acceleration is not equal to all observers. To observers on the sur- face of a rotating world, free-fall acceleration is the sum of gravitational attraction and the centrifugal acceleration due to the world’s rotation. On the other hand, free-fall accelera- tion equals gravitational attraction to an observer in inertial space. Surface-observed simu- lations (e.g. aircraft), which use non-rotating world models, may choose to model observed free fall acceleration as the “gravity” term; such a model actually combines gravitational at- traction with centrifugal acceleration due to the Earth’s rotation. However, this modeling choice invites confusion as one evolves the simulation to higher fidelity world models or adds inertial observers. Care must be taken to model gravity in concert with the world model to avoid denigrating the fidelity of modeling observed free fall.
    [Show full text]
  • Mass Properties Measurement Handbook
    SAWE PAPER No.2444 Index Category No. 3.0 Mass Properties Measurement Handbook By Richard Boynton and Kurt Wiener Space Electronics, Inc. 81 Fuller Way Berlin, Connecticut 06037 For presentation at the Annual Conference of the Society of Allied Weight Engineers Wichita, KS 18-20 May 1998 Permission to publish this paper, in full or in part, with full credit to author may be obtained by request to: S.A.W.E. Inc. The society is not responsible for statements or opinions in papers or discussions Table of Contents 1.0 Abstract .................................................................. 2 2.0 Stepsin Making a Mass Properties measurement .................................. 2 2.1 Definethe particularmass properties you need to measure and the required measurementaccuracy .............................................. 2 2.2 Choose the correcttype ofmeasuring instrument ............................ 2 2.3 Definethe coordinate system on the object to be used as the mass properties reference axes ............................................................ 2 2.4 Definethe position ofthe objecton themass properties measuring machine ....... 2 2.5 Determine the dimensional accuracy of theobject being measured .............. 2 2.6 Designthe fixturerequired to mount the object at a precise location relativeto the measuring instrument .............................................. 3 2.7 Verifythe position of theobject on the instrument. .......................... 3 2.8 Make the mass properties measurement ................................... 3 2.9
    [Show full text]
  • General Disclaimer One Or More of the Following Statements May Affect
    General Disclaimer One or more of the Following Statements may affect this Document This document has been reproduced from the best copy furnished by the organizational source. It is being released in the interest of making available as much information as possible. This document may contain data, which exceeds the sheet parameters. It was furnished in this condition by the organizational source and is the best copy available. This document may contain tone-on-tone or color graphs, charts and/or pictures, which have been reproduced in black and white. This document is paginated as submitted by the original source. Portions of this document are not fully legible due to the historical nature of some of the material. However, it is the best reproduction available from the original submission. Produced by the NASA Center for Aerospace Information (CASI) ^I NASA Technical Paper 2147 0 The NASA Geodynamics Program: An Overview Geoud ynamics Program Office Office of Space Science and Applications I Washington, DC 20546 1983 Scientific and Technical Information Branch -11 - ---- -- r NASA Technical Paper 2147 The NASA Ge odynamics Program: An Overview %, N&SA — TY — A 14 1) Ttit bA.'-A G?.QjlMARICs L, BCGhltl: Ab GVjhV14w ( NdtiO Ual AQr0GdUtlCIj 41A U 5j;dC#j AdsiaistLatjC&j) UC AU7/Mf AUI CSCL 06G uncIds 1/46 L 44119 W Geodynamics Program Office nary 1983 - Nala-al k,fo ulcs and SP re a ^I FOREWORD The study of the solid earth and its gravity and magnetic fields has advancd significantly in the past 20 years. This advancement has been greatly assisted by the emergence of methods which use space and space- derived technology for observations of the earth unattainable bit other means.
    [Show full text]
  • Gravity Measurements and the Standards Laboratory
    L..,.lldlU!> NBS TECHNICAL NOTE 491 Gravity Measurements II And the Standards Laboratory T 0F /*t* U.S. DEPARTMENT OF COMMERCE ,2.7/ Z ^ National Bureau of Standards Q \ *"»fAU 0» ' - NATIONAL BUREAU OF STANDARDS The National Bureau of Standards ' was established by an act of Congress March 3, 1901. Today, in addition to serving as the Nation's central measurement laboratory, the Bureau is a principal focal point in the Federal Government for assuring maximum application of the physical and engineering sciences to the advancement of technology in industry and commerce. To this end the Bureau conducts research and provides central national services in four broad program areas. These are: (1) basic measurements and standards, (2) materials measurements and standards, (3) technological measurements and standards, and (4) transfer of technology. The Bureau comprises the Institute for Basic Standards, the Institute for Materials Research, the Institute for Applied Technology, the Center for Radiation Research, the Center for Computer Sciences and Technology, and the Office for Information Programs. THE INSTITUTE FOR BASIC STANDARDS provides the central basis within the United States of a complete and consistent system of physical measurement; coordinates that system with measurement systems of other nations; and furnishes essential services leading to accurate and uniform physical measurements throughout the Nation's scientific community, industry, and com- merce. The Institute consists of an Office of Measurement Services and the following
    [Show full text]
  • 2.1 Gravity and the Gravity Field of the Earth
    2.1 Gravity and the gravity field of the Earth Fundamentals The mass and shape of the Earth The Geoid Satellite measurements of gravity GRACE - Gravity Recovery and Climate Experiment Fundamentals The force that attracts all matter to the Earth is the gravity field; an example of the general phenomenon of attraction between any masses. Newton's Law describes the force of attraction between two point masses, M1 and M2 separated by r: MM12 FG 2 r (2.1.1) The force per unit mass, F/M2 defines the vector gravity field or the gravitational acceleration, g when M1 is the Earth (Me) and r is the radius of the Earth, Re. So F M = g = G e r M R2 1 e (2.1.2) [Note that vectors are denoted either in boldface text or with an overlying arrow, e.g. r .] The gravitational constant G has been determined experimentally and is: -1- 6.67259 x 10-11 m3 kg-1 s-2 ( SI units ) or 6.67259 x 10-8 cm3 g-1 s-2 with an uncertainty of 100 parts per million (ppm). The gravitational acceleration is also measured experimentally and at the equator is about 9.8 meters/sec/sec (m s-2) or 980 cm s-2. (In almost all gravity surveys results are presented in c.g.s units rather than SI units. In c.g.s the unit acceleration, 1.0 cm s-2, is called the gal (short for Gallileo). A convenient subunit for surveys is the milligal, mgal, 10-3cm s-2.
    [Show full text]
  • A Subtlety of Earth's Gravity and Its Rotation
    A Subtlety of Earth’s Gravity and its Rotation David Alexander Lillis, 9 March 2019 This article provides an elementary introduction to the effect of Earth’s non-spherical shape (oblateness) and the centrifugal force on Earth’s gravity for non-Earth scientists, students and other interested readers. In this article I describe a simple model of the gravitational acceleration on Earth’s surface, taking into account both the gravitation of the Earth itself and the effect of centrifugal forces acting on bodies located on Earth’s surface. This model was developed purely for instructional purposes and is not intended to supplant more accurate and much more complex models already available in the relevant literature. Software written in the R language can be provided for those who wish to reproduce the calculations and graphs presented here. Part 1: Weight involves a combination of Gravity and Centrifugal Effect As is well known, Earth is not a perfect sphere. Because of its rotation the centrifugal effect makes it bulge slightly at the equator and slightly squashed at the poles, to produce what we refer to as an ‘oblate spheroid’. In fact, its observed diameter is approximately 43 km greater at the equator than at the poles, leading to slightly stronger gravitation at the poles than at the equator. However, the centrifugal effect also applies to a body on the Earth’s surface as it rotates about Earth’s rotation axis. This force appears to modify the effect of pure gravity and tends to reduce the quantity we call ‘weight’. Figure 1 shows both the gravitational force and the centrifugal force acting on a body at Earth’s surface at different latitudes, and their resultants (the vector sums of the gravitational forces and the centrifugal forces).
    [Show full text]
  • Typical Spacecraft Contents
    Appendix A: Typical Spacecraft This appendix contains descriptions and images of a dozen spacecraft selected from the many that are currently operating in interplanetary space or have successfully completed their missions, plus one that is now preparing for launch. Included is at least one representative of each of the eight spacecraft classifications described in Chapter 7 (see page 243). The scheme of limiting coverage of each spacecraft to a two-page spread in this appendix allows the reader to easily compare the various craft, their specifications, their missions, and their classifications, but it does not allow room to list all of a spacecraft’s activities, discoveries and questions raised; indeed entire books can and have been written on each. Complete profiles of these and other spacecraft are, how- ever, readily available at a single web site: http://nssdc.gsfc.nasa.gov/planetary. Contents: Spacecraft Classification Page Voyager Flyby 294 New Horizons Flyby 296 Spitzer Observatory 298 Chandra Observatory 300 Galileo Orbiter 302 Cassini Orbiter 304 Messenger Orbiter 306 Huygens Atmospheric 308 Phoenix Lander 310 Mars Science Laboratory Rover (launch: 2009) 312 Deep Impact Penetrator 314 Deep Space 1 Engineering 316 294 Appendix A: Typical Spacecraft The Voyager Spacecraft Fig. A.1. Each Voyager spacecraft measures about 8.5 meters from the end of the science boom across the spacecraft to the end of the RTG boom. The magnetometer boom is 13 meters long. Courtesy NASA/JPL. Classification: Flyby spacecraft Mission: Encounter giant outer planets and explore heliosphere Named: For their journeys Summary: The two similar spacecraft flew by Jupiter and Saturn.
    [Show full text]
  • HOW to CALCULATE MASS PROPERTIES (An Engineer's Practical Guide)
    HOW TO CALCULATE MASS PROPERTIES (An Engineer's Practical Guide) by Richard Boynton, President Space Electronics, Inc. Berlin, CT 06037 and Kurt Wiener, Assistant Chief Engineer Space Electronics, Inc. Berlin, CT 06037 Abstract There are numerous text books on dynamics which devote a few pages to the calculation of mass properties. However, these text books quickly jump from a very brief description of these quantities to some general mathematical formulas without giving adequate examples or explaining in enough detail how to use these formulas. The purpose of this paper is to provide a detailed procedure for the calculation of mass properties for an engineer who is inexperienced in these calculations. Hopefully this paper will also provide a convenient reference for those who are already familiar with this subject. This paper contains a number of specific examples with emphasis on units of measurement. The examples used are rockets and re-entry vehicles. The paper then describes the techniques for combining the mass properties of sub-assemblies to yield the composite mass properties of the total vehicle. Errors due to misalignment of the stages of a rocket are evaluated numerically. Methods for calculating mass property corrections are also explained. Page 1 Calculation of Mass Properties using Traditional Methods Choosing the Reference Axes The first step in calculating mass properties of an object is to assign the location of the reference axes. The center of gravity and the product of inertia of an object can have any numerical value or polarity, depending on the choice of axes that are used as a reference for the calculation.
    [Show full text]