SISMICIDAD E IMPLICACIONES ESTRUCTURALES EN EL ÁREA DE ESMERALDAS

Sismicidad e implicaciones estructurales en el área de Esmeraldas (Norte de

Ecuador): a partir de los experimentos Sublime y Esmeraldas

Kevin Manchuel1, Bernard Pontoise1, Nicole Béthoux1, Marc Régnier2, Yvonne Font1,
Vallentí Sallares3, Jordi Díaz4, Patricia Arreaga-Vargas5, Tony Monfret2, Hugo Yépes6
1
GéoAzur, Université de Nice Sophia-Antipolis, IRD, Université Pierre et Marie Curie, CNRS, Observatoire de
la Côte d’Azur, Villefranche sur mer, France

2
UGéoAzur, Université de Nice Sophia-Antipolis, IRD, Université Pierre et Marie Curie, CNRS, Observatoire
de la Côte d’Azur, Sophia-Antipolis, Valbonne, France

3
Unidad de Tecnología Marina,
Consejo Superior de Investigaciones Científicas, Passeig Maritim de la Barceloneta, 37, 08003 Barcelona,
España

4
Dpt. Geofísica y Tectónica
Institut de Ciències de la Terra ‘Jaume Almera’ IJA-CSIC, Lluis Solé i Sabaris s/n, 08028 Barcelona, España

5
Instituto Oceanográfico de la Armada (INOCAR)
Guayaquil, Ecuador

6
Instituto Geofísico de la Escuela Politécnica Nacional (IG-EPN)
Quito, Ecuador

Resumen.

La localización exacta de los terremotos costa afuera y en tierra es un desafío a largo plazo para la comunidad
científica, y este problema es de particular interés en las regiones de subducción como el margen activo de
Ecuador-Colombia, donde grandes terremotos y tsunamis han ocurrido. En el norte de Ecuador, la sismicidad
no está bien precisada debido a la ausencia de registros sísmicos entre la costa y los Andes. Por lo tanto, algunos
problemas importantes, como la geometría del slab, no están resueltos aún. Nosotros proponemos aquí estudiar
la sismicidad de la región de Esmeraldas (norte de Ecuador), usando datos de dos experimentos sismológicos
locales, SUBLIME y ESMERALDAS, llevados a cabo en esta parte de Ecuador. Estos fueron conducidos
respectivamente en 1998 y 2005, y durante ambos experimentos, estaciones en tierra y OBS (Sismómetros de
Fondo Oceánico) fueron desplegados simultáneamente. Las redes de SUBLIME y ESMERALDAS cubrieron
juntas una gran área, extendiéndose desde el margen hasta la Cordillera. Hemos localizado eventos registrados
usando diferentes modelos de velocidad de acuerdo a las áreas geológicas y topográficas. Los modelos de velocidad
usados fueron derivados a partir de perfiles sísmicos de gran ángulo a través del margen norte ecuatoriano, y de
un modelo de velocidad terrestre global. Nuestros resultados muestran que el Límite Superior pendiente arriba
(UdL) de la zona sismogénica está localizado a una profundidad superficial de ~10 Km, 40 Km al Este de la Fosa.
Utlizando resultados anteriores de sísmica activa y nuestras localizaciones, determinamos una inclinación del slab
de 10° al nivel de la fosa que aumenta a 25°-30° en la parte profunda hasta 120km de profundidad por debajo
de la Cordillera Occidental. Hemos también notado la presencia de una intensa y densa deformación en la parte
superior, desde los 40Km al Este de la Fosa hasta la Cordillera. El conocimiento de los patrones de sismicidad en
las zonas de subducción, como la localización del UdL, es un parámetro clave para la determinación de los riesgos
sísmicos y de tsunamis.

167
 GEOLOGIA Y GEOFISICA MARINA Y TERRESTRE DEL ECUADOR, PSE-001-2009

Abstract.

Accurate location of offshore and onshore earthquakes is a long-time challenge in the scientific community and
this problem is of particular interest in subduction regions like the active Ecuador-Colombian margin, where great
earthquakes and tsunamis occurred. In the north of Ecuador, the seismicity is not well constrained because of a
lack of recorded seismicity between the coast and the Andes. Therefore, some major problems, like geometry of
the slab, are not resolved yet. We propose here to study the seismicity of the Esmeraldas region (north Ecuador),
using data of two local seismological experiments, SUBLIME and ESMERALDAS, carried out in this part of
Ecuador. They were conducted respectively in 1998 and in 2005 and during both experiments, land stations and
OBSs (Ocean bottom seismometers) were deployed simultaneously. Together SUBLIME and ESMERALDAS
networks cover a large area, stretching from the margin up to the Cordillera. We locate recorded events using
different velocity models according to geological and topographical areas. The velocity models used were derived
from wide-angle seismic profiles across the north Ecuadorian margin, and from a global earth velocity model. Our
results show that the Updip Limit (UdL) of the seismogenic zone is located at a shallow ~10 km depth, 40 km east
of the trench. Using active seismic tomography of previous study and our location results we show a 10° dipping
slab in the trench area, increasing to 25°-30° down to 120km depth beneath the western Cordillera. We also note
the presence of an intensive thick deformation in the upper, from 40km east of the trench until the cordillera.
Knowledge of seismicity patterns in subduction zones, like the UdL location, is a key parameter for seismic and
tsunami risk determination.

1. Marco geodinámico y sismológico b) El segmento que se extiende desde los

0.5°N a 2.5°S está frente de Carnegie y es sometido
a un levantamiento general. No se han registrado
El margen Nor-Andino está siendo deformado grandes terremotos durante el último siglo en este
por la subducción de la placa de Nazca (5-7 cm/año) a dominio y la sismicidad está organizada en enjambres
lo largo de una dirección N80° (DeMets et al., 1990). de terremotos.
La placa de Nazca transporta la Cordillera de Carnegie,
una cordillera flotante de 200 Km de ancho (Figura Usando estas observaciones, una probabilidad
1), la cual subduce bajo el margen central ecuatoriano de 60% a 100% de ocurrencia de un gran terremoto
implicando una deformación cortical importante. superficial en la zona inter-placa fue estimada para
El flanco norte de la Cordillera de Carnegie divide el margen de Ecuador-Colombia (Nishenko, 1991).
el margen de Ecuador-Colombia en dos dominios Nuestra área de estudio está localizada en la parte
sísmica y tectónicamente opuestos (Collot et al., 2002; Norte de Ecuador, en los alrededores de la ciudad de
Collot et al., 2004; Gutscher et al., 1999; Pontoise y Esmeraldas, la cual es un sitio con alta vulnerabilidad
Monfret, 2004): ocupada por 60 000 habitantes y con una refinería de
petróleo. Por lo que, el reconocimiento de fallas activas
a) El segmento Norte que se extiende por lo y la comprensión de la sismicidad y geodinámica de
menos a lo largo del área de ruptura del terremoto esta región son cruciales.
de 1906 (Norte de 0.5°N) es subsidente, y 4 mega
terremotos ocurrieron en la zona interplaca durante el Algunos experimentos de sísmica marina
siglo 20. La zona de ruptura de 500 Km de longitud llevados a cabo durante la última década a lo
(Kelleher, 1972) del evento de 1906 (Mw=8.8) largo del margen activo de Ecuador-Colombia,
(Kanamori y Given, 1981; Kanamori y McNally, 1982) proporcionan un conocimiento detallado de la
fue parcialmente reactivada por 3 eventos compresivos estructura del margen (Calahorrano, 2005; Collot
(Mendoza y Dewey, 1984) que ocurrieron en 1942 et al., 2002; Collot et al., 2004; Gailler et al.,
(Mw=7.8) (Swenson y Beck, 1996), 1958 (Mw=7.7) 2007; Graindorge et al., 2004; Marcaillou et al.,
y 1979 (Mw=8.2) (Beck y Ruff, 1984; Herd et al., 2006; Sage et al., 2006). En el país, el Instituto
1981; Kanamori y McNally, 1982). La mayoría de Geofísico de la Escuela Politécnica Nacional (IG-
las soluciones del Tensor de Momento-Centroide del EPN) mantiene una red sísmica regional (RENSIG)
catálogo de Harvard, durante el período 1976-2008 con la mayoría de los instrumentos localizados
tienen un mecanismo de tipo compresional. en el área volcánica. La ausencia de suficiente

168
 SISMICIDAD E IMPLICACIONES ESTRUCTURALES EN EL ÁREA DE ESMERALDAS

instrumentación en la región ante-arco afecta la a los autores proponer diferentes geometrías del
resolución de la determinación hipocentral en esta slab en el norte de Ecuador. Algunos de ellos
área. Hasta ahora, los catálogos globales proveen postulan un buzamiento del slab hacia el Este
una imagen difusa de la sismicidad costera, y solo con un ángulo de 25° a 40° (Guillier et al., 2001;
pocos estudios sismológicos fueron conducidos Pontoise y Monfret, 2004; Taboada et al., 2000) y
en la región de Esmeraldas (Guillier et al., 2001; otros sugieren un slab horizontal a una profundidad
Pontoise y Monfret, 2004). La pobreza de los datos de 100Km (Gutscher et al., 1999; Gutscher et al.,
de sismicidad entre la costa y la cordillera permite 2000).

Figura 1. Croquis geodinámico del margen activo ecuatoriano. La batimetría es tomada de Michaud et al. (2006). Las líneas punteadas
representan la superficie de ruptura de los cuatro grandes terremotos de subducción que ocurrieron en el siglo 20, y las estrellas amarillas
representan sus epicentros (Collot et al., 2002). El vector de movimiento de la Placa de Nazca es tomado de Trenkamp et al. [2002].
DGM= Megafalla Dolores-Guayaquil. El Bloque Nor-Andino está empezando a desplazarse hacia el noreste a lo largo de la DGM

169
 GEOLOGIA Y GEOFISICA MARINA Y TERRESTRE DEL ECUADOR, PSE-001-2009

Con el propósito de obtener una mejor implicadas en este proyecto son el Instituto Jaume
definición de las zonas sísmicas en la región de Almera (IJA) y la Unidad Tecnológica del Mar (UTM)
Esmeraldas, adicionamos los datos a partir de dos de Barcelona, en España, el Instituto Geofísico
redes temporales ubicadas entre el margen y los de la Escuela Politécnica Nacional (IG-EPN) de
Andes. En efecto, ambos experimentos tienen Quito, en Ecuador, el Instituto Oceanográfico de
desplegadas redes sismológicas terrestres y marinas, la Armada (INOCAR) de Guayaquil, en Ecuador
que permiten una buena cobertura azimutal. Nuestro y el Instituto Colombiano de Geología y Minería
propósito es evidenciar características regionales (INGEOMINAS) de Bogotá, en Colombia.
a lo largo del margen Norte ecuatoriano, debido a
que los 3 meses de datos que hemos usado no nos Este artículo presenta los resultados a partir
permiten distinguir características a pequeña escala. de los análisis de los registros de la sismicidad
La identificación de las principales fallas y sus natural, recopilados durante ambas fases de los
mecanismos serán utilizados para caracterizar el experimentos ESMERALDAS y SUBLIME.
estilo de la deformación a través de la determinación El marco instrumental fue (ver figura 2 para la
de tensores de esfuerzos regionales del margen activo localización de las redes):
y, después, caracterizar el marco de la geodinámica
global del área. a- Las estaciones marinas:
Durante el experimento SUBLIME,
15 OBS de corto período “UTIG-ORSTOM”
fueron desplegados. Ver Pontoise y Monfret
2. Experimentos (2004) para informaciones técnicas. Dos tipos de
instrumentos fueron desplegados durante el crucero
ESMERALDAS. 7 OBS de corto periodo “UTIG-
El experimento SUBLIME (Pontoise y ORSTOM” más 19 OBS HYPOCAMPE, una nueva
Monfret, 2004) fue conducido en 1998, usando generación de OBS recientemente desarrollados
el barco de investigación ANTEA del IRD (IRD: en el IRD (UMR Geosciences Azur) (Hello et al.,
Instituto de Investigación para el Desarrollo). 1992). La duración del registro fue incrementada
Durante un período de tres semanas, 15 OBS hasta 6 meses y el registro dinámico hasta los 24
(Sismómetros de Fondo Oceánico – ver abajo bits. 9 de estos OBS pueden ser equipados con
para descripción de OBS) y 10 estaciones en tierra sensores marinos de banda ancha Guralp GMG40T
fueron desplegados. Este experimento es parte de la 3C y 10 con geófonos Mark Products 4.5Hz.
colaboración entre Francia y Ecuador. Los trabajos
fueron realizados en estrecha cooperación entre b- Las estaciones en tierra:
Géosciences Azur (CNRS, IRD, UNSA y UPMC), Durante SUBLIME, 6 estaciones LEAS
Villefranche sur mer, Francia, el Instituto Geofisico “HATHOR” y 4 estaciones telemétricas fueron
de la Escuela Politécnica Nacional (IG-EPN) de desplegadas todas equipadas con sensores de banda
Quito, Ecuador y el Instituto Oceanográfico de la corta 3C. Durante ESMERALDAS, 31 estaciones
Armada (INOCAR), Guayaquil, Ecuador. fueron desplegadas juntas. 10 estaciones son
“HATHOR” y 21 son AGECODAGIS “TITAN”.
Durante el experimento ESMERALDAS, Sus dinámicas son de 24 bits. Los sensores fueron:
conducido desde el 10 de Marzo al 14 de Junio  banda ancha (BB): 2 CMG 3, 4 CMG40
del 2005, 26 OBS y 31 estaciones en tierra fueron (Guralp)
desplegados. Dos diferentes trayectos definen el  período intermedio: 8 Lennartz-10s , 14
experimento ESMERALDAS: durante el primero, Lennartz-5s
el B/I L’Atalante investiga el área de estudio costa  corto período: 3 Mark Products L4C-1s
afuera y la red registró los disparos de cañón de aire Cada estación fue equipada con un reloj
de gran potencia a bordo del barco (18 000 disparos) GPS que proveía el tiempo absoluto continuamente
a lo largo de diferentes líneas sísmicas cruzadas y la localización exacta de la estación.
(~25). En la segunda parte, la red fue mantenida
en el lugar hasta inicios de Junio para registrar la Para la mayoría de las estaciones buscamos
sismicidad local (sismología pasiva). Este proyecto espacios tranquilos, lejos de los caminos y en las
es un trabajo de colaboración entre Francia, Ecuador, afueras de las ciudades. El uso de sensores de
Colombia y España. Las instituciones científicas banda ancha implicó instalar cuidadosamente las

170
 SISMICIDAD E IMPLICACIONES ESTRUCTURALES EN EL ÁREA DE ESMERALDAS

estaciones en el campo. Por lo tanto cavamos un barril de plástico. Escogimos una tasa de muestreo
hueco, de alrededor de un metro de profundidad, y de 125 puntos por segundo. Los datos fueron
construimos un pilote de concreto para instalar el colectados una vez al mes cuando se visitaban las
sensor, el cual fue cerrado herméticamente en un estaciones.

Figura 2. Configuración de las redes sismológicas. La Batimetría está tomada de Michaud et al. (2006). Los triángulos rojos indican
las estaciones de ESMERALDAS. Los triángulos amarillos muestran las estaciones de SUBLIME. Los triángulos verdes indican las
estaciones del IG-EPN usadas en este estudio. Las líneas rojas representan las líneas de la sísmica activa.

3. Datos localizados. Ver Pontoise y Monfret (2004) para la

determinación de la calidad de los datos. Las casi
Se detectaron los eventos sísmicos a partir de continuas fuertes lluvias debidas al gran evento
registros continuos, siguiendo dos metodologías: El Niño de 1997–1998, causó altos niveles de
Durante SUBLIME hemos analizado los datos de ruidos en las estaciones terrestres. Adicionalmente,
cada uno de los OBS y extraímos eventos sísmicos inundaciones y deslizamientos hicieron difícil
identificados a partir de los registros de todas las realizar un mantenimiento apropiado a la red. A pesar
estaciones (terrestres y marinos). 700 eventos fueron de esto, las estaciones en tierra registraron alrededor
detectados y entre ellos, 117 fueron apropiadamente de 50 sismos.

171
 GEOLOGIA Y GEOFISICA MARINA Y TERRESTRE DEL ECUADOR, PSE-001-2009

Para los datos de ESMERALDAS, hemos En el caso de la sismicidad del margen

aplicado primero un detector de eventos STA/ ecuatoriano, el medio es fuertemente heterógeneo,
LTA para el registro continuo desde las estaciones variando desde la Fosa hasta los Andes. El uso de
terrestres, fijando la condición de detecciones diferentes modelos de velocidades como una función
coincidentes en cuatro estaciones, y construyendo de la localización de la estación es fuertemente
un catálogo de eventos. Luego usamos este catálogo recomendado.
para extraer los sismógrafos de cada estación
(terrestres y marinos). También verificamos que Hemos usado la herramienta SEISAN
este catálogo contenga todos los eventos sísmicos (Havskov y Ottemöller, 2000) para crear una base
que fueron detectados anteriormente por la red de datos, leer sismogramas y localizar eventos
RENSIG. (código Hypocenter incluído en SEISAN). También
leimos polaridades con el propósito de calcular los
Adicionalmente hemos analizado mecanismos focales con FOCMEC (Snoke et al.,
visualmente los datos continuos de las dos mejores 1984). Sin embargo, los mecanismos focales no son
estaciones en tierra, y de cuatro OBS para encontrar presentados aquí debido a que tenemos todavía que
otros pequeños eventos que no fueron identificados precisar estos usando tasas de amplitud.
por el detector. Obtuvimos registros para 519
eventos correspondientes a eventos locales o Luego hemos activado las localizaciones
regionales con el primer método y otros 946 con el usando el código Hypoellipse (Lahr, 1999). Esta
segundo proceso. Finalmente 1091 eventos fueron técnica de localización permite el uso de varios
estudiados y solo 544 pudieron ser localizados, modelos de velocidad 1D, asignados a las estaciones
de los cuales 364 mostraron un RMS mejor que 1 ubicadas en diferentes regiones estructurales. El
segundo. mismo código fue usado para procesar los datos del
experimento SUBLIME.
Por otra parte, hemos extraído 53 eventos
telesísmicos localizados en el catálogo mundial de 4.2 El modelo de velocidad
Harvard, tanto para las estaciones en tierra y OBS
BB. Estos registros son usados para los estudios Cualquiera que sea el código de localización
de anisotropía y el conocimiento estructural usado, los dos principales factores que permiten una
profundo (función de transferencia). En este buena localización son las coberturas azimutales de
artículo solo nos focalizamos en la localización la red disponible, respecto al área sismogénica, y el
sísmica local. modelo de velocidad.

4. Técnicas de Localización Concerniente al experimento SUBLIME,

Monfret y Pontoise (2004) discutieron la elección
4.1 Los algoritmos de localización de sus modelos. Una interpretación preliminar
de los datos de gran ángulo obtenidos durante el
Numerosos códigos son publicados para crucero SISTEUR en el 2000, en el área de Manta,
localizar eventos sísmicos. La mayoría de ellos 100 Km al Sur de la región de Esmeraldas, indica
usan un modelo de velocidad promedio 1D, como sin ambigüedad la presencia de una inversión
el Hypo71 (Lee y Lahr, 1975), y el Hypocenter de velocidad en la estructura de la velocidad del
(Lienert et al., 1986). Algunos códigos más margen ecuatoriano. Esta inversión de velocidad
sofisticados permiten el uso de varios modelos de es inferida por el análisis de una zona de sombra
velocidad 1D, como el Hypoellipse (Lahr, 1999). de alrededor de 1s en las secciones registradas,
Los pasos clásicos son de ´´pre-localizar´´ el lo cual es interpretado como el resultado de un
evento con un código simple, usando un modelo contraste de velocidad entre la placa superior y el
apropiado 1D. tope de la placa oceánica de Nazca subductada.
Modelos de velocidad de una dimensión son
Este paso preliminar permite chequear la deducidos a partir de estos datos de gran ángulo
coherencia de los tiempos de llegada, leídos en los y son usados para la localización de terremotos
diferentes sismógrafos. Luego las localizaciones en el área de Esmeraldas (ver Figura 3). Esto
son mejoradas, usando códigos más sofisticados y la mejora ampliamente las determinaciones de los
exploración del espacio modelo. parámetros hipocentrales.

172
 SISMICIDAD E IMPLICACIONES ESTRUCTURALES EN EL ÁREA DE ESMERALDAS

Figura 3. Modelos de Velocidad deducidos a partir de los datos sísmicos de gran ángulo en el área de Manta, usados en las localizaciones
de SUBLIME. Se extrapolaron estos modelos al área de Esmeraldas para mejorar las localizaciones. La imagen es de (Pontoise y Monfret,
2004).

Con respecto a ESMERALDAS, dos modelos de SUBLIME y ESMERALDAS y los cuales exhiben
de velocidad fueron usados. El modelo marino de un RMS mejor que 1 segundo. La Figura 4b muestra la
velocidad es derivado a partir de perfiles sísmicos proyección de los hipocentros sobre un plano vertical
de gran ángulo a través del margen de Esmeraldas y N115°, normal a la Fosa. La buena cobertura de la
aplicados a estaciones OBS (Agudelo, 2005; Gailler red de ESMERALDAS permitió obtener una buena
et al., 2007). El modelo tierra adentro, deducido de precisión de la distribución de la sismicidad en la placa
CRUST2.0 (un modelo de velocidad terrestre global sobreyaciente y en el slab, a lo largo del margen Norte
especificado desde el grado2*2 (Bassin et al., 2000)), ecuatoriano, desde la superficie hasta los 120 Km de
es para las estaciones en tierra, y pone la interface profundidad hacia abajo.
del Moho a una profundidad de 32 Km. Este valor
aparece como un buen promedio entre la profundiad Inmediatamente al Oeste de la fosa, una
superficial del Moho bajo las estaciones costeras y su actividad microsísmica es observada (Figuras 4
gran profundidad bajo la Cordillera (Tabla 1). Para los y 5). Pontoise y Monfret (2004) propusieron que
modelos marinos un valor Vp/Vs de 1.81 es deducido esta actividad sísmica refleja el curvamiento de la
del estudio de Pontoise y Monfret, 2004. Para los placa cuando entra a la zona de subducción. Las
modelos tierra adentro hemos chequeado valores Vp/ localizaciones epicentrales obtenidas durante el
Vs desde 1.74 hasta 1.81 y concluimos que el mejor experimento ESMERALDAS confirman la presencia
es 1.78. de un enjambre sísmico al Oeste de la Fosa. La
profundidad del enjambre alcanza 40 Km y puede
también estar relacionado con el curvamiento del
5. Resultados y discusión slab. Algunos terremotos profundos en este enjambre
están asociados a la calidad D y no son confiables.
La Figura 4a muestra la distribución epicentral Consecuentemente no consideramos estos datos en
de 481 eventos localizados obtenidos usando los datos nuestra interpretación.

173
 GEOLOGIA Y GEOFISICA MARINA Y TERRESTRE DEL ECUADOR, PSE-001-2009

Tabla 1a. Modelo de velocidad de ondas P para la inversión de hipocentros de los datos de OBS (Agudelo, 2005 y Gailler et al., 2007).
Modelo de velocidad de las ondas S obtenido con una relación Vp/Vs de 1.81

Tabla 1b. Modelo de velocidad de ondas P para la inversión de hipocentros de los datos de las estaciones en tierra (Bassin et al., 2000).
Modelo de velocidad de las ondas S obtenido con una relación Vp/Vs de 1.78.

Desde la Fosa hasta los ~40Km hacia tierra, históricos, asociados con los terremotos de 1906 y
ninguna sismicidad es detectada. A esta distancia 1979 están en buen acuerdo con nuestra proposición.
de la Fosa, la sismicidad observada está localizada Adicionalmente a la generación co-sísmica, el tsunami
a una profundidad de ~10Km (Figuras 4b y 5). Esta puede ser amplificado por deslizamientos submarinos
observación muestra que, al este de la Fosa, hay muy activados por terremotos.
poca actividad sísmica dentro de los pocos kilómetros
más superficiales de la zona interplaca. Debido a la escasez de datos tierra adentro,
algunos autores no observaron una sismicidad cortical
Hemos interpretado esta distribución de en la placa superior en la región de Esmeraldas (Guillier
terremotos superficiales como una indicación de la et al., 2001; Pontoise y Monfret, 2004). Además, no
profundidad del Límite superior pendiente arriba(UdL) se han determinado profundidades confiables para
de la zona sismogénica (Tichelaar y Ruff, 1991). terremotos corticales dentro de la placa superior. Por
El UdL en la región de Esmeraldas está localizado lo tanto estos autores postulan que el bloque costero
a una profundidad superficial (~10 Km) comparada de Ecuador, compuesto de varios bloques oceánicos
con aquella en el Norte de Chile (~20 Km) (Hussen acrecionados (Cediel et al., 2003; Kerr et al., 2002),
et al, 1999). Como el UdL es un parámetro clave para actúa como un cuerpo no deformable. Nuestros datos
la generación de tsunamis, la profundidad superficial muestran una compleja distribución de la sismicidad
del UdL incrementa el riesgo sísmico y de tsunami en la placa superior alternando vacíos y pequeños
en el norte de la región ecuatoriana. Dos tsunamis grupos.

174
 SISMICIDAD E IMPLICACIONES ESTRUCTURALES EN EL ÁREA DE ESMERALDAS

Figura 4A. Mapa de epicentros. Los círculos amarillos son localizaciones obtenidas usando los datos de SUBLIME, y los círculos rojo,
verde, azul y negro representan la calidad de las localizaciones de ESMERALDAS. Las áreas grises muestran el ancho de la sección
transversal b. La línea roja es la traza en superficie del perfil sísmico de gran ángulo de Gailler et al. (2007).

175
 GEOLOGIA Y GEOFISICA MARINA Y TERRESTRE DEL ECUADOR, PSE-001-2009

Figura 4B. Corte de la sismicidad a lo largo de la seccion presentada en la Figura 4A. Ver leyenda en la Figura 4A.

Tabla 2. Calidad de las localizaciones y errores asociados a las localizaciones de HYPOELLIPSE.

Estas características están de acuerdo con deformación frágil. Este comportamiento mecánico
la hipótesis de un margen heterógeneo acrecionado. indica un estilo de tectónica de cobertura densa, como
Con nuestros 3 meses de datos uno de estos vacíos, se observa en el promontorio andino en Argentina
a 120 Km al este de la Fosa y extendiéndose desde (Regnier et al., 1992). Los dominios frágiles densos
la superficie hacia abajo hasta el slab parece ser en la placa superior son frecuentes en el contexto
bien resuelto. Este gap separa dos dominios de tectónico implicando colisión, fuerte acoplamiento
intensa actividad sísmica. Entre la Fosa y el vacío a lo largo de la interfase de placa y/o segmento de
la distribución de la sismicidad es superficial e subducción con slab horizontal. Esto usualmente
implica que la placa superior que está cerca de la implica un régimen de esfuerzos compresionales en la
Fosa está siendo deformada. Al Este del vacío, una corteza de la placa superior, como es probablemente
intensa sismicidad cortical dentro del bloque costero el caso en el bloque costero del Norte de Ecuador,
y la pendiente Oeste de los Andes, es organizada en hasta la pendiente Oeste de los Andes. La frecuencia
modelos con buzamiento tanto hacia el Este como al de los grandes terremotos históricos expresa también
Oeste alcanzando una profundidad de 40 Km (Figuras la presencia de un fuerte acoplamiento a lo largo del
4b y 5). Esta distribución de la sismicidad indica la límite de la placa mientras que la tasa de convergencia
extensión mínima en profundidad del dominio de sea moderada (5-7cm/año).

176
 SISMICIDAD E IMPLICACIONES ESTRUCTURALES EN EL ÁREA DE ESMERALDAS

Figura 5. Esquema de Interpretación de la sismicidad sin las localizaciones de calidad D de Esmeraldas. Todos los epicentros tienen el
mismo color para una mejor comprensión. La imagen tomográfica es de Gaillet et al. [2007] y nos permite restringir la interfase de la placa.
Ver figura 4 (A). para localización del perfil. La línea roja representa el tope del slab. Basados en esta geometría, el slab subyace una placa
superior que está significativamente deformada.

Pontoise y Monfret (2004) usando las del slab sugiere que la mayoría de los eventos del
localizaciones de SUBLIME y el catálogo ecuatoriano enjambre sísmico ocurren dentro de la placa superior,
propusieron un modelo de la geometría de la zona de la cual está bajo deformación. La hipótesis de un slab
Benioff con dos ángulos de inclinación. Conectando horizontal de 100 Km de profundidad (Gutscher et al.,
los pequeños grupos sísmicos del Este y Oeste de la 1999; Gutscher et al., 2000) aparece como no probado
Fosa ellos encontraron un primer ángulo de inclinación por la presencia de terremotos más profundos, hacia
de 10° a través de la Fosa. Hacia el Este, propusieron abajo hasta 120 Km (Figuras 4b y 5), localizados en la
un valor de ~35° para el ángulo de inclinación. A continuidad del slab.
partir de las determinaciones del experimento de
ESMERALDAS, entre los 300 y 400 Km desde la El análisis de los datos de ambos experimentos
Fosa, la inclinación del slab es de 25°-30º hacia el Este. sismológicos, SUBLIME y ESMERALDAS permitió:
Este ángulo es significativamente más bajo que los 35°
propuestos por varios autores (Guillier et al., 2001; 1) Determinar un cambio en la inclinación
Taboada et al., 2000). Cerca de la Fosa, la distribución de la placa Nazca. Al nivel de la fosa ésta
de la actividad sísmica es bastante volumétrica y no muestra una inclinación de 10º que permanece
permite delinear la zona interplaca de manera clara. constante hasta 100 km hacia el continente.
Para mejorar la localización de la interfase de la placa A partir de este punto, la inclinación aumenta
sismogénica en esta área, los hipocentros son ploteados hasta alcanzar ≈ 25º-30º hasta 120 km de
en el tope de la imagen tomográfica de Gailler et al. profundidad.
(2007) (figura 5). Gailler et al. (2007) interpreta la 2) Proponer que existe una intensa deformación
zona azul inclinada de baja velocidad como el canal en la placa superior.
de subducción. Esta capa de baja velocidad tiene una 3) Definir el límite superior (UdL) de la zona
inclinación de 10º, la cual es consistente con la idea sismogénica, pendiente arriba, a lo largo del
propuesta por Pontoise y Monfret (2004). Conectando margen norte ecuatoriano a una profundidad
esta zona de baja velocidad con el segmento profundo superficial de ~10 Km. Esta profundidad
de la WBZ definida previamente, podemos observar superficial para el UdL nos lleva a estimar un
que el ángulo de inclinación del slab aumenta de 10º alto riesgo sísmico y de tsunami para el área
a ≈ 25º - 30º en la parte más profunda. Esta posición de Esmeraldas.

177
 GEOLOGIA Y GEOFISICA MARINA Y TERRESTRE DEL ECUADOR, PSE-001-2009

Referencias Long range seismic refraction using digital OBS, Ann.

Geophys., C51, 1992.
Agudelo, W., Imagerie sismique quantitative de la marge Herd, D. G., Youd, T. L., Meyer, H., Arango C., Person, W. J., y
convergente d’Equateur-Colombie: Application des méthodes Mendoza, C., The Great Tumaco, Colombia Earthquake of 12
tomographiques aux données de sismique réflexion multitrace December 1979, Science, 211, 441-445, 1981.
et réfraction-réflexion grand-angle des campagnes SISTEUR Hussen, S., Kissling, E., Flueh, E. R., y Asch, G., Accurate
et SALIERI, Thèse de doctorat thesis, Villefranche sur mer, hypocenter determination in the seismogenic zone of the
Pierre et Marie Curie, 2005. subducting Nazca Plate in northern Chile using a combined
Bassin, C., Laske, G., y Masters, G., 2000, The Current Limits of on-/offshore network, Geophysical Journal International, 138,
Resolution for Surface Wave Tomography in North America: 687-701, 1999.
Eos Transactions American Geophysical Union, v. 81, p. 897. Kanamori, H., y Given, J. W., Use of long-epriod surface waves for
Beck, S. L., y Ruff, L. J., The rupture process of the great 1979 rapid determination of earthquake-source parameters, Physics
Colombia earthquake: evidence for the asperity model, Journal of the Earth and Planetary Interior, 27, 8-31, 1981.
of Geophysical Research, 89, 9281-9291, 1984. Kanamori, H., y McNally, K. C., Variable rupture mode of the
Calahorrano, A., Structure de la marge du Golfe de Guayaquil subduction zone along the Ecuador-Colombia coast, Bulletin
(Equateur) et propriétés physiques du chenal de subduction, of the Seismological Society of America, 72, 1241-1253, 1982.
à partir de données de sismique marine réflexion et réfraction, Kelleher, J., Rupture zones of large South American earthquakes
PhD Thesis, Université P. et M. Curie, Paris VI, 221 pp., 2005. and some predictions, Journal of Geophysical Research, 77,
Cediel, F., Shaw, R. P., y Càceres, C., Tectonic assembly of the 2087-2103, 1972.
Northern Andean Block: The Circum-Gulf of Mexico and the Kerr, A. C., Aspden, J. A., Tarney, J., y Pilatasig, L. F., The nature
Caribbean: Hydrocarbon habitats, basin formation, and plate and provenance of accreted oceanic terranes in western
tectonics, AAPG, 815-848, 2003. Ecuador: geochemical and tectonic constraints, Journal of the
Collot, J.-Y., Charvis, P., Gutscher, M. A., y Operto, S., Geological Society, 159, 577-594, 2002.
Exploring the Ecuador-Colombia active margin and interplate Lahr, J. C., Hypoellipse: A computer program for determining
seismogenic zone, EOS Transactions, American Geophysical local earthquake hypocentral parameter, magnitude and first
Union, 83, 189-190, 2002. motion pattern, U.S. Geological Sur. Open File Rep., 99(23),
Collot, J.-Y., Marcaillou, B., Sage, F., Michaud, F., Agudelo, W., 116 pp., 1999.
Charvis, P., Graindorge, D., Gutscher, M.-A., y Spence, G., Are Lee, W. H., y Lahr, J. C., Hypo71 (revised): a computer program
rupture zone limits of great subduction earthquakes controlled for determining hypocenters, magnitudes and first motion
by upper plate structures ? Evidence from multichannel seismic pattern of local earthquakes, U.S. Geological Sur. Open File
reflection data acquired across the northern Ecuador-southwest Rep., 75(311), 116pp., 1975.
Colombia margin, Journal of Geophysical Research, 109, Lienert, B. R. E., Berg, E., y Frazer, L. N., Hypocenter: An
B11103, doi:10.1029/2004JB003060, 2004. earthquake location method using centered, scaled and
DeMets, C., Gordon, R. G., Argus, D. F., y Stein, S., Current plate adaptivly least squared, BSSA, 76, 771-783, 1986.
motions, Geophysical Journal International, 101, 425-478, Marcaillou, B., Spence, G., Collot, J.-Y., y Wang, K., Thermal
1990. regime variation from bottom simulating reflector along
Gailler, A., Charvis, P., y Flueh, E. R., Segmentation of the Nazca the north Ecuador - south Colombia margin: Relation to
and South American plates along the Ecuador subduction margin segmentation and great subduction earthquakes,
zone from wide angle seismic profiles, Earth and Planetary Journal of Geophysical research, 111, B12407,
Sciences Letters, 260, 444-464, 2007. doi:10.1029/2005JB004239, 2006.
Gailler, A., Charvis, P., Sallares, V., Collot, J.Y., Graindorge., Mendoza, C., y Dewey, J. W., Seismicity associated with the great
D., y Calahorrano, Structure of the subduction channel at the Colombia-Ecuador earthquakes of 1942, 1958 and 1979:
Ecuador-Carnegie Ridge convergent margin from seismic implications for barrier models of earthquake rupture, Bulletin
modelling and inversion, Joint Assembly, Eur. Geophys. Soc., of the Seismological Society of America, 74, 577-593, 1984.
Nice, 2003. Michaud, F., Collot, J.-Y., Alvarado, A., López, E., y el personal
Graindorge, D., Calahorrano, A., Charvis, P., Collot, J.-Y., y científico y técnico del INOCAR, República del Ecuador,
Béthoux, N., Deep structures of the Ecuador convergent Batimetría y Relieve Continental, publicación IOA-CVM-02-
margin and the Carnegie Ridge, possible consequence on Post, INOCAR, Guayaquil, Ecuador, 2006.
great earthquakes recurrence interval, Geophysical Research Nishenko, S. P., 1991, Circum-Pacific seismic potential 1989-
Letters, 31, L04603, doi:10.1029/2003GL018803, 2004. 1999: Pure and Applied Geophysics, v. 135, p. 169-259.
Guillier, B., Chatelain, J.-L., Jaillard, E., Yepes, H., Poupinet, G., Pontoise, B., y Monfret, T., Shallow seismogenic zone detected
y Fels, J.-F., Seismological evidence on the geometry of the from an offshore-onshore temporary seismic network in the
orogenic system in central-northern Ecuador (South America), Esmeraldas area (northern Ecuador), Geochemistry, Geophysics,
Geophysical Research Letters, 28, 3749-3752, 2001. Geosystems, 5, Q02009, doi:10.1029/2003GC000561, 2004.
Gutscher, M. A., Malavieille, J., Lallemand, S., y Collot, J.-Y., Regnier, M., Chatelain, J.-L., Smalley, R. J., Chiu, J. M., Isacks,
Tectonic segmentation of the North Andean margin: impact B. L., y Araujo, M., Seismotectonics of the sierra pie de Palo,
of the Carnegie Ridge collision, Earth and Planetary Science a basement block uplift in the Andean foreland of Argentina,
Letters, 168, 255-270, 1999. Bulletin of the seismological Society of America, 82, 2549-
Gutscher, M. A., Spakman, W., Bijwaard, H., y Engdahl, E.R., 2571, 1992.
Geodynamics of flat subduction: Seismicity and tomographic Sage, F., Collot, J.-Y., y Ranero, C. R., Evidences for subduction-
constraints from the Andean margin, Tectonics, 19, 814-833, erosion and interplate heterogeneity at the interplate boundary
2000. along the Central Ecuador active margin, Geology, 34, 997-
Havskov, J., y Ottemöller, L., SEISAN earthquake analysis 1000, 2006.
software, Seismological Research Letters, 70, 532-534, 2000. Snoke, J. A., Munsey, J. W., Teague, A. G., y Bolliger, G. A., A
Hello, Y., Charvis, P., Pontoise, B., Nakamura, Y., y Chen, A.T., program for focal mechanism determination by combined use

178
 SISMICIDAD E IMPLICACIONES ESTRUCTURALES EN EL ÁREA DE ESMERALDAS

of polarity and SV-P amplitude ratio data, Earthquake notes, deformation (Colombia), Tectonics, 19, 787-813, 2000.
55, 15, 1984. Tichelaar, B. W., y Ruff, L., Seismic coupling along the Chilean
Swenson, J. L., y Beck, S. L., Historical 1942 Ecuador and 1942 subduction zones, Journal of geophysical research, 96, 977-
Peru subduction earthquakes, and earthquake cycles along 999, 1991.
Colombia-Equador and Peru subduction segments, Pageoph., Trenkamp, R., Kellogg, J. N., Freymueller, J. T., y Mora, H. P.,
146, 67-101, 1996. Wide plate margin deformation, southern Central America
Taboada, A., Rivera, L. A., Fuenzalida, A., Cisternas, A., Philip, and northwestern South America, CASA GPS observations,
H., Bijwaard, H., Olaya, J., y Rivera, C., 2000, Geodynamics Journal of South American Earth Sciences, 15, 157-171, 2002.
of the northern Andes: Subductions and intracontinental

179