En geodesia y navegación , un arco meridiano es la curva entre dos puntos de la superficie terrestre que tienen la misma longitud . El término puede referirse tanto a un segmento del meridiano como a su longitud .
El propósito de medir los arcos meridianos es determinar una figura de la Tierra . Se pueden utilizar una o más mediciones de arcos meridianos para inferir la forma del elipsoide de referencia que mejor se aproxima al geoide en la región de las mediciones. Las mediciones de arcos meridianos en varias latitudes a lo largo de muchos meridianos alrededor del mundo se pueden combinar para aproximarse a un elipsoide geocéntrico destinado a adaptarse a todo el mundo.
Las primeras determinaciones del tamaño de una Tierra esférica requerían un único arco. Los trabajos de topografía precisos que comenzaron en el siglo XIX requerían varias mediciones de arco en la región en la que se iba a realizar la topografía, lo que dio lugar a una proliferación de elipsoides de referencia en todo el mundo. Las determinaciones más recientes utilizan mediciones astrogeodésicas y métodos de geodesia satelital para determinar elipsoides de referencia, especialmente los elipsoides geocéntricos que ahora se utilizan para sistemas de coordenadas globales como WGS 84 (ver expresiones numéricas).
Las primeras estimaciones del tamaño de la Tierra se registraron en Grecia en el siglo IV a. C., y por parte de los eruditos de la Casa de la Sabiduría del califa en Bagdad en el siglo IX. El primer valor realista fue calculado por el científico alejandrino Eratóstenes alrededor del 240 a. C. Estimó que el meridiano tiene una longitud de 252.000 estadios , con un error en el valor real de entre el -2,4% y el +0,8% (asumiendo un valor para el estadio entre 155 y 160 metros). [1] Eratóstenes describió su técnica en un libro titulado Sobre la medida de la Tierra , que no se ha conservado. Un método similar fue utilizado por Posidonio unos 150 años después, y se calcularon resultados ligeramente mejores en 827 mediante el método de medición del arco , [2] atribuido al califa Al-Ma'mun . [ cita requerida ]
La literatura antigua utiliza el término esferoide achatado para describir una esfera "aplastada en los polos". La literatura moderna utiliza el término elipsoide de revolución en lugar de esferoide , aunque las palabras calificativas "de revolución" suelen omitirse. Un elipsoide que no es un elipsoide de revolución se denomina elipsoide triaxial. Esferoide y elipsoide se utilizan indistintamente en este artículo, con oblato implícito si no se indica.
Aunque se sabía desde la antigüedad clásica que la Tierra era esférica , en el siglo XVII se acumulaban pruebas de que no era una esfera perfecta. En 1672, Jean Richer encontró la primera evidencia de que la gravedad no era constante sobre la Tierra (como lo sería si la Tierra fuera una esfera); llevó un reloj de péndulo a Cayena , en la Guayana Francesa , y descubrió que perdía 2+1 ⁄ 2 minutos por día en comparación con su tasa en París . [3] [4] Esto indicó que la aceleración de la gravedad era menor en Cayena que en París. Los gravímetros de péndulo comenzaron a usarse en viajes a partes remotas del mundo y poco a poco se descubrió que la gravedad aumenta suavemente con el aumento de la latitud , siendo la aceleración gravitacional aproximadamente un 0,5% mayor en los polos geográficos que en el Ecuador .
En 1687, Isaac Newton había publicado en los Principia como prueba que la Tierra era un esferoide achatado de aplanamiento igual a 1/230 . [5] Esto fue disputado por algunos, pero no todos, los científicos franceses. Un arco meridiano de Jean Picard fue extendido a un arco más largo por Giovanni Domenico Cassini y su hijo Jacques Cassini durante el período 1684-1718. [6] El arco fue medido con al menos tres determinaciones de latitud, por lo que pudieron deducir curvaturas medias para las mitades norte y sur del arco, lo que permitió una determinación de la forma general. Los resultados indicaron que la Tierra era un esferoide alargado (con un radio ecuatorial menor que el radio polar). Para resolver el problema, la Academia Francesa de Ciencias (1735) emprendió expediciones a Perú ( Bouguer , Louis Godin , de La Condamine , Antonio de Ulloa , Jorge Juan ) y a Laponia ( Maupertuis , Clairaut , Camus , Le Monnier , Abbe Outhier , Anders Celsius ). Las mediciones resultantes en latitudes ecuatoriales y polares confirmaron que la Tierra se modelaba mejor mediante un esferoide achatado, lo que apoyaba a Newton. [6] Sin embargo, en 1743, el teorema de Clairaut había suplantado por completo el enfoque de Newton.
A finales de siglo, Jean Baptiste Joseph Delambre había vuelto a medir y ampliado el arco francés desde Dunkerque hasta el mar Mediterráneo (el arco meridiano de Delambre y Méchain ). Estaba dividido en cinco partes mediante cuatro determinaciones intermedias de latitud. Al combinar las mediciones junto con las del arco de Perú, se determinaron los parámetros de forma elipsoide y se calculó la distancia entre el Ecuador y el polo a lo largo del meridiano de París como5 130 762 toesas según lo especificado por la barra de toesas estándar en París. Definiendo esta distancia exactamente como10.000.000 m llevaron a la construcción de una nueva barra de metro estándar como0,513 0762 toesas. [6] : 22
En el siglo XIX, muchos astrónomos y geodestas se dedicaron a estudios detallados de la curvatura de la Tierra a lo largo de diferentes arcos meridianos. Los análisis dieron como resultado una gran cantidad de elipsoides modelo, como Plessis 1817, Airy 1830, Bessel 1841 , Everest 1830 y Clarke 1866. [ 7] En el apartado Elipsoide terrestre se ofrece una lista completa de elipsoides .
Históricamente, una milla náutica se definía como la longitud de un minuto de arco a lo largo de un meridiano de una Tierra esférica. Un modelo elipsoide conduce a una variación de la milla náutica con la latitud. Esto se resolvió definiendo la milla náutica como exactamente 1.852 metros. Sin embargo, para todos los efectos prácticos, las distancias se miden a partir de la escala de latitud de las cartas. Como dice la Royal Yachting Association en su manual para patrones de navegación diurna : "1 (minuto) de latitud = 1 milla náutica", seguido de "Para la mayoría de los fines prácticos, la distancia se mide a partir de la escala de latitud, asumiendo que un minuto de latitud equivale a una milla náutica". [8]
En una esfera, la longitud del arco meridiano es simplemente la longitud del arco circular . En un elipsoide de revolución, para arcos meridianos cortos, su longitud se puede aproximar utilizando el radio de curvatura meridional de la Tierra y la formulación del arco circular. Para arcos más largos, la longitud se obtiene restando dos distancias meridianas , la distancia desde el ecuador hasta un punto en una latitud φ . Este es un problema importante en la teoría de las proyecciones cartográficas, en particular la proyección transversal de Mercator .
Los parámetros elipsoidales principales son a , b , f , pero en el trabajo teórico resulta útil definir parámetros adicionales, en particular la excentricidad , e , y el tercer aplanamiento, n . Sólo dos de estos parámetros son independientes y existen muchas relaciones entre ellos:
Se puede demostrar que el radio de curvatura del meridiano es igual a: [9] [10]
La longitud de arco de un elemento infinitesimal del meridiano es dm = M ( φ ) dφ (con φ en radianes). Por lo tanto, la distancia del meridiano desde el ecuador hasta la latitud φ es
La fórmula de la distancia es más simple cuando se escribe en términos de latitud paramétrica ,
donde tan β = (1 − f )tan φ y e ′ 2 = y 2/1 − y 2 .
Aunque la latitud normalmente se limita al rango [− π/2 , π/2 ] , todas las fórmulas dadas aquí se aplican para medir la distancia alrededor de la elipse meridiana completa (incluido el antimeridiano). Por lo tanto, los rangos de φ , β y la latitud rectificadora μ no tienen restricciones.
La integral anterior está relacionada con un caso especial de una integral elíptica incompleta de tercer tipo . En la notación del manual en línea del NIST [11] (Sección 19.2(ii)),
También puede escribirse en términos de integrales elípticas incompletas de segundo tipo (véase la Sección 19.6(iv) del manual del NIST).
El cálculo (con precisión arbitraria) de las integrales elípticas y las aproximaciones también se analizan en el manual del NIST. Estas funciones también se implementan en programas de álgebra computacional como Mathematica [12] y Maxima. [13]
La integral anterior puede expresarse como una serie truncada infinita desarrollando el integrando en una serie de Taylor, realizando las integrales resultantes término por término y expresando el resultado como una serie trigonométrica. En 1755, Leonhard Euler derivó una expansión en la tercera excentricidad al cuadrado. [14]
Delambre en 1799 [15] derivó una expansión ampliamente utilizada en e 2 ,
dónde
Richard Rapp ofrece una derivación detallada de este resultado. [16]
Se pueden obtener series con convergencia considerablemente más rápida expandiéndolas en términos del tercer aplanamiento n en lugar de la excentricidad. Están relacionadas por
En 1837, Friedrich Bessel obtuvo una de esas series, [17] que Helmert simplificó , [18] [19]
con
Porque n cambia de signo cuando a y b se intercambian, y porque el factor inicial 1/2 ( a + b ) es constante bajo este intercambio, la mitad de los términos en las expansiones de H 2 k se desvanecen.
La serie se puede expresar con a o b como factor inicial escribiendo, por ejemplo,
y expandiendo el resultado como una serie en n . Aunque esto da como resultado series que convergen más lentamente, dichas series se utilizan en la especificación para la proyección transversal de Mercator por parte de la Agencia Nacional de Inteligencia Geoespacial [20] y el Servicio Geoespacial de Gran Bretaña [21] .
En 1825, Bessel [22] derivó una expansión de la distancia del meridiano en términos de la latitud paramétrica β en relación con su trabajo sobre geodésicas .
con
Debido a que esta serie proporciona una expansión para la integral elíptica de segundo tipo, se puede utilizar para escribir la longitud del arco en términos de la latitud geodésica como
Las series anteriores, hasta el octavo orden en excentricidad o el cuarto orden en tercer aplanamiento, proporcionan una precisión milimétrica. Con la ayuda de sistemas de álgebra simbólica, se pueden ampliar fácilmente hasta el sexto orden en tercer aplanamiento, lo que proporciona una precisión doble completa para aplicaciones terrestres.
Delambre [15] y Bessel [22] escribieron sus series de una forma que permite generalizarlas a un orden arbitrario. Los coeficientes de la serie de Bessel pueden expresarse de forma particularmente sencilla.
dónde
y k !! es el factorial doble , extendido a valores negativos mediante la relación de recursión: (−1)!! = 1 y (−3)!! = −1 .
Los coeficientes de la serie de Helmert pueden expresarse de manera similar y general mediante
Este resultado fue conjeturado por Friedrich Helmert [23] y demostrado por Kazushige Kawase. [24]
El factor extra (1 − 2 k )(1 + 2 k ) se origina a partir de la expansión adicional de que aparece en la fórmula anterior y da como resultado una convergencia más pobre de la serie en términos de φ en comparación con la de β .
La serie trigonométrica dada anteriormente se puede evaluar de manera conveniente utilizando la suma de Clenshaw . Este método evita el cálculo de la mayoría de las funciones trigonométricas y permite sumar la serie de manera rápida y precisa. La técnica también se puede utilizar para evaluar la diferencia m ( φ 1 ) − m ( φ 2 ) manteniendo al mismo tiempo una alta precisión relativa.
Sustituyendo los valores del semieje mayor y la excentricidad del elipsoide WGS84 se obtiene
donde φ ( ° ) = φ/1° es φ expresado en grados (y de manera similar para β ( ° ) ).
En el elipsoide, la distancia exacta entre los paralelos en φ 1 y φ 2 es m ( φ 1 ) − m ( φ 2 ) . Para WGS84, una expresión aproximada para la distancia Δ m entre los dos paralelos a ±0,5° del círculo en la latitud φ viene dada por
La distancia del ecuador al polo, el cuarto meridiano (análogo al cuarto de círculo ), también conocido como cuadrante terrestre , es
Fue parte de la definición histórica del metro y de la milla náutica , y se utilizó en la definición del hebdómetro .
El cuarto meridiano se puede expresar en términos de la integral elíptica completa de segundo tipo ,
¿Dónde están la primera y segunda excentricidades ?
El meridiano del cuarto de hora también viene dado por la siguiente serie generalizada:
(Para la fórmula de c 0 , consulte la sección #Series generalizadas más arriba). Este resultado fue obtenido por primera vez por James Ivory . [25]
La expresión numérica para el cuarto meridiano en el elipsoide WGS84 es
La circunferencia de la Tierra polar es simplemente cuatro veces el cuarto del meridiano:
El perímetro de una elipse meridiana también se puede reescribir en forma de perímetro de círculo rectificador, C p = 2π M r . Por lo tanto, el radio rectificador de la Tierra es:
Se puede evaluar como6 367 449 .146 metros .
En algunos problemas, necesitamos poder resolver el problema inverso: dado m , determinar φ . Esto se puede resolver mediante el método de Newton , iterando
hasta la convergencia. Una estimación inicial adecuada viene dada por φ 0 = μ donde
es la latitud rectificadora . Nótese que no es necesario diferenciar la serie para m ( φ ) , ya que se puede utilizar en su lugar la fórmula para el radio de curvatura del meridiano M ( φ ) .
Alternativamente, la serie de Helmert para la distancia del meridiano se puede revertir para dar [26] [27]
dónde
De manera similar, la serie de Bessel para m en términos de β se puede revertir para dar [28]
dónde
Adrien-Marie Legendre demostró que la distancia a lo largo de una geodésica en un esferoide es la misma que la distancia a lo largo del perímetro de una elipse. [29] Por esta razón, la expresión para m en términos de β y su inversa dada anteriormente juegan un papel clave en la solución del problema geodésico con m reemplazado por s , la distancia a lo largo de la geodésica, y β reemplazado por σ , la longitud del arco en la esfera auxiliar. [22] [30] Las series requeridas extendidas a sexto orden están dadas por Charles Karney, [31] Eqs. (17) y (21), con ε jugando el papel de n y τ jugando el papel de μ .