Geoid, Topography, and the Bouguer Plate Or Shell

Total Page:16

File Type:pdf, Size:1020Kb

Geoid, Topography, and the Bouguer Plate Or Shell Journal of Geodesy C2001) 75: 210±215 Geoid, topography, and the Bouguer plate or shell P. Vanõ cÆ ek1, P. Nova k1, Z. Martinec2 1 Department of Geodesy and Geomatics Engineering, University of New Brunswick, P.O. Box 4400, E3B 5A3 Fredericton, Canada e-mail: [email protected]; Tel.: +1-506-4535144; Fax: +1-506-4534943 2 Department of Geophysics, Charles University, V HolesÆ ovicÆ kach 2, Pragues, Czech Republic Received: 24December 1999 / Accepted: 11 December 2000 Abstract. Topography plays an important role in solv- ing many geodetic and geophysical problems. In the evaluation of a topographical eect, a planar model, a 1 Introduction spherical model or an even more sophisticated model can be used. In most applications, the planar model is Periodically, people discover that planar and spherical considered appropriate: recall the evaluation of gravity models of topography give very dierent results for reductions of the free-air, Poincare ±Prey or Bouguer Bouguer anomalies. Similarly, the results for the direct kind. For some applications, such as the evaluation of and indirect topographical eects in the Stokes±Helmert topographical eects in gravimetric geoid computations, technique for geoid computations obtained by means of it is preferable or even necessary to use at least the the planar and spherical models are found to be quite spherical model of topography. In modelling the topo- dierent. Some people claim that the planar model can graphical eect, the bulk of the eect comes from the safely be used for ``local work'' while the spherical Bouguer plate, in the case of the planar model, or from model has to be used for global work. Others still the Bouguer shell, in the case of the spherical model. The maintain that all these questions have already been dierence between the eects of the Bouguer plate and sorted out and that they do not require any more of our the Bouguer shell is studied, while the eect of the rest of attention. So what is going on? topography, the terrain, is discussed elsewhere. It is When looking into this problem CVanõ cÆ ek and Nova k argued that the classical Bouguer plate gravity reduction 1999) we discovered an interesting story, which we will should be considered as a mathematical construction try to recount here. To do so, we focus only on the with unclear physical meaning. It is shown that if the ``in®nite plate'' and the ``spherical shell'' models, leaving reduction is understood to be reducing observed gravity out the terrain eects. The corresponding dierence of onto the geoid through the Bouguer plate/shell then the planar and spherical models of terrain presents an- both models give practically identical answers, as other fascinating story which, in our opinion, requires a associated with Poincare 's and Prey's work. It is shown separate and rather more extensive paper to deal with it why only the spherical model should be used in the adequately. The main point of this ``terrain story'' is the evaluation of topographical eects in the Stokes±Helmert discovery that, contrary to popular belief, the spherical solution of Stokes' boundary-value problem. The reason model terrain eect has to be considered globally. This for this is that the Bouguer plate model does not allow for point has been already discussed by Nova k et al. C1998) a physically acceptable condensation scheme for the and by Nova k and Vanõ cÏ ek C1999). A more formal and topography. complete paper on the subject of spherical terrain model is under review CNova k et al. submitted). Once pub- lished, it will complement the present paper. Key words: Geoid ± Bouguer Reduction of Gravity ± In the present paper, we will focus our attention ®rst Stokes±Helmert's Problem on the dierence between the Bouguer plate and Bou- guer shell eects on Bouguer anomalies. Then, we will tackle another, more or less independent issue, that of the dierence between the two models in the gravimetric geoid determination by means of the Stokes±Helmert technique, where the use of the two models has pro- Correspondence to: P. Vanõ cÆ ek, P. Nova k is currently at: Depart- found consequences. In fact, the planar model, the ment of Geomatics Engineering, The University of Calgary, 2500 Bouguer plate, cannot be used at all if we wish to use a University Drive NW, Calgary, Canada. physically meaningful condensation scheme. 211 2The story of Bouguer plate reduction cally traversing the plate or shell is the same, except for higher-order terms. The situation for the potential is In order to see the pattern, let us show the gravitational naturally dierent: as the potential of the in®nite plate is potential, the gravitational attraction Cnegative ®rst in®nite, we cannot make any direct comparison between vertical derivative of the potential), and the vertical the two models. We can only observe that in the spherical gradient of gravitational attraction Cnegative second model, the potentials at the upside and on the downside of derivative of the potential) of the topographical CBou- the shell dier only by higher-order terms. guer) plate and the topographical CBouguer) shell side by What does it all mean? We wish to address here only side. To keep things simple, let us assume a constant the most interesting question of what this means in the density . Csay, the usual 2670 kg mÀ3) and the same context of the Cincomplete, i.e. without the terrain cor- thickness, H, for both the in®nite plate and the shell of rection) Bouguer gravity anomaly. The incomplete or the inner radius R. This is all shown in Fig. 1. The three simple Bouguer anomaly is computed from the follow- quantities of interest are computed at two points: one on ing formula: the top and one at the bottom of the plate/shell. In B addition, the second derivative, which is discontinuous Dg g A A À c 1 on the top Cand also at the bottom) of the plate/shell, is where g is the observed gravity on the Earth's surface Cat at this point computed in both directions: from above altitude H), A is the ``free-air reduction'' Cto the geoid) and from below. The expressions for the plate are due to the Earth's masses enclosed within the inner derived from Eqs. C3.5) and C3.7) in Heiskanen and radius of the spherical shell Cincluding the latitude and Moritz C1967) by simply extending the ®nite plate to altitude terms), AB is the ``Bouguer reduction'' Cto the in®nity. The expressions for the shell are derived directly geoid) due to the mass of the Bouguer plate, and c is the from the equations for the potential Cof a spherical shell) normal gravity at the reference ellipsoid [Heiskanen and in Wichiencharoen [1982, Eqs. C19), C24) and C25)]. Moritz 1967, Eq. C3.19)]. In this case, the so-called Now, examining Fig. 1, how dierent really are the Bouguer reduction is given by results for the planar and spherical models? Starting from the bottom, with the vertical gradient of attraction, and AB 2pG.H 2 neglecting the higher-order terms Cof the order of H=R and smaller) in the spherical model, the results are identical. For the standard value of mean topographical density of The attraction of the plate at its top is only one half of that 2670 kg mÀ3, the numerical value of the Bouguer of the shell Cat its top and neglecting the higher-order reduction is 0:1119  H mGal. terms), while the attraction at the bottom of the plate is Inspecting again Fig. 1, we can see rather easily that exactly opposite to that at the top. The attraction at the AB is not the dierence between the gravity values at the downside of the shell is zero, as it should be CKellogg top and the bottom of the in®nite plate! It is not the 1929). Note that the change in the attraction when verti- dierence between gravity values on the upside and Fig. 1. The behaviour of gravity and gravity potential induced by Bouguer shell and Bouguer plate. The symbols Oi H=R de- note higher-order terms 212 downside of the spherical shell either! This dierence, in the use of the Cphysically less satisfactory) planar model the planar case, is equal to and that we should therefore obtain the ``correct'' result by using the physically more satisfying spherical model. PP A 4pG.H 3 Rather, the Bouguer gravity reduction is the conse- known in geodesy as the ``Poincare ±Prey gravity reduc- quence of the physically unsatisfactory ``removal of the Bouguer plate'' of an in®nite mass. Thus the Bouguer tion'' [Heiskanen and Moritz 1967, Eq. C3.64)]. In the spherical case the dierence is the same except for the gravity anomaly, useful as it is in many geophysical applications, is really only an arti®cial construction. The higher-order terms, which we will neglect for simplicity. CWe note that we arrive at the stated conclusions, using dierence between using the planar and the spherical either the attraction or the vertical gradient of attraction model is of second order; it is the known spherical Â Ï formulae, as we should.) For the standard topographical correction to the Bouguer plate reduction [Vanõcek and density, the numerical value of the Poincare ±Prey Krakiwsky 1986; Eq. C21.43)] reduction is 0:2238  H mGal. H 2 The discrepancy between the planar and spherical dAB À8pG. 5 S R models was pointed out by Karl C1971). It should be mentioned that planar and spherical terrain models For the standard topographical density of 2670 kg mÀ3, display the same discrepancy.
Recommended publications
  • Analysis of Different Methodologies to Calculate Bouguer Gravity Anomalies in the Argentine Continental Margin
    Geosciences 2014, 4(2): 33-41 DOI: 10.5923/j.geo.20140402.02 Analysis of Different Methodologies to Calculate Bouguer Gravity Anomalies in the Argentine Continental Margin Ana C. Pedraza De Marchi1,2,*, Marta E. Ghidella3, Claudia N. Tocho1 1Universidad Nacional de La Plata, Facultad de Ciencias Astronómicas y Geofísicas, La Plata, 1900, Argentina 2CONICET, Concejo Nacional de Investigaciones Científicas y Técnicas, Argentina 3Instituto Antártico Argentino, Buenos Aires, C1064AAF, Argentina Abstract We have tested and used two methods to determine the Bouguer gravity anomaly in the area of the Argentine continental margin. The first method employs the relationship between the topography and gravity anomaly in the Fourier transform domain using Parker’s expression for different orders of expansion. The second method computes the complete Bouguer correction (Bullard A, B and C) with the Fortran code FA2BOUG2. The Bouguer slab correction (Bullard A), the curvature correction (Bullard B) and the terrain correction (Bullard C) are computed in several zones according to the distances between the topography and the calculation point. In each zone, different approximations of the gravitational attraction of rectangular or conic prisms are used according to the surrounding topography. Our calculations show that the anomaly generated by the fourth order in Parker’s expansion is actually compatible with the traditional Bouguer anomaly calculated with FA2BOUG, and that higher orders do not introduce significant changes. The comparison reveals that the difference between both methods in the Argentine continental margin has a quasi bimodal statistical distribution. The main disadvantage in using routines based on Parker's expansion is that an average value of the topography is needed for the calculation and, as the margin has an abrupt change of the topography in the continental slope area, it causes a bimodal distribution.
    [Show full text]
  • Geoid Determination Based on a Combination of Terrestrial and Airborne Gravity Data in South Korea
    Geoid Determination based on a Combination of Terrestrial and Airborne Gravity Data in South Korea DISSERTATION Presented in Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy in the Graduate School of The Ohio State University By Hyo Jin Yang Graduate Program in Geodetic Science and Surveying The Ohio State University 2014 Dissertation Committee: Professor Christopher Jekeli, Advisor Professor Michael Bevis Professor Ralph R.B. von Frese Copyright by Hyo Jin Yang 2014 ABSTRACT The regional gravimetric geoid model for South Korea is developed by using heterogeneous data such as gravimetric measurements, a global geopotential model, and a high resolution digital topographic model. A highly accurate gravimetric geoid model, which is a basis to support the construction of the efficient and less costly height system with GPS, requires many gravimetric observations and these are acquired by several kinds of sensors or on platforms. Especially airborne gravimetry has been widely employed to measure the earth’s gravity field in last three decades, as well as the traditional measurements on the earth’s physical surface. Therefore, it is necessary to understand the characters of each gravimetric measurement, such as the measurement surface and involved topography, and also to integrate these to a unified gravimetric data base which refers to the same gravitational field. This dissertation illustrates the methods for combining two types of available gravity data for South Korea, one is terrestrial data obtained on the earth’s surface and another is airborne data measured at altitude, and shows an accessible accuracy of the geoid model based on these data.
    [Show full text]
  • Gravity Anomalies - D
    GEOPHYSICS AND GEOCHEMISTRY – Vol.III - Gravity Anomalies - D. C. Mishra GRAVITY ANOMALIES D. C. Mishra National Geophysical Research Institute, Hyderabad, India Keywords: gravity anomalies, isostasy, Free Air and Bouguer gravity anomalies Contents 1. Introduction 2. Free Air and Bouguer Gravity Anomalies 3. Separation of Gravity Anomalies 3.1 Regional and Residual Gravity Fields 3.2 Separation Based on Surrounding Values 3.3 Polynomial Approximation 3.4 Digital Filtering 4. Analytical Operations 4.1 Continuation of the Gravity Field 4.2 Derivatives of the Gravity Field 5. Isostasy 5.1 Isostatic Regional and Residual Fields 5.2 Admittance Analysis and Effective Elastic Thickness 6. Interpretation and Modeling 6.1 Qualitative Interpretation and Some Approximate Estimates 6.2 Quantitative Modeling Due to Some Simple shapes 6.2.1 Sphere 6.2.2 Horizontal Cylinder 6.2.3 Vertical Cylinder 6.2.4 Prism 6.2.5 Contact 6.3 Gravity Anomaly Due to an Arbitrary Shaped Two-dimensional Body 6.4 Basement Relief Model 7. Applications 7.1 Bouguer Anomaly of Godavari Basin, India 7.2 Spectrum and Basement Relief 7.3 Modeling of Bouguer Anomaly of Godavari Basin Along a Profile 7.4 SomeUNESCO Special Applications – EOLSS Glossary Bibliography SAMPLE CHAPTERS Biographical Sketch Summary Gravity anomalies are defined in the form of free air and Bouguer anomalies. Various methods to separate them in the regional and the residual fields are described, and their limitations are discussed. Polynomial approximation and digital filtering for this purpose suffer from the arbitrary selection of the order of polynomial and cut off frequency, respectively. However, some constraints on the order of these anomalies can ©Encyclopedia of Life Support Systems (EOLSS) GEOPHYSICS AND GEOCHEMISTRY – Vol.III - Gravity Anomalies - D.
    [Show full text]
  • The Generalized Bouguer Anomaly
    Earth Planets Space, 58, 287–303, 2006 The generalized Bouguer anomaly Kyozo Nozaki Tsukuba Technical Research and Development Center, OYO Corporation, 43 Miyukigaoka, Tsukuba, Ibaraki 305-0841, Japan (Received November 11, 2004; Revised June 28, 2005; Accepted September 20, 2005; Online published March 10, 2006) This paper states on the new concept of the generalized Bouguer anomaly (GBA) that is defined upon the datum level of an arbitrary elevation. Discussions are particularly focused on how to realize the Bouguer anomaly that is free from the assumption of the Bouguer reduction density ρB , namely, the ρB -free Bouguer anomaly, and on what is meant by the ρB -free Bouguer anomaly in relation to the fundamental equation of physical geodesy. By introducing a new concept of the specific datum level so that GBA is not affected by the topographic masses, we show the equations of GBA upon the specific datum levels become free from ρB and/or the terrain correction. Subsequently utilizing these equations, we derive an approximate equation for estimating ρB . Finally, we show how to compute a Bouguer anomaly on the geoid by transforming the datum level of GBA from the specific datum level to the level of the geoid. These procedures yield a new method for obtaining the Bouguer anomaly in the classical sense (say, the Bouguer disturbance), which is free from the assumption of ρB . We remark that GBA upon the ρB -free datum level is the gravity disturbance and that the equation of it has a tie to the fundamental equation of physical geodesy. Key words: Generalized Bouguer anomaly, Poincare-Prey´ reduction, specific datum level of gravity reduction, free-air anomaly, Bouguer reduction density.
    [Show full text]
  • DOGAMI Open-File Report O-82-09, Gravity Anomalies in the Cascade
    STATE OF OREGON DEPARTMENT OF GEOLOGY AND MINERAL INDUSTRIES 1005 State Office Building Portland, Oregon 97201 OPEN-FILE REPORT 0-82-9 GRAVITY ANOMALIES IN THE CASCADE RANGE IN OREGON: STRUCTURAL AND THERMAL IMPLICATIONS R.W. Couch, G.S. Pitts, M. Gemperle, D.E. Braman, and C.A. Veen Geophysics Group, School of Oceanography, Oregon State University, Corvall is, OR 97331 Funds for this research were provided by the U. S. Geological Survey under Grant Number 14-08-0001 -6-393 and by the U.S. Department of Energy through the Oregon Department of Geology and Mineral Industries under Cooperative Agreement Numbers DE-FC07-79ID12044 and DE- FC07-79ET27220. NOTICE The Oregon Department of Geology and Mineral Industries is publishing this paper because the subject matter is consistent with the mission of the Department. To facilitate timely distribution of information, camera-ready copy submitted by the authors has not been edited by the staff of the Oregon Department of Geology and Mineral Industries. Tab1 e of Contents Page List of Illustrations i i Introduction Gravity Measurements Calculation of the Free-ai r and Bouguer Gravity Anomalies Free-ai r Gravity Anomal ies Bouguer Gravity Anomal ies The Separation of Regional and Residual Gravity Anomalies Regional Gravity Anornal ies Residual Gravity Anomalies of the Cascade Range in Oregon Gravity Anomaly Lineations in the Cascade Range in Oregon The Northern Oregon Crustal Cross Section Geophysical Cross Section of the Cascade Range in Southern Oregon Structural and Geothermal Implications of the Cascade Gravity Anomalies Acknowledgements References Figure Captions List of Illustrations Figures Page 1 Regional geology of the Pacific Northwest 46 2 Topographic map of the Cascade Range in northern 47 Oregon.
    [Show full text]
  • Determination of the Selected Gravity Field Functionals by the GGI Method: a Case Study of the Western Carpathians Area
    applied sciences Article Determination of the Selected Gravity Field Functionals by the GGI Method: A Case Study of the Western Carpathians Area Marek Trojanowicz 1,* , Magdalena Owczarek-Wesołowska 1, Lubomil Pospíšil 2 and Olgierd Jamroz 1 1 Institute of Geodesy and Geoinformatics, Wroclaw University of Environmental and Life Sciences, C. K. Norwida 25, 50-375 Wroclaw, Poland; [email protected] (M.O.-W.); [email protected] (O.J.) 2 Department of Geodesy and Mining Surveying, VŠB—Technical University Ostrava, 17. listopadu 2172/15, 708 00 Ostrava—Poruba, Czech Republic; [email protected] * Correspondence: [email protected] Received: 2 September 2020; Accepted: 29 October 2020; Published: 6 November 2020 Abstract: In this paper, some features of the local disturbing potential model developed by the GGI method (based on Geophysical Gravity Inversion) were analyzed. The model was developed for the area of the Western Carpathians covering the Polish–Slovak border. A detailed assessment of the model’s property was made regarding the accuracy of the disturbing potential values (height anomalies), gravity values, complete Bouguer anomalies (CBA), and differences between geoid undulations and height anomalies (N ζ). Obtained accuracies of the GGI quasigeoid model − (in terms of standard deviation of the residuals to the reference quasigeoid models) were at the level of 2.2 cm for Poland and 0.9 cm for the Slovak area. In terms of gravity, there was shown ± ± dependence of the accuracy of the GGI model on the digital elevation model (DEM) resolution, the point height, the density of gravity data used, and used reference density of topography model.
    [Show full text]
  • Analysis of Different Methodologies to Calculate Bouguer Gravity Anomalies in the Argentine Continental Margin
    Geosciences 2014, 4(2): 33-41 DOI: 10.5923/j.geo.20140402.02 Analysis of Different Methodologies to Calculate Bouguer Gravity Anomalies in the Argentine Continental Margin Ana C. Pedraza De Marchi1,2,*, Marta E. Ghidella3, Claudia N. Tocho1 1Universidad Nacional de La Plata, Facultad de Ciencias Astronómicas y Geofísicas, La Plata, 1900, Argentina 2CONICET, Concejo Nacional de Investigaciones Científicas y Técnicas, Argentina 3Instituto Antártico Argentino, Buenos Aires, C1064AAF, Argentina Abstract We have tested and used two methods to determine the Bouguer gravity anomaly in the area of the Argentine continental margin. The first method employs the relationship between the topography and gravity anomaly in the Fourier transform domain using Parker’s expression for different orders of expansion. The second method computes the complete Bouguer correction (Bullard A, B and C) with the Fortran code FA2BOUG2. The Bouguer slab correction (Bullard A), the curvature correction (Bullard B) and the terrain correction (Bullard C) are computed in several zones according to the distances between the topography and the calculation point. In each zone, different approximations of the gravitational attraction of rectangular or conic prisms are used according to the surrounding topography. Our calculations show that the anomaly generated by the fourth order in Parker’s expansion is actually compatible with the traditional Bouguer anomaly calculated with FA2BOUG, and that higher orders do not introduce significant changes. The comparison reveals that the difference between both methods in the Argentine continental margin has a quasi bimodal statistical distribution. The main disadvantage in using routines based on Parker's expansion is that an average value of the topography is needed for the calculation and, as the margin has an abrupt change of the topography in the continental slope area, it causes a bimodal distribution.
    [Show full text]
  • On Different Techniques for the Calculation of Bouguer Gravity
    University of Texas at El Paso DigitalCommons@UTEP Open Access Theses & Dissertations 2013-01-01 On Different Techniques For The alcC ulation Of Bouguer Gravity Anomalies For Joint Inversion And Model Fusion Of Geophysical Data In The Rio Grande Rift Azucena Zamora University of Texas at El Paso, [email protected] Follow this and additional works at: https://digitalcommons.utep.edu/open_etd Part of the Applied Mathematics Commons, and the Geophysics and Seismology Commons Recommended Citation Zamora, Azucena, "On Different Techniques For The alcC ulation Of Bouguer Gravity Anomalies For Joint Inversion And Model Fusion Of Geophysical Data In The Rio Grande Rift" (2013). Open Access Theses & Dissertations. 1965. https://digitalcommons.utep.edu/open_etd/1965 This is brought to you for free and open access by DigitalCommons@UTEP. It has been accepted for inclusion in Open Access Theses & Dissertations by an authorized administrator of DigitalCommons@UTEP. For more information, please contact [email protected]. ON DIFFERENT TECHNIQUES FOR THE CALCULATION OF BOUGUER GRAVITY ANOMALIES FOR JOINT INVERSION AND MODEL FUSION OF GEOPHYSICAL DATA IN THE RIO GRANDE RIFT AZUCENA ZAMORA Program in Computational Science APPROVED: Aaron Velasco, Chair, Ph.D. Rodrigo Romero, Ph.D. Vladik Kreinovich, Ph.D. Benjamin C. Flores, Ph.D. Dean of the Graduate School Para mi familia con amor Sin ustedes nada de esto seria posible ON DIFFERENT TECHNIQUES FOR THE CALCULATION OF BOUGUER GRAVITY ANOMALIES FOR JOINT INVERSION AND MODEL FUSION OF GEOPHYSICAL DATA
    [Show full text]
  • New Standards for Reducing Gravity Data: the North American Gravity Database William J
    University of South Carolina Scholar Commons Faculty Publications Earth, Ocean and Environment, School of the 2005 New standards for reducing gravity data: The North American gravity database William J. Hinze Purdue University, [email protected] Carlos Aiken University of Texas at Dallas John Brozena Bernard Coakley University of Alaska Fairbanks David Dater See next page for additional authors Follow this and additional works at: https://scholarcommons.sc.edu/geol_facpub Part of the Earth Sciences Commons Publication Info Published in Geophysics, Volume 70, Issue 4, 2005, pages J25-J32. Hinze, W. J., Aiken, C., Brozena, J., Coakley, B., Dater, D., Flanagan, G., Forsberg, R., Hildenbrand, T., Keller, G. R., Kellogg, J., Kucks, R., Li, X., Mainville, A., Morin, R., Pilkington, M., Plouff, D., Ravat, D., Roman, D., Urrutia-Fucugauchi, J., Veronneau, M., Webring, M., & Winester, D. (2005). New standards for reducing gravity data: The orN th American gravity database. Geophysics, 70 (4), J25-J32. ©Geophysics 2005, Society of Exploration Geophysicists This Article is brought to you by the Earth, Ocean and Environment, School of the at Scholar Commons. It has been accepted for inclusion in Faculty Publications by an authorized administrator of Scholar Commons. For more information, please contact [email protected]. Author(s) William J. Hinze, Carlos Aiken, John Brozena, Bernard Coakley, David Dater, Guy Flanagan, Rene Forsberg, Thomas Hildenbrand, G. Randy Keller, James N. Kellogg, Robert Kucks, Xiong Li, Andre Mainville, Robert Morin, Mark Pilkington, Donald Plouff, Dhananjay Ravat, Daniel Roman, Jamie Urrutia-Fucugauchi, Marc Veronneau, Michael Webring, and Daniel Winester This article is available at Scholar Commons: https://scholarcommons.sc.edu/geol_facpub/2 GEOPHYSICS, VOL.
    [Show full text]
  • Should Geophysicists Use the Gravity Disturbance Or the Anomaly?
    Should geophysicists use the gravity disturbance or the anomaly? Vanderlei C. Oliveira Jr∗, Leonardo Uieday, Kristoffer A. T. Hallam∗ and Val´eriaC. F. Barbosa∗ ∗Observat´orioNacional, Department of Geophysics, Rio de Janeiro, Brazil yUniversity of Hawai'i at M¯anoa, Department of Geology and Geophysics, SOEST, Honolulu, USA ABSTRACT The gravity anomaly is defined as the difference between the Earth's gravity on the geoid and the normal gravity on the reference ellipsoid. Because these quantities are not at the same point, the anomaly contains centrifugal accelerations and cannot be considered a harmonic function. The gravity disturbance is the difference between gravity and normal gravity at the same point. Consequently, the centrifugal effects can be neglected and the disturbance can be considered a harmonic function. This is the premise behind most potential-field data processing techniques (e.g., upward/downward continuation). Unlike the anomaly, the disturbance is due solely to the gravitational effects of geologic sources, making it the most appropriate for geophysical purposes. Use of the gravity anomaly in geophysics carries with it the implicit assumption that it is a good approximation for the gravity disturbance. However, bear in mind that the difference between the gravity disturbance and the free-air anomaly can be larger than 10 mGal worldwide. In fact, we argue that the assumptions made during gravity forward and inverse modeling imply that the quantity being modeled is the disturbance, not the anomaly. INTRODUCTION The gravity vector is the sum of the gravitational and centrifugal accelerations felt by a body at rest on the Earth's surface.
    [Show full text]
  • Gravity and Aeromagnetic Anomaly Maps of the Southern Peninsula Of
    The State Geological Survey collects, interprets, and Instrumentation...............................................................5 disseminates basic information on the geology and Survey Procedure ..........................................................5 mineral resources of Michigan. Reduction of Data ..........................................................5 Its activities are guided by public service available to all who are interested in the use or development of our GEOLOGICAL SIGNIFICANCE OF ANOMALY MAPS...6 resources, the protection of our environment, and sound Source of Anomalies ......................................................6 land use management. DISCUSSION OF SELECTED ANOMALIES ...................8 See page 11 for a list of other publications on geophysics and subsurface geology available from the Geological REFERENCES ................................................................10 Survey. Report of Investigation 14 ILLUSTRATIONS Gravity and Aeromagnetic Anomaly Maps of Figure 1. Normal Geomagnetic Field of the Southern Peninsula the Southern Peninsula of Michigan of Michigan..........................................................................6 Figure 2. Relief Model of the Bouguer Gravity Anomaly Field of WILLIAM J. HINZE* the Southern Peninsula of Michigan as Observed Along a Department of Geosciences Line of Sight of N 30° W and an Inclination of 45°...............7 Purdue University West Lafayette, Indiana Figure 3. Theoretical Gravity Anomalies from a Simulated Michigan Basin....................................................................8
    [Show full text]
  • LAB 4: Processing of the KM1003 Geophysics Data Due Tue, March 16 This Is a Two-Part Lab
    GG450L:Geophysical Methods Lab p. 1 LAB 4: Processing of the KM1003 Geophysics Data Due Tue, March 16 This is a two-part Lab. The purpose of this lab is to learn how to manipulate and process real shipboard geophysics data. Your goals will be to first reduce the raw gravity data to free-air anomaly, and then do a short analysis of the simple Bouguer anomaly. Second you will remove the IGRF reference field of the magnetic data and discuss the results. Part I. 1. Reducing raw gravimeter readings to free-air gravity anomaly This part of the data processing will require converting the gravimeter dial readings to mGal (using “dial_to_mGal.m”), removing the Earth’s normal field (gN), applying the Eӧtvӧs correction, and removing any short-wavelength noise that cannot be related to geology. The data are in the file “dat.txt” and can be plotted using the script “plotdata.m”. Columns 7-10 are day of the year (Julian day), hour, min, and second, respectively, for each measurement. The eastward component of the ship’s speed (i.e., v*cos(azimuth) in the Eӧtvӧs equation) can be computed by taking the difference in longitude between two locations multiplying by 60 (to get nautical miles) and dividing by the time between the two locations. Filtering noisy data can be done with, e.g., “medfilt1” or “smooth(Y,span,‘rloess’)” will probably be necessary. You should create your own plots to verify each step of the way. Please (at least) show the following plots on one page “subplot (311)” raw gravity in mGal and free-air anomaly (both having a value near or at zero at the start of the ship track), “subplot(312)” eastward component of the ship’s velocity, and “subplot(313)” bathymetry, all versus longitude.
    [Show full text]