AREQUIPA
AREQUIPA
AREQUIPA
AREQUIPA – PERÚ
2018
UNIVERSIDAD NACIONAL DE SAN AGUSTÍN DE AREQUIPA
AREQUIPA – PERÚ
2018
UNIVERSIDAD NACIONAL DE SAN AGUSTÍN DE AREQUIPA
JURADOS
PRESIDENTE: Msc. Edgar Gonzales Zenteno
VOCAL: Dr. Orlando Macedo Sánchez
SECRETARIO: Ing. Javier Ticona Paucara
AREQUIPA – PERÚ
2018
Dedicado a María y Luis
Mis padres
AGRADECIMIENTOS
En primer lugar le doy gracias a Dios por la vida. Durante el tiempo de la realización de
la presente investigación recibí el apoyo de muchas personas, cada una me dejó una
enseñanza y un recuerdo.
Un agradecimiento especial al Presidente Ejecutivo del Instituto Geofísico del Perú Dr.
Hernando Tavera Huarache quién depositó su confianza en mí, me orientó en el
desarrollo de la investigación y siempre estuvo presto a responder mis dudas. Quiero
agradecer a las personas que laboran en la Sub- Dirección de Ciencias de la Tierra
Sólida y en la Sub-Dirección de Redes Geofísicas del Instituto Geofísico del Perú.
Asimismo agradecer al Dr. Armando Minaya Lizárraga asesor interno del trabajo de
investigación.
En los años que viví en la ciudad de Lima me brindaron su mano amiga las señoras
Alicia Flores y María Elena Martínez, siempre me trataron como una hija, gracias por su
amistad y cariño a mis compañeras de piso Isela, Carmen, Janeth y Diana siempre las
recordaré.
ÍNDICE
CAPÍTULO I ................................................................................................................... 1
INTRODUCCIÓN........................................................................................................... 1
CAPÍTULO II ...................................................................................................................... 9
GEODINÁMICA............................................................................................................. 9
2.1. CONTEXTO GEODINÁMICO ................................................................................ 9
2.1.1. La Dorsal de Nazca ......................................................................................... 10
2.1.2. La Fractura de Mendaña ................................................................................. 10
2.1.3. La Fosa Perú-Chile ......................................................................................... 10
2.1.4. La Cordillera Andina ......................................................................................... 11
2.1.5. La Cadena Volcánica ...................................................................................... 13
2.1.6. Sistemas de Fallas ........................................................................................... 13
CAPÍTULO IV .............................................................................................................. 26
TERREMOTO DEL 23 DE JUNIO DEL 2001 .......................................................... 26
4.1. PARÁMETROS DEL HIPOCENTRO Y RÉPLICAS ............................................ 26
4.2. DISTRIBUCIÓN ESPACIAL DE RÉPLICAS ....................................................... 27
4.3. ORIENTACIÓN DE LA RUPTURA ...................................................................... 30
4.3.1. Sismo Principal ............................................................................................... 31
4.3.2. Réplica del 26 de Junio ...................................................................................... 32
4.3.3. Réplica del 5 de Julio ......................................................................................... 33
4.3.4. Réplica del 7 de Julio ......................................................................................... 33
4.4. TSUNAMI ............................................................................................................... 35
CAPÍTULO V ................................................................................................................ 37
DINÁMICA Y PARÁMETROS DE LA FUENTE SÍSMICA.................................. 37
5.1. MODELOS DE FUENTE SÍSMICA ...................................................................... 37
5.1.1. Modelo de Barreras ........................................................................................ 37
5.1.2. Modelo de Asperezas ..................................................................................... 38
5.2. MODELOS DE RUPTURA CIRCULAR ............................................................... 39
5.2.1. Modelo de Brune ............................................................................................ 39
5.2.2. Modelo de Madariaga ..................................................................................... 42
5.3. PARÁMETROS DE LA FUENTE SÍSMICA ........................................................ 43
5.3.1. Momento sísmico escalar ............................................................................... 43
5.3.2. Radio de fractura ............................................................................................ 44
5.3.3. Caída de esfuerzos .......................................................................................... 45
5.3.4. Energía sísmica ............................................................................................... 47
CAPÍTULO VI .............................................................................................................. 49
ANÁLISIS DE SEÑALES SÍSMICAS ........................................................................ 49
6.1. REGISTROS SÍSMICOS ........................................................................................ 49
6.2. CORRECCIONES DE LOS REGISTROS SÍSMICOS .......................................... 50
6.2.1. Corrección de línea base ................................................................................. 50
6.2.2. Corrección instrumental ................................................................................. 50
6.3. EL RUIDO SÍSMICO .............................................................................................. 53
6.3.1 Ruido cultural ................................................................................................. 53
6.3.2 Ruidos de Periodo Corto................................................................................. 54
6.3.3 Ruidos de Periodo Largo ................................................................................ 55
6.3.4 Filtros .............................................................................................................. 55
6.3.4.1 Tipos de filtros .................................................................................... 57
6.4. TRANSFORMADA RÁPIDA DE FOURIER ........................................................ 59
6.5. ANÁLISIS ESPECTRAL ........................................................................................ 61
BIBLIOGRAFÍA
ANEXOS
ANEXO 1
ANEXO 2
LISTA DE FIGURAS
Figura 1.1. Ubicación del área de estudio, epicentro del terremoto del 2001 y su serie
de réplicas.......................................................................................................................... 5
Figura 2.2. Mapa de ubicación de los principales sistemas de fallas del Perú (Macharé
et al., 2003)...................................................................................................................... 15
Figura 3.2. Mapa sísmico de Perú periodo 1960-2017. Fuente Instituto Geofísico del
Perú ................................................................................................................................. 23
Figura 4.1. Mapa de isosistas del terremoto de Arequipa (Tavera et al., 2002) ............. 27
Figura 4.2. Distribución espacial del epicentro del evento principal y réplicas
relocalizadas del terremoto del 23 de Junio del 2001. Los mecanismos focales fueron
obtenidos usando registros telesísmicos. El área rayada corresponde a las réplicas del
terremoto de Noviembre de 1996 (Tavera et al., 1998) y la sombreada, al área que no ha
experimentado ruptura durante el terremoto del 23 de Junio. Los grupos de réplicas
están delimitados por líneas punteadas. La fecha y magnitud están ubicadas encima de
cada mecanismo focal (Tavera et al., 2006) .................................................................... 29
Figura 4.5. Resultados del modelamiento de ondas P para las réplicas del terremoto de
Arequipa a) 26 de Junio, b) 5 de Julio y c) 7 de Julio (Tavera et al., 2002) .................. 34
Figura 4.6. Principales parámetros y características del tsunami que afectó la localidad
de Camaná en Arequipa. Las líneas perpendiculares a la línea de Costa indican la
máxima altura de las olas. El área sombreada corresponde al área inundada por el
tsunami (Tavera et al., 2006) ........................................................................................... 36
Figura 5.3. Sismogramas y espectros para cinco réplicas del terremoto de San Fernando
(1971). Las flechas debajo de los sismogramas indican la parte de la onda S utilizada
para el espectro de Fourier (Brune, 1973) ....................................................................... 41
Figura 6.2.a) Registro de velocidad sin corregir para una réplica del terremoto de
Arequipa del 2001, obsérvese que el registro se encuentra por debajo de la línea base; b)
Registro de velocidad, al cual se aplicó la corrección de línea base ............................... 50
Figura 6.3. Representación del cálculo de un sismograma teórico para una fuente
puntual en el campo lejano. La función temporal de la fuente sísmica (F(t)) se combina
con la atenuación del medio (Q(t)) y el instrumento (I(t)) mediante una operación de
convolución (Tavera, 1998) ............................................................................................ 51
Figura 6.5.a) Registro de la réplica del sismo de Arequipa del 2001 sin corrección
instrumental, la amplitud está expresada en cuentas/s. b) El mismo registro, al cual se
aplicó la corrección instrumental, la amplitud expresada en cm/s .................................. 53
Figura 6.6. Registros de réplicas del sismo de Arequipa con diferentes niveles de ruido
sísmico (rectángulo sombreado) ..................................................................................... 54
Figura 6.8. Respuesta frecuencial del filtro Paso-Alto, (fc) es la frecuencia de corte 57
Figura 6.13. Espectros para dos réplicas del terremoto de Arequipa, la primera
registrada el 29 de Junio en la estación Tocota y la segunda registrada el 01 de Julio en
la estación Caravelí. En el eje de las abscisas, la frecuencia está en Hz y en el eje de
ordenadas, las amplitudes en cm.se, Ω es la parte plana del espectro y fc la frecuencia
esquina............................................................................................................................. 61
Figura 7.1. Epicentros de la serie de réplicas del terremoto de Arequipa del 2001 64
Figura 7.2. Distribución de estaciones sísmicas que integran la Red Sísmica Local. La
estrella indica el epicentro del sismo principal................................................................ 65
Figura 7.3. Ejemplos de espectros de ondas S para el desplazamiento del suelo
correspondiente a réplicas registradas en la estación Tocota, N-S (Componente Norte-
Sur) y E-O (Componente Este-Oeste) ............................................................................. 69
Figura 7.8. Relación entre el logaritmo de Mo y r para las réplicas del terremoto de
Arequipa, registradas en la estación Tocota .................................................................... 80
Figura 7.9. Relación entre el logaritmo de Mo y r para las réplicas del terremoto de
Arequipa, registradas en la estación Ocoña .................................................................... 81
Figura 7.10. Relación entre el logaritmo de Mo y r para las réplicas del terremoto de
Arequipa, registradas en la estación Caravelí ................................................................. 82
Figura 7.11. Relación entre el logaritmo de Mo y r para las réplicas del terremoto de
Arequipa, registradas en la estación Bella Unión ........................................................... 83
Figura 7.12. Relación entre el logaritmo de Mo y r para las réplicas del terremoto de
Arequipa, registradas en la estación Quilca .................................................................... 84
Figura 8.1. Relación entre área de ruptura y momento sísmico. Las líneas 0.1, 1, 10 y
100 MPa muestran una caída de esfuerzos (Δσ) constante, (Kanamori y Anderson,
1975) ............................................................................................................................... 87
Figura 8.2. Relación entre el logaritmo del momento sísmico escalar y el área de la
fractura para terremotos con foco superficial, (Tavera, 1998) ........................................ 88
Figura 8.3. Relación entre el logaritmo del momento sísmico escalar (Mo) y el radio de
la fractura (r), (Abercrombie y Leary, 1993) ................................................................. 89
Figura 8.4.a) Relación entre el logaritmo del momento sísmico escalar (Mo) y el radio
de la fractura (r), Estación Tocota. Modificado Abercrombie y Leary, (1993) ............. 90
Figura 8.4.b) Relación entre el logaritmo del momento sísmico escalar (Mo) y el radio
de la fractura (r), Estación Ocoña. Modificado Abercrombie y Leary, (1993) .............. 91
Figura 8.4.c) Relación entre el logaritmo del momento sísmico escalar (Mo) y el radio
de la fractura (r), Estación Caravelí. Modificado Abercrombie y Leary, (1993) ........... 92
Figura 8.4.d) Relación entre el logaritmo del momento sísmico escalar (Mo) y el radio
de la fractura (r), Estación Bella Unión. Modificado Abercrombie y Leary, (1993) ..... 92
Figura 8.4.e) Relación entre el logaritmo del momento sísmico escalar (Mo) y el radio
de la fractura (r), Estación Quilca. Modificado Abercrombie y Leary, (1993) .............. 93
Tabla 7.1. Coordenadas de estaciones sísmicas de banda ancha que integran la red
sísmica local .................................................................................................................... 66
Los parámetros de la fuente sísmica de las réplicas del Terremoto de Arequipa del
23 de Junio del 2001 (Mw=8.2), son determinados a partir de la aplicación de la
técnica de análisis espectral. Los datos utilizados corresponden a 292 réplicas
registradas por la Red Sismológica Nacional y una Red Sísmica Local integrada por
9 estaciones a cargo del Instituto Geofísico del Perú. Todas las réplicas presentaron
magnitudes ML menores a 5.5 y focos ubicados a profundidades menores a 60Km.
Después de corregir por línea base y efecto instrumental, los espectros de Fourier
fueron obtenidos para el tren de ondas S de acuerdo a Keilis-Borok (1959) y Brune
(1970).
Los resultados obtenidos muestran que las réplicas presentan valores para el
momento sísmico escalar en el rango de 3.45x1012 a 1.63x1017N.m; el radio de
fractura entre 0.19 a 4.98km; los valores de caída de esfuerzos están comprendidos
entre 0.13 y 98.76Mpa, la Magnitud Momento Sísmico presenta valores entre 2.36
y 5.48Mw y la Magnitud Energía, valores entre 3.27 y 4.47Me. Los valores de
caída de esfuerzos para la mayoría de réplicas del terremoto de Arequipa analizadas
en este estudio se encuentran entre 0.1 a 100Mpa. Estos valores son coherentes con
los parámetros de la fuente propuestos por otros investigadores.
Asimismo, indicar que todos los eventos sísmicos presentan similares valores de
caída de esfuerzos.
The obtained results present that the aftershocks have values for the seismic
moment scalar in the range of 3.45x1012 to 1.63x1017N.m; the fracture’s radius
between 0.19 to 4.98km; the values of stress drop are between 0.13 and
98.76Mpa, the Seismic Magnitude Moment presents values between 2.36 and
5.48Mw and the Energy Magnitude, values between 3.27 and 4.47Me. The values
of stress drop for the most of aftershocks of Arequipa's earthquake analyzed in
this study are between 0.1 to 100Mpa. These values are coherent with the
parameters of the source proposed by other investigators.
Likewise, to indicate that all the seismic events present similar values of stress
drop.
INTRODUCCIÓN
1
superada la resistencia del material. La caída de esfuerzos (Δσ) es un parámetro
importante de la fuente sísmica ya que representa la parte de los esfuerzos que
produce la fractura y la radiación de energía en forma de ondas elásticas. La
mayor o menor intensidad de la energía dependerá de las dimensiones del área
afectada (A), de la rigidez del medio (µ) y del desplazamiento (ΔU), producto de
la ruptura. El producto de estos parámetros define el momento sísmico,
Mo=µ.A.ΔU.
Por otro lado, Mori y Frankel (1990), utilizando funciones de Green para estimar
los parámetros de sismos con magnitudes comprendidas entre 3.4 y 4.4ML,
ocurridos en Palm Spring (California), indican que la caída de esfuerzos se
encuentra en el rango de 3 a 80bars. De manera similar Hough et al., (1991),
estudiando las réplicas del terremoto de Loma Prieta (1989), determinan que la
caída de esfuerzos se da en el rango de 6 a 266bars. Abercrombie y Leary (1993),
analizaron en detalle la variación de los valores de caída de esfuerzos para 800
sismos de variada magnitud (1.5 a 7ML), ocurridos en California, concluyendo
2
que la caída de esfuerzos varía en un rango entre 0.1 y 100Mpa. En el caso de
Perú, un estudio realizado por Tavera (1998) analizó la caída de esfuerzos de
grandes terremotos con foco superficial, intermedio y profundo ocurridos en el
Perú durante el periodo de 1980 a 1995, determinando que todos los terremotos
con foco a profundidad superficial presentan esfuerzos medios aparentes entre
0.1Mpa y 10Mpa, similar a los definidos en estudios anteriores.
En este estudio se evalúa los parámetros de la fuente de las réplicas del terremoto
de Arequipa del 2001, obtenidos aplicando la técnica del Análisis espectral.
El desarrollo de esta investigación se ha estructurado en ocho capítulos:
3
El Capítulo VI está dedicado al análisis de las señales sísmicas, correspondientes
a las réplicas, correcciones instrumentales, tipos de ruidos y filtros, y la técnica
del Análisis Espectral.
4
Figura 1.1. Ubicación del área de estudio, epicentro del terremoto del
2001 y su serie de réplicas.
5
1.2. JUSTIFICACIÓN
En el Perú los sismos en su mayoría se originan por el proceso de subducción
de la placa de Nazca bajo la placa Sudamericana siendo la fuente sismogénica
primaria que da origen a la mayor cantidad de sismos así como los de mayor
magnitud, el terremoto de Arequipa del año 2001 de magnitud 8.2Mw se
originó en esta fuente, durante el periodo del 27 de junio al 31 de julio del
2001 el Instituto Geofísico del Perú registró 292 réplicas las cuales definen un
área de ruptura de 300x120km2, de aquí la importancia de determinar los
parámetros de la fuente sísmica de las réplicas del terremoto de Arequipa del
2001.
1.4. OBJETIVOS
1.4.1. Objetivo general
Analizar y evaluar las características espectrales de las réplicas del
terremoto de Arequipa del 2001.
6
1.5. METODOLOGÍA
La metodología empleada en esta investigación se resume en los siguientes
pasos:
Recopilación de la información bibliográfica del terremoto de Arequipa
del 2001, y de los registros sísmicos de las réplicas.
Procesamiento de las señales sísmicas en el programa Seismic Analysis
Code (SAC), elaboración de espectros sísmicos, aplicación del análisis
espectral.
Determinación de los parámetros de la fuente sísmica de las réplicas del
terremoto de Arequipa.
Análisis e interpretación de resultados.
1.6. ANTECEDENTES
En el caso de Perú, un estudio realizado por Tavera (1998) analizó la caída de
esfuerzos de grandes terremotos con foco superficial, intermedio y profundo
ocurridos en el Perú durante el periodo de 1980 a 1995, determinando que
todos los terremotos con foco a profundidad superficial presentan caída de
esfuerzos entre 0.7 y 10Mpa.
Por otro lado, Mori y Frankel (1990), utilizando funciones de Green para
estimar los parámetros de sismos con magnitudes comprendidas entre 3.4 y
4.4ML, ocurridos en Palm Spring (California), indican que la caída de
7
esfuerzos se encuentra en el rango de 3 a 80bars. Hough et al., (1991),
estudiando las réplicas del terremoto de Loma Prieta (1989), determinan que
la caída de esfuerzos se da en el rango de 6 a 266bars.
8
CAPÍTULO II
GEODINÁMICA
- La Dorsal de Nazca
- La Fractura de Mendaña
- La Fosa Peruano-Chilena
- La Cordillera Andina
- La Cadena Volcánica
- Sistemas de Fallas
9
A continuación, se describe las características principales de cada uno de estos
elementos.
10
2.1.4. La Cordillera Andina, la Cordillera Andina se formó a finales del Cretácico
tardío, por el proceso de subducción de la placa de Nazca bajo la Placa
Sudamericana. Esta cadena montañosa se extiende a lo largo del borde
occidental de Sudamérica y tiene una extensión de 7500km aproximadamente
con alturas máximas del orden de 6000msnm. En Perú de Norte a Sur la
Cordillera de los Andes está bien definida y muestra dos deflexiones:
Deflexión de Huancabamba, a la altura de 5º Sur y la Deflexión de Abancay,
a la altura de 14º Sur, estas deflexiones cambian la orientación de la
Cordillera Andina en dirección Noreste-Suroeste y Este-Oeste
respectivamente. Desde el punto de vista geomorfológico la Cordillera
Andina está constituida de una serie de unidades morfoestructurales que se
distribuyen de Oeste a Este y que se ilustran en la Figura 2.1: La Zona
Costera, la Cordillera Occidental, el Altiplano, la Cordillera Oriental y la
Zona Subandina (Dalmayarac et al., 1988; Tavera y Buforn, 1998). Estas
unidades son descritas a continuación:
Zona Costera, corresponde a una estrecha franja que se extiende a lo largo
del litoral peruano, con elevaciones que varían entre 50 y 1500msnm y ancho
entre 15 y 100km en las regiones Norte y Centro y aproximadamente 50 km
en la región Sur. La zona Norte y Centro presenta suaves plegamientos de
rocas volcánicas y sedimentarias, la región Sur se caracteriza por presentar
fuertes plegamientos de rocas basálticas.
11
Figura 2.1. Principales unidades geomorfológicas presentes en el Perú (Dalmayrac et al.,
1988; Tavera y Buforn, 1998).
12
La Cordillera Oriental, morfológicamente corresponde a un extenso
anticlinal formado por depósitos de tipo intrusivo, cuyo levantamiento fue
controlado por una serie de fallas regionales. Esta cordillera se extiende de
Norte a Sur siguiendo una orientación NO-SE, con elevaciones que varían
entre 3700 a 4000msnm y anchos entre 70 a 100km aproximadamente. En la
región Sur soporta un fuerte arqueamiento en dirección E-O a la altura de la
deflexión de Abancay; mientras que en la región Norte, desaparece debido a
la presencia de la deflexión de Huancabamba.
2.1.5. La Cadena Volcánica, se ubica en la región Sur del Perú por debajo de la
deflexión de Abancay entre los 14º hasta los 25º Sur en Chile. Esta cadena se
distribuye sobre la cordillera Occidental, siguiendo en el Perú un aparente
alineamiento con orientación NO-SE. Las características geométricas de cada
uno de los volcanes que integran esta cadena, muestran que la actividad
tectónica es contemporánea a la orogenia extensional que experimenta la
Cordillera Andina cerca del Cuaternario Medio y Reciente (Sebrier et al.,
1985).
13
describen las características de los principales sistemas de fallas según su
ubicación:
Zona Costera, en esta zona, al Norte de 5º Sur se ubica la Falla de Huaypira
con orientación NE-SO y E-O. Entre 14º y 16º Sur, se encuentra la Falla de
Marcona ubicada al NE de San Juan de Marcona (Ica). A la altura de 16.5º S
destaca la Falla de la Planchada en la provincia de Islay (Arequipa). Estas
fallas presentan longitudes del orden de 90km en promedio y son de tipo
normal.
Cordillera Oriental y Occidental, el importante proceso de levantamiento que
ha soportado la Cordillera Andina, ha permitido la formación del sistema de
Fallas de la Cordillera Blanca, siendo su ramal Norte conocido como Falla de
Quiches. Ambas fallas presentan un buzamiento SO.
El sistema de Fallas de Huambo-Cabanaconde compuesto por fallas Trigal y
Solarpampa, que se ubican en el extremo Norte del complejo volcánico
Ampato-Sabancaya Hualca-Hualca entre las localidades de Huambo y
Cabanaconde. Estas fallas son de tipo normal con rumbo en dirección E-O y
buzamiento hacia el Sur con valores entre 60° y 70° (Sébrier et al., 1985) y
presentan longitudes del orden de 10 a 12km aproximadamente.
La Falla de Ichupampa se ubica en el extremo NE del complejo volcánico
Ampato-Sabancaya-Hualca Hualca, es de tipo normal con rumbo en dirección
NO-SE y buzamiento hacia el Este. La longitud aproximada de esta falla es
de 25km.
En la latitud -12.5º se ubica la Falla de Huaytapallana al pie del nevado de la
Cordillera de Huaytapallana (Huancayo), presenta movimiento de tipo
inverso en dirección NO-SE con buzamiento de 65º hacia el NE. Esta falla se
reactivó con los sismos de 1969.
En la región Sur del Perú se encuentra el sistema de fallas de Tambomachay
ubicada entre 13º y 14.5º Sur, al Norte de la ciudad del Cusco, está compuesta
por las fallas de Tambomachay, Tamboray, Koricocha y Cusco.
En la Zona Subandina destacan los sistemas de Rioja, Satipo y Madre de Dios
con longitudes de 300-500km y en general, se orientan paralelas a la
Cordillera Andina con buzamiento hacia el SO. En la Figura 2.2 se presenta
el mapa de ubicación de los principales sistemas de fallas (Macharé et al.,
2003).
14
Figura 2.2. Mapa de ubicación de los principales sistemas de fallas del Perú (Macharé et al.,
2003).
15
CAPÍTULO III
El historiador José Toribio Polo (1904), estimó que se habían producido más
de 2500 temblores en territorio peruano desde la conquista hasta fines del
siglo XIX y advirtió que por varias causas no se anotaron muchos sismos
durante el periodo de 1600 a 1700. En ese periodo varias ciudades del Perú
fueron destruidas por sismos, Arequipa en el año 1582, 1600 y 1784; Cusco
en 1650; Trujillo en 1619 y 1725. En la región Central sobresalen los sismos
ocurridos en 1586, primer gran sismo para el cual se tiene documentación
histórica; 1687 y 1746, los mismos que destruyeron casi completamente a la
ciudad de Lima siendo la cantidad de pérdidas humanas inmensa (Silgado,
1978). Los daños materiales fueron cuantiosos debido a que las
construcciones eran inadecuadas, se construían con adobe y quincha en la
costa y con piedra en las regiones altas. Durante el siglo XIX sucedieron
varios sismos, uno de los principales por su intensidad, fue el de 1868 que
devastó Arequipa, Tacna y Arica, para luego ser seguido por un tsunami que
afectó a todo el Océano Pacífico llegando a playas de Japón, Nueva Zelanda y
Australia.
16
Según Silgado (1978), en el presente siglo fueron notables los sismos
ocurridos en Piura-Huancabamba (1912), Caravelí (1913), Chachapoyas
(1928), Lima (1940), Nazca (1942), Quiches-Ancash (1946), Satipo (1947),
Cusco (1950), Tumbes (1953), Arequipa (1958-1960), Lima (1966),
Chimbote-Callejón de Huaylas (1970) y Lima (1974).
Los últimos estudios realizados por Tavera y Buforn (2001), Bernal et al,
(2001) y Tavera y Bernal (2009), usando datos telesísmicos para el periodo
1964-1996 con mb ≥ 5.0, concluyen que la distribución espacial de los sismos
en el Perú es compleja; no obstante, se pueden identificar áreas de sismicidad
constituidas por importantes fuentes sismogénicas. En el borde Occidental del
Perú durante el periodo de 1960-2010 han ocurrido siete terremotos de
magnitudes elevadas en las regiones Centro y Sur (Figura 3.1).
17
de los más destructores, perecieron 67000 personas y 150000 heridos
(Ericksen et al., 1970). El terremoto provocó un alud como consecuencia del
desprendimiento de una parte del nevado Huascarán y sepultó las localidades
de Yungay y Ranrahirca. La máxima intensidad fue de IX (MM).
- El 23 de junio del 2001, ocurrió un sismo que afectó todo el Sur del Perú.
El terremoto produjo un tsunami con olas de 7-8m de altura. La máxima
intensidad evaluada fue de VIII (MM) y su magnitud de 8.2Mw. Este
sismo será descrito con detalle en capítulos posteriores.
18
Figura 3.1. Distribución espacial de grandes terremotos ocurridos en el borde occidental del
Perú, periodo instrumental (1960-2010). Las barras de color azul indican los sismos que
produjeron tsunami.
19
3.3. PRINCIPALES SISMOS DESTRUCTORES EN EL SUR DEL PERÚ
A continuación se describen las principales características de los sismos más
importantes ocurridos en la región Sur del Perú.
20
• Sismo del 24/08/1942
Sismo que destruyó los alrededores de Ica y Arequipa. Se produjo un
tsunami con olas que alcanzaron 3 metros de altura. Perecieron 30
personas y 25 heridos por diversas causas. La máxima intensidad evaluada
de este sismo fue de IX (MM) y su magnitud fue de 6.7Ms.
21
y la máxima intensidad evaluada VIII (MM).
22
Figura 3.2. Mapa sísmico de Perú para el periodo 1960-2017. Fuente Instituto Geofísico del Perú.
23
toda la región Sur de Perú y se caracteriza por presentar un mayor número
de sismos a una distancia de 500km de longitud desde la fosa Perú-Chile,
alcanzan profundidades de 300km aproximadamente.
24
del magma desde la placa subducente, que produce el engrosamiento de la
corteza y el plegamiento de los sedimentos hasta formar una franja de
cadenas montañosas con volcanismo activo paralelo al margen de la placa.
25
CAPÍTULO IV
El sábado 23 de junio del 2001 a las 15 horas con 33 minutos y 14.2 segundos
(Hora local) ocurrió en la zona Sur del Perú un terremoto de elevada magnitud
(Mw=8.2) con epicentro a 82km al NO del distrito de Ocoña, provincia de
Camaná en el departamento de Arequipa.
26
82km al NO de la localidad de Ocoña y foco a una profundidad focal de
29km. Este evento se caracterizó por su gran duración, puesta de manifiesto
en la complejidad de su registro que inicialmente podría evidenciar un
proceso de ruptura muy heterogéneo, así como el modo de propagación de la
energía que produjo el ondulamiento de la superficie y la generación de un
tsunami de carácter netamente local. Este sismo produjo en las primeras
veinticuatro horas 134 réplicas (ML 3), todas con sus epicentros ubicados
al SE del epicentro del evento principal, lo cual sugiere la propagación de la
ruptura en esa dirección.
Figura 4.1. Mapa de isosistas del terremoto de Arequipa (Tavera et al., 2002).
27
Las isosistas muestran una distribución elíptica. La máxima intensidad
evaluada fue de VIII (MM) en las localidades de Ocoña, Camaná y Mollendo.
La isosista IV tiene el mayor eje (1150km) orientado en dirección NO-SE
paralela a la línea de costa. En la ciudad de Arequipa muchas iglesias y
monumentos históricos fueron dañados. En los departamentos de Moquegua y
Tacna el mayor daño corresponde a las edificaciones antiguas de adobe y
quincha, en distritos que han ocupado zonas con suelos con características
geotécnicas inestables.
En la carretera Panamericana Sur aparecieron varias grietas con orientación
paralela a la línea de costa y ocurrieron varios deslizamientos. En continente
el sismo fue sentido en Cusco (Perú) y en La Paz (Bolivia), evaluándose
intensidades de III (MM) para ambas ciudades. Asimismo, la DIREMER
Región I de Chile informó que el terremoto de Arequipa produjo en la ciudad
de Arica intensidades del orden de VII, en Iquique de VI y en Calama de IV-
V en la escala Mercalli Modificada.
Los parámetros hipocentrales del evento sísmico del 23 de Junio y sus 220
réplicas han sido relocalizadas con información de la Red Sísmica Nacional
del Perú y de 8 estaciones instaladas después de ocurrido el sismo, siendo
todos los registros de estaciones de banda ancha de tres componentes. La
relocalización de réplicas se realizó teniendo como mínimo 4 estaciones de
registro y el algoritmo Hypoinverse (Klein, 1978), se usó el modelo de
velocidad de Lindo et al. (1993) y la relación Vp/Vs=1.78 calculado con casi
800 réplicas.
28
ocurrencia de dos réplicas de gran magnitud (Mw=6.6 y 7.5), definiendo la
presencia de una importante barrera.
Figura 4.2. Distribución espacial del epicentro del evento principal y réplicas relocalizadas
del terremoto del 23 de Junio del 2001. Los mecanismos focales fueron obtenidos usando
registros telesísmicos. El área rayada corresponde a las réplicas del terremoto de Noviembre
de 1996 (Tavera et al., 1998) y la sombreada, al área que no ha experimentado ruptura
durante el terremoto del 23 de Junio. Los grupos de réplicas están delimitados por líneas
punteadas. La fecha y magnitud están ubicadas encima de cada mecanismo focal
(Redibujado de Tavera et al., 2006).
29
producido en un medio asísmico (Área sombreada de la Figura 4.2). Ambas
hipótesis son válidas cuando se considera que el terremoto del 23 de Junio del
2001 tuvo un proceso de ruptura complejo a través de una superficie muy
heterogénea (Tavera et al., 2006).
30
nodales orientados hacia la misma dirección del evento principal, sin embargo
los grupos G12 y G13 tienen una pequeña componente rumbo-buzamiento.
En general el patrón de deformación sobre toda el área de réplicas fue
aproximadamente el mismo (Tavera et al., 2002).
El mecanismo focal del terremoto principal y de tres réplicas también fueron
evaluadas a partir del modelado de las amplitudes y formas de la onda P
registradas en estaciones sísmicas localizadas a distancias telesísmicas
(30°d 90°). Tavera et al. (2002) hace uso de la metodología propuesta por
Nabeleck (1984) con la finalidad de cuantificar la orientación de la fuente
(mecanismo focal), la profundidad del foco (h), la función temporal de la
fuente sísmica (STF) y el momento sísmico escalar (Mo). Mediante un
proceso iterativo se hace mínima la desviación estándar de los residuos de la
correlación existente entre los registros observados y teóricos, siendo este
evaluado a partir de la “raíz media cuadrática (RMS)” del ajuste.
31
El momento sísmico para el sismo principal es de 1.3x1021N.m
equivalente a una magnitud de 8.1Mw.
32
La STF es simple con una duración de aproximadamente 5 segundos y
un momento sísmico de 1.19x1019N.m equivalente a una magnitud de
6.6Mw.
33
Figura 4.5. Resultados del modelamiento de ondas P para las réplicas del terremoto
de Arequipa a) 26 de Junio, b) 5 de Julio y c) 7 de Julio (Redibujado de Tavera et al.,
2002).
34
4.4. TSUNAMI
El terremoto de Arequipa del 2001 fue seguido por un tsunami que se inició
15 minutos después de ocurrido el sismo principal, con el replegamiento del
mar, principalmente entre las playas de las localidades de Chala e Ilo. Por
ejemplo, en las playas de las localidades de Atico y Ocoña, el mar se replegó
hasta 200 metros y en su retorno inundó las playas entre 30 a 50 metros
principalmente por la desembocadura del río Ocoña.
Los datos obtenidos indican que los volúmenes de agua que siguieron al
tsunami entre las localidades de Chala, Atico y Ocoña formaron corrientes de
agua que discurrieron en dirección SE; mientras que en las localidades de
Matarani y Mollendo las olas siguieron una dirección NO. Esta información
indica que el tsunami se propagó de forma casi perpendicular a la línea de
costa.
35
Figura 4.6. Principales parámetros y características del tsunami que afectó la localidad de
Camaná en Arequipa. Las líneas perpendiculares a la línea de Costa indican la máxima altura
de las olas. El área sombreada corresponde al área inundada por el tsunami (Redibujado de
Tavera et al., 2006).
36
CAPÍTULO V
37
barreras o zonas de alta resistencia. El deslizamiento se produce
únicamente en dos puntos de la falla, son las áreas sin sombrear, porque
las barreras reducen considerablemente el movimiento de sus
compartimientos. En general, el modelo de Barreras describe mejor los
terremotos asociados a deformaciones corticales, fallas que se
caracterizan por producir terremotos de magnitud moderada (Tavera,
1992).
38
El nuevo terremoto iniciará en una de estas áreas o en cadena.
Explicando el comportamiento de los esfuerzos, casi la totalidad de la
falla ya liberó energía y los esfuerzos se encuentran en las zonas
resistentes. Durante los terremotos posteriores, la energía se liberará
totalmente y la situación final será de esfuerzos casi uniformes. El
deslizamiento producido es similar al Modelo de Barreras, en las áreas
con caída de esfuerzos apenas si hay deslizamiento (Tavera, 1992). El
modelo de Asperezas se aplica a los grandes terremotos ocurridos en el
borde occidental del Perú, producto del proceso de subducción.
39
instantáneamente sobre toda el área de la falla, el pulso genera una onda
de cizalla que se propaga perpendicularmente al plano de falla (Brune,
1970).
x
∆σ(x,t) =∆σH (t- ) (5.1)
β
En la ecuación 5.1, Δσ corresponde al esfuerzo de cizalla efectivo, ΔH
es la transformada de la función temporal de la fuente, x es la distancia
normal al plano de falla y β la velocidad de la onda S. El
desplazamiento de cizalla Δu, para x=0, es obtenido por integración de
la ecuación 5.1 teniendo en cuenta que σ = µ𝜕u/𝜕x es:
∆σ
∆u(t) = H(t) βt (5.2)
μ
∆σβ
∆U(ω)= - (5.3)
μω2
El espectro es:
∆σβ 1
U(ω)= (5.4)
μ ω2+b2
40
2.34β
b= (5.5)
a
Donde a es el radio de la falla. El espectro tiene una parte plana y una
caída de frecuencias, entonces la frecuencia esquina 𝜔𝑐 = b .
Figura 5.3. Sismogramas y espectros para cinco réplicas del terremoto de San
Fernando (1971). Las flechas debajo de los sismogramas indican la parte de la onda S
utilizada para el espectro de Fourier (Brune, 1973).
a = 0.37β/ωc (5.6)
41
5.2.2. Modelo de Madariaga
Este modelo considera una ruptura nucleada en un punto y se expande a
una velocidad constante (Madariaga, 1976). La falla detiene su
crecimiento instantáneamente en el tiempo:
tr=a/ VR (5.7)
Donde:
tr = tiempo de ruptura
a = radio
VR = velocidad de ruptura
42
resolver numéricamente en coordenadas cilíndricas. En este modelo se
supone que la ruptura se propaga con una velocidad constante v, hasta
alcanzar una barrera situada a un radio finito a.
Ω 4π ρ α3R
Mo = (5.8)
Rθφ F
Donde:
43
α = velocidad de la onda S
R = distancia hipocentral
r = 0.37α/fc (5.9)
Donde:
r = radio de una fuente circular
α = velocidad de la onda S
fc = frecuencia esquina
44
5.3.3. Caída de esfuerzos
En un modelo de fractura simplificada, el desplazamiento relativo Δ
de los lados de una falla es debido a los esfuerzos de cizalla, los cuales
en un momento dado exceden la resistencia del material o la fricción
que mantiene la falla.
Si los esfuerzos de cizalla que actúan en el plano de falla antes y
después de un terremoto son σ0 y σ1, se pueden definir dos nuevos
__
conceptos llamados, el “esfuerzo medio 𝜎 ” que es el valor promedio de
la suma de los esfuerzos actuantes antes y después de un terremoto y la
“caída de esfuerzos Δσ”, la cual se define como la diferencia entre el
esfuerzo inicial (σ0) y esfuerzo final (σ1); es decir, antes y después de la
ruptura (Figura 5.5).
1
Esfuerzo Medio σ = 2 (σ0+σ1) (5.10)
45
El esfuerzo inicial σ0 es el esfuerzo tectónico responsable de la
deformación en la región focal. Conocido Δσ, otro parámetro
sismológico que es posible estimar es la “energía total liberada” y para
ello se utiliza la siguiente relación:
E= σ ∆µ S (5.12)
46
5.3.4. Energía sísmica
Cuando se produce un terremoto la cantidad de energía liberada es
difícil de calcular con precisión, la mayor parte de dicha energía se
acumula en la roca y al disiparse lo hace en forma térmica y en forma
de ondas sísmicas. La energía sísmica es la única que puede ser
estimada a partir del registro de los sismos; por lo tanto, se asume que
el tamaño de los sismos puede ser conocido a partir de la cantidad de
energía liberada por ellos. En sismología, el tamaño de los sismos es
estimado a partir de las amplitudes máximas del registro de las ondas de
volumen ó superficiales expresados en valores para la magnitud mb ó
Ms, tal como lo propone Gutenberg y Richter (1956):
Tabla 5.1. Magnitudes en comparación con la energía sísmica liberada y la cantidad de TNT.
Magnitud Energía Equivalencia Referencias
Local (ergios) (TNT)
-1.5 1.58E+02 1g Rotura de roca en laboratorio
1.0 1.58E+08 170 g Explosión en construcción
3.5 1.58E+14 455 kg Explosión en mina
4.5 3.98E+16 32 t Tornado promedio
5.5 1.00E+19 500 t Terremoto de Cusco, 1986
6.0 1.58E+20 1270 t Terremoto de Lima, 1993
6.5 2.51E+21 31550 t Terremoto de Moyobamba,1991
7.0 3.98E+22 199000 t Terremoto de Chimbote, 1996
7.5 6.31E+23 1000000 t Terremoto de Nazca, 1996
8.0 1.00E+25 6270000 t Terremoto de Lima, 1940
8.5 1.58E+26 31550000 t Terremoto de Lima, 1746
9.0 2.51E+27 19999000 t Terremoto de Arequipa, 1868
10.0 6.31E+29 6.3 billón t Terremoto de Chile, 1960
47
La energía liberada por la bomba atómica es equivalente a un sismo
ocurrido en superficie con magnitud 6mb. Choy y Boatwright en 1995
definieron la magnitud energía con la siguiente relación:
48
CAPÍTULO VI
Figura 6.1. Ejemplo del registro de un sismo en el campo cercano. P y S, corresponden a los
tiempos de llegada de ambas fases.
49
que la lectura e interpretación de los mismos se convierte en una tarea
compleja.
Figura 6.2.a) Registro de velocidad sin corregir para una réplica del terremoto de Arequipa
del 2001, obsérvese que el registro se encuentra por debajo de la línea base; b) Registro de
velocidad, al cual se aplicó la corrección de línea base.
50
F(t)= función que caracteriza a la fuente
Q(t)= perturbación del medio o atenuación
I(t)= función de respuesta del instrumento
S(t)=F(t)*Q(t)*I(t) (6.1)
Figura 6.3. Representación del cálculo de un sismograma teórico para una fuente
puntual en el campo lejano. La función temporal de la fuente sísmica (F(t)) se
combina con la atenuación del medio (Q(t)) y el instrumento (I(t)) mediante una
operación de convolución (Tavera, 1998).
S(t)
= F(t)*Q(t) (6.2)
I(t)
S(t)
S(t) = real (IFF ) = S(t) = F(t)*Q(t) (6.3)
I(t)
51
Figura 6.4. Ficha técnica de un sismómetro Guralp 40T.
52
En la figura 6.5 se muestra el registro de una réplica del sismo de
Arequipa, registrada el 01/07/2001 en la estación Caravelí, sin
corrección instrumental 6.5.a, la amplitud del registro sísmico está en
cuentas y con corrección instrumental 6.5.b sus unidades están en cm/s.
Figura 6.5.a) Registro de la réplica del sismo de Arequipa del 2001 sin corrección
instrumental, la amplitud está expresada en cuentas/s. b) El mismo registro, al cual
se aplicó la corrección instrumental, la amplitud expresada en cm/s.
53
distancia y profundidad. El ruido cultural muestra fuertes variaciones
diurnas y tiene frecuencias características que dependen de la fuente de
disturbio (Cutipa y Tavera, 2006).
Figura 6.6. Registros de réplicas del sismo de Arequipa con diferentes niveles de ruido
sísmico (rectángulo sombreado).
54
temperatura durante el día y la noche. Se debe considerar como ruidos
de periodo corto, a los producidos por las explosiones en las minas, el
paso de camiones en las carreteras, o por la inducción del viento sobre
la estación sísmica con rangos de frecuencia mayores a 10Hz. Este tipo
de ruido suele ser importante si en las cercanías de la estación sísmica
hay vegetación o montañas (Cutipa y Tavera, 2006).
6.3.4 Filtros
A fin de eliminar el ruido de fondo presente en los registros sísmicos es
necesario aplicar filtros.
55
de entrada en el tiempo y amplitud, realiza un procesamiento
matemático sobre dicha señal, mediante el uso de la Transformada
rápida de Fourier.
Un filtro ideal tiene una banda pasante de corte totalmente plano y unas
zonas de transición entre ambas nulas, en la práctica esto nunca se
consigue (Figura 6.7).
𝑥𝑡 → ℎ0 → 𝑦𝑡 = ℎ0 𝑥𝑡 (6.4)
56
𝑥𝑡 → ℎ1𝑧 → 𝑦𝑡 = ℎ1𝑥𝑡−1 (6.5)
57
b) Filtro Paso-Bajo, este tipo de filtro permite el paso de las
frecuencias más bajas y atenuar las frecuencias altas (Figura
6.8).
58
Figura 6.11. Respuesta de un filtro Elimina-banda.
1 +∞
f(t)= ∫ F(ω)ejωt.dω (6.9)
2π -∞
+∞
F(ω)= ∫ f(t)e-jωt.dt (6.10)
-∞
59
señales a procesar son discretas, con un número finito de valores, entonces es
necesario modificar las ecuaciones 6.7 y 6.8 para adaptarlas a la naturaleza de
la información a procesar. De este modo, surge la Transformada de Fourier
Discreta (DFT), (Figura 6.12), la cual hace posible aplicar Fourier en el
análisis de señales digitales.
60
Es importante resaltar que la DFT reduce considerablemente los errores de
redondeo, asociados a la representación de datos numéricos en una
computadora.
Figura 6.13. Espectros para dos réplicas del terremoto de Arequipa, la primera registrada el
29 de Junio en la estación Tocota y la segunda registrada el 01 de Julio en la estación
Caravelí. En el eje de las abscisas, la frecuencia está en Hz y en el eje de ordenadas, las
amplitudes en cm.se, Ω es la parte plana del espectro y fc la frecuencia esquina.
61
Mo = 4π.ρ.α3.Ω.R / RθφF (6.11)
r =0.37α/fc (6.12)
62
CAPÍTULO VII
63
Camaná, zona donde se produjo un tsunami con olas de 8 metros de altura
aproximadamente.
Figura 7.1. Epicentros de la serie de réplicas del terremoto de Arequipa del 2001.
64
Bella Unión y Tocota fueron trasladadas e instaladas en las localidades de
Pocoma y Aplao respectivamente. Todas las estaciones sísmicas son de banda
ancha, con excepción de la estación de San Gregorio que funciona con un
sismómetro de periodo corto. Además se usaron registros de tres estaciones
adicionales de la Red Sísmica Nacional (RSN) (Conima, Huancayo y
Toquepala). En la Figura 7.2, se presenta su distribución espacial.
Figura 7.2. Distribución de estaciones sísmicas que integran la Red Sísmica Local. La estrella
indica el epicentro del sismo principal.
65
Tabla 7.1. Estaciones sísmicas de banda ancha que integran la red sísmica local.
66
7.4. METODOLOGÍA PARA LA ESTIMACIÓN DE LOS PARÁMETROS
DE LA FUENTE SÍSMICA
Para determinar los parámetros de la fuente sísmica de las réplicas se han
seguido los siguientes pasos:
1. Identificar, seleccionar y cortar los registros de las réplicas del Terremoto
de Arequipa del 2001, utilizando las componentes horizontales.
2. Aplicar la corrección de línea base.
3. Aplicación de corrección instrumental a los registros seleccionados.
4. Para eliminar el efecto de las frecuencias bajas en la señal sísmica, se
aplica un filtro pasa-altos Butterworth n2 cos 0.3. Integrar la señal para
disponer del registro en desplazamiento.
5. Calcular el espectro de desplazamiento de la onda S, aplicando la
Transformada Rápida de Fourier.
6. En el espectro de desplazamiento se aplica el Análisis espectral, se
identifica la parte plana del espectro (Ω) y la frecuencia esquina (fc) para
cada espectro (Componentes horizontales) y se promedian.
7. Se determinan los parámetros de la fuente sísmica, para calcular el
Momento sísmico escalar (Mo) se usa la relación propuesta por Keilis-
Borok (1960) y el modelo de velocidad propuesto por Grange et al.,
(1984), (Tabla 7.2). Calculado el Mo, se estima el radio de ruptura y luego
la caída de esfuerzos (Δσ). Conocido el Mo, la magnitud Mw se estima
con la relación Mw=2/3 Log(Mo)-6, (Kanamori y Anderson, 1975). La
Magnitud Energía es estimada en función de la energía radiada por el
sismo usando la relación Me = 2/3 Log Es-2.9, (Choy y Boatwright, 1995).
67
8. Elaboración de los diagramas de caída de esfuerzos (Δσ) para las
estaciones Tocota, Ocoña, Caravelí, Bella Unión y Quilca.
Todo el procesamiento de las señales sísmicas se realizó en el programa
Seismic Analysis Code (SAC) en el entorno de Linux.
7.5. RESULTADOS
Para la presentación y descripción de los resultados obtenidos en esta
investigación se realiza estación por estación:
68
Figura 7.3. Ejemplos de espectros de ondas S para el desplazamiento del suelo correspondiente a
réplicas registradas en la estación Tocota, N-S (Componente Norte- Sur) y E-O (Componente Este-
Oeste).
69
Figura 7.4. Ejemplos de espectros de ondas S para el desplazamiento del suelo correspondiente a
réplicas registradas en la estación Ocoña, N-S (Componente Norte- Sur) y E-O (Componente Este-
Oeste).
70
7.5.3. Estación Caravelí
La estación Caravelí ubicada en la provincia de Caravelí, departamento
de Arequipa, registró 191 réplicas se muestran dos espectros
representativos en la Figura 7.5. Los resultados presentan las siguientes
características espectrales:
• La parte plana del espectro (Ω), varía en un rango comprendido entre
6e-6 y 3.37e-2cm.s.
• La frecuencia esquina (fc) varía entre 0.40 y 6.75Hz.
71
Figura 7.5. Ejemplos de espectros de ondas S para el desplazamiento del suelo correspondiente a
réplicas registradas en la estación Caravelí, N-S (Componente Norte-Sur) y E-O (Componente Este-
Oeste).
72
Figura 7.6. Ejemplos de espectros de ondas S para el desplazamiento del suelo correspondiente a
réplicas registradas en la estación Bella Unión, N-S (Componente Norte-Sur) y E-O (Componente
Este-Oeste).
73
Figura 7.7. Ejemplos de espectros de ondas S para el desplazamiento del suelo correspondiente a
réplicas registradas en la estación Quilca, N-S (Componente Norte- Sur) y E-O (Componente Este-
Oeste).
74
En la Tabla 7.3 se presentan la lista de las 235 réplicas registradas en la
estación Tocota, para las cuales se ha elaborado los espectros de
desplazamiento, también se presenta sus parámetros hipocentrales como
epicentros, profundidad del foco y Magnitud Local, además los
parámetros evaluados como momento sísmico, radio, caída de
esfuerzos, Magnitud Momento (Mw) y Magnitud Energía (Me). Las
tablas con los parámetros de la fuente sísmica obtenidos para las
estaciones Ocoña, Caravelí, Bella Unión y Quilca se presentan en el
Anexo 1, en el Anexo 2 se consideran treinta espectros de
desplazamiento representativos para cada una de las estaciones
sísmicas.
75
Tabla 7.3. Continuación
76
Tabla 7.3. Continuación
77
Tabla 7.3. Continuación
78
Tabla 7.3. Continuación
79
Tabla 7.3. Continuación
Figura 7.8. Relación entre el logaritmo de Mo y r para las réplicas del terremoto de
Arequipa, registradas en la estación Tocota.
80
En la Figura 7.8 se ha representado el esfuerzo medio aparente con las
líneas de 0.1, 1, 10 y 100MPa (Kanamori y Anderson, 1975). Se
observa que el mayor número de réplicas presentan valores de caída de
esfuerzos entre 0.1 y 100Mpa, similares a valores obtenidos para
réplicas y microsismicidad (Abercrombie y Leary, 1993). Para la
estación Tocota, la variación de los valores de Momento sísmico escalar
(Mo), radio de fractura, caída de esfuerzos, Magnitud momento (Mw) y
Magnitud energía (Me) se presentan en la Tabla 7.4. Se consideran los
dos decimales para analizar la variación de los parámetros de la fuente
con detalle.
Tabla 7.4. Rangos de valores de momento sísmico, radio de fractura, caída de esfuerzos,
Magnitud momento y Magnitud energía para réplicas registradas en la estación Tocota.
Momento sísmico Radio de fractura Caída de esfuerzos Magnitud Momento Magnitud Energía
(N.m) (m) (Mpa) (Mw) (Me)
3.45x1012 - 1.63x1017 190-4122 0.13-72.59 2.36-5.48 3.27-4.47
Figura 7.9. Relación entre el logaritmo de Mo y r para las réplicas del terremoto de
Arequipa, registradas en la estación Ocoña.
81
En la Figura 7.9 se ha representado el esfuerzo medio aparente con las
líneas de 0.1, 1, 10 y 100MPa (Kanamori y Anderson, 1975). Se
observa que el mayor número de réplicas presentan valores de caída de
esfuerzos entre 0.1 y 100Mpa, similares a valores obtenidos para
réplicas, emjambres y microsismicidad (Abercrombie y Leary, 1993).
Para la estación Ocoña, la variación de los valores de Momento sísmico
escalar (Mo), radio de fractura, caída de esfuerzos, Magnitud momento
(Mw) y Magnitud energía (Me) se presentan en la Tabla 7.5.
Momento sísmico Radio de fractura Caída de esfuerzos Magnitud Momento Magnitud Energía
(N.m) (m) (Mpa) (Mw) (Me)
1.39x1013 - 1.63x1017 190-2220 0.48-86.65 2.76-5.47 3.43-4.47
Figura 7.10. Relación entre el logaritmo de Mo y r para las réplicas del terremoto de
Arequipa, registradas en la estación Caravelí.
82
En la Figura 7.10 se ha representado el esfuerzo medio aparente con las
líneas de 0.1, 1, 10 y 100MPa (Kanamori y Anderson, 1975). Se
observa que el mayor número de réplicas presentan valores de caída de
esfuerzos entre 0.1 y 100Mpa, similares a valores obtenidos para
réplicas, emjambres y microsismicidad (Abercrombie y Leary, 1993).
Para la estación Caravelí, la variación de los valores de Momento
sísmico escalar (Mo), radio de fractura, caída de esfuerzos, Magnitud
momento (Mw) y Magnitud energía (Me) se presentan en la Tabla 7.6.
Tabla 7.6. Rangos de valores de momento sísmico, radio de fractura, caída de esfuerzos,
Magnitud momento y Magnitud energía para réplicas de la estación Caravelí.
Momento sísmico Radio de fractura Caída de esfuerzos Magnitud Momento Magnitud Energía
(N.m) (m) (Mpa) (Mw) (Me)
1.35x1013 - 1.63x1017 210-3610 0.20-98.76 2.75-5.48 3.43-4.47
Figura 7.11. Relación entre el logaritmo de Mo y r para las réplicas del terremoto de
Arequipa, registradas en la estación Bella Unión.
83
En la Figura 7.11 se ha representado el esfuerzo medio aparente con las
líneas de 0.1, 1, 10 y 100MPa (Kanamori y Anderson, 1975). Se
observa que el mayor número de réplicas presentan valores de caída de
esfuerzos entre 0.1 y 100Mpa, similares a valores obtenidos para
réplicas, emjambres y microsismicidad (Abercrombie y Leary, 1993).
Para la estación Bella Unión, la variación de los valores de Momento
sísmico escalar (Mo), radio de fractura, caída de esfuerzos, Magnitud
momento (Mw) y Magnitud energía (Me) se presentan en la Tabla 7.7.
Tabla 7.7. Rangos de valores de momento sísmico, radio de fractura, caída de esfuerzos,
Magnitud momento y Magnitud energía para réplicas de la estación Bella Unión.
Momento sísmico Radio de fractura Caída de esfuerzos Magnitud Momento Magnitud Energía
(N.m) (m) (Mpa) (Mw) (Me)
1.04x1013 – 7.56x1015 210-580 0.25-90.94 2.68-4.59 3.40-4.14
Figura 7.12. Relación entre el logaritmo de Mo y r para las réplicas del terremoto de
Arequipa, registradas en la estación Quilca.
84
En la Figura 7.12 se ha representado el esfuerzo medio aparente con las
líneas de 0.1, 1, 10 y 100MPa (Kanamori y Anderson, 1975). Se
observa que el mayor número de réplicas presentan valores de caída de
esfuerzos entre 0.1 y 100Mpa.
Para la estación Quilca, la variación de los valores de Momento sísmico
escalar (Mo), radio de fractura, caída de esfuerzos, Magnitud momento
(Mw) y Magnitud energía (Me) se presentan en la Tabla 7.8.
85
CAPÍTULO VIII
ANÁLISIS E INTERPRETACIÓN
8.1. Interpretación
El objetivo principal del presente estudio es determinar los parámetros de la
fuente sísmica de las réplicas del sismo de Arequipa del 2001, aplicando la
técnica del análisis espectral. El paso más complejo en la aplicación de la
técnica de análisis espectral, es la determinación correcta de las
características del espectro sísmico; es decir, la parte plana (Ω) y la
frecuencia esquina (fc), lo cual requiere de experiencia además de un buen
conocimiento del modelo de la fuente sísmica y la escala del espectro
sísmico. A pesar de la complejidad para resolver el problema de la fuente
sísmica, el modelo de Brune para sismos con ML≤6 presenta características
simples. El momento sísmico escalar es uno de los parámetros más
importantes de la fuente sísmica (Aki, 1967), porque determina el tamaño del
sismo, de allí que se pueda definir correctamente parámetros como Mo y Δσ.
La cantidad de esfuerzos liberados durante un terremoto es una característica
fundamental del proceso de ruptura. En concreto, la caída de esfuerzos
representa la cantidad de esfuerzos actuantes que es empleado para producir
el deslizamiento de la falla (Aki y Richards, 2002).
86
manera individual mayores o menores (Madariaga, 1994) y ello se debe a que
las rupturas sísmicas son extremadamente complejas.
Figura 8.1. Relación entre área de ruptura y momento sísmico. Las líneas 0.1, 1, 10 y
100MPa muestran una caída de esfuerzos (Δσ) constante, (Kanamori y Anderson, 1975).
87
Figura 8.2. Relación entre el logaritmo del momento sísmico escalar y el área de la fractura
para terremotos con foco superficial, (Tavera, 1998).
Otro de los estudios que ha analizado con detalle la variación de los valores
de caída de esfuerzos para sismos de variada magnitud (1.5-7ML), es el de
Abercrombie y Leary (1993). Estos autores realizaron la recopilación de 800
eventos sísmicos (Terremotos de mina, sismicidad volcánica e inducida,
réplicas de grandes terremotos e hidrofracturas) con epicentros en California
y el mundo, sus magnitudes son pequeñas, moderadas y grandes, esta
recopilación considera sismos provenientes de varios estudios como se
observa en la Figura 8.3.
88
k) Réplicas de Loma Prieta (Guo et al., 1992)
l) Hidrofracturas (Fehler and Phillips, 1991)
m) Terremotos recientes (D. Wald, 1993).
Figura 8.3. Relación entre el logaritmo del momento sísmico escalar (Mo) y el radio de la
fractura (r), (Abercrombie y Leary, 1993).
Según la Figura 8.3 los valores de caída de esfuerzos para los 800 sismos
analizados varían en un rango de 0.1 a 100Mpa. En la Figura 8.4 se presenta
los valores de caída de esfuerzos obtenidos para las réplicas del terremoto de
Arequipa del 23 de Junio del 2001 analizadas en cinco estaciones sísmicas
por lo que se tiene la Figura 8.4.a) Estación Tocota, Figura 8.4.b) Estación
Ocoña, Figura 8.4.c) Estación Caravelí, Figura 8.4.d) Estación Bella Unión y
Figura 8.4.e) Estación Quilca, las réplicas del Terremoto de Arequipa
analizadas están representadas por los círculos color verde claro y se ubican
sobre los asteriscos color celeste que corresponden a las réplicas del Lago
89
Mammoth con magnitud 3.9≤ ML≤5.1 (Archuleta et al., 1982), los círculos
amarillos representan los Terremotos de California con magnitud 3≤ML≤7
(Thatcher y Hanks, 1973), triángulos verdes que representan las réplicas de
Loma Prieta 1.5≤ML≤4.5 (Guo et al., 1992) y los triángulos rosados que son
las réplicas de Palm Spring, en California con magnitud 3.4≤ML≤4.4 (Mori y
Frankel, 1990). En la Figura 8.4.a) están representadas las réplicas de la
Estación Tocota, están delimitadas por la elipse roja y la mayoría de réplicas
presentan valores de momento sísmico escalar dentro del rango de 1012 a
1015N.m equivalente a una magnitud de 2.4-4.4Mw y el radio de la fractura
varía de 190 a 500m, además se tienen dos réplicas que se encuentran
aisladas, son las réplicas número 32 y 68 para ambas réplicas el radio de la
fractura es del orden de 4000m y el momento sísmico escalar es del orden de
1017N.m equivalente a una magnitud de 5.4Mw. Es importante señalar que el
diagrama de caída de esfuerzos elaborado para la estación Tocota, se
considera el diagrama de caída de esfuerzos representativo de la
investigación debido a que los registros sísmicos de la estación Tocota tienen
mejor calidad que las estaciones Ocoña, Caravelí, Bella Unión y Quilca.
Figura 8.4.a) Relación entre el logaritmo del momento sísmico escalar (Mo) y el radio de la
fractura (r), Estación Tocota. Modificado Abercrombie y Leary, (1993).
90
La Figura 8.4.b) Estación Ocoña gráficamente muestra que los valores de
momento sísmico escalar están en el rango de 1013 a 1016N.m y el radio de la
fractura para la mayor cantidad de réplicas es de 190 a 600m
aproximadamente, existen dos réplicas que aparecen alejadas del grupo
principal, corresponden a la réplica número 32 y la 68 para estas réplicas el
momento sísmico escalar está en el orden de 1017N.m y el radio de fractura
es de 2000m.
Figura 8.4.b) Relación entre el logaritmo del momento sísmico escalar (Mo) y el radio de la
fractura (r), Estación Ocoña. Modificado Abercrombie y Leary, (1993).
91
Figura 8.4.c) Relación entre el logaritmo del momento sísmico escalar (Mo) y el radio de la
fractura (r), Estación Caravelí. Modificado Abercrombie y Leary, (1993).
Figura 8.4.d) Relación entre el logaritmo del momento sísmico escalar (Mo) y el radio de la
fractura (r), Estación Bella Unión. Modificado Abercrombie y Leary, (1993).
92
La estación Bella Unión registró 116 réplicas del terremoto de Arequipa, sin
registrar la réplica número 32 y la 68, en la Figura 8.4.d) se observa que los
valores de momento sísmico escalar están en el orden de 1013 a 1016N.m y el
radio de la fractura en el orden de 200 a 500m. En la Figura 8.4.e) se
presentan las réplicas registradas en la estación Quilca que registró 35
réplicas, se observa la réplica 32 que está fuera de la elipse roja, los valores
de momento sísmico escalar de las réplicas están en el orden de 1013 a
1017N.m y el radio de fractura está en el orden de 240 a 4900m.
Figura 8.4.e) Relación entre el logaritmo del momento sísmico escalar (Mo) y el radio de la
fractura (r), Estación Quilca. Modificado Abercrombie y Leary, (1993).
93
momento sísmico escalar se encuentran en el rango de 1012 a 1016, además se
presentan en forma alejada del grupo principal dos réplicas (réplica 32 y 68)
con un momento sísmico escalar de 1017N.m.
Figura 8.5. Relación entre el logaritmo de la frecuencia esquina, el momento sísmico (escala
inferior) y la Magnitud Momento Sísmico (escala superior). Las líneas rojas representan la
caída de esfuerzos de 0.1, 1, 10 y 100Mpa (Allman y Shearer, 2009).
94
En la Figura 8.6 se han representado los resultados de frecuencia esquina (fc),
momento sísmico escalar (M0)y Magnitud Momento (Mw) obtenidos en esta
investigación, con la finalidad de observar su comportamiento, para las
estaciones 8.6.a) Estación Tocota, 8.6.b) Estación Ocoña, 8.6.c) Estación
Caravelí, 8.6.d) Estación Bella Unión y 8.6.e) Estación Quilca, (círculos color
naranja dentro de la elipse roja), observándose correlación con los resultados
obtenidos por Allman y Shearer (2009); por lo tanto, la caída de esfuerzos
(∆𝜎) no depende de la magnitud del sismo, el valor siempre es constante. De
acuerdo a la figura de Allman y Shearer las réplicas del terremoto de
Arequipa del 2001 se ubican junto a las réplicas del Nor-este de Norte
América con magnitud 3≤ML6≤ que son las estrellas rojas (Boatwright,
1994), cerca de los rombos en azul que corresponden a 82 terremotos de
México con magnitudes de 3-7ML (Humphrey y Anderson, 1994), las
estrellas verdes son las réplicas de Loma Prieta (Hough,1996), cuadrados en
amarillo réplicas de Palm Spring (Mori y Frankel,1990), triángulos en celeste
réplicas de Lago Mammoth (Archuleta et al., 1982). En la Figura 8.6.a) las
réplicas de la estación Tocota se encuentran dentro de la elipse roja, en el eje
de ordenadas se tiene la frecuencia esquina, sus unidades son Herzt (Hz), en
la escala inferior el momento sísmico escalar (M0) sus unidades
Newton.metro y en la escala superior la Magnitud Momento (Mw). Los
valores de frecuencia esquina varían de 3.5 a 10Hz, sin embargo las réplicas
32 y 68 tienen frecuencias esquina de 0.3 y 0.4Hz, es probable que ambos
valores de frecuencia esquina se encuentran en este rango debido a la calidad
del registro sísmico, el momento sísmico escalar está en el rango de 10 12 a
1015N.m y la magnitud momento para la mayoría de réplicas es de 2.5 a 5Mw
de acuerdo a la Figura 8.6.a), esto es coherente con los resultados obtenidos
del análisis espectral, la magnitud momento determinada para la Estación
Tocota es de 2.36 a 5.48Mw, entonces el diagrama de tres escalas utilizado
permite correlacionar estos tres parámetros de la fuente sísmica, los valores
de caída de esfuerzos varían de 0.1 a 100Mpa.
La Figura 8.6.b) presenta las réplicas de la estación Ocoña los valores de
frecuencia esquina del grupo grande varían de 3 a 8Hz, en tanto que las
95
réplicas 32 y 68 tienen como frecuencia esquina 0.6 y 0.7Hz. Los valores de
momento sísmico escalar está en el rango de 1013 a 1016N.m. La magnitud
momento de las réplicas está comprendida de 2.7 a 5Mw.
96
Las réplicas de la estación Caravelí se representan en la Figura 8.6.c) la
mayor cantidad de réplicas tiene una frecuencia esquina de 3 a 7Hz, están
presentes la réplica 32 y 64 con frecuencias esquina de 0.4Hz, los valores de
momento sísmico están en el rango de 1013 a 1016N.m. La magnitud momento
varía de 2.7 a 5,4Mw.
En la Figura 8.6.d) se tienen las réplicas de la estación Bella Unión (sin las
réplicas número 32 y 68), los valores de frecuencia esquina están en el
intervalo de 3 a 7Hz, el momento sísmico escalar está en el orden de 1013 a
1015N.m y la magnitud momento está comprendida de 2.6 a 4.5Mw.
Las réplicas de la estación Quilca se representan en la Figura 8.6.e) la mayor
cantidad de réplicas tiene una frecuencia esquina de 3 a 7 Hz, está presente la
réplica 32 con frecuencia esquina de 0.4 Hz, los valores de momento sísmico
están en el rango de 1012 a 1015N.m. La magnitud momento varía de 2.5 a
4Mw.
97
Figura 8.6.d) Relación entre el logaritmo de la frecuencia esquina, el momento sísmico
(escala inferior) y la Magnitud Momento Sísmico (escala superior) Estación Bella Unión. Las
líneas rojas representan la caída de esfuerzos de 0.1, 1, 10 y 100Mpa (Modificado Allman y
Shearer, 2009).
98
Finalmente los valores de caída de esfuerzos para las réplicas del sismo de
Arequipa son similares a los obtenidos por Abercrombie y Leary (1993) para
terremotos de California y el mundo. Los parámetros de la fuente como
momento sísmico, radio de la fractura, caída de esfuerzos, magnitud momento
y magnitud energía han sido correctamente determinados, con valores
razonables para sus rangos de magnitud, como se puede demostrar en la Tabla
7.3. Estación Tocota y en el Anexo 1, donde se presentan los resultados para
las estaciones Ocoña, Caravelí, Bella Unión y Quilca.
99
CONCLUSIONES
Figura 7.8. Relación entre el logaritmo de Mo y r para las réplicas del terremoto de Arequipa,
registradas en la estación Tocota.
100
RECOMENDACIONES
101
BIBLIOGRAFÍA
Aki, K. y Richards, G., (1980). Quantitative Seismology, Freeman and Co., New
York.
Aki, K., (1966). 4. Generation and propagation of G waves from the Niigata
earthquake of June 14, 1964. Part 2. Estimation of earthquake moment, released
energy and stress-strain drop from G wave spectrum. Bulletin of the Earthquake
Research Institute 44: 73–88.
Aki, K., (1967). Scaling law of seismic spectrum, J. Geophys. Res., 72, 1217–
1231.
Allmann, B. y Shearer, P., (2009). Global variations of stress drop for moderate to
large earthquakes. J. Geophys. Res. 114, B01310, doi: 10.1029/2008JB005821.
Boatwright, J., (1980). A spectral theory for circular seismic sources: simple
estimates of source dimension, dynamic stress drop, and radiated energy. Bull.
Seism. Soc. Am, 70, 1-26.
Boatwright, J., (1980). Preliminary body-wave analysis of the St. Elias, Alaska
earthquake of february 28, 1979. BSSA, 70, 419-436.
Brigham, E. O., (1988). The Fast Fourier Transform and Its Applications,
Prentice-Hall, Englewood Cliffs, NJ.
Brune, J. N., (1970). Tectonic stress and the spectra of seismic shear waves from
earthquake, J. Geophys. Res., 75, 4997–5009.
Cahill, T. y Isacks B., (1992). Seismicity and shape of the subducted Nazca plate.
J. Geophys. Res, 97, 17.503-17.529.
Choy, G. L., y Boatwright, J., (1995). Global patterns of radiated seismic energy
and apparent stress. J. Geophys. Res., 100, B9, 18,205-18,228.
Cutipa, G., (2006). Análisis de los niveles de ruido sísmico en estaciones sísmicas
de Banda Ancha de la Red Sísmica Nacional del Perú. Tesis de grado UNSA. 138
p.
Dass, S y Aki, K., (1977). Fault planes with Barriers: A versatile earthquake
model, J. Geophys. Res., p.5658-5670.
Goodman, R.E. (1989). Introduction to Rock Mechanics. New York, John Wiley,
562 p.
Grange, F., (1984). Etude sismotectonique detaille de la subduction lithospherique
au sud Perou. Tesis Doctoral, Origm, Grenoble, Francia.
Gutenberg, B., y Richter, C., (1956). Magnitude and energy of earthquakes, Ann.
Geofis., 9, 1-15.
Hanks, T. C. y Wyss, M., (1972). The use of body wave spectra in the
determination of seismic source parameters. Bull seism. Soc. Am., 62, 561-589.
Haskell, N.A., (1964). Total energy and energy spectral density of elastic wave
radiation from propagating faults, Bull. Seism. Soc. Am., 54, 1, 811 1,841.
Heaton, T. H., (1990) Evidence for and implications of self-healing pulses of slip
in earthquake rupture, Phys. Earth Planet, Int.64, 1-20.
Hough S.E., Seeber, L., Lerner-Lam, A., Armbruster, J.C. & Guo, H., (1991).
Empirical Green’s function analysis of Loma Prieta aftershocks, Bull. Seism.
Soc. Am.,81. 1737–1753.
Klein, F.W., (1978). Hypocenter location program HYPOIN VERSE Part 1. Users
Guide to Versions 1, 2, 3 and4. US. Geol. Surv. Open-File Rept, 78-694.
Kulhanek, O., Torild Van Eck, Norris . J., Meyer. K. y Wahlstroem, R., (1983).
Spectra of the earthquake sequence February-March. 1981, in South-Cental
Sweden. Tectonophysics, 93: 337-350.
Lindo, R., (1993). Seismotectonique des andes du Perou Central: Apport des
donnes sismologiques de haute precisión. Tesis Doctoral. Universidad Louis-
Pasteur de Strasbourg, Francia.
Macharé, J. Fenton, C.H., Machette, M.N, Lavenu, A., Costa, C. y Dart,. ( 2003).
Map and Database of Quaternary Faults and Folds in Peru ans its offshore
regions: U.S. Geological Survey Open-File Report 03- 0451, 54 p., 1 plate
(1:3,000,000 scale).
Madariaga, R., (1976). Dynamics of an expanding circular fault, Bull. Seism. Soc.
Am., 66, 639–666.
Madariaga, R., (1979). On the relation between seismic moment and stress drop in
the presence of stress and strength heterogeneity. Journal of Geophysical
Research. Vol.84, N°.B5.8B0926.
Madariaga, R., (1994). Dynamic of seismic source. Física de la tierra, Núm.6.29-
72. Editorial Complutense.Madrid.
Molnar, P., B.E. Tucker, and J.N. Brune (1973): Corner frequencies of P and S
waves and models of earthquake sources; Bull. Seism. Soc. Am., 63, 2091-2104.
Mori, J. y Frankel A. (1990). Source parameters for aftershocks of the 1986 North
Palm Springs earthquake, Bull. Seism, Soc. Am. 80, p.278-295.
Muller, T. y Zurn, W., (1983). Observation of gravity changes during the passage
of cold fronts, J. Geophys., 53, 155–162.
Sebrier, M., Huaman, D., Blanc, J., Machare, J., Bonnot, D. y Cabrera, J. (1985).
Observaciones acerca de la neotectónica del Perú. IGP-LDGDI pro. Sisra, Grup.
Peruano de Neotectónica, 96 pp.
Scherbaum, F. y Stoll, D., (1983). Source parameters and scaling laws of the 1978
Swabian Jura (southwest Germany) aftershocks, Bull. Seismol. Soc. Am., 73,
1321-1343.
Silgado, E., (1978). Historia de los sismos más notables ocurridos en el Perú,
1513-1974. Instituto Geológico Minero del Perú, Lima,131 p.
Stauder, W., (1975). Subduction of the Nazca plate under Perú as evidenced by
focal mechanism and by seismicity, J. Geophys. Rs. 80, 053-1064.
Tajima, R., y F. Tajima (2007). Seismic scaling relations and aftershock activity
from the sequences of the 2004 mid Niigata and the 2005 west off Fukuoka
earthquakes (Mw 6.6) in Japan, J. Geophys. Res., 112, B10302.
Tavera, H. (1992). El Proceso de Ruptura Sísmica: Barrera o Aspereza. Boletín de
la Sociedad Geológica del Perú. Vol. 83.p. 69-73
Tavera, H., (2002). Terremoto de la región Sur de Perú del 23 de Junio del 2001.
Centro Nacional de Datos Geofísicos.
Tavera, H., Buforn, E., Bernal, I., Antayhua, Y. y Vilcapoma, L., (2002). The
Arequipa (Peru) earthquake of June 23, 2001. Journal of Seismology 6: 279-283.
Tavera, H., Fernández, E., Bernal, I., Antayhua, Y., Agüero, C., Salas, H.,
Rodriguez, S., Vilcapoma, L., Zamudio, Y., Portugal, D., Inza, A., Carpio, J.,
Callo, F. y Valdivia, I., (2006). The southern región of Peru earthquake of June
23rd, 2001. Journal of Seismology (2006). Doi: 10.1007/s10950-006-9014-2.
Thatcher, W., (1972). Regional variations and seismic source parameters in the
northern Baja California area. J. Geophys. Res., 77: 1549-1566.
Parámetros de la fuente sísmica de las réplicas del Terremoto de Arequipa del 2001
*La Tabla 7.3 con resultados de la Estación Tocota se encuentra en el Capítulo VII.
ANEXO 1
(*)Se debe tener en cuenta que en la parte superior derecha de cada uno de los espectros
aparece el nombre de la estación sísmica, la fecha y la hora de la réplica, sin embargo puede
no coincidir con la hora UTC (Universal Time Coordinated) de la réplica, por ello se debe
considerar el número de réplica de la Tabla 7.3.
a. Estación Tocota
a. Continuación
a. Continuación
a. Continuación
a. Continuación
b. Estación Ocoña
b. Continuación
b. Continuación
b. Continuación
b. Continuación
c. Estación Caravelí
c. Continuación
c. Continuación
c. Continuación
c. Continuación
d. Estación Bella Unión
d. Continuación
d. Continuación
d. Continuación
d. Continuación
e. Estación Quilca
e. Continuación
e. Continuación
e. Continuación
e. Continuación