Tema1 geodesia

Curso avanzado de posicionamiento por satélite Madrid, noviembre 2009
TEMA 1. Conceptos Fundamentales en Geodesia
1. Introducción.
En un curso en el que se desarrollan las principales características del GPS es
necesario revisar y recordar los conceptos fundamentales de la Geodesia, sobre todo lo
que se refiere a los Sistemas Geodésicos de Referencia y Sistemas de Coordenadas.
La geodesia es la ciencia que estudia la forma, dimensiones y campo gravitatorio de la
Tierra. En esencia, la Geodesia es la base de la Topografía y la Cartografía, de tal
forma que un estudio de estas dos siempre lleva implicado un conocimiento más o
menos riguroso de la Geodesia.
El campo tan amplio abarcado por la Geodesia hace que se pueda dividir en:
9 Geodesia Física: estudia el campo gravitatorio de la Tierra, partiendo de
mediciones del mismo y el modelado del mismo.
9 Geodesia Geométrica: estudia la figura de la Tierra en su aspecto geométrico,
dimensiones.
9 Astronomía Geodésica: estudia los métodos astronómicos que permiten
determinar las coordenadas geográficas sobre la superficie terrestre de una
serie de puntos fundamentales sobre los que se basarán las redes geodésicas.
9 Geodesia Espacial: utiliza las mediciones a cuerpos externos a la Tierra para
el posicionamiento y las mediciones geodésicas.
F. R. Helmert (1880) define la Geodesia como “la ciencia de la medida y
representación de la Tierra”. Esto implica en primer lugar, obtener un conocimiento de
la forma y dimensiones de la Tierra, pero por otra parte, también se podría decir que
es la ciencia que proporciona o determina coordenadas para cada punto de la
superficie terrestre.
Tema 1 - 1

2. La figura de la Tierra.
Clairaut, 1743, dice: “Théorie de la figure de la terre tireé des principles de
l´hydrostatique”. Desde esos días, la base del conocimiento de la figura de la Tierra no
ha variado, ya que sigue siendo un término moderno aplicado a la definición del
geoide, el cual es definido como la superficie (de nivel) equipotencial del campo
gravitatorio de la Tierra. Esta superficie, como media, coincide con el nivel medio del
mar (MSL) en un océano abierto sin perturbaciones o su extensión hipotética por
debajo de las masas continentales (Bomford, 1977; Clarke, 1880; Fischer, 1845;
Listing, 1873; Torge, 1991).
Un punto de la superficie terrestre está sometido casi exclusivamente a la atracción de
la Tierra y a la fuerza centrífuga derivada de su rotación, de tal manera que siendo V el
potencial gravitatorio y φ el potencial centrífugo, el potencial de la gravedad será:
W = V + φ
y por tanto será función de las coordenadas de cada punto. Los puntos para los cuales
W = cte definen una superficie equipotencial cuyo gradiente determina en cada punto
el vector de gravedad g, siendo la dirección la que define la vertical del lugar:
g = ∇W
La superficie equipotencial o de nivel materializado por los océanos cuando se
prescinde del efecto perturbador de las mareas (casi la superficie del nivel medio de
los mares) se denomina geoide y es, como todos sabemos, la superficie de referencia
para la altitud.
Veremos más adelante que la introducción del geoide como superficie de nivel cuyo
campo de gravedad verifica W = W0 = cte, tiene sentido físico, pero en Geodesia
Geométrica, y más en GPS, hemos de idealizar esta superficie en la práctica con
medidas, ya que su determinación física es uno de los problemas más intrincados de la
Geodesia.
Por todo ello, para la materialización de coordenadas planas, hemos de considerar la
Tierra como un elipsoide de revolución achatado por los polos que gira con movimiento
uniforme alrededor de uno de sus ejes principales de inercia (eje polar o menor).
Los sistemas tradicionales en los que se utiliza como Sistema de Referencia uno
bidimensional sobre un elipsoide de referencia pasan a ser en la actualidad, con las
Tema 1 - 2

nuevas técnicas de satélite a ser más complejos y sobre un sistema de coordenadas
cartesianas geocéntricas.
3. Sistemas elipsoidales de referencia.
Como ya se ha dicho, la superficie de la Tierra puede representarse con mucha
aproximación mediante un elipsoide de revolución, definiéndose este sistema con:
ƒ Superficie de referencia: dimensiones (semiejes a, b), excentricidad (e).
ƒ Ejes o líneas de referencia en la superficie.
ƒ Sentidos de medida.
A partir de las dimensiones de los semiejes del elipsoide a (semieje mayor) y b
(semieje menor) se definen los parámetros geométricos:
Achatamiento:
a
b
a
b
a
f −
=
−
= 1
Primera excentricidad:
a
b
a
e
2
2
−
= o como se suele expresar:
2
2
2
2
a
b
a
e
−
=
Segunda excentricidad:
b
b
a
e
2
2
'
−
=
Las líneas de referencia elegidas de la superficie elipsoidal son los meridianos y los
paralelos, con plano fundamental el del Ecuador (π) y como plano secundario, el plano
meridiano que contiene al meridiano origen (π').
Tema 1 - 3

Fig. 1. Coordenadas geodésicas en el elipsoide
Sobre esta superficie se definen las coordenadas geodésicas:
9 Latitud geográfica (ϕ): ángulo medido sobre el plano meridiano que contiene al
punto entre el plano ecuatorial y la normal al elipsoide en P.
9 Longitud geográfica (λ): ángulo medido sobre el plano ecuatorial entre el
meridiano origen y el plano meridiano que pasa por P.
La relación entre latitud geocéntrica (Φ), considerando el vector O (geocentro) al punto
P, y geodésica queda definida por:
φ
ϕ tan
tan
b
a
=
Si en lugar de considerar la normal al elipsoide en P consideramos la normal al geoide
(línea de la plomada), obtenemos las coordenadas astronómicas.
Las coordenadas astronómicas de un punto definen la dirección de la vertical de este
lugar con respecto al eje de rotación de la Tierra. Se obtienen por observación
astronómica y son independientes de cualquier hipótesis realizada sobre la forma de la
Tierra, ya que se refieren a un dato físico como es la vertical del lugar.
Sin embargo, las coordenadas geodésicas se obtienen mediante cálculos y
observaciones hechas sobre una superficie convencional (elipsoide de referencia),
referidos a la normal al elipsoide trazada por dicho lugar.
Una relación muy conocida que se aplicaba en las antiguas redes geodésicas para
efectuar la corrección geodésico-astronómico en los acimutes (puntos Laplace) era la
ecuación de Laplace:
Tema 1 - 4

( ) Φ
−
Λ
=
− sen
A
A g
a λ
Otra opción es conocer las componentes de la desviación de la vertical en el punto:
ϕ
ξ −
Φ
=
( ) ϕ
λ
η cos
−
Λ
=
3.1. Radios de curvatura del elipsoide de revolución.
El radio de curvatura de una curva plana y = f (x) se define como:
x
d
y
d
dx
dy
2
2
2
3
2
1
⎥
⎥
⎦
⎤
⎢
⎢
⎣
⎡
⎟
⎠
⎞
⎜
⎝
⎛
+
−
=
ρ
Operando sobre la ecuación del elipsoide y sustituyendo x e y por las coordenadas
cartesianas, obtenemos el radio de curvatura del meridano que pasa por el punto, ρ:
( )
( )2
3
2
2
2
sin
1
1
ϕ
ρ
e
e
a
−
−
=
El radio del paralelo p es:
ϕ
ϕ
2
2
sin
1
cos
e
a
p
−
=
Y la gran normal (N) es el radio de curvatura del primer vertical que pasa por el punto
(Tª Meusnier):
ϕ
2
2
sin
1 e
a
N
−
=
ya que N = p /cos ϕ.
Cualquier sección del elipsoide tendrá un
radio de curvatura comprendido entre N
(máximo) y ρ (mínimo).
Fig. 2. Radio principal de curvatura de la elipse meridiana
Tema 1 - 5

La curvatura de una sección normal arbitraria de acimut α se puede calcular según la
fórmula de Euler:
N
R
α
ρ
α
α
2
2
sin
cos
1
+
=
A partir de estas definiciones, las longitudes de arco de meridiano y paralelo se
definen:
dsm = ρ δϕ
dsp = N cos ϕ δλ
La longitud de arco de meridiano:
( )
∫ ∫
−
−
=
⋅
=
ϕ ϕ
ϕ
ϕ
δϕ
ρ
0 0 2
3
2
2
2
)
sin
1
(
1
e
d
e
a
sm
No existe una expresión cerrada para el cálculo de esta integral, utilizando diferentes
expresiones con desarrollos en serie para calcular la longitud de arco de meridiano.
Una de las más usadas puede ser (Snyder, 1979):
⎥
⎦
⎤
⎢
⎣
⎡
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
+
+
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
+
+
−
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
−
−
−
= ϕ
ϕ
ϕ 4
sin
1024
45
256
15
2
sin
1024
45
32
3
8
3
256
5
64
3
4
1
6
4
6
4
2
6
4
2
e
e
e
e
e
e
e
e
a
s
La longitud de arco de paralelo entre dos longitudes geográficas dadas λ1 y λ2:
∫ −
⋅
=
=
2
1
)
(
cos
cos 1
2
λ
λ
λ
λ
ϕ
ϕδλ N
N
sp
3.2. Reducción de medidas al elipsoide.
Al ser el elipsoide la superficie de referencia sobre la que se realizan los cálculos de
coordenadas geodésicas, es necesario reducir al mismo las mediciones realizadas sobre
el terreno.
En el caso de las medidas angulares, la reducción, aunque rigurosamente existe, es
prácticamente insignificante incluso para lados de primer orden (30 km).
La reducción clásica que se hace es en medidas de distancia:
1. Reducción al horizonte medio.
Tema 1 - 6

Suponiendo la distancia que une dos puntos de diferente altitud, la corrección a
realizar para obtener una distancia horizontal sobre el horizonte medio de ambos
puntos es:
3
4
2
8
2 D
h
D
h
c
Δ
−
Δ
−
=
siendo Δh el incremento de altitud entre ambos puntos y D, la distancia medida. Por
supuesto, la corrección es negativa (la distancia reducida es siempre inferior a la
medida):
D1 = D + c
siendo D1 la distancia reducida y D la distancia medida.
2. Reducción al nivel del mar.
El segundo paso sería la reducción de la distancia del horizonte medio al nivel del mar.
Estrictamente, tendría que hacerse al elipsoide, pero las altitudes estarán referidas
normalmente al nivel medio del mar y no al elipsoide (la diferencia es la ondulación del
geoide). Sin embargo, la ondulación del geoide en España con respecto al elipsoide de
Hayford es de unos 40 metros, lo cual no introduce variación sustancial en la
reducción. La distancia reducida al nivel del mar viene dada por:
m
1
2
h
R
R
D
D
+
⋅
=
siendo hm la altitud media de los dos extremos de la base y R puede ser el radio de la
esfera media, aunque estrictamente es el radio de la sección normal entre los dos
puntos, que viene dado por la fórmula de Euler:
z
cos
N
sen
ρ
ρ
N
R
N
α
sen
ρ
α
cos
R
1
2
2
α
2
2
α ⋅
+
⋅
⋅
=
⇒
+
=
z
3. Paso de la cuerda al arco.
La última corrección a hacer es el paso de la cuerda al arco, ya que la distancia para el
cálculo de coordenadas tiene que estar sobre el arco (del elipsoide). En este caso,
Tema 1 - 7

como se puede intuir, la corrección es positiva (la distancia sobre el arco es mayor que
la distancia sobre la cuerda):
2
3
2
2
3
24R
D
D
D +
=
En este caso, cabe hacer la misma consideración respecto a R.
Las correcciones sobre el horizonte medio y paso al nivel del mar se pueden hacer
también en un solo paso a partir de la fórmula:
2
/
1
2
1
2
2
2
1
1 ⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎣
⎡
⎟
⎠
⎞
⎜
⎝
⎛
+
⋅
⎟
⎠
⎞
⎜
⎝
⎛
+
Δ
−
=
R
h
R
h
h
D
D
3.3. Problemas directo e inverso sobre el elipsoide.
El cálculo de coordenadas a partir de una distancia y un acimut desde un punto a otro,
conocidas las coordenadas de uno de los puntos se denomina como problema directo
Tema 1 - 8

de la geodesia. El problema inverso sería la determinación del acimut y la distancia a
partir de las coordenadas de un par de puntos.
La resolución de estos problemas ha dado lugar a numerosos algoritmos propuestos
por diferentes autores (Levallois, Helmert, Coticchia, Rudoe, etc), cada uno con un
cierto grado de aproximación según los términos utilizados en cada fórmula.
Unas fórmulas bastante simples con buenos resultados hasta 30 km son las
aproximadas de Clarke:
X
X
N
p
ρ
⋅
⋅
⋅
⋅
=
2
cosA
A
sen
D2
q X
tan
tanA
p ϕ
⋅
⋅
=
q
−
⎟
⎠
⎞
⎜
⎝
⎛
−
⋅
=
Δ
Y
ρ
p
3
2
A
cos
D
ϕ
⎟
⎠
⎞
⎜
⎝
⎛
+
⋅
⎟
⎠
⎞
⎜
⎝
⎛
−
⋅
=
Δ
q
3
1
cos
N
p
3
1
A
sen
D
2
X ϕ
λ
donde x se refiere a
ρ
ϕ
ϕ
cosA
D
1
X
⋅
+
=
y se refiere a ( )
X
Y ϕ
ϕ
ϕ +
= 1
2
1
Para lados superiores a 30 km habría que utilizar otros desarrollos con más términos.
Por otra parte, el acimut recíproco del directo (la diferencia es la convergencia de
meridianos) se puede calcular a partir de la expresión:
( ) p
q
3
2
sen
λ
λ
180º
θ
θ 2
1
2
2
1
1
2 −
⎟
⎠
⎞
⎜
⎝
⎛
−
⋅
−
+
±
= ϕ
Como se ha dicho, estas fórmulas no dejan de ser aproximadas. Para un cálculo
riguroso a grandes distancias, es necesario aplicar más términos o resolver el
problema por otros métodos (punto auxiliar, método de Krasovsky, argumentos
medios con desarrollos en serie...) que en este resumen no tienen cabida.
Tema 1 - 9

El problema geodésico inverso (determinación de distancia y acimut a partir de las
coordenadas de dos puntos) es algo más complicado. Sin embargo, se pueden utilizar
también procedimientos aproximados más simples (cálculo de los arcos de meridiano y
paralelo o reducción a coordenadas UTM). Una aproximación bastante buena viene
dada por la denominada “fórmula topográfica” (errores inferiores a 0,3 m a distancias
hasta 40 km):
⎟
⎠
⎞
⎜
⎝
⎛
+
⋅
=
ΔA
2
1
A
cos
Δ
ρ
L m ϕ
o también:
⎟
⎠
⎞
⎜
⎝
⎛
+
⋅
⋅
=
ΔA
2
1
A
sen
cos
Δ
N
L m m
ϕ
λ
con ΔA = Δλ sen ϕm
3.4. Elipsoides de referencia.
La obtención de una superficie de referencia con un desarrollo matemático sencillo es
imprescindible para poder realizar la proyección de los puntos sobre la misma, de tal
forma que el geoide no puede ser esa superficie. El elipsoide de revolución que mejor
se adapte al geoide en la zona con un punto donde ambos coinciden o bien la normal a
ambos es la solución adoptada: elipsoide de referencia, constituyendo el concepto de
Datum geodésico el conjunto de ambos:
9 Una superficie de referencia (elipsoide de revolución).
9 Un punto fundamental en el que coinciden normales al geoide y al elipsoide.
En geodesia tendremos dos Datum: horizontal, visto hasta ahora, donde coincidirán las
coordenadas astronómicas y las geodésicas, y el vertical, que es la superficie de
referencia respecto a la que se definen las altitudes. En este caso, lo más normal es
que sea el geoide, con diferentes tipos de altitudes que veremos posteriormente.
A lo largo de la historia diversos elipsoides se han utilizado para definir el Sistema de
Referencia de cada país, de tal forma que se define aquel que mejor se ajuste al
geoide. Sucesivamente se han ido desarrollando esfuerzos encaminados a la adopción
de un sistema único, lo cual ha sido posible gracias al desarrollo de las técnicas
espaciales.
Tema 1 - 10

La elección del elipsoide se justifica en su sencillez matemática y el ajuste al geoide
con aproximación de primer orden, teniendo en cuenta que:
9 El centro gravitatorio terrestre debería coincidir con el centro del elipsoide.
9 El plano definido por el Ecuador terrestre debe coincidir con el del elipsoide.
9 La suma de los cuadrados de las ondulaciones del geoide debe ser mínima.
Bomford admite que aun adoptando valores arbitrarios para la latitud y longitud
geodésicas en el Datum, es necesario mantener el paralelismo entre el eje menor del
elipsoide y el de rotación de la Tierra.
Es evidente que la definición de numerosos elipsoides de referencia llevados a cabo en
los países dificulta notablemente la conexión de trabajos internacionales y cartografía,
de tal forma que desde principios del siglo XX se mantuvieron intenciones de
establecer una cartografía uniforme, siendo Hayford en 1910 el que, estableciendo un
elipsoide para la representación de EEUU, propuso en 1924 en la Asamblea
Internacional de Geodesia y Geofísica (Madrid) un Elipsoide Internacional de
Referencia, con a = 6378388 y α = 1/297.
La Asamblea General en Estocolmo adoptó para este elipsoide la fórmula de gravedad
internacional establecida por G. Cassinis correspondiente a la fórmula de gravedad
normal, esencial para la adopción de un elipsoide de nivel:
( ) 2
2
2
0 2
sin
0000059
.
0
sin
0052884
.
0
1
78049
.
9 −
⋅
−
+
= s
m
ϕ
ϕ
γ
Este elipsoide fue utilizado ampliamente por la mayoría de países, no siendo
perfeccionado hasta 1964, donde la Unión Astronómica Internacional, en Hamburgo
estableció unos nuevos valores de a = 6378160 y α = 1/298,25.
Posteriormente, en la Asamblea General de la U.G.G.I. en Lucerna (1967) el Sistema
de Referencia 1924 fue sustituido por el Sistema Geodésico de Referencia 1967 con las
constantes:
9 a = 6378160 m
9 GM = 398603 ⋅ 109
m3
s-2
9 J2 = 1082.7 ⋅ 10-6
9 w = 7.2921151467⋅ 10-5
rad s-1
En la orientación de este sistema se estipula:
Tema 1 - 11

a) El eje menor del elipsoide de referencia es paralelo a la dirección definida por
el origen internacional convencional (O.I.C.) para el movimiento del polo.
b) El meridiano de referencia es paralelo al meridiano cero adoptado por el BIH
para las longitudes (Greenwich).
Fig. 3. Principales elipsoides utilizados en cartografía
3.5. World Geodetic System 1984 (WGS84).
Desde 1987, el GPS utiliza el World Geodetic System WGS-84, que es un sistema de
referencia terrestre único para referenciar las posiciones y vectores. Se estableció este
sistema utilizando observaciones Doppler al sistema de satélites de navegación NNSS o
Transit, de tal forma que se adaptara lo mejor posible a toda la Tierra.
Se define del siguiente modo:
- Origen, centro de masas de la Tierra, incluyendo océanos y atmósfera.
Tema 1 - 12

- Eje Z paralelo a la dirección del polo CIO o polo medio definido por el BIH,
época 1984.0 con una precisión de 0,005”
- El eje X la intersección del meridiano origen, Greenwich, y el plano que pasa
por el origen y es perpendicular al eje Z, el meridiano de referencia coincide
con el meridiano cero del BIH en la época 1984.0 con una precisión de
0,005”. Realmente el meridiano origen se define como el IERS Reference
Meridian (IRM).
- El eje Y ortogonal a los anteriores, pasando por el origen.
- Terna rectangular dextrosum
Este sistema de referencia tiene asociado el elipsoide de referencia definido por los
siguientes parámetros:
Semieje mayor de la elipse a =6.378,137 km
Semieje menor de la elipse b = 6.356,752 km (parámetro derivado)
Factor de achatamiento f = 1/298,257223563
Velocidad angular de la tierra s
rad
E /
10
7292115 11
−
⋅
=
ω
Constante de gravitación
2
3
8
/
10
)
008
.
0
418
,
3986004
( s
m
GM ⋅
±
=
Fig. 4. Sistema de referencia WGS84
Tema 1 - 13

Es necesario hacer una distinción entre la definición de un sistema de coordenadas y la
realización práctica de un marco de referencia. Para conseguir una realización práctica
de un marco geodésico global de referencia se tienen que establecer una serie de
puntos con un conjunto de coordenadas. Un conjunto de puntos consistentes infieren:
• la localización de un origen
• la orientación del sistema de ejes cartesianos ortogonales
• una escala
En términos modernos, un conjunto de estaciones con coordenadas bien determinadas
constituyen o representan una realización de un Marco de Referencia Terrestre (TRF,
Terrestrial Reference Frame). El marco de referencia original de WGS84 fue establecido
en 1987 a través de un conjunto de estaciones NNSS (Navy Navigation Satellite
System) o TRANSIT (Doppler).
El principal objetivo a la hora de la realización del marco se centró en alinear lo más
preciso posible, el origen, la escala y la orientación del marco WGS84 con respecto al
Sistema Terrestre definido por el BIH (BTS, BIH Terrestrial System) en la época
1984.0.
Como se ha dicho, el marco inicial de WGS84 fue definido en enero de 1987 por la
Defense Mapping Agency (DMA) usando efemérides precisas de satélites TRANSIT con
estaciones de seguimiento Doppler. Esta realización estuvo vigente hasta 1994, ya que
varios estudios demostraron que existía un error sistemático entre las altitudes
elipsoidales derivadas de GPS y las de Doppler para un mismo punto (error en la
escala).
Se han realizado dos actualizaciones del WGS84, una de ellas denominada WGS84
(G730), en la época 1994, donde la G significa que se ha realizado por técnicas GPS y
730 es la semana GPS en la que se realizó la actualización, y otra denominada WGS84
(G873), en la época 1997.0.
Las nuevas exigencias de precisión hacen que este Sistema de Referencia, geocéntrico,
con una indeterminación en el geocentro de 1 m, con un eje Z variable, etc, hacen que
para la definición de Sistemas de Referencia más precisos sea necesario el uso de
Marcos de Referencia, los denominados ITRF, concepto que veremos posteriormente
con más detalle.
La razón proviene de cambios temporales de la corteza de la Tierra que deben ser
estimados o modelados. Los más importantes son el movimiento de las placas
Tema 1 - 14

tectónicas y los efectos de marea de la corteza de la Tierra. Como consecuencia de
estos efectos temporales, cualquier juego de coordenadas ha de referirse en el marco
ITRF vigente correspondiente a una determinada época.
Los parámetros geométricos y físicos que se derivan de los principales en WGS84 son:
Tabla 1. Parámetros de WGS84.
El elipsoide WGS84 definido como una superficie equipotencial tiene un potencial
gravitacional teórico (U). Este potencial puede ser calculado independientemente de la
distribución de masas en el interior, ya que es teórico, usando los parámetros y
constantes definidas para el elipsoide: semieje mayor (a), inversa del achatamiento
(1/f), velocidad angular de la Tierra (ω) y la constante gravitacional egocéntrica (GM):
Tema 1 - 15

- Gravedad normal en la superficie del elipsoide.
La gravedad normal teórica (γ) es el gradiente del potencial normal U y es dado en la
superficie del elipsoide por la fórmula cerrada de Somigliana:
ϕ
ϕ
γ
γ
2
2
2
sin
1
sin
1
e
k
e
−
+
=
donde 1
−
=
e
p
a
b
k
γ
γ
, siendo γe y γp la gravedad normal en el ecuador y en los polos
respectivamente.
Esta es la denominada Fórmula de la Gravedad del elipsoide WGS84, cuya superficie
equipotencial sirve como referencia horizontal y vertical, como figura geométrica de la
Tierra y como superficie de referencia para la gravedad normal de la Tierra.
- Gravedad normal en un punto por encima del elipsoide.
Se suele usar la forma de γ desarrollada en serie de Taylor con la derivada respecto a
la altura h:
2
2
2
2
1
h
h
h
h
h
∂
∂
+
∂
∂
+
=
γ
γ
γ
γ
Desarrollando la expresión se llega a:
⎥
⎦
⎤
⎢
⎣
⎡
+
−
+
+
−
= 2
2
2 3
)
sin
2
1
(
2
1 h
a
h
f
m
f
a
h ϕ
γ
γ
donde
GM
b
a
m
2
2
ω
= y γ es la gravedad normal en el elipsoide a latitud ϕ.
Esta fórmula sería válida siempre que h sea relativamente pequeño. Para altitudes muy
grandes o mayores exigencias de precisión se utilizan otros desarrollos.
- Transformaciones de Datum locales a WGS84.
Para aplicaciones poco precisas (cartografía de escalas pequeñas), la transformación
puede realizarse mediante las fórmulas de transformación estándar de Molodensky:
ϕWGS84 = ϕLOCAL + Δϕ
λWGS84 = λLOCAL + Δλ
Tema 1 - 16

hWGS84 = hLOCAL + Δh
con:
h
N
a
b
a
b
f
a
Ne
a
Z
Y
X
+
⎥
⎦
⎤
⎢
⎣
⎡
⎟
⎠
⎞
⎜
⎝
⎛
+
⎟
⎠
⎞
⎜
⎝
⎛
Δ
+
Δ
+
Δ
+
Δ
−
Δ
−
=
Δ
ϕ
ϕ
ρ
ρ
ϕ
ϕ
ϕ
λ
ϕ
λ
ϕ
ϕ
cos
sin
/
)
cos
sin
(
cos
sin
sin
cos
sin 2
ϕ
λ
λ
λ
cos
)
(
cos
sin
h
N
Y
X
+
Δ
+
Δ
−
=
Δ
ϕ
ϕ
λ
ϕ
λ
ϕ 2
sin
sin
sin
cos
cos
cos N
a
b
f
N
a
a
Z
Y
X
h ⎟
⎠
⎞
⎜
⎝
⎛
Δ
+
⎟
⎠
⎞
⎜
⎝
⎛
Δ
−
Δ
+
Δ
+
Δ
=
Δ
donde ΔX, ΔY, ΔΖ son las diferencias en coordenadas geocéntricas entre los centros del
datum local y el WGS84 y N y ρ los radios principales de curvatura.
4. Sistema de coordenadas cartesianas geocéntricas.
Para situar un punto en el espacio, se nos ocurre utilizar un sistema de coordenadas
basado en tres ejes perpendiculares entre
sí. El punto O, donde se cruzan los tres
ejes, lo colocaremos en el centro de la
Tierra. Al eje Z, lo colocamos de manera
que coincida con el eje de rotación del
planeta. Para terminar, haremos coincidir
el plano definido por OXZ con un
meridiano origen (Greenwich).
Fig. 5. Coordenadas cartesianas
Acabamos de definir un sistema geodésico de referencia. Con este sistema de
coordenadas, cualquier punto M puede ser situado con los valores de las componentes
X,Y,Z. A este tipo de coordenadas se les llama coordenadas tridimensionales
cartesianas.
Tema 1 - 17

Fig 6. Coordenadas cartesianas geocéntricas.
4.1. Paso coord. geográficas → Coord. Cartesianas
El problema planteado es pasar de coordenadas geodésica (ϕ, λ, h) al sistema
cartesiano elipsoidal, donde h es la altura del punto, pero medida sobre la normal al
elipsoide (proyección Helmert).
El radio vector OQ es:
⎟
⎟
⎟
⎠
⎞
⎜
⎜
⎜
⎝
⎛
−
=
⎟
⎟
⎟
⎠
⎞
⎜
⎜
⎜
⎝
⎛
=
ϕ
λ
ϕ
λ
ϕ
sin
)
1
(
sin
cos
cos
cos
2
e
N
N
N
ZQ
YQ
XQ
OQ
De la misma forma, h tiene la misma dirección que el vector OQ, por lo que:
⎟
⎟
⎟
⎠
⎞
⎜
⎜
⎜
⎝
⎛
=
=
ϕ
λ
ϕ
λ
ϕ
sin
sin
cos
cos
cos
h
h
h
QP
h
Tenemos que OP = OQ + QP, con lo cual tenemos finalmente:
⎟
⎟
⎟
⎠
⎞
⎜
⎜
⎜
⎝
⎛
+
−
+
+
=
⎟
⎟
⎟
⎠
⎞
⎜
⎜
⎜
⎝
⎛
=
ϕ
λ
ϕ
λ
ϕ
sin
)
)
1
(
(
sin
cos
)
(
cos
cos
)
(
2
h
e
N
h
N
h
N
Z
Y
X
OP
P
P
P
Tema 1 - 18

4.2. Paso Coordenadas cartesianas → Coord. Geográficas
En el problema inverso, a partir de coordenadas geocéntricas (X, Y, Z), obtener
coordenadas geográficas, hay que operar iterativamente, de tal forma que:
N
Y
X
h −
+
=
ϕ
cos
2
2
1
2
2
2
1
arctan
−
⎟
⎠
⎞
⎜
⎝
⎛
+
−
+
=
h
N
N
e
Y
X
Z
ϕ
X
Y
arctan
=
λ
El proceso iterativo empezaría entrando en la fórmula primera con h=0, deduciendo N.
Cálculo de N
Cálculo nuevo de ϕ
Cálculo de h
Cálculo de ϕ con h=0
La convergencia del sistema de ecuaciones es muy rápida, puesto que N >> h.
Existen también fórmulas aproximadas directas en las que no es necesaria la iteración:
θ
θ
ϕ 3
2
3
2
cos
sin
'
arctan
a
e
p
b
e
Z
−
+
=
X
Y
ar tan
=
λ
N
p
h −
=
ϕ
cos
donde θ es una cantidad auxiliar:
pb
Za
ar tan
=
θ
con
2
2
Y
X
p +
=
Tema 1 - 19

4.2. Transformación de Sistemas en coordenadas cartesianas.
Es evidente que este sistema de coordenadas cartesianas es más intuitivo y fácil de
manejar que las coordenadas geodésicas, no teniendo que tener en cuenta ninguna
superficie de referencia ni tampoco los parámetros de geometría del elipsoide.
Precisamente éste es el "sistema básico" sobre el que el GPS trabaja en coordenadas
(prescindiendo del sistema de coordenadas en el plano orbital).
También el paso de un Sistema Geodésico a otro es fácil de resolver haciendo el paso
intermedio de transformar las coordenadas geodésicas a coordenadas cartesianas en
cada sistema y hacer la transformación entre uno y otro sistema en cartesianas.
Para pasar de un
sistema geodésico a
otro, son necesarias
tres operaciones :
• una traslación
del centro de
la Tierra O
con un ΔX,
ΔY, ΔZ.
Fig. 7. Transformación entre sistemas cartesianos.
• una rotación R(si los tres ejes no son paralelos entre sí)
• un factor de escala, λ.
Dado dos sistemas geodésicos (1) y (2), y admitiendo la posibilidad que los dos
conjuntos de vectores de base no tengan la misma norma (problema de escala o de
definición de las unidades de distancia), la transformación más corriente será una
similitud (con 7 parámetros: 3 de translación, 3 de rotación y 1 de factor de escala).
Si (Xi Yi Zi) son las coordenadas de un punto P en el sistema "i" (con i=1 o 2), la
transformación será, matricialmente:
⎟
⎟
⎟
⎠
⎞
⎜
⎜
⎜
⎝
⎛
+
+
⎟
⎟
⎟
⎠
⎞
⎜
⎜
⎜
⎝
⎛
=
⎟
⎟
⎟
⎠
⎞
⎜
⎜
⎜
⎝
⎛
1
1
1
2
2
2
)
1
(
Z
Y
X
R
T
T
T
Z
Y
X
Z
Y
X
λ
Normalmente, los diferentes sistemas no son demasiado distantes, lo que
Tema 1 - 20

numéricamente, significa que:
¾ TX, TY, y TZ no exceden algunos centenares de metros.
¾ Que λ es de unas ppm (10-6
).
¾ Que R es una rotación infinitesimal de la forma:
⎟
⎟
⎟
⎠
⎞
⎜
⎜
⎜
⎝
⎛
−
−
−
=
1
1
1
X
Y
X
Z
Y
Z
R
R
R
R
R
R
R
donde RX, RY y RZ son rotaciones sobre X, Y, Z que no exceden algunos
segundos de arco.
La transformación puede ser remplazada por su modelo linealizado sin pérdida de
precisión:
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎠
⎞
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎝
⎛
⎟
⎟
⎟
⎠
⎞
⎜
⎜
⎜
⎝
⎛
−
−
−
=
⎟
⎟
⎟
⎠
⎞
⎜
⎜
⎜
⎝
⎛
Z
Y
X
Z
Y
X
Z
Y
X
R
R
R
T
T
T
X
Y
Z
X
Z
Y
Y
Z
X
D
D
D
λ
0
1
0
0
0
0
1
0
0
0
0
1
1
1
1
1
1
1
1
1
1
donde:
1
1
1
1
2
1
1
1
1
2
1
1
1
1
2
X
R
Y
R
Z
T
Z
Z
D
X
R
Z
R
Y
T
Y
Y
D
Y
R
Z
R
X
T
X
X
D
Y
X
Z
Z
Z
X
Y
Y
Z
Y
X
X
−
+
+
=
−
=
+
−
+
=
−
=
−
+
+
=
−
=
λ
λ
λ
Esta es la transformación básica general entre sistemas, si bien se estudiarán en
posteriores capítulos otras transformaciones más perfeccionadas y/o complejas de
forma más detallada (Badekas-Molodensky, DMA, NIMA, Bursa-Wolf, regresión
múltiple...) y también otras fórmulas para la transformación directamente en
coordenadas geodésicas. Se contemplarán asímismo los casos posibles que se pueden
presentar.
- Transformación de semejanza bidimensional (Helmert).
Tema 1 - 21

Se trata de una transformación bidimensional (2D) simple en el que se determinan dos
traslaciones, un giro y un factor de escala entre ambos sistemas.
i
Y
X
i
Y
X
T
T
Y
X
Y
X
⎥
⎦
⎤
⎢
⎣
⎡
⋅
⎥
⎦
⎤
⎢
⎣
⎡ −
⋅
+
⎥
⎦
⎤
⎢
⎣
⎡
+
⎥
⎦
⎤
⎢
⎣
⎡
=
⎥
⎦
⎤
⎢
⎣
⎡
1
1
1
1
2
2
cos
sin
sin
cos
α
α
α
α
λ
Si α es muy pequeño, entonces el sistema puede quedar como:
i
i
Y
X
i
Y
X
i
Y
X
Y
X
T
T
Y
X
T
T
Y
X
⎥
⎦
⎤
⎢
⎣
⎡
⋅
⎥
⎦
⎤
⎢
⎣
⎡ −
+
⎥
⎦
⎤
⎢
⎣
⎡
+
⎥
⎦
⎤
⎢
⎣
⎡
=
⎥
⎦
⎤
⎢
⎣
⎡
⋅
⎥
⎦
⎤
⎢
⎣
⎡ −
⋅
+
+
⎥
⎦
⎤
⎢
⎣
⎡
=
⎥
⎦
⎤
⎢
⎣
⎡
1
1
1
1
1
1
2
2
1
1
)
1
(
λ
α
α
λ
α
α
λ
Con dos puntos con coordenadas en ambos sistemas, se resuelve. Normalmente, se
suelen utilizar al menos cuatro puntos para resolver por mínimos cuadrados.
- Transformaciones de semejanza tridimensionales.
La transformación de semejanza tridimensional presenta problemas y condicionantes
ligeramente más complejos. La aplicación de las fórmulas en sistemas globales se hace
con coordenadas cartesianas geocéntricas, para lo cual es necesario aplicar la
expresión descrita anteriormente:
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎣
⎡
−
+
+
+
=
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎣
⎡
ϕ
ϕ
λ
ϕ
λ
ϕ
Nsen
e
h
N
sen
h
N
h
N
Z
Y
X
2
sin
)
(
cos
)
(
cos
cos
)
(
En donde h es realmente la altitud sobre el elipsoide y no sobre el geoide. El problema
se plantea aquí porque normalmente en el sistema local no se va a poder conocer la
altitud elipsoidal.
Los parámetros a determinar en este caso son 7: tres traslaciones, tres giros y un
factor de escala.
Por otro lado, al tratarse de sistemas globales, se supone que los ejes son "casi"
paralelos, en cuyo caso los cosenos del ángulo se igualan a 1 y los senos al ángulo.
El ámbito de aplicación de la transformación y su resolución da lugar a los dos modelos
más conocidos para estos propósitos.
- Transformación de Bursa-Wolf.
Este modelo es adecuado para transformación entre sistemas globales en un ámbito
geográfico amplio.
Tema 1 - 22

2
2
1
1
1
1
)
1
(
SR
x
y
x
z
y
z
SR
Z
Y
X
SR
Z
Y
X
T
T
T
Z
Y
X
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎣
⎡
⋅
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎣
⎡
Ω
−
Ω
Ω
Ω
−
Ω
−
Ω
+
+
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎣
⎡
=
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎣
⎡
λ
Preparando el sistema para su resolución por mínimos cuadrados (3 o más puntos
comunes en ambos sistemas de referencia) quedaría en la forma:
2
1
1
1
1
1
1
1
2
1
1
1
1
1
1
1
1
1
0
0
0
1
0
1
0
0
0
0
1
0
0
0
0
1
SR
i
SR
i
Z
Y
X
Z
Y
X
SR
i
i
i X
Z
Y
X
X
Z
Y
X
T
T
T
Y
Z
X
X
Y
Z
X
Z
Y
Y
Z
X
⎥
⎥
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎢
⎢
⎣
⎡
−
⎥
⎥
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎢
⎢
⎣
⎡
=
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎢
⎣
⎡
Ω
Ω
Ω
⎥
⎥
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎢
⎢
⎣
⎡
−
−
−
−
#
#
#
#
#
#
#
#
#
μ
Es decir, en la forma:
A ⋅ X = L
Y resolviendo:
X = (AT
⋅ A)-1
⋅ AT
⋅ L
- Transformación de Badekas-Molodensky.
Esta transformación es idónea para la aplicación a una zona limitada, como un país o
una región concreta. Para la resolución de este problema se divide en dos pasos, en los
que en uno se resuelven los giros y el factor de escala y en el siguiente, las tres
traslaciones. En ambos casos la solución se refiere al baricentro de los puntos, de tal
forma que se calculan en primer lugar las coordenadas de los baricentros de ambos
sistemas:
n
X
X
n
i
i
∑
=
= 1
n
Y
Y
n
i
i
∑
=
= 1
n
Z
Z
n
i
i
∑
=
= 1
Si se hace:
X
X
x −
= 1
1 Y
Y
y −
= 1
1 Z
Z
z −
= 1
1
Se plantean las ecuaciones para resolver los primeros cuatro parámetros:
Tema 1 - 23

⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎣
⎡
Ω
Ω
Ω
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎣
⎡
−
−
−
+
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎣
⎡
=
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎣
⎡
Z
Y
X
SR
SR
SR
x
y
z
x
z
y
y
z
x
z
y
x
z
y
x λ
2
1
1
1
1
1
1
1
1
1
2
1
1
1
1
1
1
1
0
0
0
#
#
#
#
#
#
Y finalmente se plantean las ecuaciones para resolver las tres traslaciones:
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎣
⎡
Ω
Ω
Ω
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎣
⎡
−
−
−
−
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎣
⎡
−
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎣
⎡
=
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎣
⎡
Z
Y
X
SR
SR
SR
Z
Y
X
X
Y
Z
X
Z
Y
Y
Z
X
Z
Y
X
Z
Y
X
T
T
T
μ
2
2
1
0
0
0
La aplicación práctica tanto en este caso como el anterior ha de hacerse con el
conocimiento de altitudes elipsoidales, lo cual en WGS84 es lo común, pero no lo es en
otros Datums locales, por lo que es necesario un proceso iterativo, en el cual en primer
lugar se aplica la transformación a los puntos sobre el geoide (en ambos sistemas). Si
se conoce la ondulación sobre WGS84, por ejemplo, se toma h=N para el conjunto de
puntos en este sistema y h aproximada para el Datum local también sobre el geoide,
resolviendo la transformación.
4.3. Sistema de coordenadas local y global.
Consideremos un sistema de coordenadas cartesianas locales con origen en el punto Pi
y las direcciones norte, este y arriba (ejes n, e, u) del plano tangente al punto.
Las coordenadas globales son iguales que las geocéntricas definidas anteriormente,
usando la notación vectorial X en lugar de las componentes X, Y, Z. De la misma
forma, los vectores i
X y j
X representan dos puntos de la superficie terrestre, Pi y Pj.
Definiendo el vector entre esos dos puntos en el sistema de coordenadas global como
i
j
ij X
X
X −
= , el vector en el sistema de coordenadas local referido al plano tangente
en P es definido como ij
x .
Tema 1 - 24

Fig. 8. Sistema de coordenadas local NEU.
En el sistema denominado NEU, los vectores n y e forman el plano tangente en P y la
tercera coordenada u es ortogonal, coincidiendo con la normal al elipsoide:
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎣
⎡
−
−
=
i
i
i
i
i
i
n
ϕ
λ
ϕ
λ
ϕ
cos
sin
sin
cos
sin
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎣
⎡−
=
0
cos
sin
i
i
i
e λ
λ
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎣
⎡
=
i
i
i
i
i
i
u
ϕ
λ
ϕ
λ
ϕ
sin
sin
cos
cos
cos
Realmente hay que tener en cuenta para ser más preciso, que el sistema local se
refiere a coordenadas astronómicas y no geodésicas y la normal al plano tangente en P
es normal al geoide, no al elipsoide, pero las diferencias para distancias pequeñas no
son apreciables.
Por consiguiente, el vector ij
x en el sistema NEU será:
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎣
⎡
⋅
⋅
⋅
=
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎣
⎡
=
ij
i
ij
i
ij
i
ij
ij
ij
ij
X
u
X
e
X
n
u
e
n
x
Uniendo los vectores i
i
i u
e
n ,
, del sistema local en la matriz i
D tenemos:
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎣
⎡
−
−
−
=
i
i
i
i
i
i
i
i
i
i
i
i
i
D
ϕ
ϕ
λ
ϕ
λ
λ
ϕ
λ
ϕ
λ
λ
ϕ
sin
0
cos
sin
cos
cos
sin
sin
cos
cos
sin
cos
sin
Tema 1 - 25

y por tanto la expresión anterior queda como:
ij
T
i
ij X
D
x ⋅
=
Se pueden también expresar los componentes de ij
x en función de la distancia espacial
sij, el azimut αij y el ángulo cenital zij mediante la relación:
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎣
⎡
=
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎣
⎡
=
ij
ij
ij
ij
ij
ij
ij
ij
ij
ij
ij
ij
z
s
z
s
z
s
u
e
n
x
cos
sin
sin
cos
sin
α
α
donde el azimut y el ángulo cenital se refieren a las medidas desde P.
Invirtiendo la expresión para obtener las medidas a partir de NEU:
2
2
2
ij
ij
ij
ij u
e
n
s +
+
=
ij
ij
ij
n
e
=
α
tan
2
2
2
cos
ij
ij
ij
ij
ij
u
e
n
u
z
+
+
=
5. Sistemas de altitudes.
5.1. El geoide como superficie de referencia.
La palabra geoide significa "forma de la Tierra" y fue introducida por Listing en 1873
como la superficie equipotencial que resulta de considerar una superficie equipotencial
de los océanos en reposo y prolongados por debajo de los continentes. Sería la
superficie de equilibrio de las masas oceánicas sometidas a la acción gravitatoria y a la
fuerza centrífuga ocasionada por la rotación de la Tierra, de manera que en cualquier
punto la dirección de la gravedad es perpendicular a esta superficie.
En este sentido, el geoide es una superficie definida por magnitudes físicas (el
potencial terrestre), teniendo en cuenta las anomalías gravimétricas (debidas a la
distribución de masas continentales y la diferente densidad de las cortezas) y el
achatamiento de los polos.
Tema 1 - 26

La diferencia en un punto determinado entre geoide y elipsoide se denomina
ondulación del geoide, concepto fundamental en la evaluación de altitudes
determinadas con GPS.
El ángulo existente en un punto determinado entre la normal al geoide y la normal al
elipsoide se denomina desviación de la vertical, siendo ésta nula en el Datum.
La desviación de la vertical se puede descomponer en dos componentes una en el
sentido de la latitud (ξ) y otra en el sentido de la longitud (η), siendo sus expresiones
en función de las coordenadas astronómicas y geodésicas:
ξ = ϕ − φ
η = (λ − Λ) cos φ
Fig. 9. Ondulación y desviación de la vertical
Las fórmulas de Vening-Meinesz son muy útiles, puesto que proporcionan las
componentes en función de la ondulación del geoide en el sentido de la longitud o
latitud:
δϕ
δ
ρ
ξ
N
⋅
−
=
1
δλ
δ
ϕ
η
N
N
⋅
−
=
cos
1
Tema 1 - 27

Y la desviación de la vertical total ε en una dirección de acimut α será:
α
η
α
ξ
ε sin
cos +
=
Estas expresiones también resultan muy útiles cuando se quiere mezclar observaciones
clásicas con GPS en una red para calcular, por ejemplo, altitudes elipsoidales, puesto
que habría que corregir los ángulos cenitales de desviación de la vertical (normal
geoide – normal elipsoide).
5.2. Determinación del geoide.
La determinación de la figura del geoide ha sido uno de los problemas fundamentales
de la Geodesia Física desde el siglo pasado, estribando sus principales dificultades en
la irregularidad del mismo sobre la superficie de la Tierra y la dificultad de obtener
valores de gravedad en la superficie terrestre con precisión. A todo esto ha de añadirse
la dificultad matemática y la complejidad del proceso.
Como se ha dicho anteriormente, el potencial total gravitatorio de la Tierra se define
como:
W = V + Φ
Siendo V el potencial gravitatorio y Φ el potencial debido a la rotación de la Tierra. Si ω
es la velocidad angular de rotación, entonces:
)
(
2
1 2
2
y
x +
= ω
φ
donde x e y son las coordenadas geocéntricas de un punto.
La función del potencial gravitacional se define como:
⎥
⎥
⎦
⎤
⎢
⎢
⎣
⎡
+
⎟
⎠
⎞
⎜
⎝
⎛
+
= ∑∑
= =
)
sin
cos
)(
'
(sin
1
max
2 0
λ
λ
φ m
S
m
C
P
r
a
r
GM
V nm
n
n
n
m
nm
nm
n
donde:
• r = distancia desde el centro de masas de la Tierra.
• a = semieje mayor elipsoide.
• n, m = grado y orden, respectivamente de los coeficientes gravitacionales.
• Φ’ = latitud geocéntrica.
Tema 1 - 28

• λ = longitud geocéntrica (o geodésica, da lo mismo).
• nm
nm S
C , = coeficientes gravitacionales normalizados = nm
C
k
n
m
n
m
n
⋅
⎥
⎦
⎤
⎢
⎣
⎡
+
−
+ 2
1
)
1
2
(
)!
(
)!
(
(lo mismo para Snm). Los coeficientes gravitacionales Cnm y Snm son para m=0
⇒ k=1 y m>1 ⇒ k=2.
• )
'
(sinφ
nm
P = Función de Legendre normalizada = )
'
(sin
)!
(
)
1
2
(
)!
( 2
1
φ
nm
P
m
n
k
n
m
n
⎥
⎦
⎤
⎢
⎣
⎡
+
+
−
• )
'
(sinφ
nm
P = Función de Legendre = [ ]
)
'
(sin
)
'
(sin
)
'
(cos φ
φ
φ n
m
m
m
P
d
d
• )
'
(sinφ
n = Polinomio de Legendre =
P n
n
n
n
d
d
n
)
1
'
(sin
)
'
(sin
!
2
1 2
−
φ
φ
La determinación de los coeficientes de los armónicos esféricos del potencial
gravitatorio se hace a partir del desarrollo de las anomalías de la gravedad Δg, las
cuales son determinadas a nivel global en toda la Tierra a partir de observaciones
gravimétricas o más concretamente, a partir de observaciones de satélites. Los
modelos geopotenciales de la Tierra nos dan estos coeficientes (normalizados) hasta
un determinado grado y orden (máximo, 360).
Es evidente que las series son válidas para r ≥ a, es decir, para la superficie de la
Tierra o para r mayor, no pueden usarse para masas en el interior de la Tierra.
Tema 1 - 29

Tabla 2. Coeficientes de EGM96.
Tema 1 - 30

Fig.10. GRACE gravity model.
Diferentes métodos han sido empleados para la determinación del geoide, siendo el
fundamental el método gravimétrico. Sin entrar en demasiados detalles sobre la
descripción del método, ya que sería muy arduo su estudio en este curso, digamos que
se basa en la evaluación en cada punto de gravedad observada de la evaluación de la
famosa integral de Stokes:
( ) ( )
∫∫ Ψ
⋅
Δ
=
σ
σ
λ
ϕ
πγ
d
S
g
R
Ngrav ,
4
donde R es el radio medio de la Tierra, γ es la gravedad normal del elipsoide, σ es la
esfera de integración, Δg son las anomalías de la gravedad reducidas al geoide y S(Ψ)
es la función de Stokes dada por:
( ) )
ln(
)
6
3
(
10
6
4
1 2
2
2
t
t
t
t
t
t
S +
−
−
+
−
−
=
Ψ con
2
sin
Ψ
=
t
siendo las anomalías de gravedad:
mod
g
g
g obs
red Δ
−
Δ
=
Δ
Δgobs son las anomalías gravimétricas observadas:
γ
δ
δ
δ −
+
−
+
=
Δ top
B
F
obs
obs g
g
Tema 1 - 31

donde gobs es la gravedad observada en el terreno, δF es la corrección aire libre o Faye,
δB es la corrección Bouguer simple, δtop es la corrección topográfica o Bouguer
compuesta (normalmente, efectuar con modelo digital del terreno) y γ es la gravedad
del elipsoide a esa latitud.
La anomalía del modelo Δgmod puede representarse en forma de armónicos esféricos
como el potencial perturbador T que resulta de restar el potencial gravitatorio terrestre
W del potencial gravitatorio del elipsoide U:
( ) Φ
+
⎪
⎭
⎪
⎬
⎫
⎪
⎩
⎪
⎨
⎧
+
⎟
⎠
⎞
⎜
⎝
⎛
+
= ∑ ∑
∞
= =
1 0
)
(cos
sin
cos
1
)
,
,
(
n
n
m
nm
nm
nm
n
P
m
S
m
C
r
a
r
GM
r
W θ
λ
λ
λ
θ
donde r, θ, λ son las coordenadas polares (radio geocéntrico, colatitud y latitud), GM es
la constante gravitacional geocéntrica, a es el semieje mayor del elipsoide, Pnm(cosθ)
es el polinomio de Legendre normalizado de primera clase, Cnm son los coeficientes
armónicos esféricos normalizados del campo gravitacional terrestre y Φ es el potencial
centrífugo de la Tierra, Φ = (1/2)ω2
r2
sin ϕ.
( ) Φ
+
⎪
⎭
⎪
⎬
⎫
⎪
⎩
⎪
⎨
⎧
⎟
⎠
⎞
⎜
⎝
⎛
+
= ∑
∞
=1
2
2
cos
1
'
)
,
(
n
n
n
P
r
a
r
GM
r
U θ
θ
siendo M' en este caso la masa del elipsoide de referencia. El potencial perturbador
viene dado por:
( )
⎪
⎭
⎪
⎬
⎫
⎪
⎩
⎪
⎨
⎧
Δ
+
Δ
⎟
⎠
⎞
⎜
⎝
⎛
=
−
= ∑ ∑
∞
= =
2 0
)
(cos
sin
cos
)
,
,
(
n
n
m
nm
nm
nm
n
P
m
S
m
C
r
a
r
GM
U
W
r
T θ
λ
λ
λ
θ
donde ΔCnm y ΔSnm son las diferencias entre los coeficientes del geopotencial y el
potencial del elipsoide. La diferencia M-M' es muy pequeña y M' puede ser reemplazada
por M.
Con esto llegamos a la expresión de la anomalía gravimétrica del modelo geopotencial:
( )
⎪
⎭
⎪
⎬
⎫
⎪
⎩
⎪
⎨
⎧
Δ
+
Δ
−
⎟
⎠
⎞
⎜
⎝
⎛
=
Δ ∑ ∑
∞
= =
2 0
2
)
(cos
sin
cos
)
1
(
)
,
,
(
n
n
m
nm
nm
nm
n
P
m
S
m
C
n
r
a
r
GM
r
g θ
λ
λ
λ
θ
Se suele utilizar la ecuación de Bruns siguiente para el cálculo de ondulaciones del
geoide:
Tema 1 - 32

U
W
T
N
δ
−
=
y por tanto:
( )
⎪
⎭
⎪
⎬
⎫
⎪
⎩
⎪
⎨
⎧
+
⎟
⎠
⎞
⎜
⎝
⎛
= ∑ ∑
∞
= =
2 0
)
(cos
sin
cos
)
,
,
(
n
n
m
nm
nm
nm
n
P
m
S
m
C
r
a
r
GM
r
N θ
λ
λ
γ
λ
θ
Con n=m=360 se llegan a obtener unos 10000 coeficientes. El más usado de todos
ellos puede ser el OSU91A o el EGM96, este último más moderno.
En cualquier caso, se trata de modelos útiles para obtener anomalías aire libre y
calcular anomalías modelo, pero para la obtención del modelo de geoide gravimétrico,
en la práctica se divide la ondulación N en tres componentes:
N = Nmod + Nind + Ngra
donde Nmod es la contribución del modelo de armónicos esféricos (EGM96, OSU91A,
etc), Nind es la contribución de la reducción al terreno (efecto indirecto o cogeoide) y
Ngra es la contribución de las observaciones del campo gravitatorio terrestre: anomalías
aire libre después de quitar el efecto del modelo geopotencial global y la topografía
(corrección topográfica con modelo digital del terreno).
Como se ve, el tema es complejo y se han de tener en cuenta más variables no
mencionadas aquí, la resolución de la integral de Stokes necesita técnicas avanzadas
(transformada de Fourier ...). Veremos posteriormente cómo el geoide resultante ha
de ser escalado, ya que se hacen una serie de suposiciones respecto a la masa del
elipsoide, densidad media de la corteza, etc que originan unos valores que han de ser
ajustados a partir de observaciones reales sobre el terreno. El GPS juega aquí un papel
fundamental en este sentido, de tal forma que por un lado la definición de un geoide
preciso necesita del GPS y por otro lado, el inverso, también se cumple: el GPS
necesita de un modelo de geoide.
Tema 1 - 33

Fig. 11. Modelo EGM96
Fig. 12. Modelo EGG97 (Europa).
Tema 1 - 34

Actualmente, los satélites de observación terrestre han facilitado la labor en cuanto a
la definición del geoide a través de sus observaciones (ver ejemplo de
TOPEX/POSEIDON).
En el caso del modelo global EGM96, la ondulación viene dada por:
)
,
(
)
,
(
)
,
,
(
)
,
( 0 λ
ϕ
γ
λ
ϕ
λ
ϕ
ξ
λ
ϕ H
g
r
N
N B
Δ
+
+
=
con:
• N0 = -0.53 m (esta cantidad es la diferencia entre el semieje mayor de WGS84
y el elipsoide que mejor se ajusta a la Tierra).
• ΔgBA (φ, λ) = anomalía Bouguer = Anomalía aire libre – 0.1119 H (ϕ, λ).
• Anomalía aire libre = g + 0.3086 H (ϕ, λ) - γ.
• γ = valor medio de gravedad normal.
• H (ϕ, λ) = definido a partir del MDT JGP95E (Joint Gravity Project 95).
EGM96 está públicamente disponible, con valores de ondulación del geoide con un
rango de error de entre 0,5 y 1,0 metros (1 sigma) a nivel mundial. Los productos de
EGM96 como modelo geoidal para WGS84 incluyen:
• Un fichero de ondulaciones del geoide sobre WGS84 de 15’ x 15’.
• Los valores de los coeficientes armónicos esféricos completos hasta grado y
orden 360.
Tema 1 - 35

Fig. 13. Software de NIMA para modelo EGM96.
5.3. Sistemas de altitudes.
Intuitivamente parece fácil calcular la altitud de un punto cualquiera de la superficie
terrestre. Bastaría partir del mareógrafo más cercano e ir determinando por cualquier
método las diferencias de altitud entre puntos. Sin embargo, un análisis más profundo
del problema nos hace llegar a la conclusión que el problema no es tan sencillo.
En principio, el nivel medio de los mares queda más o menos bien determinado por la
existencia en zonas concretas y adecuadas de mareógrafos, que registran a lo largo de
un período histórico este nivel, y a los que volveremos con más detalle más adelante.
Sin embargo, esta superficie, el geoide, presenta problemas al tratar de extenderla
matemáticamente por debajo de los continentes, pues el potencial exterior de las
masas verifica la ecuación de Laplace:
ΔU
U
x
U
y
U
z
= + + =
∂
∂
∂
∂
∂
∂
2
2
2
2
2
2
0
mientras que en el interior de las masas se verifica la ecuación de Poisson:
Tema 1 - 36

Tema 1 - 37
ΔU
U
x
U
y
U
z
G
= + + = −
∂
∂
∂
∂
∂
∂
π ρ
2
2
2
2
2
2
4
donde para calcular la superficie U = cte., es necesario integrar 4 π G ρ, siendo ρ la
densidad, con el inconveniente de que no conocemos su distribución en todos los
puntos.
donde para calcular la superficie U = cte., es necesario integrar 4 π G ρ, siendo ρ la
densidad, con el inconveniente de que no conocemos su distribución en todos los
puntos.
En la nivelación geométrica, como el nivel se equilibra con su eje óptico tangente a la
superficie equipotencial que pasa por el centro del nivel, si hiciéramos una nivelación
entre dos puntos situados en una misma superficie equipotencial, obtendríamos una
diferencia de altitud nula, y esto es evidentemente falso desde el punto de vista
geométrico, debido a que las líneas de fuerza del campo tienen una cierta curvatura.
Como las superficies equipotenciales son perpendiculares a ellas, estas no son
paralelas, y por tanto, sus puntos no equidistan de la superficie original, también
equipotencial.
En la nivelación geométrica, como el nivel se equilibra con su eje óptico tangente a la
superficie equipotencial que pasa por el centro del nivel, si hiciéramos una nivelación
entre dos puntos situados en una misma superficie equipotencial, obtendríamos una
diferencia de altitud nula, y esto es evidentemente falso desde el punto de vista
geométrico, debido a que las líneas de fuerza del campo tienen una cierta curvatura.
Como las superficies equipotenciales son perpendiculares a ellas, estas no son
paralelas, y por tanto, sus puntos no equidistan de la superficie original, también
equipotencial.
Por ello podemos tener dos conceptos para medir la altitud:
Por ello podemos tener dos conceptos para medir la altitud:
- Distancia a una superficie de referencia (altitud geométrica).
- Distancia a una superficie de referencia (altitud geométrica).
- Dos puntos están a la misma altura si se encuentran sobre la misma superficie
equipotencial (altitud dinámica).
- Dos puntos están a la misma altura si se encuentran sobre la misma superficie
equipotencial (altitud dinámica).
Cuando estacionamos un nivel y hacemos lecturas a dos miras dispuestas
verticalmente, estamos midiendo diferencias de cotas geométricas.
Cuando estacionamos un nivel y hacemos lecturas a dos miras dispuestas
verticalmente, estamos midiendo diferencias de cotas geométricas.
Suponemos dos puntos alejados A y
B de la superficie terrestre. Haciendo
la nivelación cerrada veremos que la
suma de todas las diferencias de
altitudes medidas no será 0, aún
salvando los errores inherentes a
toda medida. Esto es debido a que
matemáticamente la integral de
diferencia de altitudes a lo largo de
un circuito cerrado no es cero: Fig. 14. Superficies equipotenciales
Suponemos dos puntos alejados A y
B de la superficie terrestre. Haciendo
la nivelación cerrada veremos que la
suma de todas las diferencias de
altitudes medidas no será 0, aún
salvando los errores inherentes a
toda medida. Esto es debido a que
matemáticamente la integral de
diferencia de altitudes a lo largo de
un circuito cerrado no es cero: Fig. 14. Superficies equipotenciales
4
δn ≠
∫ 0

Tema 1 - 38
Esto es debido a que el incremento de nivelación δn es distinto al incremento δn’,
debido a la falta de paralelismo de las superficies de nivel.
Puesto que el campo gravitatorio terrestre es conservativo, la diferencia de potencial
no depende del camino. La diferencia de potencial entre los puntos A y B es:
- δW = g δn = g’ δn’
siendo g la gravedad en la estación y g’ la gravedad sobre la línea de la plomada en B.
δ δ δ
n
g
g
n n
'
'
'
= ⇒ ≠δn
Esto indica que las alturas medidas por un camino y otro son distintas, la diferencia de
altitudes medida depende del camino recorrido. Por todo ello, en trabajos geodésicos,
donde las distancias son largas y las líneas equipotenciales convergen (no son
paralelas), no tiene sentido hablar de altitudes geométricas si no añadimos el
componente físico de la gravedad.
5.3.1. Cota Geopotencial.
Considerando un punto O sobre el geoide (a nivel del mar) y un punto A, conectado a
O mediante una línea de nivelación, se define la cota geopotencial como la diferencia
de potencial entre los dos puntos:
C W W g n
A
O
A
= − = ⋅
∫
0 δ
La cota geopotencial definida pues de esta forma, es independiente del camino
recorrido. Es la misma para todos los puntos de una superficie de nivel, puede
considerarse como una medida natural de altitud, aunque no tenga dimensiones de
longitud. Se mide en unidades geopotenciales (u.g.p.): 1 u.g.p.= 1 Kgal m = 1000 g m
= 1 gal Km.
A partir de la definición del número geopotencial se pueden definir algunos tipos de
alturas, que algunos autores llaman alturas geopotenciales. La definición de estas
alturas será la relación o cociente entre el incremento del potencial y un determinado
valor de la gravedad.
H
C
g
i
i
i
=

Según los potenciales elegidos y el valor de g las alturas más usuales son: alturas
aproximadas, alturas ortométricas y alturas dinámicas.
5.3.2. El potencial normal.
Para poder determinar algunas alturas se ha de definir previamente lo que es el
potencial normal. El potencial gravitatorio terrestre es:
W G
dv
r
x y
V
= + +
∫
ρ ω2
2 2
2
( )
El conocimiento riguroso del potencial implica conocer en cada punto de la tierra el
valor de su densidad. Ante esta dificultad interesa a veces descomponer el valor del
potencial en dos términos: uno de fácil cálculo (U) que se acerque lo más posible al
valor del potencial real, y el otro la diferencia (T). Así, W = U + T, siendo U el llamado
potencial normal y T el potencial perturbador o potencial anómalo.
Para la determinación de U se eligen diversos modelos de Tierra, de forma que su
potencial normal U en cada punto sea muy similar al verdadero W. Se adopta por ello
como modelo de Tierra un elipsoide de revolución con una masa igual a la masa real
de la Tierra y con una velocidad de rotación ω igual a la terrestre.
El teorema de Stokes hace posible el cálculo del potencial normal, sin ninguna
ambigüedad, al conocerse la superficie de nivel exterior y la masa.
Como consecuencia de la definición del potencial normal surge la gravedad normal
como γ = grad U. Para tener una idea del orden de las aproximaciones, en general la
relación T/W es del orden de 10-5
.
5.3.3. Altitudes dinámicas.
La altitud dinámica de un punto introduce un factor de escala en el concepto de cota
geopotencial y hace que el resultado venga expresado en metros. La altitud dinámica
de un punto es:
H
C
din
=
γ 0
Tema 1 - 39

Tema 1 - 40
Sería el resultado de dividir C por una constante, gravedad normal a una latitud
estándar (45º), convirtiendo así la cota geopotencial en una longitud. Se usa para la
determinación de γo la fórmula de 1967:
γ1967=978.0318(1+0.0053024sen2
φ − 0.0000059 sen2
2φ )
lo que para una φ = 45º da un valor de 0,98061898.
5.3.4. Altitudes ortométricas.
Si designamos por P0 la intersección del geoide con la línea de la plomada que pasa por
el punto P y H la altitud ortométrica, esto es, la longitud del segmento de línea a lo
largo de la línea de la plomada entre P0 y P, podemos efectuar la siguiente integración
a lo largo de la línea de la plomada, puesto que el resultado es independiente de la
trayectoria:
C g d
H
= ⋅
∫
0
H
multiplicando y dividiendo por H obtenemos:
C H
H
g dH
H
= ⋅ ⋅ ⋅
∫
1
0
de modo que C g H
= ⋅ , donde:
g
H
g dH
H
= ⋅ ⋅
∫
1
0
es el valor medio de la gravedad sobre la línea de la plomada entre el geoide, (punto
P0), y el terreno, (punto P), que también puede escribirse como:
g
H
g z dz
H
= ⋅ ∫
1
0
( )
Para convertir por tanto los resultados de nivelación en altitudes ortométricas,
necesitamos conocer la gravedad dentro de la Tierra. Esta debe calcularse con los
únicos valores que realmente conocemos, que son los valores de la gravedad en la
superficie, y se hace reduciendo los valores medidos de la gravedad por el método de
Poincaré-Prey.

Tema 1 - 41
Observando la figura, Q es el punto
en que se ha de calcular g’, tal que
g’=gQ. Sea P el correspondiente
punto de la superficie de modo que
P y Q están situados sobre la
misma línea de la plomada. La
gravedad en P, se supone medida.
Para calcular gQ se usaría la
fórmula:
g g
g
h
dH
Q p
P
Q
= − ∫
∂
∂
Fig. 15. Altitudes ortométricas
suponiendo que el gradiente real de la gravedad δg/δh dentro de la Tierra fuese
conocido. Este puede obtenerse por la fórmula de Bruns:
∂
∂
π ρ ϖ
g
h
gJ G
= − + −
2 4 2 2
si se conoce la curvatura media J de las superficies geopotenciales y la densidad entre
P y Q.
La reducción de Prey supone que las curvaturas medias de las superficies
geopotenciales J coincide con la curvatura media de las superficies equipotenciales del
campo de gravedad normal J0. Por ello:
∂γ
∂
γ ϖ
h
J
= − −
2 2
0
2
Entonces:
∂
∂
∂γ
∂
π ρ
g
h h
G
= + 4
Numéricamente, despreciando la variación de ∂γ ∂
/ h con la latitud, y tomando ρ =
2,67 gr/cm3
para la densidad y G=66,7x 10-9
unidades c.g.s, resulta:
∂
∂
g
h
gal km
= − + = −
0 3086 0 2238 0 0848
, , , / ,
llegando a la expresión final de la reducción de Poincaré-Prey como:

gQ = gP + 0.0848 ( Hp - HQ )
expresando g en gales y H en kilómetros.
Esta fórmula no tiene rigor absoluto, en tanto en cuanto estamos igualando gJ = γ J0,
sin embargo se usa en la práctica casi siempre, ya que el conocimiento de la curvatura
media real J de las superficies geopotenciales es bastante imperfecto.
Con estos condicionantes, ya se puede sustituir en la fórmula de la altitud ortométrica
el término de la gravedad media a lo largo de la línea de la plomada:
( )
z
H
g
z
g −
⋅
+
= 0848
.
0
)
(
donde g es la gravedad medida en el punto P del terreno, de tal forma que ahora
puede integrarse directamente, resultando:
[ ]
H
0
2
z
H
0
2
z
H
0848
.
0
H
1
g
dz
)
z
H
(
0848
.
0
g
H
1
g ⎥
⎦
⎤
⎢
⎣
⎡
−
⋅
+
=
−
+
= ∫
g g H
= + ⋅
0 0424
. (g en gals, H en Km)
Con esta aproximación obtenemos las llamadas altitudes ortométricas Helmert:
H
g
g H
H gH C
=
+
⇒ + − =
0 0424
0 0424 0
2
.
.
Esta aproximación sustituye el terreno por una lámina Bouguer infinita de densidad
constante y altura H, siendo suficiente la precisión para redes altimétricas de pequeña
magnitud y zonas no demasiado montañosas.
5.3.5. Altitudes normales.
Este tipo de altitudes fueron introducidas por Molodensky, de tal forma que la
superficie del geoide se sustituye por otra superficie (quasi-geoide). Supone que el
campo de gravedad de la Tierra es normal, es decir, W=U, g=γ, T=0. Bajo esta
hipótesis se sustituye en la fórmula de altitud ortométrica la gravedad por la gravedad
normal del elipsoide γ y las altitudes pasan a ser normales H*
:
∫
∫ γ
=
⇒
γ
=
⇒
⋅
γ
=
−
=
C
0
*
*
H
0
*
0
C
H
dC
H
dH
W
W
C
Tema 1 - 42

Para calcular γ existen fórmulas desarrolladas, llegando a la expresión de H*
en
función de C y los parámetros del elipsoide y la latitud:
⎥
⎥
⎦
⎤
⎢
⎢
⎣
⎡
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
γ
+
γ
φ
−
+
+
+
γ
=
2
2
*
a
C
a
C
)
sin
f
2
m
f
1
(
1
C
H
donde m=ω2
a/γecuador (Teorema de Clairaut) y γ es la gravedad normal del elipsoide
para la latitud φ.
5.3.6. Altitudes elipsoidales.
En el caso del GPS, las altitudes se refieren siempre al elipsoide de referencia, GRS80
(WGS84). Este tipo de altitudes se denominan elipsoidales, ya que están medidas
desde el punto de la superficie terrestre al elipsoide, medidas sobre la normal al
elipsoide que pasa por el punto.
Es evidente que no tienen un significado geométrico ni mucho menos físico, tal como
los sistemas anteriores. Es necesario un
perfecto conocimiento del geoide para
relacionar ambos tipos de altitudes, ya que:
h = H - N
conocida relación en GPS, siendo h la altitud
elipsoidal, H la ortométrica y N la ondulación
del geoide en el punto. Por ello se hace
imprescindible el conocimiento de algún
modelo de geoide de suficiente precisión en la
nivelación con GPS.
Fig. 16. Altitud elipsoidal.
Por aspectos de desconocimiento de valores absolutos del geoide, lo normal es utilizar
la relación anterior en la práctica en modo diferencial Δh = ΔH + ΔN.
Hemos definido el geoide como la superficie equipotencial que más se aproxima
al nivel medio del mar. Esta definición tan sólo tiene un significado físico (nunca
geométrico en el sentido estricto) y por tanto el geoide no es en realidad el nivel medio
del mar. Sin embargo, en el sentido práctico geodésico, hemos de admitirlo como tal
haciendo válidas las definiciones clásicas.
Tema 1 - 43

Tema 1 - 44
Al hacer una transformación del sistema GPS a un sistema geodésico local o nacional,
hemos de tener en cuenta este aspecto, ya que una transformación global clásica
tridimensional puede introducir unos errores considerables cuando trabajamos a
distancias grandes. La superficie del geoide queda enmascarada y desconocida muchas
veces en una transformación tridimensional, de tal forma que si se conoce el geoide
con suficiente precisión, es recomendable establecer una transformación 2D para
planimetría y una 1D introduciendo la ondulación del geoide en altimetría. Estos
aspectos serán detallados más adelante.
Aunque vamos adentrándonos poco a poco en los aspectos que nos interesan del GPS
en esta introducción, vamos a ver los métodos que tendríamos de dar altitudes con
GPS. Atendiendo exclusivamente a la posibilidad o no de tener un modelo de geoide y
sin tener en cuenta nº de receptores, extensión y geometría de la zona, etc:
-1. Se dispone de carta del geoide local con la precisión necesaria.
Nivelación GPS:
H1-H2=h2-h1-N2+N1 ----> ΔH = Δh + ΔN
• Precisión requerida de al menos de 1 ppm en N, no fácilmente alcanzable en muchas
áreas.
• Permite eliminar errores en h (obs. simultánea, troposfera....)
- 2. Se dispone de puntos con H conocida en el entorno de trabajo.
Utilizándolos como puntos de control, dotándolos de altitud elipsoidal con GPS y
realizando interpolación, se obtendrían perfiles del geoide en la zona.
Se obtienen excelentes resultados en trabajos de ámbito local y si el número de puntos
de control es reducido, en zonas donde sepamos que la topografía del geoide no es
muy movida (se puede comprobar y detectar anomalías con modelo geopotencial
previo).
3. Se dispone de un modelo de geopotencial.
• Normalmente suelen ser globales y no dan la precisión necesaria, aunque sí un
conocimiento aproximado de N en términos absolutos.
)
,
(
)
,
(
)
,
,
(
)
,
( 0 λ
ϕ
γ
λ
ϕ
λ
ϕ
ξ
λ
ϕ H
g
r
N
N B
Δ
+
+
=

• Generalmente suelen ser de orden y grado 360. Para mayor precisión se necesitaría
mayor grado, orden 10000.
Estos modelos, aunque nominalmente suelen estar en precisiones de ± 1 m, es
evidente que por su carácter global (mallas 15’x15’) no pueden recoger anomalías
locales y suelen dar perfiles planos del geoide.
Sin embargo, pueden ser útiles para un posicionamiento global absoluto del trabajo.
- 4. Se dispone de un levantamiento gravimétrico.
Se está trabajando actualmente mucho en esta técnica para obtener deformaciones
locales de un geoide ya conocido con menor precisión. Para ello se combinan
observaciones GPS con las anomalías gravimétricas para obtener deformaciones
locales más detalladas.
5. No se dispone de nada de lo anterior: Transformación tridimensional (7 parámetros,
Borsa-Wolf, Badekas-Molodenski...).
Esta clasificación siempre va a depender de la amplitud del área en que estemos
trabajando y de la densidad de puntos de control que tengamos. En un área pequeña
(topografía), normalmente se hace esto último (transformación 3D) con buenos
resultados, siempre que tengamos un número suficiente de puntos con doble juego de
coordenadas.
Conclusión: la transformación WGS84-Sistema Local (en distancias > (depende)), debe
ser una Transformación 2D + 1D, y no 3-D. El conocimiento del geoide “vuelve” a ser
cuestión fundamental de la Geodesia.
Ejemplo: Dotación de altitud ortométrica a los aeropuertos de España.
En 1999 AENA-Eurocontrol solicitó al IGN la determinación de altitudes ortométricas a
puntos situados en cabecera y pie de pista de 40 aeropuertos de toda España.
Mediante la utilización de perfiles del geoide IBERGEO95 desde dos clavos NAP (red de
nivelación de alta precisión) situados a distancias de hasta 25 Km, a sendos puntos del
aeropuerto, se tenía una doble comprobación de la altitud ortométrica de cada punto.
Tema 1 - 45

Fig. 17. Esquema de trabajo.
Esquema de trabajo:
6 receptores, 2 vértices REGENTE (red geodésica GPS), 2 NAP, 2 Señales Aeropuerto.
- Transferencia de h desde vértices REGENTE a NAP.
- Transferencia de h desde NAP a Señales Aeropuerto.
Fig. 15. Esquema de trabajo
- Obtención de N obs en NAP.
- Cálculo de DN entre NAP y Señales Aeropuerto.
- Cálculo de H Señales Aeropuerto.
HAE1 = HNAP1 - ( hNAP1 - hAE1 ) + ( NNAP1 - NAE1 )
HAE1 = HNAP2 - ( hNAP2 - hAE1 ) + ( NNAP2 - NAE1 )
-Se obtuvieron así con GPS altitudes ortométricas con precisión centímetrica, precisión
que solo se hubiera podido alcanzar con técnicas clásicas de nivelación.
Tema 1 - 46

6. Sistemas geodésicos modernos: Marcos y Sistemas de
Referencia.
6.1. Introducción.
Muchos han sido los cambios conceptuales y de metodología en todas las Ciencias en
los últimos años debido al gran desarrollo tecnológico surgido. En lo que respecta a las
Ciencias de la Tierra, desde que en 1957 comenzara la carrera espacial con el
lanzamiento del Sputnik-1, los científicos dedicados tanto a la Geodesia como a la
Geofísica empezaron a intuir las posibles aplicaciones y aprovechamientos de tales
artilugios. Tanto en el estudio de la forma y dimensiones de la Tierra (Geodesia), como
en el estudio de los fenómenos físicos que afectan y condicionan dicha forma y
dimensiones (Geofísica), los cambios, sobre todo en los últimos quince años, han sido
espectaculares debido a las aplicaciones surgidas del aprovechamiento de satélites
artificiales.
Dentro de la Geodesia, los métodos de posicionamiento por satélite han cambiado por
completo la metodología de esta ciencia. El Sistema de Posicionamiento Global (GPS)
es hoy en día herramienta fundamental para todo tipo de trabajos geodésicos.
La investigación y el tratamiento de estos sistemas de posicionamiento por satélite, ha
llevado en la actualidad a que sean utilizados para fines científicos, destacando el
estudio de la atmósfera terrestre, de sus capas, fenómenos, y muy especialmente para
el estudio de la ionosfera, desconocida en muchos aspectos y con una gran influencia
sobre los distintos fenómenos que ocurren en nuestro planeta.
Pero quizá, las aplicaciones en las cuales estos sistemas han calado más hondo son la
Geodesia y la Topografía, a partir del descubrimiento de que dichos sistemas de
posicionamiento podían aportar las precisiones requeridas para el desarrollo de estas
ciencias y su aplicación en el desarrollo de infraestructuras, cartografía, sistemas de
información geográfica, estudios de movimientos y deformaciones, etc.
Actualmente se acuña una nueva definición de Geodesia: "Ciencia que determina la
forma y dimensiones de la Tierra y sus componentes dinámicas".
A finales de los años 80, la demanda de información geográfica georreferenciada en un
sistema geodésico de referencia uniforme creció de forma masiva debido a la
disponibilidad y versatilidad del GPS en aplicaciones de topografía, navegación,
transporte y logística, entre otras. En esta época, la red europea más uniforme hecha
Tema 1 - 47

hasta el momento era el European Datum ED87, calculado por la Comisión RETrig de la
IAG (International Association of Geodesy), como resultado de la combinación y el
reajuste de las redes nacionales de triangulación. Ni este ni el World Geodetic System
1984 (WGS84), realizado a través de observaciones GPS y con una precisión de unos
metros, cumplían las expectativas.
Con los nuevos conceptos de Geodesia, se establece una nueva jerarquía de Redes
Geodésicas (Work Group VIII, de CERCO), que según el tipo de observaciones de que
se disponga se puede dividir en redes geodésicas de (IAG, 1991):
- Clase A: Conjunto de puntos integrados en un determinado Marco de
Referencia con campos de velocidades, con σ < 1 cm, independiente de la
época de observación, garantizada por las estaciones permanentes GPS.
- Clase B: Redes Continentales Fundamentales sin campos de velocidades, σ <
1 cm, pero garantizada solo en una época específica, en sitios donde los
campos de velocidades no estén definidos, establecidas generalmente por
campañas GPS puntuales desde 1993.
- Clase C: Redes apoyadas en Clase B o densificación de las mismas, σ = 5 cm.
De acuerdo con estas nuevas definiciones en las que la componente dinámica debida a
los Sistemas de Referencia (que posteriormente veremos) es tenida en cuenta, las
redes de clase A quedan definidas únicamente por Estaciones Permanentes GPS con un
tiempo de observación y proceso de datos relativamente considerable (2 años).
La definición y densificación de los Marcos y de Referencia actuales en los que se
apoyan los Sistemas de Posicionamiento Espacial sólo son posibles a través de las
Estaciones Permanentes de técnicas espaciales: GPS, VLBI, DORIS y SLR/LLR, como
posteriormente se describirá.
Como ya se ha dicho anteriormente, uno de los aspectos fundamentales de la Geodesia
moderna es la definición y el mantenimiento del Marco de Referencia Terrestre. La
precisión y la realización del mismo es de fundamental importancia para estudiar las
propiedades de la Tierra local y globalmente, incluidos el rebote post-glacial, cambios
del nivel medio del mar, placas tectónicas, subsidencias regionales, deformaciones en
límites de placas, etc.
Los problemas que se abordan en cuanto a las deformaciones terrestres (Geodinámica)
tienen que ligarse a un sistema global y de ahí la importancia de un marco de
referencia perfectamente determinado.
Tema 1 - 48

Un sistema de referencia queda definido por su origen, tres direcciones fijas en el
espacio y el sistema de medida. Para relacionar los sistemas de referencia terrestres o
no inerciales con los espaciales o inerciales, ambos deben estar perfectamente
relacionados geométricamente y en el tiempo, lo que constituye una de las razones
fundamentales para explicar la necesidad, en primer lugar, de mantener una base
precisa de tiempos. Ello es debido a que en el Universo no hay un marco de referencia
que permanezca fijo, en el transcurso del tiempo, por lo que las direcciones que se
eligen como sistema de referencia, se refieren a un instante determinado en el tiempo.
La consecuencia importante de todo ello es que WGS84 no sirve para trabajos
geodésicos de precisión: indefinición centro masas de la Tierra y atmósfera,
movimiento del polo, movimiento de placas tectónicas... y es necesario definir
los Sistemas y Marcos de Referencia.
Los efectos temporales más importantes a tener en cuenta en la definición de un
Sistema Convencional de Referencia Terrestre serían cuatro:
• movimiento del Polo.
• movimiento de placas tectónicas.
• mareas de la Tierra solida.
• carga oceánica.
Estos cuatro efectos son casi perfectamente modelados en la definición de los
parámetros que intervienen en un CTRS (y por supuesto, en las observaciones GPS
derivadas), pero en observaciones GPS de gran precisión con líneas base muy grandes,
el software científico de cálculo introduce las oportunas correcciones.
6.1.1. Movimiento del Polo.
La intersección del eje instantáneo de rotación de la Tierra se mueve con el tiempo en
lo que se denomina movimiento del Polo (polar motion).
Tema 1 - 49

Figura 18. Movimiento del polo 1996 - 2000.5.
El movimiento es prácticamente periódico, con un movimiento principal de unos 434
días (periodo de Chandler) y cuya amplitud no excede de unos 10 metros. A partir de
un modelo geofísico se puede prever perfectamente el movimiento, aunque,
lógicamente, no de una manera completa.
Para evitar variaciones en longitud y latitud provenientes del movimiento del Polo, se
define el Conventional Terrestrial Pole (CTP), fijado a la corteza de la Tierra y
originalmente definido como el Polo medio entre los años 1900 y 1906. El eje de
rotación instantáneo se refiere al CTP mediante las coordenadas del Polo (xp, yp),
concepto que veremos posteriormente con más detalle.
6.1.2. Movimiento de placas tectónicas.
Para mantener una estabilidad a nivel centimétrico en un CTRS (Sistema Convencional
de Referencia Terrestre), las coordenadas de un conjunto de estaciones deben
establecerse para una determinada época y actualizarlas de movimiento de placas
tectónicas.
El movimiento observado ha sido como máximo de 7 cm/año (DoD) en las estaciones
GPS. Una forma de controlar los movimientos horizontales es estimar parámetros de
velocidad en las estaciones.
Tema 1 - 50

El modelo de movimiento de placas recomendado es conocido como NNR-NUVEL1A, en
el que se establecen 16 placas tectónicas principales, como se muestra en la figura
(IERS Conventions, 1996).
Fig. 19. Principales placas tectónicas del modelo NNR-NUVEL1A.
El modelo es una mejora de la versión original NUVEL-1 y consiste en una lista de
componentes de velocidad angular cartesiana para cada una de las 13 placas
principales, de tal forma que la matriz de rotación R es:
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎣
⎡
Ω
Ω
−
Ω
−
Ω
Ω
Ω
−
=
Ω
0
0
0
)
(
x
y
x
z
y
z
R
y la transformación entre dos épocas se define mediante (McCarthy, 1996):
[ ] )
(
)
(
)
(
10
84813681
.
4
)
( 0
0
9
t
x
t
t
R
I
t
x ⋅
−
⋅
Ω
⋅
⋅
+
= −
La definición de los ITRF se hace entonces con campos de velocidades en las
estaciones de seguimiento y para el resto, puede aplicarse el modelo.
ITRFyy
t
ITRFyy
t
ITRFyy
t v
t
t
x
x ,
0
,
, 0
0
)
( −
+
=
6.1.3. Efectos de marea.
Los fenómenos de marea son otra fuente de desplazamiento permanente de las
estaciones y son modelados según se indican en las convenciones IERS. Es evidente
Tema 1 - 51

que la mayoría de programas comerciales no tienen en cuenta estos efectos, que son
anulados en líneas base relativamente cortas en GPS diferencial, cuyos efectos, en
valor absoluto, pueden llegar a ser de hasta 42 cm. En aplicaciones rigurosas y en
cálculo de redes mundiales, debe ser introducido un modelo de corrección.
Las mareas son causadas por la variación temporal de la atracción gravitacional del Sol
y la Luna debido al movimiento orbital. Mientras que las mareas oceánicas son mucho
más evidentes en las zonas costeras y sus inmediaciones, las mareas de la Tierra
sólida (terrestres) son relativamente fáciles de calcular con modelos simples. Sus
periodicidades pueden ser derivadas directamente del movimiento de cuerpos celestes,
de forma similar a la nutación.
Las mareas terrestres generan desplazamientos periódicos de las estaciones que
dependen de la latitud. La variación puede ser como mucho de unos 30 cm en la
componente vertical y 5 cm en la horizontal.
La siguiente expresión proporciona los desplazamientos de mareas terrestres
(McCarthy, 1996):
( ) ( ) ( )
∑
=
⎥
⎦
⎤
⎢
⎣
⎡
⋅
−
⋅
+
⎟
⎠
⎞
⎜
⎝
⎛
−
=
Δ
3
2
2
2
2
4
)
(
3
2
1
2
3
j
j
j
j
j
j
E
E
j
e
e
r
r
e
r
l
e
r
e
h
r
r
GM
GM
x
donde:
• GME es la constante gravitacional geocéntrica de la Tierra.
• GMj es la constante gravitacional geocéntrica de la Luna (subíndice j=2) y el sol
(subíndice j=3).
• e es el vector unitario de la estación en el sistema de coordenadas geocéntricas
(x).
• r es el vector unitario del cuerpo celeste.
• h2 y l2 son los números de grado 2 de Love y Shida que describen las
propiedades elásticas de la Tierra.
Esta ecuación proporciona los valores de marea terrestre con una precisión de 5 mm.
6.1.4. Carga oceánica.
Se denomina carga oceánica a la deformación del suelo del mar y de la tierra costera
como consecuencia de la redistribución del agua del mar que tiene lugar durante la
marea oceánica. La corteza terrestre se deforma ante el peso del agua de marea.
Tema 1 - 52

El desplazamiento en un lugar en componente radial (Δc) en un tiempo t viene dado
por la expresión (McCarthy, 1996):
∑ Φ
−
+
+
=
Δ
j
cj
j
j
j
cj
j u
t
A
f
c )
cos( χ
ω
El sumatorio sobre j representa once ondas de marea tradicionalmente designadas
como semidiurnas M2, S2, N2, K2, componentes diurnas K1, O1, P1 y de largo periodo Mf,
Mm, Ssa.
ωj y χj las velocidades angulares y los argumentos astronómicos en el tiempo t = 0h
. χj
representa la posición del Sol y la Luna.
fj y uj dependen de la longitud del nodo de la Luna.
Acj y Φj son las amplitudes y las fases respectivamente de la estación, las cuales
pueden ser calculando usando modelos de marea oceánica.
El IERS proporciona estos valores para las estaciones ITRF y también el Astronomical
Institute de la Universidad de Berna tiene un servicio automático que proporciona los
coeficientes mediante la descarga en su servidor del fichero de coordenadas de las
estaciones que queramos.
También en http://www.oso.chalmers.se/∼loading (Observatorio de Onsala) se pueden
obtener los coeficientes a partir de diferentes modelos de carga (el más utilizado es el
GOT00.2).
Fig. 20. Tabla de coeficientes para carga oceánica en Onsala.
Tema 1 - 53

Para más información sobre este y los anteriores efectos de una forma más
pormenorizada, consultar “IERS Conventions 2003”, International Earth Rotation
Service (http://www.iers.org).
6.2. Sistema Internacional de Referencia Celeste (ICRS).
En principio, los sistemas inerciales, es decir, los que no están sometidos a
aceleraciones, son los más apropiados para definir la posición de objetos externos a la
Tierra, como las estrellas, cuásares o satélites artificiales. Sin embargo, no es posible
considerar un sistema totalmente inercial, ya que no se tendrían en cuenta los efectos
relativistas, que dependen del campo de la gravedad de la Tierra, estando sometido a
aceleraciones por el movimiento de traslación de la Tierra alrededor del Sol. Este
sistema se denomina “cuasi-inercial”.
Otro inconveniente es el movimiento del Polo, es decir la variación que se detecta
tanto en módulo como en dirección del vector velocidad angular. Este movimiento es la
suma de dos componentes:
9 la nutación astronómica con respecto al sistema de coordenadas inercial
celeste.
9 movimiento del Polo con respecto al sistema de referencia terrestre.
A su vez, el movimiento del Polo se divide en otros dos:
9 Rotación libre o período de Chandler, con una duración de 430 días y un radio
respecto a la posición media del Polo de 6 metros.
9 Movimiento diario de unos 50 cm debido al movimiento forzado por la acción
gravitatoria del Sol y la Luna.
Fig. 21. Movimiento del Polo y periodo de Chandler.
Tema 1 - 54

El Sistema de Referencia Celeste Internacional (ICRS), tiene como origen el baricentro
del sistema solar y es el conjunto de algoritmos y constantes con procedimientos y
parámetros aprobados internacionalmente: fórmulas de precesión y nutación,
coeficientes de precesión, amplitudes, velocidad de la luz, etc. La dirección del polo es
definida por la International Astronomical Union IAU, según la teoría de Precesión
(Lieske,1977) y la teoría de nutación (Seidelmann, 1982), para la época J2000.
Asociado a este Sistema existe un Marco Internacional de Referencia Celeste (ICRF),
definido por coordenadas ecuatoriales de precisión en el sistema J2000.0, de fuentes
extragalácticas determinadas con VLBI. Origen en el baricentro del sistema solar.
6.3. Sistema Convencional de Referencia Terrestre (CTRS).
Los sistemas de referencia terrestres o fijos a la Tierra se utilizan para determinar
coordenadas de puntos sobre la superficie terrestre o en sus proximidades, por lo que
están afectados por el movimiento de rotación de la tierra, apareciendo la aceleración
centrífuga y la de Coriolis, por lo cual se puede decir que son sistemas no inerciales.
Es evidente que el principal problema de estos sistemas es conseguir una dirección del
eje Z que sea fija, invariable con el tiempo, y para ello se adoptó por convenio un polo
medio.
El polo medio es el definido como el eje medio de rotación correspondiente al centro
del movimiento libre del polo entre los años 1900.0 y 1906.0. Este polo se denomina
Conventional International Origin CIO o Conventional Terrestrial Pole, CTP.
El sistema de referencia Terrestre convencional CTRS, se define del siguiente modo:
• Origen: Centro de masas terrestre o geocentro.
• Eje Z: coincide con el Polo terrestre Convencional CTP.
• Eje X: Perpendicular al anterior, que pasa por el origen y por el
meridiano de Greeenwich en la época 1903.0, adoptado como meridiano
de referencia.
• Eje Y: Perpendicular a los dos ejes anteriores y contenido en el ecuador
medio, formando un sistema de coordenadas trirrectángulo.
La diferencia entre el CTP y el CEP o polo instantáneo se obtiene de la observación
directa y el procesamiento de estaciones permanentes multitécnica y es uno de los
Tema 1 - 55

productos que proporciona el Servicio Internacional de Rotación de la Tierra
(International Earth Rotation Service, IERS), a partir del seguimiento de la variación
de coordenadas a lo largo del tiempo. Para el cálculo de esta variación de coordenadas
se utilizan diferentes métodos: Very Long Baseline Interferometry (VLBI), Satellite
Laser Ranging (SLR), Lunar Laser Ranging (LLR), GPS y DORIS/PRARE.
En definitiva, el ITRS es geocéntrico, en el centro de masas de la Tierra sólida
incluyendo oceános y atmósfera. La unidad de longitud es el metro y la precisión de la
dirección definida por el IRM (IERS Reference Meridian) es de ± 0.0005". La evolución
en el tiempo de esta orientación es calculada usando la condición de "no rotación" en
los movimientos tectónicos de toda la tierra.
El CTRS es la conexión práctica entre el ICRS y un sistema terrestre, un sistema de
coordenadas ligado firmemente al globo terráqueo. El Sistema Convencional de
Referencia Terrestre adoptado para el análisis de conjuntos de datos individuales
obtenidos mediante técnicas de observación (VLBI, SLR, LLR, GPS, DORIS, PRARE),
como para la combinación de soluciones individuales en un conjunto unificado de datos
(coordenadas de estaciones, parámetros de orientación de la Tierra, etc) sigue, de
acuerdo con la IAG, los siguientes criterios:
9 CTRS es un sistema local de coordenadas dentro del contexto de la teoría de la
relatividad, válido para la Tierra con una precisión relativa de 10-10
.
9 CTRS es un sistema cuasi-cartesiano en rotación, definido a partir de un
sistema geocéntrico sin rotación.
9 El tiempo coordinado es el Tiempo Geocéntrico Coordinado (TGC).
9 CTRS es geocéntrico, con centro de masas de la Tierra incluyendo la atmósfera.
9 Orientación por el eje Z del "BIH1984.0" dirigido hacia el polo IERS (igual que el
CIO) y los ejes X y Z forman el plano del meridiano cero.
9 Su evolución temporal en orientación no generará una rotación global residual
con respecto a la corteza terrestre (esto es la condición de NNR, no-net-
rotation).
6.4. Parámetros de Orientación de la Tierra.
Por todo lo dicho anteriormente es fundamental incluir en el Sistema de Referencia
Terrestre, en el cual hacemos todas las medidas, el concepto de los Parámetros de
Orientación de la Tierra (EOP).
Tema 1 - 56

La orientación de la Tierra es definida a partir de un sistema coordenado ligado al
geocentro de la Tierra Gxyz (el Sistema Terrestre materializado por las coordenadas
observadas de las estaciones) y un sistema coordenado ligado al espacio inercial GXYZ
(el Sistema Celeste, materizalizado por las coordenadas de estrellas, quásares y
objetos del Sistema Solar). La forma de describir la rotación de la Tierra es especificar
la matriz de rotación entre esos dos sistemas.
Si la tierra rotase uniformemente alrededor de un eje de rotación fijo (con respecto a
la corteza de la Tierra y el Sistema Celeste), los cambios de la rotación podrían ser
descritos únicamente con un parámetro: la variación del ángulo de rotación
linealmente con el tiempo o bien la escala de tiempo que se deriva de este ángulo de
rotación (tiempo universal). En realidad, el eje de rotación no se fija con respecto a la
corteza ni al Sistema Celeste. Los cambios en el vector de rotación de la Tierra son
causados por efectos gravitacionales causados por la Luna, el Sol y los planetas,
desplazamientos de materia en diferentes partes del planeta y otros mecanismos. Las
oscilaciones observadas pueden ser interpretadas en términos de elasticidad del
manto, irregularidades en la Tierra, propiedades de la frontera núcleo-manto, reología
del núcleo, aguas subterráneas, variabilidad oceánica y atmosférica en términos de
clima.
En principio, la orientación terrestre podría ser descrita a través de tres ángulos
independientes (por ejemplo los ángulos de Euler). Sin embargo, la monitorización
continua de la rotación de la Tierra considera por separado el movimiento del eje de
rotación en la Tierra y en el espacio. La determinación actual de esto se hace a través
de los cinco Parámetros de Orientación de la Tierra (Earth Orientation Parameters,
EOP). Los EOP del IERS describen las irregularidades de la rotación de la Tierra.
Técnicamente, son los parámetros que dan la rotación del ITRS con respecto al ICRS
como una función del tiempo.
Estos parámetros están referidos al Polo Celeste de Efemérides (Celestial Ephemeris
Pole, CEP), que está cercano al polo de rotación. Su movimiento en el espacio es la
precesión y la nutación. Prácticamente su movimiento es equivalente al del polo
geográfico (polo ITRF) para periodos inferiores a dos días. La posición del CEP es
perfectamente modelada con una precisión de 0.001".
¾ dPsi, dEps-. El movimiento del polo celeste, descrito por los modelos de
precesión y nutación. Las diferencias observadas reflejan la diferencia entre el
movimiento del polo celeste real con respecto al predicho por los modelos de
Tema 1 - 57

precesión y nutación de la IAU (Unión Astronómica Internacional). La mayoría
de estas variaciones son explicadas por errores en unos pocos términos:
términos seculares, términos periódicos (principalmente 18.6, 1, 0.5 años y 14
días). Estos términos están incluídos en el modelo de las Convenciones IERS
(1996), las cuales están de acuerdo con los modelos en magnitudes por debajo
de 0.001". Vienen indicados como (dPsi, dEps) y son observados por estaciones
VLBI. En definitiva, Dpsi y Depsilon dan los desplazamientos (offsets) en
longitud y oblicuidad del polo celeste con respecto a las posiciones definidas por
los modelos convencionales de la IAU.
¾ UT1-. El CEP permanece muy cerca del eje de rotación instantáneo (desviación
por debajo de 0.02") y por tanto es necesario un recálculo del ángulo de
rotación de la Tierra en el espacio. El IERS nos da el ángulo de rotación de la
Tíerra, pero en forma de escala de tiempo UT1, que es el tiempo requerido para
una rotación completa si la Tierra rotase uniformemente (360º en 86164.09891
s). Las tablas nos dan a los usuarios el desplazamiento con respecto a una
escala de tiempo uniforme TAI y UTC: UT1-TAI y UT1-UTC. El exceso del
periodo de rotación con respecto al periodo medio se denomina exceso de
longitud del día (LOD).
¾ xp , yp-. Después de aplicar la matriz de correcciones de precesión - nutación
por los desplazamientos del polo celeste y el ángulo de rotación alrededor del
CEP, queda una rotación residual con respecto al International Reference
Frame. Está asociado con dos pequeños ángulos que son las coordenadas del
polo CEP relativo al eje Gz del Marco Terrestre Internacional de Referencia
(International Terrestrial Reference Frame, ITRF), que es el eje de rotación o
polo geográfico. Esto se debe mayormente a la redistribución de masas
oceánicas y atmosféricas relacionadas con cambios estacionales. Se reflejan en
las Coordenadas del polo, x e y del CEP respecto al Polo de Referencia del IERS
(IERS Reference Pole, IRP). El CEP difiere del eje de rotación instantáneo en
unos términos quasi-diurnos con amplitudes por debajo de 0.01" (Seidelmann,
1982). El eje x está en la dirección del Meridiano de Referencia del IERS (IERS
Reference Meridian, IRM); el eje Y está perpendicularmente a este hacia el
oeste.
El IERS determina rotaciones diferenciales que son xp, yp , y que relacionan las
coordenadas del siguiente modo:
Tema 1 - 58

CTP
p
p
p
p
CEP
Z
Y
X
y
x
y
x
Z
Y
X
⎟
⎟
⎟
⎠
⎞
⎜
⎜
⎜
⎝
⎛
⋅
⎟
⎟
⎟
⎠
⎞
⎜
⎜
⎜
⎝
⎛ −
=
⎟
⎟
⎟
⎠
⎞
⎜
⎜
⎜
⎝
⎛
1
1
0
0
1
Tabla 3. Formato de los ficheros ERP del IERS.
IERS C04 POLE 13-JUN-02 21:01
---------------------------------------------------------------------------------------------------------------------------------------------------------------
CELESTIAL POLE OFFSET: NO SUBDAILY POLE MODEL: RAY
DATE TIME X-POLE Y-POLE UT1-UTC GPS-UTC RMS XP RMS YP RMS DT DE-CPO DP-CPO RMS EP RMS PS
YYYY MM DD HH MM (") (") (S) (S) REM (") (") (S) (") (") (") (")
2002 1 1 0 0 -0.17699 0.29400 -0.115826 13. C04 0.00000 0.00000 0.000000 -0.00067 -0.05316 0.00000 0.00000
2002 1 2 0 0 -0.17747 0.29761 -0.116703 13. C04 0.00000 0.00000 0.000000 -0.00074 -0.05295 0.00000 0.00000
2002 1 3 0 0 -0.17832 0.30085 -0.117711 13. C04 0.00000 0.00000 0.000000 -0.00084 -0.05287 0.00000 0.00000
2002 1 4 0 0 -0.17909 0.30370 -0.118895 13. C04 0.00000 0.00000 0.000000 -0.00095 -0.05294 0.00000 0.00000
2002 1 5 0 0 -0.17956 0.30649 -0.120257 13. C04 0.00000 0.00000 0.000000 -0.00107 -0.05315 0.00000 0.00000
2002 1 6 0 0 -0.17991 0.30956 -0.121650 13. C04 0.00000 0.00000 0.000000 -0.00113 -0.05343 0.00000 0.00000
2002 1 7 0 0 -0.18063 0.31273 -0.122914 13. C04 0.00000 0.00000 0.000000 -0.00113 -0.05372 0.00000 0.00000
2002 1 8 0 0 -0.18172 0.31556 -0.123968 13. C04 0.00000 0.00000 0.000000 -0.00108 -0.05393 0.00000 0.00000
2002 1 9 0 0 -0.18271 0.31828 -0.124808 13. C04 0.00000 0.00000 0.000000 -0.00100 -0.05398 0.00000 0.00000
2002 1 10 0 0 -0.18303 0.32118 -0.125421 13. C04 0.00000 0.00000 0.000000 -0.00096 -0.05388 0.00000 0.00000
2002 1 11 0 0 -0.18304 0.32423 -0.125791 13. C04 0.00000 0.00000 0.000000 -0.00100 -0.05368 0.00000 0.00000
2002 1 12 0 0 -0.18276 0.32740 -0.126048 13. C04 0.00000 0.00000 0.000000 -0.00116 -0.05347 0.00000 0.00000
2002 1 13 0 0 -0.18224 0.33044 -0.126131 13. C04 0.00000 0.00000 0.000000 -0.00136 -0.05328 0.00000 0.00000
... ... ... ... ... ... ... ... ... ...
6.5. International Terrestrial Reference Frame (ITRF).
Las realizaciones del Sistema Terrestre Internacional de Referencia (International
Terrestrial Reference System, ITRS) son producidas por la Sección de Marcos
Terrestres del IERS bajo el nombre de Marco Terrestre Internacional de Referencia
(International Terrestrial Reference Frame, ITRF), lo cual dicho en palabras muy
sencillas consiste en una lista de coordenadas y velocidades para una selección de
estaciones IERS (estaciones de registro continuo).
Actualmente, el ITRFyy es publicado periódicamente por el IERS, donde yy especifica
el último año cuyos datos se han usado para la formación del marco. Así por ejemplo,
ITRF97 designa el marco de coordenadas y velocidades construido en 1997 usando los
datos disponibles del IERS durante 1996.
El conjunto de puntos con coordenadas y velocidades son determinados conjuntamente
con datos de VLBI, GPS, SLR/LLR y DORIS. Su materialización se especifica en
coordenadas cartesianas ecuatoriales y en coordenadas geográficas en el elipsoide
GRS80. El Sistema Convencional de Referencia Terrestre se conoce como ITRS
(Sistema Internacional de Referencia Terrestre). Cada técnica geodésica espacial
Tema 1 - 59

Tema 1 - 60
suministra su correspondiente conjunto de datos, cuyo análisis conduce a una versión
del ITRS.
La posición de un punto situado sobre la superficie de la Tierra en una época "t" cuyas
coordenadas "x(t)" son funciones del tiempo, se expresa como:
∑Δ
+
−
+
=
x
X t
X
t
t
dt
dX
X
t
X )
(
)
(
)
( 0
0
0
donde ΔXX son las correcciones debidas a diversos efectos variables en el tiempo
(desplazamientos ocasionados por mareas terrestres, incluidos los efectos
permanentes, la carga oceánica, mareas polares, rebote postglacial, carga atmosférica,
variaciones del geocentro...), X0 y t0 son las coordenadas y el tiempo de la época
inicial. Para alcanzar precisión relativa de 10-10
, se deben determinar las correcciones
con una precisión de 1 mm.
El modelo para calcular posiciones en una época de referencia determinada dentro del
mismo ITRF es lineal, en la forma:
)
(
)
( 0
0 t
t
X
X
t
X R −
⋅
+
=
La definición del ITRS cumple las siguientes condiciones:
• Es geocéntrico, el centro de masas se define para la Tierra entera, incluyendo
océanos y atmósfera.
• La orientación fue inicialmente dada por el Bureau International de l`Heure
(BIH) en la época 1984.0.
• La evolución con el tiempo de la orientación es definida usando la condición
NNR (no-net rotation) con respecto a movimientos horizontales tectónicos sobre
la Tierra entera.
La Oficina Central del IERS tiene la sede en el Observatorio de París y más
concretamente en la Sección del Marco Terrestre, Laboratorio de Investigación en
Geodesia (LAREG) del Instituto Geográfico Nacional y se encarga de producir y publicar
los sucesivos marcos ITRF.
La estrategia general para el cálculo de ITRF consiste en combinar todas las soluciones
posibles ofrecidas por los centros de análisis de redes SLR, VLBI y GPS, quitando los

constreñimientos para luego poner los mínimos constreñimientos y en todo caso, dejar
unos vínculos localmente nada más.
Para establecer el Datum de un ITRF determinado se necesitan 4 parámetros:
• Orientación, desde el ITRF 88 hasta el ITRF92, la orientación fue
definida para que no existiese rotación entre los diferentes marcos de
referencia. La orientación del ITRF93 ha sido constreñida para que sea
consistente con los parámetros de orientación Terrestre en la época 88.0
definidos por el IERS.
• Origen, a partir del ITRS94, el origen se calculó como una media
ponderada de soluciones SLR y GPS.
• Escala, a partir del ITRF94, la escala se fijó teniendo en cuenta dos
aspectos, por un lado una media ponderada de soluciones VLBI. SLR y GPS,
y por otra parte que fuese consistente con la definición y realización del
ITRS realizada por el IUGG/IAU.
• Variaciones con el tiempo de estos parámetros.
Estos parámetros variarán de unos a otros, dependiendo de la precisión de los datos
obtenidos o de las técnicas empleadas para la obtención de los mismos. La
transformación de las coordenadas de un ITRF a otro ITRF se efectúan por medio de
matrices de transformación.
En particular, para el ITRF 2000, las soluciones que han intervenido en la definición
final una vez quitados los constreñimientos han sido: 3 de VLBI, 7 de SLR, 1 de LLR, 6
de soluciones globales de GPS, 2 de DORIS y 1 de una combinación multi-técnica. En
suma a estas, 6 soluciones GPS para densificaciones regionales que se añadirán
posteriormente al ITRF (SIRGAS, EUREF...).
¾ La escala del sistema se ha determinado a partir de las observaciones de VLBI
y SLR
¾ El origen únicamente a partir de las soluciones ponderadas de SLR, ya que se
supone la variación es lineal.
¾ La orientación ha sido calculada a partir de estaciones IERS (todas las técnicas)
con una calidad garantizada en base a:
- Observación continua al menos durante 3 años
Tema 1 - 61

- Estaciones situadas lejos de zona de contacto de placas o zonas de
deformación.
- Precisión final de velocidad mejor que 3 mm/ año.
- Residuos en las velocidades inferiores a 3 mm/año en al menos tres
soluciones diferentes.
En el siguiente gráfico se observan las 54 estaciones que cumplen todos los criterios
(azul), y en rojo, 41 estaciones que no cumplen el tercero.
Fig. 22. Estaciones ITRF que cumplen todos los criterios.
Anteriormente, los valores de orientación se sacaban minimizando las diferencias
observadas con respecto al modelo NNR-NUVEL1A, pero en esta versión esto se ha
desechado, puesto que se han encontrado inconsistencias fuertes con este modelo.
El resultado final es un fichero de las estaciones IERS con sus coordenadas y
velocidades, las cuales han de ser tenidas en cuenta en el cómputo de cualquier red
geodésica con precisiones suficientemente adecuadas para la Geodinámica. Cualquier
medida geodésica de precisión ligada a un Sistema de Referencia Terrestre tendrá que
estar referido al ITRF temporal y posteriormente al Sistema correspondiente (ITRS).
En el siguiente gráfico se ven todas las estaciones mundiales en sus diferentes técnicas
de observación para la definición del marco global y la partición de las diferentes
placas tectónicas:
Tema 1 - 62

Fig. 23. Estaciones ITRF y número de para su determinación.
Las estaciones que tienen un error medio cuadrático (RMS) de los desplazamientos
mayor que 5 cm / año o si la magnitud del vector es mayor de 10 cm / año, son
eliminadas del análisis.
6.6. El International Terrestrial Reference Frame 2000.
El ITRF2000 está dividido en dos partes:
1. Combinación principal de todas las soluciones conteniendo posiciones y velocidades
para todas las estaciones, dividido en las diferentes técnicas usadas:
• VLBI: ITRF2000_VLBI.SSC
• SLR: ITRF2000_SLR.SSC
• GPS: ITRF2000_GPS.SSC
• DORIS: ITRF2000_DORIS.SSC
2. Otras soluciones de densificaciones que serán añadidas en un segundo paso
(densificaciones regionales):
• Estaciones EUREF GPS: ITRF2000_EUR_GPS_PERM.SSC (red europea)
Tema 1 - 63

• Puntos EUREF (GPS y otros en Europa):ITRF2000_EUROPE.SSC (red europea
secundaria).
• CORS Network:ITRF2000_CORS.SSC (EEUU)
• ALASKA Network:ITRF2000_ALASKA.SSC (Alaska).
• South America Network: ITRF2000_SIRGAS.SSC (Sudamérica)
• SCAR Network:ITRF2000_SCAR.SSC (Antártida)
• REGAL Network:ITRF2000_REGAL.SSC (Alpes)
• RGP Network:ITRF2000_RGP.SSC (Francia)
(lareg.ensg.ign.fr en directorio pub/itrf/itrf2000):
La estrategia de análisis seguida para la generación de las soluciones ITRF se
siguen a partir de los siguientes pasos:
1. Eliminación de constreñimientos a las distintas soluciones aportadas.
2. Aplicación de mínimos constreñimientos globales.
3. Propagación para cada solución de las posiciones de la EMPV (Épocas de
Mínima Varianza en la Posición).
4. Combinación de todas las soluciones juntas.
En el caso del GPS, los Centros de Análisis que han aportado las soluciones finales al
ITRF2000 han sido:
Tabla 4. Centros de Análisis del IGS para el ITRF2000.
CENTRO ACRÓNIMO TIEMPO Nº ESTACIONES
Center for Orbit Determination
in Europe
(CODE) 93-00 160
GeoForschungsZentrum
Potsdam
(GFZ) 93-00 98
International GPS Service (IGS) 96 -00 179
Jet Propulsion Laboratory (JPL) 91-99 112
University of Newcastle upon
Tyne
(NCL) 95-99 90
NOAA, National Geodetic
Survey
(NOAA) 94-00 165
Tema 1 - 64

Tabla 5. Fragmento del ITRF2000 para SIRGAS.
DOMES NB. SITE NAME TECH. ID. X/Vx Y/Vy Z/Vz Sigmas SOLN
-----------------------m/m/y----------------------
--------------------------------------------------------------------------------------------------
30602M001 ASCENSION GPS ASC1 6118526.077 -1572344.698 -876451.166 .004 .002 .001
30602M001 -.0025 -.0048 .0102 .0015 .0007 .0004
41602M001 FORTALEZA GPS FORT 4985386.627 -3954998.587 -428426.482 .002 .002 .001
41602M001 -.0013 -.0044 .0121 .0004 .0004 .0003
41901M001 BOGOTA GPS BOGT 1744399.096 -6116037.813 512731.625 .003 .005 .001
41901M001 -.0057 .0277 .0112 .0015 .0027 .0006
42005M001 SANTA CRUZ GPS GALA -33796.116 -6377522.655 -82120.899 .002 .006 .001
42005M001 .0516 -.0011 .0123 .0009 .0022 .0005
42402M001 MARACAIBO GPS MARA 1976117.057 -5948895.194 1173592.127 .029 .073 .019
42402M001 .0141 .0023 .0080 .0102 .0255 .0065
Fig. 24. Velocidades ITRF2000.
6.7. Transformaciones entre Marcos de Referencia.
La relación entre el Marco Internacional de Referencia Celeste (ICRF) y el Marco
Internacional de Referencia Terrestre (ITRF) se describe mediante la fórmula:
ITRF
ICRF X
Y
X
T
N
P
X ⋅
⋅
⋅
⋅
⋅
=
donde P es la matriz de precesión, N la de nutación, T es la matriz que expresa la
rotación para un tiempo sidéreo S y XY son las matrices de transformación desde el
Tema 1 - 65

marco terrestre al marco conectado con el polo instantáneo de efemérides. XICRF es el
vector relacionado con el ICRF y XITRF es el vector relacionado con los ejes ITRF.
Lo que más nos interesa es la transformación entre dos marcos de referencia
terrestres. Bajo la hipótesis de un Sistema de Referencia Terrestre ideal definido como
un triedro cercano a la Tierra y rotando con ella, la transformación general de
coordenadas para cualquier punto de la Tierra desde TRS1 a TRS2 viene dada por una
similitud en tres dimensiones en la forma (Boucher, 1986):
1
2
,
1
2
,
1
2
,
1
2 X
R
T
X ⋅
⋅
+
= λ
donde T es un vector de translación, λ un factor de escala y R una matriz de rotación.
Como consecuencia de la aplicación conceptual anterior, el International Earth Rotation
Service (IERS) usa las fórmulas linealizadas. La transformación estándar entre dos
sistemas de referencia es una similitud euclidiana de siete parámetros: tres
traslaciones, un factor de escala y tres rotaciones, pero también sus derivadas con el
tiempo ( ) en un sentido estricto.
R
D
T

,
,
Utilizando la forma lineal:
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎣
⎡
⋅
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎣
⎡
−
−
−
+
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎣
⎡
+
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎣
⎡
=
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎣
⎡
1
1
1
1
1
1
2
2
2
1
2
1
3
2
3
z
y
x
D
R
R
R
D
R
R
R
D
T
T
T
z
y
x
z
y
x
z
y
x
donde T son las translaciones, R las rotaciones y D el factor de escala.
Si se consideran las derivadas con el tiempo de los parámetros anteriores y derivamos
la ecuación anterior en la forma:
1
1
1
2 X
R
X
D
T
X
X ⋅
+
⋅
+
+
= con λ=1+D y R=(I+R)
queda:
1
1
1
1
1
2 X
R
X
R
X
D
X
D
T
X
X ⋅
+
⋅
+
⋅
+
⋅
+
+
=

D y R son del orden de 10-5
y X
unos 10 cm por año, los términos 1
X
D
⋅ y 1
X
R
⋅
representan 0.1 mm en 100 años, luego se pueden despreciar, con lo que la ecuación
rigurosa anterior para la transformación de 14 parámetros queda:
1
1
1
2 X
R
X
D
T
X
X ⋅
+
⋅
+
+
=

Tema 1 - 66

Tabla 6. Parámetros de transformación entre diferentes realizaciones del ITRF.
6.8. Organismos Internacionales.
6.8.1. International Earth Rotation and Reference Systems Service (IERS).
El International Earth Rotation Service fue creado en 1987 por
el IAG con el objetivo principal de servir a la comunidad
astronómica, geodésica y geofísica:
• El International Celestial Reference System (ICRS) y su
realización, el International Celestial Reference Frame (ICRF).
• El International Terrestrial Reference System (ITRS) y su
realización, el International Terrestrial Reference Frame (ITRF).
• Los Earth orientation parameters (EOP) requeridos para el estudio de la variación
rotacional terrestre y la transformación entre el ICRF y el ITRF. (Coordenadas del
polo.
Tema 1 - 67

• Datos geofísicos para interpretar las variaciones temporal – espaciales en el ICRF,
ITRF o parámetros de orientación de la tierra, así como modelar sus variaciones.
• Estándares, constantes y modelos que sean adoptados internacionalmente.
http://www.iers.org
6.8.2. International GPS Service (IGS).
El IGS colecta, archiva y distribuye las observaciones GPS de la suficiente precisión
para satisfacer los objetivos de un amplio rango de aplicaciones y experimentos. Estos
datos en bruto de GPS son usados por el IGS para generar :
- Efemérides de satélites GPS.
- Parámetros de rotación de la Tierra (EOP).
- Coordenadas y velocidades de las estaciones IGS.
- Información del estado de reloj de los satélites GPS y de las estaciones IGS.
En particular las precisiones de los productos del IGS son suficientes y de hecho
aprovechados para la mejora y extensión del International Terrestrial Reference Frame
(ITRF), la monitorización de las deformaciones de la Tierra sólida, la monitorización y
las variaciones de la Tierra líquida, determinación de órbitas de satélites no GPS,
monitorización de la ionosfera y la recolección de datos informativos de la troposfera.
Operacionalmente el IGS se divide fundamentalmente en:
¾ Red de Estaciones Permanentes - Centros Operacionales
¾ Centros de Datos
¾ Centro Coordinador de Centros Análisis
¾ Coordinador de ITRF
¾ Grupos de Trabajo y de Proyectos Especiales
¾ Bureau Central
¾ Dirección y Coordinación General
¾ Centros de Análisis
o Centros de Análisis Diarios
o Centros Asociados productores de coordenadas, velocidades, información
troposférica, ionosférica, combina soluciones globales...
Tema 1 - 68

A día de hoy el la red GPS del IGS consta de un total de 273 estaciones repartidas por
todo el globo.
La densificación del IGS para Europa desde el punto de vista del GPS y el Marco de
Referencia es EUREF (EUropean REference Frame), con unos objetivos análogos a los
del IGS, estructura y funcionamiento completamente comunes.
http://igscb.jpl.nasa.gov
Fig. 25. Estaciones Permanentes GPS del IGS.
6.8.3. Sistema de Referencia Geocéntrico para América (SIRGAS).
El proyecto SIRGAS fue establecido durante la Conferencia Internacional de la
Definición de un Datum Geocéntrico para Sudamérica, en octubre de 1993 en
Asunción, Paraguay, por invitación de la Asociación Internacional de Geodesia (IAG), el
Instituto Panamericano de Geografía e Historia (PAIGH) y la National Imagery and
Mapping Agency (NIMA), en aquel entonces, United States Defense Mapping Agency
(DMA). Representatantes de cada entidad y la mayoría de países de Sudamérica
participaron y crearon el Comité SIRGAS y el Consejo Científico.
Tema 1 - 69

SIRGAS es, por tanto, la densificación del ITRF en América, conformado actualmente
por una red de más de 180 estaciones (aunque sólo algunas como estaciones
permanentes), cuya distribución es regularmente homogénea sobre el continente.
Los objetivos que se establecieron para el proyecto fueron:
9 Definir un Sistema de Referencia para Sudamérica.
9 Establecer y mantener una red de referencia para establecer y definir un datum
geocéntrico.
Estos objetivos del sistema de referencia para el continente fueron asumidos en la
Conferencia de Asunción con los condicionantes:
9 Sistema de Referencia SIRGAS: International Terrestrial Reference Frame
(ITRF).
9 Datum geocéntrico: ejes de coordenadas datos por el Sistema de Referencia
SIRGAS y parámetros del elipsoide GRS80, orientado según los ejes
coordenados equivalente al ITRF (SIRGAS 1997).
Se establecieron tres grupos de trabajo con vistas a cumplir los objetivos:
¾ Grupo I: Sistemas de Referencia
¾ Grupo II: Datum Geocéntrico
¾ Grupo III: Datum Vertical
Al principio SIRGAS fue restringido a América del Sur y el Sistema de Referencia fue
realizado con 58 estaciones observadas por una campaña de GPS en 1995. Las
coordenadas finales estaban en ITRF94 (época 1995.4) y los ajustes fueron realizados
por el Deutsches Geodädisches Forschungsinstitut (DGFI) y el NIMA, siendo aprobados
y presentados en la Asamblea del IAG en Rio de Janeiro, 1997.
El datum geocéntrico está realizado sobre las redes observadas con GPS por los
diferentes países y conectados con las redes de triangulación existentes,
densificaciones de SIRGAS.
El datum vertical está definido por estaciones mareógraficas (para cada país) y
observadas con GPS junto con otras señales con nivelación. Para este propósito se
llevó a cabo una campaña en mayo de 2000 para observar todos los mareógrafos y
algunas estaciones que conforman el marco de referencia vertical. Esta campaña
también fue la repetición de las mediciones hechas sobre el Marco de Referencia
Tema 1 - 70

Geocéntrico SIRGAS e incluyó el resto de América, desde los países de América del Sur
y nuevos puntos de América Central, Estados Unidos y Canadá (SIRGAS 2002). El
resultado fue una red conformada por 184 estaciones con coordenadas calculadas en el
ITRF2000, época 2000.4.
Una red de estaciones permanentes GPS proporciona también datos al IGS y son
procesadas semanalmente por el RNAAC-SIRGAS, que es una densificación regional de
la red global IGS. Incluye preferentemente estaciones de observación continua, pero
también incluye la ocupación periódica de las estaciones pasivas de SIRGAS,
proporcionando campos de velocidades de las estaciones.
Tema 1 - 71

Fig. 26. Estaciones GPS SIRGAS.
La red GPS permanente está compuesta por más de 40 puntos en Sudamérica y los
datos son procesados por el Deutsches Geodätisches Forschungsinstitut (DGFI) como
Centro de Procesamiento Regional o RNAAC (Regional Network Associate Analysis
Center) del IGS. Se calculan velocidades a partir de las estaciones permanentes, las
campañas anteriormente citadas de 1995 y 2000 y los proyectos geodinámicos de
reobservación, como por ejemplo CAP (Central Andes GPS Project), SAGA (South
America Geodynamics Activity), SNAPP (South America – Nazca Plate Motion Project) y
CASA (Central And South America GPS Geodynamics Project. El último campo de
Tema 1 - 72

velocidades para el continente sudamericano se muestra en la siguiente figura (Drewes
Heidbach, 2003).
Fig. 27. Campo de velocidades SIRGAS.
La ONU, en su Séptima Conferencia Cartográfica para las Américas (Nueva York, 2001)
recomendó la adopción de SIRGAS como sistema de referencia oficial para todos los países de
América (ver http://www.ibge.gov.br/home/geografia/geodesico/sirgas/principal.html).
Tema 1 - 73

Tema 1 - 74
7. Direcciones de interés.
Antes de continuar el curso conviene citar unas referencias básicas en la web acerca de
los temas explicados y otras utilidades relacionadas.
• International GPS Service: http://igscb.jpl.nasa.gov
• UNAVCO: http://www.unavco.org
• International Earth Rotation Service: http://www.iers.org
• National Geodetic Survey (EEUU): http://www.ngs.noaa.gov
• ITRF web site: http://itrf.ensg.ign.fr
• International Union of Geodesy and Geophysics (IUGG): http://www.iugg.org
• Scripps Orbit and Permanent Array Center (SOPAC): http://sopac.ucsd.edu
• European Reference Frame (EUREF): http://www.epncb.oma.be
• BKG Data Center: http://igs.ifag.de

Tema1 geodesia

Más contenido relacionado

La actualidad más candente (20)

Similar a Tema1 geodesia (20)

Último (20)

Tema1 geodesia