Introduccin a la Dinmica de la Atmsfera 1

5. Circulacin, vorticidad y divergencia

Como hemos visto la atmosfera est caracterizada por la presencia de torbellinos de todas
las escalas espaciales y temporales. En latitudes medias el torbellino mas importante es el
ciclon de escalas sinpticas caracterizado por centros de presin con escalas horizontales de
500-1000 km y vientos que giran alrededor. En este capitulo estudiaremos dos cantidades
fisicas que ayudan a cuantificar la rotacion: circulacion y vorticidad. Derivaremos una
ecuacion para la vorticidad de un fluido y veremos que cambios en la vorticidad estn
directamente asociados a cambios en la divergencia que, por continuidad, implican
movimientos verticales. Por ltimo, dado que la atmsfera en latitude medias est en
balance de viento trmico aproximado, consideraremos aproximaciones de las ecuaciones
de vorticidad y termodinmica apropiadas para flujos con Ro<<1, lo cual dar lugar al
sistema cuasi-geostrfico de ecuaciones que son fundamentales para la comprensin de la
dinmica del cicln extratropical.
5.1 Circulacin
Se define la circulacin C en una region simplemente conexa como
donde dl representa el vector desplazamiento a lo largo del borde de un elemento de fluido
como se ilustra en la figura 5.1.
Figura 5.1 Calculo de la circulacin usando un circuito de fluido.
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 2
Por convencin la integral anterior se realiza en sentido antihorario, por lo que C>0 (C<0)
corresponde a una rotacin ciclnica (anticiclnica) en el hemisferio norte. Reeescribiendo
la ecuacin en coordenadas cartesianas en la direccin horizontal
donde U y V representan las velocidades tangenciales alrededor de un elemento de fludo
en las direcciones x e y, respectivamente.
La circulacin es una medida de la rotacion. Para demostrarlo podemos considerar un anillo
de fluido de radio R con rotacin solida a una velocidad angular . En este caso U=xR
donde R es la distancia del eje de rotacion al anillo. Entonces, la circulacion alrededor del
anillo esta dada por
C=

U.dl =

0
2n
DR
2
d \=2DnR
2
En este caso la circulacion es 2 veces el momento angular del anillo de fluido alrededor de
su eje de rotacion. A su vez, C/R
2
=2 de tal forma que la circulacion dividida por el rea
encerrada por el anillo es dos veces la velocidad angular de rotacion del anillo. A diferencia
del momento angular, la circulacion puede ser calculada sin referencia de un eje de
rotacion.
En general estamos interesados en estudiar cambios en la circulacion pues estos se
manifiestan como intensificaciones de las altas y bajas. Por lo tanto, estudiaremos la
derivada temporal de C a lo largo del movimiento del fludo. Considerando el caso
tridimensional dl=(dx,dy,dz)
Por la regla de la cadena
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 3
Sustituyendo
El segundo trmino de la derecha puede escribirse como
pues es la integral cerrada alrededor de un elemento de fluido.
Ahora, si usamos la aceleracion absoluta podemos escribir
y como solamente la fuerza gradiente de presin y la gravitacional (despreciando friccion)
influencian la aceleracion absoluta
donde
en una superficie de altura constante. Puesto que la componente vertical del vector
desplazamiento, dl, es simplemente dl.k=dz, entonces
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 4
y
pues d es un diferencial perfecto y por lo tanto no se realiza trabajo contra la gravedad si
la parcela de fluido termina su recorrido en el mismo lugar que empezo,
independientemente de su recorrido. Por lo tanto
El trmino de la izquierda describe la razon de cambio lagrangiano de la rotacin del fluido,
por lo que la ecuacion anterior representa el analogo de la aceleracion angular en cuerpos
solidos pero para un fluido. El trmino de la derecha es el equivalente del torque para un
fluido y se conoce como trmino solenoidal.
Para un fluido barotropico, la ley de gas ideal implica que
por lo que
el trmino solenoidal es nulo y la circulacin absoluta se conserva siguiendo el movimiento
de las parcelas, lo cual se conoce como el teorema de circulacion de Kelvin.

La atmsfera real no es barotrpica, o sea que las superficies de presin y densidad
constante se intersectan. Consideremos el caso de la brisa marina, donde la columna de aire
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 5
es mas baja sobre el oceano por ser mas fria que sobre el continente. En este caso, aun si la
presion es muy parecida en superficie, las isobaras se inclinan hacia el mar mientras que las
isopicnals se inclinan hacia la tierra mas clida. La interseccion entre isobaras e isopicnals
forma una serie de paralelogramos llamados solenoides, de donde proviene el nombre del
trmino (Figura 5.2)
Figura 5.2 Clculo de la circulacin en la frontera mar-tierra.
Para el caso de la brisa marina mostrado es podible calcular el cambio lagrangiano de la
circulacion absoluta evaluando -dp/ alrededor del camino indicado. El movimiento en AB
es a lo largo de una isbara y por lo tanto no contribuye a la integral. En BC, dp<0 por lo
que el trmino es positivo. En CD el movimiento es nuevamente isobrico y no hay
contribucion a la integral. E DA el trmino es negativo. Como la densidad promedio en BC
es menor que en DA la contribucion positiva de BC es mayor que la contribucion negativa
de DA resultando en un aumento de la circulacion absoluta alrededor del camino. Por lo
tanto en este ejemplo aparecera una circulacion en la cual el fluido menos denso sube y el
mas denso baja. El resultado del movimiento ser inclinar las isopicnals de tal forma que se
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 6
vuelvan mas paralelas a las isobaras, o sea tiende a una situacion barotrpica en la cual el
cambio en la circulacion sera nula. Tal circulacin tambien baja el centro de masa del
fluido y reduce la energia potencial del sistema. Esta reduccion de la energia potencial es en
realidad una conversion a energa cintica del movimiento, conservando la energia total.
Hasta ahora consideramos la circulacin absoluta. Para calcular la circulacin relativa, las
relevante para nuestro estudio debemos calcular la circulacion que resulta de la rotacion de
la Tierra y luego restarla. Para ello consideremos la velocidad alrededor de un circulo de
latitud
donde R=acos . Por lo tanto el movimiento zonal hacia el oeste resultante por el
movimiento de la tierra es U=Dacos . Ahora calculamos la circulacion alrededor del
circuito latitud-longitud mostrado en la figura 5.3. En coordenadas esfricas la longitud de
un elemento en la direccion zonal es dx=acos6\ . Como la rotacion de la Tierra no
contribuye a movimientos en la direccion meridional la circulacion resultante Ce es
Usando las identidades trigonomtricas
podemos escribir
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 7
El rea de la caja es
y considerando los lmites
podemos expresar el rea como
Combinando las ecuaciones se ve que la circulacion que resulta de la rotacion terrestre es
o sea que
Combinando las ecuaciones llegamos a una expresion para la razn de cambio lagrangiana
de la circulacion relativa
la cual se conoce como el teorema de circulacion de Bjerkness.
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 8
Figura 5.3 Caja latitud-longitud para calcular la circulacin asociada a la rotacin
terrestre.
Para un fluido barotropico es posible integrar el teorema de Bjerkness siguiendo el
movimiento desde un estado inicial (sealado por 1) hasta un estado final (2)
lo cual indica que para un fluido barotropico la circulacion relativa para un circuito cerrado
cambiar si cambia el rea horizontal encerrada o cambia la latitud. Adems, se desarrollar
una circulacion absoluta negativa en el H.N slo si se advecta desde el H.S. un circuito
cerrado de fluido. Las bajas y altas anmalas del viento gradiente son ejemplos de
circulacion absoluta negativa.
5.2 Vorticidad
En el primer captulo consideramos la vorticidad como una propiedad cinemtica del
fludo. Adems, la vorticidad es una medida microscpica de la rotacin de un fludo y es
ms fcil de tratar que la circulacin. La vorticidad es una cantidad vectorial que se define
como el rotor del campo de velocidades. Como para caracterizar la atmsfera trabajamos en
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 9
un sistema de coordenadas que gira es necesario definir la vorticidad relativa y la vorticidad
planetaria, cuya suma es la vorticidad absoluta.
La vorticidad relativa del fluido es w= x u , mientras que la vorticidad planetaria se
define como 2, y la vorticidad absoluta es
w
a
= x u+2D
.
En general para estudiar los movimientos atmosfricos consideraremos nicamente la
componente vertical (local) de la vorticidad, la cual se expresa como
=k. w=k.

V=
v
x

u
y
Regiones con vorticidad relativa positiva en el H.N. estn asociadas a ciclones; ciclones en
el H.S. estn asociado a regiones con vorticidad relativa negativa. Por lo tanto la
distribucin de vorticidad relativa es muy til para el diagnstico de sistemas de latitudes
medias.
Para determinar la relacin fsica entre vorticidad y circulacin consideremos el elemento
de fludo mostrado en la figura 5.4. La velocidad en el lado A est dada por u mientras que
en el lado D est dada por v. Expandiendo u y v segun Taylor es posible obtener
expresiones para las velocidades en los lados C y B. Entonces, la circulacin es
Como el rea del elemento de fluido es xy encontramos que en el limite del rea
tendiendo a cero la vorticidad relativa es simplemente la circulacin relativa dividida el
rea del elemento.
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 10
Figura 5.4 - Circuito para vincular circulacion con vorticidad
Recordando el clculo de la circulacin terrestre dCe/dt=fdA/dt y la figura 5.3, la
componente vertical de la vorticidad planetaria es el parametro de Coriolis f, y es igual a
dos veces la razon de rotacion local de la Tierra (ver figura 5.5)
f =2Dsin .
Figura 5.5 Vorticidad planetaria (no absoluta!)
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 11
Por lo tanto la componente vertical de la vorticidad absoluta es
En forma mas general la relacin entre la vorticidad y la circulacin est dada por el
teorema de Stokes aplicado al vector velocidad
donde A es el rea encerrada por el contorno y n es un versor normal al elemento de rea
dA (positivo de acuerdo a la regla de la mano derecha). Por lo tanto, el teorema de Stokes
establece que la circulacin alrededor de cualquier contorno cerrado es igua a la integral de
la componente normal de la vorticidad sobre el rea encerrado por el contorno. As, para un
rea finita la circulacin dividida entre el rea da el promedio de la componente normal de
la vorticidad en la regin considerada.
5.3 Vorticidad en coordenadas naturales
La interpretacin de la vorticidad es mas simple considerando coordenadas naturales.
Consideremos la componente vertical de la vorticidad. Si calculamos la circulacin en el
contorno infinitesimal mostrado en la figura 5.6 obtenemos
Figura 5.6 Calculo de la circulacion en un circuito infinitesimal en coordenadas
naturales.
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 12
Como
donde 6 es el cambio angular en la direccin del viento en la distancia 6 s .
Entonces
y, como la vorticidad es circulacin por unidad de rea, cuando el rea tiende a cero se
obtiene (notar que la figura 5.6 es una foto y por lo tanto muestra las lneas de corriente)
donde Rs es el radio de curvatura de las lneas de corriente. La expresin anterior sugiere
que la vorticidad vertical es la suma de dos componentes: (1) la variacin de la velocidad
del flujo en la direccin normal, llamada vorticidad de corte, y (2) variacin de la
direccin del flujo a lo largo de una lnea de corriente, llamada vorticidad por curvatura.
Por lo tanto, an movimientos rectilneos pueden tener vorticidad si la velocidad cambia en
la direccin normal al eje del flujo. En el ejemplo mostrado en la figura 5.7 un elemento de
fluido situando al norte del maximo de velocidades tender a girar en sentido antihorario,
mientras que un elemento de fludo al sur del mximo tender a girar en sentido horario.
El caso de un flujo constante en una regin con curvatura se muestra en la figura 5.8. A
medida que las parcelas se mueven a travs de la vaguada la parte superior de una rueda
con paletas imaginaria debe recorrer una distancia mayor que la parte inferior, generando
un giro horario, consistente con una vorticidad por curvatura negativa.
Notar que el H.S. el mximo de vorticidad relativa tiende a ocurrir en una cua al norte de
la regin del mximo de vientos pues los dos trminos V/Rs y (-dV/dn) son positivos,
mientras que el mnimo de vorticidad ocurre en una vaguada al sur del mximo de vientos
ya que los trminos V/Rs y (-dV/dn) son negativos.
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 13
Figura 5.7 Ejemplo de vorticidad de corte.
Figura 5.8 Ejemplo de vorticidad por curvatura.
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 14
5.4 Vorticidad potencial
Tomando el diferencial logaritmo de la temperatura potencial y usando la 1a ley de la
termodinmica es posible mostrar que
c
p
d ln0
dt
=

Q
T
=6S
donde S es la entropa. Esta ecuacin establece que en movimientos adiabticos se conserva
la temperatura potencial. Los flujos en latitudes medias lejos de reas de precipitacin
conservan y por lo tanto es til usar superficies isentrpicas como coordenada vertical
para describir estos movimientos.
Considerando la ley de gases ideales, la definicin de la temperatura potencial puede
expresarse como una relacin entre la presin y la densidad para una superficie de
constante (o isentrpica)
o
Como cp=R-cv, vale
lo cual muestra que para un flujo en superficies isentrpicas la densidad es una funcin
solamente de la presin, o sea que el flujo se comporta como barotrpico. Por lo tanto el
trmino solenoidal es nulo y el teorema de Bjerkness puede escribirse como (C
rel
=C)
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 15
donde C es evaluado en una superficie isentrpica. Como =C/A, para una parcela
infinitesimal de aire vale
o sea que el producto
(
0
+ f ) A
es constante en un flujo adiabtico (y sin friccin) en
superficies isentrpicas.

es la componente vertical de la vorticidad relativa evaluada en

una superficie isentrpica.
Asumamos ahora que la parcela se mueve entre dos superficies de temperatura potencial
constante
0
y
0
+ que estn separadas por un -p como muestra la figura 5.9. La masa
de la parcela, M=(-p/g)A se debe conservar siguiendo el movimiento (por continuidad).
Por lo tanto
ya que M y son constantes. Sustituyendo en la ecuacion
y tomando el lmite p->0 se obtiene
La cantidad P [K/kg m
2
/s] es la forma en coordenadas isentrpicas de la vorticidad
potencial de Ertel. Se define con el signo negativo para que sea (en general) positivo en el
H.N. Por lo tanto la vorticidad potencial se conserva siguiendo el movimiento en un flujo
adiabtico y sin friccin. La vorticidad potencial se puede ver como el cociente entre la
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 16
vorticidad absoluta y la altura del vrtice. En la expresin de arriba la altura est dada por
la distancia entre dos superficies de temperatura potencial medidas en coordenadas de
presin (
0
p
) . La expresin de P indica que la vorticidad de una parcela de aire puede
cambiar solamente por variaciones en la latitud y en la estabilidad esttica vertical.
Figura 5.9 Columna cilndrica de aire movindose adiabticamente, conserva P.
En un fludo homogneo incompresible la expresin de conservacin de la vorticidad
potencial es mas simple. Como la densidad es constante el rea horizontal debe ser
inversamente proporcional a la altura h de la parcela
Sustituyendo para eliminar A obtenemos
+ f
h
=constante
donde se evala en superficies de altura constante.
5.4.1. Aplicacin: restriccin si h es constante.
Si la altura h de la columna es constante la ecuacin anterior establece que la vorticidad
absoluta se conserva siguiendo el movimiento, lo cual impone una fuerte restriccin al
flujo. Supongamos que en un cierto punto (x0,y0) del H.N. el flujo es en la direccin zonal
y la vorticidad relativa es nula, o sea que la vorticidad absoluta es f
0
. Por lo tanto, el
movimiento a lo largo de cualquier trayectoria que pase por (x0,y0) debe satisfacer +f=f
0
.
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 17
Como f aumenta hacia el norte las trayectorias que se curvan hacia el norte deben cumplir
=f
0
-f<0 y las que se curvan hacia el sur deben cumplir =f
0
-f>0. No obstante, como
muestra la figura 5.10 si el flujo es del oeste una curvatura hacia el norte aguas abajo
implica >0 y una curvatura hacia el sur aguas abajo implica <0. Por lo tanto vientos del
oeste debe permanecer puramente zonal para conservar la vorticidad absoluta. El viento del
este, al contrario, puede curvarse hacia el norte o hacia el sur y conservar vorticidad
absoluta.
Figura 5.10 Conservacin de vorticidad absoluta para trayectorias curvilneas.
5.4.2. Aplicacin: Flujo sobre cadeas montaosas
Cuando cambia la altura h de la columna se conserva la vorticidad potencial, no la
vorticidad absoluta. La conservacin de vorticidad potencial explica por qu la interaccin
del viento de los oestes con cadenas montaosas como los Andes puede inducir ondas
corriente abajo (figura 5.11). A medida que el aire fluye sobre la montaa la temperatura
potencial se conserva por lo que la superficie isentrpica de K se eleva sobre la montaa.
A mayor altura la superficie + K se eleva pero en menor medida pues la anomalia se
distribuye en una distancia mayor. Esto causa un aumento del espesor entre capas antes de
la montaa que es muy pequeo y despreciaremos, y principalmente que el espesor entre
las dos superficies isentrpicas sobre la montaa disminuya. Debido a esto ltimo, para
mantener la vorticidad potencial constante la vorticidad absoluta tambin debe disminuir.
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 18
Figura 5.11
Vemoslo con mas detalle. Si consideramos que el flujo al oeste de la montaa tiene
vorticidad relativa nula, como f<0 es aproximadamente constante, al subir la montaa el
flujo debe cumplir la siguiente relacin
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 19
f g
0
p

1
=( f +
2
) g
0
p

2
=f (

0
p

0
p

2
1)
Como la estabilidad vertical arriba de la montaa es mayor, 2>0 y el flujo debe asumir
vorticidad positiva, o sea debe girar en sentido antihorario en el H.S. Cuando la columna de
aire pas por encima de la montaa estar a latitudes menores por lo que f<0 tendr menor
magnitud. No obstante, como la magnitud de f es siempre mucho mayor que para escalas
sinopticas, se tiene que ( f +)0 . Por lo tanto para mantener la vorticidad potencial
constante el flujo debe adquirir vorticidad negativa (giro horario).
Cuando la parcela vuelva a su latitud de origen todava tendr una componente meridional
de velocidad y continuar moviendose hacia el polo adquiriendo ahora una vorticidad
positiva (giro antihorario). Cuando vuelva nuevamente a su latitud original tendr una
componente meridional que har disminuir la magnitud de f por lo que la columna deber
adquirir vorticidad negativa y as sucesivamente. La parcela seguir entonces moviendose
aguas abajo describiendo un movimiento ondulatorio, necesario para mantener la vorticidad
potencial constante. En general se observa una cua sobre la montaa y una vaguada aguas
abajo (figura 5.12).
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 20
Figura 5.12
5.5 Ecuacin de vorticidad
La ecuacion para la vorticidad absoluta puede derivarse a partir de las ecuaciones de
momento. Aqu nos restringiremos a derivar la ecuacion para la componente vertical de la
vorticidad en un fluido sin friccin
Del capitulo anterior teniamos
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 21
d u
dt
fv=
1
j
p
x
d v
dt
+ fu=
1
j
p
y
d
dt
=

t
+u

x
+v

y
+w

z
Aplicando

y
a la primer ecuacion,

x
a la segunda ecuacion y restando se eliminan
los terminos de presion y se obtiene
Como
d f
dt
=v
df
dy
la cual es la ecuacin de vorticidad en coordenadas de altura. La ecuacin establece que la
razn de cambio de la vorticidad absoluta est dada por la suma de tres trminos: (1) el
trmino de divergencia, (2) el trmino de inclinacin, (3) el trmino solenoidal.
Analicemos cada trmino por separado. Si la ecuacin de vorticidad est dominada por la
divergencia horizontal actuando en el fludo se tiene
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 22
Si se considera el caso de un fluido con vorticidad ciclnica inicial y la divergencia es
positiva, la ecuacin anterior indica que la vorticidad absoluta tender a ser cada vez mas
anticiclnica (figura 5.13). Esto es pues la circulacin debe conservarse: como la vorticidad
es circulacin por unidad de rea, si el rea encerrada por una cadena de parcelas de fludo
aumenta debido a la divergencia, entonces la vorticidad debe decrecer. Vale lo opuesto para
el caso de convergencia.
Por lo tanto, el trmino de divergencia puede pensarse como el anlogo para un fludo del
cambio en la velocidad angular que resulta de un cambio en el momento de inercia de un
rgido cuando se conserva el momento angular. Este resultado tiene una gran importancia
para los sistemas de latitudes medias: sistemas de baja presin en superficie estn
caracterizados por convergencia y por lo tanto son lugares de produccin de vorticidad
ciclnica en niveles bajos (y lo opuesto para los anticiclones).
Figura 5.13 Efecto de la divergencia en la vorticidad.
Consideremos ahora el trmino de inclinacin. En este caso se tiene
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 23
y para mayor simplificacin consideremos una situacin donde slo uno de los trminos de
la derecha es importante. Consideremos la situacin en la cual se tiene un cortante vertical
de la velocidad zonal y un gradiente meridional en la velocidad vertical positivo (figura
5.14). La ecuacin correspondiente es
d (+ f )
dt
=
w
y
u
z
Figura 5.14 Efecto del cortante vertical en la vorticidad horizontal y generacin de
vorticidad vertical.
Como ilustra la figura un cortante vertical en u induce a una paleta imaginaria en la
direccin meridional a girar en sentido antihorario, de tal forma que se puede pensar que un
tubo de aire alineado en la direccin y posee vorticidad positiva en esa direccin. Si adems
existe un gradiente de movimiento vertical en la direccin y
w
y
>0
entonces el extremo
norte del tubo ascender relativo al extremo sur. Por lo tanto aparecer una componente de
rotacin antihoraria (positiva) en la direccin z, que dar lugar a un incremento en la
vorticidad absoluta.
Por ltimo consideremos el trmino solenoidal. Se puede mostrar que este trmino es el
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 24
equivalente microscpico del trmino solenoidal en el teorema de circulacin. Para ello
aplicamos el teorema de Stokes al trmino solenoidal en la ecuacin de la circulacin
Como
cuya componente vertical es
lo cual demuestra que el trmino solenoidal en la ecuacin de vorticidad es el trmino
solenoidal en la ecuacin de circulacin dividido entre el rea del elemento de fludo.
Como vimos anteriormente el trmino solenoidal tiende a reordenar la masa del fludo de
tal forma de llegar al estado con menor energa potencial. En el proceso de reordenamiento
en general se produce vorticidad. Consideremos como ejemplo la configuracin de p y
que caracteriza adveccin de aire fro por el viento zonal geostrfico en el H.N. (figura
5.15). Si no existe divergencia ni inclinacin, la ecuacin es
Puesto que en este ejemplo
j
x
0
p
y
0

el trmino solenoidal generar vorticidad ciclnica. La vorticidad generada tender a rotar
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 25
las lneas de igual densidad hasta que estn paralelas a las isbaras en una configuracin en
la cual alta presin corresponda a alta densidad, y vice versa (o sea el estado de menor
gradiente de temperatura y menor energa potencial).
Figura 5.15 Configuracin de isbaras y lneas de igual densidad caracterizando
adveccin de aire fro por el viento geostrfico en el H.N. Las flechas indican el viento
geostrfico y el crculo la circulacin solenoidal inducida.
5.5.1. Ecuacin de vorticidad en coordenadas de presin
Para evitar la referencia a la densidad en las ecuaciones y hacer desaparecer al trmino
solenoidal (pues dp=0) derivaremos la ecuacin de vorticidad en coordenadas de presin.
Las ecuaciones de momento en coordenadas isobricas en ausencia de friccin son
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 26
De nuevo, tomando las derivadas cruzadas y haciendo cuantas se obtiene
la cual describe que cambios locales en la vorticidad absoluta se deben a (1) adveccin
horizontal, (2) adveccin vertical, (3) el trmino de divergencia, (4) el trmino de
inclinacin.
5.5.2. Anlisis de escala
A continuacin realizamos un anlisis de escala a la ecuacin de vorticidad escala para
hallar los trminos ms importantes en la descripcin de los movimientos sinpticos. Las
escalas tpicas eran
donde de nuevo hemos elegido una escala de tiempo advectiva pues los patrones de
vorticidad, como aquellos de presin, se mueven a velocidades comparables a la velocidad
horizontal del viento.
Usando estas escalas obtenemos
por lo tanto
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 27
o sea que para sistemas sinpticos en latitudes medias la vorticidad relativa es menor que la
vorticidad planetaria y por lo tanto podemos escribir el trmino de divergencia de la forma
Entonces los trminos de la ecuacin de vorticidad escalan de la forma
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 28
En los ltimos tres trminos se us una desigualdad para notar que en cada caso es posible
que los trminos individuales tiendan a cancelarse de tal forma que la magnitud del trmino
total es menor a la indicada. De hecho, esto debe ser el caso para el trmino de divergencia.
Como se ve del anlisis de escala los trminos mayores son la adveccin de vorticidad
horizontal, la tendencia local de vorticidad y la divergencia. Para que la divergencia sea
balanceada por los dems trminos debe valer que
Por lo tanto, para movimiento sinpticos en latitudes medias la ecuacin de vorticidad
puede ser aproximada por
d
h
(+ f )
dt
=f (
u
x
+
v
y
)
d
h
dt
=

t
+u

x
+v

y
la cual establece que el cambio en la vorticidad absoluta siguiendo el movimiento
horizontal en escalas sinpticas est dado aproximadamente por la concentracin o dilucin
de vorticidad planetaria causada por la convergencia o divergencia del flujo horizontal,
respectivamente.
Esta aproximacin no es vlida en el centro de ciclones muy intensos y en esos casos se
debe retener la vorticidad relativa en el trmino de divergencia:
d
h
(+ f )
dt
=(+ f )(
u
x
+
v
y
)
y es la concentracin o dilucin de vorticidad absoluta la cual da lugar a cambios siguiendo
el movimiento.
La aproximacin anterior no es vlida cerca de los frentes pues la escala horizontal de las
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 29
zonas frontales es de 100 km y la velocidad vertical es del orden de 10 cm/s. En este caso
es necesario retener los trminos de adveccin vertical, inclinacin y solenoidal ya que son
comparables con el trmino de divergencia.
5.5.3. Ecuacin de vorticidad para un fludo barotrpico
Consideremos que la atmsfera puede suponerse homognea e incompresible de altura
variable h(x,y,t)=z2-z1, donde z2 y z1 son las alturas de las fronteras superior e inferior,
respectivamente.
En este caso la ecuacin de continuidad es
y la ecuacin de vorticidad queda
d
h
(+ f )
dt
=(+ f )
w
z
Para un fludo barotrpico la velocidad geostrfica es independiente de la altura.
Asumiendo que podemos aproximar el viento y la vorticidad por sus aproximaciones
geostrficas podemos integrar entre z1 y z2 y obtenemos
h
d
h
(
g
+f )
dt
=(
g
+ f )| w( z
2
)w( z
1
)
Pero como w=dz/dt, y h=h(x,y,t)
w( z
2
)w( z
1
)=
dz
2
dt

dz
1
dt
=
d
h
h
dt
sustituyendo
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 30
1
(
g
+f )
d
h
(
g
+f )
dt
=
1
h
d
h
h
dt
d
h
ln(
g
+f )
dt
=
d
h
ln h
dt
lo cual implica que
d
h
dt
(
g
+ f )
h
=0
que es el teorema de conservacin de la vorticidad potencial (geostrfica) para un
fludo barotrpico, obtenido por primera vez por C. G. Rossby.
Si el flujo es puramente horizontal (=0), como es el caso para un flujo barotrpico de
profundidad constante el trmino de divergencia desaparece y obtenemos la ecuacin de
vorticidad barotrpica
d
h
dt
(
g
+ f )=0
la cual establece que la vorticidad absoluta se conserva siguiendo el movimiento horizontal.
En forma mas general, la vorticidad absoluta se conservar en cualquier capa del fludo en
la cual la divergencia horizontal sea nula. Para movimientos horizontales no divergentes el
flujo puede ser presentado por una funcin corriente (x,y) tal que
V

=(

y
,

x
)
y la vorticidad se escribe como
=
2

Por lo tanto la ecuacin de conservacin de la vorticidad barotrpica puede ser escrita como
una ecuacin de pronstico en la forma
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 31
Esta ecuacin establece que la tendencia local de la vorticidad relativa est dada por la
adveccin de vorticidad absoluta. La resolucin numrica de esta ecuacin predice la
evolucin de la funcin corriente y por lo tanto de la vorticidad y los vientos. La figura 5.16
muestra un perfil tipico de velocidad vertical y divergencia horizontal y muestra que el
flujo en la tropsfera media es cercano a no-divergente en escalas sinpticas. Por lo tanto la
ecuacin anterior es un buen modelo para el pronstico del flujo en 500 mb el cual
determina en buena medida el movimiento de los centros de baja presion en superficie.
Figura 5.16 Perfiles tipicos de velocidad vertical y divergencia.
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 32
5.6 La aproximacin cuasi-geostrfica
En esta seccin realizaremos simplificaciones a las ecuaciones de momento, continuidad y
energa para desarrollar un sistema de ecuaciones mas simple que nos permitir entender la
naturaleza de la circulacin en los sistemas sinpticos de latitudes medias. Consideraremos
coordenadas isobricas para que no entre la densidad en las ecuaciones.
Como vimos anteriormente los sistemas sinpticos de latitudes medias estn en balance
geostrfico aproximado y en balance hidrosttico. La combinacin de estos balances da
lugar al viento trmico, que como vimos restringe los movimientos en latitudes medias.
Para dar un paso mas debemos considerar la evolucin temporal del flujo. Si despreciamos
la friccin las ecuaciones de movimiento aproximadas son
d V
dt
=f kV
la ecuacin hidrosttica
la ecuacin de continuidad
y la ecuacin de energa termodinmica
donde
S
p
=T
ln 0
p
. En lo que sigue condensaremos estas 5 ecuaciones en 2 que
satisfacern el anlisis de escala apropiado para escalas sinpticas.
Comenzamos considerando la expansin de f en la direccin meridional alrededor de una
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 33
latitud

0
(y=0 en esta latitud)
f =2Dsin
0
+2
D
a
cos
0
y= f
0
+ y
Para escalas sinpticas
lo cual justifica asignar al parmetro de Coriolis un valor constante f
0
excepto cuando es
diferenciado en el trmino advectivo. Esta aproximacin se denomina plano- pues
fsicamente se considera que los movimientos ocurren en un plano hipottico tangente a la
Tierra en la latitud

0
y restringidos a movimientos meridionales con una extensin
L<<a.
Separamos la veolcidad horizontal en los componentes geostrfico y ageostrfico
donde
Para sistemas sinpticos
o sea que Vg >> Va.
Por lo anterior podemos aproximar la velocidad real por la velocidad geostrfica a O(Ro) y
la razn de cambio de momento o temperatura siguiendo el movimiento horizontal puede
ser aproximado al mismo orden para la razn de cambio siguiendo el viento geostrfico.
Por lo tanto V puede ser sustitudo por Vg y la adveccin vertical que aparece solamente
debido a la existencia del viento ageostrfico puede despreciarse. As
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 34
d V
dt
=
d
g
V
g
dt
d
g
dt


t
+V
g
. =

t
+u
g

x
+v
g

y
La aceleracin es el resultado de la desviacin del viento con respecto a su valor
geostrfico. Por lo tanto no es posible sustituir simplemente V por Vg en el trmino de
Coriolis. En su lugar, los trminos de la derecha de la ecuacin de momento quedan
Entonces la ecuacin de momento aproximada queda
d
g
V
g
dt
=f
0
kV
a
y kV
g
(1)
y todos los trminos son de O(Ro) comparado con la fuerza gradiente de presin.
As definido, el viento geostrfico es no divergente por lo que la ecuacin de continuidad
queda
(2)
o sea que w est definida por la componente ageostrfica del viento horizontal.
En la ecuacin termodinmica la adveccin horizontal puede ser aproximada por su valor
geostrfico. Por otro lado, la estabilidad vertical en escalas sinpticas es suficientemente
grande como para que el calentamiento y enfriamiento adiabtico debido a movimientos
verticales sea del mismo orden que la adveccin horzontal de temperatura. El trmino de
calentamiento y enfriamiento adiabtico puede simplificarse dividiendo el campo de
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 35
temperatura en un estado base T
0
(p) que slo dependa de la presin mas una desviacin
T(x,y,p,t)
Como
solo la porcin debido al estado base se debe considerar al calcular el trmino de
estabilidad. Por lo tanto, usando la ecuacin hidrosttica, la forma cuasi-geostrfica de la
ecuacin de energa termodinmica se puede expresar como
(

t
+V
g
. )(

p
)co=

Q
p
(3)
donde
y
(0 es la temperatura potencial correspondiente al estado base de temperatura T
0
) y
c2.510
6
m
2
/ Pa
2
/ s
2
en la tropsfera media. Las ecuaciones (1), (2) y (3) junto con la
definicin de viento geostrfico forman el sistema cuasi-geostrfico de ecuaciones para
las variables geopotencial, viento geostrfico, viento ageostrfico y omega. No obstante
esta forma de las ecuaciones no es adecuada para el pronstico y es mas conveniente usar la
ecuacin de vorticidad en lugar de la ecuacin de momento pues solamente aparecer la
divergencia del viento ageostrfico.
El viento geostrfico lo escribimos como
por lo que la componente geostrfica de la vorticidad vertical es
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 36
Notar que la ecuacin anterior permite calcular la vorticidad geostrfica partiendo del
geopotencial, o vice versa. Esta invertibilidad es una de las razones por las cuales
vorticidad es tan til como diagnstico de prediccin. Como el laplaciano de una funcin es
mximo donde la funcin es un mnimo, vorticidad positiva implica valores bajos de
geopotencial. La ecuacin cuasi-geostrfica para la vorticidad se puede calcular partiendo
de las ecuaciones
d
g
u
g
dt
f
0
v
a
y v
g
=0
d
g
v
g
dt
+f
0
u
a
+ yu
g
=0
y tomando el rotor se obtiene
d
g

g
dt
= f
0
(
u
a
x
+
v
a
y
)v
g
Notando que
d
g
f
dt
=v
g
y sustituyendo la divergencia del viento ageostrfico usando la
ecuacin de continuidad podemos reescribir la ecuacin de la forma
la cual establece que la razn local de cambio de la vorticidad geostrfica est dada por la
suma de la adveccin de la vorticidad absoluta por el viento geostrfico mas la
concentracin o dilucin de vorticidad por el estiramiento o achatamiento de las columnas
de fludo (el efecto divergencia).
Resumiendo, en esta seccin derivamos 2 ecuaciones simplificadas de vorticidad y energa
que usaremos para el estudio de los sistemas sinpticos de latitudes medias. Si asumimos
que el calentamiento dQ/dt=0, podemos reescribir las ecuaciones como
Notas: Prof. Marcelo Barreiro
Introduccin a la Dinmica de la Atmsfera 37

t
(

p
)=V
g
. (

p
)+co

g
t
=V
g
. (
g
+ f )+ f
0
o
p
Notar que como la temperatura est directamente relacionada con
(

p
)
en superficies
isobricas, la primera ecuacin establece que la razn de cambio de la temperatura es la
diferencia entre la tendencia advectiva y el calentamiento o enfriamiento adiabtico debido
al movimiento vertical. (La aproximacin de dQ/dt=0 es vlida para los movimientos
sinpticos en latitudes medias, pero no para el estudio del desarrollo de ciclones donde el
calor latente juega un papel importante).
Como el viento geostrfico, la vorticidad geostrfica y pueden escribirse en funcin del
geopotencial las dos ecuaciones representan un sistema de dos ecuaciones y dos incgnitas
( , o ). En los captulos que siguen usaremos estas ecuaciones para calcular la tendencia
del geopotencial
(

t
)
y los movimientos verticales dados los campos instantneos de
geopotencial en varios niveles de la atmsfera. Esto permite diagnosticar y predecir el
comportamiento de la atmsfera de latitudes medias sin medir directamente el campo de
velocidades. Estas ecuaciones son la base de la meteorologa dinmica moderna.
Referencias principales
Mid-latitude atmospheric dynamics, J. Martin
An introduction to dynamic meteorology, J. Holton
Notas: Prof. Marcelo Barreiro