Distancia al horizonte

Supongamos que la Tierra no tiene atmósfera, la luz se propaga a lo largo de una línea recta. El observador está en A a una altura h sobre la superficie de la Tierra. Trazamos la línea AB tangente a la superficie terrestre en B y que pasa por A. En el triángulo rectángulo OBA

( R+h ) 2 = R 2 + d 2 d 2 =h( 2R+h ) d 2Rh ,hR

Sabiendo que el radio de la Tierra R=6371 km.

Corrección debido a la refracción de la atmósfera

El problema no es sencillo si se tiene en cuenta el hecho de que la atmósfera es un medio cuyo índice de refracción disminuye con la altura. La trayectoria del rayo de luz ya no es rectilínea sino una curva cóncava hacia la superficie de la Tierra.

En la página titulada 'Refracción en la atmósfera', demostramos una importante relación de la óptica atmosférica

n1r1·sini1=n2r2·sini2=constante

La línea de color rojo, es la trayectoria seguida por el rayo de luz. Los ángulos que forman las tangentes a dicha trayectoria en A y B con la dirección radial son los ángulos de incidencia i1 e i2. Los índices de refracción de la atmósfera en las posiciones A y B son n1 y n2, respectivamente

En la situación que estudiamos, el rayo de luz parte de un punto A en la superficie de la Tierra y llega al observador situado en B a una altura h

Como apreciamos en la figura

Utilizamos la aproximación (1+x)-1≈1-x, para x<<1.

>> syms x;
>> taylor(1/(1+x),x) 
ans =- x^5 + x^4 - x^3 + x^2 - x + 1 

n(R+h)sin i 2 = n 0 Rsin i 1 sin i 2 = n 0 R ( n 0 +Δn )( R+h ) = ( 1+ Δn n 0 ) 1 ( 1+ h R ) 1 sin i 2 1 Δn n 0 h R

Los ángulo i2 y δ son complementarios, sinδ=cosi2. Teniendo en cuenta que el ángulo δ es pequeño, sinδδ

δ 1 sin 2 i 2 2 Δn n 0 +2 h R = 2 h R 1+ R h Δn n 0 2 h R ( 1+ R 2h Δn n 0 )

Utilizamos la aproximación (1+x)1/2≈1+x/2, para x<<1.

>> syms x;
>> taylor(sqrt(1+x),x)
ans =(7*x^5)/256 - (5*x^4)/128 + x^3/16 - x^2/8 + x/2 + 1

Cuando el efecto de la atmósfera es despreciable

δ 2 h R

tal como puede deducirse a partir de la primera figura, suponiendo que h<<R

Indice de refracción

En la página titulada 'Refracción en la atmósfera', proporcionamos una expresión para el valor del índice de refracción n(r) con la distancia radial

n(r)=1+kexp( rR H ) n(z)=1+kexp( αz )

donde z es la altura sobre el nivel del mar

Los valores de los parámetros α y k que proporciona el artículo que se menciona en las referencias son k=2.77·10-4, α=0.105 km-1

Δn(z)=n(z) n 0 =1+kexp( αz )(1+k)=k( 1exp( αz ) )

Para alturas z<1 km, αz es pequeño, utilizamos la aproximación, 1-exp(-x)≈x, para x<<1

>> syms x;
>> taylor(1-exp(-x),x)
ans =x^5/120 - x^4/24 + x^3/6 - x^2/2 + x

δ 2 h R ( 1 R 2h k( 1exp( αh ) ) 1+k ) 2 h R ( 1 R 2 kα ) δ(10.092) 2 h R

El factor de corrección para el ángulo δ es 1-0.092 cuando tenemos en cuenta la influencia de la atmósfera

Distancia al horizonte

Consideremos el rayo de luz que parte de A. En el punto P de su trayectoria de altura z se cumple que

n(R+z)sinθ=n0R

A la derecha, en la figura

cosθ= dr ds dz dx

dr=dz, y dx≈ds, ya que la distancia AP es aproximadamente igual a AM.

Hacemos las mismas aproximaciones que en el apartado anterior, sustituyendo la altura final h del punto B por la altura z del punto genérico P.

sinθ= n 0 R ( n 0 +Δn )( R+z ) = ( 1+ Δn n 0 ) 1 ( 1+ z R ) 1 1 Δn n 0 z R cosθ= 1 sin 2 θ 2 Δn n 0 +2 z R

La distancia al horizonte d(h) es

dz dx = 2 z R 1+ R z Δn n 0 dx= R 2z ( 1+ R z n(z) n 0 n 0 ) 1/2 dz d(h)= R 2 0 h { z( 1+ R z n(z) n 0 n 0 ) } 1/2 dz d(h)= R 2 0 h ( zR k 1+k ( 1exp( αz ) ) ) 1/2 dz

Esta integral se puede resolver de forma numérica, empleando la función integral de MATLAB

Otras aproximaciones

La raíz cuadrada se aproxima (véase el apartado anterior)

dx R 2z ( 1 R 2z Δn n 0 )dz d(h) R 2 0 h z 1/2 ( 1 R 2z n(z) n 0 n 0 )dz

El índice de refracción también se puede aproximar

n(z) n 0 =1+kexp( αz )(1+k)=k( 1exp( αz ) )kαz

Con esta aproximación, la integral es inmediata

d(h) R 2 ( 1+ R 2 kα ) 0 h z 1/2 dz d(h) 2Rh ( 1+ R 2 kα )= 2Rh ( 1+0.092 )=123.34 h

donde la altura h y la distancia al horizonte d (h) se miden en km.

Esta aproximación es válida cuando la altura del observador es menor que 1 km por encima del nivel del mar. Para otras alturas, el error no es muy grande.

La mayor incentidumbre proviene de los valores de los parámetros k y α que dependen de la temperatura y otros factores

Comparamos la solución más exacta y la aproximada hasta la altura de 1 km, apreciamos que son casi indistinguibles

R=6371; %radio de la Tierra, en km
alfa=0.105; %en km^-1
k=2.77e-4; %parámetro
f=@(x) 1./sqrt(x-R*k*(1-exp(-alfa*x))/(1+k));
hh=(10:10:1000)/1000;
dist=zeros(1,length(hh));
i=1;
for h=hh
    dist(i)=integral(f,0,h)*sqrt(R/2); %más exacta
    i=i+1;
end
dist_a=sqrt(2*R)*(1+R*k*alfa/2)*sqrt(hh); %aproximada
plot(hh,dist, hh,dist_a)
grid on
legend('integral','aproximada','location','best')
xlabel('h, km')
ylabel('d(h), km')
title('Distancia al horizonte')

Referencias

A. P. French. How far is the horizon?. Am. J. Phys. 50 (9) Sept. 1982, pp. 795-799