- Source: Theoretical gravity
In geodesy and geophysics, theoretical gravity or normal gravity is an approximation of Earth's gravity, on or near its surface, by means of a mathematical model. The most common theoretical model is a rotating Earth ellipsoid of revolution (i.e., a spheroid).
Other representations of gravity can be used in the study and analysis of other bodies, such as asteroids. Widely used representations of a gravity field in the context of geodesy include spherical harmonics, mascon models, and polyhedral gravity representations.
Principles
The type of gravity model used for the Earth depends upon the degree of fidelity required for a given problem. For many problems such as aircraft simulation, it may be sufficient to consider gravity to be a constant, defined as:
g
=
g
45
=
{\displaystyle g=g_{45}=}
9.80665 m/s2 (32.1740 ft/s2)
based upon data from World Geodetic System 1984 (WGS-84), where
g
{\displaystyle g}
is understood to be pointing 'down' in the local frame of reference.
If it is desirable to model an object's weight on Earth as a function of latitude, one could use the following:: 41
g
=
g
45
−
1
2
(
g
p
o
l
e
s
−
g
e
q
u
a
t
o
r
)
cos
(
2
φ
⋅
π
180
)
{\displaystyle g=g_{45}-{\tfrac {1}{2}}(g_{\mathrm {poles} }-g_{\mathrm {equator} })\cos \left(2\,\varphi \cdot {\frac {\pi }{180}}\right)}
where
g
p
o
l
e
s
{\displaystyle g_{\mathrm {poles} }}
= 9.832 m/s2 (32.26 ft/s2)
g
45
{\displaystyle g_{45}}
= 9.806 m/s2 (32.17 ft/s2)
g
e
q
u
a
t
o
r
{\displaystyle g_{\mathrm {equator} }}
= 9.780 m/s2 (32.09 ft/s2)
φ
{\displaystyle \varphi }
= latitude, between −90° and +90°
Neither of these accounts for changes in gravity with changes in altitude, but the model with the cosine function does take into account the centrifugal relief that is produced by the rotation of the Earth. On the rotating sphere, the sum of the force of the gravitational field and the centrifugal force yields an angular deviation of approximately
sin
(
2
φ
)
2
g
R
Ω
2
{\displaystyle {\frac {\sin(2\varphi )}{2g}}{R\Omega ^{2}}}
(in radians) between the direction of the gravitational field and the direction measured by a plumb line; the plumb line appears to point southwards on the northern hemisphere and northwards on the southern hemisphere.
Ω
≈
7.29
×
10
−
5
{\displaystyle \Omega \approx 7.29\times 10^{-5}}
rad/s is the diurnal angular speed of the Earth axis, and
R
≈
6370
{\displaystyle R\approx 6370}
km the radius of the reference sphere, and
R
sin
φ
{\displaystyle R\sin \varphi }
the distance of the point on the Earth crust to the Earth axis.
For the mass attraction effect by itself, the gravitational acceleration at the equator is about 0.18% less than that at the poles due to being located farther from the mass center. When the rotational component is included (as above), the gravity at the equator is about 0.53% less than that at the poles, with gravity at the poles being unaffected by the rotation. So the rotational component of change due to latitude (0.35%) is about twice as significant as the mass attraction change due to latitude (0.18%), but both reduce strength of gravity at the equator as compared to gravity at the poles.
Note that for satellites, orbits are decoupled from the rotation of the Earth so the orbital period is not necessarily one day, but also that errors can accumulate over multiple orbits so that accuracy is important. For such problems, the rotation of the Earth would be immaterial unless variations with longitude are modeled. Also, the variation in gravity with altitude becomes important, especially for highly elliptical orbits.
The Earth Gravitational Model 1996 (EGM96) contains 130,676 coefficients that refine the model of the Earth's gravitational field.: 40 The most significant correction term is about two orders of magnitude more significant than the next largest term.: 40 That coefficient is referred to as the
J
2
{\displaystyle J_{2}}
term, and accounts for the flattening of the poles, or the oblateness, of the Earth. (A shape elongated on its axis of symmetry, like an American football, would be called prolate.) A gravitational potential function can be written for the change in potential energy for a unit mass that is brought from infinity into proximity to the Earth. Taking partial derivatives of that function with respect to a coordinate system will then resolve the directional components of the gravitational acceleration vector, as a function of location. The component due to the Earth's rotation can then be included, if appropriate, based on a sidereal day relative to the stars (≈366.24 days/year) rather than on a solar day (≈365.24 days/year). That component is perpendicular to the axis of rotation rather than to the surface of the Earth.
A similar model adjusted for the geometry and gravitational field for Mars can be found in publication NASA SP-8010.
The barycentric gravitational acceleration at a point in space is given by:
g
=
−
G
M
r
2
r
^
{\displaystyle \mathbf {g} =-{GM \over r^{2}}\mathbf {\hat {r}} }
where:
M is the mass of the attracting object,
r
^
{\displaystyle \scriptstyle \mathbf {\hat {r}} }
is the unit vector from center-of-mass of the attracting object to the center-of-mass of the object being accelerated, r is the distance between the two objects, and G is the gravitational constant.
When this calculation is done for objects on the surface of the Earth, or aircraft that rotate with the Earth, one has to account for the fact that the Earth is rotating and the centrifugal acceleration has to be subtracted from this. For example, the equation above gives the acceleration at 9.820 m/s2, when GM = 3.986 × 1014 m3/s2, and R = 6.371 × 106 m. The centripetal radius is r = R cos(φ), and the centripetal time unit is approximately (day / 2π), reduces this, for r = 5 × 106 metres, to 9.79379 m/s2, which is closer to the observed value.
Basic formulas
Various, successively more refined, formulas for computing the theoretical gravity are referred to as the International Gravity Formula, the first of which was proposed in 1930 by the International Association of Geodesy. The general shape of that formula is:
g
(
ϕ
)
=
g
e
(
1
+
A
sin
2
(
ϕ
)
−
B
sin
2
(
2
ϕ
)
)
,
{\displaystyle g(\phi )=g_{e}\left(1+A\sin ^{2}(\phi )-B\sin ^{2}(2\phi )\right),}
in which g(φ) is the gravity as a function of the geographic latitude φ of the position whose gravity is to be determined,
g
e
{\displaystyle g_{e}}
denotes the gravity at the equator (as determined by measurement), and the coefficients A and B are parameters that must be selected to produce a good global fit to true gravity.
Using the values of the GRS80 reference system, a commonly used specific instantiation of the formula above is given by:
g
(
ϕ
)
=
9.780327
(
1
+
0.0053024
sin
2
(
ϕ
)
−
0.0000058
sin
2
(
2
ϕ
)
)
m
s
−
2
.
{\displaystyle g(\phi )=9.780327\left(1+0.0053024\sin ^{2}(\phi )-0.0000058\sin ^{2}(2\phi )\right)\,\mathrm {ms} ^{-2}.}
Using the appropriate double-angle formula in combination with the Pythagorean identity, this can be rewritten in the equivalent forms
g
(
ϕ
)
=
9.780327
(
1
+
0.0052792
sin
2
(
ϕ
)
+
0.0000232
sin
4
(
ϕ
)
)
m
s
−
2
,
=
9.780327
(
1.0053024
−
.0053256
cos
2
(
ϕ
)
+
.0000232
cos
4
(
ϕ
)
)
m
s
−
2
,
=
9.780327
(
1.0026454
−
0.0026512
cos
(
2
ϕ
)
+
.0000058
cos
2
(
2
ϕ
)
)
m
s
−
2
.
{\displaystyle {\begin{aligned}g(\phi )&=9.780327\left(1+0.0052792\sin ^{2}(\phi )+0.0000232\sin ^{4}(\phi )\right)\,\mathrm {ms} ^{-2},\\&=9.780327\left(1.0053024-.0053256\cos ^{2}(\phi )+.0000232\cos ^{4}(\phi )\right)\,\mathrm {ms} ^{-2},\\&=9.780327\left(1.0026454-0.0026512\cos(2\phi )+.0000058\cos ^{2}(2\phi )\right)\,\mathrm {ms} ^{-2}.\end{aligned}}\,\!}
Up to the 1960s, formulas based on the Hayford ellipsoid (1924) and of the famous German geodesist Helmert (1906) were often used. The difference between the semi-major axis (equatorial radius) of the Hayford ellipsoid and that of the modern WGS84 ellipsoid is 251 m; for Helmert's ellipsoid it is only 63 m.
Somigliana equation
A more recent theoretical formula for gravity as a function of latitude is the International Gravity Formula 1980 (IGF80), also based on the GRS80 ellipsoid but now using the Somigliana equation (after Carlo Somigliana (1860–1955)):
g
(
ϕ
)
=
g
e
[
1
+
k
sin
2
(
ϕ
)
1
−
e
2
sin
2
(
ϕ
)
]
,
{\displaystyle g(\phi )=g_{e}\left[{\frac {1+k\sin ^{2}(\phi )}{\sqrt {1-e^{2}\sin ^{2}(\phi )}}}\right],\,\!}
where,
k
=
b
g
p
−
a
g
e
a
g
e
{\displaystyle k={\frac {bg_{p}-ag_{e}}{ag_{e}}}}
(formula constant);
g
e
,
g
p
{\displaystyle g_{e},g_{p}}
is the defined gravity at the equator and poles, respectively;
a
,
b
{\displaystyle a,b}
are the equatorial and polar semi-axes, respectively;
e
2
=
a
2
−
b
2
a
2
{\displaystyle e^{2}={\frac {a^{2}-b^{2}}{a^{2}}}}
is the spheroid's squared eccentricity;
providing,
g
(
ϕ
)
=
9.7803267715
[
1
+
0.001931851353
sin
2
(
ϕ
)
1
−
0.0066943800229
sin
2
(
ϕ
)
]
m
s
−
2
.
{\displaystyle g(\phi )=9.7803267715\left[{\frac {1+0.001931851353\sin ^{2}(\phi )}{\sqrt {1-0.0066943800229\sin ^{2}(\phi )}}}\right]\,\mathrm {ms} ^{-2}.}
A later refinement, based on the WGS84 ellipsoid, is the WGS (World Geodetic System) 1984 Ellipsoidal Gravity Formula:
g
(
ϕ
)
=
9.7803253359
[
1
+
0.00193185265241
sin
2
(
ϕ
)
1
−
0.00669437999013
sin
2
(
ϕ
)
]
m
s
−
2
.
{\displaystyle g(\phi )=9.7803253359\left[{\frac {1+0.00193185265241\sin ^{2}(\phi )}{\sqrt {1-0.00669437999013\sin ^{2}(\phi )}}}\right]\,\mathrm {ms} ^{-2}.}
(where
g
p
{\displaystyle g_{p}}
= 9.8321849378 ms−2)
The difference with IGF80 is insignificant when used for geophysical purposes, but may be significant for other uses.
= Further details
=
For the normal gravity
γ
0
{\displaystyle \gamma _{0}}
of the sea level ellipsoid, i.e., elevation h = 0, this formula by Somigliana (1929) applies:
γ
0
(
φ
)
=
a
⋅
γ
a
⋅
cos
2
φ
+
b
⋅
γ
b
⋅
sin
2
φ
a
2
⋅
cos
2
φ
+
b
2
⋅
sin
2
φ
{\displaystyle \gamma _{0}(\varphi )={\frac {a\cdot \gamma _{a}\cdot \cos ^{2}\varphi +b\cdot \gamma _{b}\cdot \sin ^{2}\varphi }{\sqrt {a^{2}\cdot \cos ^{2}\varphi +b^{2}\cdot \sin ^{2}\varphi }}}}
with
γ
a
{\displaystyle \gamma _{a}}
= Normal gravity at Equator
γ
b
{\displaystyle \gamma _{b}}
= Normal gravity at poles
a = semi-major axis (Equator radius)
b = semi-minor axis (pole radius)
φ
{\displaystyle \varphi }
= latitude
Due to numerical issues, the formula is simplified to this:
γ
0
(
φ
)
=
γ
a
⋅
1
+
p
⋅
sin
2
φ
1
−
e
2
⋅
sin
2
φ
{\displaystyle \gamma _{0}(\varphi )=\gamma _{a}\cdot {\frac {1+p\cdot \sin ^{2}\varphi }{\sqrt {1-e^{2}\cdot \sin ^{2}\varphi }}}}
with
p
=
b
⋅
γ
b
a
⋅
γ
a
−
1
{\displaystyle p={\frac {b\cdot \gamma _{b}}{a\cdot \gamma _{a}}}-1}
e
2
=
1
−
b
2
a
2
;
{\displaystyle e^{2}=1-{\frac {b^{2}}{a^{2}}};\quad }
(e is the eccentricity)
For the Geodetic Reference System 1980 (GRS 80) the parameters are set to these values:
a
=
6
378
137
m
b
=
6
356
752
.
314
1
m
{\displaystyle a=6\,378\,137\,\mathrm {m} \quad \quad \quad \quad b=6\,356\,752{.}314\,1\,\mathrm {m} }
γ
a
=
9
.
780
326
771
5
m
s
2
γ
b
=
9
.
832
186
368
5
m
s
2
{\displaystyle \gamma _{a}=9{.}780\,326\,771\,5\,\mathrm {\frac {m}{s^{2}}} \quad \gamma _{b}=9{.}832\,186\,368\,5\,\mathrm {\frac {m}{s^{2}}} }
⇒
p
=
1
.
931
851
353
⋅
10
−
3
e
2
=
6
.
694
380
022
90
⋅
10
−
3
{\displaystyle \Rightarrow p=1{.}931\,851\,353\cdot 10^{-3}\quad e^{2}=6{.}694\,380\,022\,90\cdot 10^{-3}}
Approximation formula from series expansions
The Somigliana formula was approximated through different series expansions, following this scheme:
γ
0
(
φ
)
=
γ
a
⋅
(
1
+
β
⋅
sin
2
φ
+
β
1
⋅
sin
2
2
φ
+
…
)
{\displaystyle \gamma _{0}(\varphi )=\gamma _{a}\cdot (1+\beta \cdot \sin ^{2}\varphi +\beta _{1}\cdot \sin ^{2}2\varphi +\dots )}
= International gravity formula 1930
=
The normal gravity formula by Gino Cassinis was determined in 1930 by International Union of Geodesy and Geophysics as international gravity formula along with Hayford ellipsoid. The parameters are:
γ
a
=
9
.
78049
m
s
2
β
=
5
.
2884
⋅
10
−
3
β
1
=
−
5
.
9
⋅
10
−
6
{\displaystyle \gamma _{a}=9{.}78049{\frac {\mathrm {m} }{\mathrm {s} ^{2}}}\quad \beta =5{.}2884\cdot 10^{-3}\quad \beta _{1}=-5{.}9\cdot 10^{-6}}
In the course of time the values were improved again with newer knowledge and more exact measurement methods.
Harold Jeffreys improved the values in 1948 at:
γ
a
=
9
.
780373
m
s
2
β
=
5
.
2891
⋅
10
−
3
β
1
=
−
5
.
9
⋅
10
−
6
{\displaystyle \gamma _{a}=9{.}780373{\frac {\mathrm {m} }{\mathrm {s} ^{2}}}\quad \beta =5{.}2891\cdot 10^{-3}\quad \beta _{1}=-5{.}9\cdot 10^{-6}}
= International gravity formula 1967
=
The normal gravity formula of Geodetic Reference System 1967 is defined with the values:
γ
a
=
9
.
780318
m
s
2
β
=
5
.
3024
⋅
10
−
3
β
1
=
−
5
.
9
⋅
10
−
6
{\displaystyle \gamma _{a}=9{.}780318{\frac {\mathrm {m} }{\mathrm {s} ^{2}}}\quad \beta =5{.}3024\cdot 10^{-3}\quad \beta _{1}=-5{.}9\cdot 10^{-6}}
= International gravity formula 1980
=
From the parameters of GRS 80 comes the classic series expansion:
γ
a
=
9
.
780327
m
s
2
β
=
5
.
3024
⋅
10
−
3
β
1
=
−
5
.
8
⋅
10
−
6
{\displaystyle \gamma _{a}=9{.}780327{\frac {\mathrm {m} }{\mathrm {s} ^{2}}}\quad \beta =5{.}3024\cdot 10^{-3}\quad \beta _{1}=-5{.}8\cdot 10^{-6}}
The accuracy is about ±10−6 m/s2.
With GRS 80 the following series expansion is also introduced:
γ
0
(
φ
)
=
γ
a
⋅
(
1
+
c
1
⋅
sin
2
φ
+
c
2
⋅
sin
4
φ
+
c
3
⋅
sin
6
φ
+
c
4
⋅
sin
8
φ
+
…
)
{\displaystyle \gamma _{0}(\varphi )=\gamma _{a}\cdot (1+c_{1}\cdot \sin ^{2}\varphi +c_{2}\cdot \sin ^{4}\varphi +c_{3}\cdot \sin ^{6}\varphi +c_{4}\cdot \sin ^{8}\varphi +\dots )}
As such the parameters are:
c1 = 5.279 0414·10−3
c2 = 2.327 18·10−5
c3 = 1.262·10−7
c4 = 7·10−10
The accuracy is at about ±10−9 m/s2 exact. When the exactness is not required, the terms at further back can be omitted. But it is recommended to use this finalized formula.
Height dependence
Cassinis determined the height dependence, as:
g
(
φ
,
h
)
=
g
0
(
φ
)
−
(
3
.
08
⋅
10
−
6
1
s
2
−
4
.
19
⋅
10
−
7
c
m
3
g
⋅
s
2
⋅
ρ
)
⋅
h
{\displaystyle g(\varphi ,h)=g_{0}(\varphi )-\left(3{.}08\cdot 10^{-6}\,{\frac {1}{\mathrm {s} ^{2}}}-4{.}19\cdot 10^{-7}\,{\frac {\mathrm {cm} ^{3}}{\mathrm {g} \cdot \mathrm {s} ^{2}}}\cdot \rho \right)\cdot h}
The average rock density ρ is no longer considered.
Since GRS 1967 the dependence on the ellipsoidal elevation h is:
g
(
φ
,
h
)
=
g
0
(
φ
)
−
(
1
−
1
.
39
⋅
10
−
3
⋅
sin
2
(
φ
)
)
⋅
3
.
0877
⋅
10
−
6
1
s
2
⋅
h
+
7
.
2
⋅
10
−
13
1
m
⋅
s
2
⋅
h
2
=
g
0
(
φ
)
−
(
3
.
0877
⋅
10
−
6
−
4
.
3
⋅
10
−
9
⋅
sin
2
(
φ
)
)
1
s
2
⋅
h
+
7
.
2
⋅
10
−
13
1
m
⋅
s
2
⋅
h
2
{\displaystyle {\begin{aligned}g(\varphi ,h)&=g_{0}(\varphi )-\left(1-1{.}39\cdot 10^{-3}\cdot \sin ^{2}(\varphi )\right)\cdot 3{.}0877\cdot 10^{-6}\,{\frac {1}{\mathrm {s} ^{2}}}\cdot h+7{.}2\cdot 10^{-13}\,{\frac {1}{\mathrm {m} \cdot \mathrm {s} ^{2}}}\cdot h^{2}\\&=g_{0}(\varphi )-\left(3{.}0877\cdot 10^{-6}-4{.}3\cdot 10^{-9}\cdot \sin ^{2}(\varphi )\right)\,{\frac {1}{\mathrm {s} ^{2}}}\cdot h+7{.}2\cdot 10^{-13}\,{\frac {1}{\mathrm {m} \cdot \mathrm {s} ^{2}}}\cdot h^{2}\end{aligned}}}
Another expression is:
g
(
φ
,
h
)
=
g
0
(
φ
)
⋅
(
1
−
(
k
1
−
k
2
⋅
sin
2
φ
)
⋅
h
+
k
3
⋅
h
2
)
{\displaystyle g(\varphi ,h)=g_{0}(\varphi )\cdot (1-(k_{1}-k_{2}\cdot \sin ^{2}\varphi )\cdot h+k_{3}\cdot h^{2})}
with the parameters derived from GRS80:
k
1
=
2
⋅
(
1
+
f
+
m
)
/
a
=
3
.
157
04
⋅
10
−
7
m
−
1
{\displaystyle k_{1}=2\cdot (1+f+m)/a=3{.}157\,04\cdot 10^{-7}\,\mathrm {m^{-1}} }
k
2
=
4
⋅
f
/
a
=
2
.
102
69
⋅
10
−
9
m
−
1
{\displaystyle k_{2}=4\cdot f/a=2{.}102\,69\cdot 10^{-9}\,\mathrm {m^{-1}} }
k
3
=
3
/
(
a
2
)
=
7
.
374
52
⋅
10
−
14
m
−
2
{\displaystyle k_{3}=3/(a^{2})=7{.}374\,52\cdot 10^{-14}\,\mathrm {m^{-2}} }
where
m
{\displaystyle m}
with
ω
=
7.2921150
⋅
10
−
5
r
a
d
⋅
s
−
1
{\displaystyle \omega =7.2921150\cdot 10^{-5}\ rad\cdot s^{-1}}
:
m
=
ω
2
⋅
a
2
⋅
b
G
M
{\displaystyle m={\frac {\omega ^{2}\cdot a^{2}\cdot b}{GM}}}
This adjustment is about right for common heights in aviation; but for heights up to outer space (over ca. 100 kilometers) it is out of range.
WELMEC formula
In all German standards offices the free-fall acceleration g is calculated in respect to the average latitude φ and the average height above sea level h with the WELMEC–Formel:
g
(
φ
,
h
)
=
(
1
+
0
.
0053024
⋅
sin
2
(
φ
)
−
0
.
0000058
⋅
sin
2
(
2
φ
)
)
⋅
9
.
780318
m
s
2
−
0
.
000003085
1
s
2
⋅
h
{\displaystyle g(\varphi ,h)=\left(1+0{.}0053024\cdot \sin ^{2}(\varphi )-0{.}0000058\cdot \sin ^{2}(2\varphi )\right)\cdot 9{.}780318{\frac {\mathrm {m} }{\mathrm {s} ^{2}}}-0{.}000003085\,{\frac {1}{\mathrm {s} ^{2}}}\cdot h}
The formula is based on the International gravity formula from 1967.
The scale of free-fall acceleration at a certain place must be determined with precision measurement of several mechanical magnitudes. Weighing scales, the mass of which does measurement because of the weight, relies on the free-fall acceleration, thus for use they must be prepared with different constants in different places of use. Through the concept of so-called gravity zones, which are divided with the use of normal gravity, a weighing scale can be calibrated by the manufacturer before use.
Example
Free-fall acceleration in Schweinfurt:
Data:
Latitude: 50° 3′ 24″ = 50.0567°
Height above sea level: 229.7 m
Density of the rock plates: ca. 2.6 g/cm3
Measured free-fall acceleration: g = 9.8100 ± 0.0001 m/s2
Free-fall acceleration, calculated through normal gravity formulas:
Cassinis: g = 9.81038 m/s2
Jeffreys: g = 9.81027 m/s2
WELMEC: g = 9.81004 m/s2
See also
Gravity anomaly
Reference ellipsoid
EGM96 (Earth Gravitational Model 1996)
Standard gravity : 9.806 65 m/s2
References
Further reading
Karl Ledersteger: Astronomische und physikalische Geodäsie. Handbuch der Vermessungskunde Band 5, 10. Auflage. Metzler, Stuttgart 1969
B.Hofmann-Wellenhof, Helmut Moritz: Physical Geodesy, ISBN 3-211-23584-1, Springer-Verlag Wien 2006.
Wolfgang Torge: Geodäsie. 2. Auflage. Walter de Gruyter, Berlin u.a. 2003. ISBN 3-11-017545-2
Wolfgang Torge: Geodäsie. Walter de Gruyter, Berlin u.a. 1975 ISBN 3-11-004394-7
External links
Definition des Geodetic Reference System 1980 (GRS80) (pdf, engl.; 70 kB)
Gravity Information System der Physikalisch-Technischen Bundesanstalt, engl.
Online-Berechnung der Normalschwere mit verschiedenen Normalschwereformeln
Kata Kunci Pencarian:
7.166
6.8
No More Posts Available.
No more pages to load.
