Apuntes Geodesia Fisica 2005 PDF

Apuntes de Geodesia Física Curso 2005-2006
Apuntes de Geodesia Física
Autora: Yolanda Rubio

Modificado y actualizado por Martín Rodriguez Vales
Página 1
Página 2
Índice de contenido
1.- INTRODUCCIÓN A LA GEODESIA FÍSICA............................................................................. 6
1.1.- Definición y objeto de la Geodesia Física...............................................................................6
1.2.- Diferencia entre geodesia geométrica y geodesia física:.........................................................7
1.3.- Consideraciones históricas sobre la figura de la Tierra y sus métodos de estudio..................8
2.- TEORÍA DEL POTENCIAL. EL CAMPO GRAVITATORIO TERRESTRE........................... 11
2.1.- Ley de Newton de la Gravitación Universal......................................................................... 11
2.2.- Campo de fuerza newtoniano. Intensidad de campo y potencial. ........................................ 11
2.2.1.- Campo newtoniano debido a una masa puntual:........................................................... 11
2.2.2.- Campo newtoniano debido a un sistema de masas puntuales:.......................................13
2.2.3.- Campo newtoniano debido a un cuerpo extenso:.......................................................... 14
2.3.- Definición y propiedades de la superficie equipotencial: intensidad de campo y potencial: 15
2.4.- Flujo vector intensidad de campo. Teorema de Gauss.......................................................... 16
2.4.1.- Caso de una masa puntual interior m.............................................................................17
2.4.2.- Caso de una masa exterior a la superficie:.....................................................................18
2.4.3.- Caso varias masas interiores......................................................................................... 18
2.5.- Ecuaciones de Laplace y de Poisson. Funciones armónicas................................................. 18
2.6.- Campo y potencial de las fuerzas centrífugas....................................................................... 20
2.7.- Gravedad y potencial gravitatorio terrestres..........................................................................21
2.8.- Superficies de nivel y líneas de la plomada:......................................................................... 22
2.9.- Desarrollo del potencial gravitatorio terrestre. Polinomios de Legendre..............................23
2.10.- Problemas de contorno de la teoría del potencial. Solución de la ecuación de Laplace en
armónicos esféricos........................................................................................................................27
2.10.1.- Teorema de Stockes:....................................................................................................27
2.10.2.- Primer problema de contorno de la teoría del potencial o problema de Dirichlet:......27
3.- EL POTENCIAL Y LA GRAVEDAD NORMAL. SISTEMAS GEODÉSICOS DE
REFERENCIA................................................................................................................................... 31
3.1.- La función geopotencial........................................................................................................ 31
3.2.- Esferoide de Bruns: teoría del primer orden de aproximación..............................................32
3.3.- Primera fórmula de Clairaut en aproximación de primer orden............................................32
3.3.1.- 1ª método:..................................................................................................................... 33
3.3.2.- 2ª método:...................................................................................................................... 34
3.4.- Gravedad normal en el esferoide de Bruns o elipsoide de nivel:.......................................... 35
3.5.- Esferoide de Helmert:............................................................................................................37
3.6.- El elipsoide de nivel.............................................................................................................. 37
3.7.- Teoría de Pizzetti y Somigliana. Fórmula internacional de la gravedad............................... 37
3.8.- Ventajas del uso del elipsoide como figura de referencia:.................................................... 38
3.9.- Gravedad normal en el elipsoide:.......................................................................................... 38
3.10.- Sistema Geodésico de Referencia....................................................................................... 39
3.10.1.- Sistema Geodésico de Referencia 24-30 (1924-1930):............................................... 39
3.10.2.- Sistema Geodésico de Referencia 1967 ......................................................................40
3.10.3.- Sistema Geodésico de Referencia de 1980:.................................................................41
3.10.4.- Sistema Geodésico de Referencia WGS84:.................................................................41
3.10.5.- Conversión de anomalías gravimétricas por cambio de sistema geodésico de
referencia:..................................................................................................................................41
4.- EL CAMPO DE LA GRAVEDAD ANÓMALO. DETERMINACIÓN DEL GEOIDE............. 43
Página 3
4.1.- Potencial normal y potencial perturbador..............................................................................43

4.2.- Conceptos de anomalía de gravedad, ondulación del geoide y desviación de la vertical..... 43
4.3.- Fórmula de Bruns.................................................................................................................. 45
4.4.- Fórmula o ecuación fundamental de la geodesia física:........................................................ 45
4.5.- Aproximación esférica de la ecuación Fundamental de la Geodesia Física:.........................47
4.6.- Anomalía de la gravedad fuera del geoide:........................................................................... 48
4.7.- Fórmula de Stokes................................................................................................................. 49
4.8.- Evaluación práctica de la fórmula integral de Stokes........................................................... 49
4.9.- Componentes de la desviación de la vertical. Fórmulas de Vening- Meinesz...................... 51
5.- REDUCCIÓN DE LOS VALORES OBSERVADOS DE LA GRAVEDAD............................. 54
5.1.- Objeto de las reducciones. Tipos de reducciones.................................................................. 54
5.2.- Reducción aire libre o Faye...................................................................................................54
5.3.- Reducción Bouguer............................................................................................................... 55
5.3.1.- Reducción por lámina o placa Bouguer:........................................................................56
5.3.2.- Corrección topográfica:................................................................................................. 56
5.3.3.- Reducción Bouguer completa:.......................................................................................58
5.4.- Consideraciones generales sobre la isostasia........................................................................ 59
5.5.- El sistema isostático Pratt-Hayford....................................................................................... 60
5.6.- El sistema isostático Airy-Heiskanen....................................................................................61
5.7.- Compensación regional de Vening-Meinesz o Regional...................................................... 62
5.8.- Reducciones isostáticas. Anomalías isostáticas.................................................................... 63
5.8.1.- Sistema de Airy-Heiskanen:.......................................................................................... 64
6.- MEDIDAS DE GRAVEDAD.......................................................................................................65
6.1.- Las medidas de gravedad. Métodos para su determinación.................................................. 65
6.2.- Clasificación según los tipos de medida:..............................................................................65
6.3.- Clasificación según los métodos de medición:......................................................................65
6.4.- Métodos absolutos y relativos:.............................................................................................. 65
6.5.- Péndulo matemático o simple:...............................................................................................66
6.6.- Péndulo físico o compuesto:..................................................................................................67
6.7.- Péndulo reversible................................................................................................................. 68
6.8.- Métodos basados en el movimiento libre de graves:.............................................................69
6.8.1.- Caída libre:.....................................................................................................................69
6.8.2.- - Movimiento libre simétrico (elevación y caída):........................................................ 70
6.9.- Gravímetros........................................................................................................................... 70
6.9.1.- Gravímetro lineal........................................................................................................... 71
6.9.2.- Gravímetros no lineales o astáticos. Astatización......................................................... 72
6.9.2.1.- Correcciones a las medidas gravimétricas............................................................. 77
6.9.2.2.- Correcciones:......................................................................................................... 78
7.- SISTEMAS DE ALTITUDES...................................................................................................... 80
7.1.- Introducción. Números o cota geopotencial..........................................................................80
7.2.- Altitudes dinámicas............................................................................................................... 81
7.3.- Corrección dinámica:.............................................................................................................82
7.4.- Reducción de la gravedad de Poincaré y Prey.......................................................................82
7.5.- Altitudes ortométricas. Corrección ortométrica.................................................................... 83
7.6.- Corrección ortométrica:.........................................................................................................84
7.7.- Altitudes normales. Corrección normal.................................................................................85
7.8.- Comparación de diferentes sistemas de altitudes.................................................................. 87
Página 4
1.- INTRODUCCIÓN A LA GEODESIA FÍSICA
1.1.- Definición y objeto de la Geodesia Física

El objeto principal de la geodesia física es el estudio de la figura de la tierra, entendiendo por ello el
estudio de su forma y dimensiones. Este estudio puede considerarse bajo dos aspectos:
• Geométrico: geodesia geométrica. A partir de medidas de elementos geométricos (ángulos y
distancias) sobre la superficie terrestre se busca determinar la superficie matemática regular que
represente con suficiente precisión la figura de la tierra. Desde el punto de vista práctico la
geometría de la tierra considera el elipsoide de revolución que se define por dos parámetros
geométricos:
• a,b: semiejes mayor y menor
• a: semieje mayor
a−b 1
• α : aplanamiento geométrico α= ≈
a 300
Sobre este elipsoide se realizan los cálculos geodésicos que llevan al conocimiento de las
coordenadas geodésicas de puntos situados sobre la superficie terrestre que son:
ϕ g : latitud geodésica
λ g : longitud geodésica
h: altitud elipsoidal
Cuando intervienen las observaciones astronómicas tenemos las coordenadas astronómicas:

ϕ a : latitud astronómica
λ a : longitud astronómica n- normal al geoide
H: altitud ortométrica (geoide) n’- normal al elipsoide
ε - desviación de la vertical
La diferencia entre las dos altitudes es la ondulación del geoide: h=H+N

• Dinámico: geodesia física. El estudio de la figura de la tierra esta basado en el conocimiento del
campo gravitatorio terrestre lo que implica conocer g, la intensidad del .campo gravitatorio terrestre,
tanto en módulo como en dirección. Podríamos deducir el potencial gravitatorio terrestre, ya que:
→
g = gradW . El potencial, como magnitud escalar, es una función de un punto, o también podemos
utilizar las coordenadas esféricas:
π
W = W(x,y,z) =W(R, θ , λ ) θ - colatitud = − ϕ
2
Conociendo el potencial gravitatorio podemos definir las superficies equipotenciales o superficies
Página 5
de nivel en las que W(x,y,z) = cte. Del campo gravitatorio terrestre tendremos distintas superficies
de nivel entre las que distinguimos W = Wo correspondiente al geoide.
La línea de la plomada es aquella que corta perpendicularmente a todas las superficies
equipotenciales. Es una línea curva y la longitud medida sobre ella desde el punto de la superficie al
geoide es la denominada altitud ortométrica (H). Dos puntos sobre la misma superficie de nivel
tienen el mismo potencial para distinta H.
El objeto de la geodesia física por tanto, es la determinación del geoide (conocer su forma) y el
conocimiento de la función potencial del campo gravitatorio terrestre W. La superficie libre de los
océanos idealizada, que coincide con el nivel medio de los mares en reposo, coincide con el geoide.
1.2.- Diferencia entre geodesia geométrica y geodesia física:

En la geodesia geométrica tenemos en cuenta las direcciones del vector gravedad, ya que la
medición del arco de grado se basa en la fórmula s1 = ϕ 2 − ϕ1 donde las direcciones de la vertical
s
en los extremos del arco de meridiano s determinan el ángulo ϕ 2 − ϕ1 , que también esta presente
en las fórmulas necesarias para hallarla en el módulo elipsódico:
s α = R (ϕ 2 − ϕ1 ) s α ϕ 2 − ϕ1
=
s1º = R.1 º
s1º 1º
En la geodesia física tenemos en cuenta los valores del módulo del vector gravedad. En este caso la
figura exterior de la tierra debe determinarse en función de la magnitud de la fuerza de gravedad
como resultante de las fuerzas de atracción newtoniana y centrífuga.
La fuerza centrífuga varía según la distancia del punto al eje de rotación, será máxima en puntos del
Ecuador, y nula en los polos. También la fuerza de atracción newtoniana se ve afectada por la forma
de dicha superficie. Es decir, la forma de la tierra influye sobre los valores de la gravedad, por tanto,
del estudio de dichos valores podrá deducirse la figura de la tierra.
La geofísica nace cuando Newton enuncia la ley de gravitación universal en 1687. A partir de aquí
se inicia el estudio experimental para determinar la fuerza de atracción terrestre sobre la unidad de
masa que hemos llamado gravedad terrestre y es la resultante de dos fuerzas coincidentes en
módulo pero no en dimensiones (realmente son intensidades de campo o aceleraciones):
→ → → →
g = FN + FC FN - fuerza newtoniana o gravitacional
→
FC - fuerza centrífuga
Página 6
cm
Galileo midió dicha aceleración obteniendo un valor en gales: 1 gal= 1 (C.G.S.)
sg 2
Los valores de la gravedad sobre la superficie terrestre (es decir, sin tener en cuenta las variaciones
por altitud, o lo que es lo mismo: considerando el valor de la gravedad reducidos al geoide) oscilan
entre 978 y 983 gales. Esta variación de 5 gales es debida a la variación de la latitud. Como valor
m
medio se utiliza 981 gales = 9,81 2 .
sg
La gravedad varía dependiendo de la latitud y de la altura del punto.
1.3.- Consideraciones históricas sobre la figura de la Tierra y sus métodos de

estudio
Vamos a distinguir las distintas determinaciones de las dimensiones de la tierra:
1- Modelo esférico:
- Pitágoras hizo referencia sobre la redondez de la tierra.

- Eratóstenes fue el primer geógrafo y geodesta que se basó en observaciones astronómicas para
determinar las dimensiones de la tierra. En el solsticio de verano en Siena, observó el momento en
que la luz solar iluminaba el fondo de un pozo, es decir, los rayos del sol incidían
perpendicularmente sobre la superficie terrestre.
Se trasladó a Alejandría (situada aproximadamente sobre el mismo meridiano que Siena) y a la
misma hora del siguiente solsticio de verano observó y midió la proyección de la sombra sobre el
fondo de un pozo, de manera que conociendo la profundidad del pozo, y la longitud de la sombra
pudo determinar de modo aproximado el ángulo con que los rayos del sol incidían sobre la
superficie tg α = q p , ángulo equivalente a la diferencia angular, con vértice en el centro de la
tierra, entre Siena y Alejandría. Y conociendo aproximadamente la distancia entre ambos lugares
D
determinó el radio de la tierra: R=
α
- Durante la Edad Media se pierde todo interés por el conocimiento de las dimensiones de la tierra.
Este interés resurge en el siglo XVI.
- En el siglo XVII se impulsan los métodos de medida de las dimensiones de la tierra basados en la
medición de la longitud de un arco de un grado sobre un meridiano. Intentan demostrar la
esfericidad de la tierra basándose en que si las medidas del arco de un grado realizadas a diferentes
latitudes son iguales, la tierra tendrá forma esférica.
Holandeses y franceses realizan la primera medición del arco de grado interesante. La aplicación de
métodos de triangulación, tablas de logaritmos para cálculos trigonométricos y el desarrollo del
anteojo para realizar observaciones, mejoran e incrementan la fiabilidad de las mediciones.
2- Modelo elipsódico:
La idea surge después de que Newton enunciara en 1687 la ley de gravitación universal. Huyggens
y Newton tratan de encontrar en la física la explicación y los métodos para determinar la forma de
la tierra.
Newton considera que en algún momento la tierra fue fluida, y debido al movimiento de rotación y
la consiguiente fuerza centrífuga supone un alargamiento del semieje mayor correspondiente a la
posición ecuatorial. Propone el modelo de elipsoide de revolución con achatamiento polar, y llega a
esta conclusión teniendo en cuenta sus conocimientos sobre la fuerza centrifuga, y la de atracción
entre las masas. Presupone que la tierra es de densidad homogénea. Intentó determinar las
a−b 1
dimensiones de los semiejes del elipsoide, y obtuvo un aplanamiento: α= ≈
a 231
Página 7
Huyggens consideró toda la masa terrestre concentrada en su centro y calculó el aplanamiento

obteniendo el valor de 1 258 .
Seguidamente se intentó comprobar la validez del modelo newtoniano por el método geométrico
consistente en la medición del arco de un grado. Si el modelo es correcto, deben obtenerse arcos de
mayor longitud a medida que aumenta la latitud.
La-familia Cassini en Francia* (1700-1718), prolongan el arco de Picard (medido sobre el
meridiano de París) tanto hacia el sur hasta Polieu como hacia el norte de París hasta Dunquerque,
realizando, por triangulación, medida de su longitud. Los resultados obtenidos contradicen el
modelo newtoniano, ya que la longitud del arco de 1° en el sur es mayor que las obtenidas para
posiciones de latitud superior, y en base a ello proponen un modelo con achatamiento ecuatorial y.
ensanchado por los polos. Para aclarar esta controversia organizan en el siglo s.XVIII dos
expediciones con el propósito de medir el arco de grado en latitudes extremas:
- 1735-1742 Perú (zona del actual Ecuador) para medir el arco que atravesaba el ecuador. Participan
Godieu,Bouquer, La Condamíne, y marines españoles,
- 1736-1737 Laponia. Participan Clairaut, Maupestus,...
Como consecuencia de estas mediciones se comprobó y se aceptó como válido el modelo de
elipsoide de rotación con achatamiento polar de Newton.
Surge un interés por medir las dimensiones del elipsoide, tanto por métodos geométricos como
físicos (métodos gravimétricos).
3- El geoide y métodos actuales:
- Clairaut en 1743 sugiere una distribución de la masa de la tierra por capas elipsóidicas de
densidades variables. Sus estudios teóricos ponen de manifiesto el teorema para la obtención del
aplanamiento por métodos gravimétricos:
3 a−b γ polo − γ ecuador
α+β = m α= β=
2 a γ ecuador
w 2a
m=
γ ecuador
v2 m- masa
Fc = m
ρ v- velocidad
Fc = m.w 2 .ρ ρ − radio de curvatura
w- velocidad de rotación terrestre
ρ= a
Para calcular la gravedad en un punto cualquiera:
γ punto = γ ecuador (1 + β.sen 2 ϕ punto )
En esta ecuación aparecen dos incógnitas, β y γ ecuador , (la latitud es conocida por observación
astronómica), por lo que sí medimos la gravedad en dos puntos tendremos un sistema de dos
ecuaciones para darles solución. Pero con este método, para deducir α necesitamos conocer
también a (semieje mayor), que sigue siendo determinado por métodos geométricos.
- En el siglo XIX Gauss propone incorporar el método de mínimos cuadrados para aumentar las
precisiones en los cálculos de triangulación.
Se continúan las mediciones para determinar los parámetros del elipsoide, pero los errores
obtenidos son superiores a los que cabría esperar, de lo cual se dedujo que no se trataba
exactamente de un elipsoide. Es entonces cuando la desviación de la vertical deja de ser tratada
como un error de observación. Gauss propone como figura de la tierra el geoide: superficie de nivel
(o sup. equipotencial del campo gravitatorio) con la que coincide, en parle, la correspondiente al
Página 8
nivel medio de los mares en reposo.

- Aparecen otras figuras que contribuyen al desarrollo de la gravimetria. Stockes a finales del siglo-
XIX llega a una fórmula integral que permite calcular la ondulación del geoide por métodos
gravimétricos.
- Vening Meinesz en 1927 como consecuencia de estudios teóricos basados en la fórmula de
Stockes propone calcular la desviación de la vertical por métodos gravimétricos. Miden anomalías
de la gravedad en base a las cuales se pueden conocer tanto la ondulación como la desviación de la
vertical.
- Molodensky en 1945 define un sistema normal de altitudes e introduce una superficie
equipotencial de referencia el quasi-geoide.
- Actualmente conocidas las posiciones de los satélites, y utilizando la reflexión sobre la superficie
del océano de una emisión de ecorrádar podemos conocer la posición del nivel medio de los mares
con respecto al elipsoide.
Página 9
2.- TEORÍA DEL POTENCIAL. EL CAMPO GRAVITATORIO

TERRESTRE
2.1.- Ley de Newton de la Gravitación Universal

La enuncia basándose en observaciones astronómicas a planetas en su giro alrededor del sol,
aceptado ya el modelo heliocéntrico (no geocéntrico).
La forma de atracción entre dos partículas o puntos materiales m y m’, situados a una distancia L, es
proporcional al producto de las masas e inversamente al cuadrado de la distancia que las separa:
m.m'
F=G 2
L
Es una fuerza atractiva, su dirección es la de la línea imaginaria que une las dos masas y su sentido
depende de la posición de las partículas.
G es la constante de la gravitación universal y representa la fuerza que se ejercen dos masas
unitarias situadas a la unidad de distancia. Tiene dimensión de magnitud física:
[G] = F.L2 .M −2 = M −1 .L3 .T −2 F = M.L.T −2
La primera medida de G la realizó Cavendish en 1798 y su valor es:
N⋅m2
G=6.6⋅10−11 , en unidades del sistema internacional.
Kg 2
Como 1⋅N=1⋅Kg⋅m⋅s−2 entonces,
m3
G=6.6⋅10−11
Kg⋅s 2
También podremos calcularlo para el sistema cegesimal:
dy⋅cm2
G=6.6⋅10−8 , expresado en ucgs.
gr2
Como sabemos que una dina equivale a 1⋅dy=1⋅gr⋅cm⋅s−2 , por tanto:
3
−8 cm
G=6.6⋅10
gr⋅s2
Para conocer numéricamente la fuerza es necesario conocer el valor de esta constante.
Toda partícula que se desplaza alrededor de otra con una relación inversamente proporcional al
cuadrado de la distancia, lo que va a describir es una elipse o una hipérbola según la energía.
2.2.- Campo de fuerza newtoniano. Intensidad de campo y potencial.
2.2.1.- Campo newtoniano debido a una masa puntual:

Supongamos una masa m situada en un cierto punto O'( ε, η, ζ ). Situamos una masa unidad en un
punto P(x,y,z):
Página 10
→
Definimos la magnitud intensidad de campo f n en P como la fuerza atractiva que la masa m
(atrayente) ejerce sobre dicha masa unidad a la distancia l.
Según la ley de gravitación universal:
→ m
f N = f N = G. 2
l
→ m → ∂V → ∂V → ∂V → →
donde f N = −G. 2 . u = gradV = i+ j+ k = f ( x, y, z) es función vectorial del punto y
l ∂x ∂y ∂z
N Kg.m.s −2
tiene dimensiones de aceleración: = = m.s − 2
Kg Kg
Esta magnitud tiene carácter vectorial. Para describirla nos valemos de un vector unitario en la
→
dirección O'P y cuyo sentido nos determina el signo adoptado para f N , en este caso negativo.
→
→ l → → → → →
u= donde l = O' P = ( x − ε) i + ( y − η) j + (z − ζ ) k
l
l = ( x − ε) 2 + ( y − η) 2 + (z − ζ ) 2
→ x−ε→ y−η→ z−ζ → x −ε
Las componentes u = i+ j+ k son los cosenos directores cos α = ,
l l l l
y−η z−ζ
cos β = y cos γ = .
l l
Las expresiones cartesianas del vector intensidad de campo serán:
m m →
(f N ) x = −G. 2 cos α = −G. 3 .( x − ε) = − f N cos α
l l
m m →
(f N ) y = −G. 2 . cos β = −G. 3 .( y − η) = − f N cos β
l l
m m →
(f N ) z = −G. 2 cos γ = −G. 3 .(z − ζ ) = − f N cos γ
l l
Demostraremos que este campo de fuerzas es conservativo o irrotacional, es decir, deriva de una
m
función potencial. Por tanto, existe una función potencial (función de punto) V = G. = V( x, y, z)
l
→
de naturaleza tal que f N = gradV .
Para ello comprobaremos la identidad de cada una de las componentes del vector grad V con las del
vector intensidad de campo.
Página 11
Supondremos, a la hora de derivar, que la masa puntual m permanece fija en sus coordenadas
ε, η, ζ , siendo variable la posición del punto P(x,y,z):
∂l 2( x − ε)
− G.m.
∂V ∂ 1 ∂x = −G.m. 2l x−ε
(gradV) x = = G.m. .( ) = = −G.m. 3 = (f N ) x
∂x ∂x l l 2
l 2
l
∂l 2( y − η)
− G.m.
∂V ∂ 1 ∂y y−η
(gradV) y = = G.m. .( ) = = −G.m. 22l = −G.m. 3 = (f N ) y
∂y ∂y l l 2
l l
∂l 2(z − ζ )
− G.m.
∂V ∂ 1 ∂z = −G.m. 2l z−ζ
(gradV) z = = G.m. .( ) = = −G.m. 3 = (f N ) z
∂z ∂z l l 2
l 2
l
→
Se verifica así que f N = gradV y se dice entonces que es un campo de fuerzas conservativo o
irrotacional.
Si consideramos la masa atrayente situada en el origen de coordenadas: ε = η = ζ = 0 .
2.2.2.- Campo newtoniano debido a un sistema de masas puntuales:

En este caso se aplica el principio de superposición que afecta a los vectores intensidad de campo y
a los potenciales, de tal modo que:
→ → → → n →
f N = f N1 + f N 2 + ... + f Nn = ∑i =1 f N y estará dirigido hacia el centro de masas del sistema.
→
Como f N
deriva del potencial:
→ m1
f N1 = gradV1 siendo V1 = G. l
1
→ m2
f N 2 = gradV2 siendo V2 = G. l
2
→ mn
f Nn = gradVn siendo Vn = G. l
n
En base a las propiedades de la función gradiente:
→
f N = gradV1 + gradV2 + ... + gradVn = grad(V1 + V2 + ... + Vn ) = gradV
n mi
donde V = V1 + V2 + ... + Vn = G.∑i =1
li
Página 12
2.2.3.- Campo newtoniano debido a un cuerpo extenso:

Suponemos una masa distribuida de forma continua (no necesariamente homogénea, sino "sin
huecos") en un cuerpo de volumen τ .
Consideramos elementos de masa infinitesimal (dm) en torno a punto P’( ε, η, ζ )que distará l del
punto P(x.y,z).
Generalizando la fórmula del potencial correspondiente a un elemento diferencial de masa:

dm dm ρ .dτ
V = G. para todo volumen: V = G ∫∫∫τ = G ∫∫∫τ v
l l l
ρ v = ρ v (ε, η, ζ ) - función de punto denominada densidad cúbica
dτ = dε.dη.dζ - diferencial de volumen
Es la fórmula general del potencial en un cuerpo extenso.
→ →
Obtenemos ahora el vector intensidad de campo como gradiente de la función potencial: F N = f n . u
∂V ∂ dm ∂ 1 x−ε x−ε
(f N ) x = = G. ∫∫∫τ = G ∫∫∫τ  dm = − G ∫∫∫τ dm = −G ∫∫∫τ ρ v dτ
∂x ∂x l ∂x  l  l 3 l3
∂V ∂ dm ∂ 1 y−η y−η
(f N ) y = = G. ∫∫∫τ = G ∫∫∫τ  dm = − G ∫∫∫τ dm = −G ∫∫∫τ ρ v dτ
∂y ∂y l ∂y  l  l3 l3
∂V ∂ dm ∂ 1 z−ζ z−ζ
(f N ) z = = G. ∫∫∫τ = G ∫∫∫τ  dm = − G ∫∫∫τ dm = −G ∫∫∫τ ρ v dτ
∂z ∂z l ∂z  l  l 3 l3
Sumando las componentes obtendremos la función intensidad como:
→ → →
→ ( x − ε) i + ( y − η) j + ( z − ζ ) k dm →
f N = −G ∫∫∫τ dm = −G ∫∫∫τ u1
l3 l2
→
En este caso u 1 no es constante. Varía en dirección y sentido dependiendo de la posición del
diferencial de masa y del punto P.
Las fórmulas del potencial siguen inalteradas dado que a la hora de derivar consideramos que el
diferencial de masa no varía de posición, ε = η = ζ = ctes . La variación de la posición del diferencial
de masa y por tanto de sus coordenadas, son consideradas a la hora de integrar sobre el volumen
para hallar la función intensidad.
Página 13
2.3.- Definición y propiedades de la superficie equipotencial: intensidad de campo

y potencial:
Definimos superficie equipotencial como el lugar geométrico de los puntos en que el potencial es
constante. Para cada valor de V tendremos una superficie equipotencial.
Ej.: V = V( x , y, z) = V(R , θ, λ ) = cte
El vector grad V tiene la dirección de la normal a la superficie equipotencial en el punto
considerado, y dirigida en el sentido del crecimiento de la función.
En el caso de la tierra la función V crece con la profundidad.
Vamos a demostrar que en una superficie equipotencial no varía el potencial y que entre dos
superficies equipotenciales la máxima variación está en la perpendicular a ambas.
Para demostrar la perpendicularidad nos desplazamos desde el punto P(x,y,z) hasta otro
infinitamente próximo (x+dx, y+dy, z+dz):
∂V ∂V ∂V → → →
dV = dx + dy + dz ⇒ dV = gradV. ∂l = FN . ∂l
∂x ∂y ∂z
→ → → → ∂V → ∂V → ∂V →
donde ∂l = ∂x. i + ∂y. j + ∂z. k y gradV = i+ j+ k = ∇h
∂x ∂y ∂z
1. Al pasar de P a Q situado en la misma superficie equipotencial, lo haremos a lo largo del vector
→
∂l t tangencial y la variación del potencial será nula, por lo que: 0=gradV. ∂l .Al tratarse de un
t
producto escalar entre dos vectores de módulos no nulos, el coseno del ángulo que forman ha de
ser nulo; por tanto son perpendiculares.
2. Si nos desplazamos al punto Q' en la dirección del vector normal ∂l n :
→ ∂V
dV ≠ 0 ⇒ dV = gradV. ∂l n = gradV .dn. cos 0 = gradV .dn ⇒ gradV = es la derivada
∂n
→
direccional de la función V en la dirección de la normal siendo ∂l n = dn .
3. La circulación del vector intensidad de campo a lo largo de un contorno cerrado es cero; es decir,
→
el trabajo necesario para desplazar una masa bajo la acción de la fuerza f N a lo largo de un
contorno cerrado es nulo. En el caso de la "masa unidad, su expresión será:Trabajo =
→ → →
∫ f N . ∂l = ∫ gradV. ∂l = ∫ dV = 0
4. Si el contorno no es cerrado, sino que la masa se desplaza desde el punto A al B, tendremos:
→ → →
Trabajo = ∫AB f N . ∂l = ∫AB gradV. ∂l = ∫AB dV = VB − VA
, por lo que deducimos que el trabajo no depende del camino recorrido, sino exclusivamente de
la posición de los puntos extremos.
5. Estudiamos ahora el valor de la función potencial cuando actúa sobre una masa situada en un
punto del infinito. Para ello, se estudia inicialmente su influencia sobre un punto situado a gran
distancia, y se comprueba que V tiende al mismo valor que el que produciría una masa puntual
situada en el centro de masas del cuerpo y en la que se encontrase concentrada toda la masa de
dicho cuerpo.
M M
V = G. ⇒ lim R − >∞ V = lim R − >∞ G. =0 M = ∑ dm
R R
6. Se considera ahora la siguiente situación:
Página 14
M M
Sabemos que el potencial respecto a la masa tomará valores entre: G . V G .
r1 r2
Si multiplicamos ahora por un ri cualquiera situado en la alineación de la figura tendremos:
M⋅r i M⋅r i
G. V⋅r iG .
r1 r2
, tomando límites:
M⋅r i M⋅r i
lim r i ∞ G .  lim r i ∞ V⋅r ilim r i ∞ G .
r1 r2
G⋅M⋅1lim r i  ∞ V⋅r i G⋅M⋅1
V.R = G.M ⇒ lim R − >∞ V.R = G.M = cte
Se dice que la función potencial está normalizado al infinito.
2.4.- Flujo vector intensidad de campo. Teorema de Gauss

Teorema de Gauss: el flujo del vector intensidad de campo a través de una superficie cerrada
evaluado algebraicamente respecto a la normal exterior es negativo e igual al producto de la
constante 4πG por la masa encerrada en su interior.
2.4.1.- Caso de una masa puntual interior m
Considero una superficie cerrada de cualquier forma alrededor del punto O donde se encuentra m.
Tomo un elemento de superficie infinitesimal dS en torno al punto P y lo represento mediante un
→
vector ∂S de módulo igual a dS, dirección perpendicular al elemento de superficie y sentido hacia
el exterior de la superficie S.
Veremos que sólo influyen las masas interiores y no las exteriores a la superficie:
Página 15
→ →
Flujo = Φ = ∫∫S f N . dS = −4πGm int erior
→ →
Flujo = Φ = ∫∫S f N . dS = ∫∫S f N .dS. cos α donde α = π − ϕ ⇒ cos α = − cos ϕ
m
Flujo = ∫∫S − f N .dS. cos ϕ = ∫∫S − G. .dS. cos ϕ
r2
→
donde dS. cos ϕ = dS n , es decir, la proyección del vector dS sobre la dirección radial OP , siendo
dS n un elemento de casquete esférico con radio r y centro en O:
dS n dS n
Flujo = ∫∫S − G.m. 2
= −G.m ∫∫S
r r2
dS n dS. cos ϕ
dΩ = = es el ángulo sólido elemental o ángulo sólido subtendido desde O por cualquier
r2 r2
elemento de superficie dS o cualquier elemento de superficie esférica dS n . Se mide en
estereorradianes y el ángulo sólido completo es de Ω = 4π estereorradianes.
Por tanto:
dS n 4πr 2
Flujo = −G.m ∫∫S dΩ siendo Ω = ∫∫S dΩ = ∫∫S = = 4π
r2 r2
Flujo = −4πGm int erior
Un flujo elemental será:
→ → m
Flujo = f N . dS = f N dS. cos α = −f N dS. cos ϕ = −G dS. cos ϕ = −G.m.dΩ
r2
2.4.2.- Caso de una masa exterior a la superficie:
I1⋅∂ S10
I2⋅∂ S20
Por tanto:
d =d  1d  2=0
La contribución neta de una masa exterior a la superficie cerrada es cero, ya que los flujos
Página 16
elementales dos a dos se anulan.
2.4.3.- Caso varias masas interiores
Caso de tener varias masas interiores m int = ∑in=1 m i y el flujo será:

→ →
Flujo = ∫∫S f N . dS = −4πG (m1 + m 2 + ... + m n )
Si el cuerpo tuviese parte de su masa dentro del cuerpo del que se calcula el flujo y parte fuera del
mismo, tendríamos que considerar únicamente la parte que se encuentra en el interior.
2.5.- Ecuaciones de Laplace y de Poisson. Funciones armónicas

→
Definida f N como grad V , sabemos que el campo de gravitación newtoniano es conservativo o
irrotacional. Veremos que el potencial V cumple una serie de propiedades a tener en cuenta.
Ecuación de Poisson: el operador laplaciano del potencial V es igual a − 4πGρ siendo ρ la
dm
densidad cúbica en el punto considerado: dm = ρ.dτ ⇒ ρ=
dτ
∆V = −4πGρ
En un punto fuera de la masa atrayente, donde ρ = 0 , se verifica la ecuación de Laplace: ∆V = 0 .
Para demostrarlo nos valdremos de dos teoremas:
→ →
- Teorema de Gauss: ∫∫S f N . dS = −4πGM int
→ →
dm = ρ.dτ ⇒ M int = ∫∫∫τ ρ.dτ ⇒ ∫∫S f N . dS = −4πG ∫∫∫τ ρ.dτ
→ → →
- Teorema de la divergencia o de Gauss Ortogradsky: ∫∫S f N . dS = ∫∫∫S div f N .dτ
τ - volumen encerrado por la superficie S.
Igualando ambas expresiones:
→ →
− 4πG ∫∫∫τ ρ.dτ = ⇒ ∫∫∫τ − 4πGρ.dτ = ∫∫∫ div f N .dτ
∫∫∫S div f N .dτ S
→ → →
f N = gradV = ∇V ⇒ div f N = div(gradV) = div(∇V) = ∇. f N = ∇.∇V = ∇ 2 V = ∆V
Al ser dos integrales extendidas al mismo volumen, podemos agruparlas en una integral única,
→
sacando factor común dτ , y sustituyendo el valor de div f N :
∭τ −4πGρ−ΔV . dτ=0
Ecuación que se cumplirá cuando 4πGρ + ∆V = 0 ⇒ ∆V = −4πGρ , siendo la Ecuación de
Poisson que se cumple siempre que exista masa y puede comprobarse si conocemos la función V.
En caso de que no exista masa se cumplirá la ecuación de Laplace.
Una función se denomina armónica en los puntos de una región si cumple la ecuación de Laplace.
Toda función armónica es analítica, es decir, es continua y tiene derivadas continuas de cualquier
orden en aquellos puntos en que está definida y cumple la ecuación de Laplace. Por tanto las
soluciones V de la ecuación de Laplace son armónicas.
2.6.- Campo y potencial de las fuerzas centrífugas
Página 17
Consideramos un sistema de ejes cartesianos. En una determinada posición P hay una masa unidad
arrastrada por el movimiento de rotación del sistema. Sobre esta partícula actúa una fuerza
centrífuga en la dirección del radio de la circunferencia que describe dicha partícula al girar el
sistema y de sentido hacia fuera.
→ → v2 →
Fc = m.a c . u = m. .u = m.ω 2 .rp . u
r
→ → → →
→ rp → rp → →
m' = 1 ; u = = x i + y j ⇒ 2
Fc = 1.ω .rp . = ω2 (x i + y j )
rp rp rp
En el sistema de ejes cartesianos tiene las siguientes componentes: (Fc ) x = ω 2 x , (Fc ) y = ω 2 y y

→ →
(Fc ) z = 0 . Esta componente es nula porque la F es perpendicular al eje z. Tratamos F como un
c c
campo vectorial función de punto. Dimensionalmente coincide con una aceleración. Podemos
comprobar que es conservativo o que deriva de un potencial:
→ → →
i j k
→ → ∂ ∂ ∂ →
→
Fc = gradΦ c rot Fc = ∇ × Fc = =0
∂x ∂y ∂z
ω x ω2 y
2
0
∂Φ ω 2 y = ∂Φ
Debemos encontrar la función que cumple las tres ecuaciones diferenciales: ω 2 x = , y
∂x ∂y
∂Φ
0= . De la última vemos que Φ no va a depender de la variable z. De la 2ª obtenemos que
∂z
1 1
Φ = ∫ ω 2 y.dy = ω 2 y 2 + f 2 ( x ) y de la 1ª que Φ = ∫ ω 2 x.dx = ω 2 x 2 + f 1 ( y) donde la constante de
2 2
integración es una función que dependerá del resto de las variables. En este caso únicamente de y,
dado que la tercera componente de Fc nos indica que no depende de z.
Buscamos los valores de las funciones que figuran como constantes de integración. Suelen
corresponderse con los términos que aparecen en el resto de las soluciones correspondientes a sus
1 1
variables, así f1 ( y) = ω 2 y 2 y f 2 ( x ) = ω 2 x 2 , de lo que deducimos que la expresión correcta para
2 2
el potencial será:
1 2 2 1 2 2 1 2 2
Φ= ω x + ω y = ω (x + y 2 )
2 2 2
Página 18
→ → →
rp = x i + y j ⇒ rp2 = x 2 + y 2
1 2 2 → → → →
Φ= ω rp ⇒ Fc = ω 2 ( x i + y j ) = ω 2 . rp = gradΦ
2
Calculamos el operador laplaciano que como no es nulo no cumple la ecuación de Laplace, luego
∂ 2Φ c ∂ 2Φ c
no es armónica: ∆Φ c = 2
+ 2
= ω 2 + ω 2 = 2ω 2 .
∂x ∂y
El potencial centrífugo en coordenadas esféricas es:
1 2 2 1 2 2
rp = R.senθ ⇒ Φc = ω rp = ω R sen 2 θ
2 2
2.7.- Gravedad y potencial gravitatorio terrestres

Se llama gravedad terrestre o intensidad de campo gravitatorio terrestre a la fuerza resultante sobre
la unidad de masa de la fuerza newtoniana debida al cuerpo de la tierra y de la fuerza centrífuga
→ → →
debida a la rotación terrestre: g = f N + Fc .
Una masa unidad en reposo sobre la superficie de la tierra está sometida a estas dos fuerzas. No
tenemos en cuenta los efectos de atracción lunisolares o de otros astros ni las fuerzas de Coriolis.
Las variaciones que se producen debidas a estas otras causas son temporales, y habrá que realizar
las correcciones necesarias a las observaciones experimentales para hallar la g correcta.
Métodos para calcular la gravedad g:
4π 2 L
1- Experimento de Cavendish: g =
T2
2- Medición de longitudes recorridas por un cuerpo en caída libre durante un intervalo de tiempo:
1 2
S = v 0 .t + g.t
2
Como hemos visto ambas fuerzas derivan de potenciales y haciendo uso de las propiedades del
gradiente tendremos que:
→
g = gradV + gradΦ c = grad( V + Φ c ) = grad( W )
W = V + Φ c - es el potencial gravitatorio terrestre.
∂W x−ε
Las componentes del vector intensidad son: g x = = −G ∫∫∫τ 3 dm + ω 2 x ,
∂x l
∂W y−η ∂W z − ζ
gy = = −G ∫∫∫τ 3 dm + ω 2 y y g z = = −G ∫∫∫τ 3 dm .
∂y l ∂z l
Página 19
Propiedades de los potenciales:

- Como ∆V = −4πGρ y ∆Φ c = 2ω 2 ⇒ ∆W = −4πGρ + 2ω 2
- Si estamos fuera de las masas atrayentes: ρ = 0 , ∆V = 0  ∆W = 2ω 2 , no armónica.
dm 1 2 2
- Expresión cartesiana del potencial gravitatorio: W = G ∫∫∫τ + ω (x + y 2 )
l 2
dm 1
siendo V = G ∫∫∫τ y Φ c = ω2 (x 2 + y 2 )
l 2
dm 1 2 2
En esféricas: W = G ∫∫∫τ + ω R sen 2 θ
l 2
2.8.- Superficies de nivel y líneas de la plomada:

Partimos de la expresión W = V + Φ c que igualada a una constante representa una superficie de
nivel o superficie equipotencial del campo gravitatorio: W ( x, , y, z) = cte o W (R , θ, λ) = cte .De todas
las superficies equipotenciales, distinguimos como geoide aquella en la que el valor del potencial
coincide con el que toma dicho potencial en puntos correspondientes al nivel medio de los mares en
reposo, es decir, coincide en parte con el nivel medio del mar: W = W0 .
La línea de la plomada es una línea geométricamente curva, perpendicular a todas las superficies
→
equipotenciales desde el punto considerado hasta el geoide: g = gradW .
→
El vector g es tangente a la línea de la plomada en cada punto.
→ ∂W ∂W ∂W
dW = gradW. dl n = dx + dy + dz
∂x ∂y ∂z
→ → →
como g = gradW ⇒ dW = g . dl n
Podemos demostrar que gradW es perpendicular a la superficie de nivel W = cte dado que si nos
desplazamos a lo largo del vector tangente a ella obtenemos que dW = gradW.dl i . cos α = 0 ⇒
La altitud ortométrica H es la distancia medida a lo largo de la línea de la plomada desde el geoide o

superficie de referencia y el punto considerado: PP0
dW dW
gradW = dn dn = dH ⇒ g=−
dn dH
2.9.- Desarrollo del potencial gravitatorio terrestre. Polinomios de Legendre

Supongamos un punto fuera de las masas atrayentes, se cumple entonces la ecuación de Laplace.
Página 20
Cuando esto sucede, dicha función puede desarrollarse en términos de las funciones
armonicoesféricas. Al desarrollar la función 1 l nos aparecerán los polinomios de Legendre,
también denominados armónicos esféricos de superficie zonales.
dm ρ.dτ
1- Desarrollamos por el teorema del coseno: V = G ∫∫∫τ = G ∫∫∫τ
l l
1
  r'  2 r'   2  2
−
⇒ 1 = 1 1 +  r '  − 2 r ' cos γ 
l 2 = R 2 + r ' 2 −2Rr '. cos γ = R 2 1 +   − 2 cos γ 
  R  R 
l R   R  R 
r'
2- Mediante desarrollo de Taylor y haciendo los cambios cos γ = t y = h obtenemos como
R
[ ]
1
resultado que: 1 = 1 1 + h 2 − 2ht 2
−
l R
1
− 1 3 15 3 y
3-Aplicamos a nuestro desarrollo mediante Taylor que: (1 + x ) 2 =1− x + x2 − x ...
2 8 48
obtendremos lo siguiente:
1 1 1 2 3 15 2 
= 1 − (h − 2ht ) + (h 2 − 2ht ) 2 − (h − 2ht ) 3 + ...
l R 2 8 48 
como: (h 2 − 2ht ) 2 = h 4 − 4h 3 t + 4h 2 t 2
(h 2 − 2ht ) 3 = h 6 − 3h 4 2ht + 3h 2 4h 2 t 2 − 8h 3 t 3
1 1 1 2 3 15 6 
= 1 − (h 2ht ) + (h 4 − 4h 3 t + 4h 2 t 2 ) − (h − 3h 4 2ht + 3h 2 4h 2 t 2 − 8h 3 t 3 + ...)
l R 2 8 48 
4- Agrupamos en potencias de h:
1 1 0  1 3   3 15  
= 1h + ht + h 2  − + t 2  + h 3  − t + t 3  + h 4 ...
l R  2 2   2 6  
1 3 3 15
Los polinomios de legendre son: − + t 2 ; − t + t 3 ...
2 2 2 6
5- Sustituimos los valores de h y t:
n
1 1 ∞  r'  r' n
= ∑ n =0   Pn (cos γ ) = ∑∞
n =0 n +1 Pn (cos γ )
l R R R
1 dn
La expresión general de los polinomios de legendre es: Pn ( t ) = n n
( t 2 − 1) n
2 .n! dt
Para los distintos grados de n los valores de los polinomios son los siguientes:
n=0 P0 (cos γ ) = 1
n =1 P1 (cos γ ) = t = cos γ
1 3 1 3
n=2 P2 (cos γ ) = − + t 2 = − + cos 2 γ
2 2 2 2
3 5 3 3 5
n =3 P3 (cos γ ) = − t + t = − cos γ + cos 3 γ
2 2 2 2
Lo introducimos en la expresión del potencial:
dm 1  r' 
n 
V = G ∫τ = G ∫τ  ∑∞ n =0   Pn (cos γ )  dm
l  R R 
Calculamos los dos primeros términos del potencial:
Página 21
1 M
n=0 V0 = G ∫τ dm = G
R R
1
n = 1 V1 = G 2 ∫τ r '. cos γ.dm
R
V = V0 + V1 + ...
→ →
→ → R . r ' Xx '+ Yy'+ Zz' R = cte
Pero R . r ' = R.r '.cos γ ⇒ cos γ = =
R.r ' R . r' R’ = variable
V=G
1
2 ∫τ
r'
Xx '+ Yy'+ Zz'
R.r '
1 1
[
dm = G 3 ∫τ ( Xx '+ Yy'+ Zz') dm = G 3 X ∫τ x ' dm + Y ∫τ y' dm + Z∫τ z' dm ]
R R R
Recurriendo a la dinámica del sólido, sabemos que las coordenadas del centro de masas de un
∫ x '.dm ' ∫τ y'.dm ∫ z'.dm
cuerpo vienen dadas por las expresiones: x 'c = τ , yc = y z 'c = τ .
M M M
Si hacemos coincidir el origen del sistema de referencia con el centro de masas del cuerpo que crea
el potencial, en este caso la tierra: x 'c = y 'c = z 'c = 0 ⇒ ∫τ x '.dm = ∫τ y'.dm = ∫τ z'.dm = 0 . Por lo que
1 3 1
V1 = 0 y V2 = G 3 ∫τ
r ' 2  cos 2 γ − dm .
R 2 2
Como no sabemos la distribución de densidades, lo escribiremos en función de los momentos y de
los productos de inercia del cuerpo de la tierra respecto al sistema de referencia con origen en el
centro de la misma.
1 3 1 2 Xx '+ Yy'+ Zz' R 2 .r ' 2
V2 = G ∫τ  ( r '. cos γ ) 2
− r ' dm , como r '. cos γ = y r ' 2
= 
R3 2 2  R R2
1  3  Xx '+ Yy'+ Zz'  1  R 2 r ' 2 
2
V2 = G 3 ∫τ    −  dm
R  2  R  2  R 2 
V2 = G
2R
1
3 ∫τ
[ ]
3( Xx '+ Yy'+ Zz') 2 − R 2 r ' 2 dm , pero R 2 = X 2 + Y 2 + Z 2 y r ' 2 = x ' 2 + y' 2 + z' 2
V2 = G
2R
1
3 ∫τ
[ ( )(
3( Xx '+ Yy'+ Zz') 2 − X 2 + Y 2 + Z 2 x ' 2 + y' 2 + z' 2 dm )]
V2 = G
1
2R 3
(
2 2 2 2 2 2
∫τ [3 X x ' + Y y' + Z z' +2XYx ' y'+2XZx' z'+2YZy' z' − )
( )
− X 2 x ' 2 + X 2 y ' 2 + X 2 z ' 2 + Y 2 x ' 2 + Y 2 y ' 2 + Y 2 z ' 2 + Z 2 y' 2 + Z 2 x ' 2 + Z 2 y ' 2 + Z 2 z ' 2 ]
Agrupando términos tendremos los siguientes términos integrales:
G G
X2 5 ∫τ
( 2x ' 2 − y' 2 − z' 2 )dm XY ∫τ 6x ' y' dm = 0
2R 2R 5
G G
Y2 5 ∫τ
( 2 y' 2 − x ' 2 − z' 2 )dm XZ ∫τ 6 x ' z' dm = 0
2R 2R 5
G G
Z2 5 ∫τ
( 2z' 2 − y' 2 − x ' 2 )dm YZ ∫τ 6 y' z' dm = 0
2R 2R 5
Aparecen en estos tres últimos términos los productos de inercia respecto a dos planos del
sistema de referencia aplicados al cuerpo de la tierra. Si los ejes están referidos son los ejes
principales de inercia del cuerpo. Dichos productos de inercia son nulos y por tanto las integrales en
las que intervienen son nulos.
Página 22
Px 'y ' = ∫τ x ' y' dm

respecto a los planos (x=0, y=0)
Px 'z ' = ∫τ x ' z' dm respecto a los planos (x=0, z=0)
Py 'z ' = ∫τ y' z' dm respecto a los planos (y=0, z=0)
Si los ejes a los que están referidos son los ejes principales de inercia del cuerpo, dichos productos
de inercia son nulos y por tanto las integrales en que intervienen son nulas. Por eso en Geodesia
utilizamos un elipsoide de revolución para que tengamos ejes principales de inercia.
GX 2
2R 5
2 2 2
∫τ (2x ' − y' −z' )dm =
GX 2
2R 5
[ 2 2 2 2 2 2
∫τ ( x ' + y' ) + ( x ' + z' ) − 2( y' + z' ) dm = ] GX 2
2R 5
(C + B − 2 A )
GY 2
2R 5
2 2 2
∫τ (2 y' − x ' −z' )dm =
GY 2
2R 5
[ 2 2 2 2 2 2
∫τ ( x ' + y' ) + ( y' + z' ) − 2( x ' + z' ) dm = ] GY 2
2R 5
(C + A − 2B)
GZ 2
2R 5 ∫τ
(2z' 2 − y' 2 − x ' 2 )dm =
GZ 2
2R 5 ∫τ
[( x ' + z '
2 2
]
) + ( y' 2 + z' 2 ) − 2( y' 2 + x ' 2 ) dm =
GZ 2
2R 5
( A + B − 2C)
Sabemos que el momento de inercia de un sistema de partículas respecto a un eje es igual al

sumatorio de los productos de sus masas por los cuadrados de las distancias que las separan de
dicho eje I = ∑in=1 m i ri2 , así para un cuerpo extenso tendremos I = ∫c r 2 dm .
Por lo tanto, llamando:
A = ∫τ ( y' 2 + z' 2 )dm ⇒ Momento de inercia respecto del eje x.
B = ∫τ ( x ' 2 + z' 2 )dm ⇒ Momento de inercia respecto del eje y.
C = ∫τ ( y' 2 + x ' 2 )dm ⇒ Momento de inercia respecto del eje z.
Ya que ( y' + z' ) , ( x ' + z' ) y ( y' 2 + x ' 2 ) son las respectivas distancias de dm a los ejes x, y y z.
2 2 2 2
De manera que el término V2 nos quedará:

V2 =
G
2R 5
[X 2
(C + B − 2A ) + Y 2 ( A + C − 2B) + Z 2 ( A + B − 2C) ]
Recopilando estos tres términos:
V = V0 + V1 + V2 = G
M
R
+0+
G
2R 5
[
X 2 (C + B − 2A) + Y 2 (A + C − 2B) + Z 2 (A + B − 2C) ]
Para calcular el potencial gravitatorio terrestre con esta aproximación, sólo habría que añadir a esta
expresión la del potencial centrifugo:
W = V + Φc = G
M
R
+0+
G
2R 5
[ 1
X 2 (C + B − 2A ) + Y 2 ( A + C − 2B) + Z 2 ( A + B − 2C) + ω 2 (X 2 + Y 2 )
2
]
En caso de que los ejes del sistema de referencia no sean los principales del cuerpo, habría que
añadir los términos, en este caso no nulos, correspondientes a los productos de inercia.
Hemos impuesto tres condiciones:
4- Origen en el centro de masas de la tierra.
5- Ejes del sistema son los principales de inercia del cuerpo.
6- Simétrica en torno al Z(A-B) de inercia también.
Ahora imponemos una cuarta condición:
La tierra tiene simetría rotacional en torno al eje Z, es decir, en los momentos de inercia respecto a
los ejes x e y son iguales A = B ≠ C . Por otra parte haremos coincidir el eje Z con el eje medio de
rotación terrestre. De modo que A será el momento de inercia respecto a un eje ecuatorial y C
respecto al eje de rotación polar. Sustituyendo por tanto B por A tenemos:
V=G
M
+
G
R 2R 5
[
X 2 (C − A) + Y 2 (C − A) + Z 2 (2A − 2C) ]
Página 23
V=G
M
+
R 2R
G
5
[
(C − A ) X 2 + Y 2 + 2 Z 2 ]
R 2 = X 2 + Y 2 + Z2 ⇒ X 2 + Y 2 = −Z 2 + R 2 ⇒ V=G
M
+
R 2R
G
5
[
(C − A ) R 2 − 3Z 2 ]
Sabemos que existe un estrechamiento polar y un ensanchamiento ecuatorial. Los momentos de
inercia serán distintos debido a la irregular distribución de las masas en el interior de la tierra.
Lógicamente las masas estarán más alejadas del eje z que de los ejes x e y.
La expresión de V en coordenadas esféricas es:
V=G
M
+
G
R 2R 5
[
(C − A ) R 2 − 3R 2 cos 2 θ ] Z = R. cos θ
Sacando factor común R 2 y realizando un cambio de signo e introduciendo 1/2 en el paréntesis:

M 2⋅G
V =G − 3  C− A 
R 2R
3
2
cos2 θ−
[
1
2 ]
3 1 M C−A 
como P2 (cos θ) = cos 2 θ − ⇒ V=G 1− 2
P2 (cos θ) 
2 2 R  M.R 
1
La expresión de Φ c = ω 2 (X 2 + Y 2 ) en esféricas:
2
1 2 2
X = R.senθ. cos λ Y = R.senθ.senλ ⇒ Φc = ω R sen 2 θ
2
C−A
A la expresión .la llamamos J 2 o factor de forma dinámica. Es un coeficiente ya que C-A
M.a 2
tiene dimensiones de momento de inercia, y M.a 2 sería el momento de inercia de una masa M sita
a la distancia a del eje de rotación. a = cte y es el semieje de la ecuatorial asignado a modelos
esféricos, o el radio de la menor esfera que envuelva toda la masa de la tierra; tangente, por tanto, a
ésta en el ecuador. Si la tierra fuese esférica A = C ⇒ J 2 = 0 .
C−A C−A M a2 
J2 = ⇒ = J 2a 2 ⇒ V=G 1 − J 2 2 P2 (cos θ)
M.a 2 M R  R 
El coeficiente J 2 se llama también coeficiente del segundo armónico zonal en el desarrollo del
potencial o también coeficiente de elipticidad geopotencial.
Este potencial esta limitado al coeficiente n = 2. Para hallar la expresión general del potencial
newtoniano V, buscamos a partir de la ecuación de Laplace en coordenadas esféricas. En general:
M 
n
a
V=G 1 − ∑∞
n =2   J P
n n (cos θ) 
R  R 
Vamos a imponer otra condición:
Suponemos la tierra simétrica respecto al plano ecuatorial. Vemos en tal caso que los términos
impares del desarrollo son nulos, quedando exclusivamente los términos pares. Consideraremos a
G.M
partir del término n = 2, porque como hemos visto V0 = y V1 = 0 .
R
M  1
2n
a M  a2 1 ω2 R 3 
W≈G 1 − ∑∞n =1   J 2n P2 n (cos θ) + ω 2 R 2 sen 2 θ = G 1 − J 2 2 P2 (cos θ) + sen 2 θ
R  R  2 R  R 2 G.M 
 2⋅a 3 v 2 (ω.a ) 2
Definimos otra constante m= donde = = ω 2 a es la aceleración centrífuga en un
G⋅M a a
Página 24
M
punto del ecuador (radio a) y G la aceleración newtoniana en dicho punto. Por tanto m es la
a2
razón entre la fuerza centrífuga y la atracción newtoniana sobre la unidad de masa en un punto del
ecuador, y se denomina constante geodinámica si se aplica a la tierra.
M a2 2 1 R3 
Si multiplicamos y dividimos por a 3 tenemos: W ≈ G 1 − J 2 2 P2 (cos θ) + m 3 sen θ
R  R 2 a 
Expresión desarrollada de W limitada a términos hasta n = 2.
2.10.- Problemas de contorno de la teoría del potencial. Solución de la ecuación

de Laplace en armónicos esféricos
2.10.1.- Teorema de Stockes:

Una función armónica en el exterior de una superficie S queda determinada de forma única por sus
valores de S, es decir, hay una sola función armónica V que toma sobre una superficie unos valores
de contorno dados. Este teorema es importante a la hora del cálculo del potencial gravitatorio
terrestre.
2.10.2.- Primer problema de contorno de la teoría del potencial o problema de Dirichlet:

Consiste en determinar la función armónica V a partir de los valores de contorno dados.
Solución de la ecuación de Laplace en coordenadas esféricas: ∆V = 0.
Para ello se requiere un cambio de coordenadas de cartesianas a esféricas sabiendo que:
X = R.senθ.senλ Y = R.senθ. cos λ Z = R. cos θ X 2 + Y 2 = R 2 .sen 2 θ
La expresión del laplaciano de V será:
1 ∂  2 ∂V  1 ∂  ∂V  1 ∂2V
∆V =  R  +  senθ  + =0
R 2 ∂R  ∂R  R 2 senθ ∂θ  ∂θ  R 2 sen 2 θ ∂λ2
Al igualar a cero tenemos la ecuación diferencial de Laplace que intentaremos resolver para hallar
la función potencial V, que dependerá de las tres variables V = V(R , θ, λ) .
Dadas las características de esta solución de V, podemos considerarla compuesta de una función
radial f(R) y de otra función Yn (θ, λ) ⇒ V(R , θ, λ) = f (R )Yn (θ, λ) , y a su vez esta última
función Yn (θ, λ) podemos ponerla como producto de otras dos funciones de tal manera que
V (R , θ, λ) = f (R )g (θ)h (λ ) . Hemos disgregado la ecuación diferencial inicial en tres ecuaciones
diferenciales, cada una de ellas en función de una sola variable:
h (λ) = sen (mλ ) o h (λ) = cos(mλ ) P mn (cos θ)sen (mλ)

Yn (θ, λ )
g(θ) = Pnm (cos θ) P mn (cos θ) cos(mλ)
Estos son los armónicos esféricos de superficie o funciones armónico-esféricas de superficie. Son
apropiadas para representar al potencial V en superficie dado que no dependen del radio R.
Multiplicadas por la función f(R) tenemos las funciones armónicas esféricas sólidas. Las funciones
f, g y h son funciones armónico-esféricas por tanto.
1
La función radial es del tipo: f (R ) = siendo n ≥ 0 .
R n +1
Según los valores de n y m tenemos:
n ≠ m tenemos los armónicos esféricos de superficie teserales (θ, λ ) .
n = m tenemos los armónicos esféricos de superficie sectoriales (θ, λ) .
m = 0 : sen (mλ ) = sen 0 = 0 m
Pn (cos θ) = Pn (cos θ)
Página 25
cos(mλ ) = cos 0 = 1
que son las funciones de Legendre o armónicos de superficie zonales. Vemos que éstas solo
dependen de θ , o sea, son los ceros del estudio de las funciones anteriores, o los puntos de cambio
de signo de dichas funciones.
Sabemos que la combinación lineal de funciones armónicas es también una función armónica, por
tanto, la solución V es una combinación lineal de estas: V = AV1 + BV2 + ... + NVn
donde las V son las funciones y el resto son los coeficientes. Entonces:
1
V = V( R , θ, λ ) = ∑∞
n =0 K n Yn (θ, λ )
R n +1
1
Aquí los Kn son los coeficientes y las Vn = Yn (θ, λ) son las funciones armónicas. También
R n +1
n +1
a
podemos escribirlo como: V = V(R , θ, λ) = ∑∞n =0   Yn (θ, λ ) .
R
n +1 n +1 n +1
a a a
Si R>a ⇒   <1 ⇒   Yn (θ, λ ) < Yn (θ, λ ) ⇒ ∑∞n =0   Yn (θ, λ) < V(a , θ, λ ) ,
R R R
lo que quiere decir que es convergente para puntos del exterior de la esfera de radio a, por lo que la
constante a representa el radio de la esfera en el exterior de la cual la función potencial V cumple la
ecuación de Laplace.
A su vez la función Yn (θ, λ) puede ponerse como combinación lineal de las funciones armónicas
P mn (cos θ)sen (mλ) y P mn (cos θ) cos(mλ) por tanto:
Yn (θ, λ ) = ∑nm =0 A m ( m m
n cos(mλ ) + B n sen (mλ )Pn (cos θ) )
La expresión más general como solución de la ecuación de Laplace es:
n +1
V = V( R , θ, λ ) = ∑∞ n 
n =0 ∑ m = 0 
a

R
(A m
n cos(mλ ) + B m m
n sen (mλ ) Pn (cos θ) )
1- Separamos el término correspondiente a n = 0 ⇒ m = 0:
0 +1
a
 
R
(( A 0
0
a
) a
cos 0 + B 00 sen 0) P00 (cos θ) =  A 0 P0 (cos θ) = A 0
R R
P0 (cos θ) = 1
n +1
a
V ( R , θ, λ) =  A 0 + ∑∞
R
n a
n =1 ∑m =0  
R
(A m
n cos(mλ) + B m m
n sen ( mλ )Pn (cos θ) )
2- Multiplicamos a ambos lados de la igualdad por R:
n +1
R.V( R , θ, λ ) = a.A 0 +

∑∞n =1 ∑nm =0 
a

R
(A m
n cos(mλ) + B m m
n sen (mλ )Pn (cos θ) )
3- Además sabemos que la función V que buscamos cumple las condiciones:
lim R →∞ V = 0 y lim R →∞ V = 0
GM
como lim R →∞ R.V = a.A 0 ⇒ a.A 0 = G.M ⇒ A0 =
a
n +1
V ( R , θ, λ) =
G.M
R
+ ∑∞ n a
n =1 ∑ m =0  
R
(A m
n cos(mλ ) + B m m
n sen ( mλ ) Pn (cos θ) )
G.M a.A m a.B m
4- Sacamos factor común a e introducimos los coeficientes n
= −J m
n y
n
= −K m
n :
R G.M G.M
Página 26
 
( )
n
a m
V ( R , θ, λ) = G.M 1 − ∑∞ ∑ n
n =1 m = 0   J n cos( mλ ) + K m
n sen ( mλ ) Pn
m
(cos θ) 
 R 
5- Separamos ahora en una sumatoria aparte los términos en m = 0 en el que solo tendremos los
términos J n y los armónicos esféricos zonales o polinomios de Legendre porque
Pnm (cos θ) = Pn (cos θ) , cos(mλ ) = cos 0 = 1 y sen (mλ ) = sen 0 = 0 :
G.M  
( )
n n
a n  a 
V ( R , θ, λ) = 1 − ∑∞
n =1   J P
n n (cos θ) − ∑ ∞
∑
n =1 m =1   J m
n cos( mλ ) + K m
n sen ( mλ ) Pn
m
(cos θ)  Si el
R  R R 
origen de coordenadas del sistema de referencia lo centramos en el centro de masas de la tierra, el
sumando correspondiente a n = 1 de la expresión del potencial es nulo (demostrado por otro lado
por sustitución del desarrollo del inverso de la distancia en la ecuación integral del potencial
newtoniano), por lo que nuestro desarrollo empezará en n = 2:
G.M  
( )
n n
∞  a  ∞ n  a 
V (R , θ, λ) = 1 − ∑n = 2   J n Pn (cos θ) − ∑n = 2 ∑m =1   J m m m
n cos( mλ ) + K n sen ( mλ ) Pn (cos θ)  Esta
R  R R 
es la expresión general adecuada para representar el potencial newtoniano en el exterior de la menor
esfera posible (aquella de radio a) que encierre en su totalidad a la masa que genera dicho potencial,
dado que para la deducción de esta expresión, hemos supuesto que la función cumplía la ecuación
de Laplace y esto sólo ocurre en el exterior de dicha esfera: por tanto en el caso de la Tierra, el
máximo valor de a corresponderá al semieje ecuatorial. Para la deducción de esta expresión, no ha
sido necesario tener en cuenta la. Distribución interna de masas.
A medida que aumenta el grado n, el valor de los coeficientes disminuye; así:
a−b 1
J 2 = 1082,6 × 10 −6 = 1,0826 × 10 −3 es del orden del aplanamiento α = = ≈ 3 × 10 −3
a 300
J 3 = −2,5 × 10 −6 es decir es unas 1000 veces menor
J 4 = −1,6 × 10 −6
Si suponemos que el cuerpo tiene simetría de revolución en torno al eje Z todos los términos en los
que intervienen los armónicos esféricos de superficie sectoriales y teserales (aquellos multiplicados
por los coeficientes J y K) no intervendrán en la expresión, es decir, se eliminan todos los términos
dependientes de λ :
G.M  
n
a
V (R , θ) = V = 1 − ∑∞
n =2   J n Pn (cos θ)
R  R 
Si ahora imponemos 1a condición de simetría respecto al plano ecuatorial, sólo nos aparecerán los
términos con n par. Por tanto:
GM
- Si aproximamos el potencial al primer término; es decir, V = el cuerpo que lo genera será por
R
tanto: una masa puntual, una esfera homogénea o una esfera con capas esféricas homogéneas.
- Si aproximamos la expresión al término en n = 2, las superficies equipotenciales
correspondientes tendrán las simetrías anteriormente expuestas (simetría de revolución en torno al
eje z y respecto al plano ecuatorial).
- Si incluimos ahora el término en n = 3, desaparecen las simetrías respecto al plano ecuatorial,
como es el caso de la Tierra.
Cuanto más nos alejemos del cuerpo, menor será la influencia de los términos de mayor grado; de
hecho, lim R →∞ V = 0 , lo cual implica que los términos van decreciendo con la distancia, por ello en
dinámica planetaria en que R no es infinito, pero sí muy grande, se aproxima el potencial a su forma
Página 27
G.M
más sencilla del resto de los términos, V = , prescindiendo del resto de los términos.
R
Cuanto más cerca estemos del cuerpo, mayor influencia tienen los términos sectoriales y teserales,
éste es el caso de los satélites, cuyas órbitas requieren de expresiones más complejas. Si la tierra
fuese puntual, sus órbitas cumplirían las leyes de Kepler. Gracias al seguimiento de las órbitas
'reales' de los satélites, se han deducido los valores de un gran número de estos términos.
Página 28
3.- EL POTENCIAL Y LA GRAVEDAD NORMAL. SISTEMAS

GEODÉSICOS DE REFERENCIA
3.1.- La función geopotencial

Sabemos que las superficies equipotenciales son aquellas en las que el potencial gravitatorio
terrestre es constante: W=cte.
Queremos determinar la expresión de W, para, una vez conocida dicha función, poder analizar las
características de las superficies de ese campo. A medida que nos alejamos del cuerpo, la forma de
las superficies equipotenciales se va regularizando, y a gran distancia se van asemejando cada vez
más a esferas. Pero no existe una expresión analítica para dar el potencial correspondiente al geoide,
que se define como una superficie, lugar geométrico, en el que la constante es igual al valor que
tendría dicho potencial al nivel medio de los mares en calma. Ese geoide tiene puntos fuera de las
masas atrayentes y también puntos dentro de la tierra, en los que se producirán discontinuidades en
las derivadas de los potenciales, por lo que en dichos puntos no se cumplirá la ecuación de Laplace.
Función geopotencial o potencial gravitatorio terrestre W:
G.M  
( )
n n
a n  a 
W= 1 − ∑∞
n =2   J P
n n (cos θ) − ∞
∑n = 2 ∑m =1   Jm m m
n cos( mλ ) + K n sen ( mλ ) Pn (cos θ)  + φ c
R  R R 
Pero, aunque conociese todos los términos correspondientes a esta expresión, el geoide definido
utilizando esta función W = Wo, contendría toda la masa del cuerpo, dado que la expresión de W es
cierta fuera de las masas atrayentes, por tanto al geoide definido de esta manera vamos a llamarle
geoide geopotencial.
Vistos estos inconvenientes, no podemos abordar el problema de definir el geoide utilizando una
expresión de W con infinitos términos de los cuales no conocemos los coeficientes. Para nuestro
estudió vamos a suponer que:
• La tierra tiene simetría de revolución respecto al llamado eje medio de rotación terrestre.
• El origen de coordenadas está en el centro de masas de la tierra.
• existe simetría respecto al plano ecuatorial.
Por lo tanto, nos quedamos con la expresión simplificada del potencial newtoniano, en la que
aparecen sólo los términos zonales pares (n = par), y la llamamos función esferopotencial o
potencial normal U(R , θ) = Vaprox + φ c donde:
G.M  
n
a 1 1
Vaprox = V (R , θ) = V = 1 − ∑∞
n =2   J n Pn (cos θ)  y φ c = ω 2 ( X 2 + Y 2 ) = ω 2 R 2 sen 2 θ
R  R  2 2
G.M   1
n
∞  a  1
U ( R , θ) = 1 − ∑ n = 2   J n Pn (cos θ)  + ω 2 (X 2 + Y 2 ) = ω 2 R 2 sen 2 θ con n par.
R  R  2 2
Las superficies equipotenciales U = cte se aproximan a la esfera más que al geoide, y son los
llamados esferoides. A la superficie U = U 0 = W0 la llamamos esferoide de nivel, y es la superficie
Página 29
equipotencial que más se aproxima al geoide cumpliendo las condiciones impuestas. Uo es

numéricamente el mismo valor que Wo, pero son superficies distintas, las funciones que
representan son diferentes.
El potencial gravitatorio terrestre será por tanto: W (R , θ, λ) = U(R , θ) + T(R , θ, λ) donde T recibe el
nombre de potencial perturbador, y contiene los términos de V(R , θ, λ) .
3.2.- Esferoide de Bruns: teoría del primer orden de aproximación.

Si de la expresión de la función esfera potencial nos quedamos únicamente con los términos hasta n
= 2, o sea, nos limitamos al primer orden de aproximación, obtengo una función que representa al
llamado esferoide de Bruns:
G.M   a   1
2 4
a
U ( R , θ) = 1 −   J 2 P2 (cos θ) −   J 4 P4 (cos θ) − ... + ω 2 R 2 sen 2 θ
R   R  R  2
G.M   a   1
2
U ( R , θ) = 1 −   J 2 P2 (cos θ)  + ω 2 R 2 sen 2 θ
R   R   2
Analizaremos ahora la forma de la superficie equipotencial U = Uo = Wo o esferoide de nivel.
Demostraremos que, el esferoide de Bruns es un elipsoide de revolución. Para ello deduciremos la
expresión del radio vector R en polares.
GM
Sacando factor común puedo introducir el término centrífugo en el corchete:
R
G.M   a  
2
1 ω2 R 3
U ( R , θ) = 1 −   J 2 P2 (cos θ) + sen 2 θ
R   R  2 GM 
 2⋅a 3
Introducimos ahora la constante geodinámica m= :
G⋅M
G.M   a  
2
1 R3
U ( R , θ) = 1 −   J 2 P2 (cos θ) + m 3 sen 2 θ
R   R  2 a 
3.3.- Primera fórmula de Clairaut en aproximación de primer orden

Determinación de Uo: primera fórmula de Clairaut.
Página 30
3.3.1.- 1ª método:
La función U(R , θ) representa infinitas superficies equipotenciales. Sabemos que la función tomará
el valor de Uo = Wo (esferoide de nivel) por ejemplo en los puntos ecuatoriales, es decir, para R = a
π
y θ= en los puntos polares donde R = b y θ = 0 , θ = π .
2
Hallemos el valor Uo en función de constantes en un punto polar (R = b, θ = 0 ):
3 1 3 1
Sabiendo que P2 (cos θ) = cos 2 θ − ⇒ P2 (cos θ) = − = 1 y que senθ = 0
2 2 2 2
G.M   a  
2
U0 = −
   J2 
1
b   b  
π
Hallamos el valor Uo en función de constantes en un punto ecuatorial (R = a, θ = )
2
3 1 π 1
Sabiendo que P2 (cos θ) = cos 2 θ − ⇒ P2 (cos ) = − :
2 2 2 2
G.M  J 2 m 
U0 = 1+ +  que depende de la masa del cuerpo (M), de la velocidad angular de rotación
a  2 2
(incluida en m) y de datos geométricos (a, J 2 ).
4- Igualando este último valor de Uo a U(R , θ) tenemos:
G.M  J 2 m  G.M   a  
2
1 R3 2
1 + +  = 1 −   J 2 P2 (cos θ) + m 3 sen θ
a  2 2 R   R  2 a 
a2
5- En la teoría de la primera aproximación ≈1:
R2
[  
][
3 1 1
cos 2 θ −  msen2 θ
] 
1−J 2 −1
R 2 2 2 3 1 1 J m
= = 1−J 2  cos 2 θ−  msen2 θ 1 2 
a J2 m 2 2 2 2 2
1 
2 2
−1
J m  J m
6- Llamando x = 2 + tenemos que el término 1 + 2 +  es de la forma (1 + x ) n . Los
2 2  2 2
términos J 2 y m son del orden de 10 −3 , es decir, mucho menores que l, x será también un valor
pequeño, por lo que podemos aproximar el valor del paréntesis al de su desarrollo en serie de
Taylor:
m( n − 1) 2
(1 + x ) n = 1 + nx + x + ... pero, en teoría de primer orden de aproximación los términos
2!
J 22 , m 2 , J 2 m , se desprecian, por lo tanto los términos de x 2 en adelante serán despreciables
−1
 J m  J2 m 
1 + 2 +  ≈ (1 − x ) = 1 − − 
 2 2  2 2
R  3 1 1  J m
= 1 − J 2 ( cos 2 θ − ) + msen 2 θ 1 − 2 − 
a  2 2 2  2 2
7- Operando y despreciando los términos mencionados obtenemos:
R  3 1 1 J m
= 1 − J 2 ( cos 2 θ − ) + msen 2 θ − 2 − 
a  2 2 2 2 2
Página 31
como sen 2 θ = (1 − cos 2 θ) ⇒ R

a [ 3 1 J2 J2 m m
= 1− J 2 cos 2 θ− mcos 2 θ −  −
2 2 2 2 2 2 ]
R 3 1 3 m
=1− J 2 cos 2 θ− m cos 2 θ=1−cos 2 θ  J 2  
a 2 2 2 2
[ 2 3
8- El radio vector de esferoide en polares es: R=a 1−cos θ  J 2 
2
m
2 ]
R = f (θ) Queda demostrado que el esferoide de Bruns es un elipsoide de revolución.
3.3.2.- 2ª método:
x 2 + y2 z2
1- La ecuación del elipsoide de revolución en cartesianas es + = 1.
a2 b2
a−b
α= ⇒ b = a (1 − α )
a
Usamos la relación entre cartesianas y esféricas:
x=r⋅sen⋅cos , y=r⋅sen⋅sen y z=r⋅cos 
x2y2=r2⋅sen2  sen2 cos2 =r2⋅sen2 
r 2⋅sen2  r 2⋅cos 2  r 2⋅sen2  r 2⋅cos 2 
 =1 ⇒ 2
 2
=1 , Aproximamos 1−2≈1−2⋅ :
a2 b2 a a⋅1−
r2
a2  sen2

cos2 
1−2⋅ 
=1 [1], extrayendo lo que tenemos entre paréntesis:
cos2  1−2⋅sen2 cos2 

sen2  = ⇒
1−2⋅ 1−2⋅
2 cos2  −2⋅⋅sen2 sen2 cos 2 
sen  = ⇒
1−2⋅ 1−2⋅
cos2 
sen2  =−2⋅⋅sen2 11−2⋅−1 , quitando el término 2 y usando el desarrollo
1−2⋅
en serie 1xn≈1n⋅x en primera aproximación ⇒ tenemos:
2
cos 
sen2  ≈−2⋅⋅sen2 12⋅ ⇒
1−2⋅
2 cos2 
sen  ≈2⋅⋅1−sen2 1≈2⋅⋅cos 2 1
1−2⋅
sustituyendo en [1] y hallando su raíz cuadrada:
r

1
−
2⋅⋅cos 2 1=1 ⇒ 2
r=a12⋅⋅cos  2
a
Como 1xn≈1n⋅x en el desarrollo en serie en primera aproximación ⇒
r=a1−⋅cos2 
3 m
Por tanto por lo deducido en la 1ª forma tenemos que = ⋅J2 
2 2
El radio-vector de un elipsoide de revolución en teoría de 2º orden de aproximación tiene por
expresión:
3
R = a (1 − α cos 2 θ − α 2 sen 2 2θ)
8
Página 32
3.4.- Gravedad normal en el esferoide de Bruns o elipsoide de nivel:

→
Definimos la gravedad normal γ = gradU como aquella que deriva del potencial normal U:
G.M   a   1
2
U ( R , θ) = 1 −   J 2 P2 (cos θ)  + ω 2 R 2 sen 2 θ
R   R   2
→ ∂U ∂U
γ =− ≈− donde el signo es debido a que la dirección del vector normal es contraria a la del
∂n ∂R
vector gravedad.
→ ∂U ∂U
Esta aproximación es válida en teoría de primer orden. En realidad γ = − ≈− cos(ϕ'−ϕ) pero
∂n ∂R
el máximo valor de (ϕ'−ϕ) = 11. ´6 cuando ϕ = 45º . El error que cometemos al hacer dicha
aproximación es el causado al aproximar cos 11,6 '=0. 9999995≈1 . También hay que tener en
cuenta que al realizar esta aproximación, estamos despreciando las diferencias entre la utilización
de la latitud egocéntrica y la astronómica en las expresiones que vamos a deducir.
Hallamos cuanto vale la gravedad normal partiendo de:
[ ]
2
G.M a 3 1 1
U  R , θ = −GM 3 J 2  cos 2 θ −  ω 2 R 2 sen 2 θ
R R 2 2 2
Derivando respecto a R:
∂U GM GMa 2  3 1
=− 2 +3 J 2  cos 2 θ −  + ω 2 Rsen 2 θ
∂R R R 4
2 2
γ =−
∂U GM
∂R
a2
[
= 2 1−3 2 J 2
R R
3
2 
1
cos 2 θ− −
2
ω2 R 3
GM
sen 2 θ ] ; y como m =
ω2 a 3
GM
GM  a2 3 2 1 R3 2  a2 a3
γ= 1 − 3 J 2  cos θ −  − m sen θ  ; aproximando ≈ 1, ≈1
R 2  R2  2 2 a3  R2 R3
GM  3 1 
γ= 2 
1 − 3J 2  cos 2 θ −  − msen 2 θ = γ (θ)
R   2 2  
π
Para un punto en el ecuador R = a y θ = , en base a lo que resolvimos con anterioridad nos vamos
2
GM  3 
a encontrar que la gravedad normal en el ecuador es: γ E = 2 1 + J 2 − m  .
a  2 
GM
Para un punto en el polo R = b y θ = 0 , θ = π , será: γ P = [1−3J 2 ] .
b2
Deducción de la expresión del aplanamiento gravimétrico γ = γ E (1 + β cos 2 θ) :
Partimos de la fórmula:
Página 33
GM  3 1 
γ= 2 
1 − 3J 2  cos 2 θ −  − msen 2 θ
R  2 2 
GM  3 2 1 2 
Hacemos R=a  1−α cos θ 
2
⇒ γ= 2 2 1 − 3J 2  2 cos θ − 2  − msen θ
2
a (1 − α cos θ)    
GM  3 1 
γ= 2 
1 − 3J 2  cos 2 θ −  − msen 2 θ (1 − α cos 2 θ) −2
a  2 2 
Aplicando Taylor hacemos que (1 − α cos 2 θ) −2 = (1 + 2α cos 2 θ) :
GM  3 1 
γ= 2 
1 − 3J 2  cos 2 θ −  − msen 2 θ (1 + 2α cos 2 θ)
a  2 2 
GM  9 3 
γ = 2 1 − J 2 cos 2 θ + J 2 − msen 2 θ (1 + 2α cos 2 θ)
a  2 2 
Al hacer el producto despreciamos todos aquellos términos en α 2 , J 2 m, J 2 α, mα,... :
GM  9 3 
γ= 2 
1 − J 2 cos 2 θ + J 2 − msen 2 θ + 2α cos 2 θ
a  2 2 
Poniendo el seno en función del coseno tendremos:
GM  9 3 
γ= 2 
1 − J 2 cos 2 θ + J 2 − m + m cos 2 θ + 2α cos 2 θ
a  2 2 
GM  3  9  
γ = 2 1 + J 2 − m +  − J 2 + m + 2α  cos 2 θ
a  2  2  
 3  9  2 
1 + 2 J 2 − m +  − 2 J 2 + m + 2α  cos θ −1
γ      3  9  2  3 
= = 1 + J 2 − m +  − J 2 + m + 2α  cos θ1 + J 2 − m 
γE 3  2  2   2 
1+ J2 − m
2
−1 −1 1
 3   3   3 
Mediante Taylor: 1 + J 2 − m  = (1 + x ) n
⇒ 1 + J 2 − m  ≈ 1 − J 2 + m 
 2   2   2 
Despreciamos todos aquellos términos en α 2 , J 2 m, J 2 α, mα,... :
γ  3  9  3 
= 1 + J 2 − m +  − J 2 + m + 2α  cos 2 θ − J 2 + m 
γ E  2  2  2 
9 γ 2
Haciendo β = − J 2 + m + 2α obtenemos que γ = 1 + β cos θ
2 E
Segunda fórmula de Clairaut o fórmula de distribución de la gravedad normal:
9
β = − J 2 + m + 2α siendo β el aplanamiento gravimétrico.
2
γP − γE
Si θ = 0 ⇒ γ c = γ P y γ P = γ E (1 + β) ⇒ β= expresión parecida a la del
γE
aplanamiento geométrico.
Como consecuencia de la 1ª y 2ª fórmulas de Clairaut obtenemos la tercera fórmula de Clairaut o
teorema de Clairaut en teoría de primer orden de aproximación.:
3 m 9 3 3 9
α= J2 + ⇒ 3α = J2 + m ⇒ m − 3α = − J 2
2 2 2 2 2 2
9 3 5 5
β = − J 2 + m + 2α ⇒ β = m − 3α + 2α + m = m − α ⇒ β+α= m
2 2 2 2
Página 34
3.5.- Esferoide de Helmert:

La superficie equipotencial correspondiente a este esferoide se obtiene igualando la función del
potencial normal a la constante Uo = Wo, pero considerando la función potencial en aproximación
de segundo orden, es decir, tomando hasta el término en J4.
G.M   a   1
2 4
a
U ( R , θ) = 1 −   J 2 P2 (cos θ) −   J 4 P4 (cos θ) − ... + ω 2 R 2 sen 2 θ
R   R  R  2
Si introducimos en la expresión los valores empíricos de J2 y J4 obtenidos gracias a observación del
movimiento de los satélites, vemos que este esferoide no es exactamente elipsoide de revolución,
aunque sí muy próximo.
→
El valor deducido en este caso para la gravedad normal a partir de γ = gradU será:
γ = γ E (1 + βsen 2 ϕ − β1sen 2 2ϕ)
Al estar trabajando ya con expresiones correspondientes a la teoría de segundo orden hemos de
especificar que la latitud empleada es la geocéntrica.
3.6.- El elipsoide de nivel

Durante algún tiempo se estuvieron proponiendo como figura normal de la tierra diferentes
esferoides.
En la reunión de Estocolmo de 1930 se adopta la fórmula internacional de la gravedad normal.
Cassini (basándose en el elipsoide de Hayford) propone los valores de β y de β1 . Acuerdan que
siendo despreciables las diferencias entre el elipsoide y los esferoides, es preferible adoptar una
superficie que tenga un tratamiento matemático sencillo, y que sea la misma que ha venido
utilizando la geodesia pata cálculos geométricos.
A partir de entonces, y en la actualidad, se adopta como figura normal de la tierra la de un elipsoide
de revolución de semiejes a y b, que tenga la misma masa que la tierra, que rote en torno a su eje
menor con la misma velocidad angular que ésta, y que sea una superficie equipotencial del campo
de la gravedad normal. Se desprecian las pequeñas diferencias encontradas entre la figura real y este
elipsoide.
Este elipsoide de nivel U = Uo = Wo sirve tanto de superficie de referencia geométrica (para los
cálculos de la geodesia geométrica) como de superficie de referencia física (y las fórmulas de la
gravedad normal se desarrollan sobre ese elipsoide).
3.7.- Teoría de Pizzetti y Somigliana. Fórmula internacional de la gravedad

Utilizando teoría de segundo orden, a partir de la expresión general del potencial normal y
→
sabiendo que γ = gradU , se deducía la fórmula trinomia de la gravedad normal en aproximación de
2° orden, válida para el elipsoide y los esferoides:
γ = γ E (1 + βsen 2 ϕ − β1sen 2 2ϕ) donde ϕ = ϕ geocéntrica
En el caso de los distintos elipsoides que han ido adoptándose como mejor aproximación a la figura
de la tierra, el valor del coeficiente β1 se deduce del adoptado para α y para β es un valor forzado
para que la figura representada para la función U sea un elipsoide de revolución. No es éste el caso
del esferoide en el que, aun cumpliendo unas condiciones de simetría, no es un elipsoide.
A medida que se adoptan nuevos elipsoides (se mejora la aproximación, por lo que se adoptan
valores más precisos para los parámetros; cambia la fórmula de γ ), hay que descartar las fórmulas
anteriores y utilizar la nueva. Pero existe mucha información cartográfica, muchos datos y
anomalías de gravedad calculados con fórmulas de gravedad normal de sistemas anteriores, y hay
Página 35
que convertirlos.
Somigliana, Pizzetti y distintos geodestas desarrollaron toda la teoría del elipsoide y dedujeron
teóricamente la fórmula de la gravedad normal sobre un elipsoide de nivel. La fórmula del potencial
normal es análoga a la empleada para los esferoides. La diferencia está en que en el caso del
elipsoide, los J no pueden tomar cualquier valor, son consecuencia de desarrollar la teoría del
elipsoide (son valores deducidos teóricamente); definidos en función de J , y se conocen una vez
fijado éste. Son valores forzados para que U sea un elipsoide perfecto de revolución. Son valores
ligeramente distintos de los empíricos obtenidos a partir de! conocimiento de las perturbaciones de
las órbitas de los satélites, que eran los usados para los esferoides.
Por tanto W representa al geoide, y U representa al elipsoide de nivel. Si el valor de U es conocido,
el elipsoide estará perfectamente determinado. Este U dependerá de G, de ω y de las dimensiones
geométricas que tenga dicho elipsoide.
Son razones históricas y prácticas las que imponen la adopción del elipsoide de revolución como
superficie de referencia tanto geométrica como física. Se adopta el elipsoide como figura normal de
la tierra. La superficie de ese elipsoide es una superficie equipotencial del campo de la gravedad
normal.
3.8.- Ventajas del uso del elipsoide como figura de referencia:

• El elipsoide se venía utilizando desde hacia muchos años (más de 300) como superficie de
referencia para todas las operaciones geodésicas.
• La geometría del elipsoide y su trigonometría han llevado a fórmulas prácticas para los
problemas directo e inverso de cambio de coordenadas. Desde hace tiempo utilizadas por los
servicios geodésicos de todo el mundo.
• A igualdad de semiejes, las diferencias entre el elipsoide de revolución y los diferentes
esferoides de revolución es muy pequeña.
• La teoría del campo gravitatorio del elipsoide desarrollada por Pizzetti y Somigliana (entre
otros) da la función U y el valor de la constante U en función de la masa del elipsoide, la
velocidad angular de rotación, y los dos parámetros geométricos para definir su forma.
Igualando M a la masa de la tierra, y ω a la velocidad angular de rotación, habremos de fijar los
parámetros geométricos del elipsoide para que la diferencia entre W y U sea mínima; es decir, el
potencial perturbador T = W sea mínimo.
3.9.- Gravedad normal en el elipsoide:

La fórmula rigurosa para definir la gravedad normal en puntos del elipsoide es la fórmula cerrada
aγ E cos 2 ϕ + bγ P sen 2 ϕ
propuesta por Somigliana en 1929: γ = γ (ϕ) =
a 2 cos 2 ϕ + b 2 sen 2 ϕ
ϕ - latitud rigurosamente egocéntrica, aproximadamente geográfica o geodésica.
Conociendo los valores de γ P y γ E (definidos a su vez en función de parámetros del elipsoide, es
decir, se fijan) y el de la latitud en el punto ϕ , podemos obtener directamente el valor de la
gravedad normal en dicho punto.
Otra forma de expresar γ es haciendo cos 2 ϕ = 1 − sen 2 ϕ , dividiendo el numerador y el denominador
por el coeficiente a y sacando factor común γ E :
Página 36
1
 bγ P −aγ E sen 2 ϕ
aγ E
γ =γ E
bγ − aγ E
k= P
 1−
a 2 −b 2 
a 2
sen 2 ϕ
aγ E 1ksen ϕ
2
γ =γ E
e2 =
a 2 − b2
=
E2 1−e 2 sen 2 ϕ
a2 a2
Desde el punto de vista práctico, es costumbre utilizar fórmulas desarrolladas deducidas
1
matemáticamente (basándonos en el desarrollo de la expresión ) y a partir de la anterior en
1− x
series convencionales, con dos términos (fórmula binomia en la aproximación de primer orden), con
tres términos (fórmula trinomia en aproximación de 2º orden), ...
γ = γ E (1 + βsen 2 ϕ − β1sen 2 2ϕ − β 2 sen 2 ϕsen 2 2ϕ...)
1 1
siendo β1 = α(α + 2β) , β 2 = α 2 (α + 3β + ...) expresiones ya deducidas a partir de α y de β .
8 8
3.10.- Sistema Geodésico de Referencia

Un sistema geodésico de referencia está formado por cuatro parámetros por los que queda
determinado el elipsoide equipotencial de referencia y su campo de gravedad normal, visto
anteriormente, en el que operando adecuadamente llegábamos a la fórmula cerrada de Somigliana.
El campo de este elipsoide equipotencial queda definido al fijar cuatro parámetros. Estos
parámetros son diferentes para los distintos sistemas de referencia adoptados.
3.10.1.- Sistema Geodésico de Referencia 24-30 (1924-1930):

Basado en el elipsoide de Hayford, que ya se utilizaba desde 1909, y fue adoptado en la Asamblea
de Madrid de 1924 como oficial para los trabajos geodésicos geométricos.
Se fijaron dos parámetros para definir dimensionalmente el elipsoide:
a = 6378388 m Semieje mayor del elipsoide.
1
α= Aplanamiento geométrico.
297
Y otros dos parámetros para determinar el campo de gravedad normal asociado a dicho elipsoide:
γ E = 978,049 gales (Heiskanen, 1928). Gravedad normal en el ecuador.
rad
ω = 7,292115l x 10 −5 Velocidad angular de rotación terrestre.
sg
Como consecuencia, la gravedad normal en función de la latitud, limitando el desarrollo hasta
términos de 2º orden: γ = γ E (1 + βsen 2 ϕ − β1sen 2 2ϕ)
En la Asamblea de Estocolmo de l930 de la Unión Geofísica y Geodésica Internacional se fijaron
los valores de β y β1 , y se acordó adoptar la siguiente fórmula para el calculo de la gravedad
normal; fórmula internacional de la gravedad conocida también como fórmula de Cassini-Silva. Es
el único sistema en que se fija el valor de la gravedad normal, en los demás se fija el GM.
γ 30 = 978,049(1 + 0,0052884sen 2 ϕ − 0,0000059sen 2 2ϕ)
Se fijó el valor de la gravedad en el ecuador, que fue deducido mediante un proceso de ajuste a
partir de medidas de gravedad convenientemente reducidas y distribuidas por distintos lugares de la
Página 37
superficie terrestre.
Los valores de β y β1 fueron deducidos a partir del teorema de Clairaut modificado por Somigliana
al 2° orden de aproximación en que:
5 17 1 2 13
α+β= mχ( α ) donde χ(α) = 1 − α− α − α 3 ... es la llamada función de Clairaut.
2 35 245 18865
Fijado el valor de α tendremos determinado el valor de χ( α ) . Limitando al 2° orden de
17 5  17 
aproximación, χ( α ) = 1 − α , ya que del producto m1 − α  sale un término en m α , ya de 2°
35 2  35 
orden. Se aproxima, en la práctica, al segundo orden porque nos da la suficiente precisión para la
mayoría de los trabajos que necesitemos abordar. Estaría de más utilizar una fórmula con precisión
del miligal si vamos a medir y utilizar anomalías de gravedad de menor grado de precisión.
ω2 a
El valor de m utilizado para calcular estos coeficientes es: m =
γE
Por tanto, fijados los valores de a, α , γ E y ω obtenemos los valores numéricos de β y β1
siguientes:
χ( α ) = 0.998364552
α + β = 0,008655388 ⇒ β = 0,005288384
1
β1 = α(α + 2β) ⇒ β1 = 0,000005869
8
En tercer orden de aproximación:
γ = γ E (1 + βsen 2 ϕ − β1sen 2 2ϕ − β 2 sen 2 ϕsen 2 2ϕ...)
1 1
β 2 = α 2 (2α + 3β) − α 2 (3α + 4β) ⇒ β 2 = 0,000000022
8 32
3.10.2.- Sistema Geodésico de Referencia 1967

Basado también en la teoría del elipsoide de Hayford, fue adoptado en la Asamblea de Lucerna de
1967. Se fijaron los siguientes parámetros:
a=6378160 metros
m3
GM = 398603 × 10 9 Constante de gravitación egocéntrica
sg 2
J 2 = 10827 × 10 −7 Factor de elipticidad geopotencial o coeficiente del 2º armónico zonal
rad
ω = 7,2921151467 × 10 −5 Velocidad angular de rotación terrestre
sg
Una vez fijados estos valores, y basándose en la teoría, se pueden deducir las valores, tanto de las
constantes geométricas y físicas, como de los coeficientes de la fórmula de la gravedad normal,
gravedad en el ecuador, ...
Los valores de J4, J6, J8,... en este caso no son los empíricos deducidos de las perturbaciones de las
órbitas de los satélites, dado que si introdujésemos dichos valores en la fórmula del potencial no
obtendríamos un elipsoide de revolución sino esferoides; por tanto, son valores forzados deducidos
por teoría a partir de J2 para obtener el elipsoide equipotencial
Fórmula práctica:
γ 67 = 978,0318(1 + 0,0053024sen 2 ϕ − 0,0000059sen 2 2ϕ) gales
Se acuerda disminuir el valor de la gravedad absoluta en Postdam (origen de la red gravimétrica) en
14 miligales, pasando de 98.1,274 a 981,260 gales. Como consecuencia, han de ser modificados en
Página 38
la misma medida todos los valores de 1ª gravedad absoluta deducidos por diferencias de gravedad a
partir de dicho origen; es decir, arrastrados desde Postdam. Esto afectará a las anomalías, y se
utilizan fórmulas de conversión.
3.10.3.- Sistema Geodésico de Referencia de 1980:

En la 19° Asamblea celebrada en Gamberra en 1979, se acuerda cambiar los valores de los
parámetros aceptados en el 67, dado que no representa a la figura de la tierra con la precisión
requerida.
a = 6378137 metros
J 2 = 108263 x 108
GM = 3986005 × 10 8 m 3 sg 2
 = 7292115 x 10 −11 rad sg
3.10.4.- Sistema Geodésico de Referencia WGS84:

Conserva exactamente los valores anteriores, pero en este caso en lugar de J2 se da C2 = C2,0(n =
2,m = 0) que es el coeficiente armónico zonal normalizado de 2º grado:
J2
C 2=−
5
3.10.5.- Conversión de anomalías gravimétricas por cambio de sistema geodésico de
referencia:
Anomalía gravimétrica es ∆g = g P − γ Q siendo g P la gravedad real reducida al geoide y γ Q la
gravedad normal correspondiente a la proyección de Q en la dirección de la normal al elipsoide de
nivel del punto P. Es la diferencia entre la gravedad real reducida al geoide menos la gravedad
normal en el elipsoide.
Hemos de deducir la gravedad reducida al geoide que será igual a una gravedad observada en
superficie terrestre y corregida: g reducida = g observada + correcciones
Página 39
n: normal al geoide
n': normal al elipsoide
ε̂ : desviación de la vertical
N: ondulación del geoide
Anomalía de la gravedad:
∆g = g reducida − γ normal
Debido a que la fórmula de la gravedad normal varía de un sistema de referencia a otro, también
variarán los valores de las anomalías de gravedad.
Las gravedades observadas no proceden de una observación de gravedad absoluta, sino de un enlace
por incrementos de gravedad con la absoluta conocida de un punto de la red gravimétrica de
referencia (datum Postdam).
Página 40
4.- EL CAMPO DE LA GRAVEDAD ANÓMALO. DETERMINACIÓN DEL

GEOIDE
4.1.- Potencial normal y potencial perturbador

W (r, θ, λ ) = U(r, θ) + T (r, θ, λ )
W ( x, y, z) = U( x, y, z) + T ( x , y, z)
Esta misma fórmula podríamos haberla escrito en cartesianas refiriéndonos a la dependencia de los
tres tipos de potenciales a x,y,z. En cualquier caso, vemos que son funciones de punto; es decir,
podemos establecer dicha relación particularizada para cada punto del espacio.
4.2.- Conceptos de anomalía de gravedad, ondulación del geoide y desviación de

la vertical
También llamadas magnitudes del campo anómalo:

→ → →
- Vector anomalía de gravedad: ∆g = g P − γ Q que difieren en el módulo y en dirección formando
un ángulo ε̂ .
- Anomalía de la gravedad: ∆g = g P − γ Q estas cantidades escalares son las que vamos a tratar de
deducir a partir de los valores de gravedad observados. Es función de los módulos de g P y γ Q . No
→
es el módulo del ∆g , dado que g P y γ Q no están en la misma dirección, pero la diferencia entre
→
∆g y ∆g es despreciable. Las anomalías de gravedad son debidas a la "anómala" distribución de la
masa interna terrestre, principalmente en la corteza, dado que en capas mas profundas la masa está
distribuida de forma más regular, no homogénea, pero si en capas de diversas densidades.
Página 41
- Desviación de la vertical: ε̂ Es la diferencia entre las direcciones de las normales al geoide y al

elipsoide en el punto P. Es el ángulo que se desvía la normal al geoide n respecto a la normal
elipsódica n’ en ese punto, y sus dos componentes NS y EW nos determinan la orientación del
geoide respecto al elipsoide en el punto P. Puede ser conocido geométricamente y a partir de
anomalías de gravedad.
- Ondulación del geoide: N Es la distancia entre P y Q medida a lo largo de la normal elipsoidal.
Es la altitud h del geoide sobre el elipsoide en cada punto. Puede ser-conocido a partir de anomalías
de gravedad.
Para situar el geoide respecto al elipsoide necesitaremos conocer su separación N y su orientación
respecto a éste. La prospección gravimétrica (que junto con la prospección sísmica y magnética
forma parte de la geofísica aplicada) consiste en calcular anomalías para deducir estos parámetros.
Pero para la orientación son necesarias una componente vertical NS y una componente horizontal
EW, por ello la desviación de la vertical ε̂ viene definida por dichas componentes:
NS=ξ= 90− G − 90− A = astronómica − geodésica

EW = η = (λ astronómica − λ geodésica ) cos ϕ
· Aplicando el teorema de los senos en el triángulo PN CZ G :
senη ≈ η
sen90 sen (λ A − λ G )
= ⇒ sen (λ A − λ G ) ≈ λ A − λ G η = (λ A − λ G ) cos ϕ G
sen (90 − ϕ G ) senη
sen (90 − ϕ G ) = cos ϕ G
Consideramos  A =G
Página 42
· Aplicando el teorema de los senos en la anterior figura obtenemos:

sen (λ A − λ G ) ≈ λ A − λ G
π
sen 1 λ − λG
2 = sen (λ A − λ G ) ⇒ sen (A A − A G ) ≈ A A − A G = A
senϕ A A − A G
senϕ G sen (A A − A G )
senϕ G ≈ senϕ A
1 cos ϕ
λA − λG = (A A − A G ) ⇒ η = (A A − A G ) ⇒ η = (A A − A G ) cot gϕ
senϕ senϕ
El método anterior es geométrico. A partir de coordenadas calculamos las componentes de la
desviación de la vertical.
Partiendo de anomalías de gravedad como datos observables, las fórmulas de Vening Meinesz nos
dan otro procedimiento para calcular componentes de la desviación de la vertical, y la fórmula de
Stokes nos permitirá conocer el valor de N. Esta es la finalidad de la gravimetría.
4.3.- Fórmula de Bruns.

La fórmula de Bruns expresa la relación entre la ondulación del geoide y el potencial perturbador.
En el punto P: Wp = U p + Tp La variación de la función U desde el punto Q al P será
∂U ∂U
Up = UQ + dn ' pero la separación entre Q y P es dn’ = N, por tanto: U p = U Q + N.
∂n ' ∂n '
→ →
Sabemos que: g = gradW y γ = gradU
Sus componentes escalares en las direcciones normales n y n’ respectivamente son:
∂W ∂U
g=− y γ=−
∂n ∂n '
Sustituyendo: p = U Q − γ Q N
U ⇒ Wp = U Q − γ Q N + Tp
Tp
Como Wp = Wo y U Q = U o = Wo ⇒ W0 = W0 − γ Q N + Tp ⇒ γ Q N = Tp ⇒ N =
γQ
Fórmula de Bruns que nos permite calcular N una vez conocido el potencial perturbador T en el
punto P del geoide y calculada la gravedad normal en su proyección Q.
4.4.- Fórmula o ecuación fundamental de la geodesia física:

Expresa la relación entre la anomalía y el potencial perturbador.
Llamamos:
→ → →
• ∂g = g p − γ p Vector perturbación de gravedad. Está definido sólo en P.
• ∂g = g p − γ p Perturbación de gravedad. Cantidad escalar.

→
∂g = gradW − gradU = grad( W − U ) = gradT
→ ∂T
∂g = gradT y su componente escalar en la dirección de la normal será: ∂g = − ∂n Prescindimos de la
∂ ∂
pequeña diferencia en dirección entre ambas normales: ≈ .
∂n ∂n '
El valor de la gravedad normal en P será el valor de la gravedad en Q más un incremento de la
función en la dirección QP que viene expresada por la derivada de la función respecto a la normal n,
multiplicado por la distancia entre Q y P (dada por N):
Página 43
∂γ ∂γ
γp = γQ + N ≈ γQ + N
∂n ' ∂n
Sustituyendo los valores de ∂g y γ p obtenemos:
∂T  ∂γ  ∂γ
∂g = g p − γ p ⇒ − = gp − γQ + N  = (g p − γ Q ) − N
∂n  ∂n  ∂n
∆g = g p − γ Q
∂T 1 ∂γ
Tp − = ∆g − T
N= ∂n γ ∂n
γQ
∂T 1 ∂γ
Podemos transformar la expresión de la siguiente manera: ∆g + − T=0
∂n γ ∂n
Esta expresión se llama ecuación fundamental de la Geodesia Física, porque relaciona la cantidad
medida ∆g con el potencial anómalo desconocido T.
Tiene la forma de una ecuación en derivadas parciales. Si ∆g fuera conocido en todo el espacio,
entonces esta ecuación podría ser discutida y resuelta como una auténtica ecuación en derivadas
parciales. No obstante, como ∆g es conocida sólo en la superficie (el geoide), esta ecuación
fundamental sólo puede usarse como una condición de contorno que por sí sola no es suficiente para
calcular T. Por consiguiente el nombre de la ecuación diferencial de la geodesia física que a veces
se da a esta ecuación es erróneo. No es una ecuación diferencial, es una condición de contorno.
Ordinariamente se supone que no hay masas fuera del geoide. Está claro que esto no es realmente
cierto. Las observaciones no las hacemos directamente sobre el geoide, sino sobre la superficie
física de la Tierra. El efecto de las masas fuera del geoide es eliminado al reducir las medidas de la
gravedad al geoide por cálculo, de manera que podemos suponer lógicamente que todas las masas
quedan encerradas por el geoide.
En este caso, puesto que la densidad ρ es cero en todas partes fuera del geoide, el potencial
anómalo T es allí armónico y satisface la ecuación de Laplace: T = 0.
∂ 2T ∂ 2T ∂ 2T
∆T = + + =0
∂x 2 ∂y 2 ∂z 2
Esta es claramente una verdadera ecuación en derivadas parciales que si se suplementa con la
condición de contorno, basta para determinar T en todo punto fuera del geoide.
∂T 1 ∂γ
Si escribimos la condición de contorno en la forma: − + T = ∆g
∂n γ ∂n
donde ∆g se supone conocida en todo punto de! geoide, vemos que sobre esta superficie se conoce
∂T
una combinación lineal de T y . Por tanto la determinación de T es entonces un tercer problema
∂n
de contorno de la teoría del potencial. Si se resuelve en T, la altitud del geoide, que es la cantidad
geométrica más importante en geodesia física, puede calcularse por la fórmula de Bruns.
Aparentemente es una ecuación diferencial en derivadas parciales. Decimos aparentemente porque
tomada como ecuación diferencial implica que las ∆g podrían conocerse no solo en el geoide sino
también fuera de él. Pero el valor de ∆g hasta el momento solo lo conocemos para puntos del
geoide. Luego no es una verdadera ecuación diferencial, es lo que llamamos condición de contorno
porque es una condición que cumple los valores de ∆g sobre el geoide.
De la resolución de una ecuación diferencial obtenemos una función del tipo y = f(x). A partir de
este momento conoceríamos, para cualquier valor de una de las variables, el de la otra en todo
punto. Sin embargo, en nuestro caso, al sustituir los valores discretos de las anomalías de gravedad
medidas sobre el geoide y resolviendo la ecuación diferencial obtenemos valores concretos que
Página 44
toma la función y = f(x), pero no obtenemos la expresión de la función en sí, por ello decimos que
no se trata de una verdadera ecuación diferencial, sino de una condición de contorno.
Sabemos que T = W-U, es la diferencia entre el potencial real y e! potencial normal, por tanto para
hallar los potenciales normales y para calcular las anomalías de gravedad se reducen al geoide las
gravedades observadas en superficie, eliminando así los efectos de cualquier masa exterior al geoide
que influya sobre el valor de la gravedad como pueden ser las interpuestas entre el punto en
superficie y éste. Es decir, estamos suponiendo todas las masas en el interior del geoide por lo que
la función del potencial perturbador cumple la Ecuación de Laplace, por tanto es una función
armónica, que puede desarrollarse mediante armónicos esféricos.
Teniendo en cuenta la condición de contorno, y sabiendo que T es armónica, podría lucirse la
expresión gracias a la cual conoceríamos el valor del potencial perturbador para todo punto sobre y
fuera del geoide.
Conociendo T en puntos sobre la superficie del geoide, y teniendo en cuenta la fórmula de Bruns,
podríamos calculada la gravedad normal en cada punto, conocer la ondulación en cada punto.
Se puede desarrollar en función de los armónicos esféricos. Podemos deducir la expresión del T
para puntos en el interior y en el exterior del geoide. Conocido este valor obtenemos la ondulación
del geoide por la fórmula de Bruns.
La fórmula de Stokes trata de dar a la función del potencial perturbador, considerada como una
función armónica, una expresión en función de anomalías de gravedad, para poder conocer los
valores de T.
4.5.- Aproximación esférica de la ecuación Fundamental de la Geodesia Física:

Si la tierra fuese una esfera de radio medio y la γ en un punto cualquiera, tanto en la superficie
m
como en el exterior, sería: γ = G .
r2
Veamos el error que cometemos al aceptar esta aproximación esférica. Al considerar la superficie
como esfera, se prescinde del aplanamiento geométrico α , que es del orden de 3*10-3 , luego este es
el error relativo que vamos a cometer en el cálculo de la ondulación. El error absoluto será del
orden 3*10-3 *N donde los valores de N, respecto a elipsoides generales (centrados en el centro de
masas de la. tierra: N en España positivas. El de Hayford es tangente en Postdam N = 0, pero no
está centrado en el centro de masas de la tierra: N en España es negativa; son pequeños, no superan
los 100 metros, por tanto el error máximo de N seria de 3*10-3*100 = 0,3m. Aun aceptando que el
error en la práctica fuese mayor, nunca supera el valor de 1 metro.
Tomaremos como radio medio terrestre el de aquella esfera que tenga el mismo volumen que el
elipsoide de referencia:
4 2 4
πa b = πR m ⇒ R m = 3 a 2 b ≈ 6371km
3 3
M
La gravedad normal media sobre esa esfera será: γ m = G ≈ 979,8gales o cm/sg
R 2m
∂ ∂ ∂T 1 ∂γ
En este caso = por lo que la ecuación será ahora: ∆g + − T=0
∂n ∂r ∂r γ ∂r
m ∂γ M 1 ∂γ 2
Como γ = G 2 y = −2G 3 ⇒ =−
r ∂r r γ ∂r r
Sustituyendo en la ecuación fundamental obtenemos la ecuación fundamental de la geodesia física
∂T 2
en aproximación esférica: ∆g + + T=0
∂r r
Sólo será una ecuación diferencial si conocemos ∆g en cualquier punto fuera de la superficie.
Página 45
∂T 2
Si r = Rm = R ⇒ ∆g + + T = 0 , tendremos la condición de contorno de la superficie de ∆g
∂r r
conocida.
4.6.- Anomalía de la gravedad fuera del geoide:

Al igual que la ondulación N era la distancia medida sobre la normal al elipsoide de nivel entre el
punto sobre el geoide Po y su proyección Qo sobre el elipsoide, ahora llamaremos anomalía en
altitud Wp a la distancia medida sobre la normal elipsódica entre el punto P situado en la superficie
geopotencial Wp y su proyección Q sobre la esferopotencial U Q = Wp .
Anomalía en el geoide: ∆g p = g p0 − γ Q0
Anomalía fuera del geoide: ∆g p = g p − γ Q
Cuasigeoide es la superficie obtenida al llevar sobre cada punto Q del elipsoide de nivel, la
anomalía en altitud Np en la dirección de la normal al elipsoide. Es una superficie muy parecida al
geoide. Coincidiría con éste en los puntos situados en el nivel medio de los mares en reposo.
Para el cálculo de las anomalías de gravedad fuera de la Tierra a partir de las anomalías de la
gravedad en superficie, existe la llamada fórmula:
R2  r 2 − R 2 1 3R 
∆g p = ∫∫σ  − − 2 cos ψ ∆g.dσ

4πr  l 3 r r 
ψ -distancia angular desde donde conozco la anomalía en el geoide hasta donde lo quiero conocer.
Si una función armónica H está dada en la superficie de la Tierra, como aproximación esférica, los
valores de H fuera de la Tierra pueden calcularse por la fórmula integral de Poisson:
R  r2 − R2 
Hp = ∫∫σ  H.dσ

4π  l 3 
El símbolo ∫∫σ es la abreviatura usual de una integral extendida sobre toda la esfera unidad, o
sobre el ángulo sólido completo. dσ designa el elemento de ángulo sólido definido como el
elemento da superficie de la esfera unidad. Por eso el elemento de superficie de la esfera terrestre r
= R es R2 dσ . El valor de la función armónica en el elemento variable de superficie R2 dσ se
designa simplemente por H mientras que Hp se refiere al punto fijo P. Obviamente, entonces:
l = r 2 + R 2 + 2Rr cos ψ .
La función armónica H puede desarrollarse en serie de armónicos esféricos:
2 n +1
R R R
H =   H 0 +   H 1 + ∑∞
n =2   Hn
r r r
2 n +1
R R R
H ' = H −  H 0 −   H1 = ∑∞
n =2   Hn
r r r
1 3
siendo: H 0 = ∫∫σ H.dσ y H1 = ∫∫σ H. cos ψ.dσ
4π 4π
Sustituyendo Ho y H1 por sus valores, y expresando H por la integral de Poisson:
R  r 2 − R 2 1 3R 
H' p = ∫∫σ  − − cos ψ H.dσ

4π  l 3 r r2 
Esta es la llamada Integral de Poisson modificada. Aplicando esta fórmula a r ∆g (función
armónica en todo el espacio), tenemos:
Página 46
R  r 2 − R 2 1 3R 
r.∆g p = ∫∫σ  − − cos ψ R.∆g.dσ

4π  l 3 r r2 
R2  r 2 − R 2 1 3R 
∆g p = ∫∫σ  − − cos ψ ∆g.dσ

4πr  l 3 r r2 
4.7.- Fórmula de Stokes

∂T 2 ∂T 2
La ecuación ∆g + + T=0 ⇒ ∆g = −
− T puede considerarse como una condición de
∂r r ∂r r
contorno solamente en tanto qué las anomalías de la gravedad ∆g son conocidas sólo en la
superficie de la tierra. Sin embargo, por la integral de prolongación ascendente, podemos calcular
las anomalías de la gravedad fuera de la tierra. La ecuación cambia de significado radicalmente
convirtiéndose en una auténtica ecuación diferencial que puede integrarse con respecto a r.
El potencial perturbador es una función de punto T = T(r, θ, λ) que particularizada para puntos de la
esfera de radio R se convertirá en T = T(R , θ, λ) . Con estas condiciones podemos tener T en puntos
de la superficie del geoide, y fuera de él:
R
T= ∫∫ ∆g.S(ψ).dσ
4π σ
1 ψ
− 6sen + 1 − 5 cos ψ − 3 cos ψ. ln sen ψ + sen 2 ψ 
siendo S(ψ) = ψ 2  2 2  la función de Stokes que
sen
2
ψ
se obtiene haciendo r = R y l = 2R.sen . Viene tabulada en función de la distancia angular ψ .
2
T
Por el teorema de Bruns tendremos la ondulación del geoide: N = y con ella obtenemos la
γm
R
fórmula integral de Stokes: N = ∫∫ ∆g.S(ψ).dσ
4πγ m σ
R 1
Como dσ = dS R 2 , obtenemos N = 2 ∫∫S
∆g.S(ψ ).dS = ∫∫ ∆g.S(ψ).dS
4πγ m R S
4πγ m R
4.8.- Evaluación práctica de la fórmula integral de Stokes

La fórmula integral de Stokes ha de ser evaluada aproximadamente mediante sumatorios. En lugar
de considerar infinitos términos infinitesimales de área dS, consideraremos un número finito k de
elementos de superficie qk que se obtienen por una adecuada subdivisión de la tierra, de manera que
la integral podrá ser evaluada como el sumatorio de todos estos términos.
Para la subdivisión del terreno alrededor del punto P, en el cual queremos calcular la ondulación del
geoide, construiremos sobre material transparente y a la escala del mapa gravimétrico a emplear,
una plantilla que centrada en el punto P muestre la subdivisión zonal del terreno.
Se emplean dos métodos diferentes de subdivisión:
- Cuadrícula: La subdivisión se realiza por las líneas coordenadas de algún sistema de coordenadas
fijado, en particular las coordenadas geográficas ϕ, λ . Estas forman bloques rectangulares (p. ej;
10'x 10'), llamados cuadrados, aunque normalmente no son cuadrados, tal como se definen en
geometría plana.
- Plantillas: La subdivisión se efectúa por círculos concéntricos de diferentes radios centrados en el
punto P, en el cual queremos calcular el valor de la ondulación. Las zonas comprendidas entre
radios contiguos se subdividen a su vez en compartimentos todos ellos de igual incremento angular.
Página 47
En cada compartimento las anomalías de la gravedad se reemplazan por un valor medio de dicho
compartimento, de modo que la ecuación integral será:
R
N= ∫∫ ∆g.S(ψ).dσ
4πγ m σ
R R
N= ∑k ∫∫q k ∆g k .S(ψ).dσ = ∑ ∆g k ∫∫q k .S(ψ).dσ = ∑k ∆g k c k
4πγ m 4πγ m k
los coeficientes ck pueden integrarse sobre cada compartimento dado que ahora la integral no
depende de las anomalías ∆g .
R
ck = ∫∫ S(ψ).dσ
4πγ m q k
Si el integrando, en este caso la función de Stokes S(ψ) se mantiene razonablemente constante sobre
el compartimento qk, puede sustituirse por su valor S(ψ k ) en el centro del compartimento:
R
ck = S(ψ k ) ∫∫q dσ
4πγ m k
1 1 A S(ψ k )
ck = S(ψ k ) ∫∫q R 2 .dσ = S(ψ k ) ∫∫q dσ = k
4πγ m R k 4πγ m R k 4πRγ m
A k - área del compartimento q k considerado.
Si los compartimentos están comprendidos entre líneas del tipo ψ = cte y α = cte , otra manera de
calcular el coeficiente c k será:
R
ck = S(ψ k ) ∫∫q dσ
4πγ m k
poniendo el diferencial de ángulo sólido en coordenadas polares: dσ = senψ.dψ.dα

R α2 ψ2
ck = ∫ ∫ S(ψ).senψ.dψ.dα
4πγ m α1 ψ1
R  α 2 −α 1  ψ2
ck =
4 πγ m
∫ψ 1
S  ψ . sen ψ . dψ
1 ψ ψ
J (ψ ) = ∫ S(ψ).senψ.dψ = ∫0 F(ψ).dψ
2 0
A esta función la llamamos de esta manera y como ha sido tabulada por Lambert y Darling
podemos obtener un valor para c k :
R (α 2 − α 1 )
ck = .2[ J(ψ 2 ) − J(ψ 1 )] es para cada punto
4πγ m
Por trigonometría esférica podemos calcular la distancia entre los puntos de anomalía conocida (
ϕ i' , λ'i ) y el punto P( ϕ, λ ):
cos ψ = senϕ.senϕ i' + cos ϕ cos ϕ i' cos(λ'i − λ)
Una vez calculado ψ podemos conocer la distancia sobre la esfera de radio medio:
r = Rψ rad
Página 48
4.9.- Componentes de la desviación de la vertical. Fórmulas de Vening- Meinesz
El ángulo formado en P por el geoide y una superficie paralela al elipsoide nos determina la
desviación de la vertical. Si tomamos el corte de dichas superficies por un plano que contenga a la
normal n' del elipsoide en P y a n y n’’ obtendremos la figura anterior. Si tomamos ahora el plano
normal en P que contenga a la dirección NS, tendremos que, la componente EW = η de ε es nula,
y en este caso el diferencial de longitud ds en la dirección NS es un diferencial de meridiano, por lo
que escribiremos ds = ds ϕ = dϕ . Al tomar un plano paralelo al ecuador por P tendremos que la
componente NS = ξ es nula, y ds será un diferencial de paralelo por lo que escribiremos
ds = ds λ = dλ
Al tener un ángulo muy pequeño podemos suponer que =tg =sen 
dN
Tendremos que: tgε = − ⇒ dN =−⋅ds
ds
En la dirección NS la componente EW = η es nula y por tanto:
dN =−⋅ds 
Como sabemos además que ds  =R⋅d 
dN 1 ∂N
− =− =
ds ϕ R ∂ϕ
La componente EW es cero por tanto  es el diferencial de meridianos.
En la dirección EW la componente NS = ξ es nula y por tanto:
dN =−⋅ds
Como sabemos además que ds  =R⋅cos ⋅d 
dN 1 ∂N
− =− =
ds λ R. cos ϕ ∂λ
La componente NS es cero por tanto  es el diferencial de paralelos.
Definición de los componentes de la desviación de la vertical a partir de la ondulación del geoide:
R
N= ∫∫ ∆g.S(ψ).dσ
4πγ m σ
Calculemos las derivadas de N respecto de ϕ y λ , La única componente de la expresión de la
fórmula de la ondulación dependiente de estas variables es la función de Stokes S( ψ ). Son
derivadas parciales ya que N es función de ambas variables, y en cada caso sólo derivamos en una
de ellas.
Página 49
∂N R ∂S(ψ ) ∂S(ψ ) ∂S ∂ψ ∂N R ∂S ∂ψ
= ∫∫σ ∆g. .dσ pero = por tanto = ∫∫σ ∆g. .dσ
∂ϕ 4πγ m ∂ϕ ∂ϕ ∂ψ ∂ϕ ∂ϕ 4πγ m ∂ψ ∂ϕ
∂N R ∂S(ψ ) ∂S(ψ ) ∂S ∂ψ ∂N R ∂S ∂ψ
= ∫∫σ ∆g. .dσ pero = por tanto = ∫∫σ ∆g. .dσ
∂λ 4πγ m ∂λ ∂λ ∂ψ ∂λ ∂λ 4πγ m ∂ψ ∂λ
∂ψ ∂ψ
Demostramos ahora las igualdades = − cos α y = − cos ϕ.senα
∂ϕ ∂λ
Por trigonometría esférica fijándonos en el triángulo de la figura:

Demostración 1:
cos a= cos b cos c + sen b sen c cos A
cos=cos90−⋅cos90− 'sen90−⋅sen90− '⋅cos'−
cos =sen⋅sen'cos⋅cos'⋅cos '−
Derivamos ahora esta expresión respecto a  :

−sen =cos ⋅sen '−sen⋅cos  '⋅cos'− [1]

sen a cos B= cos b sen c- sen b cos c cos A
sen⋅cos=sen'⋅cos−cos'⋅sen⋅cos '− [2]
Dividiendo [1] entre [2]:

=−cos 

Demostración 2:
cos=cos90−⋅cos90− 'sen90−⋅sen90− '⋅cos'−
cos =sen⋅sen'cos⋅cos'⋅cos '−
Derivamos ahora esta expresión respecto a  :

−sen =cos ⋅cos '⋅sen'− [3]

Por el teorema del seno:
sen⋅sen=cos '⋅sen '− ⇒
cos⋅sen⋅sen=cos⋅cos '⋅sen '− [4]
Dividiendo [3] entre [4] tenemos:
Página 50

=−cos  sen 

∂ψ ∂ψ
Sustituyendo los valores de = − cos α y − = cos ϕ.senα tenemos las expresiones finales de las
∂ϕ ∂λ
fórmulas de Vening-Meinesz:
1 dS
NS: ξ = ∫∫σ ∆g cos α.dσ
4πγ m dψ
1 dS
EW: η = ∫∫σ ∆g senα.dσ
4πγ m dψ
dS(ψ )
Por tanto, para poder calcularlas necesitamos calcular primeramente y el acimut α .
dψ
Derivando la función de Stokes con respecto a ψ obtenemos la función de Vening-Meinesz:
ψ ψ
cos  1 − sen 
dS ψ
 2  + 8senψ − 6 cos  − 3  2  + 3senψ. ln sen ψ  + sen 2  ψ 
=−    2  
dψ 2 ψ   2 senψ     2 
2sen  
2
El acimut α se deduce dividiendo:
senψ. cos α = cos ϕ.senϕ'−senϕ. cos ϕ' cos(λ '−λ) y senψ.senα = sen (λ'−λ). cos ϕ'
o
senψ.senα = sen (λ'−λ ).sen (90 − ϕ' ) y
senψ. cos α = cos(90 − ϕ' ).sen (90 − ϕ) − cos(90 − ϕ).sen (90 − ϕ' ) cos(λ '−λ )
sen (λ '−λ ) cos ϕ'
obteniendo tgα =
φ cos ϕsenϕ'−senϕ cos ϕ' cos(λ'−λ )
Página 51
5.- REDUCCIÓN DE LOS VALORES OBSERVADOS DE LA GRAVEDAD
5.1.- Objeto de las reducciones. Tipos de reducciones

La gravedad es un valor observado sobre la superficie terrestre. Antiguamente, se deducía gracias a
un péndulo, y actualmente, se utilizan gravímetros que nos dan incrementos de gravedad. Si
tenemos tres puntos situados a la misma latitud, según la figura, y comparamos los valores de la
gravedad en los mismos, veremos que la gravedad en P1 es menor que la gravedad en P2, y ésta a
su vez, mayor que la gravedad en P3, dado que la densidad del agua situada bajo P3 es menor que la
de la tierra situada bajo P2.
Vemos que la gravedad variará dependiendo de la altitud del punto, de la profundidad del océano,
etc.; es decir, los valores directamente observados no son comparables con la gravedad normal. Para
poder utilizar los valores de la gravedad en las fórmulas, hemos de reducir estos valores observados
al geoide. De este modo, podemos calcular las anomalías sobre el geoide, en base a las cuales se
deducen gran cantidad de parámetros, también útiles para deducir la distribución de densidades o
irregularidades de masa de la corteza terrestre causantes de "anomalías" en los valores de la
gravedad.
La gravedad es medida y reducida en un número discreto de puntos de la superficie, en base a los
cuales se trazan mapas de curvas isoanómalas, a partir de los cuales extrapolaremos el valor de la
anomalía y deduciremos el de la gravedad reducida en cualquier punto, y a partir de ésta
conoceremos la ondulación en dicho punto...
Para reducir estos valores de gravedad observados al geoide, hemos de aplicarles una serie de
correcciones o reducciones.
Tipos de reducciones:
 Reducciones no isostáticas:
o Reducción de aire libre o Faye
o Reducción de Bouguer:
 Corrección por lámina o placa de Bouguer
 Corrección topográfica
 Reducciones isostáticas:
o Reducción de Pratt-Hayford
o Reducción de Airy-Heiskaner
o Reducción de Venning-Meinesz
5.2.- Reducción aire libre o Faye

En la reducción aire libre o Faye se tiene en cuenta la altitud de la estación en que se observa el
valor de la gravedad.
g = gravedad observada en P
go = gravedad en Po , reducida al geoide
Página 52
∂g
g0 = g − H
∂h
∂g
Vemos que la corrección aplicada a la gravedad observada es: F = − H = δF
∂h
∂g ∂γ miligal
En la práctica se acepta: ≈ = −0.3086
∂h ∂h m
Por tanto: δ F = F = 0.3086H mgal con H en metros
g 0 = g + F = g + δ F = g + 0.3086H
Un método de obtenerla es:
m m
g≈G 2 si g 0 = G siendo R el radio medio terrestre
( R + H) R2
−2
m m 1 m 1 m  H
g≈G =G =G =G 1 + 
R 2 + H 2 + 2RH R2  2  R2  H
2 2
R  R
1 + 2 H + H  1 + 
 R R 2   R

−2
 H
Donde H R << 1 y 1 +  es de la forma (1 + x ) n cuyo desarrollo en serie de Taylor es nuestro caso
 R
que nos quedamos sólo con el término correspondiente a la aproximación lineal o de primer orden
 H
1 − 2  :
 R
m  H  H  H
g ≈ G 2 1 − 2  = g 0 1 − 2  =  g 0 − 2g 0 
R  R  R  R
g g 0 g
g 0 = g + 2 0 H donde 2 H = F = δ F ≈ 2 m H siendo estos los valores medios de la gravedad y
R R Rm
del radio terrestre respectivamente.
g 0 = g + δ F = g + 0.3086H siendo F positivo
Sólo para casos excepcionales de altitudes superiores a 2000 m se consideran términos de siguiente
orden en el desarrollo en serie de Taylor.
Se llama reducción aire libre porque sólo se tiene en cuenta la altitud de la estación y no la masa
interpuesta.
Para los casos en que H ≥ 2000 m se utiliza el segundo orden de aproximación:
m  H H 2 
g≈G 1 − 2 + 3
R 2  R R 2 
 H H2  g0  3 H
g ≈  g 0 − 2g 0 + 3g 0 2  ⇒ g0 = g + 2 H 1 − 
 R R  R  2 R

Salvo en casos excepcionales, la fórmula práctica utilizada es:
g 0 = g obs + 0.3086H = g F
A esta gravedad reducida únicamente utilizando la corrección aire libre se la llama también
gravedad Faye gF. La anomalía de gravedad Faye será:
∆g F = g F − γ = g obs + 0.3086H − γ .
5.3.- Reducción Bouguer.

Tiene en cuenta las masas interpuestas entre el nivel de la estación y el geoide. La descomponemos
en tres partes:
Página 53
5.3.1.- Reducción por lámina o placa Bouguer:

En este caso se acepta que la estación donde medimos la gravedad descansa sobre una "placa
infinita" o meseta plana e indefinida, de densidad constante e igual a la de la corteza (
ρ 0 = ρ media.de.la.corteza = 2.67 gr ) y de un espesor constante e igual a la altitud ortométrica H de la
cm 3
estación.
Al aplicar es la reducción se elimina el efecto que ha tenido, sobre el valor de la gravedad
observado en la estación, la atracción ejercida por dicha masa laminar indefinida: es decir, se tiene
en cuenta el efecto de las masas interpuestas entre el nivel de la estación y el Geoide.
Gravedad de Bouguer Simple g BS : Llamamos así al valor de gravedad obtenido al aplicar
únicamente la corrección por placa al de la gravedad Faye: g BS = g obs + δ F + δ B
Adelantamos que por el teorema de Gauss la intensidad de campo o atracción sobre la unidad de
masa ejercida por una placa indefinida sobre un punto situado en su superficie es − 2πGρ 0 H . Por lo
tanto la corrección Bouguer Simple será: δ BS = δ B ≡ B = −2πGρ 0 H
Aplicando valores: δ BS =δ B ≡B≈−0. 1119 H miligales con H en metros
En aproximación lineal: g BS = g obs + 0.3086H − 0.1119H = g obs + 0.1367H
La anomalía de Bouguer simple será: ∆g BS = g BS − γ
5.3.2.- Corrección topográfica:

El fin de esta corrección es llevar la estación a situación de meseta.
En esta corrección tenemos en cuenta la situación real de la estación, se tienen en cuenta los
accidentes topográficos entorno a la estación de medida. Tras eliminar los efectos de los excesos y
defectos de masa que se encuentran entorno a la estación, nos encontraremos con que el valor
obtenido para g es el que le correspondería a un punto situado sobre una meseta plana e indefinida,
o sea, las condiciones necesarias para que la aplicación de la corrección por placa sea correcta.
Recordemos que, al aplicar la corrección Bouguer Simple, estábamos suponiendo que la estación se
encontraba sobre la superficie de una placa indefinida y de densidad constante; y tras hacer esta
suposición, restábamos al valor observado de la gravedad, los efectos de dicha placa, obteniendo así
la g BS< g obs . Pero esta corrección δ B aplicada a la g obs es excesiva si existe una zona vacía ( ∆m < 0
, un defecto de masa) alrededor de la estación de medida, ya que al aplicar la corrección por placa
dábamos por supuesto que dicha zona estaba rellena de materia, y para compensarlo habremos de
sumarle el efecto que produciría sobre el valor de la gravedad una masa del mismo volumen y en la
misma posición que la zona vacía. Para ello habremos de sumar la componente vertical que dicha
fuerza atractiva provocarla en la estación P.
Si por el contrario existe una zona sobre el nivel de P, con ∆m > 0 , hemos de tener en cuenta que Fz
(la componente vertical de la fuerza de atracción F que dicha masa ejerce en P) es contraria a la
dirección de la gravedad, y actuó disminuyendo el valor de ésta en el momento en que realizamos
su medida, por lo tanto hemos de contrarrestar ésta pérdida de gravedad sumando al valor de g obs el
valor del módulo de dicha componente vertical Fz.
Vemos que en cualquier caso la corrección topográfica es siempre positiva.
Para el cálculo de la corrección topográfica hemos de elegir un modelo de subdivisión del terreno
en torno a la estación. Generalmente se utiliza el modelo de zonas circulares y compartimentos
según Hayford. Cada zona (corona cilíndrica) vendrá definida por sus radios interior y exterior,
siendo el radio interior de una zona el exterior de la zona interna contigua. Cada zona se divide en
un número n de compartimentos o bloques iguales dentro de dicha zona . El ángulo
Página 54
π
correspondiente a cada bloque vendrá determinado por α = 2 . La altura H1 (HA,HB,...) es la
n
altura media correspondiente al bloque considerado, deducida en función de las curvas de nivel que
aparecen en el sector de la zona ocupado por dicho bloque.
En la práctica hemos construir a la escala del mapa topográfico (con curvas de nivel) a utilizar, una
plantilla sobre material transparente en la que estén representadas las zonas y los correspondientes
compartimentos.
En el modelo de Hayford vemos que existen:
• Zonas Literales (A,B,...) - nos dan los valores del radio exterior de cada zona hasta 166,735
Kms alrededor de la estación.
• Zonas Numéricas (18, 17, .... 1) - nos dan el valor del radio exterior en forma angular, de
modo que rext = R m .α siendo Rm el radio medio terrestre.
Cuando se hace la corrección topográfica exclusivamente para calcular correcciones Bouguer,
normalmente sólo se tienen en cuenta las zonas literales. Si además van a aplicarse correcciones
isostáticas, habrá que calcular la corrección topográfica teniendo en cuenta también las zonas
numéricas; es decir, se tiene en cuenta toda la superficie de la tierra ya que el valor angular
correspondiente a la última zona numérica (zona 1) es de 180º.
Si el terreno es muy accidentado, quizás convendría aumentar el número de compartimentos,
y podría disminuirse en caso de que el terreno sea llano.
La fórmula básica que nos permite el cálculo de dichas correcciones corresponde a la componente
vertical de la fuerza ejercida por un anillo cilíndrico o zona sobre un punto de su eje.
Fz=2⋅⋅G⋅  r 21h 2−  r 21h 21− r 22h2  r 22h 21
Hs- Altura de estación

HP- Altura del compartimento
h- distancia del punto P a la base inferior del cilindro
h1- distancia del punto P a la base superior del cilindro
Se sustituye aire por corteza por tanto =aire =2.67 g/cm 3=2670 Kg/m3
r1 - radio interior de la zona
r2 - radio exterior de la zona
Este tipo de corrección así formulada sólo es válida para un compartimento terrestre. En el caso de
tener un compartimento oceánico la correción será la suma de dos correcciones:
• Sustitución de aire por corteza. Es similar a la anterior.
Página 55
• Sustitución de agua por corteza. Tendremos =corteza−agua = 2670−1027 Kg/m3
Obtenemos un resultado positivo siempre como definitivo.

La corrección topográfica correspondiente a un bloque será el cociente entre el la corrección
anterior entre el número de compartimentos de la zona.
La corrección topográfica es la suma de las correspondientes a todos los compartimentos:
δ r = ∑in=1 ( K Z ) i
Estas fórmulas corresponden a un modelo plano. Cuando se aplican a gran distancia de P y cuando
utilizamos las zonas numéricas, hay que efectuar rectificaciones por curvatura.
La corrección topográfica podemos calcularla aplicando directamente estas fórmulas, o haciendo
uso de tablas basadas en las mismas. Utilizar las tablas es algo más preciso, ya que incorporan la
corrección por esfericidad. Son válidas también para el cálculo de la corrección de Bouguer. Las
tablas de Cassini-Doré-Ballarín dan la componente vertical correspondiente a la atracción de una
zona del modelo Hayford suponiendo que la densidad es uno y están calculadas para distintas
alturas del compartimento (h) y suponiendo una de las bases del prisma a nivel de estación. El signo
dependerá de la situación del bloque. En ellas aparecen tabulados:
Denominación de la Zona, Radio Exterior de la Zona, Número de Compartimentos, h que es el
argumento necesario para entrar en las tablas y f que es el número de Cassini, su unidad es la
centésima del miligal y representa la atracción de toda una zona suponiendo que la densidad es uno.
5.3.3.- Reducción Bouguer completa:

La corrección Bouguer completa comprenderá las anteriores:
g BC = g obs + δ F + δ B + δ1
Pero la corrección por placa δ B = −2πGρH = −0.1119H es siempre negativa por lo que:
Página 56
g BC = g obs + δ F − δ B + δ1
Vemos que − δ B + δ1 < 0 dado que la corrección topográfica δ1 es pequeña frente al efecto de la
lámina.
La anomalía de Bouguer completa será: ∆g BC = g obs + δ F − δ B + δ1 − γ = g BC − γ
La anomalía Bouguer Completa se suele comportar de forma que:
• En zona continental o de meseta es negativa, tanto más negativa cuanto más alto es el
terreno.
• En zonas oceánicas es positiva.
• En zonas litorales o costeras será aproximadamente cero.
Los organismos oficiales emiten mapas de isoanómalas basados en anomalías Faye, Bouguer
Simple y Bouguer Completa. En la leyenda se incluyen:
• La densidad con que están calculadas dichas anomalías.
• La fórmula de la gravedad normal utilizada. Según el S.G. de Referencia adoptado.
• El nombre de la última zona del modelo de Hayford considerada para los cálculos.
5.4.- Consideraciones generales sobre la isostasia

El término de isostasia = equilibrio fue acuñado a finales del siglo XIX por Dutton. El concepto de
isostasia supone que todo exceso o defecto de masa por encima o debajo del nivel del geoide está
compensado de manera que a una cierta profundidad las masas se encuentran en un equilibrio
análogo al hidrostático, llamado equilibrio isostático.
A cierta profundidad existe una superficie de compensación. Esa compensación se puede formular
en términos de igualdad de masa o en términos de igualdad de presión.
Los hechos que avalaron la idea del equilibrio isostático fueron dos:
• Los análisis de las anomalías de Bouguer:
Al correlacionar mediante un esquema sencillo y simplificado, las anomalías de Bouguer con las
altitudes vemos que son nulas en la zona al nivel del mar, más negativas con la elevación del
terreno y más positivas cuanto más profunda es la zona oceánica. Una anomalía de Bouguer
negativa indica una deficiencia de masa por debajo del nivel del geoide. Se corresponde con las
llamadas raíces de las montañas en el modelo Airy-Heiskanen. Una anomalía positiva indica un
exceso de masa por debajo del nivel del geoide, también llamadas antiraíces de los océanos.
La sismología confirma que la corteza no es constante ni en densidad ni en espesor. Los perfiles
sísmicos ponen de manifiesto que el espesor de la corteza es mayor debajo de las zonas de montaña
que de las oceánicas.
• El cálculo de la desviación de la vertical:
Los valores obtenidos para la desviación de la vertical en puntos cercanos a las grandes montañas es
menor del que debiera obtenerse dado lo abrupto del relieve ya que una gran montaña produciría
sobre un péndulo una componente horizontal de atracción. Cuanto mayor sea la masa montañosa,
mayor debería ser el ángulo ϕ , pero no es así. Al determinar por métodos astro-geodésicos y
posteriormente por gravimétricos, la desviación de la vertical en zonas donde es de esperar una gran
desviación, el ángulo obtenido es mucho menor de lo que parece sugerir el relieve.
Bouguer tuvo ocasión de comprobarlo al medir el ángulo de la desviación junto a los Andes, en la
expedición al Perú de 1735.
Las determinaciones más importantes del ángulo de la desviación de la vertical fueron las realizadas
por Everest y Pratt (entre el 1820 y 1850) en zonas de la India próximas al Himalaya, donde el
relieve es muy accidentado, no sólo por las grandes masas montañosas, sino también porque al
borde de la placa continental en el océano Indico, existen grandes fosas oceánicas donde según sus
cálculos teóricos deberían obtenerse valores para ε̂ comprendidos entre 30 y 40 segundos (los
Página 57
obtenidos estaban comprendidos entre 5 y 8 segundos), poniendo así de manifiesto que debía existir
en el interior de la tierra algo que compensaba los efectos de aquellas irregularidades de masa.
5.5.- El sistema isostático Pratt-Hayford

La idea inicial es de Pratt, y el cálculo y uso de Hayford. Combina ideas de Pratt con ideas de
Hayford.
Suponemos que existe una "superficie de compensación" o "nivel de compensación", por debajo del
cual las masas se encuentran en equilibrio. Dicha superficie se halla a una profundidad del orden de
100 km por debajo del geoide.
En la capa más externa considerada desde la superficie del terreno hasta la superficie de
compensación, hablaremos de una compensación local o por columnas verticales.
Comparando columnas de igual sección correspondientes a zonas del terreno de diversas altitudes
sobre el geoide tendremos diferentes casos:
gr
- La columna de altitud cero tendría una profundidad D = 100 km, y una densidad ρ 0 = 2.67 .
cm 3
Considerando una sección con área unidad, tendrá una masa de
M = ρ 0 .V = ρ 0 .D
- La masa de una columna de altura h sobre el geoide será: M = ρ(D + h )
- La masa de una columna a una profundidad h' bajo el geoide será: M = ρ(D − h ' ) + ρ W h '
Suponemos que la masa de estas columnas de igual sección, por encima de la superficie de
compensación es la misma, por lo tanto sus densidades deben ser distintas ya que sus dimensiones
lo son. De modo que igualando las masas:
• Para una columna terrestre:
D
ρ 0 .D = ρ.(D + h ) ⇒ ρ= .ρ 0 , como h > 0  ρ < ρ 0
D+h
 D  h h
∆ρ = ρ − ρ 0 =  − 1.ρ 0 = ρ 0 ⇒ ∆ρ = ρ0
D+h  D+h D+h
• Para una columna que subyace bajo el agua:
D h'
D.ρ 0 = (D − h ' )ρ + h ' ρ W ⇒ ρ0 − ρ =ρ
D−h ' D−h ' w
Página 58
 D  h' h' h'

∆ρ = ρ − ρ 0 =  − 1.ρ 0 − ρW = (ρ 0 − ρ W ) ⇒ ∆ρ = (ρ 0 − ρ W )
 D − h'  D − h' D − h' D − h'
gr h'
La densidad media del agua del mar ρ W = 1.027 3 ⇒ ∆ρ = 1.643 siendo h' la
cm D − h'
profundidad a que se encuentra la sección superior de la columna.
Vemos, por tanto, que estamos considerando que la corteza no es de densidad constante, no es
homogénea. Nuestra compensación isostática será debida a esta diferencia de densidades.
Por debajo de las zonas de montaña la densidad de la columna es menor que 2,67, por ello las
anomalías de Bouguer son negativas y hablamos de que existe una deficiencia de masa. Por el
contrario debajo delos océanos la densidad es mayor, de ahí que las anomalías de Bouguer sean
positivas y hablamos de exceso de masa. Es decir, los excesos de masa quedan compensados por
defectos de masa internos, y bajo los océanos existe un exceso.
Si la tierra estuviese compensada isostáticamente al 100%, la masa que sobresale sería totalmente la
que falta en el interior. No es así, aunque sí lo está en un porcentaje muy elevado.
Esta corteza no se corresponde con la del modelo corteza-manto-núcleo, sino que se aproxima a lo
que denominamos litosfera (la profundidad de las placas litosféricas varia de 80 a 120 kms), que
comprende la corteza y parte superior del manto, y descansa sobre la astenosfera (zona en que
existe el equilibrio).
Vemos que las fórmulas de este sistema están deducidas en términos de igualdad de masa. Podría
haberse formulado en términos de igualdad de presión.
5.6.- El sistema isostático Airy-Heiskanen

Actualmente es el más utilizado, y el que debemos utilizar si queremos realizar
compensaciones isostáticas.
Hipótesis de Airy (1855): suponía que la capa más externa de la tierra (corteza) está flotando sobre
una capa inferior de densidad mayor, sobre una especie de magma denso.
Al igual que en el caso anterior, estamos considerando un modelo plano, y para el cálculo de los
valores de la raíz y antirraíz, realizamos también una compensación local o por columnas verticales
de igual sección. Igualmente se formulará en términos de igualdad de masas.
Página 59
Suponemos que la densidad de la corteza se mantiene constante y que la superficie de

compensación se encuentra situada aproximadamente a T = 30 km de profundidad bajo el geoide.
gr
La densidad correspondiente a dicha superficie será: ρ1 = 3.27
cm 3
A una columna que sobresale h por encima del geoide le corresponderá una raíz de profundidad t
por debajo del nivel de compensación. Para calcular el valor de t comparamos el empuje ejercido
por el exceso de masa de altura h, con el de compensación en sentido contrario, correspondiente al
defecto de masa (le llamamos defecto porque el espacio ocupado por la raíz de densidad ρ 0 < ρ1
debería estar ocupado por el magma de densidad ρ1 , por tanto falta una cantidad de masa igual a
t (ρ1 − ρ 0 ) ) de profundidad t bajo el magma denso.
Es semejante al equilibrio mantenido por un iceberg, pero sin tener en cuenta la zona T
comprendida entre el geoide y el nivel de compensación.
El empuje ejercido por la corteza sobre el manto será el ejercido por la parte que sobresale h más el
ejercido por la raíz t. No tenemos en cuenta la parte de la columna correspondiente a T dado que si
h no existe, la columna no tendría una raíz t, y en tal paso sólo existe la parte T, que está en
equilibrio.
Empuje de h y t: (h.ρ 0 + t.ρ 0 )g
(h.ρ 0 + t.ρ 0 )g = t.ρ1 .g
Empuje compensador: t.ρ1 .g
ρ0 2.67
h.ρ 0 = t.(ρ1 − ρ 0 ) ⇒ t= h= h = 4.45h
ρ1 − ρ 0 3.27 − 2.67
Decimos que está formulado en igualdad de masas, ya que está fórmula equivale a igualar el exceso
de masa h.ρ 0 con el defecto de masa t.(ρ1 − ρ 0 ) . Vemos que t > h .
Por lo tanto, cuanto mayor sea h, mayor será la profundidad t de la raíz.
El espesor de la corteza en este caso será: T + h + t
Debajo de una columna oceánica existirá una antirraíz. Donde debiera haber material de densidad
ρ 0 existe material de ρ 0 < ρ1 , por eso hablamos de la existencia de un exceso de masa.
La deficiencia de masa, debida a la diferencia entre las densidades de la corteza y del agua ρ 0 − ρ W
(donde debería haber corteza hay agua), queda compensada en el interior por un exceso de masa,
debido ahora a la diferencia entre las densidades del magma denso y de la corteza ρ1 − ρ 0 (donde
debería haber corteza hay magma).
ρ0 − ρ W 2.67 − 1.027
h ' (ρ 0 − ρ W ) = t ' (ρ1 − ρ 0 ) ⇒ t' = h' = h ' = 2.74h '
ρ1 − ρ 0 0.6
Vemos que t ' > h ' .
Por lo tanto, el espesor de la corteza por debajo del agua será: T – h’ – t’
Vemos que en este sistema varia la profundidad de la corteza, manteniéndose constante su densidad.
5.7.- Compensación regional de Vening-Meinesz o Regional

Es una modificación del modelo de Airy - Heiskanen, ya que hace una compensación regional en
vez de local. Pero en la práctica, este modelo apenas ha sido utilizado para el cálculo de
compensaciones isostáticas.
Suponemos que la corteza descansa sobre una superficie inferior cuyo comportamiento es de tipo
elástico. La superficie de discontinuidad que separa la corteza y el manto (Discontinuidad de
Mohorovicic) se comporta como una membrana elástica, con lo cual la deformación producida para
compensar un exceso de masa h no queda localizada en forma prismática justo debajo de dicho
Página 60
exceso, como se supone en el modelo de Airy-Heiskanen, sino repartida sobre dicha superficie.
5.8.- Reducciones isostáticas. Anomalías isostáticas

El objeto de la reducción isostática de los valores observados de la gravedad es la regularización de
la corteza1 basándonos en un determinado modelo de isostasia. Puede decirse que estamos
transportando los excesos de masa h al interior.
En el modelo de Pratt-Hayford las densidades de las columnas variaban de manera que la masa de
todas y cada una de ellas por encima de la superficie de compensación era idéntica. La columna de
altura D tenía la misma masa que la D + h. Si comprimiésemos la columna de altura D + h (sin
variar su sección) hasta reducirla a la medida D, o lo que es lo mismo, si la masa correspondiente a
la zona h (lo que llamamos exceso de masa) de dicha columna la introdujésemos en la zona D, su
densidad aumentarla de ρ a ρ 0 , por ello decimos que estamos transportando los excesos de masa al
interior, y de esta manera, se regulariza la corteza, se regulariza su densidad al valor de ρ 0 . En
definitiva, ha sido lo mismo que rellenar la corteza con la masa que sobresale.
En el modelo de Airy-Heiskanen, lo que varía no son las densidades de los bloque, sino sus
dimensiones, por lo que la profundidad o espesor de la corteza es diferente bajo cada bloque.
Para regularizar la corteza hemos de trasladar los excesos de masa de cada columna hasta su raíz de
manera que la densidad de ésta aumentará de ρ 0 hasta el valor de ρ1 , ya que la masa del exceso de
h es precisamente la que le falta a la raíz t para tener la misma densidad que el manto. De este
modo, quedará regularizado el espesor o profundidad de la corteza a un valor T = cte. La masa de la
antirraíz t’ la llevaremos a la zona oceánica h’ donde hay defecto de masa.
Llamamos δ1 a la corrección isostática o efecto de compensación isostática total.
La gravedad reducida con correcciones isostática será: g1 = g obs + δ F − δ BS + δ t + δ1
La anomalía en este caso será: ∆g1 = g red − γ = g 1 − γ
Vimos que la corrección Bouguer − δ BS + δ t = − δ T crecía negativamente con la altura. Este efecto
queda amortiguado (no extinguido) al incluir las correcciones isostáticas, debido a que si el modelo
utilizado para el cálculo de las reducciones isostáticas fuese perfecto y el método utilizado el
correcto, las correcciones isostáticas compensan la Bouguer δ I = − δ T y en este caso:
∆g1 = g1 − γ = g obs + δ F − δ BS + δ t + δ1 = g obs + δ F = ∆g F
Es decir, el valor delas anomalías tras aplicar todas las reducciones se asemeja al de las anomalías
aire-libre.
Efectivamente, al introducir correcciones isostáticas vemos que se suavizan los valores de las
anomalías de gravedad. Se comprueba que los valores de la anomalía ∆g1 (anomalía con correciones
isostáticas incluidas) son muy pequeños, y están distribuidos de una manera aleatoria.
Cómo realizar las reducciones isostáticas en la práctica:
A la hora de calcular la reducción isostática, también llamado efecto de compensación isostática,
hemos de elegir un modelo de subdivisión de la Tierra en torno a la estación del tipo zonas
circulares y compartimentos (p. ej. el modelo de Hayford) que se extienda a toda la superficie
terrestre. Hay que considerar tanto las zonas literales como las zonas numéricas.
Este efecto de compensación será igual al sumatorio de los efectos cada uno de los compartimentos
de las diferentes zonas: δ1 = ∑ (k 1 ) q
Tanto para el cálculo de reducciones isostáticas en el sistema Pratt como en el Airy se utiliza el
modelo de división de Hayford.
Nunca se realizan correcciones isostáticas sin realizar las topográficas de las zonas numéricas. Por
1 Aquí el término corteza se refiere a la zona comprendida entre el geoide y la superficie de compensación.
Página 61
ello muchas veces se habla de reducción topo-isostática como suma de ambas.
5.8.1.- Sistema de Airy-Heiskanen:

Las fórmulas para calcular el efecto de la compensación isostática son:
• Zonas continentales
Fz=2⋅⋅G⋅   r 21h 2−  r 21h 21− r 22h 2  r 22h 21
h1=THs , h=THst ,  =1− 0=3.27−2.67 g/cm3=600 Kg/m3
• Zonas oceánicas
Fz=−2⋅⋅G⋅   r 21h 2− r 21h21−  r 22h 2 r 22h21 
h1=T−t'H s , h=THs
Con  =1−0 =3.27−2.67 g/cm 3=600 Kg/m 3 en los dos casos
Página 62
6.- MEDIDAS DE GRAVEDAD
6.1.- Las medidas de gravedad. Métodos para su determinación

Estas medidas son necesarias para conocer como está distribuida la gravedad sobre la superficie
terrestre. Necesario para conocer el campo gravitatorio terrestre. Llamamos unidad gravimétrica
u.g. a la décima del miligal. Es una unidad de uso fecuente en prospección gravimétrica.
6.2.- Clasificación según los tipos de medida:

• Absolutas: como resultado de una medida absoluta se obtiene el valor de la gravedad en el
punto de observación. Por tanto el instrumento a utilizar ha de darnos el valor de la gravedad
en un punto, con independencia del valor de la gravedad en cualquier otro punto. Este tipo
de mediciones son las más laboriosas y delicadas. Se realizan sólo en determinados "puntos
base".
• Relativas: como resultado de una medición relativa, se deduce la diferencia de gravedad
entre dos puntos. No obtenemos la gravedad en el punto de observación. Si con el mismo
instrumento obtenemos lecturas en varios puntos, deduciremos diferencias de gravedad entre
ellos. Se utilizan para enlazar puntos base en los que se conoce la gravedad absoluta. De este
modo, podremos conocer la gravedad en todos los puntos.
6.3.- Clasificación según los métodos de medición:

• Dinámicos: se observa el movimiento de un cuerpo bajo la acción de gravedad
midiéndose el tiempo empleado por éste para pasar de una posición a otra. Entre los
diferentes instrumentos dinámicos varían el cuerpo y el movimiento a observar. Se
emplean en determinaciones absolutas de gravedad. Inicialmente fueron empleados en
determinaciones relativas.
• Estáticos: se observa el cambio en la posición de equilibrio de un cuerpo sometido a la
acción de la gravedad y a otra fuerza antagonista que equilibra a ésta, generalmente
ejercida por un elemento elástico. El cambio en la posición de equilibrio se traduce en un
desplazamiento (cantidad medible) lineal o angular del cuerpo. Se utiliza para medidas
relativas.
6.4.- Métodos absolutos y relativos:

7. Métodos absolutos (son todos métodos dinámicos):
• Péndulos: fue el primer instrumento que se utilizó para medir la gravedad. El que fue
fundamentalmente utilizado es el péndulo reversible.
• Péndulo matemático o simple
• Péndulo físico o compuesto
• Péndulo reversible
Página 63
•Basados en el movimiento libre de graves: cuerpo pequeño o partícula material que se

mueven bajo la acción de la gravedad.
• Caída libre
• Movimiento libre simétrico (elevación y caída)
Como resultado de una medida absoluta de la gravedad obtenemos un valor de la gravedad en el
punto. Actualmente los métodos pendulares están en desuso, se emplean los basados en el
movimiento libre de graves.
8. Métodos relativos:
• Péndulos de relativas: son métodos dinámicos. En desuso.
• Gravímetros: métodos estáticos. Todos los utilizados en la actualidad.
6.5.- Péndulo matemático o simple:
Consiste en un punto material de masa m suspendido mediante un hilo inextensible de longitud l

fijo por el extremo opuesto oscilando alrededor de un eje horizontal que pasa por el punto fijo O.
La masa del hilo es despreciable frente a la masa del punto material. Si se desplaza un ángulo α de
su posición de equilibrio y se deja libre, se inicia un movimiento oscilatorio en torno a la posición
de equilibrio.
Para saber cuál es la ecuación que rige el movimiento del péndulo, hemos de conocer la fuerza
peso, que se descompone en dos componentes, siendo la tangencial la que nos interesa. Esta
componente es tangente a la trayectoria del péndulo, el cual, describe una circunferencia de radio l.
Ft = m.g.senα es la fuerza que provoca el movimiento.
d 2α
Ft = m.a t siendo a t = l la aceleración tangencial = aceleración angular x radio.
dt 2
La ecuación del movimiento es:
d 2α d 2α g
m.g.senα = −m 2
l ⇒ 2
= − senα
dt dt l
Ecuación diferencial que integrada adecuadamente nos lleva a deducir el periodo de oscilación que
es el intervalo de tiempo transcurrido entre dos posiciones consecutivas idénticas (incluso en la
dirección de su movimiento) o tiempo empleado en una oscilación completa:
l  1 
2 2
α  1.3  4 α
T = 2π 1 +   sen 2 +   sen + ...
g   2  2  2.4  2 
α
, siendo el valor máximo que aparezca en la oscilación, se corresponde con la amplitud.
Este es el periodo de oscilación para oscilaciones de cualquier amplitud. Para oscilaciones de
pequeña amplitud se puede confundir el seno de un ángulo con el valor del ángulo en radianes:
α α
sen ≈ .
2 2
Página 64
l 1 
La fórmula práctica para oscilaciones de pequeña amplitud es: T = 2π 1 + α2  .
g  16 
No obstante el período depende de la amplitud por lo que todos los métodos pendulares implicaban
unas correcciones por amplitud, dado que, por muy pequeña que quisiera hacerse esa amplitud,
siempre tenia un valor apreciable. Aplicando estas correcciones se deducía un valor T0:
l
lim α →0 T = T0 = 2π
g
que se acepta como fórmula correcta para oscilaciones de muy pequeña amplitud. En este caso, si
conseguimos medir el tiempo T0 con precisión (para lo cual medimos el que tarda en oscilar n veces
y promediamos) y conocemos la longitud del péndulo con precisión, se puede obtener valor de la
4π 2 l
gravedad: g = T02
El péndulo matemático es un ideal. En la realidad se producen deformaciones que afectan al hilo, e
incluso al soporte, por ello en la práctica se han utilizado más los péndulos físicos.
6.6.- Péndulo físico o compuesto:
Es un sólido rígido oscilando en torno a un eje horizontal, siendo esta oscilación debida solamente a
la acción de la gravedad. El cuerpo que se toma como sólido es un cuerpo con un cierto eje de
simetría, como lo puede ser una barra prismática.
La fuerza peso se manifestará actuando sobre el centro de gravedad. El punto de suspensión O es la
intersección del eje de oscilación con el plano perpendicular a él y que contiene a la fuerza peso.
La ecuación de su movimiento de oscilación vendrá dada por la igualdad entre el momento (fuerza
x distancia) respecto al eje de giro, de la fuerza peso, y el de las fuerzas exteriores:
d 2α d 2α m.g.h
m.g.h.senα = −I 2
⇒ 2
=− senα
dt dt I
2
siendo M e = −I d α el momento de las fuerzas exteriores respecto al eje de giro; igual al
dt 2
d 2α
momento de inercia respecto al eje de giro I por la aceleración angular .
dt 2
Comparando esta ecuación de movimiento con la del péndulo simple, vemos que existe gran
d 2α g
semejanza, ya que obtendremos: 2
=− senα
dt Ie
I
Sin más que igualar l e = , que dimensionalmente es una longitud dado que m.h = kg × m y
m.h
I = kg × m 2 . Le es la denominada longitud equivalente del péndulo físico. Representa la longitud que
debiera tener un péndulo matemático para que, colocado en el mismo lugar que el físico, oscilase
Página 65
con el mismo periodo que este.

Para obtener la gravedad y el período podemos emplear las mismas expresiones:
4π 2 l e l l
g= T0 = 2π = 2π e
T02 m.g.h g
Vemos que, la utilización de este tipo de péndulos requiere igualmente de la medición de tiempos y
de longitudes. Es necesario conocer bien la Le.
6.7.- Péndulo reversible

I
Si ahora transformamos la expresión l e = aplicando el teorema de Steiner que dice: "El
m.h
momento de inercia de un cuerpo respecto a un eje cualquiera es igual al momento de inercia
respecto a un eje paralelo a él más el producto de la masa del cuerpo por el cuadrado de la distancia
entre dichos ejes". Por lo tanto, poniendo el momento de inercia respecto a O en función del
correspondiente a un eje paralelo que pase por el centro de masas G del cuerpo:
I G + m.h 2 IG
I = I G + m.h 2 ⇒ le = ⇒ m.h 2 − l e .m.h + I G = 0 ⇒ h 2 − l e .h + =0
m.h m
De esta ecuación de segundo grado sacamos dos soluciones: h1 y h2. Por la propiedad de las raíces
b
de este tipo de ecuaciones: ah2 + bh + c = 0 siendo h1 y h2 soluciones, entonces h1 + h2= tenemos
a
que: h1 + h2= le (una longitud no puede ser negativa).
El péndulo físico tiene el mismo periodo de oscilación en torno a un eje horizontal que pasa por O1,
situado a una distancia h1 del centro de gravedad, que en tomo a un eje O2, situado a la distancia h2
del mismo y alineado con O1 y con dicho centro de gravedad G. Se dice entonces que el péndulo es
reversible respecto a dichos ejes. Para todo péndulo físico existirá otro eje respecto al cual será
reversible.
le
Como el periodo viene dado por T0 = 2π donde h1 + h2= le tanto si el cuerpo está suspendido
g
del eje O1 como si lo está del O2, concluimos que el cuerpo tendrá el mismo período de oscilación
tanto si lo hacemos oscilar respecto al eje horizontal que pasa por O1 como si lo hacemos oscilar
respecto al eje horizontal que pasa por O2.
I
Pero conocer le si tengo que emplear su expresión l e = requiere conocer I, m, la posición del
m.h
punto G, medir la h... ; es decir, es muy complicado a la hora de realizar las medidas absolutas de
gravedad, por ello se ideó el péndulo reversible, cuya ventaja es que nos permite conocer le
fácilmente.
Catre en el año 1818 realizó una medida de gravedad en Londres utilizando un péndulo que fue
denominado péndulo reversible de Catre. En este péndulo los ejes O1 y O2 son fijos, por lo que
también lo será le, pero existen dos masas que pueden fijarse al cuerpo principal en posiciones
distintas, de manera que variando las masas de posición ha de conseguirse que el periodo de
oscilación del sistema sea el mismo respecto a ambos ejes, de este modo sabremos que el sistema es
reversible y puesto que la distancia le es conocida, no será necesario medir la masa del cuerpo, ni
conocer la posición de su centro de gravedad, ni su momento de inercia,...
Página 66
le 4π 2 l e
T1 = T2 T = 0 = 2π donde g =
g T02
Se utilizaron mucho. Por ejemplo para dar la gravedad en Postdam.
La primera medida de gravedad absoluta que se realizó en Madrid en 1877, fue hecha por Barraquer
utilizando un péndulo reversible de Repsold, basado en el anterior, pero con las modificaciones
propuestas por Bessel en 1864 y Repsold.
Actualmente los métodos pendulares hace ya tiempo que se encuentran en desuso.
6.8.- Métodos basados en el movimiento libre de graves:

Son los empleados en la actualidad. Se basan en la medición de tiempos y distancias, actualmente
posible con las precisiones necesarias.
6.8.1.- Caída libre:
Se utiliza la ecuación del movimiento uniformemente acelerado: z = z 0 + v 0 t + 1 2 gt .

2
El origen de espacios no ha de coincidir necesariamente con el de tiempos. Cuando t0 = O

1
generalmente z 0 ≠ 0 . En esta ecuación t0 y z 0 son constantes, y también lo será g.
2
Se deja caer el cuerpo, y se miden los tiempos empleados (con relojes de cuarzo,...) y las distancias
recorridas (por métodos de medida interferométricos) por el grave respecto a un instante de
referencia. La precisión en la medida de g avanzará en la medida en que lo hagan los sistemas de
medición de tiempo y de distancias. En la práctica se miden n distancias y n tiempos, y se realiza un
ajuste de observaciones:
z i = z 0 + v 0 t i + 1 gt i2
2
El número mínimo de observaciones será de tres: z1 = z 0 + v 0 t 1 + 1 gt 12
2
z 2 = z 0 + v 0 t 2 + 1 gt 22
2
z3 = z0 + v0 t 3 + 1 gt 2
2 3
2( (z1 − z 2 )( t 1 − t 3 ) − (z1 − z 3 )(t 1 − t 3 ) )
Despejando: g =
( t 1 − t 2 )( t 1 − t 3 )( t 2 − t 3 )
Página 67
6.8.2.- - Movimiento libre simétrico (elevación y caída):
Observamos los instantes de paso del móvil por el nivel I y por el nivel II tanto en el ascenso como
en el descenso. Los niveles están en posiciones intermedias de la trayectoria. No deben elegirse ni
en el punto inicial ni en el punto más alto de la trayectoria.
Utilizando la ecuación del movimiento y considerándolo como un solo movimiento uniformemente
decelerado, y no como un decelerado hasta la posición más alta y uno acelerado desde ésta,
podremos escribir:
Ascenso: z1 = z 0 + v 0 t 1 − 1 2 gt 1 z 2 = z 0 + v 0 t 2 − 1 gt 22
2
2
Descenso: z 3 = z 0 + v 0 t 3 − 1 2 gt 3 = z 2 z 4 = z 0 + v 0 t 4 − 1 gt 24 = z1
2
2
8(z1 − z 2 ) 8H
g= 2 2 ⇒ g= siendo H = z − z
(t 4 − t1 ) − (t 3 − t 2 ) T22 − T12 2 1
T1 = t 3 − t 2
T2 = t 4 − t 1
Actualmente se fabrican estaciones portátiles para la medida absoluta de la gravedad que consisten
en cámaras de vacío en las que cae el cuerpo y utilizan métodos interferométricos para medir la
distancia. Son caras y grandes.
6.9.- Gravímetros
Un gravímetro en general consiste en un sistema elástico que se deforma bajo la acción de la
gravedad. Los métodos estáticos se utilizan en medidas relativas. El resultado da la medición es
una lectura expresada en divisiones de la escala que utilice el instrumento. La gravedad está
relacionada con la lectura G = f (L) donde f (L) es una función que debemos conocer. El calibrado
consiste en determinar esta función; es la llamada función de calibrado, que será facilitada por el
constructor del gravímetro.
Estacionando el gravímetro en dos puntos:
Punto 1: g1L1 ∆g 21 = g 2 − g1
Punto 2: 2 2g L ∆L 21 = L 2 − L1
Para pasar de diferencias de lecturas a diferencias de gravedad hemos de conocer la constante del
gravímetro:
∆g 21 = cte × ∆L 21 = C × ∆L 21 las unidades de C son miligales / división de la escala del gravímetro.
g 2 − g 1 = C( L 2 − L 1 ) (L 2 − L1 ) es el número de divisiones
Existen gravímetros con una C única para todo el rango de observación que se pueda realizar con él
Página 68
y hay en los que la constante varía según el rango de lectura.

Existen, en general, dos modelos de gravímetros estáticos:
• Geodésicos de amplio rango: deben de ser capaces de medir diferencias de gravedad entre
dos puntos cualesquiera de la superficie terrestre, aunque existan entre ellos grandes
diferencias de altitud y de latitud (amplio rango). Son los más caros. Ej: LaCoste-Romberg".
• Prospector: son modelos de pequeño rango. Se utilizan en prospecciones gravimétricas. Por
ejemplo: no pueden medir diferencias de más de 200 miligales. Por lo que, si la gravedad
entre dos puntos tiene una variación mayor que el rango de medición del instrumento, para
conocer su diferencia de gravedad habremos de realizar medidas intermedias reinicializando
la escala del instrumento. Son más baratos y sencillos.
6.9.1.- Gravímetro lineal
• Basados en resortes o muelles: son los más utilizados. Existen también de gas, fibras,...
• Lineales o balanza de resorte vertical: llamado Pesón por los franceses. En desuso. Basados
en el mismo principio que los sismómetros de componente vertical.
Está constituido por un resorte (hilo enrollado en hélice) con longitud inicial L0 del cual se
suspenderá un cuerpo de masa m. En su posición de equilibrio la longitud del resorte será L.
L0- longitud inicial del resorte sin la deformación debida a la masa m.
m -suponemos que la masa está concentrada en el centro de gravedad del cuerpo.
g - es la gravedad que queremos conocer. No es posible conocerla mediante una única observación.
Viendo la variación longitudinal del resorte al situar el instrumento en dos puntos distintos,
podremos conocer el incremento de la gravedad entre ellos.
Cuando el instrumento se encuentra en equilibrio (posición 1) actuarán las siguientes fuerzas:
• Peso = m.g
• Tensión correspondiente a la fuerza elástica: T = K(L-L0) donde K es la constante
recuperadora o elástica del resorte. Esta constante depende de su naturaleza, es decir del
material del cual está hecho, de su coeficiente de elasticidad,...; en definitiva, del
comportamiento elástico del mismo.
Ecuación del equilibrio:
K
P=T ⇒ m.g = K(L- L0) ⇒ g= (L − L 0 ) ⇒ g = f(L)
m
Si cambio el instrumento de lugar, cambiarán g y L.
Si desplazamos la masa una distancia z desde su posición de equilibrio (posición 2), el sistema
oscila y trabaja como un péndulo vertical.
d2z d2z
La ecuación del movimiento del sistema será: T = m.a = m = −kz ⇒ m + kz = 0
dt 2 dt 2
Resolviendo esta ecuación: z = z 0 . cos(ω0 t + ϕ 0 )
Página 69
2
ω0 =
2π K  2π 
es la frecuencia angular característica del resorte: ω02 = =  
T0 m  T0 
2π
T0 = = 2π m es el periodo característico del sistema. Depende de la constante recuperadora
ω K
0
del resorte y de la masa. El usuario no puede modificar este T0. Es algo fijo por construcción.
Sensibilidad mecánica del instrumento: depende de la construcción y las características elásticas del
instrumento; de la masa utilizada m y de la constante recuperadora del resorte K y es directamente
proporcional al cuadrado del periodo característico T0.
∆L m T02
m.∆g = K.∆L ⇒ = =
∆g K 4π 2
Se procura construir gravímetros con grandes T0 para poder apreciar pequeñas variaciones o
incrementos de la gravedad.
Pero la sensibilidad, en la práctica, esta limitada. Por ejemplo: sí queremos que nuestro instrumento
aprecie diferencias de gravedad de una décima de miligal y aproximando el valor de la gravedad g =
lOOO gales
∆g ∆L
, dividimos la expresión m.∆g = K.∆L por m.g = K(L - Lo) término a término: =
g L − L0
∆g
⇒ ∆L = ( L − L 0 )
g
10 −1 mgal
Si L − L 0 = 0.5m ⇒ ∆L = 0.5m 6
= 0.5 × 10 −7 metros = 0.05micras
10 mgal
Es decir, para poder medir variaciones de una décima de miligal, deberíamos poder apreciar en
nuestro instrumento un incremento de longitud ∆L de media décima de micra.
Podríamos construir un instrumento en el que L – L0 fuese muy grande, de manera que el ∆L no
resultase un valor tan pequeño, y por lo tanto tan difícil de apreciar, pero un instrumento de tales
proporciones no sería manejable.
Concluimos que con esta disposición o configuración del instrumento, la sensibilidad tiene un
límite. Debido a esto aparecen los llamados gravímetros no lineales.
6.9.2.- Gravímetros no lineales o astáticos. Astatización
Dentro de éstos están los llamados gravímetros astáticos, también llamados inestables en algunos
Página 70
libros.
En física, la palabra astático es sinónimo de equilibrio indiferente. Los llamados gravímetros
astáticos trabajan en unas condiciones tales que las zonas de equilibrio estable e inestable son muy
próximas y se dirá que las condiciones de astatismo perfecto se habrán alcanzado cuando el
equilibrio sea indiferente.
Esquema básico del gravímetro astático:
La masa está en el extremo de un brazo o balancín o palanca capaz de girar en torno a un eje
horizontal que pasa por su extremo. El sistema mecánico es capaz de girar en torno a un eje
horizontal que pase por O. Se mantiene en equilibrio bajo la acción de la fuerza peso y de la tensión
o fuerza elástica que se desarrolla en el resorte.
Es un sistema con un sólo grado de libertad. En este caso la variación es angular. Ahora la ecuación
de equilibrio vendrá dada en función de las fuerzas que se ejercen respecto al eje de giro. El
equilibrio viene dado por tanto como igualdad de momentos en torno al eje de giro.
Ecuación de equilibrio:
Momento de la fuerza peso: M g = OG ∧ P y su módulo M g = a.m.g.sen (α + δ)
Momento del resorte de módulo M R =T⋅
Deducción del valor de λ :
Tenemos dos expresiones para el área del triángulo OCD:
1
b.d.senα 1 1 b.d
2 b.d.senα = Lλ  senα = λ
1 2 2 L
Lλ
2
bd
Tendremos por tanto M R = K (L − L 0 )λ = K (L − L 0 ) senα
L
El momento resultante será: M = M g + M R
Estos momentos ejercen giros opuestos sobre el eje de giro y por lo tanto:
M = Mg − MR
La condición de equilibrio será M = 0. Tendremos: M g − M R = 0
bd
a.m.g.sen (α + δ) − K ( L − L 0 ) senα = 0
L
Página 71
dα
Ahora la sensibilidad se define como ya que α es ahora la magnitud que principalmente varía
dg
al variar la gravedad.
Estudios realizados sobre la sensibilidad pusieron de manifiesto que si L0 = 0 pueden conseguirse
mayores sensibilidades. Los resortes que cumplen esta condición son los llamados resortes de
longitud cero. Esto no quiere decir que en algún momento la longitud del resorte sea nula, quiere
decir que, dentro del rango de actuación del instrumento, el resorte se comporta de tal manera que la
tensión que ejerce es directamente proporcional a su longitud: T = KL, de manera que esta clase de
resortes se ajusta a una gráfica del tipo:
Vemos que el origen está contenido en la prolongación de la gráfica de la función, lo cual quiere
decir que si la tensión fuese nula, la longitud del resorte sería nula también; de ahí su denominación.
Los límites de actuación o rango del resorte están comprendidos entre L1...Ln.
La fabricación de resortes de longitud cero es complicada. Los resortes de este tipo fueron una
aportación muy importante para los gravímetros.
Los resortes generalmente se ajustan a una gráfica de este otro tipo, en la que podemos ver que
existe una longitud de que la tensión se anule:
Este tipo de gravímetros no lineales, está basado en el sismómetro de componente vertical de

LaCoste con la incorporación de este tipo de resortes. Son los utilizados en la actualidad. Ej:
gravímetro LaCoste-Romberg.
bd
Si L 0 = 0 ⇒ M R = Tλ = KL senα = K.b.d.senα ⇒ a.m.g.sen (α + δ) − Kbdsenα = 0
L
Para aumentar la sensibilidad consiguiendo así un sistema astático, el montaje debe hacerse de tal
modo que: α + δ = 90 siendo α mucho mayor que δ . El ángulo δ suele ser del orden de pocos
segundos. Y la máxima sensibilidad se conseguirá si δ = 0 en cuyo caso la sensibilidad será infinita,
ya que el sistema tendría un periodo propio de oscilación infinito. Con estas condiciones decimos
que la astatización es perfecta.
(a.m.g − Kbd)senα = 0
Se cumplirá para cualquier valor de α si: a.m.g − Kbd = 0 por tanto esta es la condición de
Página 72
a.m
equilibrio (LaCoste-Romberg) de donde: d = g
K.b
a.m
Vemos que d es proporcional a g ya que en todos los términos son constantes. Los cambios de
K.b
gravedad están relacionados con la variación necesaria en d para mantener el equilibrio del sistema.
De manera que, si g varía, hemos de modificar el valor de d para que el sistema recupere de nuevo
el equilibrio.
Si en el punto D situamos un tornillo micrométrico graduado, en vez de una unión fija, podremos
controlar la variación necesaria de la distancia d para mantener el brazo horizontal en cada punto,
de manera que podremos relacionar la variación de la gravedad existente entre dos puntos con las
lecturas realizadas sobre el tornillo micrométrico necesarias para mantener en ambos puntos el
instrumento en la situación de equilibrio.
Sistemas astáticos:
Teníamos que M g = M g (g, α) y M R = M R (α)
La condición de equilibrio del sistema consiste en que. la suma vectorial de ambos momentos sea
cero: M = M g + M R = 0
Ambos vectores actúan en la misma dirección y sentidos contrarios. Considerando sus módulos:
M = Mg − MR = 0 Mg = MR
Representando las funciones Mg y MR respecto a la variable α , existirá equilibrio en aquellos

puntos de corte tales como α 1 y α 2, en los que se cumple que M g = M R .
Para analizar el tipo de equilibrio existente entre estos puntos de corte, hemos de estudiar las
variaciones de Mg y MR, en definitiva, la variación del momento resultante M respecto de α :
∂M
Equilibrio inestable >0
∂α
∂M
Equilibrio estable <0
∂α
∂M
Equilibrio indiferente =0
∂α
Analicemos el tipo de equilibrio existente en el punto de equilibrio correspondiente a α 1, para ello
veamos el tipo de equilibrio existente en sus proximidades, en puntos tales como α1 + dα y α1 − dα
:
∂M
En α1 + dα con dα > 0 vemos que M g > M R ⇒ M>0 ⇒ >0
∂α
En α1 − dα con α1 + dα y dα < 0 vemos que M g < M R ⇒ M<0 ⇒ pero en este
Página 73
∂M
caso también nos da >0.
∂α
En ambos casos da positivo y por ello decimos que α 1 es un punto en el que el equilibrio es
inestable.
Analicemos ¡as proximidades del punto α 2 para lo que consideramos puntos tales como α 2 + dα y
α 2 − dα :
∂M
En α 2 + dα con dα > 0 vemos que M g < M R ⇒ M<0 ⇒ <0
∂α
En α 2 − dα con α 2 + dα y dα < 0 vemos que M g > M R ⇒ M>0 ⇒ pero en este
∂M
caso también nos da <0.
∂α
Vemos que en ambos casos nos da negativo, por ello, decimos que α 2 es un punto en el que el
equilibrio es estable.
Se entiende por astatización el conseguir incrementar la sensibilidad mecánica de un sistema
rotacional actuando sobre los momentos de la fuerza peso y la fuerza antagonista.
Para llegar a la astatización perfecta, la sensibilidad del gravímetro debe ser infinita, hay que
∂M ∂M g ∂M R ∂M g ∂M R dα
conseguir que: =0 ⇒ − =0 ⇒ ≈ ⇒ ≈∞
∂α ∂α ∂α ∂α ∂α dg
∂M g
es el llamado momento característico. Se estudia qué tipo de funciones han de ser las MR( α )
∂α
para que los puntos de corte estén lo más próximos posible. Si se consigue que ambos puntos se
confundan en uno, habremos alcanzado la astatización perfecta.
En los sistemas astáticos, las posiciones de equilibrio estable e inestable deben estar muy próximas.
∂M g ∂M R
Actuando sobre el sistema podemos conseguir que = llevando el sistema a la posición
∂α ∂α
de equilibrio indiferente.
En este caso, las curvas son tangentes y se consigue la astatización perfecta. Esto no ocurre
exclusivamente en un punto, sino en un pequeño entorno.
Demostración de las condiciones necesarias para la astatización:
∂M g ∂M g ∂M R ∂M g  ∂M g ∂M R 
M = Mg − MR ⇒ dg + dα − dα = 0 ⇒ − dg =  − dα

∂g ∂α ∂α ∂g  ∂α ∂α 
∂M g ∂M g
− −
dα ∂α ∂g ∂g
definiendo la sensibilidad como : ∂g = ∂M =
∂M
dg  g ∂M R 
 
 ∂α − ∂α  ∂α
 
dα ∂M g ∂M R
en el equilibrio indiferente: =0 ⇒ = astatización perfecta
dg ∂α ∂α
Hemos de buscar el momento recuperador cuya dependencia con la variable α sea no lineal (de ahí
el nombre de este tipo de gravímetros); es decir, que su representación no sea una recta porque en
tal caso es más difícil conseguir la astatización perfecta
El gravímetro LaCoste - Romberg utiliza una masa de lOgr. Al ser de tipo no lineal, el ángulo α + δ
= 90 en su posición de lectura (brazo horizontal), siendo δ ≈ 0 (del orden de segundos).
Vimos que si δ = O y L0 = O los momentos eran:
Página 74
∂M g
M g = m.g.a.senα ⇒ = m.g.a. cos α
∂α
∂M R
M R = k.b.d.senα ⇒ = k.b.d. cos α
∂α
vemos que su relación con α no es lineal, ya que aparece la función seno.
∂M g ∂M R
La sensibilidad será máxima s¡ − = ( m.g.a − k.b.d ) cos α = 0 para cualquier valor de α , y
∂α ∂α
dα
esto ocurrirá si (m.g.a − k.b.d ) cos α = 0 ⇒ = ∞ . Por tanto esta es la condición de
dg
astatización.
6.9.2.1.- Correcciones a las medidas gravimétricas

Factores a controlar en un gravímetro:
1- Temperatura: al variar la temperatura varían las propiedades elásticas del sistema (constantes
elásticas de los resortes, longitud de los brazos y de las palancas del mismo). Esto, produciría una
alteración en la posición de equilibrio del sistema.
Elegir materiales que no se vean afectados por la temperatura es utópico. Para evitar la influencia de
la temperatura, hay dos posibilidades:
- Colocar el sistema elástico en el interior de un termostato, que haga que la temperatura se
mantenga constante. Son sistemas termostatizados, llevan una bateria para mantener la temperatura
de régimen constante. Los gravímetros termostatizados tienen que colocarse con tiempo
(aproximadamente 24 horas) en condiciones de funcionamiento y controlar la temperatura
constantemente ya que si no están a la temperatura de régimen no se comportan adecuadamente.
- Cuando no hay termostato, se diseñan montajes termocompensados en los que si una
parte del sistema sufre dilatación, esta deformación queda compensada con la deformación de otra
zona.
2- Presión atmosférica: para eliminar su influencia el sistema elástico del gravímetro se coloca en
recipientes herméticamente cerrados y de baja presión, a los que no puede acceder directamente el
usuario.
Para controlar la presión, en los sistemas a rotación se colocan compensadores barométricos,
disponiéndose el sistema elástico según la figura:
Hay que conseguir que el llamado momento volumétrico del sistema respecto al eje de giro sea
nulo, para ello se sitúa en el extremo opuesto del brazo una masa k (qué actúa como compensador
barométrico) extremadamente ligera, y cuyo volumen es mayor que el de la masa m, de manera que:
v k rk = rm v m donde rk < rm y v k < v m
Se trata de evitar que aparezcan presiones o fuerzas que pudieran ejercer algún momento sobre el
eje de giro, alterando así el equilibrio del sistema. Se diseñan por tanto, sistemas que actúen como
compensadores barométricos.
3- Inclinación del eje de giro del gravímetro: la barra y el eje de giro del gravímetro tienen
posiciones perpendiculares y hay que controlarlas. El gravímetro LaCoste-Romberg tiene un nivel
para controlar la horizontalidad de la barra a la hora de realizar la lectura y otro para controlar la
horizontalidad del eje de giro.
Si el eje de giro no está horizontal, sino que está inclinado, la componente horizontal del momento
respecto al eje de giro estará influenciada por β si este ángulo es distinto de cero. El gravímetro nos
da por tanto una alteración en el valor de la gravedad. Aparece la llamada variación aparente de la
gravedad, que responde a una expresión del tipo:
Página 75
β2
∆g ≈ −g
2
Esta variación llevada a una gráfica nos representa una parábola, cuyo mínimo se encuentra en
β = 0:
dM
dM = M 0 senβdβ ⇒ = M 0 senβ
dβ
dM
= mínima ⇒ senβ = 0 ⇒ β=0
dβ
El sistema debe estar perfectamente nivelado y el eje debe ser rigurosamente horizontal.
4- Campo magnético exterior: para evitar la influencia del campo magnético externo, los
sistemas elásticos del gravímetro se fabrican con materiales que no sean magnético (cuarzo, sílice
fundida). Al fabricar el gravímetro, los muelles se someten a un proceso de desmagnetización y se
usan pantallas que lo aíslan del campo magnético externo.
6.9.2.2.- Correcciones:
• Lectura:
El resultado de la medida con un gravímetro es una lectura en divisiones de una escala. La lectura
es función de la gravedad: L =f(g) o g = F(L) siendo F(L) la función de calibrado del gravímetro
que es un dato facilitado por el constructor (muchas veces no es una función matemática, sino una
tabla de calibración). Es necesario conocerla para poder utilizar el gravímetro. Esta función puede
variar para diferentes rangos de lecturas. El gravímetro es calibrado por el constructor, pero el
usuario puede hacer una comprobación estacionando en puntos de gravedad conocida, existiendo
para ello líneas y polígonos de calibración. De la relación entre las lecturas del gravímetro y ¡os
valores de la gravedad en los puntos de estación, deducimos las constantes unitarias
correspondientes para los diferentes rangos de lecturas del instrumento.
• Deriva instrumental o deriva del punto 0:
La mayoría de los gravímetros utilizan un sistema de lectura que se llama de vuelta a cero o
posición de cero. Consiste en que, bajo la acción de una fuerza, suplementaria (la que realiza el
tornillo sobre el muelle), se lleva el sistema elástico a la misma posición de referencia. En la lectura
se controla la modificación de esa acción suplementaria; pero no se observa directamente, sino que
vemos la variación en el tornillo micrométrico calibrado para que la barra horizontal se mantenga
en dicha posición.
En el gravímetro LaCoste - Romberg el tornillo modifica el punto de unión del resorte de longitud
cero a la carcasa. Se varía la distancia d, de tal manera que barra del gravímetro se mantenga
siempre horizontal. La mordaza del brazo se libera exclusivamente a la hora de ir a realizar las
lecturas; es decir, una vez estacionado y nivelado. Sabremos que el balancín se encuentra horizontal
por métodos ópticos, bien porque se observa la trayectoria de un rayo reflejado sobre un espejo
situado en la barra, o bien mediante la observación de la imagen de un hilo adherido al brazo, de
manera que dicha imagen ocupa una cierta posición sobre la escala cuando el brazo se encuentra en
posición horizontal. El cambio de posición de dicha imagen se controla mediante la actuación sobre
el tomillo micrométrico. La fracción o número de vueltas ejercidas sobre dicho tornillo nos dará la
lectura.
Se define deriva como la variación de las lecturas del gravímetro con el tiempo en un determinado
lugar. La causa de la deriva de los gravímetros de tipo elástico es principalmente la fatiga que
experimentan los elementos elásticos del gravímetro. Estas variaciones en los valores de la
gravedad observados a diferentes intervalos de tiempo se deben a que el sistema elástico se deforma
y se recupera según el sistema. También a la influencia de la atracción lunisolar o marea
Página 76
gravimétrica, y no a que la gravedad esté variando, como pudiera parecer.

Si llevamos todas esas lecturas a una gráfica, vemos que se ajustan generalmente a una línea recta
Lf − Li
dependiente C D = . El valor de esta pendiente es la deriva unitaria, o variación de lectura
tf − ti
por unidad de tiempo. Antes de llevar estas observaciones a una gráfica y calcular esta corrección
por deriva, deberíamos eliminar los efectos de la marea gravimétrica. De cualquier manera, en caso
de no hacerse, este efecto de la atracción lunisolar queda en parte compensado al realizar la
corrección por deriva de este modo.
Por ello si el intervalo de tiempo es grande, suponemos que los efectos por marea gravimétrica han
variado considerablemente entre la primera y última observación, y no vamos a hacer corrección
lunisolar, la mejor manera de tener compensadas las lecturas consistirá en calcular y realizar
correcciones por deriva independientes para todas aquellas observaciones comprendidas dentro de
un rango temporal que establezcamos, por ejemplo de dos horas. Para poder hacer esto habremos de
repetir una misma observación cada dos horas.
Los efectos de la marea gravimétrica son del orden del microgal. Dependen principalmente le la
latitud del punto, de la época del año y de la hora a la que haya sido realizada la observación. El
gravímetro LaCoste - Romberg incorpora software para el cálculo de la corrección compensadora
en base a estos parámetros. Esta compensación suele realizarse para medidas de gravedad de alta
precisión. No se utiliza en prospección gravimétrica en trabajos de poca precisión, o cuando las
estaciones se encuentran próximas y el intervalo el tiempo de las observaciones es pequeño, en cuyo
caso suponemos que todas están afectadas del mismo modo.
La lectura adecuada es la primera lectura observada. Para obtener la lectura correcta, es decir, la
primera Li, a partir de cualquier otra Lf, no tenemos más que hacer:
Li = Lf- CD(tf- ti).
Podemos controlar la corrección por deriva repitiendo al menos la observación en una estación.
• Gravímetros astáticos de tipo Worden
Es también astático. El resorte es vertical y no está directamente unido al brazo principal del
gravímetro, sino a un brazo menor complementario.
• Gravímetros astáticos tipo LaCoste- Romberg
El resorte es metálico. Está termostatizado, siendo necesario para trabajar con él encontrándose a
una temperatura constante llamada temperatura de régimen o temperatura funcionamiento (53 ºC).
Para mantener esta temperatura es necesario un termostato y una batería. Como es un gravímetro
metálico, debe estar desmagnetizado y apantallado respecto campos magnéticos externos (para que
no le influya el campo magnético terrestre ni otros que puedan ser generados por centrales
eléctricas, cables de alta tensión, etc.). En el interior existe un alto grado de vacío para controlar la
presión.
Página 77
7.- SISTEMAS DE ALTITUDES
7.1.- Introducción. Números o cota geopotencial
Al utilizar un nivel, mediante visuales horizontales efectuamos lectura sobre dos miras verticales
(perpendiculares a la superficie equipotencial que pasa por dichos puntos), situadas en dos puntos,
A y B próximos entre sí.
Esta diferencia de lecturas la interpretamos como un desnivel geométrico y decimos que
corresponde a una diferencia de altitudes entre A y B.
Si los puntos se encuentran muy separados, el desnivel geométrico que existe entre ellos se calcula
mediante una serie de niveladas ∆n AB = ∑ δn , nos indicaría una diferencia entre las altitudes
ortométricas de los puntos A y B; pero no es así, ya que, en este caso las superficies equipotenciales
que pasan por dichos puntos no son equidistantes, por lo que: ∆n AB ≠ ∆H AB .
Si consideramos que las variaciones o incrementos de nivel que medimos δn son infinitesimales,
es decir, δn ≈ dn , podemos considerar el sumatorio como una integral:
B
∆n AB = ∑BA δn = ∫A dn
Cambiando el itinerario de nivelación entre los puntos, se puede demostrar que los desniveles
obtenidos ∆n AB son distintos.
Igualmente si hacemos un anillo cerrado en el punto A, ∆n AA = ∑AA δn ≈ ∫ dn ≠ 0 el desnivel obtenido
no es nulo; es un error apreciable, no sólo debido a los errores propios de las mediciones y depende
del itinerario. Esto es debido a que la dirección de la vertical es perpendicular a las superficies
equipotenciales y éstas no son equidistantes.
En este caso δn es geométricamente distinto δH B dado que los valores de la gravedad en los puntos
n y B’ sin distintos.
En el punto n: δW = −g.dn
Página 78
En el punto B’: δW = −g '.dH B

Como ambos puntos están sobre la misma superficie equipotencial W:
g
− g.dn = −g'.dH B ⇒ dH B = δn y como g ≠ g' ⇒ dH B ≠ δn
g'
Si pudiésemos medir diferencias infinitesimales podríamos considerar la igualdad
δW = −g.δn ≡ dW = −g.dn . Son diferencias de potenciales en base a los cuales podremos definir
sistemas de altitudes y calcular diferencias de altitudes entre puntos.
Vemos que el valor de la integral, que representa un itinerario entre A y B en el cual se han
realizado medidas conjuntas de gravedad y nivelación,. no depende del itinerario elegido, dado que
la función potencial depende exclusivamente de sus valores en los puntos inicial y final (puntos
extremos); y por lo tanto, el valor correspondiente a un itinerario cerrado será:
B B B
∫A dW = − ∫A g.dn ⇒ − ∫A g.dn = WB − WA
A
− ∫A g.dn = WA − WA
Se adopta una magnitud 'que represente a la diferencia entre los potenciales de un punto P (Wp) y un
punto de referencia O (W0) situado sobre el geoide, y que por lo tanto, tendrá su mismo potencial.
Esta magnitud es la llamada cota o número geopotencial.
P
Número geopotencial o cota geopotencial: C = W0 - WP = ∫0 g.dn
Por tanto la cota geopotencial de un punto sobre el geoide será cero.
Para conocer su valor, habría que realizar un itinerario de nivelación combinado con medidas de la
gravedad desde el punto O a nivel del mar hasta el punto considerado P. A cualquier punto
podremos asociarle una cota o número geopotencial.
La cota geopotencial viene expresada en las llamadas unidades geopotenciales:
1 u.g.p. = 1 kgal x m = 1 gal x km
A partir de la cota geopotencial pueden definirse diferentes sistemas de altitudes.
7.2.- Altitudes dinámicas.

din cA
La altitud dinámica correspondiente al punto A será H A = . con resultado en Km si γ 0 está en
γ0
gales y cA en u.g.p.
din
c
Genéricamente escribiremos: H = donde γ 0 = 980,6294 gales es una constante
γ0
correspondiente al valor de la gravedad normal a 45 grados de latitud, calculada utilizando la
fórmula internacional de la gravedad normal del sistema geodésico de referencia utilizado; en este
caso el del año 80. (Si no se especifica presupondremos que se trata del S.G.R. del 80).
Tiene dimensiones de longitud, pero no tiene de por sí un significado geométrico claro. Si nos
fijamos en la expresión vemos que es como realizar un cambio de escala y de unidades de las cotas
geopotenciales. Tiene la ventaja de venir expresado en unidades de longitud.
Si el incremento de altitudes dinámicas entre dos puntos A y B es nulo,
∆H din din din
AB = H B − H A = 0 podemos asegurar que ambos puntos se encuentran sobre la misma
superficie de nivel.
Las altitudes dinámicas no han sido utilizadas en Geodesia para dar la coordenada altimétrica de un
punto. Han sido utilizadas las llamadas altitudes ortométricas, que son aquellas a las que
comúnmente hacemos referencia al hablar de la altitud geométrica de un punto.
Página 79
7.3.- Corrección dinámica:

Pone de manifiesto que el incremento geométrico de la nivelación a lo largo de un itinerario entre
los puntos A y B: ∆n AB = ∑BA δn ≠ ∆H din
AB es distinto que el incremento de sus altitudes dinámicas.
cB cA 1
∆H din din din
AB = H B − H A = − = (c B − c A )
γ0 γ0 γ0
B B
c = ∫A g.dn = ∫A dc = c B − c A
1 B 1 B 1 B 1 B
∆H din
AB = ∫A g.dn = ∫A (g − γ 0 + γ 0 )dn = ∫A (g − γ 0 )dn + ∫A γ 0 .dn
γ0 γ0 γ0 γ0
B B B
∆H din
AB = ∫A dn + ∫A (g − γ 0 )dn donde ∫A dn ≈ ∆n AB
B
∆H din
AB = ∆n AB + ∫A (g − γ 0 )dn
∆H din
AB = ∆n AB + CD AB
B (g − γ 0 )
La corrección dinámica: CD AB = ∫A dn es la corrección que hemos de aplicarle al desnivel
γ0
geométrico para obtener la diferencia de altitudes dinámicas entre dos puntos. Cuando los puntos
están próximos, estamos considerando que g no varia entre A y B; o lo que es lo mismo, estamos
considerando que este término es nulo.
En la práctica, se aproxima la integral a un sumatorio, dado que no podemos realizar itinerarios de
nivelación con un número infinito de niveladas:
B (g − γ0 ) (g − γ 0 )
CD AB = ∫A dn ≈ ∑BA δn
γ0 γ0
Vemos que para calcularlo hemos de ir midiendo en cada punto del itinerario el valor de la
gravedad g, y su incremento geométrico respecto del punto anterior δ n.
7.4.- Reducción de la gravedad de Poincaré y Prey

El objetivo de esta reducción de la gravedad observada es poder deducir la gravedad en puntos del
interior de la tierra, situados sobre la línea de la plomada que pasa por un punto P de la superficie, a
partir del valor de la gravedad medido en éste.
Al igual que cuando queríamos conocer el valor de la gravedad reducida al geoide, ahora queremos
deducirlo reducido a la superficie equipotencial que pasa por Q.
Consideramos que el punto P descansa sobre una lámina Bouguer, plana e indefinida, de densidad
gr
constante ρ 0 = 2,67 y cuyo espesor es la diferencia de altitudes entre ambos puntos:
cm 3
h = HP − HQ
Seguidamente eliminamos (por lo tanto restamos) el efecto de atracción que dicha lámina ejerce
sobre la gravedad observada en P: 2πGρ 0 h . El término será: -0,1119h.
Aplicamos la corrección Aire-Libre para bajar de P a Q: 0,3086h.
Tendremos tras aplicar estas reducciones, el valor que tendría la gravedad observada en el punto Q
si no existiese terreno sobre él. Pero no es así, de manera que hemos de reducir nuevamente este
valor.
Ahora le sumaremos el efecto de atracción que ejerce la lámina Bouguer de espesor h = H P − H Q ,
pero que en este caso, está situada por encima del punto y, por lo tanto, el término para aplicar esta
reducción es de nuevo: -0,1119h. Reduce el valor de la gQ dado que su efecto atractivo es ahora
opuesto a la dirección de la gravedad.
Página 80
No hemos de realizar corrección alguna debida a los efectos de aquellas masas que se encontrasen
por encima de la lámina Bouguer dado que el valor de la gravedad observada en P, del cual
partíamos para deducir el de la gravedad en Q, ya estaba influido por los efectos de dichas masas.
La reducción será: gQ = gp + (0,3086-0,1119-0,1119)h
g Q = g P + 0,0848(H P − H Q )
Si se requieren grandes precisiones en el cálculo de la altitud ortométrica, hay que incluir
correcciones topográficas a la hora de deducir los valores de gravedad en puntos del interior.
Para generalizar esta fórmula ponemos la altura del punto genérico Q (HQ), como término variable,
es decir, sustituimos HQ por z y dH por dz. La gravedad es una función g(z)
g Q = g(z) = g P + 0,0848(H P − z) con H en km, g en gales
7.5.- Altitudes ortométricas. Corrección ortométrica

Se define la altitud ortométrica de un punto como la distancia medida desde el geoide hasta dicho
punto sobre la línea de la plomada que pasa por el mismo.
P P
Definíamos CP = W0 - WP = ∫0 g.dn = ∫P0 g.dn pero como, tanto el punto O como el Po están sobre el
geoide, en ambos el valor del potencial es el del geoide; y por lo tanto, la diferencia de potencial del
punto P respecto a cualquiera de ellos será la misma. Pero en el caso de querer hallar el potencial en
P partiendo del punto Po utilizando la línea de la plomada como itinerario, habremos de medir
gravedad en puntos Q del interior de la tierra. Al no ser esto posible, los valores de la gravedad
correspondientes a dichos puntos serán deducidos a partir del observado en el punto P. Por otro
lado, vemos que los incrementos de nivel dn son ahora incrementos de altitud ortométrica dH.
P P P P
CP = W0 - WP = ∫0 g.dn = ∫P0 g.dn = ∫P0 g.dH = ∫0 g.dH
despejando H:
dC C dC
g.dH = dC ⇒ dH = ⇒ H = ∫0 ⇒ Altitud ortométrica en función del
g g
geopotencial
dW W dW
g.dH = −dW ⇒ dH = − ⇒ H = − ∫W ⇒ Altitud ortométrica en función del
g 0 g
potencial
El signo negativo es debido a que el sentido de crecimiento de W es contrario al de crecimiento de
H.
Estas fórmulas no son manejables en la práctica. Fijémonos en que g no es el valor de gp en el punto
P, sino que es una variable.
H 1 H
C = ∫0 g.dH multiplicando y dividiendo por H tendremos C = H. ∫ g.dH
H 0
Sabemos que el valor medio de una función viene expresado por la fórmula
1 b 1 H
f (x) = ∫a f ( x )dx por tanto, g = ∫0 g.dH es el valor medio de la función g entre los valores P0 y
b−a H
P.
C CP
C = Hg ⇒ H= particularizando para un punto P: H P =
g g
Esta es la fórmula práctica, pero hay que saber cómo calcular la gravedad media, dado que no es
una cantidad constante sino que depende del punto elegido P.
Para el cálculo de esta g bastaría con conocer la gravedad en la superficie terrestre y la gravedad
en el geoide (reducida al geoide por los métodos de reducción vistos anteriormente) y hacemos la
Página 81
media aritmética entre las dos:

g Pg 0
g=
2
También podemos tener está g mediante la reducción Prey simplificada:
g  z =g 0, 0848 H −z  ,
donde g es la gravedad medida en el punto P en el terreno. Ahora la integración puede efectuarse
inmediatamente:
[ ]
H H
1 1 z2
g= ∫ [ g0,0848  H −z  ] dz ⇒ g= g⋅z 0,0848  H⋅z −  ⇒
H 0 H 2 0
g= [
1
H
g⋅H 0,0848 H⋅H −
H2
2
 ] ⇒ g=g 0,0848 H −  ⇒
H
2
H
g=g 0,0848  ⇒ g=g 0,0424⋅H
2
Tendremos por tanto:
C C
H= ⇒ H= , con C en u.g.p., g en gals y H en km.
g g0,0424⋅H
Para resolver este sistema y clacular el valor de H podremos realizar un proceso iterativo, pero lo
normal es resolverlo normalmente, puesto que es una ecuación de segundo grado. Tiene dos
soluciones y sólo una de ellas nos sirve.
7.6.- Corrección ortométrica:

La diferencia de altitudes ortométricas entre dos puntos A y B, será:
∆H AB = H B − H A
Como ∆n AB ≠ ∆H AB , habremos de añadir un término correctivo COAB de manera que:
∆H AB = ∆n AB + CO AB
Para el cálculo de la corrección ortométrica, partimos de ∆H AB = H B − H A . Sumamos y restamos
las altitudes dinámicas de A y B:
∆H AB = H B − H A + H din din din din
B − HB + HA − HA
Asociando términos por conveniencia:
∆H AB = (H B − H din din din din
B ) − ( H A − H A ) + (H B − H A )
∆H AB = (H B − H din din
B ) − ( H A − H A ) + ∆H AB
Recordemos que ∆H AB = ∆n AB + CD AB , de modo que si sustituimos este término, ya tenemos la
relación entre ∆H AB y ∆n AB , de manera que los términos entre paréntesis junto con la CD AB
constituirán la corrección ortométrica.
Calculamos los valores de los términos entre paréntesis para sustituirlos posteriormente en esta
din din din
expresión: H B − H B0 = ∆n B0B + CD B0 B donde H B0 = 0 .
Si imaginamos un itinerario de nivelación ficticio entre los puntos Bo y B que siga la línea de la
plomada que pasa por ambos, en tal caso ∆n B0 B = H B , y los diferenciales dn serán dH.
Por tanto: H din
B = H B + CD B0 B ⇒ H B − H din
B = −CD B0 B
B
Lo demostramos: CD B0 B = ∫B0
g − γ0
γ0
dH =
1 B
γ0
[ B
∫B0 g.dH − γ 0 ∫B0 dH ]
Poniendo los límites correspondientes en función de la variable de integración:
Página 82
CD B0 B =
1 HB
γ0
[ H
∫0 g.dH − γ 0 ∫0 B dH ]
1 H H H
Recordemos que g = ∫0 g.dH ⇒ g.H = ∫0 g.dH y por similitud g B0B H B = ∫0 B g.dH
H
Por tanto: CD B0 B =
1
γ0
[ H
g B0 B H B − γ 0 ∫0 B dH =
1
]
g H B − γ 0 .H B
γ 0 B0 B
[ ]
CD B0 B =
1
γ0
[ ]
gB B − γ0
g B0 B − γ 0 H B = 0
γ0
HB
B
g B0 B − γ 0 B
g B0 B − γ 0
CD B0 B = ∫B dH ≈ ∫B H
0 γ0 0 γ0 B
HB − H din
B = −CD B0 B
Por analogía: H A − H din

A = −CD A 0 A
Por otra parte: ∆H din

AB = ∆n AB + CD AB
Sustituyendo: ∆H din
AB = ∆n AB + CD AB − CD B0 B + CD A 0 A
∆H din
AB = ∆n AB + CO AB siendo CO AB = CD AB − CD B0 B + CD A 0 A
Es decir, el cálculo de esta corrección se basa en calcular correcciones dinámicas.
Para la práctica, poniendo g = g A 0 A y g = g B0 B :
(g − γ 0 ) (g − γ 0 ) (g − γ 0 )
CO AB ≈ ∑ BA δn − B HB + A H A siendo γ 0 una constante.
γ0 γ0 γ0
Vemos que para calcular esta corrección necesitamos conocer previamente unos valores
aproximados para las altitudes ortométricas de A y B, y utilizar métodos interactivos hasta obtener
la solución.
7.7.- Altitudes normales. Corrección normal

Altitud normal (H*) se define en el campo de la gravedad normal de manera análoga a como lo está
la altitud ortométrica en el campo de la gravedad real. Muy utilizada en los países del este.
En este caso la cota geopotencial vendrá dada por:
H*
C = U 0 − U H = ∫0 γdH *
Multiplicando y dividiendo por H* (altitud normal):
1 H*
C= H* ∫ γdH * = γH
H* 0
C
Fórmula práctica: H = siendo γ el valor medio de la gravedad normal entre Q0 y Q.
γ
Trazamos por el punto P la llamada plomada normal, que es la normal N’ al elipsoide, que pasa por
P. Sobre dicha normal N’ existirá un punto Q tal que el valor de su potencial normal UQ será igual a
Wp ; valor del potencial correspondiente a la gravedad real en el punto P.
Igualmente podríamos tomar otros puntos P',P",... y sobre sus normales obtendríamos igualmente
puntos Q', Q",... tales que UQ’ = Wp’, UQ’’ = Wp’’,...
Llamamos teluroide a la superficie formada por todos estos puntos Q, tiene una forma equivalente a
la de la superficie terrestre.
Página 83
La altitud normal asignada a un punto P, situado sobre la superficie terrestre, es la distancia medida
sobre la normal elipsoidal desde Qo hasta Q: H *P = Q 0 Q
Vemos que la altitud normal no es igual a la elipsoidal: H *P ≠ h p
h p = Q0P = Q0Q + ζ , h p ≈ H P + N P
siendo N ≈ N P y ζ la anomalía en altitud.
La superficie resultante de llevar, sobre la normal elipsoidal N’ de cada punto P de la superficie
terrestre, su altura normal H* hacia el interior de la tierra, es el denominado cuasigeoide. No es una
superficie física real como el geoide, aunque es muy próxima a él. .El cuasigeoide también puede
obtenerse llevando las anomalías en altitud sobre dicha normal a partir del elipsoide.
La variación de γ con la altura es:

 2 
γ = γ 0 1 − (1 + α + m − 2αsen 2 ϕ)h + ...
 a 
α = aplanamiento geométrico
h = altitud elipsoidal
γ 0 = gravedad normal en el elipsoide
ω2 a 2 b
m= es la constante geodinámica
GM
 2 2 
Genéricamente: γ = γ (z) = γ 0 1 − (1 + α + m − 2αsen ϕ)z + ...
 a 
y la gravedad media:
Página 84
1 z 1 z  2 2 
γ= ∫0 γ (z)dz = ∫0 γ 0 1 − (1 + α + m − 2αsen ϕ)z  dz
z z  a 
γ z 2 z 
γ = 0  ∫0 dz − (1 + α + m − 2αsen 2 ϕ) ∫0 zdz 
z  a 
γ  2 z  γ 
2
1 
γ = 0 z − (1 + α + m − 2αsen 2 ϕ)  = 0 z − (1 + α + m − 2αsen 2 ϕ) z 2 
z  a 2  z  a 
γ0  2 2 z2  γ0  1 
γ= z − (1 + α + m − 2αsen ϕ)  = 1 − (1 + α + m − 2αsen 2 ϕ) z 
z  a 2  z  a 
La gravedad media en este caso será:
1 H*
γ= * ∫0
γ.dH *
HP
 1 
γ = γ 0 1 − (1 + α + m − 2αsen 2 ϕ)H *P + ...
 a 
Corrección normal:
Por analogía con la corrección ortométrica pero en el campo de la gravedad normal:
∆H *AB = ∆n AB + CN AB
(g − γ ' 0 ) (γ − γ'0 ) (γ − γ'0 )
Se justifica que: CN AB ≈ ∑BA δn − B H *B + A H *A
γ'0 γ'0 γ'0
siendo γ ' 0 una constante correspondiente a la gravedad normal sobre el elipsoide ( γ ' 0 = γ 45 ) y g una
gravedad real medida.
7.8.- Comparación de diferentes sistemas de altitudes

La altitud será el número geopotencial dividido por una gravedad media:
c c
Altitud = =
gravedad g
Dependiendo del valor-tomado para la gravedad tendremos un sistema de altitudes u otro. Unos con
significado geométrico claro y otros no.
c c
H din = = Dinámica: no tiene significado geométrico claro; pero si dos puntos tienen la
γ 0 γ 45
misma altitud dinámica, tienen la misma cota geopotencial y, por tanto, están sobre la misma
superficie equipotencial.
c
H= Ortométrica: dos puntos que tengan la misma altitud ortométrica, no
g
necesariamente están en la misma superficie equipotencial, ya que la cota geopotencial será
diferente de uno a otro.
c
H* = Normal: la ventaja de la altitud normal H * frente a la ortométrica H está en que si las
γ
determinamos por estas fórmulas, para hallar la g en las altitudes ortométricas interviene la
densidad de la corteza (2,67 gr/cm3), pero los cambios de densidad por discontinuidades o fallas se
ven reflejados en los cambios de gravedad, cosa que no ocurre con la altitudes normales.
Página 85

Apuntes Geodesia Fisica 2005 PDF

Cargado por

Copyright:

Formatos disponibles

Apuntes Geodesia Fisica 2005 PDF

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Apuntes Geodesia Fisica 2005 PDF

Cargado por

Copyright:

Formatos disponibles

Apuntes de Geodesia Física Curso 2005-2006

Apuntes de Geodesia Física

Autora: Yolanda Rubio

4.1.- Potencial normal y potencial perturbador..............................................................................43

1.- INTRODUCCIÓN A LA GEODESIA FÍSICA

1.1.- Definición y objeto de la Geodesia Física

Cuando intervienen las observaciones astronómicas tenemos las coordenadas astronómicas:

La diferencia entre las dos altitudes es la ondulación del geoide: h=H+N

1.2.- Diferencia entre geodesia geométrica y geodesia física:

1.3.- Consideraciones históricas sobre la figura de la Tierra y sus métodos de

- Pitágoras hizo referencia sobre la redondez de la tierra.

Huyggens consideró toda la masa terrestre concentrada en su centro y calculó el aplanamiento

nivel medio de los mares en reposo.

2.- TEORÍA DEL POTENCIAL. EL CAMPO GRAVITATORIO

2.1.- Ley de Newton de la Gravitación Universal

2.2.- Campo de fuerza newtoniano. Intensidad de campo y potencial.

2.2.1.- Campo newtoniano debido a una masa puntual:

2.2.2.- Campo newtoniano debido a un sistema de masas puntuales:

2.2.3.- Campo newtoniano debido a un cuerpo extenso:

Generalizando la fórmula del potencial correspondiente a un elemento diferencial de masa:

2.3.- Definición y propiedades de la superficie equipotencial: intensidad de campo

2.4.- Flujo vector intensidad de campo. Teorema de Gauss

2.4.1.- Caso de una masa puntual interior m

2.4.2.- Caso de una masa exterior a la superficie:

elementales dos a dos se anulan.

2.4.3.- Caso varias masas interiores

Caso de tener varias masas interiores m int = ∑in=1 m i y el flujo será:

2.5.- Ecuaciones de Laplace y de Poisson. Funciones armónicas

2.6.- Campo y potencial de las fuerzas centrífugas

En el sistema de ejes cartesianos tiene las siguientes componentes: (Fc ) x = ω 2 x , (Fc ) y = ω 2 y y

2.7.- Gravedad y potencial gravitatorio terrestres

Propiedades de los potenciales:

2.8.- Superficies de nivel y líneas de la plomada:

La altitud ortométrica H es la distancia medida a lo largo de la línea de la plomada desde el geoide o

2.9.- Desarrollo del potencial gravitatorio terrestre. Polinomios de Legendre

Px 'y ' = ∫τ x ' y' dm

Sabemos que el momento de inercia de un sistema de partículas respecto a un eje es igual al

De manera que el término V2 nos quedará:

Sacando factor común R 2 y realizando un cambio de signo e introduciendo 1/2 en el paréntesis:

2.10.- Problemas de contorno de la teoría del potencial. Solución de la ecuación

2.10.1.- Teorema de Stockes:

2.10.2.- Primer problema de contorno de la teoría del potencial o problema de Dirichlet:

h (λ) = sen (mλ ) o h (λ) = cos(mλ ) P mn (cos θ)sen (mλ)

3.- EL POTENCIAL Y LA GRAVEDAD NORMAL. SISTEMAS

3.1.- La función geopotencial

equipotencial que más se aproxima al geoide cumpliendo las condiciones impuestas. Uo es

3.2.- Esferoide de Bruns: teoría del primer orden de aproximación.

3.3.- Primera fórmula de Clairaut en aproximación de primer orden

como sen 2 θ = (1 − cos 2 θ) ⇒ R

cos2  1−2⋅sen2 cos2 

3.4.- Gravedad normal en el esferoide de Bruns o elipsoide de nivel:

3.5.- Esferoide de Helmert:

3.6.- El elipsoide de nivel

3.7.- Teoría de Pizzetti y Somigliana. Fórmula internacional de la gravedad

3.8.- Ventajas del uso del elipsoide como figura de referencia:

3.9.- Gravedad normal en el elipsoide:

3.10.- Sistema Geodésico de Referencia

3.10.1.- Sistema Geodésico de Referencia 24-30 (1924-1930):