Adaros Rodrigo 0255m

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 90

UNIVERSIDAD DE CHILE

FACULTAD DE CIENCIAS FÍSICAS Y MATEMÁTICAS


DEPARTAMENTO DE GEOFÍSICA

SISMICIDAD Y TECTÓNICA DEL EXTREMO SUR


DE CHILE

RODRJGO EDUARDO ADAROS CÁRCAMO

PROFESOR GUÍA EMILIO VERA S.


PROFESORES DE COMISIÓN MARIO PARDO
DIANA COMTE
PROFESOR INVITADO DOUGLAS WIENS

TESIS PARA OPTAR AL GRADO DE


MAGÍSTER EN CIENCIAS, MENCIÓN GEOFÍSICA

SANTIAGO DE CHILE
Agosto 2003
Resumen
Durante los años 1997 y 1998, 320 eventos sísmicos fueron registrados por una red de
cinco estaciones de banda ancha en el extremo Sur de Chile, con el objetivo de determinar el
grado de microsismicidad presente en la región de Magallanes.
Con la información adquirida se determinó en forma simultánea una estructura óptima de
velocidades unidimensional (velocidad versus profundidad), y una relocalización de los eventos
(hipocentros) acorde con este modelo de velocidades. Como resultado de este procedimiento,
un total de 263 eventos pudo ser relocalizado.
El modelo de velocidades muestra una capa superficial de 32 km. de espesor con una velo-
cidad de 5.7 km/s, sobre una capa más veloz de 8.4 km/s. Esta estructura se interpreta como
la corteza continental que yace directamente sobre el manto superior, situación que se refleja
en las curvas camino-tiempo, donde las fases del manto (Pn) forman las primeras llegadas
para distancias epicentrales mayores a 180 km. Los hipocentros muestran profundidades que
en su gt:an mayoría no sobrepasan los 80 km. La integración de la información hipocentral
con perfiles sísmicos regionales muestran una correlaciÓn entre la sismicidad y sistemas de
fallas normales con profundidades que no sobrepasan los 10 km. lo que lleva a plantear la
existencia de fallas superficiales activas para la región de Magahanes. Sin embargo no es
posible determinar los mecanismos de foco, por lo cual se mantiene en duda la naturaleza de
dichos eventos.
La distribución epicentral de la sismicidad muestra una concentración de eventos al Norte
del Estrecho de Magallanes (Falla de Magallanes), situación que puede ser explicada por la
baja cobertura de estaciones al Sur de dicha falla. A pesar de lo anterior, la sismicidad a
lo largo de ella es tenue, contraviniendo la información entregada por los registros históricos
aportados por la red mundial que, apoyados en soluciones de mecanismos de foco, verifican la
presencia de tres placas (Scotia-Antártica-Sudamérica) interactuando en el Sur de Chile. El
punto específico de la triple unión entre estas placas no queda particularmente evidenciado
en este estudio. No obstante lo anterior, la SJsmicidad obtenida permite ver lineamientos
concurrentes y paralelos a la falla de Magallanes, lo que estaría de acuerdo no con la idea de
una falla concentrada en una sola línea, sino que en una zona de falla más ancha con una
liberación de energía más extendida.
Resulta interesante señalar que un grupo importante de eventos registrados entre los meses
de Abril y Mayo de 1998 muestran hipocentros en las cercanías del volcán Reclus, que al igual
que los demás volcanes de la región se presumen extintos. La presente observación sería una
de las primeras evidencias directas de la existencia de actividad actual asociada al volcán
Reclus.
Tanto por la cobertura de las estaciones como por la baja energía de los eventos, no
fue posible determinar los ejes principales de esfuerzo para la región, obteniendo sólo los
mecanismos de foco de los dos eventos más grandes del periodo que muestran soluciones de
fallamiento inverso, tal cual lo presenta una solución de 1976 obtenida por la red mundial ,
evidenciando de esta forma los procesos de subducción que ocurren entre la placa Antártica
y Sudamericana.
Agradecimientos

Quisiera agradecer al profesor Emilio Vera por haber ayudado a hacer posible esta tesis y
por haberme permitido participar en tan hermoso proyecto, descubriendo zonas tan salvajes
y bellas a la vez.

A los profesores Douglas Wiens y Patrick Shore con quienes pasé largas campañas por los
mares antárticos y por las estepas patagónicas. A la Washington University, St. Louis-Missouri
y en especial a Stacey Robertson, por su cálida hospitalidad.

A los profesores de la comisión y en forma especial al profesor Mario Pardo, quien siempre
me dio el apoyo y el empuje para terminar esta tesis. Al departamento de Exploración de
ENAP Magllanes por permitirme utilizar sus figuras , asi tambien a su grupo de geólogos y
en particular a Sebastián Trejo, por sus comentarios acerca de la geología del Sur de Chile.
En forma especial agradezco al Departamento de Geofísica de la Universida de Chile, a mis
compañeros y funcionarios, con quienes compartí tantas jornadas haciendo más placentera
mi permanencia en esas aulas. A llicardo Zapata por sus consejos de GMT y J5IEX. A Sergio
Aguilera por su ayuda computacional.

Mi trabajo de campo no hubiera sido posible sin el soporte del Instituto Antártico Chileno
y de National Science Foundation (USA). Mi especial agradecimiento a John Evans, logístico
del NSF y a todas aquellas personas que ayudaron a manter en funcionamiento la red temporal
de estaciones en Patagonia y Península Antártica.

Finalmente me gustaría agradecer a mi familia que siempre me ha apoyado, a mi esposa


Gabriela por su constante ayuda y a mi hija Emilia que vino a alegrar la última etapa de este
trabajo.

*
Esta tesis fue financiada por la Comisión Nacional de Investigación Científica y Tecnológica
(CONICYT) con la Beca de Magíster en Ciencias, año 1999.
,
Indice general

l. Introducción 1
1.1. Antecedentes 1
1.2. Objetivos . . 4

2. Geología y Tectónica del Extremo Austral Sudamericano 6


2.1. Introducción .. . . . . . . . . . . . . . . .. . . . . . . 6
2.1.1. Evolución Geotectónica del Margen Chileno .. 8
2.1.2. Evolución Tectónica para el sur de Sudamérica . 12
2.2. Situación Actual .. . .. . . . . . .. . . .. . . . . . 19
2.2.1. Telesismicidad de la región austral de Sudamérica 22
2.2.2. Terremotos ocurridos en la región . . . . . . . . 23

3. Descripción y Procesamiento de Datos Sismológicos 25


3.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . 25
3.1.1. Operación de campo . . . . . . . . . . .. . . . 25
3.1.2. Secuencia de trabajo con la información recolectada 30
3.2. Sismicidad local . .. . . . . . .. . .. . 34
3.2.1. Elección del modelo de velocidad 34
3.2.2. Resultados hipocentrales 36
3.2.3. Estructura Cortical .. . 45

4. Resultados:Análisis y Discusión 49
4.1. Sisrnicidad local durante el periodo del Proyecto SEPA 49
4.2. Mecanismos de foco . . . . . . . . . 51
4.3. Sisrnicidad asociado a volcanismo . . . . . . 59
4.4. Breve reseña del volcanismo austral . . . . . 61
4.4.1. Sismicidad asociada al Volcán Reclus 63
4.5. Perfiles sísmicos regionales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66

5. Conclusiones 70

A. Función de transferencia del sensor STS-2 73

B. Determinación de los parámetros hipocentrales y modelos de velocidad 75


B.l. Tiempo de viaje de un rayo simple . . . . . . . . . . . . . . . . . . . . . . . 75
B.l.l. El problema acoplado: Hipocentro y modelo de velocidad . . . . . . . 76
B.l.2. Separación de los Parámetros Hipocentrales y del modelo de velocidad 78
B.l.3. Solución Numérica al Problema Inverso Separado . . . . . . . . . . . . 78

II
,
lndice de figuras

1.1. Distribución de estaciones que conforman el proyecto SEPA. 3

2.1. Tectónica actual propuesta para la zona de estudio . . . . 7


2.2. Evolución geotectónica del extremo austral de Sudamérica 9
2.3. Mapa geológico y rasgos morfoestructurales . . . . . . . . 10
2.4. Estructura actualmente existente en el extremo austral chileno 11
2.5. Faja corrida y plegada .. . . . . . . . . . . . . . . . . . 12
2.6. Régimen extensivo dominante hacia el Jurásico Superior . . . 13
2.7. Geotectonismo imperante en Cretácico Inferior . . . . . . . . . 14
2.8. Tectónica de las placas Farallon y Phoenix y su posterior conversión a las
placas Nazca y Antártica . . . . . . . . . . . . . . . . . . . . . . . . 16
2.9. Migración de la dorsal chilena desde los 18 Ma hasta la actualidad. 17
2.10. Evolución de los modelos tectónicos para la región de Scotia .. . . 21
2.11. Telesismos registrados en los bordes de la placa de Scotia y alrededor del triple
punto Antártica-Scotia-Sudamericana . . . . . . . . . . . . . . . . . . . . . . . 22
2.12. Epicentros estimados para los dos grandes eventos del día 17 de Diciembre de
1949 . . . . . . . . . . . . . . . . . . . . 24

3.1. Respuesta en velocidad del sensor STS-2 26


3.2. Cronograma de funcionamiento de las estaciones SEPA 27
3.3. Sismogramas sin filtro para un evento local .. . . . . 28
3.4. Sismogramas con filtro pasa alto para un evento local . 28
3.5. Sismogramas sin filtro para un evento telesísmico . . . 29
3.6. Sismogramas con filtro pasa bajo para un evento telesísmico 29
3.7. Secuencia de trabajo con los datos del experimento SEPA. . 30
3.8. Epicentro del mayor evento registrado en 1997 por las estaciones del proyecto
SEPA . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32

III
3.9. Sismograma asociado al mayor evento del año 1997 . . . . . . . . . . . . . . . 32
3.10. Epicentro del mayor evento registrado en 1998 por las estaciones del proyecto
SEPA . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
3.11. Sismograma asociado al mayor evento del año 1998 . . 33
3.12. Últimos modelos de velocidad entregado por VELEST . 35
3.13. Sismicidad total registrada por la red temporal SEPA . 37
3.14. Espesores de corteza aproximado, mostrando espesores del orden de 30-40 km.
para el sur de Sudamérica . . . . . . . . . . . . . . . . 45
3.15. Modelo de dos capas. . . . . . . . . . . . . . . . . . . . . 46
3.16. Distancia epicentral vjs tiempo. Componente vertical. .. 47
3.17. Distancia epicentral vjs tiempo. Componente horizontal. 48

4.1. Principales lineamientos generados por SEPA . . . . . . 50


4.2. Inversión para la componente radial en estación FELL 52
4.3. Inversión para la componente transversal en estación FELL . 53
4.4. Inversión para la componente vertical en estación FELL . . . 54
4.5. Inversión para la componente radial en estación MILO . . . 55
4.6. Inversión para la componente transversal en estación MILO . 56
4.7. Inversión para la componente vertical estación MILO . . . . 57
4.8. Soluciones de mecanismos de foco existentes para la Región de Magallanes 58
4.9. Cadena de volcanes de Los Andes australes . . . . . . . . . . . . . . . . . . 60
4.10. Sismicidad superficial ( <3 km profundidad) y su relación con el volcanismo 64
4.11. Sismograma típico de los eventos superficiales del sector del volcán Reclus . 65
4.12. Ubicación de cuatro secciones regionales. . . . . . . . . . . . . . . . . . 67
4.13. Interpretación estructural para tres perfiles de la región de Magallanes . 68
4.14. Detalle de la sismicidad en la sección B-B' . . . . 68
4.15. Sección sísmica interpretada del sector Atlántico . . . . . . . . . . . . . 69

IV
,
Indice de cuadros

1.1. Coordenadas de estaciones del proyecto SEPA . . . . . . . . . . . . . . . . . . 2

2.1. Coordenadas estimadas de los dos mayores eventos registrados durante el siglo
XX. Evento 1 posee una elipse de confianza (eje mayor:l12 km , eje menor:57
km, azimut:N78°W). Evento 2 posee una elipse menor a la primera (eje ma-
yor:62 km, eje menor:31 km, azimut:N75°W). Evento 3 es el mismo al evento
2, pero calculado por otro autor. . . . . . . . . . . . . . . . . . . . . . . . . . . 24

3.1. Modelo de velocidad obtenido de experimentos de refracción sísmica en el Es-


trecho de Magallanes por Ludwing et al. [1965] . . . . . . . . . . . . . . . . . . 34
3.2. Modelo final de velocidad, obtenido por sucesivas iteraciones de minimización
de los residuales utilizando el programa VELEST [Kissling, 1988] 35
3.3. Archivo final de localizaciones . . . . . . . . . . . . . . . . . . . . . . . . . . . 38

4.1. Mecanismos de foco obtenido del catálogo Harvard-CMT [2002] para los ma-
yores eventos registrados desde el año 1979. . . . . . . . . . . . . . . . . . . . 59

V
Capítulo 1

Introducción

1.1. Antecedentes
En los últimos dos siglos se ha avanzado en la descripción y estudio de la actividad sísmica,
logrando reconstruir la historia geológica del continente con un marcado acento en el norte y centro
de nuestro país. En cambio, de la sismología y tectónica del sur de Chile es muy poco lo que se sabe,
no tan sólo por su aislamiento, sino por su abrupta geografía, marcada por la presencia de profundos
fiordos , canales y lagos. Adentrarse por estas regiones es, sin duda, una suerte de azar y fortuna.
La actividad sísmica del margen sur de Sudamérica se conoce sólo a través de los registros
generados a partir de la red mundial de estaciones sismológicas, donde es posible apreciar una baja
sismicidad. Diversos estudios reconocen que la mayor parte de la actividad sísmica del margen
costero de Chile corresponde a un proceso activo, en que la subducción y fusión parcial de la placa
oceánica ha dado lugar a una gran actividad magmática desde el Paleozoico a la actualidad [Ramos,
1989]. Sin embargo, en el sur de Chile (sur de la península de Taitao) el proceso difiere de las otras
zonas estudiadas: la subducción se produce por el encuentro de tres placas tectónicas, originando
de esta forma un complejo modelo geodinámko.
Es así como surge la idea de desarrollar un registro específico para esta zona y el proyecto SEPA se
materializa como un modo de estudio sismológico más acabado para la región. Durante 1997 y 1998,
cinco sismógrafos de banda ancha y de tres componentes fueron instalados en la Patagonia Chilena,
junto con ocho sismógrafos de las mismas características en Islas Shetland del Sur y Península
Antártica. Esta red formó parte del Experimento Sísmico en la Patagonia y Antártica (SEPA),
experimento que fue un esfuerzo conjunto entre la Universidad de Chile y Washington University,
Saint Louis (MO) USA.
La red local de estaciones en Patagonia cubrió parte de la región de Magallanes (figura 1.1) , lo

1
l . Introducción 1.1 Antecedentes

que significa un experimento único y pionero para una zona muy poco estudiada.
A finales de 1998 se retiraron las cinco estaciones de la zona de estudio, reinstalando solamente
una, ubicada en las cercanías de Punta Arenas y que se encuentra operativa hasta esta fecha.

Estación Latitud Longitud Elevación


ELEF -61.17 -55. 10 0.0
DECP -62.98 -60.67 0.0
FREI -62.19 -58.98 0.0
ome -63.33 -57.90 0.0
PRAT -62.47 -59.66 0.0
LOWI -63.25 -62. 18 0.0
SPPT -64.30 -61.05 0.0
MILO -51.57 -72.62 0.0
FELL -51.57 -72.62 0.0
SALM -52.55 -72.03 0.0
HAMB -53.61 -70.93 0.0
VTDF -54.14 -68.71 0.0

Cuadro 1.1: Coordenadas de estaciones del proyecto SEPA

2
l. Introducción 1.1 Antecedentes

·52"

Ore* e PCl!is~

·52"

Figura 1.1: DU;tribuci6n de r.t;tacionetl lrismológir.as que formaron parte del proyecto SEPA

3
l . Introducción 1.2 Objetivos

1.2. Objetivos
El objetivo principal de esta tesis es determinar la microsismicidad de la región de Magallanes
para poder contribuir al conocimiento de la tectónica del extremo sur de Sudamérica y, a su vez,
establecer sus posibles relaciones con el riesgo sísmico de la región. Para ello se utilizaron datos
sismológicos de una red local y mundial, conjuntamente con datos geológicos.
Tal como se indica en el título, este trabajo está centrado en el extremo sur de Sudamérica,
específicamente en la Región de Magallanes, donde los parámetros que controlan el proceso de
subducción de las placas Antártica y Sudamericana permanecen prácticamente constantes a lo largo
del margen, a excepción de la oblicuidad. Además, factores tales como la baja actividad sísmica
y volcánica registrada, así como los pocos estudios previos realizados en la zona han sido una
constante en el desarrollo de esta tesis. De los objetivos generales se desprenden los siguientes
objetivos específicos:

• Determinación de estructura de velocidades de propagación de ondas sísmicas {elásticas) pri-


marias (P) y secundarias (S) para la zona de estudio: En este punto se utilizarán algoritmos
de inversión simultáneo para velocidaes en 1-D (Vp y Vs) y parámetros hipocentrales.

• Determinación de soluciones bipocentrales: Una vez determinado el mejor modelo de velocida-


des para la zona, se procederá a la determinación de los hipocentros para los eventos detectados
por la red.

• Correspondencia de las soluciones bipocentrales a la actividad volcánica de la región.

• Determinación de los mecanismos de foco: A partir de los objetivos 1 y 2 se determinarán las


posibles soluciones focales sólo de aquellos eventos de magnitudes {mb) mayores a 4.5.

Todo lo anterior permite un desarrollo de un total de cinco capítulos que se detallan a continuación:
Capítulo 1: Se trata de una introducción general del tema que incluye los Objetivos Generales y
específicos de la tesis. Junto con lo anterior se explica la materialización del Proyecto SEPA como
una necesidad de estudiar la zona a nivel local.
Capítulo 2 : Reúne los antecedentes generales de la zona estudiada que incluyen el marco sismo-
tectónico, los modelos geodinámicos para la región y su respectivo marco Geológico.
Capítulo 3 : Describe la secuencia de trabajo desde la adquisición hasta el procesamiento de
datos sismólogicos, tales como los sistemas de operación de los instrumentos utilizados y los aspectos
teóricos de las localizaciones hipocentrales.
Capítulo 4 : En este punto se reúnen todos los resultados obtenidos en el proceso para su análisis
y discusión.

4
l. Introducción 1.2 Objetivos

Capítulo 5 : Se concluye los aspectos más relevantes originados por los resultados del experimento.

5
Capítulo 2

Geología y Tectónica del Extremo


Austral Sudamericano

2.1. Introducción
En este capítulo se presenta un resumen de la evolución geodinámica del margen sur de Su-
damérica, específicamente entre los 50 os y 57 °S.
Hacia el margen occidental de la placa Sudamericana existen diferentes terrenos tectónicos orien-
tados en dirección N-S, definiendo en el extremo austral de Sudamérica tres placas: Sudamericana,
Antártica y la placa de Scotia [Ramos, 1989] (Figura 2.1). En sus bordes Norte y Sur, la placa de
Scotia está en contacto con las placas Sudamericana y Antártica a través de las dorsales Norte y Sur
de Scotia respectivamente (Pelayo y Wiens, 1989] . Estos bordes presentan un movimiento relativo
predominantemente transcurrente sinestral, con una pequeña componente de compresión. El límite
Oeste de la placa de Scotia se ubica en la zona del paso Drake, entre Sudamérica y la Península
Antártica y presenta movimientos relativos de convergencia y transcurrencia. En el extremo Sur de
Sudamérica, la placa Antártica se subducta por debajo de la de Scotia. Hacia el Sur, sin embargo,
esta convergencia se hace cada vez más oblicua terminando en un movimiento relativo de transcu-
rrencia a lo largo de la zona de fractura de Schakleton, la que finalmente intercepta la dorsal Sur de
Scotia cerca de Isla Elefante (islas Shetland del Sur) .
El extremo norte de la placa de Scotia queda bien definido con una orientación Este-Oeste
que sigue la dorsal norte de la placa de Scotia. Evidencias geológicas [Winslow, 1982] indican que
hacia el Oeste sigue la Falla de Magallanes, aproximadamente a lo largo del Lago Fagnano, Seno
Almirantazgo, y brazo Oeste del Estrecho de Magallanes, hasta su intersección con la fosa Chilena
(52 os, 76 °W) donde definiría un punto triple entre las placas Sudamericana, Antártica y de Scotia

6
2. Goolog{a y Tectónica del Extremo Austral Sudamericano 2.1 Introducción

{Figura 2.1). Cabe señalar que hasta el día de hoy dicho punto triple no está del todo defirúdo, y su
ubicación actual es sólo una prolongación de la Ial.la de Magallanes hacia la placa Antártica.
El área de estudio {recuadro Figura 2.1) involucra parte del territorio chileno y otra parte del
territorio argentino, sin embargo, se utilizaron sólo datos registrados íntegramente en territorio
nacional.
-80°
-10° _j

1
, ~

NAZCA

-60°¡ ¡-600
-80° -600
Figura 2.1: Tectónica actual de la zona Sur de Sudamérica. Recuadro indica área de
estudio.

7
2. Geología y Tectónica del Extremo Austral Sudamericano 2.1 Introducción

2.1.1. Evolución Geotectónica del Margen Chileno

Para todo el margen chileno, la subducción de corteza oceánica bajo la placa Sudamericana dio
lugar a una actividad magmática importante, afectando de manera considerable su margen occi-
dental [Ramos, 1989]. Según Mpodozis y Ramos [1990], durante el paleozoico inferior el margen
suroccidental de Gondwana estaba cerca del límite oeste de las actuales sierras pampeanas argen-
tinas. Al oeste del continente de Gondwana existe una plataforma carbonatada que limitaba con el
dominio oceánico, que a su vez, separaba el continente de Gondwana del rnicrocontinente de Chile-
nía. La subducción del dominio oceánico y la posterior acreción de Chilenia a Gondwana se produjo
en el Devónico Superior, conformando un complejo acrecional que actualmente se conoce con el
nombre de Cordillera de la Costa, en donde se desarrolló un importante arco volcán ico asociado a
una cuenca de tras arco. A finales del Paleozoico finalizó este proceso de acreción lateral, el cual
resurgiría hacia el Jurásico a lo largo de todo el país, definiendo según Mpodozis y Ramos [1990]
cinco segmentos de acuerdo a la evolución tectónica particular de cada uno de ellos. Para los autores
antes mencionados, el segmento más austral de ellos (49°- 56°S ) evolucionó como una cuenca de
antearco con basamento de tipo oceánico, conocido con el nombre de cuenca marginal de Rocas
Verdes, formada al Este del arco, en el Cretácico Inferior. Se reconoce la formación de un complejo
metamórfico producto de una inversión y cierre de la cuenca en el Cretácico Medio (Figura 2.2).
Importantes eventos deformativos afectaron a la región andina durante el Cretácico Medio y
Superior, proceso en el que los distintos segmentos respondieron de forma particular en cuanto
desarrollaron djstintos sistemas de fallas y pliegues de un segmento a otro. Durante el Mioceno, la
región Andina experimentó una nueva reorganización que condujo a la actual división morfotectórúca
controlada aparentemente por factores tales como la velocidad de convergencia entre las placas y la
inclinación de la zona de Benioff. Como consecuencia de estos factores se explica la distribución del
volcanismo Cenozoico, el levantamiento del altiplano y las sierras pampeanas de Argentina y, en la
misma medida, todos los cambios tectónicos posteriores que se pueden reconocer en la actualidad,
en particular el sistema de Faja Corrida y Plegada y la Cuenca de Antepaís del extremo a ustral
chileno (Figuras 2.3 y 2.4).

8
2. Geología y Tectónica del Extremo Austral Sudamericano 2.1 Introducción

JUilÁSICO MEDIO
( AalenianoCaiOIIionol

t·¡::
1

~
Q) JURÁSICO SUPERIOR
"O ( Oxfordlano -ntonlono
oo_
.2
L1J

t
·¡::
CRETÁCICO INFERIOR
( Beniasiono - Aptiano )
-t;
oo_
Q)
"O
oo_
.2
L1J

CRETÁCICO SUPERIOR
( Alblono 1 - Cenoman\on4>

, 5!_?
o
o_
CRETÁCICO SUPERIOR
2e ( Turoniano 1 - Maaslrichtlono

<{
Q)
"O
oo_
.2
UJ
TERCIAAIO
( PoleocenoPIIoceno )

200 Km.

Figura 2.2: Evolución geotect ónica del extremo Sur de Sudamérica, desde J urásico Medio
hasta el Terciario

9
2. Goolog{a y Thctónica del Extremo Austral Sudamericano 2.1 Introducción

-50°

D Cretóclco

D Jurósico

D Basamento ¿Paleozoico?

• Batollto Patagónico
(Meso-Cenozoico)

-70°

Figura 2.3: Mapa Geológico y rasgos moríoestructurales del extremo Sur de Sudamérica
(Diraison et al., 1997; Klepeis y Austin, 1996].

10
2. Geolog{a y Tectónica del Extremo Austral Sudamericano 2. 1 Introducción

1+- Antefosa ----- ---------------------------· Anlepais


Cinturón
Plegado
o E

so N

Figura 2.4: Estructura actual existente en el extremo austral de Chile

11
2. Geolog{a y Thctónica del Extremo Austral Sudamericano 2. I Introducción

2.1.2. Evolución Tectónica para el sur de Sudamérica

El desarrollo tectónico del margen sur de los Andes tiene su inicio en el Paleozoico Superior, o
tal vez antes, pues se reconoce la formación de un complejo acrecional en el margen occidental de
Gondwana (lado Pacífico) que es lo que actualmente constituye el basamento de la zona en estudio.
Hervé et al. [1981] , Klepeis y Austin (1996] dataron rocas que afloran en el sector de la faja plegada y
corrida del extremo austral, y que cabalgan sobre rocas sedimentarias de edades jurásicas superiores
{Figura 2.5), las que sugieren una edad de 224 ± 38 Ma, ubicando al basamento entre el Paleozoico
Superior y el Mesozoico Inferior.

20 Km. de acortamiento
horizontal mínimo

LEYENDA
O Oligoceno
O Eoceno
• Paleoceno
O Cret . sup .
• Cret . inf.

................._...__._....... lO km

Figura 2.5: Faja corrida y plegada del ~>.xtremo Sur-OccidP.ntal de Sudamérica

Un importante cambio paleogeográ.fico a lo largo de todo el margen Pacífico se produjo tras un


periodo de levantamiento y erosión. Por una parte, el arco magmá.tico se trasladó hacia el oeste,
intruyendo el complejo acrecional Paleozoico. Como consecuencia de lo anterior se reanuda la sub-
ducción del dominio oceánico entre el Jurásico Superior y el Terciario [Hervé et al., 1984; Bruce et
al.., 1991], constituyendo el batolito calcoalcal.ino de la Patagonia o batolito patagónico {Figura 2.3)
que se extiende desde el Cabo de Hornos (56°8) hasta el norte (50°S) con un ancho promedio de
60 km. Dicho cuerpo granodiorítico se menciona como referencia de la raíz de la cadena volcánica
andesítica iniciada en el Jurásico Superior. Por otra parte, la ruptura de Gondwana inicia el desa-
rrollo de un régimen tectónico extensivo [Dalziel, 1981; Klepeis y Austin, 1996) en el que se produce
volcanismo silicico y el adelgazamiento del cratón Sudamericano por medio de fallas normales, junto
con el desarrollo de la cuenca de Magallanes [Wilson, 1991; Klepeis y Austin, 1996), figura 2.7. Este
régimen acompañado de un importante pulso magmático durante el Jurásico Superior y Cretácico

12
2. Geolog(a y Thctóoica del Extremo Aust.ral Sudamericano 2.1 Introducción

Inferior propicia la intrusión de magma basáltico en la corteza continental al Este del complejo
volcánico que, como consecuencia, fue desplazado hacia el pacífico origina.ndo una cuenca. marginal
con basamento oceánico ubicada entre el arco volcánico y el continente Sudamericano [Cunningham,
1995; Da.lziel et aL, 1974].

JURASICO
SUPERIOR
Oxfordlano
(150 m.a.)

u.-•costa
Sprlnghlll

Dorslll

Arco
MagrMtlco

Figura 2.6: Regimen extensivo dominante hacia el Jurásico Superior

La cuenca. marginal de tras arco conocida con el nombre de Rocas Verdes es rellenada con
sedimentitas del arco activo, rocas volcánjca.s silicias, sedjmentos volca.noclásticos del Jurásico Medio-
Superior [Da.lziel et al., 1974] , y depósitos del Cretácico Inferior. En la actualidad, la. localización
de las rocas asociadas al complejo de Rocas Verdes se ubica. desde los 50° S hasta el Cabo de Hornos,
abarcando unos 700 km, a los cuales hay que añadir unos 50 km sobre la plataforma continental de
Sudamérica al este del Cabo de Hornos, junto con 350 km más que se encuentran sobre la plataforma
sumergida. de Georgia del Sur. Cabe recordar que aunque esta isla. se encuentra a 1700 km de
distancia formaba parte del bloque de tras arco de Sudamérica. antes de producirse el movimiento
transformante de la dorsal norte del mar de Scotia [Dalziel, 1981].
A lo largo del margen occidental se produjo un aumento de la convergencia Este-Oeste después
de la apertura definitiva del Atlántico sur en el Cretácico Inferior. El surgimiento de un componente
de movimiento sinestral Este-Oeste de las placas de la región sur se produjo junto al desarrollo de la
placa Sudamericana, separada de la placa. Antártica. y con una deriva hacia el Oeste mayor que esta
última con respecto a la placa. Africana [Cunningham et al., 1991; Cunningham, 1995], en la medida

13
2. Geologfa y Tectónica del Extremo Austral Sudamericano 2.1 Introducción

CRETACICO
INFERIOR
Barren lano
(120 m.a)

Unudecosta
Sprlnghlll

1 Dorsal

0 Arco
Magmitlco

Figura 2.7: Geotectonismo imperante en el Cretácico Inferior

en que empezó a separarse de la futura península antártica. Cunningham explica la existencia de


ambos fenómenos como:
a) el inicio del cierre de la cuenca marginal de las Rocas Verdes en su parte meridional, b) el
inicio de la subducción del basamento oceánico de la cuenca por debajo del margen occidental y e)
el inicio de la formación de la Cordillera Darwin [Cunningham, 1995].
El inicio del desarrollo de la cuenca de antepaís se produce con el depósito de los primeros
sedimentos provenientes del nuevo orógeno en la cuenca de Magallanes [Klepeis y Austin, 1996], y
cuya subducción origina un importante plegamiento y metamorfismo de las rocas que conforman
la cubierta de la cuenca. Tal proceso se acrecentó en los inicios del Cretácico Superior en la zona
de la cuenca de Rocas Verdes y en la actual cordillera de Darwin, proceso que fue acompañado de
una migración de la compresión hacia la zona de antepaís. La subducción alcanza su punto crítico a
mediados del Cretácico superior cuando subduce la raíz oceánica de la cuenca; se produce un bloqueo
de la entrada de la corteza continental en la zona de subducción, lo que causa la ruptura de la raíz
oceánica que subducía, produciendo el rápido levantamiento de la Cordillera Darwin como un modo
de restablecer el equilibrio isostático [Cunningham, 1995]. Tal levantamiento fue a su vez favorecido
por la rotación en sentido antihorario de la región más meridional de los Andes relacionada con el
orógeno patagónico.

14
2. Geología y Tectónica del Extremo Austral Sudamericano 2.1 Introducción

Según lo antes visto, el plegamiento de la cordillera estaría ligado a un sistema de fallas direc-
cionales con movimiento siniestra!, tanto en el extremo sur de Sudamérica como en la península
Antártica. Tal sistema de fallas responde a la zona de fractura de Schakleton, que poseen una orien-
tación WNW en que se habría sintetizado un movimiento de componente EW de tipo direccional
siniestra} y el anterior N-S divergente. La zona de fractura de Schakleton habría actuado como límite
tectónico entre la zona deformada del sur de Sudamérica y la parte norte de la península Antártica
no deformada. Sin embargo, Cunningham señala que el proceso de separación entre ambas zonas
y que dio lugar a la formación de la corteza oceánica en la parte Oeste del mar de Scotia se ace-
leró entre los 55 y 40 ma cuando el movimiento relativo entre Sudamérica y África pasó a ser de
dirección WNE mientras que el mismo movimiento entre Antártica y África era de dirección WSE.
Cunningham et al. señala que desde los 50 ma hasta la apertura del paso del Drake (aprox 30 ma)
el desplazamiento EW de tipo direccional siniestra} entre el sur de Sudamérica y el norte de la
Antártica es de unos 290 km, mientras que su desplazamiento divergente NS es de 150 km.
Por otro lado Cande y Leslie [1986] explican la situación tectónica del extremo sur de Sudamérica
y península Antártica como producto de una colisión entre las placas de Phoenix-Farallon (prede-
cesora de la placa de Nazca) y la placa Sudamericana durante el Cenozoico Inferior (Figura 2.8).
Como consecuencia de la colisión, el punto triple resultante migró rápidamente hacia el sur durante
el Eoceno, localizándose en la zona de Tierra del Fuego hace unos 40 Ma.
La apertura del paso Drake fue originada, de esta forma, por la separación gradual del sur de
Sudamérica con la Península Antártica [Cunningbam, 1995]. El fondo oceánico se expande hacia
la zona occidental del mar de Scotia. En este punto, el movimiento siniestra! con orientación Este-
Oeste junto con el divergente Norte-Sur fue reacomodado en la región central y occidental del mar de
Scotia, cuyo límite occidental de orientación Noreste lo constituye la zona de fractura de Schakleton.
Cunningham explica tal fenómeno por medio de una rotación en sentido antiborario del extremo
meridional del continente Sudamericano.
La subducción de la dorsal chilena bajo la placa sudamericana se estima que comenzó hace 14
Ma. Cande y Leslie [1986] sugieren que la colisión entre la dorsal chilena y la fosa chilena se inicia al
Oeste de Tierra del Fuego (Figura 2.9}. El triple punto resultante de esta colisión tuvo una migración
rápida hacia el Norte, situándose al sur del Golfo de Penas hace unos 10 Ma. Entre 6 Ma y 3 Ma
secciones más cortas de la dorsal subdujeron a lo largo del margen continental. Posteriormente, hace
unos 100.000 años se produce la colisión de la dorsal chilena situada entre las zonas de fractura
de Taitao y Darwin con la fosa chilena. Interpolando este mismo tipo de movimiento al futuro,
podrá estimarse que al Norte de la fractura de Darwin se generará una colisión similar dentro de
un millón de años. Cande y Leslie señalan como consecuencia de estas colisiones la formación de
centros de expansión y fallas transformantes de la dorsal chilena.

15
2. Geología y Tectónica del Extremo Austral Sudamericano 2.1 Introducción

FARALLON
(NAZCA)

ANTARTICA

-40Ma

Figura 2.8: Tectónica de 40 Ma donde se observa la subducción y posterior desaparición


de las Placas de Farallon y Phoenix originando en forma posterior la actuales
placas de Nazca y Antártica (modificado de Cande y Leslie [1986]).

16
2. Geología y Tectónica del Extrem o A ustra/ Sudamericano 2.1 Introducción

Figura 2.9: Migración de la dorsal chilena desde 18 Ma hasta su ubicación actual [Cande
y Leslie, 1986].

17
2. Geolog{a y Tectónica del Extremo Austral Sudamericano 2.1 Introducción

En Resumen, la convergencia de placas en el periodo Cenozoico se puede reducir a tres fases:


a) La subducción de la placa de Phoenix bajo la placa Sudamericana y Antártica que colisiona en
forma oblicua a la dorsal Phoenix-Farallon (Nazca).
b) La subducción de la placa de Nazca bajo la Sudamericana que termina con la colisión de la
actual dorsal chilena (Nazca-Antártica) con el margen continental a unos 90 mm/ a (Cande y Leslie
1986; De Mets et al, 1990; De Mets et al , 1994).
e) La subducción de la placa Antártica bajo la placa Sudamericana a una razón de 15 a 20 mm/ a
al norte de 52°S y oblicua bajo la placa de Scotia a unos 10-15 mm/a al sur de 52°S [Cande y Leslie,
1986; DeMets et al., 1990, 1994]

18
2. Geología y Tectónica del Extremo Austral Sudamericano 2.2 Situación Actual

2.2. Situación Actual


La situación tectónica actual está determinada por la interacción de las tres placas: la Sudameri-
cana, la Antártica y la de Scotia. Uno de los primeros autores en formular un modelo tectónico para
la zona del mar de Scotia fue Barker [1970], quien identificó el límite occidental del mar de Scotia
con la zona de fractura de Shackleton y el límite oriental con las islas Sandwich del Sur (Figura
2.10a). Para el autor, dichas islas corresponderían a un arco constituido por un centro de expansión
de trasarco hacia el oeste y por una zona de subducción hacia el Este. Los correspondientes límites
norte y sur también logran ser definidos: el límite norte se caracteriza por un conjunto de centros de
expansión y fallas transformantes que comienzan alrededor del punto medio de la zona de fractura
de Shackleton que va desde los 56°S de latitud y 52°W de longitud , desde este punto en adelante el
límite se prolonga hacia el sudeste como una zona de fractura.

Forsyth [1975] es quien propone y define la placa de Scotia como tal. El autor señala que las
placas Sudamericana y Antártica tienen un movimiento sinestral a través de la dorsal Norte y Sur
del mar de Scotia las que posteriormente logran definir una dorsal de orientación S-W que se pro-
pondría como el límite Oeste de una placa menor que serán los límites Norte y Sur de la misma
placa, sin embargo, no se propone un límite claro para el borde Oeste. Al sudoeste de la Zona de
Fractura de Shackleton se reconoce una posible dorsal de orientación S-W que se propondría como el
límite Oeste de una placa menor que será luego subducida por debajo de la Peninsula Antártica en
su borde Oeste (Figura 2.10b). El limite Norte presentaría una falla sinestral de orientación E-W a
lo largo de la Dorsal Norte del mar de Scotia que se prolongaría hasta el margen chileno, definiendo
el triple punto creado por la intersección de las placas Sudamericana, Antártica y Scotia.

Actualmente, autores como Pelayo y Wiens [1989] han propuesto un movimiento sinestral de
strike-slip con una componente de compresión que aumenta hacia el Este a lo largo de la dorsal
(Figura 2.10c), basándose en mecanismos focales de terremotos ocurridos en la zona , igualmente,
los autores señalan que en su límite Este.,las placas de Scotia y Sandwich están separadas por un
centro de expansión de trasarco.
Las placas de Scotia y Antártica mantienen un contacto de 800 km con una orientación SE
siguiendo la Zona de Fractura de Shackleton, pasando de una subducción oblicua [Cande y Leslie,
1986] a un movimiento de falla de rumbo [Forsyth, 1975; Pelayo y Wiens, 1989]. Al norte de la
Zona de Fractura de Shackleton se encontraría el difuso punto triple de intersección de las placas
Sudamericana, Antártica y Scotia.
Tomando la placa Sudamericana, la placa Antártica, la placa Scotia y la placa Sandwich del Sur

19
2. Geología y Tectónica del Extremo Austral Sudamericano 2.2 Situación Actual

y fijando el polo de Euler entre las placas Sudamericana y Antártica tal como lo dice el modelo
NUVEL-1 [DeMets et al., 1990, 1994], Pelayo y Wiens [1989] obtienen un velocidad relativa de 5
mm/año con movimiento sinestral (componente compresiva) a lo largo de la Dorsal Norte de Scotia
y de 10 mm/año con movimiento sinestral (componente extensiva). Para el sector de la Zona de
Fractura de Shackleton se tiene velocidades de convergencia relativa entre las placas Antártica y
Scotia de 11 y 13 mm/año.
Al norte del punto triple de las placas Sudamericana, Scotia y Antártica se encuentra que la
velocidad de convergencia entre las placas Antártica y Scotia es de aproximadamente 20 mm/ año
[DeMets et al. , 1990, 1994].
Lo anterior es equivalente a ver las velocidades en un sistema absoluto, lo que da una velocidad
de 14 mm/año para la placa de Scotia. Por otro lado se tiene que la velocidad de convergencia
absoluta de la placa Antártica es de 5 mm/año con dirección hacia el SW.
Estudios geodésicos efectuados en el área indican que la velocidad relativa entre las placas Su-
damericana y Scotia es de ...... 6.5 mm/año [Smalley et al. , 2003], de igual forma, logran verificar la
velocidad de 12 mm/año entre la placa Sudamericana y la placa Antártica.

20
2. Geología y Tectónica del Extremo Austral Sudamericano 2.2 Situación Actual

Sudamé~ca

....

-so o

Sudamérica
~ C) 'Voo

Figuza 2.10: Evolución de los modelos tectónicos para la zona de estudio.A)Propuesto


por Barker [1970]; B) Modelo tectónico propuesto por Forsyth [1975] y C)
Modelo actual del extremo Sur de Sudamérica, Placa de Scotia y parte
norte de la Pe.n ínsula Antártica. DNS: Dorsal Norte de Scotia; DSS: Dorsal
Sur de Scotia; ZFS: Zona de Fractura de Shackleton; A indica zonas de
subducción. F lechas negras indican velocidad (mm/ año) absoluta de las
respectivas placas [Pelayo y Wiens, 1989] y flechas grises indican velocidad
relativa entre placas [Gripp y Gordon, 1990].

21
2. Geologfa y Thct<Snica del Extremo Austral Sudamericano 2.2 Situación Actual

2.2.1. Telesismicidad de la región austral de Sudamérica

A nivel mundial, el extremo sur de Chlle ha sido considerado una zona de baja sisrnicidad, debido
principalmente a la escasa información registrada por la red mundial de estaciones sismológicas. Tal
como se observa en la figura (Figura 2.11), los mayores eventos asociados al área de estudio se
concentran en la dorsal norte de la placa de Scotia, zona de Fractura de Schackleton y en mayor
número a la rnicroplaca de Sandwich del Sur. Por tratarse de eventos telesísmicos es necesario que
la energía liberada sea mayor o igual a 4 mb para ser detectado, indicando de esta forma que sólo
una pequeña fracción de los eventos que ocurren en la Patagonia son realmente registrados por las
estaciones de la red mundial.

. ·.


•'

....
....
"' '
,. . · "aW:·
···-·~....
. ..
. . ·,_#';,·· ·r-:.

Figura 2.11: Telesismos registrados en los bordes de la placa de Scotia y alrededor del
triple punto: Aotártica...Nazca..Sudamérica (Fuente: NEIC (2002]).

22
2. Geología y Tectónica del Extremo Austral Sudamericano 2.2 Situación Actual

2.2.2. Terremotos ocurridos en la región

La baja tasa de movimiento entre las placas de Scotia y Suda.mérica origina un bajo nivel de
sismicidad a lo largo del límite de dichas placas [Pelayo y Wiens, 1989], sin embargo durante el pasado
siglo se han reportado grandes eventos que indicarían actividad en la falla de Magallanes-Fagnano.
Los reportes má.s antiguos de actividad sísmica en Tierra del Fuego proceden de la llegada de los
primeros colonos europeos, quienes en Febrero de 1879 reportaron un terremoto que fue percibido
con mayor fuerza en la ciudad de Ushuaia (Sur de Tierra del Fuego).
El 17 de Diciembre de 1949, fueron percibidos dos grandes eventos tanto en la Patagonia Sur
como en la isla Tierra del Fuego, ambos con una diferencia de 8 horas [Lomnitz, 1970]. Sólo el primer
evento tuvo una magnitud reportada por ISC (7.5). La localización hecha por ISC, y relocalización
hecha por Pelayo y Wiens [1989] son pobres debido a la poca cantidad y calidad limitada de registros
sísmicos, fijando dicho epicentro fuera de la denominada falla de Magallanes-Fagnano, sin embargo
estudios efectuados por los mismos autores involucran la falla a una corta distancia de la elipse
de error calcuJada (Tabla 2.1), por lo que se puede asumir el epicentro en las cercanías de la boca
occidental del estrecho de Magallanes sobre la mencionada falla. Reportes de la época señalan que se
produjo licuefacción a 170 km al norte de Río Grande (Tierra del Fuego) y a 50 km al sur de Punta
Arenas [Lomnitz, 1970]. Desplazamientos horizontales de 5 metros fueron reportados por algunos
trabajadores de estancias (granjas) a lo largo de la costa Atlántica donde la falla Magallanes-Fagnano
cae al mar.
El segundo evento tuvo una intensidad VI en Punta Arenas sin registro de magnitud. La mayoría
de los daños fueron como consecuencia de la caída de techumbre y muros de edificaciones. Informes
indican que se produjeron unas 24 réplicas durante el primer día [Lomnitz, 1970]. Dicho documento
también da cuenta de la muerte de tres personas en el sector de San Nicolás a unos 70 km al sur
de Punta Arenas. Grandes olas fueron observadas en el sector de Porvenir y Seno Almirantazgo. El
epicentro (Figura 2.12) fue estimado al Oeste del Estrecho de Magallanes, en la continuación del
Seno Almirantazgo. Aún permanece en duda la ubicación de dicho terremoto, principalmente por la
escasa información de registros sísmicos ya que en aquellos años no existía una buena cobertura ni
estaciones de buena calidad como ocurre en nuestros días.
Smalley et al. [2003] indican que si aceptamos el reporte de 5 metros de desplazamiento, el que
resulta ser razonable para la magnitud del primer evento y una velocidad GPS de 6.5 mm/año, el
periodo de ocurrencia de grandes terremotos asociados a la falla Magallanes-Fagnano sería del orden
de 800 anos, sin embargo la incerteza de dicha estimación es muy grande.

23
2. Geologfa y Tectónica del Extremo Austral Sudamericano 2.2 Situación Actual

Evento Fecha T Origen Lat os Loo ow Nro arr. Mag Int Fuente
1 Dic. 17, 1949 065337 53.40 69.22 39 7.75 (Pelayo y Wiens, 1989]
2 Dic. 17, 1949 150754 53.99 68.77 40 (Pelayo y Wiens, 1989]
3 Dic. 17, 1949 150600 54.00 71.00 VII (Lo!IliÜtz, 1970)

Cuadro 2.1: Coordenadas estimadas de los dos mayores eventos registrados durante el
siglo XX. Evento 1 posee una elipse de confianza (eje mayor:ll2 km, eje
menor:57 km, azimut:N78°W). Evento 2 posee una elipse menor a la primera
(eje mayor:62 km, eje menor:31 km, azimut:N75°W). Evento 3 es el mismo
al evento 2, pero calculado por otro autor.

Figura 2.12: EpicentrO!\ estimados para los dos grandes eventOR del día 17 de Diciem-
bre de 1949. Los epicentros 1 y 2 y sus respectivas elipses corresponden a
Pelayo y Wiens [1989] mientras que la ubicación del epicentro 3, fue estima-
do por Lomnitz [1970). Línea. punteada representa la ubicación de la falla
Magallanes-Fagnano.

24
Capítulo 3

Descripción y Procesamiento de Datos


Sismológicos

3. 1. Introducción
Como se describió en el capítulo I, el proyecto SEPA mantuvo una red de 5 estaciones de banda
ancha en la Patagonia chilena, desde Enero de 1997 hasta Noviembre de 1998; las estaciones FELL,
MILO, SALM, HAMB estuvieron localizadas en el sector continental de la región de Magallanes
mientras que la estación VTDF estuvo ubicada al sur de la isla Tierra del Fuego.
Los datos recolectados por esta red son la materia prima para el desarrollo de esta tesis y
particularmente en la parte correspondiente a la microsismicidad local.
Los equipos instalados se componen de un sensor de velocidad de movimiento del suelo de tres
componentes Streckeisen, modelo STS-2 de banda ancha, con respuesta plana entre 0.01 y 70 Hz
(Figura 3.1), una estación Reftek de control y registro, programable con digitalizador de 24 bits,
un receptor GPS para establecer el tiempo universal de referencia, paneles solares y j o conexión a
red/generador y un disco duro SCSI de 1.2, 2, o 4 Gbytes para almacenamiento de la información.

3.1.1. Operación de campo

La red temporal de estaciones de banda ancha permaneció operativa durante los meses de Enero
de 1997 hasta Noviembre del año 1998, con un sistema de chequeo semestral que consistió principal-
mente en mantener las unidades de poder (baterías de 12 volt) , drenar las posibles acumulaciones
de agua cerca del sensor, recambiar los discos duros de 1.2 gb y 4.0 gb los que registraban a una
razón de 24 muestras por segundo para las estaciones FELL, MILO, HAMB, VTDF y a una razón
de 40 muestras por segundo para la estación SALM.

25
3. Descripción y Procesamiento de Datos Sismológicos 3.1 Introducción

Por ser estaciones ubicadas en lugares remotos, existió una alta probabilidad de fallas de distinto
origen y que sólo podrían ser reparados en las visitas de mantenimiento semestral; es así como la
figura 3.2 muestra en forma aproximada el tiempo que las unidades estuvieron registrando en forma
continua. En general se puede apreciar que en su conjunto, las cinco estaciones trabajaron en forma
ininterrumpida la mayor parte del tiempo y sólo un pequeño porcentaje del tiempo total estando
apagadas.

RESPUESTA VELOCIDAD SENSOR STS-2

'

S
"':::J
Q)
1
o.
"'
Q) V
0:: 10
'

tOo
10"" to "1 10 o 10 ' 10'

Frecuencia (Hz)

Figura 3.1: Respuesta del sensor STS-2 usado en las estaciones sismológicas del proyecto SEPA.

Resulta fundamental para la labor de recolección de datos contar con un sistema de registro
continuo del movimiento del suelo, pues en esta zona se cuenta con la presencia de agentes externos
tales como el fuerte viento, la lluvia y cercanía de vías urbanas que, unido a un bajo nivel de
sismicidad, dificultan la adquisición de registros sismológicos vía "trigger". Si bien el ancho de banda
nos permite observar eventos locales como eventos registrados a distancias telesísmicas, este trabajo
se orienta en los aspectos locales de la sismicidad de la región de Magallanes y sólo como casos
aislados aquellos eventos lejanos que pudieron ser registrados por las estaciones de la Patagonia.
Como ejemplo de lo anterior se muestran los sismogramas originales y filtrados para dos eventos:
uno local {Figura 3.3) y otro registrado en el centro de Chile {Figura 3.5). En ellos es fácil ver que
sin la presencia de filtros, la información registrada no es procesable, ya que es imposible distinguir

26
3. Descripción y Procesamiento de Datos Sismológicos 3.1 Introducción

1997 1998

FELL
MILO

-
SALM
HAMB
VTDF

146 292 146 292


Día Juliano

Figura 3.2: Cronograma de funcionamiento de las distintas estaciones sismólogicas


durante los años 1997 y 1998.

las distintas fases presentes en los sismogramas.


Para el evento local, sólo en la estación MILO es posible ver la forma de onda y reconocer
algunas fases , lo que resulta imposible para las otras dos estaciones (HAMB, SALM). Esta situación
es reparable mediante el uso adecuado de filtros, que por lo general para los eventos locales (Figura
3.4) consitió en filtros pasa alto con frecuencia de corte de 3 Hz, siendo de esta forma identificable
las distintas fases de llegadas para las distintas estaciones. Una situación similar ocurre para los
eventos telesísmicos, donde las fases no son identificables sin el uso de un filtro apropiado, en este
caso un filtro eficiente es aquel que deja pasar las bajas frecuencias (pasa bajo), teniendo como
frecuencias de corte 0.015 Hz y 0.05 Hz {Figura 3.6) identificando, de esta forma las fases de llegada
para el evento telesísmico que significó un terremoto con daños considerables para los alrededores
de Punitaqui, IV Región.

27
3. Descripción y Procesamiento de Datos Sismológicos 3.1 Introducción

1997:017:08:32, No fdle<

¡:

S(C from 08-32

Figura 3.3: Sismogramas sin filtro de las estaciones MILO, HAMB, SALM para un evento
local con sus respectivas componentes vertical (Z}, norte-sur (N-S) y este-
oeste (E-W).

1997:017·08:32 HP 3 Hz

HAMBZ
_HAMB_l\J-S
HAMB E-W
""r ,..
¡¡...

MILOZ
..IJ.. i.lli..
r· l"'r· MILO N-S
j""'
MILO E-W
....
1
11"
SALMZ
I.L
l'f ~
~
¡
~ . . '1". ~

S(C from
..
08:32
.. .. ..
SALM E-W
-
Figura 3.4: Sismogramas con filtro pasa alto con frecuencia de corte de 3 Hz para el
mismo evento anterior.

28
3. Descripción y Procesamiento de Datos Sismológicos 3.1 Introducción

. . .8, < Octt5 >, Est Potoqonio, No rater


PUN11AOO 1997:288:01:03-30
1 -,... Mil L
.....
..,..
.... -............
...-
Mil N:ii

! ...... ./ '\....,. M_I~

!
! .,.
.d.. .....
.,. '"1 t-l:l N:>

! "T J±J .l:W_;

! ..,..Jl ~
.JiL dio.
!
l .....
.... ~, '' ~

! ''11'
. ...l:1l:illll.t

,,.
' T.,.
.,. -'lillE.L __¡;
....
. ,.
.Jt- N::.r.

•• ..,. fiOI- t:W.


SEC lrom O1:00

Figura 3.5: Sismogramas sin filtro para un evento t elesísmico ocurrido en el centro de
Chile (t erremoto de P unitaqui) y registrado por todas las estaciones del
proyecto SEPA.

PUNIIAQUI, 1997:288:01:03:30.8, < Ocl.15 >, [si. Poloqonio, 8P 0.015-0.05 Hz


y .,..,,.
lllill.
Mil ..L
1 11111111
! IVIIL 11'\l:::.r.
-VV\[V' '""'' \A" ,A. ..f\.1<.
! '" •1\t- •v
,(l ....
! ·vy t- t- l

! '"A
'V
•1\11 n
yl/Vy
....tt. ....t:.W_(;
,u«&u,
! v· 1''''' -'l8LM..L
! ,.,.
.n. "'' ::, . MN::ío
.•f\11/'1, Al
! ::.AJ.M I:VIr.
.Ulllll•
! -'"llll!J' .t1BMI1_
.UIII•
! V1'111111' ' ..J:I.'::llill.~

! '/1/'Jyy¡
" '' HJ\MH t-VI.fó

! VI V .Jt-
.n1
! •t .lt_j\.j:i';

! 'fv.JivA¡, _jJJ lt....t.'l!.l:ó


SEC lrom 01:00

Figura 3.6: Sismogramas con filtro pasa bajo con frecuencias de corte de 0.015 y 0.05 Hz
para el mismo evento telesísmico anterior.

29
3. Descripción y Procesamiento de Datos Sismológicos 3.11ntroducción

3.1.2. Secuencia de trabajo con la información recolectada

La secuencia de trabajo comienza como en toda campaña sismológica con la adquisición de los
datos crudos desde las unidades de registro para posteriormente elaborar las bases de datos con la
información útil. Por tratarse de estaciones con registros continuos, la labor de depurar la información
cruda no resultó tarea fácil, debido principalmente al excesivo ruido ambiente presente en dichos
registros, situación que llevó a desechar la idea de utilizar un programa de lectura automática, que
marca los tiempos de fases de posibles eventos sísmicos. Es por esto que se procedió a examinar
en forma visual el total de la información registrada en los dos años que estuvieron trabajando
las estaciones de la Patagonia Chilena. Dicha tarea involucró un gasto de tiempo de alrededor de
medio año y la utilización de programas computacionales diseñados para mostrar en pantalla los
sismogramas y manipularlos mediante la utilización de filtros.

Nuevo modelo de velocodades


o partir de FOCMI!C

Figura 3.7: Secuencia de trabajo realizado con Jos datos del proyecto SEPA.

La figura 3. 7 muestra en forma esquemática los principales puntos que se tomaron en cuenta
para procesar la información desde principio a fin. Es así como es posible ver que luego de la etapa
de la búsqueda visual de los eventos locales se realizó en forma conjunta una localización hipocentral
preliminar utilizando un modelo de velocidad de ondas general (IASP91 [Kennett y Engdahl, 1991]}
y en forma paralela se fue confeccionando una base de datos con aquellos candidatos a eventos loca-
les. Dicha operación consitió en cortar Jos sismogramas originales en pequeños sismogramas de 200
seg de largo, donde se presente en forma completa el registro sísmico. Para tal efecto, se utilizó el
programa de análisis sísmico SAC (Goldstein, 1996], que es un programa interactivo multipropósito

30
3. Descripción y Procesamiento de Datos Sismológicos 3.1 Introducción

diseñado para estudiar datos en series de tiempo. Además incluye soportes para análisis de operacio-
nes numéricas en general, transformadas de Fourier y muchas otras alternativas que la constituyen
en una herramienta importante en la investigación sismológica.
Como resultado de este primer paso, se logró obtener un total de 300 eventos, tomando en cuenta
sólo la calidad del registro sísmico y dejando de lado los aspectos de localización hipocentral y los
errores involucrados en dichos procesos. Por lo tanto el próximo punto a solucionar era determinar
la calidad de los registros, es decir, ver cuáles eventos tienen buenas residuales en tiempo y cuáles
no, lo que significó deshechar el modelo IASP91 y determinar el mejor modelo para la región. En
este modelo se determinaron cuáles eventos tiene buena calidad en función del valor rms.
Una vez construido el modelo de velocidad y descartado aquellos eventos con mala calidad se
construyó una base de datos con 250 eventos locales, de los cuales se leyeron los tiempos de arribo
de las fases P y S. De este total, sólo dos tuvieron magnitudes mayores a 4. El primero fue el evento
ocurrido el día 26 de Abril de 1997, a las 23:46 hrs (GMT) con magnitud 4.3 mb, y el segundo fue
el que ocurrió el día 23 de Mayo de 1998, a las 20:50 hrs (GMT) cuya magnitud fue de 4.5 mb.
Tal como se muestra en las figuras 3.8 y 3.10 se puede observar que ambos eventos (indicados
con una estrella) están separados por una distancia de alrededor de 100 km.

31
3. Descripción y Procesamiento de Datos Sismológicos 3.1 Introducción

-66•

·--

Figura 3.8: Epicentro del mayor evento registrado durante el año 1997 el que fue locali-
zado en la cercanías de Seno Skyring (estrella) a una profundidad de 19 km
con magnitud 4.3 mb.

!
!
---
1998·116·2H6 ff' 3 Hz ftl mio

.,.;.
.,
....
SllÍ11 .tdf
r.•

·-·
-
~ '1'1"
...- .
!'
¡_
.... -·
¡_
•.L
..... w.·...
! ,..
........
·--·
.,- .,,,
!
r
..
!
¡-
.....
"1' -·
r.•
i

¡
r .·

Figura 3.9: Sismograma asociado al mayor ew.nto del año 1997 (26 de Abril), registrado
por las cuatro estaciones que se encontraban operativas para dicha fecha.

32
3. Descripción y Procesam.iento de Datos Sismológicos 3.1 Introducción

Figura. 3.10: Epicentro del mayor evento registrado durante el año 1998 el que fue loca-
lizado en la cercanías de Puerto Natales (estrella.) a una profundidad de 25
km. con magnitud 4.5 mb

m
-r:.·-·..... .
1998:1 47'20:50. HP J Hz. Id, mio.

i
.L
,.. .,..,..
llll..ll ,

........
¡ ...,....

i
.J.....a..
"'ll''f'"l"' -
~-~

j, .......
"""......-
=-\.t.-
JitJ. eo..•... - ·
l'P -
-·-
1
....

1

1 -
--·
~-~

1
- - - - - -
Figura 3.11: SiRmogra.ma. ru;ociado al mayor evento de 1998 (27 de Mayo) con un filtro
pasa alto de 3 Hz

33
3. Descripción y Procesamiento de Datos Sismológicos 3.2 Sismicidad local

3.2. Sismicidad local


A continuación se muestra el modelo de velocidades utilizado para la región y el modo en que
fue obtenido. De igual forma, se muestran los resultados hipocentrales obtenidos a partir del modelo
antes mencionado.

3.2.1. Elección del modelo de velocidad

El problema a resolver con nuestros datos es el de un modelo acoplado, o sea, resolver las variables
hipocentrales y la determinación del modelo de velocidad.
Para ello se utilizó un programa computacional (VELEST [Kissling, 1988]), cuyo objetivo es
minimizar los errores entre los tiempos observados y calculados dado un modelo de velocidad inicial
y de esta forma hacer las respectivas correcciones de velocidad e hipocentrales (tiempo origen,
latitud, longitud, profundidad). El primer paso fue estimar un modelo inicial de velocidad originado
en sísmica de refracción en el estrecho de Magallanes (tabla 3.1), el que fue obtenido a partir de los
trabajos de Ludwing et al. [1965].

Capa Inte rvalo Ve locidad


capa 1 0-3 3.37
capa 2 3-6 4.1
capa 3 6-16 5.87
capa 4 16-32 6.7
capa 5 32-80 7.95

Cuadro 3.1: Modelo de velocidad obtenido de experimentos de refracción sísmica en el


Estrecho de Magallanes por Ludwing et al. [1965)

El control de calidad de las sucesivas ~aciones se fijó en el valor de la raíz cuadrática media
rms, el que no debía ser mayor a 0.5.
Luego de sucesivas ~aciones se procedió a depurar la base de datos, eliminando aquellos eventos
con errores rms muy elevados y que mostraban conductas anómalas, es decir, posibles errores en
los tiempos de llegadas de las ondas P y S, lo que significó reducir el catálogo final de eventos y
trabajar sólo con los que presentaban un buen índice rms, para llegar finalmente a los modelos en
que se apreciaba cierta convergencia (Figura 3.12).
De esta última figura se observa que a partir de un modelo de 5 capas, se llega a un modelo
final de 3 capas y un semiespacio que representa básicamemte la solución al problema acoplado, es
decir, la minimización del error hipocentral. El cuadro 3.2 muestra en más detalle la solución final
del modelo de velocidades para la zona de estudio, donde es posible ver cómo las primeras dos capas

34
3. Descripción y Procesamiento de Datos Sismológicos 3.2 Sismicidad local

propuestas por Ludwing et al. [1965] convergen en una sola, al igual que las capas 3 y 4.
La solución propuesta por el programa VELEST es muy sensible a las condiciones de entrada,
o sea, el modelo inicial de velocidades, ya que para cada modelo inicial existirá un modelo de
convergencia distinto, es por esto que el resultado no es único. En nuestro caso se probaron varias
corridas con distintos modelos de entrada existiendo algunos casos en que no era posible llegar a la
convergencia por lo que se procedió a eliminar dichos modelos.

L..-
~
...
~
~
-20

-30
.-

/
Modelo Final

!
-oo


PWao.te Velocity (llrrñ)

Figura 3.12: Modelos de velocidad generados a partir de las últimas iteraciones del pro-
grama de relocalización conjunta (VELEST)

Capa Intervalo Velocidad


capa 1 0-6 4.67
capa 2 6-32 5.73
capa 3 32-80 8.43
capa 4 80-oo 8.93

Cuadro 3.2: Modelo final de velocidad, obtenido por sucesivas iteraciones de minimiza-
ción de los residuales utilizando el programa VELEST (Kissling, 1988]

35
3. Descripción y Procesamiento de Datos Sismológicos 3.2 Sismicidad local

3.2.2. R esultad os hipocent r a les

Durante los dos años que estuvo operando la red sismológica del proyecto SEPA, se logró registrar
una cantidad de 320 eventos con magnitudes inferiores a 4.5 mb, siendo en su mayoría imperceptibles
para la población de la región. Luego del proceso de selección que dependió del valor rms para cada
evento, se redujo el catálogo a 250 eventos, es decir, eventos con valores rms < 0.5.
La figura 3.13 muestra la sismicidad total registrada durante los años 1997 y 1998 por las
cinco estaciones que conformaron el proyecto SEPA. La distribución epicentral muestra una mayor
concentración de eventos para el área Norte y Centro de la región y un menor número en el sector
de la isla Tierra del Fuego.
Los principales núcleos de sismicidad se registraron en los alrededores de Seno Skyring y en el
sector de Cueva del Milodón y es justo en estos núcleos donde se registran los mayores eventos
para el periodo 1997-1998 {Figuras 3.8 y 3.10). De igual forma se tiene que las profundidades
varían desde los 80 km. de(profundid~Jhasta valores muy superficiales. En Chile no es común
encontrar profundidades eplCentrales cercanas a la superficie, y es por esta razón que se revisó con
mucho cuidado este factor, lo que se tradujo en una variación de ciertos parámetros tales como
fijar profundidades mayores y luego determinar los errores asociados, llegando a concluir que dichas
profundidades someras respondían bien a los modelos de entrada.
Es posible distinguir sólo dos eventos ubicados costa afuera, específicamente al Oeste del Estrecho
de Magallanes, esto indicaría la presencia de actividad, producto de la subducción oblicua de la placa
Antártica bajo la placa de Scotia. Resulta raro a su vez la presencia de actividad en sectores tan
orientales como en los alrededores de la estación FELL donde supuestamente no existe un tectonismo
activo capaz de originar dichos eventos y más bien sólo quedan resabios de un pasado volcánico
activo.
La estaciones de la Red Mundial de Estaciones Sismológicas lograron registrar sólo los dos
grandes eventos del periodo, dejando de manifiesto que para los ojos del resto del mundo es una
zona con baja sismicidad. Lo anterior viene a justificar la realización del proyecto SEPA y cambiar
de esta forma nuestra visión del extremo sur de Chile.

36
3. Descripción y Procesamiento de Dat08 SismológiCOB 3.2 Sismicidad local

-72" -700

-500

LEGEND

• SEPA Statlon

-5ZJ -52"

ocfb.(\
\(\G
,.G
-540 -540



-7ff' -74° -7ZJ -700 -6SO -660

Figura 3.13: Sismicidad total de los años 1997 y 1998 obtenidas a partir del experimento
SEPA. La sismicidad se concentra principalmente en el norte y centro de la
región y las magnitudes no superan los 4.5 mb (Thbla 3.3)

37
3. Descripción y Procesamiento de Datos Sismológicos 3.2 Sismicidad local

Cuadro 3.3: Archivos final de localizaciones sísmicas del proyecto SEPA,


años 1997-1998.

Año Mes Día Tiempo Origen Latitud Longitud Profundidad RMS


97 02 11 063329.90 53.18568 73.1663W 12.19 0.546
97 02 11 063634.30 52.31178 72.2464W 0.19 0.462
97 02 12 065419.95 53.62038 72.4286W 2.92 0.120
97 02 12 110624.35 52.45528 71.6435W 16.38 0.263
97 02 13 072144.25 53.63118 72.6939W 2.98 0.276
97 02 14 002131.74 53.32518 72.5734W 23.17 0.435
97 02 15 175038.66 50.67168 72.6968W 12.02 0.663
97 02 15 201325.90 52.72718 72.5919W 2.01 0.493
97 02 16 042006.69 52.65828 72.2969W 24.26 0.132
97 02 18 084638.22 50.76598 72.3994W 44.59 0.134
97 02 18 152519.12 53.12658 71.3608W 46.02 0.357
97 02 20 094951.51 54.21638 71.7285W 13.72 0.318
97 02 21 025203.98 52.45068 71.6688W 17.80 0.046
97 02 23 085803.55 52.65528 72.0020W 15.54 0.261
97 02 28 080009.23 54.07648 73.2635W 99.90 0.541
97 03 01 105343.77 51.87268 74.2109W 61.93 0.411
97 03 02 135247.94 53.34178 71.6583W 13.65 0.265
97 03 09 030423.37 52.31668 72.1117W 15.20 0.295
97 03 09 032728.18 52.35428 71.9427W 23.99 0.076
97 03 10 222816.32 52.51788 71.8817W 23.10 0.151
97 03 12 225210.12 54.30308 74.2742W 65.27 0.528
97 03 13 043908.17 53.07898 72.7170W 12.16 0.249
97 03 14 140502.30 52.37618 72.2218W 22.24 0.177
97 03 14 152240.39 52.95708 71 .1734W 18.35 0.233
97 03 15 165342.65 52.59078 71.7043W 25.77 0.190
97 03 16 071814.32 51.25208 73.4311W 50.08 0.038
97 03 18 081729.23 51.98138 74.7743W 49.01 0.036
97 04 07 005340.04 53.15248 73.4849W 2.92 0.141
97 04 13 231148.49 51.88588 71.0601W 32.43 0.011
97 04 19 055903.93 52.97428 73.3819W 2.79 0.243
97 04 26 234652.53 52.49988 71.6278W 19.35 0.121
97 04 26 234912.51 51.68748 71.3045W 106.94 0.262
97 04 27 003311 .54 51.79478 71.2734W 75.44 0.423
97 04 27 005758.42 51.64908 71.1570W 99.89 0.078
97 04 27 035640.07 51.64588 71.2911W 112.55 0.272
97 04 28 071027.68 52.39188 71.6452W 2.53 0.085
97 04 29 154712.67 52.53418 71.7019W 22.29 0.175
97 04 30 032500.66 52.84398 74.1208W 32.68 0.298
97 05 01 061809.53 52.58418 71.6725W 25.99 0.147

38
3. Descripción y Procesamiento de Datos Sismológicos 3.2 Sismicidad local

Cuadro 3.3: Continuación

Año Mes Dia Tiempo Origen Latitud Longitud Profundidad RMS


97 05 01 064731.14 52.31628 71.9700W 2.91 0.423
97 05 04 083707.21 52.56968 71.7076W 24.90 0.151
97 05 04 084716.05 55.86248 71.5571W 25.37 0.313
97 05 05 202723.97 53.0507S 73.1133W 10.58 0.142
97 05 06 173438.49 53.4966S 71.9793W 2.81 0.301
97 05 08 221229.88 52.7591S 72.9052W 16.60 0.027
97 05 09 020051.92 52.5134S 71.8814W 23.37 0.155
97 05 09 093119.55 51.2513S 72.6769W 43.14 0.486
97 05 10 135312.91 51.24648 72.6785W 42.83 0.485
97 05 10 180240.14 52.38468 71.3780W 1.98 0.792
97 05 10 182910.25 52.6193S 72.1723W 26.38 0.131
97 05 10 193429.24 52.4710S 71.6859W 20.69 0.148
97 05 10 195702.24 52.4515S 71.6231W 17.06 0.203
97 05 11 095800.66 52.43688 71.5835W 14.91 0.170
97 05 11 103921.90 52.4435S 71.6092W 16.55 0.126
97 05 11 111921.78 52.43818 71.5695W 13.70 0.250
97 05 11 201005.53 52.89868 72.3723W 14.21 0.198
97 05 12 041144.85 52.50698 71.6499W 21.13 0.037
97 05 12 045505.93 52.74058 72.3050W 24.26 0.101
97 05 12 231455.14 52.4555S 71.6373W 18.17 0.158
97 05 13 061442.11 52.63678 72.2044W 25.67 0.163
97 05 13 100744.55 52.44268 71.6587W 15.75 0.163
97 05 13 160954.23 52.44668 71.5925W 15.44 0.189
97 05 14 120204.14 52.65108 72.6870W 19.78 0.445
97 05 15 044624.98 52.45488 71.6791W 18.19 0.163
97 05 15 070209.75 52.47148 71.6997W 17.61 0.119
97 05 16 032426.49 52.58758 71.6149W 3.00 0.243
97 05 17 043008.53 52.48208 71.7058W 21 .45 0.242
97 05 17 214640.23 52.44988 71.5844W 15.40 0.206
97 05 17 220025.04 52.45648 71.6176W 17.63 0.209
97 05 17 220229.21 52.44878 71.6216W 16.84 0.114
97 05 17 220333.93 52.41158 71.7038W 15.07 0.170
97 05 17 224445.14 52.46338 71.6445W 18.24 0.131
97 05 17 224643.78 52.44488 71.5888W 16.38 0.153
97 05 18 005443.79 52.46338 71.8191W 22.03 0.454
97 05 18 050816.25 52.46408 71.6550W 19.41 0.146
97 05 18 050734.91 52.45298 71.6081W 16.20 0.150
97 05 18 051406.21 52.46428 71.6574W 19.10 0.098
97 05 18 055832.28 52.97538 73.2017W 14.56 0.064
97 05 21 042952.06 52.45988 71.6654W 83.43 0.167

39
3. Descripción y Procesamiento de Datos Sismológicos 3.2 Sismicidad local

Cuadro 3.3: Continuación

Año Mes Dia Tiempo Origen Latitud Longitud Profundidad RMS


97 05 23 061948.72 52.96078 70.7233W 33.72 0.253
97 05 26 193639.92 52.42908 71.6553W 87.38 0.134
97 05 30 072017.13 52.44455 71.5947W 19.53 0.075
97 06 07 163004.76 52.60365 73.0251W 20.54 0.226
97 06 10 032507.09 52.46755 71.6999W 27.53 0.079
97 06 11 093300.01 52.02725 71.3910W 128.44 0.092
97 06 17 071724.18 53.42625 72.5977W 32.04 0.184
97 07 04 095558.49 51.10928 70.3850W 0.44 0.020
97 07 09 022733.03 52.49208 71.7173W 19.92 0.262
97 08 06 034744.25 53.65865 71.1225W 24.11 0.237
97 09 10 041843.51 53.07155 71.7216W 21.03 0.081
97 09 23 224255.36 51.67305 73.9221W 107.65 0.073
97 09 26 095606.48 52.68958 71.3636W 15.86 0.275
97 10 04 105850.36 51.00198 71.1796W 21.58 0.014
97 10 10 100955.34 51.85588 74.4976W 55.89 0.532
97 10 24 202133.09 51.27148 72.5219W 23.30 0.273
97 10 26 075309.45 50.49328 72.5550W 28.72 0.175
97 11 02 182316.92 50.89478 72.6343W 79.09 0.341
97 11 04 063225.21 52.98875 71.6665W 19.74 0.199
97 12 08 052721.07 53.14048 72.6092W 7.71 0.424
97 12 10 154857.91 52.73838 71.8861W 18.09 0.154
97 12 12 055345.32 51.22588 73.2663W 25.09 0.358
97 12 12 084423.54 52.43288 71.6362W 14.72 0.210
97 12 16 093505.07 53.31028 72.0164W 15.79 0.191
97 12 17 051556.00 53.13688 72.5546W 10.09 0.238
97 12 25 212314.99 51.09995 71.0189W 21.68 0.012
97 12 28 030928.34 52.94305 73.4136W 3.75 0.740
97 12 30 034543.38 51.49155 72.3451W 20.77 0.331
97 12 30 074501.48 51.50398 72.3688W 21.11 0.468
97 12 30 074914.73 51.54378 72.3183W 17.35 0.258
97 12 30 163106.73 54.55375 73.9882W 32.49 0.504
97 12 31 034909.37 51.57985 70.6232W 32.15 0.123
98 01 05 092008.37 52.85625 72.0526W 23.99 0.181
98 01 12 103502.94 53.82045 71.0403W 17.32 0.501
98 01 12 104458.32 53.85375 71.0357W 14.66 0.076
98 01 15 005108.58 52.48785 71.6573W 20.22 0.095
98 01 17 033507.09 52.45755 72.7373W 20.69 0.082
98 01 21 150455.44 51.95575 74.0059W 10.78 0.746
98 01 21 231200.60 52.23085 73.2068W 1.17 0.200
98 01 22 232926.01 52.73555 70.2264W 32.17 0.254

40
3. Descripción y Procesamiento de Datos Sismológicos 3.2 Sismicidad local

Cuadro 3.3: Continuación

Año Mes Dia Tiempo Origen Latitud Longitud Profundidad RMS


98 01 30 130832.87 54.97038 66.5484W 60.70 0.349
98 02 09 021137.75 53.51098 71.9353W 2.29 0.314
98 02 13 100719.81 53.24738 73.6127W 0.80 0.172
98 02 16 190405.84 52.05938 70.0197W 2.33 0.000
98 02 17 093457.16 51.73378 74.0315W 88.15 0.532
98 02 17 130320.82 52.74698 72.3065W 21.03 0.078
98 02 21 200450.38 51.26968 73.0358W 135.45 0.033
98 03 02 171122.00 53.76758 71.7085W 30.57 0.142
98 03 04 030307.80 53.08538 72.7384W 10.03 0.238
98 03 06 031031.61 53.13718 72.6598W 2.86 0.497
98 03 09 012858.84 54.47958 69.2750W 2.78 0.698
98 03 18 235211.96 51.63528 73.6786W 35.42 0.174
98 03 19 212310.19 52.50538 73.3717W 7.16 0.606
98 03 21 081629.36 51.35208 73.5958W 44.82 0.087
98 03 29 112848.45 50.95708 73.5569W 1.52 0.143
98 03 29 161007.10 50.95538 73.4714W 2.24 0.493
98 03 29 205045.78 50.98398 73.4367W 2.05 0.246
98 03 29 235908.98 53.41478 71.8629W 2.90 0.618
98 03 31 214702.40 50.93628 73.3126W 7.62 0.235
98 04 01 051010.54 50.90668 73.3262W 0.30 0.306
98 04 03 093612.20 49.67408 74.0209W 37.24 0.097
98 04 03 153359.41 50.96028 73.3866W 0.15 0.676
98 04 03 175609.16 51.49228 72.3170W 18.41 0.122
98 04 03 183137.02 51.22248 73.7029W 0.34 0.444
98 04 03 183135.51 51.02588 73.6819W 1.88 0.250
98 04 07 030512.32 52.33068 69.5987W 0.19 0.023
98 04 07 045350.01 49.75958 74.3672W 38.94 0.139
98 04 15 013733.08 52.63718 72.2274W 22.67 0.051
98 04 15 165544.23 52.82168 73.5938W 0.09 0.451
98 04 15 224429.37 51.51578 72.4005W 22.16 0.500
98 04 18 234409.26 53.06568 72.2220W 19.23 0.380
98 04 20 163702.06 50.94268 73.3850W 0.73 0.095
98 04 22 035056.47 50.79578 73.1243W 6.18 0.134
98 04 22 042910.62 50.88318 73.4014W 3.51 0.702
98 04 22 070311.30 51.02038 73.5695W 0.48 0.282
98 04 23 100126.65 51.58108 72.2885W 5.27 0.031
98 04 25 014813.21 53.56368 72.4994W 2.93 0.199
98 04 28 221429.78 53.27868 72.1631W 0.90 0.216
98 05 01 040344.43 52.24538 72.2889W 16.41 0.181
98 05 03 201917.67 52.34968 73.5986W 8.28 0.027

41
3. Descripción y Procesamiento de Datos Sismológicos 3.2 Sismicidad local

Cuadro 3.3: Continuación

Año Mes Día Tiempo Origen Latitud Longitud Profundidad RMS


98 05 04 073234.74 50.73588 73.4066W 1.97 0.322
98 05 04 081257.71 50.98878 73.5369W 0.22 0.244
98 05 06 064630.03 50.87628 73.4568W 1.20 0.101
98 05 06 133740.82 51.01418 73.5363W 0.14 0.317
98 05 07 070008.31 50.90088 73.4123W 1.19 0.090
98 05 13 150600.44 50.78348 72.1868W 13.28 0.218
98 05 14 231859.48 50.79248 72.1969W 12.74 0.217
98 05 15 194559.39 50.96118 73.5220W 0.88 0.313
98 05 15 200327.00 51.03548 73.4803W 0.73 0.177
98 05 16 053021.89 50.96458 73.4527W 0.19 0.140
98 05 17 095754.89 50.83068 72.2863W 15.77 0.056
98 05 21 093933.27 52.51418 72.3321W 15.18 0.234
98 05 27 205041.33 51.46518 72.5442W 25.49 0.154
98 05 27 205624.52 51.50938 72.4287W 22.13 0.045
98 05 28 065412.66 51.52948 72.2885W 16.33 0.087
98 05 28 142119.58 51.48978 72.4594W 23.53 0.054
98 05 29 040025.07 51.54028 72.3728W 20.13 0.007
98 05 29 103506.85 52.33598 71.9491W 44.44 0.093
98 05 29 114429.61 51.48718 72.4658W 23.95 0.050
98 05 31 063255.53 51.60068 72.4736W 19.15 0.284
98 05 31 074433.70 51.54858 72.2929W 15.72 0.145
98 06 06 175714.82 53.36428 71.4723W 9.72 0.782
98 06 07 152945.39 52.84478 73.5278W 5.50 0.222
98 06 08 170731.44 50.53648 72.7879W 52.14 0.354
98 06 11 083435.33 50.46825 72.6159W 26.39 0.340
98 06 16 090904.25 52.82108 73.4681W 5.52 0.089
98 06 18 011001.62 51.54518 72.4643W 21.88 0.160
98 06 20 044334.54 51.50268 72.2722W 16.98 0.077
98 06 20 081404.81 51.53078 72.3360W 18.35 0.180
98 06 21 122514.17 52.49008 71.6997W 23.70 0.013
98 06 25 064622.71 51.56578 72.4686W 21.73 0.275
98 06 25 103239.70 51.50308 72.3790W 20.91 0.148
98 06 27 065623.53 53.04688 73.2785W 8.35 0.264
98 06 29 014802.36 51.55748 72.4580W 21.57 0.205
98 06 29 194424.13 53.13588 71.8725W 16.90 0.350
98 07 02 080246.71 53.03808 71.6364W 8.04 0.319
98 07 03 200508.57 52.81678 73.5912W 0.68 0.329
98 07 05 201257.86 51.46668 72.4649W 25.04 0.187
98 07 06 192224.31 51.48515 72.4370W 23.55 0.225
98 07 07 203740.36 51.48818 72.3903W 22.35 0.255

42
3. Descripción y Procesamiento de Datos Sismológicos 3.2 Sismicidad local

Cuadro 3.3: Continuación

Año Mes Dia Tiempo Origen Latitud Longitud Profundidad RMS


98 07 09 173618.44 52.52488 71.8589W 25.65 0.613
98 07 14 191600.62 52.85528 73.4613W 7.37 0.192
98 07 16 124233.53 52.45538 71.6221W 17.55 0.204
98 07 20 113148.24 52.44388 72.2700W 27.23 0.035
98 07 22 102405.16 51.54908 72.4569W 21.66 0.159
98 07 28 112600.88 51.99348 72.4928W 16.55 0.260
98 07 31 014552.43 53.06348 72.8206W 17.94 0.121
98 07 31 220059.50 51.43638 73.1990W 129.84 0.047
98 07 31 231410.72 51.94228 73.8908W 2.21 0.296
98 08 01 043608.16 51.53118 72.3694W 20.09 0.006
98 08 02 232810.74 51 .84638 74.3077W 32.10 0.605
98 08 09 005126.05 51.94348 70.1707W 13.11 0.008
98 08 10 072317.11 52.43808 71.6644W 17.08 0.169
98 08 13 054817.94 53.27258 71.7891W 1.48 0.260
98 08 14 064822.53 52.48648 71.6794W 18.71 0.148
98 08 14 065034.09 52.60978 72.1925W 25 .04 0.194
98 08 27 154408.49 52.87138 73.8702W 0.31 0.549
98 08 29 215920.22 52.72758 72.3611W 21.80 0.100
98 09 01 013425.03 53.16818 71.7341W 16.58 0.116
98 09 03 224627.00 54.77318 71.3161W 70.72 0.355
98 09 06 093326.58 52.03648 71.4117W 32.73 0.002
98 09 16 102850.03 53.21045 72.4143W 19.38 0.121
98 09 19 002216.39 53.17678 73.1203W 6.95 0.054
98 09 22 205957.45 52.47295 71.4083W 2.51 0.212
98 09 23 113943.49 52.72655 72.3170W 25.76 0.053
98 09 24 225731.44 51.47078 72.5394W 27.08 0.477
98 09 24 233641.00 51.54305 72.4820W 22.06 0.174
98 09 25 030212.75 51.50358 72.2208W 12.91 0.165
98 09 25 040503.06 51.45585 72.4529W 24.77 0.457
98 09 25 055041.64 51.57205 72.3546W 17.75 0.305
98 09 26 025300.70 53.93438 72.1607W 17.91 0.621
98 09 27 123446.21 51.68038 70.3333W 17.16 0.118
98 09 27 202856.60 52.04908 69.9780W 0.29 0.000
98 09 29 110009.84 50.60718 73.2328W 56.46 0.110
98 09 30 151340.89 50.96568 73.2089W 110.22 0.795
98 10 02 121443.63 51.51198 72.3373W 19.37 0.202
98 10 03 095548.18 53.07408 72.8654W 13.24 0.406
98 10 05 224306.72 50.70568 71.5046W 25.98 0.181
98 10 09 102524.42 50.99298 72.6021W 27.29 0.596
98 10 09 160740.20 51.64048 70.7059W 30.17 0.011

43
3. Descripción y Procesamiento de Datos Sismológicos 3.2 Sismicidad local

Cuadro 3.3: Continuación

Año Mes Di a Tiempo O rigen Latitud Longitud Profundidad RMS


98 10 09 164634.23 51.14928 73.1819W 36.05 0.101
98 10 11 071642.32 51.53518 72.3374W 17.50 0.006
98 10 11 203804.85 51.08028 73.3285W 81.18 0.075
98 10 11 224114.92 51.08768 73.6830W 33.64 0.387
98 10 15 215406.74 50.80998 72.2765W 14.65 0.118
98 10 16 081608.67 50.79068 72.3846W 29.03 0.056
98 10 21 120151.87 50.79068 72.3846W 29.03 0.056
98 10 27 022105.28 53.48228 71.9276W 3.10 0.263
98 10 29 131609.46 55.50858 71.8008W 60.81 0.568
98 10 31 052156.07 51.51748 72.3854W 21.53 0.369
98 11 01 061147.23 52.28468 73.4801W 0.44 0.277
98 11 02 053529.30 52.64488 72.5052W 22.40 0.783
98 11 02 070233.14 51.66848 72.7358W 3.06 0.654
98 11 04 060843.44 54.14998 72.2176W 11.45 0.164
98 11 05 054046.29 51.49818 72.3670W 20.28 0.098
98 11 05 210324.56 51.34938 72.3324W 23.52 0.126
98 11 07 184952.60 53.81588 70.5781W 4.12 0.392
98 11 08 011944.03 52.45808 73.8035W 2.78 0.360
98 11 08 131219.92 51.72958 74.0799W 76.24 0.575
98 11 09 051249.23 52.64128 71.9346W 12.96 0.329
98 11 10 085331.40 54.44488 70.1980W 8.58 0.702
98 11 11 144345.91 50.75728 74.8726W 0.30 0.127
98 11 12 112030.58 50.72628 72.4480W 60.33 0.393
98 11 13 112318.31 51.38688 73.5529W 36.24 0.145

44
3. Descripción y Procesamiento de Datos Sismológicos 3.2 Sismicidad local

3.2.3. Estructura Cortical

A partir de la solución del modelo de velocidad generado por el programa VELEST se puede
inferir que para profundidades mayores de 32 km se estaría en presencia del manto. En virtud de
las variaciones de velocidad, se puede concluir que el espesor de la corteza es de aproximadamente
32 km, consistente con los valores de espesor de corteza asignados para la Patagonia, donde se le
asignan valores de alrededor de 30 km. (Figura 3.14).
La determinación con más exactitud del espesor de la corteza nos permite inferir acerca del
comportamiento de las ondas sísmicas, al ser el manto y la corteza dos unidades estructurales muy
bien definidas tanto por composición como por el contraste de las velocidades de propagación de las
ondas P y S. Es claro asumir además que el espesor de la corteza para el Norte de Chile es mayor
que en el lado Sur, ratificado principalmete por el espesor de la raíz de la cordillera de los Andes.

Figura 3.14: Espesores de corteza a nivel mundial, asumiendo espesores de 30-40 km.
para el Sur de Sudamérica (http:/ / mahi.ucsd.edu/Gabi/rem.dirfcrust/crust2.html).

45
3. Descripción y Procesamiento de Datos Sismológicos 3.2 Sismicidad local

AnaHzando la sismicidad superficial, es posible ver en forma más gráfica la interacción entre cor-
teza y manto. La figura 3.15 muestra en forma esquemática cómo se comportan los rayos originados
en una fuente sísmica en la zona de la corteza que, por tener una velocidad de propagación menor
a la zona del manto, se produce en la estación receptora un primer arribo (onda Pn) para después
llegar la onda primaria P.
Por tal motivo se diseñaron gráficos de distancia epicentral v /s tiempo para la sismicidad regis-
trada a profundidades menores de 3 km, con el objetivo de poder observar de mejor forma la llegada
de las distintas fases (P, Pn, S).
La figura 3.16 muestra los sismogramas registrados en las componentes verticales para aquelJos
eventos con profundidad menor a 3 km. En dicha figura es posible ver en forma clara los tiempos
de arribo para la onda primaria mientras que en forma más o menos difusa las ondas generadas por
la refracción del manto.
La figura 3.17 muestra los sismogramas registrados por las componentes horizontales para la
misma profundidad anterior, observándose de igual forma el carácter refractor para la onda S.

Pn

Corteza ____.

Manto ----.

Figura 3.15: Modelo de dos capas donde se muestra la generación de las fases refractadas
(Pn, Sn) así como las fases directas (P , S).

46
3. Descripción y Procesamiento de Datos Sismológicos 3.2 Sismicidad local

100 200 400


(Km)
300

50

(Seg)

Figura 3.16: Distancia epicentral vjs tiempo (componente vertical) para eventos locali-
zados a profundidades menores de 10 km.

47
3. Descripción y Procesamiento de Datos Sismológicos 3.2 Sismicidad local

100 200 300 400


(Km)

10

(Seg)

Figura 3.17: Distancia epicentral v/s tiempo (componentes horizontales) para eventos
localizados a profundidades menores de 10 km.

48
Capítulo 4

Resultados:Análisis y Discusión

4.1. Sismicidad local durante el periodo del Proyecto SEPA


La sismicidad registrada durante los años 1997 y 1998 por los instrumentos del proyecto SEPA
permiten ver por primera vez cómo es la distribución de microsismicidad para el extremo Sur de
Sudamérica, situación que no era reflejado en las estaciones de la Red Mundial. Dicha distribución
epicentral nos muestra que la mayor parte de ésta, o en otras palabras las zonas más sísmicas de
la región, se encuentran en la zona de la faja corrida y plegada (Figura 2.3) sobresaliendo la alta
concentración en la zona de Puerto Natales y Seno Skyring.
Lo anterior nos da una percepción errada de la existencia de la Falla de Magallanes ya que no es
posible ver en forma clara un lineamiento en torno a ésta y de esa manera fortalecer las hipótesis que
plantean a esta falla como uno de los principales agentes tectónicos, siendo causante de los mayores
eventos para la región (Terremoto de 1949, M=7.5).
La baja sismicidad asociada a la falla de Magallanes puede ser explicado en términos de que
dicha falla sólo es activada por grandes terremotos, dejando disipar la energía acumulada a través
de pequeños eventos en fallas laterales a la falla principal (Figura 4.1). Dicha figura muestra ciertos
lineamientos en forma oblicua a la falla principal (lineas azules) que, en términos estructurales,
dichas alineaciones podrían ser explicadas como un sistema de fallas Riedel Shear (R shear) ya que
por definición, estas fallas se desarrollan en ángulos pequeños (10°- 20°) de la falla principal.
Fortaleciendo dicha hipótesis se debe recordar que las fracturas R shear son sintéticas a la falla
principal, es decir, son subparalelas y tiene el mismo sentido de cizalle que la falla principal.
Vera et al. [1997] postulan otro sistema de fallas para explicar la liberación de energía de la falla
principal. Ellos propusieron un sistema como el descrito en la figura 4.1 (líneas verdes) . Si bien el
sistema es coherente con los esfuerzos involucrados en el área, la sismicidad registradda no avala

49
4. Resultados:Análisis y Discusión 4.1 Sismicidad local durante el periodo del Proyecto SEPA

-7'P -70° -680 -640

-SOO

-5'P -5'P

-70° -6SO

Figura 4.1: Lineamientos de aismicidad: Líneaa azulea repreaentan los principalea linea-
mientos generados por SEPA. Líneas verdea, muestran el modelo de distri-
bución de fallas generados por Vera et al . (1997]

50
4. Resultados:Análisis y Discusión 4.2 Mecanismos de foco

dicho modelo o quizás una registración más prolongada en el tiempo podría dar más valor a dicho
modelo.
Por otro lado, un hecho importante de destacar es la presencia de una sismicidad difusa localizada
al oriente de la zona de estudio, Provincia de Santa Cruz, Argentina, siendo que estas zonas se
encuentran alejadas de las principales estructuras (placas) que serían las responsables de la actividad
sísmica de la región. Revisando la actividad exploratoria del sur de Argentina, queda descartado que
dichos eventos tengan UD origen en la exploración de hidrocarburos, aunque es probable que dicha
actividad sea asociada a la presencia de fallas geológicas superficiales.
Los núcleos más importantes de sismicidad se encuentran en los alrededores de Puerto Natales y
es en ese lugar donde se produjeron los mayores eventos registrados durante los dos años que duró el
proyecto.

4.2. Mecanismos de foco

La sismicidad local aportada por el proyecto SEPA es principalmente de muy baja energía
con magnitudes menores a 4.5 mb. Por tratarse de un pequeño número de estaciones resulta casi
imposible determinar las soluciones de los mecanismos focales a través de las polaridades de la
primera llegada de la onda P. Es por esto que se utilizó un metodo de inversión de las formas de
ondas P, SV y SH para los dos mayores eventos registrados por la red de estaciones del proyecto
SEPA.
La principal característica de esta inversión es que a partir de la variación de los parámetros
de velocidad y/o párametros de la fuente sísmica se pueden generar sismogramas sintéticos los que
pueden ser comparados con los datos reales.
Las figuras 4.2-4.7 muestran los resultados de la inversión para el evento del 27 de Mayo de
1998 cuya magnitud fue 4.5 mb. Dicha inversión muestra un mecanismo de falla tipo inverso con
- t_:na componente de rumbo. De igual forma, se generó el mecanismo de foco para el evento del 17
de Abril de 1997 (4.3 mb) , mostrando un mecanismo de falla inversa con UD rumbo de dirección
Norte (figura 4.8). Este último mecanismo estaría indicando que en el área cercana a Seno Skiring,
la subducción de la Placa Antártica bajo la Sudamericana sería la causante de la actividad sísmica
del área.

51
4. Resultados:Análisis y Discusión 4.2 Mecanismos de foco

íELL BBR
íREOUENCY LIMITS O. 0.200 TIME PTS 128 DT 2.
NUMBER OF LPYERS IN MODEL 22 TOTPL THICKNESS 449
OEPTH 7. RPNCE J. ~ZIMUTH 72.5
STRIKE 5 OIP 55 SLIP 90 DURPTION 2.
-!ME SHIFT -J. MOMENT J074.6J * JB NM ERROR 0.250 WEICHT J.losalmfel
toce\domod op25

Pg Sg

\
1 \
f\
1 \
1 \
1
1 \
1 1
1
1 1
1 1
1 1
1
1 1
\ 1
\

o 20 40 50 80 100
TIME (sl

Figura 4.2: Componente radial. Estación FELL. Evento 27-05-98. La línea continua re-
presenta la onda teórica mientras que la línea segmentada muestra la forma
real de la onda.

52
4. Resultados:Análisis y Discusión 4.2 Mecanismos de foco

FELL BBT
FREOUENCY LIMITS O. 0.200 TIME PTS 128 OT 2.
NUMBER OF LRYERS IN MODEL 22 TOTRL THICKNESS 449
DEPTH 7. RRNGE 1. RZIMUTH 72.5
STRJKE 5 DJP 55 SLIP 90 DURRTION 2.
TIME SHIFT 0.50 MOMENT 1074.61 **18 NM ERROR 0.169 WEIGHT 1.
toceyomod op25

,...
1 \
1 1
1
1
1 1
1 1
\
\
\
\
\
\ 1
1 1
11
\..

o 20 40 50 80 100
TIME (sl
Figura 4.3: Componente t ransversal. Estación FELL. Evento 27-05-98. La lfnea continua
representa la onda teórica mientras que la línea segmentada muestra la forma
real de la onda.

53
4. Resultados:Análisis y Discusión 4.2 Mecanismos de foco

FELL BBZ
FREOUENCY LIMI TS O. 0.200 TIME PTS 128 OT 2.
NUMBER OF L~YERS IN MOOEL 22 TOT~L THICKNESS 449
OEPTH 7. R~NGE 1. AZIMUTH 72.5
STRIKE 5 DIP 55 SLIP 90 DURATION 2.
TIME SHIFT O. MOMENT 1074.61 ~*18 NM ERROR 0 .1 01 WEIGHT 1.
1aceyorry6.d ap26

1 "' '\
h \ ./
\
\.1
1 '
1
1
1
1
1
1
\ 1
\1

o 20 40 60 80 100
TI ME (sl
Figura 4.4: Componente vertical. Estación FELL. Evento 27-05-98. La línea cont inua
representa la onda teórica mientras que la línea segmentada muestra la forma
real de la onda.

54
4. Resultados:Análisis y Discusión 4.2 Mecanismos de foco

MILO BBR
FREOUENCY LIMITS O. 0.200 TIME PTS 128 OT 2.
NUMBER OF L~YERS IN MODEL 22 TOT~L THICKNESS 449
DEPTH 7. R~NGE 1. ~ZIMUTH 322.2
STRIKE 5 DIP 55 SLIP 90 OURRTION 2.
TIME SHIFT -2. MOMENT 1074.61 **18 NM ERROR 0.264 ~ EIGHT 1.
loceyomod op26

P9 S9

'\
1
1
1
1
\
o 1
1
1 '\
7 1
1
1 1 1
\./ 1
1
1
1 1 1
1 1 1 1
'1 \ 1
./

o 20 40 60 80 100
TIME (sl
Figura 4.5: Componente radial. Estación MILO. Evento 27-05-98 . La línea continua
representa la onda teórica mientras que la línea segmentada muestra la forma
real de la onda.

55
4. Resultados:Análisis y Discusión 4.2 Mecanismos de foco

MILO BBT
FREOUENCY LIMITS O. 0.200 TIME PTS 128 OT 2.
NUMBER OF L~YERS IN MOOEL 22 TOT~L THICKNESS 449
OEPTH 7. R~NGE 1. ~Z I MUTH 322.2
STRlKE 5 DIP 55 SLIP 90 OUR~T!ON 2.
TIME SHIFT -1.50 MOMENT 1074.61 **1~ NM ERROR 0.118 ~EIGHT J.
toce y o mod op26 ,' 1

Pg Sg

1 ""

IÍ \
\ 1
7 \ 1
\/

o 60 80 100
( sl
\1

Figura 4.6: Componente transversal. Estación MILO. Evento 27-05-98. La línea continua
representa la onda teórica mientras que la línea segmentada muestra la forma
real de la onda.

56
4. Resultados:Análisis y Discusión 4.2 Mecanismos de foco

MILO BBZ
FREOUENCY LIMITS Q. 0.200 TIME PTS 128 DT 2.
NUMBER OF L~YERS IN MODEL 22 TOT~L THICKNESS 449
DEPTH 7. R~NGE l. AZIMUTH 322.2
STRIKE 5 DIP 55 SLIP 90 OUR~TION 2.
TIME SHIFT -1. MOMENT 1074,61 **18 NM ERROR 0.121 WEIGHT 1.
oceld o m b~ op26
l
1

Oo2
,...
r
1 \ 1 \ .......
/, l \ 1 \
\ 1
o \
/
\ /

~ 1
1 \_.,
'
,1 1
1
-Oo2 1
1
1
1 1
11
1..
-Oo4

o 20 40 60 80 100
TIME (sl
Figura 4.7: Componente vertical. Estación MILO. Evento 27-05-98. La línea continua
representa la onda teórica mientras que la línea segmentada muestra la forma
real de la. onda..

57
4. R.esultados:Análisis y Discusión 4.2 Mecanismos de foco

-7SO

-50° -50°

LEGEND

• SEPA Statlon

-52'1

-70°

Figura 4.8: Soluciones de mecanismos de foco para la región de Magallanes. Mecanismos


de foco color negro fueron determinados a partir de los datos del proyecto
SEPA , mientras que los de color gris, son mecanismos determinados a partir
de la red mundial (Harvard-CMT, 2002].

58
4. Resultados:Análisis y Discusión 4.3 Sismicidad asociado a volcanismo

El catálogo Harvard [Harvard-CMT, 2002) de mecanismos de foco para distancias telesísmicas,


muestra sólo tres soluciones focales para la Región de Magallanes (tabla 4.1) con párametros similares
a los determinados con los datos locales.
Fe cha M agnitud ( mb) Latit ud Longitud Profundidad St rike Dip Slip
06/06/1979 5.6 -52.93 -74.65 15 7 8 104
30/07/1992 5.4 -51.12 -72.59 15 31 25 114
31/08/1996 4.9 -53.43 -73.12 30.2 54 77 -148

Cuadro 4.1: Mecanismos de foco obtenido del catálogo Harvard-CMT [2002] para los
mayores eventos registrados desde el año 1979.

Mostrando una singular solución ubicada en los márgenes de la falla de Magallanes que muestra
un mecanismo de falla de rumbo, haciendo notar la actividad de la falla de Magallanes, la que puede
ser robustecida si se muestra la sismicidad superficial asociada a la parte Oeste del Estrecho de Ma-
gallanes, donde es posible ver lineamientos en la continuación del Seno Almirantazgo y lineamientos
secundarios en torno a la falla de Magallanes.

4.3. Sismicidad asociado a volcanismo

La téctonica de la región austral de Sudamérica difiere en forma notable de la zona Norte y


Centro-Sur de Chile. Esto se refleja por la escasa actividad volcánica generada desde el P leistoceno
basta la fecha, relacionada a la baja velocidad de subducción, es decir, 3 cm por año de la Placa
Antártica bajo la placa Sudamericana [DeMets et al., 1994]. De esta forma se puede ver (figura
4.9) que al sur del punto triple de las placas de Nazca-Antártica-Sudamericana existe un "vacío",
es decir, ausencia de centros volcánicos, siendo de esta forma, el Volcán Hudson el último volcán
asociado a la placa de Nazca mientras que el volcán Lautaro es el primero en asociarse a la placa
Antártica.
El cinturón volcánico relacionado a la acción de la P laca Antártica sobre la placa Sudamericana
comenzaría con el volcán Laut aro por el Norte y 700 km al Sur con el volcán Cook. A pesar de la
baja actividad volcánica registrada en la zona.

59
4. R esultados:Análisís y Discusión 4.3 Sísmícidad asociado a volcanísmo

Figura 4.9: Cadena volcánica de los andes australes según Gonzalez-Ferran [1995] dentro
del contexto tectónico de la región.

60
4. Resultados:Análisis y Discusión 4.4 Breve reseña del volcanismo austral

4.4. Breve reseña del volcanismo austral


A continuación se mencionan las principales características de la cadena volcánica de la zona de
estudio, junto con algunos reportes de actividad en épocas pasadas.

Volcán Lau.taro: Este volcán ubicado en pleno Campo de Hielo Sur (49°01' S, 73°33' W, 3380
m) corresponde a un complejo estrato-volcánico afectado en forma parcial por la erosión glaciar, ya
que cerca de un 90% de éste se encuentra cubierto por una gruesa capa de hielo. El cráter principal
se emplaza en la base de una antigua caldera de aproximadamente un kilómetro de diámetro.
Los primeros informes de actividad eruptiva datan de mediados del siglo XIX , sin embargo, en
el año 1959 fue la primera vez que se divisó una erupción en forma directa, gracias a los reportes
de un piloto civil que sobrevolaba la región. Dicha actividad produjo un pequeño flujo de lavas y
piroclastos. Luego en el año 1960 se vuelven a observar abundantes fumarolas en un cráter activo a
100 m bajo la cumbre.

Volcán Viedma: El volcán Viedma (49°22' S, 73°19' W, 1500 m) está ubicado a unos 45 km al sur-
sureste del Volcán Lautaro. Dicho volcán está formado por cuatro cráteres volcánicos con diámetros
que varían de 1.5 a 3 km. Actualmente sólo el cráter sur permanece activo [Gonzalez-Ferran, 1995].
Fotos aéreas tomadas el año 1946 muestran tres franjas con depósitos de piroclastos sobre el
glaciar, mostrando de esta forma la ocurrencia de erupciones explosivas.

Volcán Aguilera: Ubicado en la costa noroccidental del fiordo Peel (50°20' S, 73°45' W, 2546
m). La composición de este volcán lo forman lavas y piroclastos de composición dacítica [Gonzalez-
Ferran, 1995].

Volcán Reclu.s: Corresponde a un edificio volcánico ubicado en el centro del glaciar Amalia (50°57'
50" S, 73°35' 05" W 1000 m). El cráter presenta un diámetro de 1 km, abierto hacia el noroeste y
cubierto por hielo. La ubicación del volcán fue descubierta por primera vez por el geólogo Salvador
Harambour el año 1987, quien señaló que se trataba de un cono volcánico con actividad post glacial.

Volcán Burney: Este volcán se encuentra localizado en la parte noroeste de la Península Muñoz
Camero específicamente a 52°20' S 73°24' W y con una altura que alcanza los 1700 m. Corresponde
a un complejo volcánico, compuesto por un estrato-caldera basal [Gonzalez-Ferran, 1995]. Un 70 %
del complejo volcánico se encuentra cubierto de hielo. Sus ciclos evolutivos habrían ocurrido durante
el Pleistoceno, Holoceno y Reciente.

61
4. Resultados:Análisis y Discusión 4.4 Breve reseña del volcanismo austral

Volcan es de Pali-A ike : Ubicado en el sector limítrofe de Chile y Argentina (51 °30' a 52°30' S,
70°00' W, 282 m) , corresponde a un campo de erupciones monogenéticas que cubre un área aproxi-
mada de 3000 km 2 . En dicho campo es posible encontrar maares, flujos lávicos y conos piroclásticos,
controlados por un sistema de fallas de rumbo N 45°W [Gonzalez-Ferran, 1995]. La composición
de lavas y piroclastos corresponde a basaltos olivínicos. Las estructuras volcánicas de Pali-Aike no
presentan erosión glacial por lo que se supone una edad menor que 1 Ma. El centro eruptivo más
reciente está asociado al volcán Cerro del Diablo cuya edad es menor a 15000 años.

Volcán Cook: Se encuentra ubicado en el sector oriental de la isla Cook (54°57' S, 70° 15' W , 150
m). En términos generales se trata de un conjunto de al menos ocho conos piroclásticos en donde
uno de los mayores cráteres presentan un diámetro de 150m. Estas estructuras no presentan erosión
glaciar, por lo que se considera una edad post-glacial menor a 11000 años.

62
4. Resultados:Análisis y Discusión 4.4 Breve reseña del volcanismo austral

4.4.1. Sismicidad asociada al Volcán Reclus

La sismicidad superficial {figura 4.10) asociada con profundidades menores a 3 km muestran


cierta tendencia hacia la zona del volcán Reclus mientras que para el resto de los volcanes no sucede
lo mismo.
Los eventos asociados a la zona del volcán Reclus ocurren en un período de tiempo determinado:
Abril-Mayo de 1998 mostrando todos ellos sismogramas de similares características (figura 4.11)
donde las fases de las ondas P y S se tornan muy difusas a diferencia de otros eventos de similares
profundidades epicentrales pero localizadas fuera del área del Reclus, donde sí es posible ver en
forma clara, las fases de llegada.
A pesar de lo difuso de algunas fases, se puede ver que para dichos eventos, las cinco estaciones
fueron capaces de registrar tal actividad, inclusive estaciones tan lejana como VTDF que estaba en
la parte Sur de Tierra del Fuego a aproximadamente 500 km del Área del volcán Reclus.
En el detalle {Figura 4.10) se puede observar la concentración de eventos superficiales alrededor
del volcán Reclus, lo que lleva a conjeturar en forma más rigurosa que dicha actividad fue originado
por el volcán durante los meses de Abril y Mayo, no manifest~dose de nuevo durante el resto del
experimento SEPA.

63
4. Resultados:Análisis y Discusi6n 4.4 Breve reseña del volcan.ismo austral

-74°

-500

LEGEND

• SEPA Sútlon

-52"

-76° -74° -72" -70° -68° -66"

Figura 4.10: Sismicidad superficial ( <3km) registrado durante los años 1997 y 1998 por
la. red temporal SEPA junto a. la cadena. de volcanes de la. región. El acer-
camiento muestra. la actividad sismológica ocurrida en los alrededores del
volcán Reclus 64
4. Resultados:Análisis y Discusión 4.4 Breve reseña del volcanismo austral

1998:088:11 :28. HP 3 Hz, fell, milo, solm, vtdf

J.
"'" .tii.IJiü>.
~r .,......,...,.. -·
J.
...... .... -·
.!r.'
. ' l"T '!'[ ""
wl .
~

...,.,.
...l.tlllll ...... ~
~-·
·-.
J.
&...

l. ..,, .........
.i&J
,.. -·
J.
lJ.

J.
LM,
r"" ' -·
::.·.....·
~
u...
.,..., ... ~-·

--" r.·

.....-_.
7.•
J.
.
-~·
...IL ¡.{
.,.. .

..L. ..1. .....


-·~ .
.., !"

Figura 4.11: Sismograma típico de los eventos superficiales loca.lizado en los a.lrededores
del volcán Reclus durante el período Abril-Mayo de 1998.

65
4. Resultados:Análisis y Discusión 4.5 Perfiles sfsmicos regionales

4.5. Perfiles sísmicos regionales

La industria del petróleo ha explorado una amplia zona de la región de Magallanes con motivos
principalmente comerciales, es por eso que hasta la fecha se cuenta con miles de km de sísmica 2-D
así como un volumen considerable de sísmica 3-D tanto para la Isla Tierra del Fuego como para el
continente. Desde los años 50 que se vienen registrando campañas exploratorias principalmente de
sísmica de reflexión y algunos pocos perfiles orientados a la sísmica de refracción.
Tanto la Empresa Nacional de Petróleo (ENAP) como investigadores externos se han dedicado
a recopilar información de perfiles sísmicos y así junto a abservaciones geológicas de superficie han
podido recrear las distintas estructuras que componen la cuenca de Magallanes, es decir estructuras
que no se encuentran a profundidades mayores a 3 km y en casos excepcionales hasta los 6 km de
profundidad.
Autores como Diraison et al. [2000] han propuesto tres secciones estructurales para la región
(perfiles A, By C en la figuras 4.5 y 4.13) donde el fallamiento tipo inverso pareciera ser el dominante.
El poder determinar si dicho sistema idealizado de fallas se encuetran activos en la actualidad es
casi responsabilidad única de las soluciones de mecanismos de foco. Lamentablemente las soluciones
obtenidas en este estudio no son lo suficientemente grandes como para poder hacer una población
de mecanismos y así poder determinar con mayor precisión los ejes P , T. Sin embargo, haciendo
correlaciones entre los modelos de Diraison et al. [2000] y las soluciones hipocentrales, es posible ver
ciertas analogías (Figura 4.14). A partir de la premisa anterior se podría concluir que los sistemas
de fallas del perfil B-B ' podrían ser los causantes de la actividad sísmica originado en las cercanías
de dicho perfil.
No es posible decir lo mismo de los otros dos perfiles (A-A ' y C-C') ya que no cruzan zonas con
densidades hlpocentrales lo suficientemente grandes como para hacer el mismo ejercicio del perfil
B-B '.
'ITabajos hechos por Lodolo et al. [2002) dan cuenta de una campaña sísmica en el extremo
Sur-Oriente de la isla Tierra del Fuego, específicamente sobre la falla Magallanes-Fagnano. Para tal
efecto registraron sísmica de reflexión utilizando receptores de 96 canales y dos cañones aire.
La presencia de la falla Magallanes-Fagnano queda de manifiesto en los perfiles obtenidos, que
logran cortar dicha falla, es así como la figura 4.15 nos muestra una cuenca asimétrica asociada a una
falla principal. Dichos perfiles nos muestran la falla, pero son incapaces de mostrar en qué grado
de actividad se encuentran y si es que realmente han sido los causantes de los grandes eventos
registrados en la región.
De igual forma Lodolo et al. [2003] se atreve a postular la Falla Magallanes-Fagnano como un
sistema de fallas de rumbo que se entrecorta en el sector de Tierra del Fuego para luego seguir como

66
4. R.esultados:Análisis y Discusión 4.5 Perfiles slsmicos regionales

una falla única hasta el triple punto de las placas Antártica-Scotia-Sudamérica 4.15. Dicho sistema
de fallas también fue postuldao por Diraison et al. (1997] y Klepeis y Austin (1996).
El incremento de perfiles sísmicos sobre la falla de Magallanes ayudará. en futuro a una mayor
comprensión de la forma que sigue esta falla.

Figura 4.12: Ubicación de cuatro secciones regionales obtenidas a partir de sfsmica e


intrepretación geológica superficial.

67
4. Resu.ltados:Análisis y Discusión 4.5 Perfiles slsmicOB regionales

...

--
1\
- Tcrcomo
- e,..;.,.,~
CJ Cn.úc:ica l.úcrior
SOicm 0 -..-o

PtY f»}t([i
8 8"

~ 1;
t
SC.I km ---- .............................. ·-------------

Figura 4.13: Interpretación estructural para tres secciones de la Región de Magallanes


(Diraison et al., 2000).

B'

·.

·.
·.
·.

10

Figura 4.14: Detalle de la sismicidad superficial en los alrededores de la sección B-B'

68
4. Resultados:Análisis y Discusión 4.5 Perfiles sfsmicos regionales

~. _

Figura 4.15: Detalle de una sección sísmica (96 canales) sobre el borde SE de Tierra
del Fuego (Sección D-D'). La interpretación muestra una cuenca asimétrica
cortada por una falla principal (Magallanes-Fagnano) [Lodolo et al. , 2002].
La figura inferior muestra el modelo propuesto por Lodolo et al. [2003) hecha
en base a 3 perfiles sísmicos que cortan la falla de Magallanes.

69
Capítulo 5

Conclusiones

La tectónica propuesta para la región austral de Sudamérica ha sido ampliamente fortalecida


durante los últimos veinte años, definiendo así zonas de contacto entre placas (Scotia-Antártica-
Sudamérica) tales como zonas de subducción y jo movimientos transformantes que son avalados
principalmente por estudios geológicos y geofísicos. Hasta la fecha del presente estudio, sólo se
contaba con datos de telesismicidad, representando de esta forma una zona de baja actividad sísmica,
a diferencia de lo que ocurre en el centro y norte de Chile, donde la interacción de las placas de Nazca
y Sudamérica son las causantes de un gran número de eventos. Los principales rasgos geomorfológicos
para la zona de estudio muestran una subducción entre la placa Antártica y la placa Sudamericana
y una subducción oblicua entre las placas Antártica y Scotia. De esta misma forma, una de las
principales características que puede señalarse para la región es la presencia de un sistema de fallas
de rumbo entre las placas Sudamericana y de Scotia. A dicho sistema de fallas se le reconoce
como falla de Magallanes o Magallanes-Fagnano, la que se orienta siguiendo la boca occidental del
Estrecho de Magallanes, continuando a través del Seno Almirantazgo y Lago Fagnano en el sector de
Tierra del Fuego. Lo anterior estaría actuando como responsable de los principales eventos sísmicos
(terremotos) de la región.
Resulta fundamental contar con datos de sismología a nivel local que permitan estudiar en
forma profunda y detallada los rasgos dinámicos de las tres placas involucradas. Es por esto que
la red de estaciones sismológicas del proyecto SEPA registró, durante un periodo de dos años, la
microsismicidad local dentro de la región de Magallanes, obteniendo por primera vez un mapa con la
distribución epicentral que no era posible lograr desde la Red Mundial de Estaciones Sismológicas.
La secuencia de trabajo incluyó desde la adquisición cruda hasta la creación de una base de datos
con registros digitales, involucrando en el camino una serie de otros procesos, tales como:
- la selección visual de los eventos, por tratarse de registros continuos

70
5. Conclusiones

- uso de un filtro adecuado


- elección de modelos de velocidades, entre otros.
La selección visual se realizó debido a la presencia de distintos ruidos ambientales, tales como el
fuerte viento y la cercanía de vías urbanas, siendo ineficiente el uso de un sistema de Trigger.
Por tratarse de estaciones de banda ancha, la información recolectada osciló en un amplio rango
del espectro, registrando de esta forma tanto eventos locales como telesísmicos, situación que sólo era
descubierta mediante el uso de filtros de alta frecuencia para los eventos locales y de baja frecuencia
para los lejanos.
Por su parte, el modelo de velocidad deducido en esta tesis se obtuvo utilizando un programa
iterativo, donde la función objetivo se concentra en minimizar los errores residuales, tomando como
función de entrada un modelo de velocidad inicial y los tiempos de llegada de las fases P y S para
las distintas estaciones.
El modelo final da como resultado un modelo de tierra plana de tres capas, encontrándose la
discontinuidad manto-corteza a los 32 km. de profundidad. La distribución hipocentral obtenida
a través de este modelo final de velocidades muestra una región en que la sismicidad se agrupa
principalmente en la faja plegada y corrida de Magallanes, siendo más notorio en las zonas aledañas
al área de Puerto Na tales y Seno Skiring y una sismicidad menor asociado a los fiordos y canales
cercanos a la boca occidental del Estrecho de Magallanes. De igual forma, se observa una sismicidad
difusa en las zonas llanas de la Provincia de Santa Cruz (Argentina), lo que resulta un hecho curioso
debido a que esta zona se encuentra alejada de los estructuras tectónicas propuestas para la región.
A lo largo de la Falla de Magallanes es posible observar lineamientos que siguen la orientación
de la misma, junto con otro grupo periférico que la interceptan en forma oblicua. Sin embargo, no es
posible concluir que dicha falla sea el agente causante de los eventos sísmicos de la región, tal como
lo sugiere la literatura. Tal vez esta falla sólo se manifiesta en aquellos eventos sísmicos de mayor
magnitud.
Uno de los rasgos a destacar es la profundidad hipocentral en que se encontró un grupo de even-
tos, variando entre los O y los 5km. La aparición de estos eventos someros presentan características
análogas, tanto por la similaridad en sus formas de ondas, por el periodo de tiempo en que fueron
registrados (Abril-Julio, 1998) y por la cercanía de su solución epicentral. Estas soluciones epicen-
trales coinciden geográficamente con la presencia del volcán Reclus, ubicado en la zona terminal
de Campo de hielo Sur que, por lo remoto de su ubicación, se desconoce su actividad. Por todo lo
anterior, se postula que dicha actividad se origina por cierta actividad generada en los alrededores
del volcán.
En cuanto a la magnitud de los eventos totales registrados en la zona de estudio, se puede ver que
la mayoría de los eventos fueron imperceptibles para la población siendo sus magnitudes menores a

71
5. Conclusiones

4. Sólo dos eventos registrados durante el periodo del experimento alcanzaron magnitudes mayores
a 4, localizándose dichos epicentros en la zona de Puerto Natales y al Norte de la ciudad de Punta
Arenas. Es notable observar que sólo tres de los trescientos eventos producidos en la región fueron
registrados por la Red Mundial de Estaciones Sismológicas.
Debido a los puntos antes mencionado se puede postular a la región austral de Sudamérica
como una zona sísmicamente activa, donde los eventos sísmicos son producidos principalmete por
la subducción de la placa Antártica bajo la placa Sudamericana. Esto es ratificado por los dos
mecanismos de foco obtenidos durante el experimento, los que poseen componentes compresionales
del tipo falla inversa. Asimismo, los tres mecanismos focales obtenidos a través de la red mundial
muestran el carácter de falla de rumbo producido en el sector de la falla de Magallanes.
En términos de tectónica se postula a la falla de Magallanes como eje principal para una secuencia
de fallas laterales (Riedel shear) que serían las encargadas de disipar la energía acumulada en la falla
principal, siendo necesario ratificar lo último mediante una solución más acabada de los mecanismos
de foco involucrados en tales lineamientos.

72
Apéndice A

Función de transferencia del sensor


STS-2

A continuación, se deriva la función de transferencia del sismógrafo Streckeisen STS-2 y las


respuestas en fase y amplitud del instrumento. Dicho sensor contiene 3 sensores ortogonales inde-
pendientes orientado 45°de la vertical. La respuesta vertical y horizontal son sumas vectoriales de
estas tres señales. El STS-2 posee una amplia respuesta en frecuencia alcanzando hasta 100 segundos
de período. Esto es posible gracias a la suspensión mecánica usada para la masa que es cercana a
cero Hertz.
La función de transferencia entre el movimiento del suelo y la salida de voltage es escrita en
términos de la variable de Laplace s = a + iw. El problema es dividido en dos partes: primero el
movimiento de la masa M con respecto al instrumento y la señal resultante nos da la función de trans-
ferencia V(s)/Xh(s); luego el movimiento del sistema inercial con respecto al suelo Xh(s)fX 9 (s).
La función total de transferencia es el producto de estas dos funciones:

V(s) = Xh(s) V(s)


(A.1)
JC9(s) JC9(s) )(h(s)

Ecuación del Movimiento de Masa:


.
'(t)
dv
= C- v
+-
dt Rt
1
+R2r
- ¡ vdt (A.2)

(A.3)

k¡i = -Mxm (A.4)

73
A. Función de transferencia del sensor STS-2

v (s) -M s3 ¡ k1 e
(A.5)
X()
h S
M
k 1 k 2 Cs
3+ S 2+ RlCS
1 + R2Cr
1

Ecuación del sistema inercial:

(A.6)

(A.7)

Xh(s) k+ bs wg + 2wo(s
(A.8)
X 9 (s) = Mhs 2 + bs +k = s 2 + swo(s + wg
Aquí, wfi =k/Mh y 2w0 ( =b/Mh . Quedando la función de transferencia con respecto a la
velocidad del terreno X9 (s) = sX9 (s) como:

V(s) ( wg+2wo(s ) ( -Ms 2 /k 1C ) (A.9)


X9 (s) = 2
s + 2wo(s + wfi . k1t;cs 3 + s2 + R~cs + rdc:r
El manual Streckkeisen especifica las siguientes constantes:

M= 0.3 Kg
e= 7.8 micrones
k1 = 50 Newtons/amp
R 1 = 1.7 meg ohm
R2 = 0.6 ohm

Resultando un factor de escala para el instrumento de:

k~c=769.2 volt secfm

74
Apéndice B

Determinación de los parámetros


hipocentrales y modelos de velocidad

Para determinar los parámetros hipocentrales se hace necesario explicar el modelo matemático
que sostiene la teoría de rayos y su posterior aplicación al modelo acoplado: hipocentro-modelo de
velocidad.

B .l. Tiempo de viaje de un rayo simple


La teoría de rayos es uno de los métodos más eficientes para calcular los tiempos de viaje entre
la fuente sísmica y un receptor para un campo de velocidades. En tal teoría, los tiempos de arribo de
las ondas sísmicas generadas por un terremoto y que luego son observados por una estación sísmica
son una función de los cuatro parámetros hipocentrales (tiempo origen y coordenadas espaciales del
hipocentro) y el campo de velocidades usado para cada forma de rayo:

tobs = F(to,xo, Yo ,zo, V( z)) (B.l)

Donde
tobs tiempo de arribo observado;
to tiempo origen;
xo, Yo, zo coordenadas del hipocentro;
V(z) campo de velocidades.

En general, no se conocen ni los parámetros hipocentrales ni los campos de velocidades, siendo


sólo medible los tiempo de arribo en B.l. Usando un simple modelo de velocidades promedio, es

75
B. Determinación de los parámetros hipocentrales y modelos de velocidad B.l Tiempo de viaje de un rayo simple

posible trazar rayos desde la fuente hacia los receptores y de esta forma calcular los tiempos teóricos
de llegadas:

tcalc = F(t~,x 0 ,y0 ,z0 , V*(z)) = F(hj,mjJ (B.2)


Donde tcatc es el tiempo teórico de viaje, y los parámetros estimados hipocentrales son:
hj = t0, x 0,y 0, z0 y mk = V*(z) como el parámetro de velocidad estimado. La diferencia entre los
tiempos de arribo observado y calculado se denomina tiempo de viaje residual. Obviamente este
tiempo de viaje residual es función de los parámetros estimados y verdaderos de la velocidad y del
hipocentro.

tres= toba- tcalc = F(dto , dxo , dzo , dV(z)) (B.3)


Donde tres es el tiempo de viaje residual y d es la diferencia entre los parámetros desconocidos
verdaderos y estimados.
Para calcular los ajustes a los parámetros hipocentrales o de velocidades {d), necesitamos conocer
la dependencia de todos estos parámetros sobre los tiempos de viaje observados. Para todos estos
parámetros, excepto el tiempo origen, se tiene una dependencia no lineal . Usando series de expansión
de Taylor , se puede escribir la ecuación B.l como una aproximación lineal:

tobs = tcalc + F(hj, mk) =


4
8F n 8F
=L 8hh. +L 8mk · ó.mk (B.4)
j=l k=l

Gf¡ (para hj = hj) es la derivada parcial del tiempo de viaje con respecto al ith parámetro
hipocentral y J!" (para mk = mk) es la derivada parcial del tiempo de viaje con respecto al kth
parámetro del modelo de velocidad. Usando la ecuación B.3 podemos reescribir el tiempo de viaje
residual (tres) de una observación con ajustes (d) a los parámetros hipocentrales y de velocidad:

4
8F
tres = tobs - tcalc = .L: 8 · ó.mk {B.5)
j=l k
El cálculo de tcalc y las derivadas parciales es conocido como el problema directo, mientras que
el cálculo de D.hj y ó.mk se le denomina el problema inverso.

B .l.l. El problema acoplado: Hipocentro y modelo de velocidad

El problema acoplado de los parámetros hlpocentrales y de velocidad se refiere a la determinación


de los parámetros hlpocentrales desconocidos y los parámetros de velocidad desde un conjunto de

76
B. Determinación de los parámetros hipocentrales y modelos de velocidad B.l Tiempo de viaje de un rayo simple

tiempos de arribos de varios eventos sísmicos. El acoplamiento de dos problemas inversos es visible
en cada camino de rayo unitario. Dado un conjunto de tiempos de arribos desde un evento y usando
una estimación para la ubicación hipocentral y tomando un modelo inicial de velocidades, podríamos
obtener un sistema de ecuaciones de la forma de la ecuación B.5:

8F 8F 8F 8F 8F
· ~zo + L -- ·~mk
n
ti = - · ~to +- · ~xo + - · ~Yo+ - (B.6)
8t ax 8y 8z k=l 8mk

donde ti es el tiempo residual; i= 1, . . . ,nobs, con npbs es el número total de observaciones para
este evento. Este sistema de (nobs) ecuaciones por un evento puede ser escrito en notación matricial
como:
t= Ad (B.7)

donde t es el vector de tiempos de viaje residuales, A la matriz con derivadas parciales , y d el


vector de ajuste para los parámetros hipocentrales y modelos de velocidad.
Si asumimos que el modelo de velocidad inicial es correcto y lo utilizamos en la la ecuación B.6,
consideraremos el problema que es ampliamente resuelto por procedimientos de localización de
eventos sísmicos tal como los programas HYP071 [Lee y Lahr, 1985], HYPOINVERSE [Klein, 1978],
o HYPOCENTER [Lienert et al., 1988]. Con cuatro incógnitas y un gran número de observaciones
(NOBS > 4), este sistema lineal de ecuaciones es transformado usualmente en un problema de
mínimos cuadrados que puede ser resuelto por varios métodos numéricos. El programa VELEST
usa el método de descomposición de Cholesky. Con la introducción de los parámetros de velocidad
desconocidos en la ecuación B.7, se necesita resolver el sistema de ecuaciones en forma simultánea
para un número de eventos. Esta formulación conduce a un número total de observaciones (nobsto)
más grande que el número total de incógnitas, las que han aumentado a 4 veces el número de eventos
más el número de parámetros del modelo de velocidad. Agregando un vector de errores que contenga
la parte restante de los residuales de los tiempos de viaje que no pudieran haber sido desarrollados
por ajustes de los parámetros desconocidos, obtenemos:

t = Ad +e (B.8)

d= (ahl,j! a h2,j, ... ahnev,j, am¡, a m2, ... amn) j = l. .. 4


nev número total de eventos sísmicos
n número total de los parámetros del modelo de velocidad desconocido
índice para las observaciones (residuales del tiempo de viaje)
e vector que contiene los errores de linealización y de tiempo.

77
B. Determinación de los parámetros hipocentrales y modelos de velocidad B.l Tiempo de viaje de un rayo simple

_!lli_ ...!lli_
8hl,j
o o _!!k_
8m1 8mn

o
A=
o _!!!;__ _!!k_ ...!lli_
8hnev,j 8m1 8mn

B.1.2. Separación de los Parámetros Hipocentrales y del modelo de velocidad

La matriz de derivadas parciales definido anteriormente se puede separar en dos matrices:

A=H+M (B.9)

con d como el vector de ajustes a todas las variables desconocidas, h como el el vector de ajuste
a los parámetros hipocentrales, y finalmente m como el vector de ajuste a los parámetros del modelo
de velocidades.

B.1.3. Solución Numérica al Problema Inverso Separado

Por la separación de los parámetros hipocentrales del modelo de velocidad, obtenemos del pro-
blema inverso un sistema condicional de ecuadones, las que contienen los parámetros del modelo de
velocidad como incógnitas:

(B.lO)

donde ea es el vector de error que podría ser incluido en el vector transformado del tiempo
de viaje residual t 11 • M 8 es la matriz transformada por la separación del modelo de velocidad. La
separación se lleva a cabo por la multiplicaión de la ecuación B.8 por una matriz ortogonal a H (ver
ecuación B.9).
Se resuelve la aproximación linealizada de este problema no lineal por un procedimiento interactivo,
donde en cada iteración la t area es resolver un sistema de ecuaciones lineales (ecuación B.lO) y
aplicar los ajustes de velocidad resultante al campo de velocidad. mk es la solución de la kth
iteración, obteniendo finalmente el modelo de velocidades:

m f inal = meatimada + m¡ + m2 + ... + m k + ... (B.ll)

78
B. Determinación de los parámetros hipocentrales y modelos de velocidad B.l Tiempo de viaje de un rayo simple

La matriz de derivadas parciales (Mk) y el vector de residuales de tiempo de viaje (tk) tiene
que ser calculada para cada paso de la iteración usando los parámetros actualizados del modelo de
velocidad (mk-1)· Para un sistema de ecuaciones lineales (B.lO) la solución de mínimos cuadrados
es dado por la ecuación normal

(B.l2)

y usando la notación de inversa generalizada [Aki et al., 1976], se obtiene la solución por la
inversión de la matriz simétrica (MTM).

(B.l3)

79
Bibliografía

Barker, P. F. Plate tectonics in the scotia sea region. Nature, 228, págs. 489-503, 1970.

Bruce, R. M., E. P. Nelson, S. G. Weaver y D. R. Lux. Temporal and spatial variations in the
southern patagonian batholith: Contrains on magmatic are development. Geol. Soc. Am. Spec.
Paper, 265, págs. 1-12, 1991.

Cande, S. C. y R. B. Leslie. Late cenozoic tectonics of southern Chile trench. Journal oj Geophysical
Research, 91(B1), págs. 471- 496, 1986.

Cunningham, W . D. Orogenesis at the southern tip of the americas: the structural evolution of the
cordilleradarwin methamorphic complex, southernmost chile. Tectonophysics, 244, págs. 197- 229,
1995.

Cunningham, W. D. , K. A. Klepeis, W. A. Gose y l. W . D. Dalziel. The patagonian orocline: New


paleomagnetic data from the andean magmatic are in tierra del fuego, chile. J. Geophys. R es.,
196, págs. 16 ,061- 16,069, 1991.

Dalziel, l. W. D. Back-arc extension in the southern andes: A review and critica} reappraisal. Phil.
TI-ans. R . S oc. Lond. A 300, págs. 319-335, 1981.

Dalziel, l. W. D. , M. F. D. Wit y K. F. Palmer. Fossil marginal basin in the southern andes. Nature,
250, págs. 291- 294, 1974.

DeMets, C., R. G. Gordon, D. F. Argus y S. Stein. Current plate motions. Geophys. J. Int. , 101,
págs. 425- 478, 1990.

DeMets, C., R. G. Gordon, D. F. Argus y S. Stein. Effect of recent revisions to the geomagnetic
reversa} time scale on estimates of current plate motions. Geophys. R es. Let., 21(20), págs. 2191-
2194, 1994.

Diraison, M. , P. R. Cobbold, D. Gapais y E. A. Rossello. Magellan Strait: Part of a Neogene rift


system . Geology , 25(8), págs. 703- 708, 1997.

80
Diraison, M., P. R. Cobbold, D. Gapais, E. A. Rossello y C. L. Corre. Cenozoic crustal thickening,
wrenching and rifting in the foothills of the southernmost Andes. Tectonophysics, 316, págs. 91-
119, 2000.

Forsyth, D. W. Fault plane solutions and tectonics of the south atlantic and scotia sea. J. Geophys.
Res., 80, págs. 1429- 1443, 1975.

Goldstein, P. Sac: Seismic analysis code. Inf. téc., Lawrence Livermore National Laboratory, 1996.

Gonzalez-Ferran, O. Volcanes de Chile. Instituto Geográfico Militar, Santiago, Chile, 1995.

Gripp, A. E. y R. G. Gordon. Current plate velocities relative to the hotspots incorporating the
nuvel-1 global plate motion model. Geophys. Res. Lett, 117, págs. 1109- 1112, 1990.

Harvard-CMT. Seismology, CMT Catalog Search. http:/ /www.seismology.harvard.edu/CMTsearch.html,


2002.

Hervé, F., E.Ñelson, K. Kawashita y M. Suárez. New isotopic ages and the timing of orogenic events
in the Cordillera Darwin, southernmost Chilean Andes. Earth Planet. Sci. Lett., 55, págs. 257- 265,
1981.

Hervé, F., M. Suárez y A. Puig. The patagonian batholith s of tierra del fuego, chile: Timing and
tectonic implication. J. Geol. Soc. London, 141, págs. 909- 917, 1984.

Kennett, B. L.Ñ. y E. R. Engdahl. Traveltimes for global earthquake location and phase identifica-
tion. J. Geophys. Res. , 122, págs. 429- 465, 1991.

Kissling, E. Geotomography with local earthquake data. Reviews of Geophysics , 26, págs. 659-698,
1988.

Klein, F. W. Hypocenter location program hypoinverse. Inf. téc., Open File Report , U. S. Geological
Survey, 1978.

Klepeis, K. A. y J. A. Austin. Contrasting styles of superposed deformation in the southernmost


andes. Tectonics , 101, págs. 20,211- 20,231, 1996.

Lee, W. H. K. y C. M. Lahr. Hyp071 : A personal computer version of the hypo71 earthquake


location program. Inf. téc., Open File Report, U. S. Geological Survey, 1985.

Lienert, B. R. E. , E. Berg y L.Ñ. Frazer. Hypocenter: An earthquake location method using centered,
scaled, and adaptively least squares. Bull. Seism. Soc. Am., 76, págs. 771- 783, 1988.

81
Lodolo, E., M. Menichetti, A. Tassone, R. Geletti y H. Lippai. Morphostructure of the south
america-scotia transform margin in the tierra del fuego region, 2003.

Lodolo, E., M.Menichetti, A. Tassone, R. Geletti, P. Sterzai, H. Lippai y J. Hormachea. Researchers


target a continental transform fault in tierra del fuego. EOS, 83(1) , 2002.

Lomnitz, C. Geologische rundschau. En S. Springer Verlag, ed., Major Earthquakes and Tsunamis
in Chile, págs. 951- 955. 1970.

Ludwing, W . J. , J. l. Ewing y M. Ewing. Seismic-refraction measurments in the magellan straits.


Joumal of Geophysical R esearch, 70(8), págs. 1855- 1975, 1965.

Mpodozis, C. y V. Ramos. The andes of chile and argentina. En Geology of the Andes and
its relation to Hydrocarbon and Mineral resources, ed., Circum-Pacific Council for Energy and
Mineral resources Earth Science Series , tomo 11, págs. 59- 90. G. E. Ericksen, M. T. Cañas and
J. A. Reinemund, 1990.

NEIC. National Earthquake Information Center. http:/ jwwwneic.cr.usgs.govjneis/epic/epic.html,


2002.

Pelayo, A. M. y D. A. Wiens. Seismotectonics and relative plate motion in the scotia sea region. J.
Geophys. Res. , 94, págs. 7293- 7320, 1989.

Ramos, V. A. The Birth of Southern America. Amer. Scien, 77, págs. 444-450, 1989.

Smalley, R. , E. Kendrick, M. Bevis, l. Dalziel, F. Taylor, E. Lauría, R. Barriga, G. Casassa, E. Olivero


y E. Piana. Geodetic determination of relative plate motion crustal deformation across the scotia-
southh america plate boundary in eastern tierra del fuego, 2003.

Vera, E. E. , G. A. Yáñez, A. Díaz y M. Torne. Tectonic framework of the southernmost chilean


margin: Gravimetric and seismic constraints. En Simposio Geología Marina de Chile José Corvalán
Díaz. VIII congreso Geológico Chileno, tomo 111, págs. 1849-1853. Universidad Católica del Norte.
Departamento de Ciencias Geológicas, 1997.

Wilson, T. J. Transition from back-arc to foreland basin developement in the southernmost Andes:
Stratigraphic record from Ultima Esperanza District, Chile. Geol. Soc. Am. Bull, 103, págs. 98-
111, 1991.

Winslow, M. A. The structural evolution of the magallanes basin and neotectonics in the southern-
most andes. En Antartic Geoscience. 1982.

82

También podría gustarte