11.

Dinámica Meteorológica: Ecuaciones del Movimiento Meteorología

DINÁMICA METEOROLÓGICA: ECUACIONES DEL

MOVIMIENTO
11.1 Ecuación del Movimiento Absoluto

Para derivar los diversos términos en la ecuación de un fluido en movimiento, es necesario

distinguir entre el movimiento referido a un sistema de coordenadas fijo en el espacio (Sistema
de Referencia Inercial) y el que se refiere a unas coordenadas solidario a la Tierra y que forman
parte de su rotación (Sistema de Referencia no Inercial). Al primero se le llama absoluto y al
segundo relativo. Para identificar estos movimientos usaremos los subíndices a y r.
Los principios para establecer la ecuación que rige dichos movimientos son: la ley de
conservación de la masa y la segunda ley de Newton según la cual: "La variación por unidad
de tiempo de la cantidad de movimiento de un cuerpo es igual a la suma de fuerzas que actúan
sobre él".

La cantidad de movimiento, p , de una porción infinitesimal de fluido de masa dm=ρ·dv,
es decir el producto de la masa por la velocidad absoluta viene dado por:

  
· a = ρdvVa
p = dmV

Y su variación por unidad de tiempo:

 d p d 
∑i i dt dt
=F = ρdvV(
a )
Pero como la masa permanece invariable, quedará:


 d p  dVa 
∑i Fi = dt = ρdv  dt  (11.1)
 

Las fuerzas que actúan sobre esa porción infinitesimal de fluido son:
 
p Fuerza de la presión atmosférica: Fp = −∇p·dv
 
Φ Geopotencial (gravedad “g”) Fg = −∇Φ a ρdv
 
FR Fuerzas de fricción y viscosidad FR = f R ·dv

Analicemos cada uno de ellas con algo más de detalle.

11-1
E.P.S.I. Náutica, Máquinas y Radioelectrónica Naval Universidad de La Laguna

11.1.1 La presión

  
Viene dada por: Fp = −∇p·dv , donde ∇p es el gradiente de presión por unidad de volumen,

y el conjunto la fuerza total. Para conocer de dónde procede dicha ecuación consideremos la
presión sobre un paralelepípedo elemental de fluido de volumen dx·dy·dz.
Sea en la figura el elemento de fluido ABCD, rectangular y de área dy dz.

z
C H

B
N
p dz

dy D I
A dx M
y x

A uno y otro lado actúan las presiones:

p (su sentido coincide con el eje x),
–p (su sentido es opuesto con el eje x).
∂p
Lo que aumenta la presión por unidad de distancia a lo largo del eje x es: ,
∂x
∂p
y para una distancia dx, será: dx
∂x
Por lo tanto, sobre el elemento de fluido MNHI las presiones son:

∂p ∂p
p+ dx , a la izquierda y −  p + dx  a la derecha.
∂x  ∂x 

Sobre el paralelepípedo elemental de fluido de volumen dv=dx·dy·dz las fuerzas que

actuarán sobre las caras elementales serán (Fuerza=Presión × Área):

 ∂p 
FABCD = pdydz y FMNHI =
− p + dx  dydz
 ∂x 

La fuerza neta sobre la parcela en la dirección positiva del eje x será la suma de ambas:

 ∂p   ∂p 
Fp ( x ) =
−   dxdydz =
−   dv
 ∂x   ∂x 

11-2
11. Dinámica Meteorológica: Ecuaciones del Movimiento Meteorología

Obteniendo resultados análogos sobre los ejes oy, oz.

En definitiva, la fuerza total será:

  ∂p  ∂p  ∂p   
Fp = −  i + j + k  dv = −∇p·dv
 ∂x ∂y ∂z 

Y por unidad de volumen (dividiendo por el volumen)


∑F p

= −∇p
dv

Que es el gradiente tridimensional de la presión.

 
Quedando para la unidad de masa:
∑F = −
∇p
p

dm ρ

1
Y haciendo en la práctica meteorológica = α , volumen específico o el ocupado por la
ρ
unidad de masa, la fuerza de la presión por unidad de masa:


∑ p
F 
= −α∇p
dm

EJERCICIO 11.1
Calcular la fuerza de la presión por unidad de volumen y unidad de masa en el punto M de
la figura, en donde se muestran los siguientes datos:

1004 mb

120 NM
1000 mb
● M

1002 mb

Presión: 1002 mb; diferencia de presión entre isobaras (dp): 4 mb

Distancia entre isobaras (dh): 120 millas. Sea la temperatura: 10 ℃.

11-3
E.P.S.I. Náutica, Máquinas y Radioelectrónica Naval Universidad de La Laguna


Hemos escogido para calcular el gradiente máximo el eje x en la dirección del −∇p y el
eje “y” estará lógicamente dirigido en la dirección de la isobara (coordenadas naturales),
∂p
por lo que = 0.
∂y
∂p
De este modo, todo el gradiente estará dirigido según (escogidos los ejes de esta
∂x
∂p
manera, solemos escribir ∇p = ).
∂h
La fuerza de la presión por unidad de volumen será, en módulo:

∑F p ∂p
=∇p = = 
4  mb 
=
400  N / m 2 
=1.80·10−3 [N/m3 ]
dv  
∂h 120  NM  120 × 1852  m  

p 100200
Y por unidad de masa siendo la densidad
= ρ =
*
= 1.23 [kg/m3 ]
R ⋅ T 287 ⋅ 283

Nota: Recordemos que R* es la constante de los gases ideales expresada en [Pa·m3/(kg·K)]

o [N·m/(kg·K)], por lo que, para el aire seco es aprox. 287. En [J/(mol·K)] o [N·m/(mol·K)]
es igual a 8.314472 dado que la masa molar del aire seco es 28.96·10-3 kg/mol.

N 
∑ F= p ∇p 1.80·10−3  m3 
= = 1.4633·10−3 [N/kg]
dm ρ 1.23  kg 
 m3 

11.1.2 El geopotencial (potencial de la atracción de la gravedad)

  
Fg = −∇Φ a ρdv , en donde el vector ∇Φ a es el gradiente del geopotencial y el conjunto

nos indica la fuerza de la gravedad que actúa sobre la porción infinitesimal de fluido.
Con el fin de asegurarnos de que dos parcelas de aire situadas sobre el nivel del mar a la
misma altura geométrica, corresponden a la misma superficie equipotencial, y por tanto tienen
la misma energía potencial, se utiliza, para especificar las alturas a la que está una superficie
isobárica, en vez del metro ordinario o geométrico el metro geopotencial. Por tanto, debemos
tener siempre presente que las alturas geopotenciales en un mapa de espesores o contornos,
significan realmente “energía potencial gravitatoria”.

11-4
11. Dinámica Meteorológica: Ecuaciones del Movimiento Meteorología

El geopotencial es el trabajo ejecutado para desplazar a la unidad de masa desde el nivel

del mar hasta la altura z , o la diferencia de potencial gravitatorio entre un punto cualquiera de
la atmósfera y el nivel del mar:
dΦ=gdz

Las unidades serán [m2/seg2] o, [J/kg], que es la unidad de trabajo para la unidad de masa.
Por tanto:

dΦ
g= = ∇Φ
dz

Y vectorialmente:

 
g = −∇Φ

Si integramos la expresión anterior del dΦ suponiendo g constante y dado que el

geopotencial Φ es nulo a nivel del mar:

Φ z z
∫=
0
dΦ ∫=
0
gdz g ∫ dz
0

0 g ( z − 0)
Φ −=

Haciendo Φ =1 y g=9.80665 (gravedad estándar para ϕ=45°), resulta el espesor de un

estrato desde el nivel del mar de potencial igual a la unidad:

1 = 9.80665·z [m 2 /s 2 ]

que nos da el valor del denominado decímetro dinámico: z = 0.10197 m.

Sin embargo, g no es constante por lo que, realmente, sólo coincide con el metro ordinario
a la latitud ϕ=45° y al nivel del mar.
Por acuerdo internacional se considera la siguiente expresión para el cálculo del mgp:

Φ z g
Z [mgp]= =∫ dz
9.80665 0 9.80665

Nótese que g cambia en función de la latitud ϕ y de la altura z según la siguiente expresión

(http://glossary.ametsoc.org/wiki/Acceleration_of_gravity):

11-5
E.P.S.I. Náutica, Máquinas y Radioelectrónica Naval Universidad de La Laguna

g ϕ − (3.085462·10−6 + 2.27·10−9 cos(2ϕ))·z

g= †
2
g ϕ = 9.80616·(1 − 0.0026373cos (2ϕ) + 0.000 0059 cos (2ϕ))

Si bien se puede encontrar una expresión aún más precisa para gϕ según el Ministerio de
Defensa de EE.UU. (http://earth-info.nga.mil/GandG/publications/tr8350.2/wgs84fin.pdf):

1 + 0.00193185265241sin 2 (ϕ)
g ϕ = 9.7803253359
1 − 0.00669437999013sin 2 (ϕ)

Dando lugar al siguiente comportamiento:

Figura 11.1. Variación de g con la latitud y la altura.

EJERCICIO 11.2
Calcular la altura geopotencial de un punto de la topografía de los 500 mb, en la latitud
ϕ= 30°, siendo la altura geométrica 5800 m.

†
gϕ da cuenta de la constante gravitacional corregida por la fuerza centrífuga de la Tierra.

11-6
11. Dinámica Meteorológica: Ecuaciones del Movimiento Meteorología

Aplicando la ecuación:
−6
Φ 5800 g 5800 g ϕ − (3.085462·10 + 2.27·10−9 cos(2ϕ))·z
[mgp ] = ∫
Z= = dz ∫ dz
9.80665 0 9.80665 0 9.80665
que: g ϕ 9.80616·(1 − 0.0026373 ⋅ 0.5 + 0.000 0059 ⋅ 0.25 ) =
Para ϕ=30° se tiene= 9.79324

9.79324 − (3.085462·10−6 + 2.27·10−9 ·0.5)·z

5800
Z [mgp ] ∫= dz
0 9.80665
 9.79324 5800 3.087 ⋅10−6 z 2 5800  58002
=  ⋅z 0 − ⋅ = 0.9986 ⋅ 5800 − 3.148 ⋅10−7 = 5786.78 [mgp ]
 9.80665 9.80665 2 0  2
 

Para el mismo caso, en la latitud ϕ = 45°, g ϕ = 9.80616 por lo que Z = 5794.42 [mgp].

Un mapa de 500 mb (=500 hPa) es una representación a presión constante y nos muestra
hasta que altura en mgp debemos ascender para que la presión disminuya hasta ese valor
(aproximadamente unos 5500 m.s.n.m). Por tanto, los cambios en los niveles de mgp que se
observan se deben a la presencia de regiones anticiclónicas y ciclónicas tal que, las zonas que
aparezcan con mayor altura geopotencial corresponden a collados y anticiclones, mientras que
las de menos a vaguadas y borrascas.

Figura 11.2. Mapa de superficie de 500 hPa y temperaturas a 850 mb generado por el
European Centre for Medium-Range Weather Forecasts (ECMWF). Las unidades para el
geopotencial están en decímetros dinámicos (dam) y la temperatura en grados centígrados.

11-7
E.P.S.I. Náutica, Máquinas y Radioelectrónica Naval Universidad de La Laguna

11.1.3 La fricción


Si el vector f R es la fuerza de la fricción por unidad de volumen, la fuerza total de la
 
fricción será: FR = f R ·dv

11.1.4 Ecuación del movimiento absoluto

La suma de las tres fuerzas anteriores da:

   
∑ F = − ∇p·dv − ∇Φ ρdv + f
i
i a R ·dv

Y como dicho sumatorio es igual a la variación de la cantidad de movimiento:


 dVa    
ρdv  = ( −∇p − ∇Φ ρ + f )·dv
a R
 dt 
ó bien:

 dVa    
ρ = ( −∇p − ∇Φ ρ + f )a R
 dt 

Que expresa la ecuación del movimiento absoluto, referidas las fuerzas por unidad de
volumen.
Si dividimos por la densidad ρ y haciendo 1/ρ = α (volumen específico) o el ocupado por
la unidad de masa, tendremos la ecuación del movimiento absoluto referidas las fuerzas por
unidad de masa:

 dVa    
  = −∇Φ a − α∇p + αf R (11.2)
 dt 

11.2 Ecuación del Movimiento Relativo ‡

Para obtener la relación entre el movimiento absoluto y el relativo consideremos la

siguiente figura. El punto M es un punto arbitrario que mantiene su posición relativa a la Tierra
(gira de manera solidaria con ella).

‡
Para una explicación más detallada ver el Apéndice 11-A

11-8
11. Dinámica Meteorológica: Ecuaciones del Movimiento Meteorología

y
δ

O ϕ

Figura 11.3. Velocidad de una parcela de aire que se desplaza a una velocidad Vr respecto del
punto M.

Ahora, si una porción de aire se desplaza desde M (Figura 11.3), con una velocidad relativa

Vr a la Tierra en rotación, su velocidad absoluta será:

     
Va= Vr + VM= Vr + Ω × R (11.3)

Puesto que R es el vector de posición de la porción, tenemos, derivando según el
movimiento absoluto y relativo:

 
  dR    dR 
Va =   , y, Vr =  
 dt a  dt  r

La ecuación (11.3) puede por tanto escribirse:

 
 dR   dR    §
=    +Ω×R
 dt  a  dt  r


§
Se demuestra que, en general, para una situación como la indicada, para cualquier vector b se tiene que (Teorema
 
 db   db   
de Poisson en sólidos):  =    + Ω×b
 dt a  dt r

11-9
E.P.S.I. Náutica, Máquinas y Radioelectrónica Naval Universidad de La Laguna

Luego tenemos el siguiente operador para relacionar un vector de nuestro sistema relativo
d  d  
al absoluto:  =    + Ω × .Si se lo aplicamos ahora a la velocidad absoluta tendremos:
 dt  a  dt  r

 
 dVa   dVa   
=    + Ω × Va
 dt  a  dt  r

Y sustituyendo el valor de la Va de la ecuación (11.3) en la anterior, tenemos:

   
 dVa   d (Vr + Ω × R)     
=
    + Ω × (Vr + Ω × R)
 dt a  dt r
de donde
   
 dVa   dVr   d (Ω × R )      
=    +  + Ω × Vr + Ω × (Ω × R)
 dt a  dt r  dt r

   
 d (Ω × R )   d Ω     dR   
Teniendo en cuenta que:   =   × R + Ω×  = Ω × Vr
 dt  r  dt  r  dt  r

 
 dVa   dVr      
=    + 2Ω × Vr + Ω × (Ω × R )
 dt a  dt  r

Y si efectuamos el triple producto vectorial:

   
Ω × (Ω × R ) = −Ω 2·r (ver final de este apartado)

En donde r es el vector de posición (radio vector) perpendicular al eje de la Tierra a la
latitud ϕ, y despejando la aceleración relativa:

 
 dVr   dVa    2
=    − 2Ω × Vr + Ω r
 dt  r  dt a

De donde si sustituimos la aceleración absoluta según (11.2):


 dVr       2
  = −∇Φ a − α∇p + αf R − 2Ω × Vr + Ω r (11.4)
 dt r

Puede verse que, debido a la rotación de la tierra aparecen dos fuerzas “ficticias”:

  
−2Ω × Vr y + Ω 2 r

11-10
11. Dinámica Meteorológica: Ecuaciones del Movimiento Meteorología

La primera es la aceleración de Coriolis y de la que se hablará en el siguiente apartado,

mientras que la segunda corresponde a la fuerza centrífuga debida a la rotación de la Tierra,
perpendicular al eje y dirigida hacia fuera de ella.
Al igual que la fuerza real de la atracción de la gravedad (absoluta), la fuerza centrífuga es
constante en cada lugar e independiente del movimiento atmosférico.
Ambas fuerzas podemos combinarlas en una sola, obteniéndose la gravedad relativa (ó
aparente), que es la que se observa en la superficie de la Tierra:

  
−∇Φ r = −∇Φ a + Ω 2 r

Ω2 r
ga
gr

Introduciendo la gravedad relativa en (11.4), obtenemos la ecuación del movimiento

relativo:


 dVr      
  = −α∇p − ∇Φ r + αf R − 2Ω × Vr (11.5)
 dt r
  
El término Ω × (Ω × R ) es la aceleración centrípeta del punto M (Figura 11.3) debida
 
solamente a la rotación terrestre. El vector (Ω × R ) es un vector cuyo módulo es Ω·R·sin(δ),

siendo δ el ángulo formado entre ellos y cuya dirección es perpendicular a ambos vectores. Es
fácil ver que dicho módulo se puede escribir también como Ωr, siendo r el módulo del radio

vector r a la latitud ϕ.
 
Si aplicamos nuevamente el producto escalar entre Ω y el vector anterior Ωr el resultado
 
será un vector que se dirige hacia el eje de la Tierra de módulo Ω2r. Nótese que Ω y Ωr son
perpendiculares por lo que el seno del ángulo formado entre ambos es igual a la unidad.

  
Ω × (Ωr ) = −Ω 2 r

11-11
E.P.S.I. Náutica, Máquinas y Radioelectrónica Naval Universidad de La Laguna

11.3 La Aceleración de Coriolis

 
Observando la expresión para la aceleración de Coriolis: −2Ω × Vr vemos que, por

definición de producto vectorial, la fuerza de Coriolis es perpendicular al vector Vr , actuando

hacia la derecha en el hemisferio norte y hacia la izquierda en el sur.

Supongamos que estamos fuera de la Tierra, en reposo absoluto, entonces la veríamos girar
de W-E alrededor de su eje. Si observamos una masa de aire que se traslada horizontalmente y
en línea recta desde el punto P, en cualquier latitud en el hemisferio norte y fuera del ecuador;
la masa de aire avanzará hacia el punto A.

A*

P
A
Plano
del
horizonte

Pero como el plano del horizonte gira alrededor de su vertical de W a E, mientras el aire
avanza hacia el punto que inicialmente está en su trayectoria, como el punto A va girando de
W a E, cuando el aire llegue a él ya no está allí, sino a la izquierda de la masa de aire, en A*.
Desde el espacio exterior veríamos eso, pero para quien está en la tierra que supone está
quieta, cree ver que es el aire el que se ha desviado hacia la derecha de la trayectoria.
Esta desviación aparente, debida a la rotación del plano del horizonte, es la desviación de
Coriolis. En el Hemisferio Sur ocurriría igual, pero la desviación, por razones obvias sería hacia
la izquierda.
Observemos ahora en la siguiente figura el plano del horizonte de un polo: lo veríamos
girar alrededor de su vertical que precisamente en ese lugar coincide con el eje de la Tierra, de
manera que en ese plano y sin que su vertical sufra ningún traslado, gira una vuelta completa
en veinticuatro horas.
Observemos así mismo el plano del horizonte de un punto situado en el ecuador, allí su
vertical es obviamente perpendicular al eje de la tierra, y el plano del ecuador y su vertical dan
una vuelta completa alrededor del eje de la tierra en veinticuatro horas, pero sin que el plano

11-12
11. Dinámica Meteorológica: Ecuaciones del Movimiento Meteorología

gire lo más mínimo alrededor de su vertical. Es decir que en el polo el horizonte gira alrededor
de su vertical sin traslación de ésta.

Ω
Ω

Ωz
Ωy

En el ecuador solo hay traslado de la vertical y el horizonte no gira en absoluto.

Imaginemos, por ejemplo, que estamos viendo el sol o las estrellas. En el Polo Norte giran
alrededor nuestro, por lo que para seguirlas sólo debemos girarnos, pero sin subir ni bajar la
mirada. En el ecuador sólo se trasladan por lo que no debemos girar la cabeza en horizontal,
sólo en vertical. En cualquier punto entre 0 y 90 grados de latitud el plano del horizonte y su
vertical dan vueltas alrededor del eje de la Tierra, pero simultáneamente el plano del horizonte
gira alrededor de su vertical en el mismo sentido que la Tierra, pero más lentamente. En este
caso, para seguir al sol o a las estrellas deberemos tanto girarnos como cambiar el ángulo
vertical. Como vemos a continuación para ϕ = 90°, y en [grados/hora] seria 15°·sin(ϕ), o lo que
es lo mismo, gira con la misma velocidad que la Tierra.
No gira nada para ϕ = 0°.
Y gira a la mitad de la velocidad de la Tierra para ϕ = 30°
Como podremos comprobar en la figura la componente perpendicular al plano del

horizonte de un observador, Ω z , vale:
 
Ω z =Ω sin(ϕ)k

Que es realmente la componente que interviene en el movimiento horizontal, por lo que la

aceleración de Coriolis (o la fuerza de Coriolis por unidad de masa) valdrá para el movimiento
horizontal:
      
FCo (horizontal) =−2Ω z × V =−2Ω sin(ϕ)k × V =− fk × V (11.6)

A la parte 2Ωsin(ϕ), se le llama parámetro de Coriolis y como ya hemos visto se le asigna

la letra “ f ”, por lo que al conjunto se le designa como “fV”. La acción de Coriolis es

11-13
E.P.S.I. Náutica, Máquinas y Radioelectrónica Naval Universidad de La Laguna

despreciable por tanto siempre en las proximidades del ecuador, así como en cualquier lugar si
el recorrido del viento es corto.

Figura 11.4. Componente de Coriolis para un viento Vr que se desplace hacia el norte según
un sistema de referencia situado en un punto de la superficie del planeta en el HN y en el HS.

11.3.1 Las componentes de las ecuaciones

Escogemos el eje de las X hacia el Este, el de las Y hacia el Norte y el de las Z hacia el
cénit. La velocidad angular de la rotación de la Tierra tiene, como hemos visto en el apartado
anterior las componentes:

Ωx=0 Ωy=Ωcos(ϕ) Ωz=Ωsin(ϕ)

Ahora, si u , v y w son las componentes de la velocidad del viento según los ejes X, Y, Z
respectivamente, las componentes de la fuerza de Coriolis son:

  
i j k
  
−2Ω × V =−2 0 Ω cos(ϕ) Ω sin(ϕ) =−2Ω[ w cos(ϕ) − v sin(ϕ)]i +
u v w
 
+2Ω[−u sin(ϕ)] j − 2Ω[−u cos(ϕ)]k

Con lo que la ecuación vectorial (11.5) puede ahora sustituirse por sus respectivas
ecuaciones, componentes de la aceleración:

11-14
11. Dinámica Meteorológica: Ecuaciones del Movimiento Meteorología

 du ∂p
 dt = −α ∂x − 2Ω[ w cos(ϕ) − v sin(ϕ)] + αf x

 dv ∂p
 = −α − 2Ω[u sin(ϕ)] + αf y (11.7)
 dt ∂y
 dw ∂p ∂Φ
 = −α − 2Ω[−u cos(ϕ)] − + αf z
 dt ∂z ∂z

Y si consideramos las velocidades verticales despreciables y no incluimos el rozamiento:

 du ∂p
 dt = −α ∂x − 2Ω sin(ϕ)( −v )
 dv (11.8)
 = −α ∂p − 2Ω sin(ϕ)( +u )
 dt ∂y

Esta ecuación puede ser escrita en forma vectorial como:


dVH   
= −α∇p − fk × VH (11.9)
dt

11-15
E.P.S.I. Náutica, Máquinas y Radioelectrónica Naval Universidad de La Laguna

Apéndice A

ω
z z′ P

 
r r′ y′

 O′
O R
y

x′
x

Figura A.1. Sistema de referencia inercial (S.R.I. - Oxyz ) y no inercial (S.R.N.I - O′x′y′z ′ )

Si suponemos un sistema de referencia inercial Oxyz y otro no inercial O′x′y′z ′ , tendremos

que el vector posición de un punto situado en P puede ser calculado como:

   
r = r ′ + R = x′iˆ′ + y′ˆj ′ + z ′kˆ′ + R

Si suponemos ahora que el punto P se mueve a cierta velocidad v′ con respecto al S.R.N.I.,

  dR
el cual a su vez rota con una velocidad angular ω y se traslada con una cierta velocidad V =
dt
podemos calcular la velocidad de P respecto del S.R.I. de la siguiente forma:

   
 dr dr ′ dR dx′ ˆ dy′ ˆ dz ′ ˆ diˆ′ djˆ′ dkˆ′ dR
v= = + = i′ + j′ + k ′ + x′ + y ′ + z ′ +
dt dt dt dt dt dt dt dt dt dt

dado que para el observador situado en el sistema inercial los vectores unitarios de O′x′y′z ′ , es

decir iˆ′, ˆj ′, kˆ′ , rotan y por tanto cambian con el tiempo, aplicando el teorema de Poisson
tendremos que:

diˆ′  ˆ djˆ′  ˆ dkˆ′  ˆ

= ω× i′; = ω× j ′; = ω× k ′
dt dt dt 

 dr   ˆ ˆ ˆ dR    dR
v= = v′ + ω× ( x′i′ + y′j ′ + z ′k ′) + = v′ + ω× r ′ +
dt dt dt

Mientras que para la aceleración tendremos que:

11-16
11. Dinámica Meteorológica: Ecuaciones del Movimiento Meteorología

  
 d 2r d 2r′ d 2 R d 2 x′ ˆ d 2 y′ ˆ d 2 z ′ ˆ dx′ diˆ′ dy′ djˆ′ dz ′ dkˆ′
a= 2 = + = i′ + 2 j ′ + 2 k ′ + + + +
dt dt 2 dt 2 dt 2 dt dt dt dt dt dt dt dt

dx′ diˆ′ dy′ djˆ′ dz ′ dkˆ′ d 2iˆ′ d 2 ˆj ′ d 2 kˆ′ d 2 R
+ + + + x′ 2 + y ′ 2 + z ′ 2 + 2
dt dt dt dt dt dt dt dt dt dt

Como:

 
d 2iˆ′ d ω ˆ  diˆ′ d ω ˆ   ˆ
=
dt dt
× i′ + ω× =
dt dt
(
× i′ + ω× ω× i′ )
 
d 2 ˆj ′ d ω ˆ  djˆ′ d ω ˆ   ˆ
dt
=
dt
× j ′ + ω× =
dt dt
(
× j ′ + ω× ω× j ′ )
 
d 2 kˆ′ d ω ˆ  dkˆ′ d ω ˆ   ˆ
dt
=
dt
′
× k + ω× =
dt dt
(
× k ′ + ω× ω× k ′ )

    dx′ ˆ dy′ ˆ dz ′ ˆ  d ω
a =a′ + 2ω×  i′ + j′ + k′ + × ( x′iˆ′ + y′ˆj ′ + z ′kˆ′) +
 dt dt dt  dt

     d 2R
+ ω× (ω× ( x′i′ + y′j ′ + z ′k ′)) + 2
 dt
 2
       dω  d R
a= a′ + 2ω× v′ + ω× (ω× r ′) + × r′ + 2
dt dt

Para el caso particular de un sistema situado en la superficie de la Tierra y que rote con

ella tenemos que, de una parte, la velocidad angular de la Tierra, ω , no cambia, y por tanto

dω
= 0 . Además, para el movimiento de traslación del sistema O′ situado sobre la superficie
dt
de la Tierra:

   
dR  d 2R d   d ω   dR   
dt
= ω× R y por tanto
dt 2
=
dt
ω× R =
dt
( dt
)
× R + ω× = 0 + ω× (ω× R)

Con esto, substituyendo obtenemos para la velocidad y aceleración relativas:

     
v′= v − ω× r ′ − ω× R
         
a′= a − 2ω× v′ − ω× (ω× r ′) − ω× (ω× R)

    
donde −2ω× v′ es la aceleración de Coriolis y −ω× (ω× R) es la centrífuga.
La Figura A.2 muestra las componentes del viento en un ejemplo en el que nuestro
elemento de atmósfera parte de cada uno de los puntos cardinales del S.R.N.I., a una distancia
 
r ′ del centro O′ , con una velocidad relativa v′ respecto de dicho centro, y según una dirección
hacia el Este.

11-17
E.P.S.I. Náutica, Máquinas y Radioelectrónica Naval Universidad de La Laguna

Figura A.2. Velocidades instantáneas para un elemento de atmósfera que se desplaza a una

velocidad relativa v′ hacia el E partiendo desde cada punto cardinal del S.R.N.I. En verde las
    
velocidades relativas v′ , en rojo ω× R , en naranja las ω× r ′ y en negro el resultado final.

En el caso más común en el que el movimiento del elemento de atmósfera parta desde O’,

es decir, colocamos a nuestro S.R.N.I. en el punto P, no existe, lógicamente, la distancia r ′
hasta el centro, por lo que tendremos para la velocidad y aceleración relativas:

   
v′= v − ω× R
      
a′= a − 2ω× v′ − ω× (ω× R)

11-18