ENSO, ROTACIÓN TERRESTRE, VOLCANES Y SISMICIDAD

OSCAR MESA
Escuela de Geociencias y Medio Ambiente, Universidad Nacional de Colombia, Medellín.

Recibido para revisar 23 de Septiembre de 2001, aceptado Junio 2002, versión final Agosto de 2002.

RESUMEN: El fenómeno de El Niño- Oscilación del Sur (ENSO) integra paradigmáticamente todas las
geociencias. La interacción océano-atmósfera se acepta comúnmente como el actor central de esta oscilación
interanual del sistema climático. Algo más desconocida es la influencia del ENSO en la rotación terrestre,
aunque hay excelente evidencia empírica y sólida explicación teórica. Con respecto a la influencia del ENSO
en la sismicidad la situación es bastante polémica, las observaciones son discutibles y las posibles
explicaciones teóricas son apenas ideas sugestivas. Hay también autores que postulan que el ENSO obedece a
los efectos radiativos de aerosoles volcánicos en la estratosfera, involucrando de nuevo a la tierra sólida. En
este trabajo se hace una revisión crítica de la literatura científica relativa a estas interacciones y se aportan
algunos cálculos sobre la influencia en la sismicidad. La conclusión general es que tanto para la relación de
ENSO con la sismicidad como para el efecto del volcanismo en ENSO, las preguntas permanecen abiertas.
PALABRAS CLAVES: El Niño, Rotación Terrestre, Sismicidad, Vulcanismo.
ABSTRACT: The study of El Niño- Southern Oscillation (ENSO) is a paradigm for the necessary dialogue
of all the geosciences. Ocean-atmosphere interaction is readily accepted as the central actor of the interannual
oscillation of the climate system. The influence of ENSO in earth's rotation is not widely acknowledged
outside very specialized audiences, nevertheless empirical evidence in support is strong, and theory to explain
it is very well founded. ENSO influence on sismicity is even more polemic: observations are not clear and
theory is still under development. Some authors postulate that ENSO is a response to the effect of
stratospheric volcanic aerosol. This works is a review of the literature related to this issues and some
computations related to possible ENSO effect on seismicity are reported. The main conclusion is that issues
remain open both with respect to the relation of ENSO with seismicity and with respect to the relation of
volcanic aerosols and ENSO
KEYWORDS: El Niño, Earth rotation, Sismicidad, Volcanism

(ENSO por sus siglas en Inglés) significó la

1. INTRODUCCIÓN aceptación de la importancia de la interacción
entre la atmósfera y el océano. De manera
Por un extraño azar, el trabajo moderno sobre el tradicional, los meteorólogos toman las
clima implica el rompimiento de varias temperaturas superficiales del mar como
fronteras: Entre la hidrología, las ciencias de la condición de borde inferior, lo cual es
atmósfera y las del océano; entre las ciencias del justificado para escalas de tiempo muy cortas
clima y las ciencias de la tierra sólida; entre la pero ignora que los vientos juegan un papel
investigación aplicada y la investigación pura; importante en la determinación de éstas. Los
entre lo local y lo global; entre lo determinístico oceanógrafos de manera similar, toman los
y lo estocástico. La anterior lista incluye sólo las vientos como una condición de frontera superior
fronteras más importantes. Esta circunstancia impuesta externamente, sin tener en cuenta que
afortunada parece ser condición necesaria para los gradientes de temperatura son mecanismo
el avance y sin lugar a dudas es un reto difícil. fundamental para la producción de los vientos.
El fenómeno de El Niño Oscilación del Sur Esta interacción mutua es más clara e
 Mesa
42

importante en los trópicos y la zona ecuatorial y periódicamente entre dos y diez años. EN es una
a escalas de tiempo mensuales o mayores. amplificación del ciclo anual. Los cambios
Nuestra hipótesis es que sobre los continentes la oceánicos van acompañados de cambios
humedad del suelo juega un papel similar al de atmosféricos, la llamada Oscilación del Sur
la temperatura del mar y que su variabilidad es (SO), que corresponde a cambios en la masa
un factor no despreciable de la dinámica atmosférica entre el Este y el Oeste de la cuenca
climática. En estas condiciones tendríamos la Pacífica, con el correspondiente cambio en la
necesidad de integrar las miradas de las ciencias presión atmosférica, intensidad y dirección de
atmosféricas con las oceánicas y las los vientos y ubicación de las lluvias
hidrológicas. En este punto es conveniente predominantes. Estos cambios se reflejan a su
proporcionar un resumen condensado del vez en cambios en el nivel del mar y la
fenómeno ENSO. Para esto es necesario primero profundidad de la termoclina. Esta visión
explicar la condición normal del sistema océano- simplificada debe tomarse únicamente como una
atmósfera en el Pacífico tropical. Existe sobre el invitación a consultar estudios a mayor
Pacífico tropical una circulación térmica, profundidad (Diaz and Markgraff, 1992; Glantz,
llamada la circulación de Walker, con alta Katz and Nicholls, 1991; Hastenrath, 1991).
presión atmosférica en el Sudeste (Tahití) y baja Otros estudios confirman la relación causal
presión en Indonesia y el Norte de Australia entre la hidrología colombiana y los eventos
(Darwin). La diferencia de presiones produce ENSO (Poveda and Mesa, 1997). Tanto las
vientos superficiales del Este, que refuerzan los observaciones disponibles como el análisis
alisios del Noreste en el hemisferio Norte y del meteorológico y climático soportan esta
Sudeste en el hemisferio Sur y que se encuentran conclusión. Sin embargo, la relación no es lineal
en la zona de convergencia intertropical. Dichos y simple; tal vez una afirmación más correcta
vientos mueven el agua oceánica superficial sea decir que la variabilidad del clima
mediante esfuerzos de Reynolds. Como colombiano es parte integral de un ciclo
resultado se obtiene la corriente ecuatorial hacia aperiódico en la interacción océano-atmósfera-
el Oeste y las diferencias en la temperatura continente en la zona ecuatorial, caracterizada
superficial, nivel del mar y profundidad de la por una escala espacial global y una baja
termoclina entre el Este y el Oeste frecuencia (2 a 5 años) y cuya principal
(respectivamente 4ºC mayor, 40cm más bajo y manifestación es el ENSO. Aunque es bien
150m más profunda en el Oeste). Esta diferencia reconocido el efecto de la oscilación ENSO
de temperatura es la que produce la circulación sobre las lluvias, no se reconoce el hecho de que
de Walker en la atmósfera. El aire asciende esta relación es recíproca y que los efectos
sobre las aguas cálidas del Pacífico Oeste, fluye energéticos de la evapotranspiración, la
hacia el Este en la alta troposfera, y desciende condensación de vapor y la precipitación
sobre las aguas frías del Pacífico Este. La tropical son a su vez retroalimentadas y tienen
circulación de Walker es acompañada por efectos no despreciables sobre el resto del clima.
lluvias torrenciales en el Oeste y clima Pero la integración no termina allí. Los cambios
predominantemente seco en el Este. La en el momento angular atmosférico global
retroalimentación positiva, que acabamos de asociado con el fenómeno ENSO, aunque
describir, entre el océano y la atmósfera, puede pequeños en relación con el momento angular
eventualmente reversarse, lo cual sucede durante terrestre, producen cambios observables en la
los eventos ENSO. El Niño (EN) representa la rotación de la tierra y por tanto en la duración
componente oceánica de dicha oscilación, que del día. La razón está en que el momento
corresponde a un calentamiento anómalo de las angular total debe conservarse, excepto por el
aguas superficiales del Este de la cuenca del pequeño efecto de la fricción de las fuerzas de
Pacífico tropical. Normalmente las temperaturas mareas. La manera específica por medio de la
allí se incrementan durante febrero a marzo en cual se logran estos ajustes no está totalmente
aproximadamente 4ºC. Cuando hay un evento elucidada. Se sabe que la fricción, los cambios
EN este incremento es mayor y permanece por en el nivel del mar y las fuerzas sobre las
más tiempo. Estos eventos recurren cadenas de montañas asociadas a las diferencias
 Dyna 136, 2002 43

en la presión atmosférica son parte integrante del considerar la conservación del momento angular
balance. Hay elementos para pensar que estos como aplicable. A escalas mayores de tiempo el
ajustes pueden tener incidencias en la efecto de la fricción de marea no es despreciable.
sismicidad. Las ciencias de la tierra sólida Por ejemplo, una extrapolación de las
también entran en la integración. Además hay condiciones actuales hasta unos 1500 millones de
cierta evidencia del efecto de las emisiones años atrás, nos llevaría a que la luna estaría tan
volcánicas de aerosoles a la estratosfera tropical cerca de la tierra como para producir cataclismos
como uno de los factores en la dinámica de la que no se observan en el registro geológico
oscilación ENSO. (Williams, 2000). La reconstrucción de la historia
de la desaceleración de la tierra por fricción de
marea es todavía un gran desafío para la
2. ENSO Y ROTACIÓN TERRESTRE geología. A escala no tan grande, pero
definitivamente mayor de los 10 años, hay otros
2.1 Enso y rotación terrestre fenómenos geofísicos que afectan la rotación de
la tierra, por ejemplo acoplamiento en el manto,
La velocidad de rotación de la tierra "sólida" transferencias de masa por ciclos de glaciación y
exhibe pequeñas pero complicadas variaciones de rebote elástico por deglaciación (Peltier, 1998).
hasta varias partes en 100 millones, que Las componentes estacionales e
corresponden a variaciones de varios intraestacionales de la variación en la duración
milisegundos (ms) en la duración de un día e del día están bien explicadas en términos de la
incluso mayores cambios en la dirección del eje conservación de momento angular y en
de rotación (movimiento polar). Estos cambios consecuencia, los intercambios de momento
ocurren a varias escalas de tiempo que van desde angular entre la atmósfera y la tierra sólida
días hasta siglos, milenios y mayores. Lo que (Lambeck, 1988; Hide and Dickey, 1991). Los
refleja la gran variedad de fenómenos geofísicos factores meteorológicos difícilmente explican las
y astronómicos involucrados. En principio, para variaciones decadales o mayores. La pregunta
el sistema tierra-atmósfera-océano se conserva el interesante se refiere a la causa de la variabilidad
momento angular y por tanto los cambios en la interanual. La amplitud de esta componente
velocidad de rotación deben ir acompañados de (hasta 0,5ms) es comparable con la del ciclo
cambios compensatorios en alguna de las otras estacional. Se puede mostrar que los efectos
componentes del sistema. Estrictamente atmosféricos del ENSO son fundamentalmente
hablando, la fricción por fuerzas de marea los responsables de los cambios a esta escala de
produce un pequeño retardo en la rotación de la tiempo.
tierra que se manifiesta en un incremento en la La influencia de ENSO en la rotación terrestre
duración del día de aproximadamente

 por ha sido estudiada de tiempo atrás, la siguiente
año (-6x10-²²rad s-², Stacey, 1977), lo cual sólo es descripción sigue el trabajo de (Dickey et al.,
apreciable a escalas de tiempo milenarias o 1994). En efecto, El Niño de 1982-83 estuvo
mayores y que efectivamente se puede detectar acompañado por la mayor variabilidad interanual
del registro de eclipses por antiguos astrónomos registrada en la tasa de rotación de la tierra. La
(Kant, 1754; MacDonald, 1964; Munk, 1968; causa dominante de esta anomalía fue la
Williams, 2000). El origen de estas fuerzas está atmósfera: el momento angular atmosférico
en el torque que la luna y el sol ejercen sobre la (AAM) integrado desde 1000 hasta 1 hPa
protuberancia de la marea. Por efecto de la (troposfera y estratosfera) explica hasta el 92%
fricción, la tierra responde retrasada y la de la varianza interanual en la serie de longitud
protuberancia antecede al eje tierra-luna en un del día (LOD). Los vientos estratosféricos entre
ángulo pequeño, que da origen al torque que 100 y 1 hPa contribuyen a explicar casi al 20%
además acelera el movimiento de rotación de la de la varianza, lo que indica que la estratosfera
luna, lo que causa su alejamiento de la tierra. En juega un papel significativo en el balance de
lo que sigue nos concentraremos en fenómenos a momento angular terrestre a escala interanual. El
escala interanual (10 años o menos) y por lo tanto análisis de las series de momento angular
podemos descontar este efecto sin mayor error y atmosférico (AAM, las componentes ecuatoriales
 Mesa
44

H1 y H2, y la componente axial H3), la duración y también de su distribución con la latitud. El

del día (LOD,
   índice modificado de otro término tiene en cuenta los cambios en la
oscilación del Sur (MSOI, s) para un período de distribución de la masa, y es generalmente
15 años centrado en 1982-83 sugiere que esta menor.
tremenda anomalía en la rotación terrestre resultó La serie de duración del día (
  
de interferencia constructiva entre la componente 1 (Dickey et al., 1994) es del Jet Propulsión Lab.,
de baja frecuencia de ENSO (4 a 6 años) y la está filtrada (Kalman) y combina resultados
componente cuasi-bienal (2-3 años) del ENSO y basados en mediciones astronómicas y ópticas
de la QBO (siglas en inglés de la Oscilación (VLBI, interferometría de línea base muy larga, y
Cuasi-Bienal en los vientos zonales LLR, rangos lunares con láser). Después de 1985
estratosféricos tropicales, que cambian de la calidad es muy buena. La incertidumbre en
dirección entre el Este y el Oeste con un período 1976 es del orden de 0,09ms y de allí decrece
cercano a los dos años). Las fuentes de linealmente hasta 0,03ms en 1985. Las
discrepancia restantes para la serie de longitud del observaciones base para la construcción de la
día pueden estar en ruidos y errores sistemáticos serie están espaciadas cada 5 días. La
en las series de AAM y/o deberse a incertidumbre en la variabilidad interanual está
contribuciones oceánicas. estimada en 
La asociación entre los cambios en la
duración del día y los cambios en el momento
angular atmosférico ha sido objeto de varios
estudios. El avance más importante proviene de
los adelantos en la capacidad de observación
tanto atmosférica como astronómica. Para ilustrar
todas estas ideas consideremos las funciones
    i, i = 1, 2, 3, introducidas en
(Barnes et al., 1983) y que están relacionadas con
las dos componentes ecuatoriales y la
componente axial del momento angular, y que
permiten explicar estas relaciones. La
componente axial se define mediante la ecuación
1. Figura 1. Serie de tiempo de la duración del día y
sus distintas componentes (Dickey et al. 1994.)
R3 ps
2

X3 


u cos 2 dddp
2

gCm  0 
0 Las series de momento angular atmosférico

2
(Hi) tienen dos orígenes: El Centro Nacional de
0.7 R 4 2




ps cos 3 dd Meteorología Americano (NMC) y el Centro
2
(1)
gCm 0 
Europeo para Predicción Climática de Mediano
2
Plazo (ECMWF). En el primer caso van hasta los
Donde ( 1) son latitud y longitud 100 hPa y en el segundo hasta los 50 hPa. El
respectivamente, ps es la presión atmosférica NMC incorpora correcciones basadas en la
superficial y u es la componente de velocidad del hipótesis del barómetro invertido. También se
viento hacia el Este, R = 6,37x106 m es el radio usaron cálculos desde los 100 hasta los 1 hPa
medio de la tierra sólida, g = 9,81 m²/s es la haciendo suposiciones geostróficas a partir de
aceleración media de la gravedad, el coeficiente sondeos.
0,7 es el número de Love que tiene en cuenta las El índice usado para ENSO (s) es una versión
correcciones al tensor de inercia por la modificada del SOI (MSOI): diferencia entre las
"imperfección" de la rigidez terrestre. La presiones atmosféricas en Darwin y Tahití por tal
contribución dominante viene por el primer motivo la correlación con la serie de LOD es
término debido al viento, que depende de la positiva (durante los eventos El Niño el día es
intensidad del viento zonal, integrado en longitud más largo).
 Dyna 136, 2002 45

Hay excelente coincidencia de los componentes

interanuales ver Figura 2, (Dickey et al., 1994).
Durante un evento ENSO cálido se tiene un
máximo del MSOI (mínimo del SOI) y un
máximo de LOD. En estos casos hay un
debilitamiento generalizado de los vientos
ecuatoriales del Este y además fortalecimiento de
las corrientes chorro subtropicales del Oeste. Las
anomalías de los vientos tienen por lo tanto signo
positivo (hacia el Este). La conservación del
momento angular terrestre obliga a que la tierra
sólida tenga anomalías negativas (hacia el Oeste)
lo que se refleja en disminución de la velocidad
de rotación de la tierra y por tanto aumento de la
duración del día. Como se aprecia en la Figura 3
que muestra la Tierra vista desde un punto por
encima del polo Norte, señalando la dirección de
la rotación. El sistema de coordenadas tiene el eje
x a lo largo del ecuador hacía el Este. Las Figura 3. Esquema de la tierra vista desde arriba
anomalías positivas de momento angular del Polo Norte.
atmosférico (hacía el Este) corresponden a
debilitamiento en los vientos predominantes del Dickey et al., (1994) reportan que la máxima
Este y/o incrementos en vientos predominantes correlación (0,67) se da para un mes de rezago en
del Oeste, van acompañadas de aumento de la la serie de LOD con respecto al MSOI. Para el
duración del día, típico de el fenómeno de El evento 82-83, 0,50ms, con una raíz del
Niño y de la estación de invierno en el hemisferio error cuadrático medio (rms) igual a . La
Norte contribución estratosférica es importante, usando
los datos del Centro Europeo (viento y presión
hasta 100 hPa) da una reducción de rms de 193 a
  (73% de la varianza), incluyendo hasta 50
hPa da 84,6% de la varianza explicada y hasta 1
hPa permite explicar hasta el 91,6%, dejando un
residuo de . La estratosfera contribuye por lo
tanto con aproximadamente el 20% de la
varianza. Esto se explica teniendo en cuenta que
los vientos estratosféricos son mucho mayores
que los troposféricos (doble), al mayor brazo para
el momento y debido al cambio total de dirección
entre las distintas fases de la oscilación cuasi-
bienal. Es posible mejorar la correlación
Figura 2. Series desde 1964 hasta 1989 de la introduciendo factores de escala. Así se llega
componente interannual de los cambios en la hasta el 97,3% con un residuo de
. Aunque
duración del día y el índice modificado de la la estratosfera tiene sólo el 10% de la masa
oscilación del Sur, parte superior y Series desde atmosférica, contribuye al cambio de momento
1976 hasta 1991 de los cambios en la duración angular atmosférico interanual con un porcentaje
del día y las anomalías de momento angular significativo de la varianza. Sin embargo, las
atmosférico, parte inferior. (Dickey et al., 1994) conexiones entre QBO estratosférica y la QBO
troposférica permanecen elusivas.
 Mesa
46

2.2 Propagación de Anomalías de Momento fuertes de esos esfuerzos sobre el océano. De allí
Angular Atmosférico se deriva que el océano puede ser un
intermediario importante.
Dickey et al., (1992) estudian las observaciones y
demuestran que hay una propagación lenta,
coherente, a escala global e interanual, dirigida
del ecuador hacía los polos, de las anomalías de
momento angular atmosférico (AAMA). Estas
anomalías se originan en regiones ecuatoriales,
donde anteceden a las principales anomalías
atmosféricas del ENSO en casi dos años. De allí
se propagan hacía latitudes mayores de 60° en
ambos hemisferios, donde siguen rezagadas al
ciclo ENSO en aproximadamente cuatro años. Es
posible identificar dos componentes principales
en esta propagación, con patrones latitudinales
diferentes y períodos de 4,2 (baja frecuencia) y
2,4 (cuasi bienal) años, respectivamente. Los
episodios intensos de ENSO, en sus fases caliente
o fría corresponden a los casos en los que hay
interferencia constructiva. A escalas intranuales
se observa la importancia de oscilaciones de 40-
50 días.
La Figura 4 (Dickey et al., 1992) muestran la
coherencia cuadrada y el diagrama de fase -
coherencia entre las series de AAMA tropicales
de un lado y las de AAMA del hemisferio Norte y
el Sur del otro. La coherencia es significativa en Figura 4. (a) Coherencia cuadrada entre la serie
períodos centrados alrededor de los 4,2 años y tropical (de 7.5 ºN a 7.5ºS) y la subtropical (de
con las anomalías en los subtrópicos tendiendo a 15º a 31.4º de cada hemisferio) de anomalías de
ir rezagada de la de los trópicos en momento angular atmosférico en la parte superior
aproximadamente un cuarto de ciclo para y (b) el diagrama fase coherencia de las mismas
períodos mayores de aproximadamente 2 años. series en la parte inferior. (Dickey et al., 1992).
También se nota que los subtrópicos de los dos
hemisferios se separan a escalas del orden de 1,5
años, donde la coherencia es alta para el Esto también puede estar asociado a
hemisferio Norte que lidera en aproximadamente redistribución de masa     
5 meses. cambios en nivel del mar también pueden ser
muy significativos (si el aumento se debe a
2.3 Papel de los océanos expansión térmica no hay cambios rotacionales).
Dickey et al. (1994), apoyados en (Eubanks,
1993,) usando un modelo del océano en dos
El intercambio de AAM entre la atmósfera y capas estiman aproximadamente la contribución
la tierra sólida parece muy rápido a juzgar por los de los cambios en circulación, nivel del mar y
datos, pues el rezago es de apenas unos 10 a 20 profundidad de la termoclina. Se ha observado
días. Una porción significativa de este que en la zona ecuatorial el aumento de nivel del
intercambio puede ser por fricción y torque en las mar va acompañado de profundización de la
montañas. Hay alguna evidencia que a escala termoclina. Se estima que el efecto neto de un
interanual la contribución del esfuerzo cortante incremento de nivel del orden de 20cm, con la
sobre el océano puede ser significativa. ENSO en correspondiente profundización de la termoclina,
particular está asociado con fluctuaciones muy representa una disminución en la presión sobre el
 Dyna 136, 2002 47

fondo del mar del orden de 1,2 hPa. Esto del barómetro invertido. Esto justifica la
corresponde aproximadamente a lo ocurrido en necesidad de un modelo dinámico.
un evento El Niño como en 1982. Esto El océano puede considerarse estratificado en
corresponde a un incremento en la carga de dos capas, la más superficial de apenas unas
presión debido a la disminución de la masa cuantas centenas de metros y más caliente. Por
ecuatorial del océano del orden de 3x1014 kg, lo debajo hay una capa más fría, mucho más
que implica un incremento de LOD de 
. En profunda, de varios kilómetros de espesor. La
consecuencia, es muy probable que el resto de la termoclina, la frontera entre ambas capas
variación en LOD no explicada por cambios en representa un frente, con un cambio fuerte de
AAM venga del océano. temperatura, densidad y otros parámetros
Eubanks (1993) presenta un modelo dinámico asociados. Las oscilaciones libres de tal modelo
de las ondas oceánicas barotrópicas y baroclínicas se separan entre barotrópicas (densidad sólo es
para estudiar las posibles excitaciones oceánicas a función de la presión) y baroclínicas (densidad
la rotación terrestre. Las mediciones de la presión función de temperatura y/o salinidad, además de
en el fondo del mar sugieren fuertemente que en la presión). En la primera las superficies por
su variación está la causa de los movimientos encima de las dos capas oscilan en fase y en la
polares hasta ahora no explicados en cortos segunda fuera de fase. La velocidad de fase de las
períodos (semanas a meses). Estos cambios ondas baroclínicas es del orden de unos cuantos
obedecen a oscilaciones barotrópicas, forzadas a m-¹ y tardan centenas de días para cruzar el
su vez por esfuerzos superficiales de los vientos. Pacífico, mientras que las ondas barotrópicas se
A períodos mayores, las ondas baroclínicas que propagan a velocidades 100 veces mayor y sólo
están asociadas a la oscilación ENSO causan tardan unos cuantos días en cruzar el Pacífico. Un
cambios significativos en la tasa de rotación. Por simple modelo lineal de dos capas puede servir
contraste, la alta correlación entre la tasa de para representar la dinámica.
rotación y el momento angular atmosférico Sea z la coordenada vertical, con origen en la
indican que hay un rápido mecanismo de superficie no perturbada entre el océano y la
transferencia de las perturbaciones de los vientos atmósfera. Sean H1 y H2 las profundidades no
al océano y de allí a la tierra sólida. La evidencia perturbadas de cada capa. El fondo del mar está
indica que esto se logra mediante rápida por tanto a una profundidad H = H1 + H2. Sean
propagación zonal de ondas barotrópicas h1 y h2 las perturbaciones en H1 y H2 y suponga
oceánicas hasta las márgenes de los continentes. que cada capa tiene una densidad constante e
Como no hay registros completos y igual a 1 y 2 respectivamente. Sea 2 –1
suficientemente largos de presiones en el fondo 1 . El eje x tiene dirección hacía el Este, ui y
del mar es necesario recurrir a modelos. Una vi son las velocidades correspondientes en la capa
posibilidad es trabajar directamente con la i. La velocidad angular de rotación de la tierra se
presión atmosférica a nivel del mar e ignorar lo representa por [Link] ecuaciones de momento
que sucede en el fondo del océano, otra linearizadas son las siguientes:
posibilidad es el uso de la hipótesis del barómetro
invertido, que usa únicamente la presión 
atmosférica. En este modelo la superficie del mar u1  v1   x
v   20  u  
y 1
permanece como una superficie equipotencial gh (2)
bajo los efectos de la variación en la presión t  1  1   
 
atmosférica. La presión en el fondo se compensa  
con cambios, de signo contrario, en el nivel del
mar. La hipótesis del barómetro invertido no tiene

en cuenta los efectos dinámicos y la influencia de u2   v2   x
los esfuerzos del viento en la superficie. v  2 0  u 2   g
h1 1 h2  /  2 
t 2   y
Eubanks et. al., (1988) han mostrado que
   
ninguno de los métodos anteriores puede explicar  
satisfactoriamente las observaciones del (3)
movimiento polar, siendo menos mala la hipótesis
 Mesa
48

Y las ecuaciones de continuidad p  1 gh1 gh2 con p f   2 gh1 (7)

f

h1 h  u v  p f   1 gh1
 2  H1  1  1 (4)
t t  x y 
H1 1  H (8)
donde  
H 2 2
h2  u v 
 H 2  2  2 (5)
t  x y 
Con estas presiones es posible calcular la
contribución a los cambios en el momento
Como las ecuaciones son lineales es posible angular axial, según el último término de la
expresar la respuesta a forzamientos externos expresión para 3 dada por la ecuación (1). Con
mediante superposición con los modos valores realistas de los parámetros,  ,
principales y también es posible hacer separación H  , H1   se obtiene -0,06, y por
de variables. Los modos principales dan: lo tanto un incremento baroclínico de 10 cm en el
nivel del mar implica un decrecimiento de la

 Ondas de gravedad que son muy rápidas en presión en el fondo del mar de 0,6 hPa.
comparación con la rotación terrestre, para Hay ausencia de mediciones y registros
ellas el término de Coriolis se puede sistemáticos de la presión en el fondo del mar.
despreciar. Algunas campañas como el experimento

 Ondas Rossby, que son muy lentas en BEMPEX (Barotropic Electromagnetic and
comparación con la rotación terrestre y por Pressure Experiment) por 10 meses en el centro
tanto para ellas los términos
i/
t y
i/
 se del Pacífico Norte, reportan cambios de 5 hPa,
pueden despreciar. que son comparables a los cambios en la presión

 Ondas Kelvin, atrapadas en el ecuador, para atmosférica. También hay algunos datos de
las cuales vi es cero. PMEL (Laboratorio ambiental marino del
Pacífico).
Todas estas ondas tienen la misma estructura Se había señalado que a escala interanual, los
vertical, aunque sus velocidades de propagación y cambios en el momento angular atmosférico
estructura horizontal sean muy diferentes. La asociados a los eventos ENSO explican buena
solución a la estructura vertical se puede parte de la variabilidad en la duración del día,
encontrar introduciendo una función propia quedando un residuo de
 y que los cambios
escalar , tal que oceánicos podrían explicar este residuo. En
efecto, durante El Niño hay cambios en el nivel
h2 h1 (6) del mar de decenas de centímetros en el Pacífico
ecuatorial. Estos cambios son producidos por los
Teniendo en cuenta que H1<<H2 y que cambios en los vientos superficiales.
i es posible encontrar dos autofunciones Generalmente se ha supuesto que los cambios en
para solución del problema planteado, con el nivel del mar son estéricos, es decir por
soluciones expansión térmica. Si este fuera el caso, no habría
lugar a ningún efecto rotacional. En realidad,
H2 H 2 estos cambios en el nivel del mar se producen por
 1 y   redistribución de aguas cálidas a gran escala, que
H H 2  pueden modelarse como ondas baroclinas en la
parte superior del océano, y por lo tanto pueden
que corresponden al modo barotrópico y excitar variaciones en la rotación de la tierra.
baroclínico respectivamente. En cada caso se Wyrtki (1985) mostró la alta correlación
puede calcular la perturbación de la presión sobre entre niveles del mar y profundidad de la
el fondo del océano pf debida a los cambios en el termoclina en el Pacífico tropical. Con base en
nivel de la superficie y la profundidad de la esto estimó el cambio en el volumen total de agua
termoclina en el Pacífico tropical, entre 15° de latitud Norte
 Dyna 136, 2002 49

y Sur, en función del tiempo, para el evento de criticalidad, al borde de la falla, pequeños
1982-83. Su estimación indica que entre perturbaciones se amplifican y producen grandes
mediados del 82 y mediados del 83 el volumen cambios. La conexión se apoya en una mirada a
disminuyó aproximadamente en 5•10¹² m³, los mecanismos de transferencia de momento
reflejando pérdida de agua que se propagó a los angular entre la atmósfera, los océanos y la tierra
polos por las costas. Aunque posteriormente ha sólida. El numeral anterior mostró como la
habido revalorizaciones de estos cálculos, siguen variabilidad interanual del clima explica bastante
siendo correctos en lo fundamental. bien los cambios a escala interanual en la
Usando el modelo de dos capas,   -0,06 y rotación de la tierra, representados por la serie de
los datos de Wyrtki se obtiene un incremento en duración del día. No se entró en detalles, pero a
la carga de presión sobre el fondo, por la escala intra-anual sucede algo semejante (Oort,
disminución de la masa ecuatorial de 3•1011 kg, 1989; Peixoto and Oort, 1992). La atmósfera rota
que equivale a un incremento en la duración del más rápido que la tierra sólida durante el invierno
día de , si esta agua viene del Norte y Sur de del hemisferio Norte y viceversa en el verano. La
los 45°. Aunque estas estimaciones son sólo del corriente de chorro subtropical que se intensifica
orden del 30% del residuo no explicado, son del en invierno es la principal responsable de este
signo correcto y ocurren en el momento correcto. ciclo. En consecuencia el día es más largo
Note que este efecto es integrado sobre todo el durante el invierno del hemisferio Norte, por
Pacífico tropical, y que hay diferencia de signo conservación de momento angular. La atmósfera
entre los cambios que ocurren en el Este y en el que tiene aproximadamente una masa apenas 10-6
Oeste. veces la masa de la tierra sólida, es sin embargo
Por último, aunque por falta de espacio no se responsable por estos cambios en la duración del
entró en detalles sobre las celeridades de las día. Esto se logra mediante intercambios de
ondas que se obtienen del modelo de dos capas, momento angular, torques entre la atmósfera y la
esta idealización permite obtener la información superficie. La pregunta es si las fuerzas generadas
básica sobre estos fenómenos oceánicos, de en este intercambio juegan algún papel en la
fundamental importancia para entender la generación de algunos sismos.
dinámica de ENSO, que evidentemente requiere Los mecanismos por medio de los cuales
incorporar la parte atmosférica y la intercambian momento angular la atmósfera y la
retroalimentación entre ambas componentes. tierra sólida son básicamente fuerzas de fricción
sobre la superficie (océano y continente), torques
en las montañas por diferencias de presión entre
3. CLIMA Y VARIACIONES EN sus laderas Este y Oeste, torques sobre los
continentes por diferencias en el nivel del mar
SISMICIDAD
entre la costa Este y Oeste, que son ocasionados
por los esfuerzos superficiales de los vientos
La teoría de la tectónica de placas es aceptada sobre el océano. Además hay que considerar que
hoy para explicar la sismicidad, los mecanismos los vientos superficiales son predominantemente
para la generación de sismos y las diferencias de del Este en los trópicos y del Oeste en los
actividad sísmica entre regiones (England, 1992). extratrópicos, por lo que hay un patrón
Sin embargo, conviene tener en cuenta que determinado de propagación de momento entre
trabajos muy recientes (Bostrom, 2000) postulan los trópicos y los extratrópicos.
que la actividad tectónica de la tierra es modulada La fricción tiende a frenar tanto los vientos
por la rotación de la tierra y las mareas, que del Este como los del Oeste, de tal manera que la
organizan el proceso de convección en el manto. superficie en los trópicos es una fuente de
La reflexión que sigue usa la tectónica como momento angular del Oeste (hacía el Este) y en
punto de partida para postular algunas posibles las latitudes medias es un sumidero. El momento
relaciones de la variabilidad de la sismicidad en se transporta a la parte alta de la atmósfera
el tiempo, con la variabilidad climática, mediante las circulaciones meridionales (Hadley
especialmente a escala interanual. La teoría y Ferrel), pero al nivel de la corriente de chorro
moderna del caos ha señalado cómo cerca de la en las latitudes medias el momento se transporta
 Mesa
50

por medio de las vaguadas y riscos de presión que que pueden influir en la variabilidad temporal de
tienen una inclinación SW-NE en el hemisferio la sismicidad.
Norte y SE-NW en el hemisferio Sur. De todas
maneras, la transferencia de momento en cada
banda de latitud ocurre por fricción y torque en
las montañas. En el largo plazo la diferencia de
niveles del mar entre las costas Este y Oeste de
los continentes es también importante.
Un posible escenario es que los continentes
liberen esfuerzos intermitentemente mediante
movimientos internos, tales como los que ocurren
a lo largo de las fallas. En tal caso es de esperar
que las líneas de falla tengan una orientación
acorde, tanto en las latitudes medias de cada
hemisferio como en los trópicos, lo cual se puede Figura 5. Series del Indice de sismicidad en la
verificar en un mapa. Por ejemplo la falla de San parte noroccidental de América del sur y del
Andrés en California, o las fallas en Nueva índice de Oscilación del Sur.
Zelanda, o en el borde Pacífico de Las Américas.
Además se espera una variabilidad en la Sí es mucho más clara la existencia de un ciclo
sismicidad debido a la ocurrencia de eventos anual en la sismicidad, siguiendo el ciclo
como El Niño, que tienen asociado cambios correspondiente en el momento angular
fuertes en el momento angular en la dirección del atmosférico (Figura 6). Allí se representa para el
eje de rotación. Incluso se deben dar cambios con mismo catálogo de temblores, el análisis de la
la época del año. distribución del número de sismos según la
Algunos autores han calculado el efecto de magnitud, separando en cuatro grupos, según la
grandes sismos en el momento angular en las época del año de su ocurrencia. Para todas las
direcciones perpendiculares al eje de rotación y magnitudes es mayor el número de sismos en
que afectan la posición geográfica del Polo Norte. diciembre-enero-febrero y menor en septiembre-
Incluso se ha argumentado que los sismos son octubre-noviembre, siendo intermedia para los
mecanismo para excitar el llamado bamboleo de meses restantes. Para las magnitudes altas la
Chandler (Mansinha y Smile, 1967; Shimazaki diferencia es mayor, llegando a ser más del doble
and Takeuchi, 1972; Preisig, 1992; Chinnery and el número de temblores de magnitud mayor que
Wells, 1972). 6.5 en diciembre-enero-febrero que en
Con el ánimo de explorar la relación a partir septiembre-octubre-noviembre.
de algunos datos, la Figura 5 tiene dos series, una
con la variación en el siglo XX de un índice de la
sismicidad en la parte occidental y más ecuatorial
de Centro América y América del Sur,
específicamente en la región comprendida entre
latitudes 10°N y 10°S, y longitudes 70°O y 80°O.
El índice se calcula con base en la energía
liberada por los sismos de magnitud mayor a 5
(datos tomados del U.S. Geological Survey,
National Earthquake Information Center). La otra
serie corresponde al Índice de Oscilación del Sur,
que representa bien la ocurrencia (valores
negativos) del fenómeno de El Niño o de la Niña
(valores positivos). Aunque pueden haber algunas
épocas de buena coincidencia, no se puede Figura 6. Ciclo anual de la distribución del
apreciar unívocamente una relación directa. Esto número de sismos para una determinada
era de esperarse pues hay muchas otras causas magnitud en la misma región.
 Dyna 136, 2002 51

A escala global, la Figura 7 tiene tres series, eventos El Niño. Walker sugiere que hay
en la primera se presenta el número de sismos precursores sísmicos de los eventos El Niño y
superficiales (profundidad menor a 100 km) con para justificar su afirmación da algún soporte con
magnitud mayor o igual a 5, en el Catálogo del observaciones en los más recientes eventos. Ya
Centro Sismológico Internacional. En la segunda que no hay una identificación clara de los
se muestra el Índice de Oscilación del Sur mecanismos que desencadenan El Niño, y con
Modificado (negativo del SOI) y en la tercera se base a la anterior observación, Walker postula
presenta la variación en la duración del día que el incremento en la sismicidad está
(LOD). Las bandas verticales muestran la probablemente asociado a la intensificación del
ocurrencia de eventos El Niño. (Chinnery, M, flujo de calor en la dorsal, lo que explicaría
1991). además el inicio de El Niño.
En contra de la hipótesis de Walker hay
varios argumentos. En primer lugar, el papel de
las anomalías en los vientos atmosféricos como
detonante de El Niño está bastante bien aceptado.
En segundo lugar, parece tomar cuerpo un
mecanismo extratropical en la dinámica
atmosférica y oceánica para desencadenar las
anomalías atmosféricas que producen El Niño
(Yu and Rienecker, 1998; Hunt, 2000). Por lo
tanto la correlación geográfica de Walker entre el
lugar donde ocurre la sismicidad y eventual
aumento de liberación de calor en la dorsal, y los
cambios en la superficie del mar tienen que
recorrer un camino muy largo para ir hasta el
Pacífico Occidental y/o los extratrópicos y
regresar.
Subsiste la pregunta sobre las observaciones
de Walker: ¿Son estas correlaciones producto de
la casualidad, o están asociadas a un mecanismo
Figura 7. Serie del número de sismos físico, debido a aumento en los esfuerzos en el
superficiales (profundidad menor a 100 km) con fondo del mar por los cambios tanto en presión
magnitud mayor o igual a 5. En el II cuadro se atmosférica como en presión "hidráulica" por los
muestra el Índice de Oscilación del Sur cambios de nivel y profundidad de la termoclina?
Modificado (negativo del SOI) y en el III cuadro El tema parece interesante. Para mirarlo con
se presenta la variación en la duración del día algún rigor habría que desarrollar un modelo
(LOD). Las bandas verticales muestran la elastoplástico de la corteza terrestre y someter las
ocurrencia de eventos El Niño. (Chinnery, M., placas a los esfuerzos que se derivan de los
1991) modelos atmosféricos y oceánicos. Tampoco es
conveniente dar lugar a especulaciones poco
El análisis de estas tres figuras da lugar a científicas, basadas en temores y en la idea
especulaciones acerca de la existencia de alguna errónea de que todo tiene que ver con todo, sin
relación. cuantificar y sin demostrar la importancia.
Walker (1988, 1995 y 1999) ha presentado Además, no hay que ignorar que los
alguna evidencia de aumento de sismicidad en el catálogos pueden tener cambios en sismicidad en
Este del Pacífico tropical, en la zona de más el tiempo debido a cambios en la manera de
rápida expansión del suelo marino en la Dorsal observar y no a la sismicidad misma. Habermann
del Pacífico Este, durante los eventos El Niño. (1991) presenta una manera interesante de
Esta zona está ubicada directamente debajo de la verificar la existencia de estas perturbaciones
zona de alta presión atmosférica predominante, artificiales usando la ley de Ritcher de
que se altera fuertemente con la ocurrencia de los distribución de magnitudes. Obtiene en particular
 Mesa
52

que para las Américas hay efectos importantes de de la atmósfera y la superficie terrestre
esta naturaleza aproximadamente hacía 1970. De experimentan una pérdida neta de radiación solar
allí en adelante el registro es fundamentalmente por dispersión hacía atrás (lo que aumenta el
homogéneo. albedo) y además una redistribución de la
absorción de la radiación solar. Mucha de la
radiación solar se dispersa hacia delante,
4. VOLCANES, CLIMA Y ENSO produciendo un incremento en la radiación difusa
hacia abajo, lo que compensa parcialmente la
4.1 Generalidades disminución en el rayo solar directo. La Figura 8
ilustra esquemáticamente los distintos impactos.
Plutarco señaló que la erupción del Etna en el A simple vista estos efectos se aprecian en
año 44 a.C. opacó el sol, y atribuyó al cambios de color en el cielo. Por ejemplo el azul
enfriamiento resultante, la falla en las cosechas y celeste normal se vuelve lechoso por el aumento
la hambruna en Roma y Egipto. Benjamin de la radiación difusa. La reflexión del sol
Franklin sugirió que la erupción del Lakagigar en poniente en las capas inferiores del aerosol
Islandia en 1783 fue la responsable de que en estratosférico producen los característicos
Europa, en ese año, el verano y el subsecuente atardeceres rojizos (el famoso cuadro "El Grito"
invierno fueran anormalmente fríos. Algunos le de Munch, 1892, muestra un atardecer en Oslo
atribuyen un papel a la megaerupción del volcán después de la erupción del volcán Awu). La
Toba, 71000 años a.C. en la glaciación instrumentación moderna ha permitido
subsiguiente. Revisiones detalladas sobre el determinar con mayor detalle el efecto radiativo
efecto climático de las erupciones volcánicas se de los aerosoles a diferentes alturas y longitudes
encuentran en Lamb (1970), Toon and Pollack de onda. Por ejemplo, en la Figura 9 se muestra el
(1980), Toon (1982), Ellsaesser (1983), calentamiento radiativo para tres longitudes de
Kondratyev (1988), Robock (1989, 1991 y 2000), onda, en función de la latitud y la altura, para el
y Kondratyev and Galindo (1997). caso del Pinatubo. Se observa un calentamiento
El material de las erupciones volcánicas está en la parte superior de la nube de aerosol, por
compuesto principalmente por cenizas y gases de absorción en el infrarrojo cercano. Este efecto
dióxido de azufre, ácido clorhídrico, vapor de domina al incremento en el enfriamiento en el
agua y H2S. El efecto de la ceniza en el clima infrarrojo debido al aumento en la emisividad por
local puede ser importante, pero no en el clima la presencia de aerosoles. En la estratosfera baja,
global; porque la ceniza dura muy poco en la la atmósfera se calienta por absorción de
atmósfera y no se extiende a grandes áreas. radiación de onda larga de la troposfera y la
Si la erupción de gases penetra solamente la superficie. También se producen cambios en la
troposfera su efecto climático también es menor, microfísica de las nubes, la química de la lluvia,
pues se lavarán en cuestión de meses las tormentas eléctricas y en la convección
(Ramanathan, 1988). Si la fuerza de la erupción profunda como resultado de la abundancia de
es suficiente para penetrar a la estratosfera la núcleos de condensación. Estos aerosoles tienen a
situación es completamente diferente. Aunque la su vez un impacto severo sobre la química de la
ceniza se precipita en los primeros meses, la estratosfera, en particular sobre el ozono. De
estabilidad vertical de la atmósfera en esta capa hecho, varios autores (Crutzen, 1976; Tie et al.,
trae un efecto mucho más duradero de los gases, 1996; Solomon et al., 1996 y Solomon, 1999) han
que se dispersan horizontalmente debido a señalado que el agujero en la capa de ozono está
cortante. El dióxido de azufre reacciona con el fuertemente relacionado con la actividad
agua y produce ácido sulfúrico, que a la volcánica natural, además de la actividad
temperatura prevaleciente (-50° C a -80° C) se antrópica. Sobra reconocer que a su vez, la
condensa para formar aerosol de tamaños disminución de la concentración del ozono trae
inferiores al micrómetro (). Estas partículas efectos sobre el equilibrio radiativo, por su
son muy eficientes para dispersar la luz solar que capacidad de absorber radiación ultravioleta.
tiene longitudes de onda de ese tamaño, En el pasado se creía que la cantidad de gases
aproximadamente. En consecuencia, la parte baja liberados por una erupción volcánica estaba
 Dyna 136, 2002 53

directamente relacionada con la cantidad de Algunos volcanes de escudo pueden producir

ceniza. Pero recientemente se ha encontrado que grandes cantidades de ceniza y relativamente
está más relacionada con la química del magma y poco gas. Mientras que hay casos como el de
que la relación entre cantidades de ceniza y gas Nyamiragira, Africa, en 1981, que emitió
puede variar por factores del orden de 100 a 1000 enormes cantidades de gas y prácticamente
(Palais and Sigurdson, 1989). ninguna ceniza (Smithsonian, 1981). El gas puede
ser invisible lo que dificulta la interpretación de
los registros históricos. Además, la altura de la
pluma de ceniza puede ser muy diferente de la
altura de la pluma de gas.
Los datos para El Chichón, México, 1982,
indican que la retro-dispersión de la radiación
solar fue en las primeras semanas entre el 5 y el
10%. En el caso del Pinatubo, se estima que la
disminución de radiación incidente en la
superficie en zonas tropicales fue de 8 Wm-2.
Con el tiempo los aerosoles se dispersan
hacía los polos, con mayor transporte en el
invierno en cada hemisferio. Este proceso puede
tomar hasta 2 años. Es posible que haya cruce del
Figura 8. Esquema de las emisiones volcánicas a ecuador con el cambio de las estaciones, todo
la atmósfera y sus efectos, (Robock, 2000). dependiendo de la latitud de la erupción. El
efecto climático de una erupción también es por
lo tanto función de su latitud.

4.2 Procedimientos de Observación

La siguiente es una enumeración de diversos

procedimientos de observación que pueden
emplearse para estudiar el efecto climático de las
erupciones volcánicas:

 Medición del rayo solar directo. Desde 1883

hay observaciones en superficie del rayo solar
directo mediante radiómetro. En 1924 había 8
estaciones, todas al Norte de la latitud 40° N,
con calidad variable y no estandarizada. La
falta de observaciones tropicales y en el
hemisferio Sur, sumada a las dificultades de
calidad de la información le restan
aplicabilidad a estas observaciones. Sin
embargo hay rasgos claros de algunas
erupciones en este registro.

 El radiómetro TOMS (Total Ozone Mapping
Spectrometer, Espectrómetro para mapear el
Figura 9. Calentamiento radiativo para tres ozono total) a bordo del satélite Nimbus-7
longitudes de onda, en función de la latitud y la lanzado en 1979, puede ser usado para
altura, para el caso del Pinatubo, (Robock, 2000). detectar el gas SO2 emitido por los volcanes,
aunque su diseño original era las bandas de
 Mesa
54

absorción ultravioleta del ozono. Por ejemplo, depositados y la datación se hace con base en
el Pinatubo en 1991 puede haber inyectado O18/O16.
hasta 20•109 kg de SO2 a la estratosfera (Bluth
et al., 1992; Read, et al., 1983). 4.1 Efectos Observados en el Clima

 El radiómetro de alta resolución AVHRR Regional y Global
(Advanced Very High Resolution
Radiometer) a bordo de los satélites En términos generales se ha argumentado que las
NOAA/11 y siguientes mide la dispersión de erupciones volcánicas producen reducción del
la luz, pero sólo sobre mar (Stowe et al., ciclo diurno, enfriamiento en verano, reducción
1992). de precipitación en los trópicos, calentamiento

 Otra técnica satelital usa el instrumento estratosférico y calentamiento en invierno,
SAGE II (Stratospheric Aerosol and Gas fenómeno ENSO y la pequeña edad de hielo de la
Experiment) a bordo del satélite ERBS (Earth edad media.
Radiation Budget Satellite). El instrumento La evidencia más clara es la disminución de
observa la atenuación de la luz solar en los temperatura. Lamb (1970) discute la evidencia
atardeceres o amaneceres, es decir desde 1500 hasta 1968. En general, las mayores
observando el sol en el borde de la atmósfera. erupciones han significado un enfriamiento
El resultado es la profundidad óptica de la global de la baja atmósfera. En la Figura 10
atmósfera. Se ha encontrado clara evidencia tomada de Handler (1993) que contiene la
de las mayores erupciones volcánicas evolución de la temperatura en el período
(McCormick and Veiga, 1992; Veiga, 1993). histórico se puede observar este efecto, en

 Otra técnica de medición usa el lidar (láser- particular para los volcanes Kraratoa, Agung y
radar). Se usa un pulso de láser dirigido a la Pinatubo. Note que menor radiación no
estratosfera, se mide la amplitud y el tiempo necesariamente produce enfriamiento local. La
de llegada del rayo reflejado para determinar energía se conserva, pero no la temperatura. La
la dispersión producida por el aerosol energía se distribuye en el cambio de temperatura
estratosférico (DeFoor et al., 1992; Jäger, del hielo, agua, vapor, aire, tierra, vegetación, en
1992). energía cinética del aire y el océano.

 Registro de erupciones. Mediante testimonio
directo se ha reconstruido un catálogo de las
mayores erupciones en época reciente. Debe
tenerse en cuenta que es posible que haya
erupciones con aerosoles importantes que
pasen inadvertidas, por no estar cerca de
sitios poblados o por no estar acompañadas de
emisiones significativas de ceniza. Los
índices de las erupciones volcánicas más
empleados son: el DVI (Dust Veil Index),
índice de velo por polvo, basado en
observaciones de los atardeceres (Lamb,
1970); el VEI (Volcanic Explosivity Index),
índice de explosividad volcánico, basado en
registros geológicos e históricos (Newhall and Figura 10. Evolución de la temperatura global
Self, 1982) y el IVI (Icecore Volcanic Index), promedio desde 1880 hasta 1990, suavizada con
índice de vulcanismo de núcleos de hielo un promedio móvil de 5 años. Se puede observar
(Robock and Free, 1995). el efecto de enfriamiento por erupción volcánica,

 Otro estudio para deducir la historia del en particular para los volcanes Kraratoa, Agung y
vulcanismo es mediante núcleos de hielo de Pinatubo. (Handler, 1993).
glaciares en Groenlandia, Antártica u otros
glaciares en nevados tropicales. En el hielo se
determina la concentración de aerosoles
 Dyna 136, 2002 55

4.1.1 Relación con lluvias del monzón hindú. calentamiento en Norte América, Europa y
Siberia. Este patrón de temperaturas está
Handler (1993) encuentra alta correlación entre la directamente relacionado con la circulación
ocurrencia de erupciones volcánicas con troposférica, dominada por el chorro polar
aerosoles en latitudes bajas y las lluvias en el nocturno, que ocurre simultáneamente con un
subcontinente. Observa lluvias inferiores a lo patrón de baja presión sobre Groenlandia y de
normal después de eventos volcánicos con alta presión sobre Norte América, Europa y el
inyección de aerosoles en latitudes bajas, y Este de Asia. Este patrón parece ser un modo de
superiores a lo normal después de erupciones en vibrar natural de la atmósfera en invierno que se
latitudes altas. puede forzar por las perturbaciones
estratosféricas asociadas a las erupciones
4.1.2 Relación con ENSO. volcánicas. La explicación teórica es que el fuerte
chorro polar nocturno atrapa las ondas planetarias
Handler (1993) dice que las erupciones que se propagan verticalmente, lo que interfiere
volcánicas con alta inyección de aerosol constructivamente para producir el patrón
estratosférico en latitudes bajas desencadenan mencionado (Robock, 2000). El chorro
eventos ENSO cálidos. Para esto presenta alguna intensificado evita los súbitos calentamientos
evidencia (muy débil) desde el punto de vista estratosféricos tarde en el invierno, por lo que se
estadístico. Ante la falta de predecesores de perpetua él mismo, manteniendo la baja
algunos eventos ENSO muy intensos (1899-1900, atmósfera polar fría en el vórtice polar aislado.
1925-26, 1939-41 y 1957-1959) hace una revisión Mientras que cuando este patrón se presenta de
en los registros de hielo en Groenlandia, manera natural, sin forzamiento externo, es
aportando algo de evidencia. frecuente que ceda a los súbitos calentamientos
estratósfericos. Este patrón ha seguido a todas las
4.1.3 Correlación con espesor de anillos de erupciones tropicales desde que se tienen
árboles. observaciones con radiosondas. Debe sin
embargo notarse que otros forzamientos pueden
La Marche y Hircshboeck (1984) reportan alta producir el mismo patrón, entre ellos el ciclo
correlación entre ocurrencia de erupciones y solar de 11 años, el calentamiento global y el
registros de anillos en los Estados Unidos. agujero en la capa de ozono. El ciclo solar tiene
efecto en la temperatura de la alta estratosfera, el
Varios científicos, entre ellos Robock (1991) y flujo de radiación ultravioleta y la cantidad de
Robock & Mao (1992) han simulado el efecto de ozono que se combinan para producir
perturbaciones radiativas producidas por calentamiento y circulación anómala en la
aerosoles debidos a erupciones volcánicas. Los estratosfera (Kodera et al., 1991; Robock, 1996;
modelos indican que es posible que algunas Haigh, 1996; Shindell et al., 1999; Labitzke &
regiones se calienten y otras se enfríen. Por van Loon, 1988; van Loon & Labitzke, 1990). El
ejemplo en 1991-92, después de Pinatubo, el calentamiento global aumenta más el espesor de
invierno en Norte América fue muy suave casi en la troposfera cerca al ecuador que en los extra-
todos los lugares excepto en el medio Oeste. trópicos e intensifica el gradiente de temperatura
Robock (2000) presenta evidencia en el sentido del ecuador a los polos (Graf et al., 1995). El
que este calentamiento en invierno en ciertas agujero en la capa de ozono produce mayor
regiones del hemisferio norte está asociado al enfriamiento cerca de los polos (Ramaswamy et
patrón de la oscilación del Atlántico Norte al., 1995). Los estudios mediante modelos
(NAO). Hay observaciones de inviernos suaves globales tienden a confirmar esta teoría (Robock,
en Rusia asociados a las erupciones volcánicas 2000). Cabe anotar que en otras investigaciones
(Groisman, 1985). En Norte América la evidencia hemos evidenciado la influencia de la oscilación
no es tan clara, pues hay mezcla con el efecto de NAO (Oscilación del Atlántico Norte) sobre el
ENSO y tal vez alguna diferencia regional. De clima de Colombia (Mesa, Poveda y Carvajal,
hecho el patrón incluye enfriamiento en Alaska, 1997).
Groenlandia, el Medio Este y China, y
 Mesa
56

Rind et al. (1992) usando el modelo global de Los vientos alisios son máximos en Julio y
clima de la NASA encontraron que un incremento mínimos en Diciembre. Es claro que las
del aerosol estratosférico tropical produce temperaturas del mar están inversamente ligadas
calentamiento de la baja estratosfera, lo que a los vientos con un rezago de dos a tres meses.
disminuye el gradiente vertical de temperatura La razón para esto está en los procesos de
tropical, lo que a su vez debilita la circulación de evaporación, mezcla con aguas más profundas y
Hadley (del orden de 1 ms-1). Esto puede producir más frías, y surgencia de aguas profundas; todos
eventos ENSO. La perturbación usada es parecida estos fenómenos son función directa de la
a la del volcán Chichón en 1982. velocidad del viento y tienen un claro efecto de
Tres estudios recientes (Crowley & Kim, enfriamiento de las aguas superficiales. Los
1999; Free and Robock, 1999 y D'Arrigo et al., vientos a su vez están controlados por gradientes
1999) de la evolución climática en los últimos de presión atmosférica a nivel del mar (SLP, Sea
600 años, usando la variación solar y los índices Level Pressure), que son función de la cantidad
de las erupciones volcánicas mediante modelos de masa de aire en la columna atmosférica. Esta
han encontrado buena explicación de los registros transferencia de masa de aire entre los
directos e indirectos, excepto en los últimos continentes y los océanos está controlada por el
tiempos, cuando el efecto antrópico es evidente. calentamiento diferencial de mar y tierra, la
distribución geográfica de los continentes y el
ciclo anual de radiación, en particular, es
4.2 Hipótesis de Relación ENSO-
fundamental recordar que el hemisferio Norte
Volcanismo tiene mucha más masa emergida que el Sur. La
Figura 11 presenta esquemáticamente las
variaciones en presión atmosférica en los
Handler (1993) postula que pequeñas anticiclones oceánicos y los principales centros
modulaciones en el ciclo anual de presión de acción continentales, entre los veranos de
atmosférica a nivel del mar pueden producir el ambos hemisferios. Se nota que todos los
ENSO. A partir del análisis del monzón que anticiclones oceánicos incrementan la presión en
obedece a un calentamiento diferencial entre mar julio, en aproximadamente 6 hPa, incluso en el
y tierra plantea que se desarrolla un gradiente de hemisferio Norte, cuando es más caliente allí.
presión (calentamiento asociado a baja presión) el Esto ocurre básicamente por el calentamiento
cual produce vientos que traen humedad, ascenso continental y la circulación monzónica. La
de aire, convección y lluvia convectiva, todo esto transferencia de masa desde el sur del continente
asociado a transferencia de masa atmosférica. euro-asiático hacía los océanos es
Como la lluvia monzónica también disminuye particularmente significativa. Los anticiclones
durante los eventos ENSO, sería interesante oceánicos pueden acomodar más masa al
explicar ambas circunstancias. desplazarse hacía los polos sobre aguas más frías.
El ciclo anual de temperatura superficial del Realmente todo el Este del Pacífico subtropical
mar en el Pacífico ecuatorial del Este presenta un tiene presión atmosférica mayor durante este
máximo en febrero y un mínimo en octubre, con período. Los vientos responden a estos cambios
una amplitud de entre 3 y 4 °C. Durante los consistentemente. El índice de oscilación del Sur
eventos ENSO, fase cálida, la amplitud se es una medida del gradiente de presión a lo largo
incrementa hasta 5 o 7 °C, por incremento del del ecuador.
máximo, siguiendo un debilitamiento de los La presencia de aerosol estratosférico en
alisios. El ciclo anual de temperatura en la zona latitudes bajas ocasionada por una erupción
ecuatorial no obedece directamente al ciclo volcánica atenúa el ciclo anual de la radiación y
semianual de radiación solar. Esto indica que el puede alterar la transferencia de masa entre los
clima en esta zona está dominado por regiones en continentes y los anticiclones oceánicos,
las que el ciclo anual es dominante. Un factor particularmente en Eurasia. Esto produce centros
determinante es la transferencia de masas de aire anticiclónicos sobre el Pacífico con presión
desde los continentes hacia los anticiclones atmosférica más baja de lo normal, desplazados
oceánicos en ambos hemisferios (Gordon, 1953). hacía el ecuador, por lo tanto vientos alisios más
 Dyna 136, 2002 57

débiles y temperaturas mayores en la superficie de los eventos El Niño. Sin embargo, el análisis
del mar. Esto puede explicar que las anomalías en estadístico de otros autores no soporta esta
enero sean más intensas que en julio. hipótesis. En general, la comunidad
El aerosol inyectado a altas latitudes meteorológica y oceanográfica tiene dificultad
normalmente no alcanza a los trópicos o los con estas ideas, probablemente por sesgos
subtrópicos, por lo tanto estas zonas reciben la profesionales. Por ejemplo, los modelos que
cantidad normal de radiación, mientras que los tienen en cuenta los aerosoles y las mediciones
extratrópicos estarán más fríos de lo normal. En apenas empiezan a aparecer.
este caso se tienen gradientes más fuertes de
temperatura del ecuador al polo correspondiente y
se tendrá una circulación de Hadley intensificada
y no hay lugar a eventos El Niño.

Figura 11. Redistribución de la masa atmosférica Figura 12. Serie de tiempo de las temperaturas
de la tierra a los anticiclones oceánicos a medida del mar en el Pacífico tropical del Este desde
que el sol se “mueve” del hemisferio Sur en enero 1950 y ocurrencia de las erupciones precursoras
(cuadro superior) hasta el hemisferio Norte en de los evento El Niño (Handler, 1990).
julio. La presión atmosférica en HPa en cada uno
de los principales anticiclones y centros de acción
continental está indicada, así como los vientos 5. CONCLUSIONES
alisios. (Handler, 1990).
El patrón de variabilidad climática interanual
Los efectos en las lluvias monzónicas y la conocido como ENSO, por sus siglas en Inglés
alta correlación entre éstas y ENSO estarían representa un caso paradigmático para la
explicadas mediante este mecanismo propuesto interrelación entre diferentes componentes de las
por Handler. Las observaciones soportan ciencias geofísicas.
débilmente esta hipótesis. La Figura 12 muestra La atmósfera rota más rápido que la tierra
una serie de tiempo de las temperaturas del mar sólida durante el invierno del hemisferio Norte y
en el Pacífico tropical del Este desde 1950 y viceversa en el verano. La corriente de chorro
señala la ocurrencia de las erupciones precursoras subtropical que se intensifica en invierno es la
 Mesa
58

principal responsable de este ciclo. En ciclo correspondiente en el momento angular

consecuencia por conservación de momento atmosférico. Para todas las magnitudes es mayor
angular, el día es más largo durante el invierno el número de sismos en diciembre-enero-febrero
del hemisferio Norte (del orden de 1,5ms). De y menor en septiembre-octubre-noviembre,
manera semejante, durante un evento ENSO siendo intermedia para los meses restantes. A
cálido, ocurre un aumento de la duración del día, escala interanual, el estudio de estos registros en
por el debilitamiento generalizado de los vientos comparación con índices del ENSO, muestra que
ecuatoriales del Este y el fortalecimiento de las aunque pueden haber algunas épocas de buena
corrientes chorro subtropicales del Oeste. El coincidencia, no se puede apreciar unívocamente
cambio observado corresponde a un incremento una relación directa.
de 0,50 ms en la duración del día. Hoy se dispone Para mirar el tema de las relaciones clima y
de excelentes mecanismos para observar tanto la sismicidad con algún rigor habría que desarrollar
duración del día como el momento angular un modelo elastoplástico de la corteza terrestre y
atmosférico. Los cambios en el momento angular someter las placas a los esfuerzos que se derivan
atmosférico, incluyendo la troposfera y la de los modelos atmosféricos y oceánicos. No es
estratosfera hasta 1 hPa, explican hasta el 91,6% conveniente dar lugar a especulaciones poco
de la varianza interanual en la serie de duración científicas, basadas en temores y en la idea
del día, dejando un residuo de 55s. La errónea de que todo tiene que ver con todo, sin
estratosfera contribuye con aproximadamente el cuantificar y sin demostrar la importancia.
20% de la varianza. Los aerosoles producto de las erupciones
Generalmente se ha asumido que los cambios volcánicas pueden tener efectos importantes en el
en el nivel del mar asociados al ENSO son clima terrestre. La latitud de una erupción
estéricos, es decir por expansión térmica. Si este determina substancialmente su efecto climático.
fuera el caso, no habría lugar a ningún efecto Varios autores presentan evidencia en el sentido
rotacional. En realidad, estos cambios se de que las erupciones volcánicas tienden a
producen por redistribución de aguas cálidas a reforzar el patrón de la oscilación del Atlántico
gran escala, que pueden modelarse como ondas Norte y algunos también presentan algo de
baroclinas en la parte superior del océano, y que evidencia en el sentido de que hay relación con
pueden excitar variaciones en la rotación de la ENSO. A pesar de que estos autores presentan
tierra. Usando un modelo de dos capas se obtiene explicaciones teóricas para soportar estas
un incremento en la carga de presión sobre el sugerencias, todavía siguen siendo
fondo, por incremento de la masa ecuatorial de 3 especulaciones que necesitan más observación y
x 1014 kg, que equivale a un incremento en la análisis.
duración del día de  , si esta agua viene del
Norte y Sur de los 45°. Aunque estas
estimaciones son sólo del orden del 30% del REFERENCIAS
residuo no explicado, son del signo correcto y
ocurren en el momento correcto.
Sin perder de vista la tectónica de placas, es Barnes, R. et al., Atmospheric Angular
posible examinar la posibilidad de que existan Momentum fluctuations, length of day and
relaciones entre la variabilidad de la sismicidad polar motion. Proc. Royal Soc. London Ser. A.
en el tiempo, con la variabilidad climática, 387, 31-73, 1983
especialmente a escala interanual. Pequeñas CBluth, G., et al., Global tracking of the SO2
perturbaciones pueden amplificarse y producir clouds from the June 1991 Mount Pinatubo
grandes cambios. La posible conexión se apoya eruptions, Geophys. Res. Lett., 19, 151-154,
en los mecanismos de transferencia de momento 1992
angular entre la atmósfera, los océanos y la tierra Bostrom, R., Tectonic Consequences of the
sólida. El estudio de registros de temblores en la Earth's Rotation, Oxford Univ Press, 2000.
zona tropical del occidente de Sudamérica y Chinnery, M., comunicación personal, 1991
Centroamérica muestra claramente la existencia Chinnery, M. and Wells, F., On the correlation
de un ciclo anual en la sismicidad, siguiendo el between earthquake occurrence and
 Dyna 136, 2002 59

disturbances in the path of the rotation pole. Gordon, A., Seasonal changes in the mean
Rotation of the Earth. Edited by Melchior and pressure distribution over the world and some
Yumi. IAU, 1972 inferences about the general circulation. Bull.
Crowley, T. and Kim, K.,Modeling the Am. Meteorrol. Soc. 34, 357-367, 1953.
temperature response to forced climate change Graf, H., et al., On the interrelationship between
over the last six centuries. Geophys. Res. Lett. recent climatic trends, ozone change, and
26, 1901-1904, 1999. increasing greenhouse forcing. In:
Crutzen, P., The possible importance of CSO for Atmospheric Ozone as a Climate Gas. Edited
the sulfate layer of the stratosphere. Geophys. by Wang and Isaksen, Springer-Verlag, 1995.
Res. Lett. 3, 73-76., 1976 Groisman, P. Y. Regional climatic consequences
D'Arrigo, R., et al., Northern hemisphere of volcanic eruptions. Meteorol. Hydrol. 4, 39-
tempera-ture variability for the past three 45, 1985.
centuries: tree ring and model estimates, Clim. Habermann, R., Seismicity rate variations and
Change, 42, 663-675, 1999. systematic changes in magnitudes in
DeFoor, T., Robinson, E. and Ryan, S., Early teleseismic catalogs. Tectonophysics. 193, 277-
lidar observations of the June 1991 Pinatubo 289, 1991.
eruption plume at Mauna Loa Observatory, Haigh, J., The impact of solar variability on
Hawai. Geophysical Res. Lett. 19,187-190, climate. Science. 272. 981-984., 1996.
1992. Handler, P., The Response of Global and
Diaz, H. and Markgraf, V., El Niño, Historical Regional Climates to Volcanic Aerosols.
and Paleoclimatic Aspects of the Southern Submitted to Reviews of Geophysics, 1993.
Oscillation, Cambridge University Press, 1192.. Hastenrath, S., Climate Dynamics of the Tropics,
Dickey, J., et al., Global Propagation of Kluwer Academic Publishers. Boston, 1991.
Interannual Fluctuations in Atmospheric Hide, R., and Dickey, J., Earth's variable
Angular Momentum. Nature,357, 484-488, Rotation. Science. 253, 629-637, 1991.
1992. Hunt, A., A physical interpretation of the
Dickey, J., S. L. Marcus, R. Hide, T. M. Eubanks, correlation between El Niño and global
and D. H. Boggs. Angular momentum warming. Proceedings of the 1999 NOAA
exchange among the solid Earth, atmosphere, workshop on climate diagnostics and
and oceans: A case study of the 1982-1983 prediction. NOAA, 2000.
event, J. Geoph. Res., 99, B12, 23921-23937, Jäger, H., The Pinatubo eruption cloud observed
1994. by lidar at Garmisch-Patenkirchen. Geophys.
Ellsaesser, H., Isolating the climatological effects Res. Lett. 191-194, 1992.
of volcanoes, Rep. UCRL-89161, Lawrence Kant, Examination of the question whether the
Livermore Natl. Lab., 1983. Earth has undergone an alteration of its axial
England, [Link] of the Continental Crust. rotation, 1754.
In: Brown, G. et al.. Editors: Understanding the Kodera, K., Chiba, M., and Shibata, K., A
Earth. Cambridge University Press, 1992. general circulation model study of the solar and
Eubanks, T.M., et al., Causes of Rapidmotions of QBO modulation of the stratospheric
the Earth's pole. Nature, 334, 115-119, 1988. circulation during the Northern Hemisphere
Eubanks, T., Interactions between the winter. Geophys. Res. Lett. 18, 1209-1212,
atmosphere, oceans, and crust: possible oceanic 1991.
signals in earth rotation. Adv. Space Res.13, Kondratyev, K., Volcanoes and Climate, WCP-
(11),291-300, 1993. 54, WMO/TD-166, 1998.
Free, M. and Robock, A. Global warming in the Kondratyev, K. and Galindo, I., Volcanic
context of Little Ice Age. J. Geophys. Res. Activity and [Link], 1997
104,19057-19070,1999. Labitzke, K. and van Loon, H., Association
Glantz, M., et al., Teleconnections linking between the 11-year solar cycle, the QBO and
worldwide climate anomalies. Scientific basis the atmosphere. J. Atmosph. Terr. Phys., 50,
and societal impacts. Cambridge University 197-207, 1988.
Press, 1991.
 Mesa
60

Lamb, H., Volcanic dust in the atmosphere, with Peltier, W., Postglacial Variations inthe Level of
a cronology and assesment of its the Sea: Implications for Climate Dynamics
meteorological significance, Phylos. Trans. and Solid Earth Geophysics. Rev. of Geophy.,
Royal Soc. London, Ser. A. 266. 425-533, 1970. 36,(4),603-689, 1998.
Lamb, P. and Peppler R. A.,West Africa. In: Poveda, G. and Mesa, O.,"Feedbacks between
Teleconnections Linking Worldwide Climate hydrological processes in tropical South
Anomalies. Ed. M.H. Glantz, R.W. Katz y N. America and large scale oceanic-atmospheric
Nicholls. Cambridge University Press, 1991. phenomena", J. Climate. 10, 2690-2702, 1997.
Lambeck, K., Geophysical geodesy, the slow Preisig, J. R. Polar motion, atmospheric angular
deformation of the earth. Clarendon Oxford, momentum excitation and earthquakes-
1988. correlation and significance. Geophys. J.,
Lambeck, K., The Earth's variable rotation, Intern. 108, 161-178, 1992.
Cambridge University Press, New York, 1980. Ramanathan, V., The greenhouse theory of
MacDonald, Tidal Friction. Rev. Geophys. 2,467- climate change:A test by an inadvertant global
541, 1964. experiment. Science. 240, 293-299,1988.
Mansinha, L., and Smyle, D., Effect of Ramaswamy, V., Schuarzkopf, M. D., and
earthquakes on the Chandler Wobble and the Randel W. [Link] of ozone depletion in
secular polar shift. Jour. Geophys. Res. 72, the spatial and temporal pattern of recent
(18), 4731-4743,1967. lower-stratospheric cooling, Nature, 382, 616-
McCormick, P., and Veiga, R., SAGE II 618, 1995
measurements of early Pinatubo aerosols. Read, W., Frodeveax, L. and Waters, W.,
Geophys. Res. Lett. 19, 155-158, 1992. Microwave limb sounder measurements of SO
Mesa, O., Poveda, G. y Carvajal, L.,Introducción from Mt Pinatubo volcano, Submitted to
al Clima de Colombia, Universidad Nacional de Geophys. Rres. Lett., 1983.
Colombia, 1997. Rind, D., Balachandran, N. K., and Suozzo, R.,
Munk, W., Once again: Tidal friction. Q. J. R. Climate change and the middle atmosphere,
Astron. Soc. 9,352-375, 1968. Part II, The impact of volcanic aerosols. J.
Munk, W. and Mac Donald, G., The Rotation of Climate. 5, 189-208, 1992.
the Earth, Cambridge University Press. Robock, A., Volcanoes and Climate, in Climate
Newhall, C. and Self, S. The volcanic explosivity and Geosciences: A challenge for science and
index, VEI: an estimate of explosive magnitude society in the 21st Century. Nato Asi Ser. C.
for historical volcanism. J. Geophys. Res. 87. 285, 309-314, 1989.
1231-1238, 1982. Robock, A.,The volcanic contribution to climate
Oort, A., Angular Momentum Cycle in the change of the past 100 years. In: Greehouse
Atmosphere-Ocean-Solid Earth System. Gas Induced Climatic Change: A critical
Bulletin American Meteorological Society. 70, appraisal of simulations and observations,
(10), 1231-1242, 1989. edited by Schelesinger. Elsevier 429-444,1991.
Palais, J. and Sigurdson, H., Petrological Robock, A., Stratospheric control of climate.
evidence of volatile emissions from major Science. 272, 972-972, 1996.
historic and prehistoric volcanic eruptions. In: Robock, A. and Mao, J., Winter warming for
Berger, Dickinson and Kidson, editors of large volcanic eruptions. Geophys. Res. Lett.
Understanding Climate Change. Geophysical 19, 2405-2408, 1992.
Monographs 52, 1989. Robock, A.,and Free, M., Ice cores as an index of
Patterno, L., et al., Angular momentum transport global volcanism from 1850 to the present. J.
by Reynolds stresses determined from the Geophys. Res. 100, 11549-11567, 1995.
analysis of 100-year sunspot motions and its Robock, A., Volcanic eruptions and climate.
variations with solar cycle. Astronomy and Reviews of Geophysics. 38, 191-219, 2000.
Astrophysics. 252, 337-342, 1991. Shimazaki, K., and Takeuchi, H. Excitation of the
Peixoto, J. and Oort, A., Physics of Climate. Chandler wobble by large earthquakes. In:
American Institute of Physics, 1992. Rotation of the Earth. Edited by Melchior and
Yumi, 1972.
 Dyna 136, 2002 61

Shindell, D., Rind, D., Balachandran, N., Lean, J. Moon's Orbit. Rev. of Geophy., 38,(1), 37-59,
and Lonergan, P., Solar cycle variability, ozone 2000.
and climate. Science. 284, 305-308, 1999. Wyrtki, K., 1985. Water displacements in the
Smithsonian, 1981. Bulletin of the global Pacific and the genesis of El Niño cycles. J.
volcanism network, 7, (12), 1981 Geophysical Res., 90, 7129-32, 1985.
Solomon, S., Stratospheric ozone depletion: A Yu, L. and Rienecker, M., 1998. Evidence of an
review of concepts and history. Rev. Geophys. extratropical atmospheric influence during the
37, 275-316, 1999. onset of the 1997-98 El Niño. Geo. Res. Lett.,
Solomon, S., Portman, R. W. Garcia, R. R., 25, 3537, 1998.
Thomason, L. W., Poole, L. R. and
McCormick, M., The role of aerosol variations
in antropogenic ozone depletion at northern
midlati-tudes. J. Geophys. Res. 101, 6713-6727,
1996.
Stacey, F. D., Physics of the Earth. John Wiley &
Sons, New York, 1977.
Stowe, L., Carey, R. and Pellegrino, P.,
Monitoring the Mount Pinatubo aerosol layer
with NOAA 11 AVHRR data. Geophys. Res.
Lett. 19, 159-162, 1992.
Tie, X.., Brasseur, G., Granier, C., De Reudder,
A. and Larsen, N., Model study of polar
stratospheric clouds and their effect on
stratospheric ozone, 2, Model results. J.
Geophys. Res. 101, 12575-12584, 1996.
Toon, O., Volcanoes and climate. In:
Atmospheric effects and potential climatic
impact of the 1980 eruptions of Mount St.
Helens. Edited by Deepak, Nasa Conf. Publ.
2240, 15-36, 1982.
Toon, O. and Pollack, J., Atmospheric aerosol
and climate. Am. Sci. 68, 268-278, 1980
Van Loon, H. and Labitzke, K., Association
between the 11 years solar cycle, QBO and the
atmosphere, part IV, The stratosphere not
grouped by the phase of the QBO. J. Climate.
3, 827-837, 1990.
Veiga, R., SAGE II measurements of aerosols. In:
Optical Remote Sensing of the Atmosphere
Technical Digest. Optical Soc. Am., 5, 399-402,
1993.
Walker, D., 1988. Seismicity of the East Pacific
rise: Correlations with the Southern Oscillation
Index? Eos, Trans., 69, 857, 1988.
Walker, D., More evidence indicates link between
El Niño and seismicity. Eos, Trans., 76,
33,1995.
Walker, D., Seismic predictors of El Niño
revisited, Eos, Trans., 80, 281, 1999.
Williams, G., Geological Constraints on the
Precambrian History of Earth's Rotation and