MEAN GRAVITY ALONG PLUMBLINE. University of New Brunswick, Department of Geodesy and Geomatics Engineering, Fredericton, N.B.

Similar documents
GEOID-QUASIGEOID CORRECTION IN FORMULATION OF THE FUNDAMENTAL FORMULA OF PHYSICAL GEODESY. Robert Tenzer 1 Petr Vaníček 2

11. HAFAT İş-Enerji Power of a force: Power in the ability of a force to do work

Visco-elastic Layers

On Some Numerical Aspects of Primary Indirect Topographical Effect Computation in the Stokes-Helmert Theory of Geoid Determination

Part 2 KINEMATICS Motion in One and Two Dimensions Projectile Motion Circular Motion Kinematics Problems

Heat transfer between shell and rigid body through the thin heat-conducting layer taking into account mechanical contact

Computer Propagation Analysis Tools

ROLE OF NO TOPOGRAPHY SPACE IN STOKES-HELMERT SCHEME FOR GEOID DETERMINATION

Lecture 17: Kinetics of Phase Growth in a Two-component System:

Maxwell Equations. Dr. Ray Kwok sjsu

Transient Radial Flow Toward a Well Aquifer Equation, based on assumptions becomes a 1D PDE for h(r,t) : Transient Radial Flow Toward a Well

On Some Aspects of Accuracy of the Geopotential Model EMG96 at the Territory of Canada

7 Wave Equation in Higher Dimensions

2. The units in which the rate of a chemical reaction in solution is measured are (could be); 4rate. sec L.sec

Maximum Cross Section Reduction Ratio of Billet in a Single Wire Forming Pass Based on Unified Strength Theory. Xiaowei Li1,2, a

Fig. 1S. The antenna construction: (a) main geometrical parameters, (b) the wire support pillar and (c) the console link between wire and coaxial

Kinematics Review Outline

Lecture 18: Kinetics of Phase Growth in a Two-component System: general kinetics analysis based on the dilute-solution approximation

Lecture 3. Electrostatics

Lecture 22 Electromagnetic Waves

The 37th International Physics Olympiad Singapore. Experimental Competition. Wednesday, 12 July, Sample Solution

CHAPTER 24 GAUSS LAW


Relative and Circular Motion

Journal of Theoretics

, on the power of the transmitter P t fed to it, and on the distance R between the antenna and the observation point as. r r t

NSEP EXAMINATION

Mechanics and Special Relativity (MAPH10030) Assignment 3

Sharif University of Technology - CEDRA By: Professor Ali Meghdari

Automatic Measuring of English Language Proficiency using MT Evaluation Technology

Strees Analysis in Elastic Half Space Due To a Thermoelastic Strain

( ) + θ θ. ω rotation rate. θ g geographic latitude - - θ geocentric latitude - - Reference Earth Model - WGS84 (Copyright 2002, David T.

The Components of Vector B. The Components of Vector B. Vector Components. Component Method of Vector Addition. Vector Components

Neutron Slowing Down Distances and Times in Hydrogenous Materials. Erin Boyd May 10, 2005

Lecture 3: Resistive forces, and Energy

EN221 - Fall HW # 7 Solutions

Determination of a geoid model for Ghana using the Stokes-Helmert method. Michael Adjei Klu. BSc Geodetic Engineering

Physics Courseware Physics I Constant Acceleration

( ) ( ) CHAPTER 8. Lai et al, Introduction to Continuum Mechanics. Copyright 2010, Elsevier Inc 8-1

= 4 3 π( m) 3 (5480 kg m 3 ) = kg.

MATHEMATICAL FOUNDATIONS FOR APPROXIMATING PARTICLE BEHAVIOUR AT RADIUS OF THE PLANCK LENGTH

Incompressible Viscous Flows

The Flatness Problem as A Natural Cosmological Phenomenon

Review: Electrostatics and Magnetostatics

CHAPTER GAUSS'S LAW

WYSE Academic Challenge Sectional Mathematics 2006 Solution Set

STUDY OF THE BIOHEAT EQUATION WITH A SPHERICAL HEAT SOURCE FOR LOCAL MAGNETIC HYPERTHERMIA

r r r r r EE334 Electromagnetic Theory I Todd Kaiser

Chapter 13 Gravitation

STUDY OF THE STRESS-STRENGTH RELIABILITY AMONG THE PARAMETERS OF GENERALIZED INVERSE WEIBULL DISTRIBUTION

Low-complexity Algorithms for MIMO Multiplexing Systems

KINEMATICS OF RIGID BODIES

Design Guideline for Buried Hume Pipe Subject to Coupling Forces

Milos Pick Geophysical Institute of the Czech Academy of Sciences, Czech Republic.

5.1 Angles and Their Measure

Example

Newton s Laws, Kepler s Laws, and Planetary Orbits

[ ] [ ] 3.3 Given: turning corner radius, r ε = 0 mm lead angle, ψ r = 15 back rake angle, γ p = 5 side rake angle, γ f = 5

Prerna Tower, Road No 2, Contractors Area, Bistupur, Jamshedpur , Tel (0657) , PART A PHYSICS

Lecture 5. Chapter 3. Electromagnetic Theory, Photons, and Light

The Gradient and Applications This unit is based on Sections 9.5 and 9.6, Chapter 9. All assigned readings and exercises are from the textbook

( t) Steady Shear Flow Material Functions. Material function definitions. How do we predict material functions?

Solving Problems of Advance of Mercury s Perihelion and Deflection of. Photon Around the Sun with New Newton s Formula of Gravity

Application of Fractional Calculus in Food Rheology E. Vozáry, Gy. Csima, L. Csapó, F. Mohos

Work, Energy, and Power. AP Physics C

Coherent PSK. The functional model of passband data transmission system is. Signal transmission encoder. x Signal. decoder.

Example 11: The man shown in Figure (a) pulls on the cord with a force of 70

EECE 301 Signals & Systems Prof. Mark Fowler

PHYS PRACTICE EXAM 2

THE LAPLACE EQUATION. The Laplace (or potential) equation is the equation. u = 0. = 2 x 2. x y 2 in R 2

Two-Pion Exchange Currents in Photodisintegration of the Deuteron

General Non-Arbitrage Model. I. Partial Differential Equation for Pricing A. Traded Underlying Security

Chapter. s r. check whether your calculator is in all other parts of the body. When a rigid body rotates through a given angle, all

Chapter 31 Faraday s Law

Three dimensional flow analysis in Axial Flow Compressors

ÖRNEK 1: THE LINEAR IMPULSE-MOMENTUM RELATION Calculate the linear momentum of a particle of mass m=10 kg which has a. kg m s

Lecture 4 ( ) Some points of vertical motion: Here we assumed t 0 =0 and the y axis to be vertical.

Total Deformation and its Role in Heavy Precipitation Events Associated with Deformation-Dominant Flow Patterns

AST1100 Lecture Notes

21.9 Magnetic Materials

ON 3-DIMENSIONAL CONTACT METRIC MANIFOLDS

Electric Charge. Electric charge is quantized. Electric charge is conserved

Physics 207 Lecture 13

AP Physics 1 MC Practice Kinematics 1D

Motion Along a Straight Line

MEEN 617 Handout #11 MODAL ANALYSIS OF MDOF Systems with VISCOUS DAMPING

Circular Motion. Radians. One revolution is equivalent to which is also equivalent to 2π radians. Therefore we can.

Chapter 13: Gravitation

A) N B) 0.0 N C) N D) N E) N

ENGI 4430 Advanced Calculus for Engineering Faculty of Engineering and Applied Science Problem Set 9 Solutions [Theorems of Gauss and Stokes]

Combinatorial Approach to M/M/1 Queues. Using Hypergeometric Functions

ANTENNAS. Vector and Scalar Potentials. Maxwell's Equations. D = εe. For a linear, homogeneous, isotropic medium µ and ε are contant.

Today - Lecture 13. Today s lecture continue with rotations, torque, Note that chapters 11, 12, 13 all involve rotations

Qualifying Examination Electricity and Magnetism Solutions January 12, 2006

The shortest path between two truths in the real domain passes through the complex domain. J. Hadamard

10.7 Temperature-dependent Viscoelastic Materials

Field emission of Electrons from Negatively Charged Cylindrical Particles with Nonlinear Screening in a Dusty Plasma

Note: Please use the actual date you accessed this material in your citation.

Central Force Problem. Central Force Motion. Two Body Problem: Center of Mass Coordinates. Reduction of Two Body Problem 8.01 W14D1. + m 2. m 2.

CHAPTER 5: Circular Motion; Gravitation

Transcription:

MEA GRAVITY ALG PLUMBLIE Beh-Anne Main 1, Chis MacPhee, Rbe Tenze 1, Pe Vaníek 1 and Macel Sans 1 1. Inducin 1 Univesiy f ew Bunswick, Depamen f Gedesy and Gemaics Engineeing, Fedeicn,.B., E3B 5A3, Canada Advanced Cmpuainal Reseach Labay, Faculy f Cmpue Science, Univesiy f ew Bunswick, Fedeicn,.B., E3B 5A3, Canada The mean value f gaviy alng he plumbline beween he geid and he eah suface depends n he mass densiy disibuin wihin he Eah and he shape f he Eah. Since he acual values f gaviy alng he plumbline cann be measued, he mean value f gaviy alng he plumbline has be cmpued fm he gaviy bseved a he suface f he Eah. This can be dne by educing he bseved gaviy accding sme acceped physical mdel.. Mean gaviy alng plumbline Accding he heem f mean inegal values, he mean value g f gaviy alng he plumbline beween he geid and he eah suface eads (eiskanen and Miz, 1967, Eq. 4-0) 1 g g g, Ω d, Ω = g (1) g is he gaviy a a pin f which gecenic psiin is Ω ) φ, λ), and he gecenic Ω f he eah suface is given (wih an accuacy f a ms a few millimees) by he whee adius gecenic adius f he geid g plus he hmeic heigh. The gaviy g can be expessed as (Tenze and Vaníek, 003) g φ ) δg V γ, () whee γ (,φ ) is he nmal gaviy f he gecenic efeence ellipsid, g δ is he gaviy disubance in he Tpgaphy gaviy space (i.e., he gaviy disubance geneaed by he mass wihin he geid iself, Vaníek e al., 003), and V / is he gaviainal aacin f pgaphical masses. g Ω akes he fllwing fm Subsiuing Eq. () back in Eq. (1), he mean gaviy g 1 g V φ ) δg = g 3. Mean nmal gaviy alng plumbline γ d. (3) γ f nmal gaviy alng he plumbline beween he geid and he eah suface 1 g γ γ φ ) d. (4) Ω = g The fis em n he igh-hand-side f Eq. (3) defines he mean value 1

eglecing he deflecin f plumbline and he cecin caused by diffeen lenghs f he plumbline (beween he geid and he eah suface) and he ellipsidal nmal (beween he gecenic efeence ellipsid and he elluid), he mean nmal gaviy alng he plumbline can be ewien as (Tenze and Vaníek, 003) γ 1, (5) = ( φ ) n ( φ ) γ φ ) γ, φ d whee ( φ ) is he gecenic adius f he gecenic efeence ellipsid, (Mldensky, 1945), and is he geidal heigh. = ( φ ) is he nmal heigh The mean value f nmal gaviy alng he nmal beween he efeence (gecenic) ellipsid and he elluid is evaluaed by he fllwing fmula (eiskanen and Miz, 1967, Eq. 4-4) 1 Ω a b Ω Ω φ φ ) ( φ ) γ = γ d γ 1 1 f f sin ϕ, (6) = ( φ ) GM a a whee a, b ae he semi-axes and f = ( a - b) / a is he fis numeical flaening f he gecenic efeence ellipsid, is he mean angula velciy f Eah s ain, ϕ is he gedeic laiude, GM is he gecenic gaviainal cnsan, and γ ( φ ) is he nmal gaviy a he suface f he gecenic efeence ellipsid. Fig. 1: Relain beween he mean value γ f nmal gaviy (efeence ellipsid GRS-80) and he nmal heigh.

55 981.4 [gal] 981.3 [gal] 981. [gal] 5 981.1 [gal] 981.0 [gal] 980.9 [gal] 980.8 [gal] 35 36 37 38 39 980.7 [gal] The cecin nmal gaviy Fig. : Mean values f nmal gaviy. ε due he geid undulain, i.e., cecin f a shif f he γ inegain ineval fm he efeence ellipsid he geid, is given by (Tenze and Vaníek, 003) γ ( φ ) ( φ ) γ, γ φ ε γ γ g =. (7) n a = ( φ ) Fig. 3: Relain beween he cecin γ and he geidal heigh. ε nmal gaviy due he geid undulain 3

55 5. [mgal] 5.1 [mgal] 5 35 36 37 38 39 Fig. 4: Cecin nmal gaviy due he geid undulain. 5.0 [mgal] 4.9 [mgal] 4.8 [mgal] 4.7 [mgal] 4.6 [mgal] 4.5 [mgal] 4.4 [mgal] 4.3 [mgal] 4. [mgal] 4.1 [mgal] 4.0 [mgal] 3.9 [mgal] 3.8 [mgal] 3.7 [mgal] 3.6 [mgal] 4. Mean geid-geneaed gaviy disubance The mean value g δ f geid-geneaed gaviy disubance alng he plumbline beween he geid and he eah suface is given by he secnd em n he igh-hand-side f Eq. (3), i.e., R g = R Since he geid-geneaed gaviy disubance g 1 δ δ d. (8) g δ muliplied by a gecenic adial disance is hamnic abve he geid in he Tpgaphy gaviy space, i.e., saisfies he Laplace equain, he gaviy disubance δ g is evaluaed by slving Diichle s bunday value pblem (Kellgg, 199) 1 R δ g = K[, ψ ( Ω, R] δg ( R, Ω ) dω, (9) 4 Ω Ω whee he gaviy disubance δg ( R, δg g is efeed n he c-geid, R is he mean adius f he Eah which appximaes he gecenic adius f he geid ( c-geid) suface, and Pissn s K, ψ Ω, Ω, R is given by inegal kenel [ ] R K[, ψ ( Ω, R ] = R. (10) 3 l [, ψ ( Ω, R] Subsiuing Eqns. (10) and (9) back in Eq. (8) and pefming he adial inegain wih espec, Eq. (8) akes he fllwing fm (Tenze and Vaníek, 003) 4

δg = 1 4 R l[, ψ ( Ω, R] Ω Ω R agsinh R sinψ ( Ω, Ω ) anψ ( Ω, Ω ) 1 R 1 = R δg ( R, Ω ) dω. (11) 55 10 [mgal] 100 [mgal] 80 [mgal] 60 [mgal] 40 [mgal] 0 [mgal] 5 0 [mgal] -0 [mgal] -40 [mgal] -60 [mgal] -80 [mgal] -100 [mgal] -10 [mgal] -140 [mgal] 35 36 37 38 39 Fig. 5: Mean geid-geneaed gaviy disubances g δ. -160 [mgal] 5. Mean pgaphy-geneaed gaviainal aacin The mean value f pgaphy-geneaed gaviainal aacin alng he plumbline beween he geid and he eah suface (given by he hid em n he igh-hand-side f Eq. 3) can be deived as (Tenze and Vaníek, 003) whee V ( R, and ( ) 1 g V 1 R V d = = Ω g = R V ( R, V ( ) = d, (1) V ae he gaviainal penials f pgaphical masses as a efeed n he geid and a he eah suface. V,Ω f he pgaphical masses f a pin Accding Mainec (1998), he gaviainal penial inside he pgaphical mass R R eads 5

V ( = G R R 3 R 1 [ ], 3 3 G Ω Ω G Ω Ω R = R R = R ( Ω ) ( Ω ) l 1 [, ψ ( Ω, ] 1 (, Ω ) l [, ψ ( Ω, ] d dω δρ d dω, (13) whee G is ewn s gaviainal cnsan, l [,ψ ( Ω, ] is he spaial disance beween (, (,Ω ), and ψ ( Ω,Ω ) is he spheical disance beween he gecenic diecin Ω and Ω. and The fis em n he igh-hand-side f Eq. (13) is he gaviainal penial f he spheical Bugue shell (f he mean pgaphical densiy and hickness f ) inside he pgaphical masses (Wichienchaen, 198). The secnd em sands f he gaviainal penial f he eain ughness em f densiy (Mainec, 1998), and he hid em epesens he effec f anmalus pgaphical densiy δρ disibuin n he gaviainal penial. Subsiuing Eq. (13) back in Eq. (1), he mean value f pgaphy-geneaed gaviainal aacin alng he adial diecin (which appximaes he plumbline) becmes (Tenze and Vaníek, 003) 1 R = R G V G Ω Ω Ω Ω R d = = R R = R ( Ω ) ( Ω ) 3R G Ω Ω 1 1 ( l [ R, ψ ( Ω, ] l [, ψ ( Ω, ] ) 3 1 1 (, Ω ) l [ R, ψ ( Ω, ] l [, ψ ( Ω, ] d dω δρ d dω. (14) b Fig. 6: Relain beween he mean gaviy hmeic heigh g geneaed by he spheical Bugue shell and he. 6

55 40 [mgal] 0 [mgal] 00 [mgal] 180 [mgal] 160 [mgal] 5 140 [mgal] 10 [mgal] 100 [mgal] 80 [mgal] 60 [mgal] 35 36 37 38 39 40 [mgal] 0 [mgal] 0 [mgal] Fig. 7: Mean values f he gaviainal aacin caused by he spheical Bugue shell f densiy. 55 30 [mgal] 0 [mgal] 5 10 [mgal] 0 [mgal] -10 [mgal] -0 [mgal] -30 [mgal] -40 [mgal] 35 36 37 38 39 - [mgal] Fig. 8: Mean values f he gaviainal aacin caused by he eain ughness em. 7

δρ Fig. 9: Relain beween he mean gaviainal aacin pgaphical densiy disibuin and he hmeic heigh g f he anmalus (laeal). 55 5 [mgal] 5 0 [mgal] 15 [mgal] 10 [mgal] 5 [mgal] 0 [mgal] -5 [mgal] 35 36 37 38 39 6. Cnclusins Fig. 10: Mean values f he gaviainal aacin caused by he anmalus (laeally vaying) pgaphical densiy disibuin. -10 [mgal] The mean value f gaviy (geneaed by he slid Eah wihu he amsphee) alng he plumbline beween he geid and he eah suface (Fig. 11) is descibed as a sum f he mean nmal gaviy beween he ellipsid and he elluid (Eq. 6), he cecin nmal gaviy f geid undulain (Eq. 7), he mean geid-geneaed gaviy disubance (Eq. 11), and he mean pgaphy-geneaed gaviainal aacin (Eq. 14). The mean pgaphy-geneaed gaviainal aacin cnsiss f he 8

mean gaviainal aacins caused by he spheical Bugue shell f densiy (fis em in Eq. 14), he eain ughness em (secnd em in Eq. 14), and anmalus pgaphical densiy disibuin (hid em in Eq. 14). Minimum, maximum, and aveage values f hmeic heighs (Fig. 1), geidal heighs (Fig. 13) and he laeal vaiain f he pgaphical densiy (Fig. 14) a he esing aea ϕ, 55, λ 35, 39 in pa f he Canadian Rcky Munains ae in Tab. 1. The esuls f he numeical invesigain f all cmpnens he mean gaviy ae summaized in Tab.. 55 981.5 [gal] 981.5 [gal] 981.4 [gal] 981.3 [gal] 981.3 [gal] 5 981.3 [gal] 981. [gal] 981. [gal] 981.1 [gal] 981.0 [gal] 981.0 [gal] 981.0 [gal] 980.9 [gal] 35 36 37 38 39 980.8 [gal] Fig. 11: Mean values f gaviy alng he plumbline beween he geid and he eah suface a a pa f he Canadian Rcky Munains. Tab. 1: hmeic eighs [m] (Fig. 1) Geidal eighs [m] (Fig. 13) Anmalus Laeal Densiy [g.cm -3 ] (Fig. 14) Min. Max. Aveage 4 736 1166-17.17-11.90-14.63-0.18 0.31 0.04 9

Tab. : Mean mal Gaviy (Fig. ) Cecin mal Gaviy f Geid Undulain (Fig. 4) Mean Geid-Geneaed Gaviy Disubance (Fig. 5) Mean Gaviain Aacin f Spheical Bugue Shell (Fig. 7) Mean Gaviainal Aacin f Teain Rughness Tem (Fig. 8) Mean Gaviainal Aacin f Anmalus Tpgaphical Densiy (Fig. 10) Mean Gaviy (Fig. 11) Min. [mgal] Max. [mgal] Aveage [mgal] 9807.5 98139.1 981111.0 3.7 5.3 4.6-169. 18. -14.6 0.0 86.1 130.4-64.3 35.0 -.1-10.1 9.0.3 980858.5 9817.1 9813. 55 400 [m] 00 [m] 000 [m] 1800 [m] 5 1600 [m] 1400 [m] 100 [m] 1000 [m] 800 [m] 600 [m] 400 [m] 00 [m] 35 36 37 38 39 Fig. 1: Teain a a pa f he Canadian Rcky Munains. 0 [m] 10

55-1 [m] -13 [m] 5-14 [m] -15 [m] -16 [m] -17 [m] 35 36 37 38 39 Fig. 13: Geid a a pa f he Canadian Rcky Munains. -18 [m] 55 0.3 [g.cm-3] 5 0. [g.cm-3] 0.1 [g.cm-3] 0.0 [g.cm-3] -0.1 [g.cm-3] 35 36 37 38 39-0. [g.cm-3] Fig: 14: Laeal vaiain f pgaphical densiy a a pa f he Canadian Rcky Munains. 11

Refeences: eiskanen W. A., Miz., 1967: Physical gedesy. W.. Feeman and C., San Fancisc. Kellgg.D., 199: Fundains f penial hey. Spinge. Belin. Mainec Z., 1998: Bunday value pblems f gavimeic deeminain f a pecise. Lecue nes in eah sciences, Vl. 73, Spinge. Mldensky M.S., 1945: Fundamenal pblems f Gedeic Gavimey (in Russian). TRUDY Ts IIGAIK 4, Gedezizda, Mscw. Tenze R., Vaníek P., 003: Rigus definiin f mean gaviy in hey f heigh. Junal f Gedesy, Spinge. (in pepaain) Vaníek P., Tenze R., Sjöbeg L.E., Mainec Z., Feahesne W.E., 003: ew views f he spheical Bugue gaviy anmaly. Gephysical Junal Inenainal. (submied) Wichienchaen C., 198: The indiec effecs n he cmpuain f geid undulains. Dep. f Ged. Sci. Rep.336, hi Sae Univesiy. 1