Aplicaciones de La Geodesia
Aplicaciones de La Geodesia
Aplicaciones de La Geodesia
Objetivo y Aplicaciones
CONTENIDO PÁG.
2
Definición, división y objetivo de la Geodesia
1. ¿Qué es la Geodesia?
La palabra Geodesia literalmente expresa división de la Tierra, sin embargo, diversos
autores notables establecen distintas definiciones de este concepto. Para unos existe una clara
diferencia entre la Geodesia Teórica y la Geodesia Práctica, indicando que la primera estudia
la forma y dimensiones de la Tierra, en cambio la segunda establece los procedimientos para
la medida de porciones terrestres. Para otros autores esta diferencia no es tan clara, por ello se
refieren a la Geodesia como una ciencia cuyo objetivo es el de proporcionar un armazón o
estructura geométrica precisa para el apoyo de los levantamientos topográficos.
Podemos concluir este apartado diciendo que la Geodesia es una ciencia, que desde la
antigüedad, se ha dedicado al estudio de la medida y forma del globo terráqueo, adatándose a
las necesidades de la época para aplicarse a problemas prácticos, como son básicamente la
confección de mapas nacionales e internacionales, así como la preparación de cartas para
aplicaciones específicas como las geológicas e hidrográficas, entre otras. Pudiendo afirmar
que la Geodesia se ha necesitado y seguirá siendo necesaria mientras se proyecten obras
humanas que requieran precisiones cada vez mayores.
2. Problemas de la Geodesia
3
•3 Determinaciones astronómicas de las latitudes y longitudes terrestres.
3. División de la Geodesia
En esta situación de aislamiento, sucede que la esfera es una figura de equilibrio para
una masa homogénea en reposo, siendo el único movimiento posible para una masa
homogénea que se mueve como un sólido, el de una rotación uniforme alrededor de uno de
sus ejes principales de inercia. Ambas conclusiones, unidas al hecho de que una pequeña
rotación produce un achatamiento de la forma esférica, nos lleva a considerar que la Tierra
tiene una figura de equilibrio dada por un elipsoide de revolución, ligeramente achatado en los
polos, que gira alrededor de un eje que pasa aproximadamente por los polos. Conviene sin
embargo recordar, que la Tierra no es un cuerpo rígido y homogéneo, por tanto este modelo
teórico está muy simplificado. Otra simplificación importante que se ha indicado consiste en
ignorar los efectos gravitatorios de los demás cuerpos del sistema solar, incluido el Sol, sobre
la superficie terrestre. Es un fenómeno bien conocido la existencia de mareas tanto oceánicas
como terrestres y atmosféricas, todas éstas son debidas a la atracción gravitatoria que ejercen
4
sobre la superficie terrestre, principalmente el Sol por su gran masa y la Luna por su cercanía
a la Tierra.
•9 La Geodesia Regional, que es practicada por cada país con el fin de resolver los
problemas planteados por la Cartografía y la Geografía, entre otras.
•10 La Geodesia Topográfica, que trata de precisar detalles de una cierta superficie de
pequeñas dimensiones, para ello la considera una superficie plana o esférica según sean
sus dimensiones.
5
obteniendo el valor 1/50 parte de un círculo (7º 12'). Entonces, dado que la distancia entre
Alejandría
y Siene era bien conocida en la época
6
En lo que se refiere a nuestro país, desde la mitad del siglo XIX hasta la mitad del
siglo XX, se fueron desarrollando los trabajos destinados a extender una red geodésica
formada por cadenas de triángulos, los cuales corren sensiblemente a lo largo de los
meridianos, paralelos y costas (ver figura 2). Forman parte de esta red las cadenas de enlace
del Archipiélago Balear con la Península, así como los grandes cuadriláteros que ligan nuestra
costa con Argelia y Marruecos. El Archipiélago Canario posee su propia red que liga las islas
entre sí, y éstas con el continente africano (ver figura 2). Los trabajos geodésicos cuyo
objetivo fue la formación del mapa nacional en escala 1:50.000, se iniciaron en 1858
finalizándose en su conjunto hacia 1930. El punto fundamental de la red geodésica nacional o
dátum está en el Observatorio Astronómico de Madrid, cuyo meridiano ha sido utilizado
también como origen de longitudes. Sin embargo, la red geodésica de las Canarias por ser
independiente tiene su propio dátum, el cual se halla en el vértice denominado Pico de las
Nieves, en la isla de Gran Canaria.
Más recientemente y gracias al desarrollo de nuevas tecnologías como GPS, todas las
redes europeas se han podido unificar para formar una sola red de precisión. Esta unificación
ha permitido unir toda el área ibérica con Europa, mediante una red de alta precisión que nos
proporciona actualmente medidas de gran precisión para ciertos puntos, que pueden usarse
como vértices de referencia para realizar tareas científicas o técnicas. Esta unificación ha sido
posible gracias al desarrollo de la red EUREF (ver figura 3), la cual está formada por
estaciones GPS permanentes. Esta red de estaciones permanentes proporciona coordenadas de
alta precisión para diversas aplicaciones tanto científicas como técnicas.
Fig. 3. Mapa de la red de estaciones GPS permanentes que forman la red EUREF.
7
con una precisión centimétrica. Esta red de nivelación que recibe el nombre de UELN (United
European Levelling Network) está representada en la figura 4(a).
(b)
(a)
Fig. 4. (a) Mapa de la red de UELN formada por la unificación de todas las redes de
nivelación europeas. (b) Mapa de la red vertical europea EVRS (European Vertical Reference
System).
En el área ibérica esta red posee los vértices de precisión mostrados en la figura 5. Las
coordenadas y alturas de estos puntos, junto con los datos de los otros puntos situados por
toda Europa, pueden descargarse desde internet en la dirección:
http://crs.bkg.bund.de/evrs/tabelle_neu.html
8
Aplicaciones de la Geodesia Geométrica
En las campañas geodésicas, como las mencionadas anteriormente, se suelen determinar los
ángulos sobre la superficie terrestre, tras lo cual el problema consiste en calcular las longitudes de
los lados, para cada triángulo de la red. Según este procedimiento, si conocemos un lado de la
triangulación y determinamos todos los ángulos existentes en dicha red, es posible hallar las
dimensiones de la red completa (ver figura 1(a)). Por otra parte, la triangulación es también necesaria
en muchos casos, para pasar de una base medida a una base mayor o ampliada, debido a las
dificultades que puede presentar el terreno, para efectuar medidas de distancias sobre el mismo,
medidas que serían necesarias para determinar la base ampliada directamente (ver figura 1(b)).
Fig. 1. (a) Red geodésica entre dos puntos A y B. (b) Red geodésica de ampliación de base.
En consecuencia, es necesario llevar a cabo algún procedimiento que nos permita resolver los
triángulos planteados en una red geodésica. Para ello, existen diversos métodos de los cuales vamos
a comentar los más relevantes, como son el método de Legendre y el método de los aditamentos.
Comenzando por el primero, hay que decir que su nombre proviene del teorema que se aplica,
teorema de Legendre, este teorema permite asociar a cada triángulo geodésico como el mostrado en
la figura 2(a), un triángulo plano como el mostrado en la figura 2(b), de tal forma que las distancias
(a, b, c) en ambos triángulos son las mismas.
(a) (b)
9
En consecuencia, nos planteamos el problema de calcular la diferencia que existe entre
los ángulos (A, B, C) del triángulo geodésico, y los ángulos del triángulo plano asociado (A 1,
B1, C1), porque si conocemos estas diferencias, podemos establecer una relación entre ambos
triángulos, de tal forma, que es posible aplicar las conocidas y sencillas fórmulas de
trigonometría plana, a los triángulos geodésicos de la red considerada.
Por otra parte, hay que tener en cuenta que los triángulos de una red geodésica, poseen
unos lados cuyas longitudes rara vez superan los 100 km, en este caso podemos considerar los
triángulos geodésicos como triángulos esféricos, para los cuales es válida la ley de los
cosenos. Aplicando esta ley al triángulo de la figura 2(a) tenemos
a bc bc
cos = cos cos + sen sen cosA
R RR RR
si aplicamos ahora el desarrollo en serie de Taylor de las funciones trigonométricas, que aparecen
en esta expresión, se tiene
cosA 4 4
+ + +
2R2 24R4 2R2 24R4yV 2R2 24R4y
b c c3
^
R 3
6R y R
ab habiendo
despreciado los términos de quinto grado y superiores, en las cantidades , y
c RR
R Así, después de multiplicar y agrupar términos obtendremos
cosA 2 2 2 4 4 4 2 2 22 2 2
+
b +c -a a +b + c - 2ab - 2a c - 2cb
2bc 24R2bc
Ahora bien, para el triángulo plano de la figura 2(b), aplicando la ley de los cosenos, podemos
escribir
a2 = b2 + c2 - 2bc cos A1
de donde
cos A1 b2 + c2 - a2 2bc
10
A-A
22 2A - a4 - b4 - c4 + 2a2b2 + 2a2c22A+ 2c2b2
bcsen
cosA - cosA
sen A1 = 11-sen
=-2 cos A11 = -(A - A1)sen A1 = - 22 2
1
4b c
Entonces, combinando las expresiones obtenidas, con la fórmula anteriormente escrita para el
cos A, deducimos que
bcsen2 A
1
cosA = cosA1 - 2
6R
y sabiendo que
A-A A+A 1 1
cosA
1
- cosA = - 2 s e n s e n
22
dado que la cantidad A-A1 es muy pequeña comparada con A, podemos escribir dentro de la
precisión aceptada
1 11
2 2 6R2
concretamente, la diferencia de ángulos buscada A-A1 podrá escribirse como
bcsen A1
A - A1= 2 = (área del triángulo plano)/3R2
6R
recordando que el exceso esférico ε de un triángulo esférico, es precisamente el área del
triángulo dividida por R2, notamos que la corrección a considerar es la tercera parte del exceso
esférico. Análogamente, para los demás ángulos podremos escribir expresiones similares, de
tal forma que la corrección buscada para los ángulos del triángulo de la figura 2(a),
determinados sobre el terreno, será en cada caso restar un tercio del exceso esférico, que como
sabemos puede calcularse mediante
ε = A + B + C - 180º
En consecuencia, las ecuaciones
εε ε
A1-A=- B1-B=- C1-C=-
3 3 3
expresan el contenido del teorema de Legendre, para pequeños triángulos geodésicos (cuyos
lados no sobrepasen los 100 km).
No obstante, debemos notar que las ecuaciones anteriores, permiten resolver triángulos
geodésicos, es decir, calcular todas los lados de los triángulos geodésicos que forman una red,
aplicando correcciones a los ángulos determinados sobre el terreno. Ésta es una forma válida
de operar, pero también es posible aplicar correcciones a los lados, de tal forma que el
resultado sea, directamente, las distancias buscadas, es decir, los lados del triángulo
geodésico. Este método operacional, recibe el nombre de método de los aditamentos, siendo
el objetivo final del mismo, determinar las correcciones o aditamentos, que es necesario
aplicar a los lados de un triángulo, cuyos ángulos sean los determinados sobre el terreno (A,
B, C), con el objeto de convertirlo en un triángulo plano, al que podemos aplicar la
trigonometría plana para determinar sus lados.
1
1
sen b - sentí senB
R ^R^senA
donde aplicamos, como hicimos antes, los desarrollos en serie de las funciones
trigonométricas, considerando despreciables los términos de grado cuarto y superiores, en las
cantidades — y —, teniendo en consecuencia
R R
ab 3 3
RA
sen B 6R2J"R1 6R2JsenA
si designamos
a3 b3
— = a - A = a' b - — = b - A = b'
6R 2 a 6R 2 b
podemos escribir
b' = a'senB
sen A
análogamente para el lado c, tendremos
, sen C
c =a---------
sen A
De esta forma, introducimos una cantidad As (donde s designa el lado considerado), llamada
aditamento del lado s. Notando que las expresiones anteriores nos permiten escribir ahora
Para llevar a cabo la resolución de una red de triángulos, siguiendo el método de los
aditamentos, operamos como sigue:
15. Con el valor hallado a’ se calculan b y c’, empleando para ello la ley de los senos.
b'3 = b3 c'3 = c3
debido a que la diferencia b - b es muy pequeña frente a b.
4. Con los valores de los aditamentos obtenidos, podemos escribir finalmente los valores de
los lados buscados (b, c), en la forma
12
b = b′ + Ab c = c′ + Ab
6. Determinación de distancias
En cualquier campaña de triangulación, como las mencionadas anteriormente, se
suelen medir bases sobre el terreno, siguiendo para ello, un método que consiste en medir con
precisión una pequeña distancia, y a continuación, mediante una red de ampliación de base, se
obtiene la distancia entre los dos vértices correspondientes. Este procedimiento fue
mencionado en el apartado anterior, como una aplicación básica de la geodesia geométrica.
Llegados a este punto es conveniente especificar cómo es posible hallar una distancia sobre la
superficie de la Tierra, considerada en este caso como la superficie del elipsoide de referencia
utilizado durante esa campaña.
En primer lugar, hay que decir que la forma de llevar a cabo las medidas directas que
conlleva el cálculo de una base, ha cambiado mucho desde la antigüedad hasta nuestros días.
Esta evolución es continua, y por tanto, sigue sucediendo en actualmente, de tal forma que los
métodos a seguir estarán siempre sujetos a una continua revisión, para incorporar los recientes
avances habidos en este tema.
Para ser conscientes de esta evolución, basta recordar que la medida de bases con
precisión no fue realizada hasta el siglo XIX, fecha en la que se comenzaron a utilizar las
reglas bimetálicas, que fueron inventadas por Borda en 1792. El método que se utilizaba con
estas reglas, consistía en colocar cada regla, una y otra vez, a continuación de la otra. De esta
forma, podían medirse distancias de 10 km, como mucho, con una precisión de 10 -6. El
problema de este método está en lo laborioso del mismo y en la necesidad de conocer con
precisión la longitud de la regla empleada. Con respecto a esto último, la precisión con la que
se conoce la longitud de la regla empleada y cómo cambia esta longitud durante el proceso de
medida de bases, hay que recordar que esta longitud puede conocerse de forma precisa,
gracias a que la regla bimetálica está constituida por dos reglas de diferentes metales, sujetas
por un extremo, siendo el otro extremo libre. De esta forma, puede medirse el desplazamiento
relativo de una de ellas respecto de la otra, gracias a una regleta graduada colocada en la regla
inferior. Entonces, la alteración de las longitudes habidas por la distinta dilatación térmica de
estas reglas, puede medirse con precisión, hallándose también la verdadera longitud de la
regla bimetálica empleada mediante la fórmula
b = A + Bδ
donde A y B son constantes conocidas, b es la longitud de la regla bimetálica en metros
(aproximadamente 4 m) y δ es la lectura tomada en la regleta en mm.
En una versión posterior de este método, empleada por primera vez en 1885 por el
sueco Jáderin, las reglas bimetálicas fueron sustituidas por los hilos de invar, estos hilos se
colocan entre dos soportes y se tensan mediante dos pesas. Los hilos de invar están
compuestos de un material llamado invar, cuya composición es 36 % de níquel y 64 % de
acero; su longitud es de unos 24 m, con una sección recta de 1.65 mm de diámetro. Debido a
su composición, estos hilos presentan la ventaja de variar muy poco de longitud por efecto de
los cambios de temperatura, así las correcciones a efectuar por este efecto son siempre muy
pequeñas. Otra ventaja que presenta el uso de estos hilos, reside en que con ellos medimos las
distancias de 24 en 24 metros, mientras que con las reglas bimetálicas las distancias se miden
de 4 en 4 metros. Con todo, el método sigue siendo muy laborioso, ya que, es necesario que
no se produzca una desalineación, cuando se coloca el patrón una y otra vez, delante de la
13
posición anterior. Esto ha hecho que estos métodos antiguos hayan sido gradualmente
descartados, en favor de métodos más modernos de medida basados en la interferometría, o
también en la medida de los tiempos de retraso de una señal enviada y recibida por una misma
estación de registro.
Fig. 3. (a) Representación de las secciones normales recíprocas correspondientes a los puntos
A y B. (b) Representación de la línea geodésica que une los puntos A y B.
La figura 3(a) representa este problema, en ella vemos que para dos puntos dados
sobre la superficie de referencia, la cuerda es la línea representada por un trazo recto
discontinuo. En cambio, la verdadera distancia entre tales puntos considerados, debe
calcularse siguiendo la línea geodésica (curva de mínima distancia entre dos puntos, medida
sobre la superficie) representada en la figura 3(b), como puede verse esta línea está situada
dentro de las secciones normales recíprocas, representadas por las curvas AaB y BbA. La
sección normal AaB es una curva que se define mediante la intersección de la superficie de
referencia considerada (un elipsoide de revolución achatado por los polos), con el plano
definido por la dirección normal a dicha superficie de referencia en el punto A (dirección na) y
el punto B. Así, el corte de este plano y la superficie de referencia, es la curva AaB. De forma
similar definimos BbA intercambiando los papeles que juegan los puntos A y B.
14
Cuando consideramos distancias AB menores de 500 km, podemos aceptar que la
medida del arco sAB sobre la línea geodésica, es aproximadamente igual a la media de las
medidas sobre las secciones normales recíprocas, es decir
sAB = (sA + sB) / 2
obteniendo las cantidades sA y sB mediante el desarrollo en serie (Cid y Ferrer, 1997)
Así, realizadas todas las operaciones la expresión anterior nos dará la distancia
buscada. Debemos notar que cuando los puntos A y B sean muy próximos, la cuerda y el arco
serán prácticamente la misma cantidad. Esta situación suele ser frecuente en la medida directa
de bases sobre el terreno, donde las distancias consideradas no superan los 15 o 20 km.
7. Determinación de ángulos
Tal como puede verse en la figura 3(a), cuando nos situamos con un instrumento en el
punto A, observando el punto B, o al contrario, cuando observamos el punto A desde el B;
estamos determinando secciones normales recíprocas, de tal forma, que las medidas de los
ángulos, llevadas a cabo en una triangulación estarán falseadas en una cantidad 8, como la
representada en la figura 4, en la que podemos ver un triángulo geodésico ABC y las secciones
normales recíprocas, correspondientes a las observaciones realizadas en sus vértices.
Para corregir estas medidas angulares se debe obtener, en primer lugar, el ángulo A que
forman dos secciones normales recíprocas en un vértice, por ejemplo, el ángulo que
15
forman las curvas AaB y BbA, en el punto A o en el B. En segundo lugar, debemos relacionar
este ángulo con la pequeña cantidad ¿que constituye la corrección.
Para hallar el ángulo A tenemos en cuenta que éste será siempre una cantidad muy
pequeña, entonces, obteniendo su tangente mediante la expresión (Cid y Ferrer, 1997)
'l+. s2
KT + s|-K'T-KT'
tgA 2
donde r denota en este caso la torsión geodésica (que es distinta de cero). De esta fórmula y
tras despreciar términos de grado s3 y superiores, podemos deducir (Cid y Ferrer, 1997)
2 2
A= 4a cos 24)sen 2A
donde <f> es la latitud del vértice para el que se calcula A, e es la excentricidad del elipsoide de
referencia, a es su semieje mayor y A es el acimut de la sección normal directa.
donde k es un factor de escala cuyo valor suele ser del orden de 10- 6 y R es la matriz de
rotación que define la transformación de coordenadas que relaciona ambos sistemas, dada por
(Torge, 1991)
(8.1) + (1 + k)
J
Sy -Sx
\zJ vzoy
vz J
donde los ángulos (sx, Sy, Sz) expresan las rotaciones indicadas en la figura 5(a). Estos ángulos
suelen tener un valor muy pequeño, por eso se expresan habitualmente en segundos de arco.
Es importante notar que la ecuación (8.1) expresa una relación entre coordenadas cartesianas.
No obstante, en muchas ocasiones no trabajamos con las coordenadas cartesianas (x, y, z) de
un punto P, sino con sus coordenadas geodésicas (<|>, X, h), donde § es la latitud geodésica
(medida en grados norte), X es la longitud geodésica (medida en grados este) y h es la altura
elipsoidal (medida en metros) ilustrada en la figura 5(b).
16
(a) (b)
Fig. 5. (a) Representación del sistema de referencia no geocéntrico xyz, con respecto al
sistema de referencia geocéntrico XYZ. (b) Posición de un punto P medida sobre un sistema
de coordenadas cartesiano y sobre el elipsoide de referencia.
Entonces, nos damos cuenta de que es necesario poder convertir las coordenadas
geodésicas (<|>, X, h) en coordenadas cartesianas (x, y, z) y viceversa. Para ello, podemos
recurrir de nuevo a una suma vectorial, como la representada en la figura 5(b)
r = rQ + hN
1
7
Llegados a este punto, tenemos ya las herramientas necesarias para abordar el
problema de transformación de coordenadas planteado al principio de este apartado, así como,
algunos otros problemas implicados en el mismo. Entonces, dependiendo de los problemas de
transformación de coordenadas que tengamos planteados, podemos optar por aplicar una u
otra fórmula, de las anteriormente obtenidas. Concretamente, los problemas que podemos
resolver con las fórmulas (8.1), (8.2) o sus inversas; son los siguientes:
17. Transformación de coordenadas geodésicas (é, A h) en coordenadas cartesianas (x y,
z), dadas en el mismo dátum. Para resolver este problema recurrimos a aplicar directamente
la fórmula (8.2).
18. Transformación de coordenadas cartesianas (x y, z) en coordenadas geodésicas (é, A
h), dadas en el mismo dátum. Para resolver este problema aplicamos la fórmula inversa de
(8.2), que puede ser computada fácilmente de forma analítica y numérica.
19. Transformación de coordenadas cartesianas (x’, y’, z’) dadas en el dátum1 a
coordenadas cartesianas (x’, y’, z’) dadas en el dátum2. Para resolver este problema
operamos en distintas etapas:
a) Aplicamos la fórmula (8.1) para convertir coordenadas cartesianas (x’, y’, z’), dadas en
el dátum 1, en coordenadas cartesianas geocéntricas (x, y, z). Para ello, necesitamos
establecer los valores de los parámetros (k, xo, yo, zo, sx, Sy, Sz), correspondientes al
datum1. Estos valores pueden obtenerse desde varias fuentes. Por ejemplo, Torge (1991) ha
establecido los valores de estos parámetros, para transformar coordenadas desde algunos
de los dátum más utilizados al sistema geocéntrico WGS84 (el que se utiliza con GPS).
b) Aplicamos la inversa de la fórmula (8.1) para convertir las coordenadas cartesianas
geocéntricas (x, y, z), obtenidas en el paso anterior, en las coordenadas cartesianas (x’, y’,
z’) dadas en el dátum2. Para ello, necesitamos establecer los valores de los parámetros (k,
xo, yo, zo, sx, Sy, Sz), correspondientes ahora al datum2. La inversa de la fórmula (8.1) puede
realizarse fácilmente, programando las ecuaciones (8.1) de forma matricial en un
ordenador, en ese caso todo el problema se reduce a invertir una matriz. Este problema de
invertir una matriz es muy conocido y es fácil de resolver.
4) Transformación de coordenadas geodésicas (é, A h) dadas en el dátum1 a coordenadas
geodésicas (é, A h) dadas en el dátum2. Para resolver este problema podemos dividirlo en
tres problemas más simples que ya hemos resuelto antes:
20.Transformar primero las coordenadas geodésicas (<|>, X, h) en coordenadas cartesianas
(x’, y’, z’), para el dátum1 en el que están dadas. Para resolver este problema aplicamos
directamente la fórmula (8.2), tal como en el problema (1) anteriormente resuelto.
21.Transformarlas coordenadas cartesianas (x’, y’, z’), obtenidas en el paso anterior, a
coordenadas cartesianas (x’, y’, z’) dadas en el dátum2. Para resolver este problema
operamos como se ha indicado antes en el problema (3).
22.Finalmente, transformamos las coordenadas cartesianas (x’, y’, z’), obtenidas en el paso
anterior, en coordenadas geodésicas (<|>, X, h). Para resolver este problema aplicamos la
fórmula inversa de (8.2), tal como se ha indicado en el problema (2) anteriormente
resuelto.
18
Aplicaciones de la Geodesia Física
En este sentido, el importante papel que juega el campo de gravedad terrestre, queda
todavía más claro, si tenemos en cuenta que cuando el topógrafo representa la altura de un
punto de la superficie de la Tierra, no son las alturas referidas al elipsoide de referencia las
que tienen sentido, sino las alturas referidas al geoide (alturas ortométricas H). Así, el
conocimiento de las alturas de los puntos de la superficie terrestre, determina el relieve que el
topógrafo debe representar en los mapas topográficos. En consecuencia, altura ortométrica H
de los puntos de la superficie terrestre, es una de las coordenadas que determinan la figura de
la Tierra, junto con las coordenadas geodésicas latitud y longitud, siendo esta altura H una
cantidad que requiere el conocimiento del campo de gravedad terrestre para su determinación.
19
10. El campo de gravedad terrestre. Superficies de nivel. El geoide
El geoide como superficie de referencia fundamental, es una superficie cuya forma
refleja la distribución de masas en la Tierra. Por ello, debemos empezar recordando cuál es la
fuerza gravitatoria producida por una cierta distribución de masa y su potencial asociado, para
entender cómo surge de estos conceptos y del concepto de superficie equipotencial, una
definición de superficie de referencia, tal como es la definición de geoide.
Así, debemos comenzar recordando que cuando consideramos dos masas puntuales,
separadas una distancia s, la fuerza con la que cada una atrae a la otra es
km1m2
(Ley de gravitación de Newton)
(b)
(a)
Fig. 1. (a) Fuerza de atracción gravitatoria debida a una masa puntual. (b) Fuerza centrífuga
debida a la rotación de la Tierra en punto P de la misma.
20
También debemos recordar que el campo gravitatorio es un campo conservativo, por
ello puede obtenerse a través del gradiente de un campo escalar, es decir
grad V = —,—,— = F
5>Zkm V- kmi
obviamente, para una distribución continua de masa M, el sumatorio debe ser reemplazado
por la integral, escribiendo
donde p denota la densidad de masa y V denota el volumen del cuerpo considerado. En el caso
de que este cuerpo sea la Tierra, debemos tener en cuenta que la misma no está en reposo, sino
rotando con una velocidad angular de 0.7297x10-4 seg-1. En este caso para poder considerar que
la Tierra está en reposo y aplicar la ecuación anterior, debemos incluir una fuerza de inercia
conocida como fuerza centrífuga, de tal forma que para cualquier punto P sobre la superficie
de la Tierra debemos considerar una fuerza f, como la representada en la figura 1(b). Dicha
fuerza puede escribirse como
f =ca2 p = (ca2 x,ca2 y, 0) = grad^=í^,^,^l ,, * = 12 (x 2 + y 2 )
\dx dy dz) 2
donde/; es un vector distancia, paralelo al plano ecuatorial, que va desde el eje de rotación de
la Tierra hasta el punto P. Esto nos lleva a una definición de potencial gravitatorio del tipo
grad W =,,— = g
l dx dy dz)
El módulo de este vector, o magnitud de la gravedad, se suele expresar en gales siendo 1 gal =
1 cm/seg2. La magnitud de este vector se incrementa sobre la superficie de la Tierra cuando
vamos desde el ecuador hacia los polos, debido a que disminuye la fuerza centrífuga.
Concretamente, la gravedad en el ecuador es 978 gal y en los polos 983 gal,
aproximadamente.
21
En lo referente a W, el potencial de la gravedad, hay que recordar que la superficie
para la cual éste es constante, es decir, W(x, y, z) = W0, se denomina superficie equipotencial
o superficie de nivel, siendo el geoide la superficie de equipotencial del campo de gravedad
terrestre que coincide con el nivel medio de los océanos. La figura 2(a) nos muestra la
relación que existe entre las superficies de nivel y la línea de la plomada, notando que estas
superficies son normales a la línea de la plomada.
(a)
Fig. 2. (a) Superficies equipotenciales del campo de gravedad terrestre y línea de la plomada
en un punto P. (b) Gravedad sobre una superficie equipotencial.
También hay que notar que esta definición de altura, en la que la superficie de
referencia es el geoide, nos lleva a establecer un sistema de alturas en el que obtenemos la
altura de un punto P en la forma
22
C dC
H=
W 0 g
donde
Wg C=W
- W0 es
=J denomi
nado
númer
o
geopot
encial.
L
le
ga
d
os
a
es
te
p
u
nt
o
ca
be
pr
eg
u
nt
ar
se
q
ué
re
la
ci
ó
n
ex
ist
e
en
tr
e
es
te
si
st
e
m
a
de esta h = H
al cuestió (b)
tu n
ra puede Fig. 3.
s verse (a)
y bien Altura
la ilustrad ortomé
s a en la trica y
al figura elipsoi
tu 3(a), en dal en
ra la que un
s represe punto
de ntamos P. (b)
te para un Elipsoi
r punto de de
m P, sus referen
in alturas cia.
ad asociad
as as, E
so tanto la n
br elipsoi consec
e dal h uencia,
el como determ
el la inada h
ip ortomé o H,
so trica H, obtener
id existie la otra
e ndo la altura
de sencilla es
re relació sencill
fe n o si se
entr
re e conoce
nc amb N, la
ia as ondula
. altur ción
L as. del
a geoide.
re Éste es
sp uno de
ue los
st proble
a mas
a princip
ales de más cercano mi
la a la figura
geodes de la Tierra, ca
ia tomando por de
práctic ello un
a o elipsoide de un
aplicad revolución pu
a, la ligeramente
nt
determ aplanado
inación por los o
de las polos. Este P.
alturas elipsoide
del vendrá (b)
geoide definido De
sobre entonces por
un dos svi
elipsoi parámetros, aci
de de su semieje
ón
referen mayor a
cia (aproximada de
dado. mente el la
Para radio
resolve ecuatorial ve
r este terrestre) y
proble su rti
ma aplanamient cal
hemos o f definido
de por en
fijar, f = a el
-b
en
primer a pu
lugar, ( nt
un o
elipsoi b) Fig.
de de 4. (a) P.
referen
Latitud
cia que
23
sea lo astronó
donde b es el semieje menor. La figura 3(b) muestra las coordenadas elipsoidales de un punto
P, sobre el elipsoide de referencia elegido. Estas coordenadas (<|>, X, h) no son obviamente las
coordenadas que determinamos en una campaña geodésica. En una campaña geodésica
realizada con instrumentos tradicionales (es decir, sin GPS) determinamos las coordenadas
naturales (O, A, H). Esto es debido a que cuando situamos un instrumento sobre el terreno, la
vertical es para nosotros la dirección de la línea de la plomada, tal como puede verse en la
figura 4(a), siendo la latitud y longitud del lugar, la latitud y longitud astronómicas. La
separación que existe entre la línea de la plomada y la normal al elipsoide de referencia en el
punto P, punto en el que situamos el instrumento, está representada en la figura 4(b). A partir
de esta figura pueden obtenerse las relaciones
£ = <D-<|>
TI = (R cos <|> )(A -X) = (A-X) cos §
donde R es el radio de la esfera situada en el punto P, cuyo valor es la unidad. Las cantidades
(£,, TI) reciben el nombre de componentes de la desviación de la vertical. Así, las medidas
realizadas sobre el terreno en un vértice de una red geodésica (OP, AP, HP), deben ser reducidas
a las cantidades (&P, hP, hP), medidas sobre el elipsoide de referencia utilizado en la campaña de
triangulación, mediante las fórmulas
<|> = <!>-£ „ X = A- cos „ h=H+N
donde si son conocidas la ondulación del geoide y las componentes de la desviación de la
vertical, será fácil pasar del sistema de referencia físico al sistema de referencia matemático o
viceversa.
24
N= (Fórmula de Bruns) ,, y =
y grad U
donde y es denominada la
gravedad normal. Este vector
tendrá la dirección de la
normal al elipsoide de
referencia, tal como se indica
en las figuras 3(b) y 4(b),
siendo su diferencia con el
vector gravedad denominada
anomalía de la gravedad Ag.
Esta última cantidad definida
nos permite obtener el
potencial perturbador y como
consecuencia la ondulación
del geoide, si llevamos a cabo
un desarrollo en armónicos en
esféricos de este potencial
perturbador T (Heiskanen y
Moritz, 1985), teniendo como
resultado final las importantes
relaciones
íí
sen (v|//2)
R
Ag SO) da (Fórmula de Stokes)
R T Ag S(»da ,, N = se
1
= 4n
2 TTG - 6 sen V
(v|//2) + 1 - 5
2)]
::
\\iln[ sen
(
(v|//2) +
v
|/
)
donde ores er la fórmula de Stokes. (b)
R esexten Relación entre las
el didas coordenadas polares (y, a) y las
radio a toda coordenadas geográficas (<|>,
medio la X), en la esfera de radio
terrest esfera unidad.
re unida
(6371 d, tal La fórmula de Stokes
km), como publicada por él en 1849, es la
G esse expresión matemática más
el indica importante de la Geodesia
valor en laFísica, debido a que hace
medio figura posible conocer la figura de la
de la5(a). Tierra (es decir, la ondulación
grave En del geoide) a partir de datos
dad esta gravimétricos. Debemos notar
sobre figura también que la ecuación
la y enintegral de Stokes es el
Tierra la vínculo fundamental entre el
(979.8figura campo de gravedad terrestre y
gales) 5(b), la figura de la Tierra,
y a espode confirmando la afirmación
el mos realizada al principio de este
capítulo, cuando se decía que
área ver
para avanzar hacia una figura
de laque el
de la Tierra más precisa que
esfera papel
un elipsoide de revolución, era
(
25
determinadas por toda la Tierra, entre otros datos) podemos obtener la ondulación del geoide,
tal como se indica en las figuras 6(a) y 6(b).
(a) (b)
Fig. 6. (a) Mapa de la ondulación del geoide sobre un elipsoide de revolución. (b) Imagen del
geoide en 3D mostrando los detalles de la forma del geoide, a una escala en la que se han
exagerado los valores de la ondulación para ver mejor dichos detalles.
23.La fórmula presentada para la ondulación del geoide ha sido obtenida considerando
que el elipsoide de referencia (superficie U(x, y, z) = U0), tiene el mismo potencial que el
geoide (U0 = W0). Además, esta superficie debe encerrar en su interior una masa que sea
numéricamente igual a la de toda la Tierra.
25. El potencial anómalo o potencial perturbador debe ser armónico fuera del geoide, es
decir, T(x, y, z) puede ser desarrollado en armónicos esféricos. Esto significa que el efecto
de las masas por encima del geoide, debe ser eliminado aplicando adecuadas reducciones
a los valores de la gravedad determinados. De lo contrario, no será posible aplicar la
expresión integral de Stokes.
26
de anomalías de la gravedad distribuidas por toda la Tierra. Estas fórmulas son debidas a
Vening Meinesz y fueron obtenidas por este geodesta en 1928, su expresión es (Heiskanen y
Moritz, 1985)
Ll 1
4*G
2n
a=0 v=0
n dS( ) rcosa]
Ag(\|/,a)^N
dW [sena
senydyda
donde a y y son las coordenadas polares sobre la esfera de radio unidad, tal como se indica en
la figura 5(b), siendo da = sen y d\|/ da y S(\|/) la función de Stokes, cuya expresión ha sido
dada anteriormente.
En efecto, puede hallarse con precisión la altura ortométrica H de un punto sin tener
conocimiento del geoide, o al menos sin un conocimiento detallado o muy preciso del mismo.
El procedimiento que se sigue para calcular las alturas ortométricas directamente sobre el
terreno se llama nivelación de precisión. A continuación vamos a describir este proceso con
cierto detalle.
(b)
Fig. 7.
(a)
Esque
ma del
proced
imient
o de
nivelac
ión
geomé
trica.
(b)
Difere
ncias
de
alturas
ortomé
tricas
debida
s al no
paralel
ismo
de las
superfi
cies de
nivel.
P
ara
entend
er lo
que
sucede
basta
consid
erar
dos plomad
puntos a que
A y B, pasa
suficie por B.
ntemen Ello es
te debido
alejado al no
s, de paraleli
tal smo de
manera las
que el superfi
proces cies de
o nivel,
mostra tal
do en como
la se ve
figura en la
7(a), figura,
debe en
ser consec
aplicad uencia,
o si
repetid design
as amos
veces. por 5W
Esta el
situaci increm
ón es ento
la que del
está potenci
represe al de la
ntada graved
en la ad,
figura tenemo
7(b), s
en esta
-8W =
figura
g5n =
la
g’5HB
diferen
cia de donde
altitud g es la
nivelad graved
a 8n, es ad en
distinta la
al estació
increm n de
ento nivelac
5HB en ión y
la línea g’ es la
de la graved
ad la
sobre combi
la línea nación
de la de
ploma medida
da de s
B en geomét
5HB. ricas
con la
La
graved
relació
ad, nos
n propor
anterio ciona
r la
implic diferen
a que cia de
5HB = potenci
g 5n * al 5W,
5
n de tal
forma
no que
habien B WB - WA =
do -^ g5n
relació A
n A
geomét sí, la
rica nivelac
directa ión
entre el combi
resulta nada
do de con
la medida
nivelac s de la
ión graved
geomét ad nos
rica y propor
la ciona
diferen diferen
cia del cias de
alturas potenci
ortomé al, esto
trica. es,
Esta cantida
consec des
uencia físicas
es no
inmedi geomét
ata si ricas
notam con
os que una
ventaja = A
añadid 0 g
a, la = Ad
diferen n
cia de donde
potenci hemos
al entre cambia
dos do el
puntos sumato
es rio por
indepe la
ndiente integra
del l por
camino ser
que se teórica
siga, mente
obtenié más
ndose riguros
el o. Este
mismo resulta
resulta do es
do con imposi
indepe ble, en
ndenci general
a de , si
los consid
itinerar eramos
ios que sólo
se medida
sigan, s
por geomét
ello, si ricas,
medim ya que,
os a lo la
largo suma
de un de
circuit increm
o de entos
nivelac de
ión nivelac
(segui ión
mos un depend
trayect e de la
o trayect
cerrad oria,
o) esto es
sucede
que
W -W
28
A A
AnAB = B5n ,,
A A dn = error de cierre
Fig. 8. Altura ortométrica de un punto P
sobre la superficie de la Tierra.
Como
consecuencia de
todo lo dicho, las
diferencias de
potencial son
básicas en toda
teoría de altitudes,
concluyendo que la
nivelación sin
medidas de
gravedad, que en la
práctica puede
aplicarse a
pequeñas distancias
(para las que las
superficies de nivel
son prácticamente
paralelas), no tiene
significado desde
un punto de vista
riguroso, llevando a
contradicciones
como es el error de
cierre. Todo esto
hace que se
considere, en
trabajos geodésicos
de gran precisión
como es la
nivelación de
precisión, la
diferencia de
altura ortométrica
δHAB en lugar de la
diferencia de altura
geométrica ∆nAB,
ya que, tal como
puede verse en la
figura 7(b), las
alturas de los
puntos A y B sobre el
geoide HA y HB, son H
siempre las mismas g= g dH
con independencia 0
del trayecto de es el valor medio
nivelación que se siga de la gravedad
entre ambos, debido a sobre la línea de la
que la definición de plomada entre P y
altura ortométrica P0. La expresión
es consecuencia de anterior nos
principios físicos, no permite obtener la
una cantidad obtenida altura ortométrica
de un sumatorio. del punto P en la
forma
Quedando
claro qué sistema de g
alturas debe usarse H=
para obtener una gran
precisión, el esta ecuación no es
problema planteado una tautología, es
es la realización de realmente una
esta medida de altura, fórmula de utilidad
es decir, el práctica, debido a
procedimiento que la gravedad
operativo que debe
media g no
seguirse para calcular
depende
la altura ortométrica
de un punto P. Para fuertemente con la
ello, vamos a altura ortométrica.
considerar un punto P En consecuencia,
sobre la superficie de sin un
la Tierra, tal como se conocimiento
indica en la figura 8, preciso de la altura
designando por P0 la del punto P sobre el
intersección con el geoide, es decir,
geoide de la línea de con un valor inicial
la plomada que pasa de H
por P. Entonces, suficientemente
recordando la cercano a la
definición de número
verdadera altura
geopotencial,
ortométrica,
podemos escribir
podemos abordar el
cálculo de la
gravedad g . En
d H H
W=
0 g
I este sentido,
H 0 podemos encontrar
expresiones
- d
WH aproximadas de g
en diversos libros
== de geodesia
H (Heiskanen y
C 0 Moritz, 1985).
Por otra parte,
si se quiere aplicar la
sencilla fórmula
anterior es necesario
hallar, además de g ,
el número
geopotencial de P
designado por C. Para
ello, un
procedimiento que
29
puede seguirse es elegir un punto O al nivel del mar (es decir, sobre el geoide), eligiendo un
punto apropiado sobre la costa. Entonces, puede determinarse el número geopotencial de P
aplicando la integral
C = W0 – W =∫ 0
P
g dn
siendo como ya sabemos independiente del itinerario de nivelación utilizado para relacionar
los puntos P y O. El número geopotencial se mide en unidades geopotenciales (u. g. p.). Estas
unidades están derivadas de las unidades de la gravedad y de la altura, mediante la relación
1 ugp = 1 kgal.m = 1000 gal.m
como g es aproximadamente 0.98 kgal, tenemos que C es aproximadamente igual a 0.98H.
Esta técnica es muy empleada en la actualidad para obtener las alturas del geoide sobre
el elipsoide de referencia, en trabajos geodésicos de gran precisión, basta mencionar como
ejemplo el proyecto ALGESTAR (Bürki and Marti, 1990), cuyo objetivo fue llevar a cabo la
computación de un nuevo geoide para Suiza, comprobando al mismo tiempo la fiabilidad de
las medidas GPS en una zona muy montañosa. Para ello, se obtuvieron las ondulaciones del
geoide aplicando la ecuación
N=h-H
donde h fue obtenida a partir de observaciones GPS y H a partir de la computación de los
números geopotenciales (nivelación de precisión), siguiendo las líneas de nivelación de la Red
Nacional de Nivelación que une los 40 vértices geodésicos de la Red Suiza de Triangulación.
El resultado de este proyecto fue la obtención de las ondulaciones geoide para Suiza con una
gran precisión.
30
la latitud y longitud del punto considerado. En este aspecto se centra una de las principales
aplicaciones de la geodesia física, puesto que permite obtener la altura del geoide para
cualquier punto de la Tierra. Estas alturas del geoide se refieren, generalmente, al Sistema
Geodésico de Referencia de 1980 (GRS80).
Los modelos geopotenciales pueden servir para calcular la ondulación del geoide. No
obstante, las alturas del geoide provenientes de un modelo geopotencial no son, en general,
suficientemente precisas. Por ello, es necesario determinar una corrección a las mismas. Esta
corrección se determina a partir de las anomalías de la gravedad medidas sobre el terreno,
para computar así un modelo de geoide más preciso. Debido a que las anomalías observadas
sobre el terreno, están dadas en un área limitada y no pueden por tanto utilizarse para resolver
la longitud de onda larga del campo de gravedad terrestre, es necesario considerar un modelo
geopotencial para la computación de tales longitudes de onda. Como veremos a continuación,
la contribución de esta longitud de onda larga corresponde a una aproximación suave del
geoide, para la región bajo estudio.
El cálculo gravimétrico de una altura del geoide muy precisa, para una región dada,
tendrá una gran utilidad práctica directa e inmediata para toda la comunidad de usuarios de
GPS (Kaula, 1987). Estos resultados permitirán incorporar modelos más realísticos a los
algoritmos de cálculo de las alturas ortométricas, reportando una mejora en la precisión de las
alturas mencionadas, lo que conllevará evidentes beneficios para la geo-referenciación, la
fotogrametría o el control geométrico de obras. En suma, serán beneficiarios de tales
resultados aquellos científicos y técnicos, interesados en una localización más precisa sobre el
terreno. Puesto que, usando medidas de posición con el sistema GPS obtendrán a la vez una
precisión suficiente en planimetría, junto con una altura ortométrica precisa, hallada a partir
de un modelo de geoide (cuya ondulación sea precisa) y la altura elipsoidal dada por GPS,
pudiendo realizar medidas precisas de alturas ortométricas sin nivelación (Schwarz et al.,
1987).
Fotogrametría. Tener un modelo de geoide, supone poder tomar puntos de apoyo sin
necesidad de apoyarse en las redes geodésicas, que es lo que se hace actualmente para
calibrar la altura elipsoidal (la obtenida con GPS), a altura ortométrica en la zona de
trabajo (la que se necesita para la orientación absoluta del modelo estereoscópico).
Para ello, suponemos que elipsoide y geoide son paralelos en esa zona, cosa que no es
rigurosamente cierta. La limitación del trabajo de apoyo en campo en Fotogrametría,
ya sólo estaría limitada por la capacidad de los equipos empleados y se elevaría el
rendimiento de éstos, puesto que no habría que inmovilizar equipos GPS para las
calibraciones. Además de lo mencionado anteriormente, si se dispusiese de un marco
de referencia común entre los puntos sobre el terreno y el centro de proyección de las
cámaras aéreas, se simplificaría todavía más el trabajo de apoyo de campo.
31
Caracterización agronómica y medioambiental. La posibilidad de determinar la altitud
sobre el nivel del mar de forma rápida y eficaz, permitiría añadir la altitud como
variable de estudio en los SIG (Sistemas de Información Geográfica), para su inclusión
en los posteriores análisis con las demás variables geográficas.
Vista la importancia que tiene usar un modelo de geoide preciso, vamos a ver aquí cuál
es el procedimiento operacional para obtenerlo. Para comenzar, hay que decir que las técnicas
que vamos a describir a continuación, se han mostrado herramientas eficaces a la hora de
indagar la ondulación del geoide, en distintas regiones de la Tierra (Kearsley, 1988; Forsberg
and Kearsley, 1989; Torge et al., 1989; Yun, 1994). Estas técnicas pretenden evitar el
problema que conlleva la obtención de anomalías de la gravedad en un área limitada, puesto
que con estos datos no podemos resolver las largas longitudes de onda del campo gravitatorio
terrestre. Por ello, necesitamos considerar un modelo geopotencial con el que podemos
obtener estas longitudes de onda larga. Así, en un paso previo de computación, obtendremos
las anomalías aire-libre referidas al modelo geopotencial, computando la diferencia con las
anomalías predichas por este modelo, estas diferencias serán usadas en un procedimiento de
integración con la FFT (Schwarz et al., 1990), para obtener finalmente la ondulación del
geoide, la cual en un proceso de computación posterior será añadida a la ondulación del
geoide predicha por el modelo geopotencial considerado. En suma, el proceso a seguir consta
de las etapas que a continuación serán descritas (Corchete et al., 2005).
Interpolación de los datos de gravedad. En primer lugar, hay que saber que si queremos
aplicar técnicas basadas en la Transformada de Fourier Rápida (FFT), es necesario pasar de la
distribución irregular de puntos, dada por la localización geográfica de nuestras medidas de
anomalías aire-libre de la gravedad, a una distribución regular de puntos en forma de una
rejilla de puntos equidistantes. Como ejemplo, podemos ver en la figura 9(a) los datos de
gravedad disponibles para calcular un modelo de geoide para el área ibérica. A la vista de esta
figura notamos que estos puntos no están distribuidos de forma regular, hay zonas con muchos
puntos y otras áreas que tienen muy pocos. Con esta distribución de valores no podemos
aplicar las técnicas basadas en la FFT, pues estas técnicas requieren una distribución regular de
puntos en forma de una rejilla de puntos equidistantes, como muestra la figura 9(b). Por otra
parte, antes del proceso de interpolación es necesario reducir los valores de estas anomalías de
la gravedad, con el objeto de trabajar con una colección de puntos cuyos valores sean lo más
parecidos entre ellos que sea posible, para facilitar el trabajo de la interpolación y obtener
entonces tras este proceso la mejor rejilla de valores posible. Para ello, antes de la
interpolación es muy conveniente eliminar los efectos de corta longitud de onda debidos a la
topografía y la batimetría. Esta corrección se llama corrección RTM (Forsberg and
Tscherning, 1997) y nos permite obtener un campo de anomalías aire-libre suavizado, muy
adecuado para la interpolación a una rejilla regular. Esta corrección suele ser aplicada junto
con una corrección de onda larga, obtenida restando las anomalías de la gravedad predichas
por un modelo geopotencial, a los valores de las anomalías aire-libre observadas. En suma, la
corrección total que debemos aplicar a los valores de anomalías aire-libre mostrados en la
figura 9(a) está dada por
32
Agpretds = Agfprtese-27ikp(h-href)pts + cpts-AgGptMs (13.1)
donde el superíndice pts denota cada punto distribuido sobre el área de estudio de forma
aleatoria, Agfree es la anomalía aire-libre determinada en ese punto (figura 9(a)), k es la
constante de gravitación de Newton, p es la densidad de la topografía (2.67 g/cm3) para la
corrección RTM en tierra o la densidad de la topografía menos la densidad del agua marina
(2.67 - 1.03 = 1.64 g/cm3) para la corrección RTM marina, h es la elevación, href denota la
elevación de la superficie de referencia, c es la corrección del terreno computada para puntos
terrestres y marinos y AgGM es la anomalía de la gravedad computada desde un modelo
geopotencial (esta computación será luego descrita).
(a)
(b)
Fig. 9. (a) Distribución geográfica de los valores de anomalías aire-libre de la gravedad,
medidos para el área ibérica. (b) Mapa de las anomalías de la gravedad interpoladas tras ser
reducidas mediante la fórmula (13.1).
33
(a) (b)
Fig. 10. (a) Modelo de elevaciones del terreno proporcionado por el modelo digital del terreno
SRTM, completado con la batimetría dada por ETOPO2. (b) Superficie de referencia obtenida
mediante el filtrado 2D del campo de elevaciones representado en la figura 10(a).
La corrección del terreno puede ser positiva o negativa cuando es computada en áreas
marinas, dependiendo de la forma de la superficie del fondo marino (Torge, 1989). Esto no es
posible en áreas continentales en las que la corrección del terreno es siempre positiva.
34
siguiendo el segundo método de condensación de Helmert. A continuación se indicará cómo son
computados los tres términos de la fórmula (13.3).
Contribución de onda larga (término N1). Está contribución se obtiene mediante un modelo
geopotencial a través del desarrollo en serie de armónicos esféricos dado por (Heiskanen y
Moritz, 1985)
kry
M / / Pnm(sin(p) [5Jnmcosm^ + 5Knmsin m^] + N0 (13.4)
Á—¡\ r I <¿—/
N1 m=0
donde (5Jnm, SKnm) son los coeficientes geopotenciales totalmente normalizados del potencial
anómalo de a gravedad, a es el semieje mayor del elipsoide de referencia, kM es la constante
gravitacional geocéntrica de la Tierra, Pnm son las funciones de Legendre completamente
normalizadas, (r , cp, X) son la posición geocéntrica de los puntos sobre la superficie del
geoide, y es la gravedad normal sobre el elipsoide de referencia, nmax es el grado máximo de el
modelo geopotencial considerado y N0 el término de grado cero del desarrollo.
(a) (b)
Fig. 13. (a) Ondulación del geoide obtenida a través de un modelo geopotencial para el área
de estudio. Intervalo entre isolíneas: 1m. (b) Anomalía de la gravedad obtenida a través de un
modelo geopotencial para el área de estudio.
35
En la computación de la ecuación (13.4), es necesario obtener la posición geocéntrica
r del punto considerado sobre el geoide, para ello, podemos seguir el procedimiento en dos
pasos indicado por Kuroishi (1995). Según este procedimiento, calculamos en primer lugar el
radio geocéntrico elipsoidal rE asociado a las coordenadas (cp, X) del punto considerado. En un
segundo paso y por medio de la ecuación (13.4), poniendo r = rE obtenemos la anomalía en
altura £E, la cual puede ser usada como una primera aproximación de la ondulación del geoide
en el punto considerado. Finalmente, obtendremos el valor de N(cp, X) en la posición
geocéntrica del punto considerado r, empleando la fórmula (13.4) con el valor r = rE + £E. Hay
que notar que según este método de computación, r es considerado como la posición
geocéntrica de un punto P sobre el geoide, en lugar de considerar la posición geocéntrica del
punto P sobre la superficie de la Tierra. Por esta razón, siguiendo este método de computación
en dos pasos, lo que obtenemos son ondulaciones del geoide en lugar de anomalías en altura.
En la figura 13(a) podemos ver los resultados obtenidos mediante la aplicación de la fórmula
(13.4), a los puntos de la rejilla regular ubicada en el área de estudio.
n=2 m=0
donde Ag0 es el término de grado cero del desarrollo. En la computación de la formula (13.5), hemos
usados el valor de la posición geocéntrica r obtenida previamente para la fórmula (13.4). La figura
13(b) muestra los resultados obtenidos mediante la aplicación de la fórmula (13.5), a los puntos de
la rejilla regular ubicada en el área de estudio.
Efecto indirecto (término N2). El efecto indirecto en el geoide debido a la aplicación del segundo
método de condensación de Helmert se puede calcular en aproximación plana mediante los dos
términos (Sideris, 1990)
N2=-^h 2 ( x , y ) - kP ff h (x,y)~h (xp,yp ) dxdy (13.6)
3/2
y 6y L _x)2+( _ )2
pp 6y JJ
ríf
f(x,y)= (x2+y2)3/2 g= f(x,y)dxdy
A
36
Así, mediante la aplicación de la formula de convolución (13.7), podemos calcular el
efecto indirecto en el geoide, debido a la aplicación del segundo método de condensación de
Helmert, obteniendo para el área de estudio considerada aquí el resultado ilustrado en la
figura 14(a).
(a) (b)
Fig. 14. (a) Efecto indirecto en la ondulación del geoide producido en la eliminación de las
masas exteriores al geoide, según el segundo método de condensación de Helmert. (b)
Corrección del terreno asociada al segundo método de condensación de Helmert.
37
(a) (b)
Fig. 15. (a) Anomalías residuales obtenidas mediante la aplicación del segundo método de
condensación de Helmert, fórmula (13.9), a las anomalías calculadas mediante la fórmula
(13.2). (b) Ondulación del geoide obtenida integrando las anomalías residuales dadas por
(13.9), mediante la fórmula (13.10).
Cuando los términos de la formula (13.9) están todos calculados, estamos en situación
de realizar la integral (13.8). Para ello, podemos expresar esta integración en forma de
convolución usando la FFT en 1D (Haagmans et al., 1993), obteniendo el término N 3 en la
forma
RA4)A?
n F1(Agcos) F1(S)
iF
N13
4TTY (13.10)
>'H1
donde S es la función de Stokes, F1 denota la transformada de Fourier discreta y periódica y
F1-1 es su inversa. Haagmans et al. (1993) han descrito la forma correcta de llevar a cabo la
computación dada por (13.10), por ello, siguiendo estas indicaciones podemos obtener el
valor de la ondulación residual del geoide mostrada en la figura 15(b).
Validación del geoide. Cuando sumamos los tres términos dados por la formula (13.3),
obtenemos el modelo de geoide buscado, en este caso es el geoide del área ibérica mostrado
en la figura 16(a). Para chequear la validez del modelo obtenido, debemos comparar las
ondulaciones del geoide predichas por este modelo, con las ondulaciones del geoide obtenidas
de forma independiente a este modelo. Estas medidas independientes deben haber sido
calculadas con gran precisión, para poder establecer entonces que las diferencias, entre los
valores predichos por el modelo y los valores observados o medidos, son debidas únicas y
exclusivamente al error que tiene el modelo de geoide obtenido. Para ello, lógicamente
necesitamos utilizar una colección de medidas de gran precisión, pues de lo contrario todo el
proceso de validación sería dudoso. En el área ibérica disponemos de tal colección de medidas
gracias a la red GPS-nivelación que está desplegada para toda Europa, tal como se explicó en
el apartado 4. Esta red cuyo nombre es EVRS (European Vertical Reference System),
constituye un marco de referencia fundamental para la medición de alturas con precisión
centimétrica. Concretamente, en el área ibérica esta red posee los vértices de precisión
mostrados en la figura 16(b). Cuando comprobamos el modelo obtenido y representado en la
figura 16(a), en los puntos de precisión mostrados en la figura 16(b), obtenemos un error de
27.8 cm (Corchete et al., 2005). Este modelo no ha sido mejorado hasta 2008, fecha en la que
han sido publicados los modelos del norte (Corchete, 2008) y sur de Iberia (Corchete et al.,
2008), cuyos errores son de 12.6 y 17.0 cm, respectivamente.
38
(a) (b)
Fig. 16. (a) Modelo de geoide obtenido para el área ibérica, mediante la suma de los tres
términos dados por (13.3). (b) Red de puntos para los cuales es conocida la ondulación del
geoide con gran precisión.
39
nmáx n
esto implica que los valores de (£, T)) dados por (15.1) deberán ser corregidos mediante
terreno = geoide " 50 „ Tferreno = Tgeoide – 5Acos <|>
-2
(b) Fig. 17. (a) Componente N-S de la
desviación de la vertical (£). (b) Componente E-W de la
desviación de la vertical (TI). Ambas componentes han sido obtenidas mediante (15.2) a partir
del geoide NIBGEO.
41
Bibliografía
Bibliografía básica
Bomford G. (1971). Geodesy. Oxford University Press.
Bürki B. and Marti U. (1990). The Swiss Geoid Computation: A Status Report. International
Association of Geodesy Simposia. Symposia nº 106: Determination of the Geoid.
Springer Verlag.
Cid R. y Ferrer S. (1997). Geodesia Geométrica, Física y por Satélites. Instituto Geográfico
Nacional. Ministerio de Fomento.
Corchete V., Chourak M. and Khattach D. (2005). The high-resolution gravimetric geoid of
Iberia: IGG2005. Geophys. J. Int., 162, 676–684.
Corchete V. (2008). The high-resolution gravimetric geoid of North Iberia: NIBGEO. Terra
Nova, 20, 489–493.
Corchete V., Chourak M., Khattach D. and Benaim E. H. (2008). A new high-resolution
gravimetric geoid for South Spain and the Gibraltar Strait area: SOSGIS. Journal
of African Earth Sciences, 51, 145–150.
Forsberg R. and Kearsley A. H. W. (1989). Precise Gravimetric Geoid Computations over
Large Regions. Developments in Four-Dimensional Geodesy. Lecture Notes in Earth
Sciences Ser.: Vol. 29. Springer-Verlag.
Forsberg R. and Tscherning C. C. (1997). Topographic effects in gravity field modelling. In:
Geodetic boundary value problems in view of the one centimeter geoid, F. Sanso and R.
Rummel (Eds.), pp. 241-272, Lecture Notes in Earth Sciences, Springer.
Haagmans R., deMin E. and van Gelderen M. (1993). Fast evaluation of convolution
integrals on the sphere using 1D-FFT and a comparison with existing methods for
Stokes’ integral. Manuscripta Geodaetica, 18, 82-90.
Heiskanen W. A. y Moritz H. (1985). Geodesia Física. Instituto Geográfico Nacional e
Instituto de Astronomía y Geodesia.
Hofmann-Wellenhof B. and Lichtenegger H. (1994). Global Positioning System. Theory and
Practice. Springer-Verlag, Berlín.
Kaula W. M. (1987). The need for vertical control. Surveying and Mapping, 47, pp. 57-64.
Kearsley A. H. W. (1988). Tests on the recovery of the precise geoid height differences from
gravimetry. Journal of Geophysical Research, 93, pp. 6559-6570.
Kuroishi I. (1995). Precise gravimetric determination of geoid in the vicinity of Japan. Bull.
Geographical Surv. Inst., 41, 1-94.
Leick A. (1995). GPS satellite surveying. John Wiley & Sons.
Martín Asín F. (1983). Geodesia y Cartografía Matemática. Instituto Geográfico Nacional.
Schwarz K. P., Sideris M. G. and Forsberg R. (1987). Orthometric Heights without Levelling.
Journal of Surveying Engineering, 113, pp. 28-40.
42
Schwarz K. P., Sideris M. G. and Forsberg R. (1990). The use of FFT techniques in physical
geodesy. Geophys. J. Int., 100, pp. 485-514.
Sideris M. G. (1990). Rigorous gravimetric terrain modelling using Molodensky’s operator.
Manuscripta Geodaetica, 15, 97-106.
Smith D. A. and Roman D. R. (2001). GEOID99 and G99SSS: 1-arc-minute geoid models for
the United States. Journal of Geodesy, 75, 469-490.
Torge W., Basic T., Denker H., Doliff J., Wenzel H.-G. (1989). Long Range Geoid Control
through the European GPS Traverse. Deutsche Geodätische Kommission, Reihe B, nº
290.
Torge W. (1989). Gravimetry. Walter de Gruyter. Berlin-New York.
Torge W. (1991). Geodesy, 2nd Edition. Walter de Gruyter. Berlin-New York.
Vanicek P. and Krakiwsky E. (1986). Geodesy: The Concepts. Elsevier Science Publishers.
Yun H. S. (1994). Determination of a Gravimetric Geoid Solution for South Korea.
International Association of Geodesy Simposia. Symposia nº 113: Gravity and Geoid.
Springer-Verlag.
Zakatov P. S. (1981). Curso de Geodesia Superior. MIR.
43