En geodesia y navegación , un arco de meridiano es la curva entre dos puntos de la superficie terrestre que tienen la misma longitud . El término puede referirse a un segmento del meridiano o a su longitud .
El propósito de medir los arcos de meridianos es determinar una figura de la Tierra . Se pueden usar una o más mediciones de arcos de 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 de meridianos en varias latitudes a lo largo de muchos meridianos de todo el mundo se pueden combinar para aproximar un elipsoide geocéntrico destinado a adaptarse al mundo entero.
Las primeras determinaciones del tamaño de una Tierra esférica requirieron un solo arco. Los trabajos topográficos precisos que comenzaron en el siglo XIX requirieron varias mediciones de arco en la región donde se iba a realizar el estudio, lo que llevó a una proliferación de elipsoides de referencia en todo el mundo. Las últimas determinaciones 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 registran en Grecia en el siglo IV a. C., y por 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 año 240 a.C. Estimó que el meridiano tiene una longitud de 252.000 estadios , con un error sobre el valor real entre -2,4% y +0,8% (suponiendo 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. Posidonio utilizó un método similar unos 150 años después, y en 827 se calcularon resultados ligeramente mejores mediante el método de medición del arco , [2] atribuido al califa Al-Ma'mun . [ cita necesaria ]
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" generalmente se omiten. Un elipsoide que no es un elipsoide de revolución se llama elipsoide triaxial. Esferoide y elipsoide se usan indistintamente en este artículo, y se implica achatado 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 acumulaba evidencia 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 , Guayana Francesa y descubrió que perdía 2+1 ⁄ 2 minutos por día en comparación con su tarifa en París . [3] [4] Esto indicó que la aceleración de la gravedad era menor en Cayena que en París. Se empezaron a llevar gravímetros de péndulo en viajes a partes remotas del mundo, y poco a poco se descubrió que la gravedad aumenta suavemente al aumentar 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 una prueba de que la Tierra era un esferoide achatado de aplanamiento igual a 1/230 . [5] Esto fue cuestionado por algunos, pero no todos, los científicos franceses. Giovanni Domenico Cassini y su hijo Jacques Cassini ampliaron un arco meridiano de Jean Picard a un arco más largo durante el período 1684-1718. [6] El arco se midió 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ó determinar la forma general. Los resultados indicaron que la Tierra era un esferoide alargado (con un radio ecuatorial menor que el radio polar). Para resolver la cuestión, la Academia Francesa de Ciencias (1735) emprendió expediciones al 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 mejor manera de modelar la Tierra era mediante un esferoide achatado que sustentara 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 por cuatro determinaciones intermedias de latitud. Combinando las mediciones con las del arco de Perú, se determinaron los parámetros de forma del elipsoide y se calculó la distancia entre el Ecuador y el polo a lo largo del Meridiano de París como5.130.762 toises según lo especificado por la barra estándar de toises de París. Definiendo esta distancia como exactamente10.000.000 m llevaron a la construcción de una nueva barra de metros estándar como0,513 0762 toesas. [6] : 22
En el siglo XIX, muchos astrónomos y geodesistas 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] Se proporciona una lista completa de elipsoides en Elipsoide terrestre .
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 en exactamente 1.852 metros. Sin embargo, a 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 navegantes diurnos : "1 (minuto) de latitud = 1 milla marina", seguido de "Para la mayoría de los fines prácticos, la distancia se mide a partir de la escala de latitud, suponiendo que un minuto de latitud equivale a un minuto náutico". milla". [8]
En una esfera, la longitud del arco del 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 de la resta de dos distancias de meridianos , 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, particularmente la proyección transversal de Mercator .
Los principales parámetros elipsoidales son a , b , f , pero en el trabajo teórico es útil definir parámetros adicionales, particularmente 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 del 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 = mi 2/1 - mi 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 del meridiano completo (incluido el antimeridiano). Por lo tanto, los rangos de φ , β y la latitud rectificadora μ no están restringidos.
La integral anterior está relacionada con un caso especial de 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 se puede escribir en términos de integrales elípticas incompletas de segundo tipo (consulte 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 informáticos de álgebra como Mathematica [12] y Maxima. [13]
La integral anterior se puede expresar como una serie truncada infinita expandiendo 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 dedujo 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 una convergencia considerablemente más rápida expandiendo en términos del tercer aplanamiento n en lugar de la excentricidad. Están relacionados por
En 1837, Friedrich Bessel obtuvo una de esas series, [17] que Helmert puso en una forma más simple , [18] [19]
con
Porque n cambia de signo cuando se intercambian a y b , 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 desaparecen.
La serie se puede expresar con a o b como factor inicial escribiendo, por ejemplo,
y expandir 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 de la Agencia Nacional de Inteligencia Geoespacial [20] y el Ordnance Survey de Gran Bretaña . [21]
En 1825, Bessel [22] dedujo 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 usar para escribir la longitud del arco en términos de la latitud geodésica como
La serie anterior, de octavo orden en excentricidad o de cuarto orden en tercer aplanamiento, proporciona una precisión milimétrica. Con la ayuda de sistemas de álgebra simbólica, se pueden extender fácilmente al sexto orden en el tercer aplanamiento, lo que proporciona precisión doble total para aplicaciones terrestres.
Delambre [15] y Bessel [22] escribieron sus series de una forma que les permite generalizarse en un orden arbitrario. Los coeficientes de la serie de Bessel se pueden expresar de forma particularmente sencilla
dónde
yk !! es el factorial doble , extendido a valores negativos mediante la relación de recursividad: (−1)!! = 1 y (−3)!! = −1 .
Los coeficientes de la serie de Helmert pueden expresarse de manera similar en general mediante
Este resultado fue conjeturado por Friedrich Helmert [23] y demostrado por Kazushige Kawase. [24]
El factor adicional (1 − 2 k )(1 + 2 k ) se origina en la expansión adicional 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 convenientemente utilizando la suma de Clenshaw . Este método evita el cálculo de la mayoría de las funciones trigonométricas y permite sumar las series de forma rápida y precisa. La técnica también se puede utilizar para evaluar la diferencia m ( φ 1 ) − m ( φ 2 ) manteniendo 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 de la Tierra , es
Formó parte de la definición histórica del metro y de la milla náutica , y se utilizó en la definición del hebdomómetro .
El cuarto de meridiano se puede expresar en términos de la integral elíptica completa de segundo tipo ,
¿Dónde están la primera y la segunda excentricidad ?
El cuarto de meridiano también viene dado por la siguiente serie generalizada:
(Para conocer la fórmula de c 0 , consulte la sección #Serie generalizada más arriba). Este resultado fue obtenido por primera vez por James Ivory . [25]
La expresión numérica para el cuarto de meridiano en el elipsoide WGS84 es
La circunferencia polar de la Tierra es simplemente cuatro veces un cuarto de meridiano:
El perímetro de una elipse meridiana también se puede reescribir en forma de perímetro de un círculo rectificador, C p = 2π M r . Por tanto, el radio rectificador de la Tierra es:
Se puede evaluar como6 367 449 , 146 m .
En algunos problemas, debemos poder resolver el problema inverso: dado m , determine φ . Esto puede resolverse mediante el método de Newton , iterando
hasta la convergencia. Una estimación inicial adecuada viene dada por φ 0 = μ donde
es la latitud rectificadora . Tenga en cuenta que no es necesario diferenciar la serie para m ( φ ) , ya que en su lugar se puede utilizar 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 β puede revertirse para dar [28]
dónde
Adrien-Marie Legendre demostró que la distancia a lo largo de una geodésica sobre 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 juega 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 al sexto orden las proporciona Charles Karney, [31] Ecs. (17) y (21), con ε desempeñando el papel de n y τ desempeñando el papel de μ .