Tavera Condori 2011

UNIVERSIDAD NACIONAL DE SAN AGUSTN DE AREQUIPA
FACULTAD DE GEOLOGA, GEOFSICA Y MINAS

ESCUELA PROFESIONAL DE INGENIERA GEOFSICA
REAS PROBABLES DE RUPTURA SSMICA E EL BORDE

OCCIDETAL DEL PER, A PARTIR DE LA VARIACI DEL
PARMETRO b
Tesis Presentada por el Bachiller en
Ciencias Geofsicas:
Cristobal Condori Quispe
Para optar el Titulo Profesional de
INGENIERO GEOFSICO
AREQUIPA - PER
2011
Dedicada a toda
Mi Familia
Quienes me apoyaron en mi formacin profesional
Los amo mucho.
AGRADECIMIETOS
Deseo expresar mi gratitud a todas aquellas personas que con sus sugerencias
contribuyeron a la realizacin de esta tesis.
En primer lugar, debo agradecer a Dios por bendecirme y guiarme en el buen
camino, por darme la vida y los amigos que tengo. Deseo tambin agradecer a mi
padre que seguramente esta en alguna parte muy cerca de Dios. Te extrao mucho.
Doy mi agradecimiento eterno a mi Asesor de Tesis, Dr. Hernando Tavera
Director del rea de Sismologa del Instituto Geofsico del Per (IGP), por su valiosa
asesora y por brindarme la oportunidad de formarme profesionalmente en el campo
de la investigacin en Sismologa.
Mi agradecimiento al Instituto Geofsico del Per, por haberme brindado una
beca y facilitado un ambiente agradable de trabajo y la logstica necesaria para
concluir esta tesis.
Quiero agradecer especialmente a los Magsteres Juan Carlos Villegas, Sheila
Yauri, Isabel Bernal y Bhila Herrera, por todos los conocimientos que compartieron
conmigo y por sus valiosos tiempos dedicados a ayudarme a lograr esta tesis.
A todo el personal que labora en el rea de Sismologa (IGP). De igual manera,
a todos mis compaeros tesistas.
A mis profesores de la Universidad Nacional de San Agustn, M.Sc. Armando
Minaya, Jorge Soto, Javier Ticona, Edgar Gonzales, Dr. Orlando Macedo e Ing.
Sebastin Ziga, muchas gracias por sus enseanzas. Asimismo, a los Ingenieros
Vctor Aguilar, Robert Salazar y Jos Campano del Instituto Geofsico de Characato,
lugar donde inicie mis primeras prcticas pre-profesionales en el campo de la
Sismologa.
A Gloria, por darme su amor, apoyo, confianza, y compartir nuevos e
inolvidables momentos en mi vida, te amo y espero contar siempre contigo.
DICE
RESUME
CAPTULO I: INTRODUCCIN
1.1 REA DE ESTUDIO
1.2 OBJETIVOS
1.3 TRABAJOS ANTERIORES
CAPTULO II: CONTEXTO GEODINMICO DEL PER

2.1 PRINCIPALES RASGOS TECTNICOS
2.1.1 Dorsal de Nazca
2.1.2 La Fractura de Mendaa
11
2.1.3 La Fosa Marina
11
2.1.4 La Cordillera Andina
11
2.1.5 La Cadena Volcnica
12
2.1.6 Los Sistemas de Fallas
12
2.2 ZONIFICACIN MORFOLGICA
13
2.2.1 La Cordillera Costanera
13
2.2.2 La Zona Costera
13
2.2.3 La Cordillera Occidental
14
2.2.4 Zona Intercordillerana
14
2.2.5 La Cordillera oriental
16
2.2.6 Zona Subandina
16
2.2.7 Llanura Amaznica
16
CAPTULO III: CARACTERSTICAS GENERALES DE LA SISMICIDAD EN

EL PER
3.1 PERIODO HISTRICO
20
3.2 PERIODO INSTRUMENTAL
23
3.3 DISTRIBUCIN ESPACIAL
27
3.3.1 Sismicidad con foco Superficial
27
3.3.2 Sismicidad con foco Intermedio
28
3.3.3 Sismicidad con foco Profundo
28
3.4 PERFILES SSMICOS
28
3.4.1 Regin Norte
29
3.4.2 Regin Centro
29
3.4.3 Regin Sur
29
3.5 GEOMETRA DE LA SUBDUCCIN
31
3.6 DISTRIBUCIN DE REAS DE RUPTURA
35
CAPTULO IV: MODELOS DE PROCESOS DE RUPTURA

4.1 MODELO DE FUENTE SISMICA
41
4.1.1 MODELO DE BARRERAS
41
4.1.2 MODELO DE ASPEREZAS
42
CAPTULO V: METODOLOGA PARA EL CLCULO DEL PARAMETRO b

5.1 DISTRIBUCIN FRECUENCIA-MAGNITUD
47
5.1.1 Los Parmetros a y b
50
5.1.1.1 El Parmetro a
50
5.1.1.2 El Parmetro b
50
5.2 ASPEREZAS Y VALORES MNIMOS DEL PARAMETRO b
53
5.3 VARIACIONES DEL PARMETRO b
54
5.4 MTODOS PARA EL CALCULO DE a y b
56
5.4.1 Mtodo de Mnimos Cuadrados
56
5.4.2 Mtodo de Mximo Verosimilitud
57
5.5 METODOLOGIA PARA EL CALCULO DE a y b
59
5.6 PERIODOS DE RECURRENCIA
62
5.7 PROBABILIDAD DE RECURRENCIA
64
5.8 APLICACIONES DEL PARMETRO b
65
5.8.1 Estudios en las Zonas de Subduccin
66
5.8.1.1 Zona de Subduccin de Mxico
66
5.8.1.2 Zona de Subduccin Alaska y Nueva

Zelanda
67
5.8.1.3 Zona de Subduccin de Per
68
5.8.2 Sistemas de Fallas
70
5.8.3 Aplicacin en Volcanes
71
CAPTULO VI: ANLISIS Y SELECCIN DE DATOS (APLICACIN AL

BORDE OCCIDENTAL DEL PER)
6.1 BASE DE DATOS SSMICO
6.1.1 Catlogo Ssmico del Per
75
75
6.2 HOMOGENEIZACIN DE LA MAGNITUD
79
6.3 MAGNITUD DE COMPLETITUD (Mc)
82
6.4 ELIMINACIN DE RPLICAS
82
6.4.1 Ley de Omori-Utsu
82
6.4.2 Algoritmo de Reasenberg
83
6.5 BASE FINAL DE DATOS SSMICOS
87
CAPTULO VII: DISTRIBUCIN DEL PARMETRO b (APLICACIN AL

BORDE OCCIDENTAL DEL PER)
7.1 DISTRIBUCION DEL PARMETRO b
91
7.2 RESULTADOS DEL PERIODO DE RECURRENCIA
94
7.3 RESULTADOS DE LA PROBABILIDAD DE RECURRENCIA
96
7.4 ESTIMACIN DE LA MAGNITUD A PARTIR DEL AREA DE

ASPEREZAS
CAPTULO VIII: DISCUSIN E INTERPRETACIN
98
101
COCLUSIOES
RECOMEDACIOES
BIBLIOGRAFA
AEXOS
Anexo A.
Resumen extendido de los resultados preliminares de esta

investigacin presentado al Congreso peruano de Geologa Cusco2010.
RESUME
Se ha evaluado la presencia de asperezas para el borde occidental del Per,
propensas a producir, en el futuro, sismos de gran magnitud, todos asociados al
proceso de subduccin de la placa de Nazca bajo la Sudamericana. Para tal fin, se
analiza la variacin espacio-tiempo del parmetro b siguiendo la metodologa
propuesta por Wiemer y Wyss. (1997). La data utilizada considera 1367 eventos
ssmicos del catlogo ssmico del Instituto Geofsico del Per (periodo 1970 a 2010),
todos con focos ubicados a profundidades menores a 60 Km. Con una magnitud de
completitud de 3.8Ms. Los resultados obtenidos permiten identificar la existencia de
al menos 5 asperezas con valores anmalos del parmetro b en las cuales existe
mayor acumulacin de energa elstica. Histricamente, estas asperezas produjeron
en el pasado terremotos de gran magnitud. La primera aspereza (A1) se encuentra
frente a la costa del extremo sur del departamento de Arequipa, Moquegua y Tacna
(terremoto de 1868); la segunda (A2), frente a la costa del extremo
norte del
departamento de Arequipa (terremoto de 1913); la tercera (A3) y cuarta (A4), frente

a la costa de los departamentos de Lima y Ancash (terremoto de 1746) y la quinta
(A5), frente a la costa de los departamentos de Lambayeque y Piura (terremoto de
1619).
Los resultados del estudio muestran que las asperezas A1, A2 y A3 presentan
probabilidades de un 75% de producir en los prximos 50 aos, terremotos con
magnitud mayores a 7.0Ms.
Captulo I: Introduccin
CAPTULO I
ITRODUCCI
El Per forma parte de la regin ssmica ms activa del mundo debido a que se
encuentra ubicado dentro del denominado Cinturn de Fuego del Pacfico, donde se
libera ms del 80% de la energa ssmica a nivel mundial. En Per la actividad
ssmica est directamente asociada al proceso de subduccin de la placa de Nazca
bajo la Sudamericana, el cual se da a una velocidad relativa de 6 a 7 cm/ao (DeMets
et al., 1990; Norabuena et al., 1999; Kendrick, et al., 2003, Villegas, 2009).
Este proceso dio origen a los terremotos mas grandes ocurridos en el Per,
como los del 17 de octubre de 1966 (7.7Mw Lima), 31 de mayo de 1970 (7.9Mw
Ancash), 3 de octubre de 1974 (7.9Mw Lima), 12 noviembre de 1996 (7.7Mw
Nazca), 23 de junio de 2001 (8.2Mw Arequipa), y 15 de agosto de 2007 (7.9Mw
Pisco). Todos estos terremotos produjeron cuantiosos daos materiales y prdidas de
vidas humanas.
En la actualidad se tiene conocimiento amplio acerca de este peligro, ms no se
sabe cuando ha de ocurrir, es por eso que durante los ltimos aos varios
investigadores han desarrollado diversos mtodos de prediccin, sin lograr el
objetivo de conocer la fecha de recurrencia de los grandes sismos (Tavera, 2003). Sin
embargo, s ha sido posible caracterizar y localizar las reas probables donde podra
generarse un futuro terremoto de magnitud elevada (7.0Ms).
1
reas probables de ruptura ssmica en el borde occidental del Per, a partir de la
variacin del parmetro b.
Una de las metodologas propuestas para identificar estas reas fue

desarrollada por Wiemer y Wyss (1997), y se basa en la distribucin frecuenciamagnitud de los sismos Log!=a-bM (Gutenberg y Richter, 1944). Estos
investigadores, sealan que las zonas con alta concentracin de esfuerzos
(Asperezas) se caracterizan por presentar valores bajos del parmetro de b, siendo
sus dimensiones proporcionales a las del rea de ruptura y a la magnitud del sismo.
As mismo; este mtodo permite calcular el periodo de retorno de los sismos
asociados a dichas asperezas. El trmino Aspereza define a un rea presente en la
superficie de friccin de placas en la cual no ha ocurrido un sismo de gran magnitud
durante un largo periodo de tiempo, lo que ocasiona que se acumule la energa
elstica suficiente como para producir un gran sismo.
El mtodo propuesto por Wiemer y Wyss (1997), fue aplicado en diversas

zonas de subduccin en el mundo, tal es el caso de Ziga y Wyss (2001) que
buscaron localizar las grandes asperezas a lo largo de la Costa occidental de Mxico
y Katsumata, (2006), en el estudio de la variacin del valor de b en profundidad en la
misma regin. En ambos casos, el objetivo fue identificar la presencia de asperezas o
reas con mayor probabilidad de ocurrencia de sismos de magnitud 7.2Ms.
En el Per, el mtodo fue aplicado por Heras y Tavera (2002) al borde

occidental, llegndose a identificar asperezas frente a los departamentos de Arequipa
e Ica. La primera dio origen el sismo el 23 de junio del 2001 y el segundo, el sismo
de Pisco del 15 de agosto del 2007.
En el presente estudio, se emplea el mtodo de Wiemer y Wyss (1997),

para analizar el valor de la distribucin espacial del parmetro b en el borde
occidental del Per y a partir de los resultados identificar la ubicacin de nuevas
Asperezas.
El desarrollo de este trabajo de investigacin se ha estructurado en nueve

captulos, que permiten entender de manera general el contenido del estudio. A
continuacin se describe cada uno de estos captulos.
2
En el Captulo I, se realiza la descripcin y se detalla la importancia de aplicar

el mtodo de Wiemer y Wyss (1997) al borde Occidental del Per. Se describe los
objetivos del estudio, la definicin del rea de estudio y finalmente se hace una
referencia a los trabajos anteriores aplicados en el Per y en otros lugares.
En el Captulo II, se describe el contexto geodinmico del Per considerando

sus principales rasgos tectnicos y las unidades estructurales que son responsables de
la actual geodinmica superficial.
En el Captulo III, se discute sobre las caractersticas generales de la sismicidad

histrica e instrumental del Per, as como su distribucin espacial y el anlisis de las
reas de ruptura de los grandes terremotos ocurridos en su borde occidental.
En el Captulo IV, se explica y discute los modelos propuestos para describir las
caractersticas de las reas de ruptura asociados a grandes terremotos con origen en la
zona de subduccin y/o cortical.
En el Captulo V, se realiza el anlisis y descripcin del mtodo desarrollado

por Wiemer y Wyss (1997), para el clculo del parmetro b. Asimismo, se discute
los fundamentos sobre el cual se basa la relacin frecuencia-magnitud definido por
Gutenberg y Richter en 1954. Seguidamente, se describen algunos estudios
realizados en diversas zonas de subduccin, sistemas de fallas y regiones volcnicas
aplicando esta metodologa.
En el Captulo VI, se realiza el anlisis y seleccin de la base de datos ssmicos

a utilizarse en el presente estudio para el clculo del parmetro b.
En el Captulo VII,
se presenta los resultados obtenidos para el borde
occidental del Per a partir de la aplicacin del mtodo de Wiemer y Wyss (1997).
Se evala los valores del parmetro b, periodos de retorno y probabilidad de
ocurrencia para terremotos con magnitudes 7.0Ms. As mismo, se realiza la
estimacin de las magnitudes Mw para las asperezas identificadas.
3
En el Captulo VIII, se realiza la discusin de los resultados obtenidos en este

estudio en relacin a los obtenidos por otros autores.
Finalmente, se presenta las conclusiones y recomendaciones propuestas a partir

de los resultados obtenidos en este estudio. Asimismo, se presenta un anexo
conteniendo una copia del artculo de investigacin realizado por el autor durante la
elaboracin del presente estudio, el mismo ha sido publicado en congresos y revistas
especializadas.
1.1 REA DE ESTUDIO

La sismicidad y la tectnica del Per son controladas por el proceso de
subduccin entre dos placas, la de Nazca (ocenica) bajo la Sudamericana
(continental), y cuyo contacto se inicia en la fosa Per-Chile. Este proceso genera un
gran nmero de sismos de diversas magnitudes, siendo los mayores, los que
producirn daos en las ciudades prximas al epicentro. Dentro de este contexto es
necesario conocer la ubicacin de las reas posibles que podran generar un sismo en
el futuro.
En este estudio se aplica la metodologa propuesta por Wiemer y Wyss (1997)

al borde occidental de Per (Figura 1.1), delimitado por las siguientes coordenadas
geogrficas, de -2,5 a -20,0 de Latitud Sur; abarcando todos los departamentos que
se ubican a lo largo de la lnea de costa: Tacna, Moquegua, Arequipa, Ica, Lima,
Ancash, la Libertad, Lambayeque, Piura y Tumbes.
4
Figura 1.1: Mapa topogrfico del Per mostrando la ubicacin del rea de estudio
delimitado con el polgono de color azul, y que comprende todo el borde Occidental del
Per. Las flechas rojas indican la direccin del movimiento de placas.
5
1.2 OBJETIVOS
Los objetivos del presente estudio son:
Analizar la distribucin espacial del parmetro de b en todo el borde
occidental del Per.
Correlacionar las reas de valores bajos del parmetro de b con reas de
ruptura asociados a grandes terremotos ocurridos en todo el borde
occidental del Per.
Calcular los periodos de recurrencia y densidad de probabilidad para
sismos con magnitud mayor/igual a 7.0Ms en todo el borde occidental del
Per.
Analizar las Asperezas y estimar su magnitud Mw en base a las relaciones
propuestas por Kanamori y Anderson (1975).
1.3 TRABAJOS ATERIORES
La primera aplicacin del mtodo de Wiemer y Wyss (1997), para el borde

occidental del Per fue realizado por Heras y Tavera (2002), permitiendo obtener
anomalas del parmetro b bien definidas, las mismas se detallan en el Capitulo V.
En este estudio se vuelve aplicar esta metodologa con el objetivo de localizar las
asperezas despus de los terremotos del 2001 y 2007, para ello se cuenta con una
base de datos ms consistentes e homogneos del catalogo ssmico del IGP para el
periodo de 1970-2010.
6
Captulo II: Contexto Geodinmico del Per
CAPTULO II
COTEXTO GEODIMICO DEL PER

El territorio peruano est sometido a una fuerte actividad dinmica debido a la
interaccin de las Placas Sudamericana y Nazca, ambas convergen de Este a Oeste
con una velocidad promedio de 6 a 7 cm/ao (DeMets et al., 1990; Norabuena et al.,
1999; Kendrick, et al., 2003, Villegas, 2009). Durante el desarrollo de este proceso,
en los ltimos 65 millones de aos (Ma), el desplazamiento de las placas produjo en
el continente el levantamiento y resquebrajamiento de rocas sedimentarias
permitiendo la formacin de la cordillera Andina (Soler y Sbrier, 1990; Mergard,
1978). Esta cordillera se extiende paralela al margen costera, desde Venezuela hasta
el sur de Chile cubriendo una longitud de 7000 Km, y con un ancho que oscila entre
250 Km en la regin central del Per y de 500 Km en el lmite entre Per, Chile y
Bolivia. Las fuerzas tectnicas que producen el desplazamiento de las placas generan
una constante acumulacin de energa que se libera en forma de terremotos (Tavera y
Buforn, 1998). En el presente captulo se describe las caractersticas geodinmicas
que controlan la ocurrencia de la sismicidad en el Per.
7
2.1 PRICIPALES RASGOS TECTICOS

La formacin de los rasgos tectnicos caractersticos del Per son debidos al
proceso de la subduccin entre las placas de Nazca y Sudamericana. En general, la
evolucin geodinmica de Per est controlada por los siguientes rasgos tectnicos:
la Dorsal de Nazca, la Fractura de Mendaa, la Fosa Peruano-Chilena, la Cordillera
Andina, la Cadena Volcnica y los diferentes Sistemas de Fallas distribuidas en el
interior del continente. La ubicacin geogrfica de estos rasgos tectnicos se muestra
en la Figura 2.1.
2.1.1 Dorsal de azca: es una cordillera ocenica que se localiza en el
extremo NO de la regin Sur de Per frente al departamento de Ica. Esta
cordillera sigue una orientacin NE-SO perpendicular a la lnea de la fosa
peruanochilena (entre 15 y 24 Sur), de tal modo que su extremo NE se
ubica frente al departamento de Ica en donde presenta un ancho de
aproximadamente 220 km sobre la cota de 2000 metros. Sin embargo; su
ancho y altitud disminuye gradualmente hacia su extremo SO. Segn la Figura
2.1, la cota de 2000 metros de esta dorsal, se localiza a 50 km de distancia
aproximadamente de la lnea de fosa; mientras que, las cotas menores ya
subducieron bajo la placa continental (Sebrier et al, 1985). La Dorsal de Nazca
presenta una forma asimtrica, siendo probablemente esta caracterstica
determinante para los diferentes procesos geodinmicos que se producen en
esta regin.
Estudios recientes, sobre anomalas magnticas, permiten considerar la

hiptesis de que la Dorsal de Nazca debe su origen a una antigua zona de
creacin de corteza que ces su actividad hace 5 a 10 millones de aos
aproximadamente (Marocco, 1980; Sebrier et al, 1985).
Las principales consecuencias de la migracin lateral de la dorsal de Nazca

son:
8
1. Cambio en la geometra de la Costa. Se observa que la distancia entre la

fosa y la costa es mayor entre 6 y 14 Sur como consecuencia del paso
de la dorsal de Nazca.
2. Dio origen a dos tipos de subduccin: subhorizontal en la regin norte y

centro, y normal en la regin sur del Per.
3. Modific el patrn de la sismicidad en el Per. En la regin norte los
sismos son menos frecuentes, en la regin centro predominan sismos
hasta magnitudes de 8Ms y en la regin sur con magnitudes mayores a
8Ms
4. Ocasiono el cese de la actividad volcnica en el norte y centro del Per.
9
Figura 2.1: Principales rasgos tectnicos presentes en el Per, La Fosa Per-Chile,

La Dorsal de 'azca, La Fractura de Mendaa, La Cordillera de los Andes, La
Cadena Volcnica (tringulos) y los principales sistemas de Fallas (AM=Amazonas,
CB=Codillera Blanca, SA=Satipo, HUA=Huaytapallana, TM=Tambomachay,
MD=Madre de Dios) (Machar et al., 2003).
10
2.1.2 La Fractura de Mendaa: la fractura de Mendaa corresponde a una

discontinuidad de la corteza ocenica que se localiza en el extremo NO de la
regin central de Per, entre 10 y 12 de Latitud Sur frente al departamento de
Ancash. Esta fractura es perpendicular a la lnea de la fosa peruano-chilena,
con un ancho de 150 Km y una orientacin N62E. En la Figura 2.1 se observa
que la fractura de Mendaa se encuentra aproximadamente a una distancia de
180 Km de la lnea de fosa.
La factura de Mendaa sera una antigua zona de divergencia de placa. En
la actualidad, no se dispone de mayor informacin sobre las caractersticas
fsicas de esta fractura.
2.1.3 La Fosa Marina: marca el inicio de la subduccin de la placa de Nazca

bajo la de Sudamrica sobre una longitud de ms de 5000 Km desde Colombia
hasta la Tierra del Fuego en Chile. Se le denomina fosa Peruana-Chilena,
porque abarca toda la costa de Per y Chile. Litolgicamente, la fosa est
formada por sedimentos de diferente potencia depositados sobre rocas preexistentes, la fosa peruana-chilena alcanza una profundidad mxima de 6000
metros y su orientacin es NNO-SSE en la regin Norte y Centro del Per, y
NO-SE en la regin Sur de Per (Heras y Tavera, 2002). El cambio en la
orientacin de la fosa se produce frente a la Dorsal de Nazca (Figura 2.1). La
distancia entre la fosa y la lnea de costa en la regin Norte es de 80 Km, en la
regin central vara entre 60 y 150 Km y en la regin Sur es de 120 Km
aproximadamente.
2.1.4 La Cordillera Andina: la Cordillera de los Andes Peruanos comprende

un conjunto de cordilleras emplazadas entre la fosa Peruano-Chilena y el Llano
Amaznico. Es una cadena montaosa que se extiende a lo largo del borde
Occidental de Sudamrica, desde Venezuela hasta la Tierra del fuego en Chile,
sobre una longitud de 7000 km con espesores que fluctan entre 50 a 70 km
(James, 1978) y alturas mxima de 6000 metros sobre el nivel del mar. Los
Andes peruanos pertenecen al llamado Andes Centrales. En Per, y de norte
a sur, la Cordillera de los Andes se presenta bien definida; sin embargo, es
notoria la presencia de dos inflexiones, la primera a la altura de 5 Sur,
11
denominada deflexin de Huancabamba (D. H.) y la segunda a los 14 Sur,

denomina la deflexin de Abancay (D. A.) (Figura 2.1). La Cordillera Andina
presenta diversas unidades morfo estructurales y anchos que oscilan entre 250
km. en la regin norte y centro de Per hasta 500 km en la frontera entre Per,
Chile y Bolivia (Marocco, 1980; Tavera y Buform, 1998).
2.1.5 La Cadena Volcnica: se encuentra ubicada en la regin sur del Per,
dentro de los 14 hasta los 25 sur en Chile, y se distribuye sobre la cordillera
andina siguiendo un aparente alineamiento con orientaron NO-SE (Figura 2.1).
Las caractersticas geomtricas de cada uno de los volcanes que integran

esta cadena, muestran que la actividad tectnica es contempornea a la
orogenia extensional que experimenta la Cordillera Andina cerca del
Cuaternario Medio y Reciente (Sebrier et al., 1985). Los principales volcanes
presentes en la regin Sur de Per son: Coropuna (6425 msnm), Sabancaya
(5795 msnm), Misti (5825 msnm), Ubinas (5672 msnm), Chachani
(3745msnm), Huaynaputina (4800 msnm), Yucamane (5508 msnm). Es
importante remarcar que en la regin Norte y Centro de Per, la actividad
Volcnica disminuy o desapareci hace 8 millones de aos (Marocco, 1980)
probablemente debido a la subduccin de la Dorsal de Nazca.
2.1.6 Los Sistemas de Fallas: son resultados del constante proceso de

deformacin que sufre la corteza continental a causa del proceso de
subduccin, estos sistemas estn presentes en mayor nmero, de Norte a Sur
sobre la zona Sub-Andina al pie del borde Oriental de la cordillera Andina
(Figura 2.1), afectando los principales plegamientos del cratn brasileo;
dentro de los cuales destacan, los sistemas de fallas de Alto Mayo (AM),
Satipo (SA), Madre de Dios (MD), etc, y a la vez se puede distinguir algunos
sistemas de fallas sobre la Alta Cordillera y el Altiplano destacando entre ellos
los sistemas de fallas de Cordillera Blanca (CB), Huaytapallana (HU),
Tambomachay (TM), etc. En general, todos los sistemas de fallas deben su
origen
una
distribucin
heterognea
de
esfuerzos
tensionales
compresionales en el interior del continente (James, 1978).
12
2.2 ZOIFICACI MORFOLGICA

Este tipo de zonificacin permite distinguir para el Per, cinco zonas
morfolgicas referidas a la estructura y caractersticas del relieve terrestre (Figura
2.2). A continuacin se realiza una descripcin general de cada una de estas zonas
segn Tavera y Buforn (1998).
2.2.1 La Cordillera Costanera: la Cordillera Costanera (Figura 2.2), se origin hace

4500 aos durante las fases tectnicas del Pre-Cmbrico, habiendo sido plegada
durante las orogenias Hercnica (Devnico) y Andina (Cretcico al Plioceno);
adems de ser atravesada por numerosos diques magmticos de composicin variada.
Estos afloramientos estn parcialmente cubiertos por discordancias, rocas
sedimentaras y volcnicas del Terciario. La cordillera de la costa alcanza
elevaciones entre 900 y 1200 metros sobre el nivel del mar. Se presenta segmentada
a lo largo y prximo al litoral siguiendo una direccin NO-SE. El segmento Sur de
esta cordillera, est formada por el macizo de Arequipa ubicado entre los 14 y 18
sur como parte de la deflexin de Abancay y el segmento Norte entre 2 y 7 Sur, el
mismo que forma parte de la deflexin de Huancabamba que se proyecta hacia el
territorio ecuatoriano. Entre las latitudes de 6 y 14 Sur, la cordillera desaparece
posiblemente debido a una posible subsidencia de una parte del margen continental
por la subduccin y migracin lateral de la dorsal de Nazca.
2.2.2 Zona Costera: la Zona Costanera se encuentra entre la Cordillera de la Costa y
la Cordillera Occidental (Figura 2.2), y habindose formado en el Cenozoico. Esta
unidad recorre el territorio peruano de Norte a Sur, es estrecha con elevaciones de 50
a 1500 metros sobre el nivel del mar. Los anchos varan en la regin Norte y Centro,
de 16 a 100 Km.; mientras que, en la regin Sur es ms angosta de 40 a 50 Km.
Durante el Cenozoico, los Andes soportaron gran actividad glaciar que durante la
glaciacin fueron acarreados por los ros de la Vertiente Occidental, rellenando y
formando las pampas costaneras. Despus de esta fase, hasta la actualidad, la
morfologa de la Costa no ha variado en mayor grado.
13
2.2.3 La Cordillera Occidental: la Cordillera Occidental est localizada entre la

zona costera por el Oeste y la Cordillera Oriental (entre 4 y 9 Sur) y con el
Altiplano (entre 9 y 18 Sur) por el Este (Figura 2.2). Formada por elevaciones que
sobrepasan los 5000 metros que van decreciendo conforme se desciende al Sur, con
anchos que varan entre 90 y 100 Km. Esta unidad sigue una direccin NO-SE,
aunque cerca de la deflexin de Huancabamba cambia a direccin hacia NE y en la
regin Sur entre (13 y 14) hacia al Oeste cerca de la deflexin de Abancay.
La Cordillera Occidental est compuesta por rocas volcnicas y plutnicas, del

Mesozoico y Cenozoico, fuertemente plegados y con la presencia de fallas inversas
con grandes escurrimientos. En la regin Sur presenta un alineamiento de conos
volcnicos abajo de 15 latitud Sur, siguiendo una direccin NO-SE. Esta cadena la
integran los volcanes Coropuna, Ampato, Chachani, Misti y Ubinas como los ms
activos.
2.2.4 Zona Intercordillerana: el Altiplano est limitado por Oeste con la

Cordillera Occidental y al Este con la Cordillera Oriental, con elevaciones de 3000
metros sobre el nivel del mar, y anchos variables, as para la regin Norte y Centro
entre 10 a 50 Km y en la regin Sur, es entre 140 a 200 Km. Esta unidad se
extiende de 9 Sur hasta cubrir todo el altiplano peruano-boliviano siguiendo una
orientacin NO a SE (Figura 2.2). Durante el Cenozoico el Altiplano ha sido una
cuenca deposicional de sedimentos clsticos continentales derivados de las
Cordilleras Occidental y Oriental (James, 1978).
14
Figura 2.2: Principales Unidades geomorfolgicos presentes en el Per (Dalmayarac

et al., 1981; Tavera y Buforn, 1998). Los tringulos en blanco indican la presencia de
la cadena volcnica.
15
2.2.5 La Cordillera Oriental: la Cordillera Oriental est formada por rocas

precmbricas y paleozoicas, presenta elevaciones desde 3700 a 4000 metros sobre el
nivel del mar y anchos que varan entre 70 y 100 Km. Aproximadamente, esta unidad
sigue una orientacin NO-SE; sin embargo, soporta un fuerte arqueamiento a 14 Sur
en direccin E-O (Deflexin de Abancay) en la regin Sur. En la Figura 2.2 se
observa que en la regin Norte a ~5 Sur, la Cordillera Oriental desaparece por causa
de la Deflexin de Huancabamba.
La formacin de esta cordillera se inicia durante el tectonismo Hercnico

(Devnico) sobre un basamento o ncleo Precmbrico. El levantamiento fue
controlado por fallas regionales longitudinales. La morfologa de la Cordillera
Oriental ha sido igualmente labrada por actividad glaciar durante el Pleistoceno.
2.2.6 Zona Subandina: ccorresponde a la zona de transicin entre la Cordillera

Oriental y la Llanura Amaznica (Figura 2.2). La unidad presenta plegamientos en
las formaciones sedimentarias Mesozoicas y Terciarias ocurridas durante el
Plioceno las dimensiones de estos segmentos disminuyen en intensidad a medida
que se aproxima a la Llanura Amaznica. Esta unidad se caracteriza por la presencia
de numerosas fallas inversas, sobre escurrimientos y plegamientos de estratos. La
traza de las fallas y los ejes de los pliegues siguen una orientacin de NO-SE.
2.2.7
Llanura Amaznica: la llanura Amaznica se extiende desde la zona
Subandina sobre todo el escudo brasileo, es una amplia llanura compuesta de

importantes secuencias de sedimentos. Las series sedimentarias Mesozoicas y
Cenozoicas son en general planas y slo presentan ligero arqueamiento en el borde
de la zona Subandina, (Figura 2.2).
16
Captulo III: Caractersticas Generales de la Sismicidad en el Per
CAPTULO III
CARACTERISTICAS GEERALES DE LA SISMICIDAD E

EL PER
La constante friccin entre las Placas de Nazca y Sudamericana da origen a
los sismos ms violentos conocidos a lo largo de la historia del Per (Silgado, 1978;
Dorbath et al., 1990). Esta friccin se realiza a lo largo de toda su superficie de
contacto a una profundidad aproximada de 60 Km (Guardia y Tavera, 2010) y sobre
una longitud de 7000 Km aproximadamente, paralelo a la lnea de costa. Sin
embargo los deslizamientos que dan origen a la actividad ssmica solamente
involucran segmentos de longitud menor, permitiendo considerar reas pequeas de
ruptura. Estas reas son conocidas como fuentes sismognicas y sus superficies de
friccin no son visibles en superficie debido a que se encuentran por debajo del
nivel del mar; y sus dimensiones pueden ser estimadas a partir de la distribucin
espacial de las rplicas de los eventos ms importantes (Tavera y Bernal, 2005).
En el Per han ocurrido numerosos eventos ssmicos de gran magnitud, todas

relacionados con el proceso de subduccin entre la placa de Nazca y Sudamericana.
(Tavera y Buforn, 1998; Bernal y Tavera, 2002), siendo la fuente sismognica ms
importante en el Per. La informacin ms completa sobre la actividad ssmica del
periodo histrico (Tabla 3.1) en Per, fue recopilada por Silgado (1978) en su obra,
Historia de los sismos ms notables ocurridos en el Per (1513-1974).
Posteriormente, otros autores como Dorbath et al (1990) contribuyeron a este
estudio. Segn estos investigadores los ms grandes terremotos han producido daos
17
en las ciudades de Tumbes, Trujillo, Lima, Ica, Arequipa, Moquegua y Tacna.

Muchos de estos ocasionaron tsunamis con olas de 15-20 metros de altura que
causaron destruccin en Callao, Pisco, Ilo, Mollendo, Tacna, Arica e Iquique en
Chile; dichos terremotos generaron longitudes de ruptura del orden de 100 a 500 Km.
La informacin ms completa sobre los parmetros hipocentrales de un sismo,

data desde el ao de 1961 aproximadamente, fecha a partir del cual se considera
como periodo Instrumental (Tabla 3.2). La informacin ssmica detallada de este
periodo ha permitido realizar los primeros anlisis y evaluaciones de la sismicidad
que ocurre en Per y as, identificar la presencia de fuentes sismognicas en
superficie y en profundidad. La distribucin espacial de los sismos en el Per es por
dems compleja; sin embargo, es posible identificar la presencia de diferentes
fuentes sismognicas (Bernal, 2000; Tavera y Buforn, 2001).
18
TABLA 3.1
Sismos de gran magnitud e intensidad ocurridos durante el periodo histrico en el
borde occidental del Per.

COORDEADAS
FECHA
MAGITUD
LOGITUD DE
RUPTURA
ITESIDAD
LOCALIDADES
AFECTADAS
#
1
2
3
4
5
6
7
8
9
10
Lon. O
-71,60
-77,00
-70,90
-79,00
-77,60
-76,80
-77,80
-77,50
-77,20
-72,00
Lat. S
-16,60
-12,10
-17,80
-7,90
-12,30
-14,00
-12,30
-13,00
-11,90
-16,50
ao mes da
1582
1
22
1586
7
9
1604 11 24
1619
2
14
1655 11 13
1664
5
12
1678
6
16
1687
9
20
1746 10 20
1784
5
13
Mw
7.5
8.1
8.7
7.7
7.7
7.5
8.0
8.4
8.6
8.4
Ms
7.9
8.1
8.4
7.8
7.7
7.8
8.2
8.4
8.4
Km.
80
175
450
100-150
75
100-150
300
350
300
M.M.
X
IX
IX
IX
IX
X
IX
IX
X
X
Ciudades
Arequipa
Lima
Sur Per
Trujillo
Lima
Ica
Lima Callao
Lima
Lima Callao
Arequipa
11
-72,90
-16,00
1821
10
7.9
7.9
VIII
Arequipa
12
13
14
15
-70,60
-77,70
-74,70
-80,50
-18,30
-11,20
-15,50
-3,80
1868
1940
1942
1953
8
5
8
12
13
24
24
12
8.6
8.1
8.2
6.7
8.2
8.4
6.7
500
180
200
-
X
VIII
IX
VIII
Sur Per
Lima
Ica Arequipa
Tumbes
Mw: Magnitud Momento obtenido por Dorbath et al (1990); Ms: Magnitud de Ondas
Superficiales, obtenido por Silgado (1978), MM; Intensidad en escala de Mercalli
Modificada.
TABLA 3.2
Sismos de gran magnitud ocurridos durante el periodo instrumental en el borde
occidental del Per.
#
1
2
3
4
5
6
7
8
COORDEADAS
Lon. O
-78,60
-78,80
-80,60
-77,80
-79.84
-76,40
-73,75
-76,85
Lat. S
-10,70
-9,20
-4,60
-12,30
-9.80
-15,30
-16,20
-13,49
FECHA
ao
1966
1970
1970
1974
1996
1996
2001
2007
mes da
17 10
31
5
12 10
3
10
2
21
11 12
23
6
15
8
MAGITUD
Mw
7.7
8.0
Ms
7.2
7.9
7.7*
7.7*
8.2*
7.9*
7.5
7.9
8.4
LOGITUD
DE
RUPTURA
ITESIDAD
Km.
100
180
75
140
15
125
350
370
M.M.
VIII
VII-VIII
VIII
VII
VII
VII-VIII
VII
VI-IX
LOCALIDADES
AFECTADAS
Ciudades
3orte de Lima
Huaraz
Tumbes
Lima
Chimbote
3azca
Arequipa
Ica Pisco
Mw: Magnitud Momento obtenido por Dorbath et al (1990); Ms; Magnitud de Ondas
superficiales; Mw* (IGP); MM; Intensidad en escala de Mercalli Modificada.
19
3.1 PERIODO HISTRICO

La sismicidad histrica considera a los sismos ocurridos entre los aos 1500 y
1959, esta informacin data del tiempo de la conquista espaola y colonizacin del
pas, y los datos que se disponen en la actualidad son informaciones recopiladas e
interpretada por Polo (1904), Barriga y Vctor (1951), Silgado (1978), Dorbath et al,
(1990). En general, estos autores revelan que el sismo ms antiguo para el cual se
dispone informacin data del ao 1513. La principal fuente de informacin fue a
travs de crnicas, diversas obras inditas o poco conocidas, crnicas de los
religiosos, narraciones de los viajeros ilustres que visitaron esta parte del continente,
informes administrativos por parte de los gobernantes de la poca. La calidad de
estos datos depender de la distribucin y densidad de la poblacin en las regiones
afectadas por estos movimientos ssmicos; por lo tanto, cabe la posibilidad de que
hayan ocurrido sismos importantes en reas no pobladas o en localidades con las
cuales era difcil establecer comunicacin (Bernal y Tavera, 2002). A continuacin,
se describe los aportes ms resaltantes de los investigadores que se dedicaron a
estudiar la sismicidad histrica de Per.
T. Polo (1904). El historiador Toribio Polo recalca que se haban
producido ms de 2500 temblores en el territorio peruano, desde la
conquista hasta los finales del siglo XIX. Asimismo, indica que los
grandes sismos produjeron cuantiosos daos materiales debido a que las
construcciones de aquellos aos eran inadecuadas para resistir los
movimientos violentos del suelo.
P. Barriga (1939). El padre Barriga describe las caractersticas de todos
los sismos de mayor magnitud ocurridos en la regin Sur de Per hasta
1868. Dicha informacin servira como base para otros investigadores.
E. Silgado (1978). Recopila la informacin de los sismos importantes
ocurridos en Per entre los aos 1513 y 1974. Describe las caractersticas
de los sismos detalladamente consignando los valores de intensidad local
y regional producidos por estos sismos. Asimismo, el autor estima la
magnitud de los sismos a partir del rea de intensidad mxima a fin de
poder compararlos con sismos ms recientes. El trabajo de este autor se
20
constituye como la recopilacin ms completa de sismos importantes

ocurridos en Per entre 1513 y 1974.
Dorbath et al (1990). Recientemente estos autores realizaron una revisin
detallada de toda la informacin existente sobre sismos histricos a fin de
correlacionar las reas de intensidad mxima y las longitudes de ruptura
de los grandes sismos y a partir de esta, estimar su magnitud (Mw). Los
resultados obtenidos por estos autores, son comparados con los
respectivos valores obtenidos en otras investigaciones aplicando
diferentes mtodos. Asimismo, los autores sugieren que los grandes
sismos en la regin sur del Per tienen un periodo de recurrencia del
orden de una centuria (100 aos).
Entre los sismos ms importantes ocurridos durante el periodo histrico en la
regin Norte fueron los de 1619 y 1953 (VII MM), los mismos que produjeron
muerte y destruccin en las ciudades de Trujillo y Tumbes. En la regin Central
sobresalen los sismos ocurridos en los aos 1586 (IX MM), 1687 (VIII MM) y 1746
(X MM), los mismos que destruyeron casi completamente la ciudad de Lima. El
sismo de 1746 gener un tsunami con olas de 15-20 metros de altura que inundo
totalmente al puerto del Callao. En la regin Sur, ocurrieron sismos importantes en
los aos de 1604 (IX MM), 1748 (X MM) y 1868 (X MM) que destruyeron
principalmente a las ciudades de Arequipa, Moquegua, Tacna, Puno y Norte de
Chile. El terremoto de 1868 habra producido un tsunami con olas de 12 a 16 metros
de altura. La profundidad focal del total de los sismos histricos no ha sido
determinada con precisin; sin embargo, debido a su ubicacin geogrfica entre la
fosa y la lnea de costa y por los daos causados en superficie, estos fueron
considerados en su mayora como sismos superficiales.
En la Figura 3.1 se muestra la localizacin epicentral de los sismos histricos

ocurridos entre los aos 1582 y 1959, en el borde occidental del Per. Estos sismos
han presentado magnitudes mayores a 6,6Mw y muchas de ellas produjeron tsunamis
con olas de diferentes alturas. En esta figura se observa que los sismos se distribuyen
principalmente entre la lnea de fosa y la costa, localizndose en mayor nmero en la
regin Centro y Sur de Per.
21
Figura 3.1: Distribucin espacial de grandes terremotos ocurridos en el borde occidental

del Per, dentro el periodo histrico (1471-1959) segn Silgado (1978). Las barras de
color azul indican los sismos que produjeron tsunami.
22
3.2 PERIODO ISTRUMETAL

La informacin sismolgica instrumental del Per se encuentra recopilada en
distintos Catlogos Ssmicos elaborados por el Instituto Geofsico del Per (Ocola,
1984; Huaco, 1986; Tavera y Agero 2000). El clculo de los parmetros
hipocentrales (Latitud, Longitud, Profundidad y Magnitud) de los sismos ocurridos
durante este periodo son ms precisos, y las de mayor magnitud fueron estudiados
fueron estudiados por diversos investigadores como: Stauder, (1975); Barazangi y
Isacks, (1976); Hasegawa y Sacks, (1981); Grange et al., (1984); Bevis y Isacks,
(1984); Schneider y Sacks (1987); Rodriguez y Tavera, (1991), Cahill y Isacks,
(1992); Tavera y Buforn, (1998); Bernal, (2000) utilizando diferentes bases de datos;
telessmicos y locales.
Los resultados obtenidos por estos investigadores han permitido conocer las
caractersticas espaciales, en superficie y profundidad de la sismicidad asociadas a la
deformacin cortical y al proceso de subduccin de la placa de Nazca bajo la
Sudamericana en el Per. Recientemente, Tavera y Buforn, (2001), as como Bernal
y Tavera, (2002), concluyen en que la distribucin espacial de los sismos en el Per
es por dems compleja; no obstante, se pueden identificar reas y/o agrupaciones de
sismicidad constituidas por importantes fuentes sismognicas. Durante el periodo de
1960-2010 en el borde Occidental del Per han ocurrido siete terremotos de
magnitudes elevadas en las regiones Centro y Sur (Figura 3.2).
En 1966 se produce un terremoto con magnitud de 7.5Mw con epicentro en
la parte Central del Per, y present una longitud de ruptura de 100 Km y a su
vez gener un tsunami con olas de altura de 1.6m.
Las intensidades
evaluadas en la escala de Mercalli Modificada (MM) fue de VIII en la regin

de Huacho, y VII en Lima.
El 31 de mayo de 1970 vuelve a ocurrir un terremoto con epicentro frente a la

costa del departamento de Ancash con magnitud 7.8Mw. Este sismo fue uno
de los ms destructores y lleg a producir causando aproximadamente 70,000
muertos y 50,000 desaparecidos (Ericksen et al., 1970).
23
Figura 3.2: Distribucin 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.
24
El terremoto provoc un alud como consecuencia del desprendimiento de

una parte del nevado de Huascaran y enterr las localidades de Yungay y
Ranrahirca. El tamao fue similar al terremoto de 1966, pero ms destructivo
probablemente porque su localizacin epicentral estuvo
dentro del
continente.
El 3 de octubre de 1974 ocurre otro terremoto de magnitud 7.5Mw con

epicentro a 80 Km al oeste de Lima. Este terremoto caus severos daos en
Lima y localidades cercanas con intensidades mximas de IX (MM) en
algunos lugares de Lima. En distritos al sur de Lima la intensidad no excedi
los VIII (MM). El sismo produjo un tsunami local de con olas de 3 metros.
El 21 de febrero de 1996, se produjo un sismo cuyo epicentro fue localizado a

210 Km. al SO de la ciudad de Chimbote, con magnitud 7.5Ms. Este sismo
produjo un tsunami que caus daos materiales y prdidas de 15 vidas
humanas en el departamento de Chimbote. En Salaverry causando daos
materiales de poca consideracin.
En la regin central del Per, el 12 de noviembre de 1996 ocurre un sismo de

magnitud 7.7Mw, con epicentro a 98 Km al SO de la ciudad de Nazca,
Departamento de Ica (IGP, 1996). Este sismo present su epicentro en la
misma rea donde ocurri el sismo del 24 de agosto de 1942 (Figura 3.1), y
pondra en evidencia el rol importante que cumple la Dorsal de Nazca en la
generacin de sismos en esta regin. El sismo produjo intensidades mximas
de VII (MM) sobre un rea elptica de 4000 km2 en torno a la ciudad de
Nazca.
El da 23 de Junio de 2001, a horas 15 y 33 minutos (Hora Local), la regin

Sur del Per fue nuevamente testigo de la fuerza con que la naturaleza nos
indica que la Tierra, sobre la cual habitamos, es dinmicamente activa. Ese
da se predijo la ocurrencia de un terremoto de magnitud 8.0 Mw que produjo
muerte y destruccin principalmente en los departamentos de Arequipa,
Moquegua y Tacna en Per, Arica e Iquique en Chile y La Paz en Bolivia.
25
Las intensidades mximas en algunas localidades y ciudades, prximas a la

costa, llegaron a VII-VIII en la escala Mercalli Modificada. Este terremoto se
constituye como el ms grande ocurrido en esta regin desde el ao de 1868.
El 15 de Agosto de 2007 la zona sur de la regin central de Per fue afectada

por un terremoto de magnitud 8.0Mw que en algunos segundos produjo
muerte y destruccin en las ciudades de Pisco, Ica y Chincha en donde se
evaluaron intensidades mximas de VII-VIII en la escala de Mercalli
Modificada (MM). El terremoto produjo intensidades de V (MM) en la
ciudad de Lima y fue sentido hasta las ciudades de Piura, Arequipa, Cusco y
Pucallpa con intensidades del orden de II-III (MM). Este sismo tuvo como
caractersticas principales su gran duracin y el aparente proceso complejo de
ruptura que experimento, para luego ser seguido por un gran nmero de
rplicas con magnitudes que no sobrepasaron el grado 6.5Mw. El terremoto
del 15 de Agosto se constituye como el de mayor magnitud y duracin
ocurrida en esta regin en los ltimos 50 aos. El rea de ruptura fue de
170x120 Km con el eje mayor paralelo a la lnea de costa.
26
3.3 DISTRIBUCI ESPACIAL

Teniendo en cuenta la profundidad de los focos, los terremotos pueden
clasificarse en:
Terremotos de foco Superficial : (h 60 Km.)

Terremotos de foco Intermedio : (61 < h 300 Km.)
Terremotos de foco Profundo
: (h >300 Km.).
3.3.1 Sismicidad con foco Superficial

Los epicentros de estos sismos se localizan principalmente entre la fosa y la
lnea de costa, (Figura 3.3a, crculos rojos), y en ella se considera la ocurrencia de
terremotos de magnitud elevada, como los del 17 de octubre de 1966 (7.7Mw), 31 de
mayo de 1970 (7.8Mw), 3 de octubre de 1974 (7.9Mw), 23 de junio de 2001
(8.0Mw), 21 de febrero de 1996 (7.5Mw), 12 de noviembre de 1996 (7.7Mw) y 15 de
agosto de 2007 (8.0Mw).
En el interior del continente se tiene otro grupo importante de sismos con foco
superficial, se distribuyen, de norte a sur, a lo largo de la Cordillera Oriental y en el
margen Occidental de la zona Subandina (3-13 S). Estos sismos pueden ser
asociados a las deformaciones corticales que dieron origen a los principales sistemas
de fallas presentes en el Per, tales como: la falla de Moyobamba, Huytapallana,
Ayacucho, Tambomachay, Cordillera Blanca, Marcota, etc. (Tavera y Buforn, 1998).
En la zona del altiplano, tambin existen sismos superficiales, pero en menor nmero
y ms dispersos.
27
3.3.2 Sismicidad con foco Intermedio

Estos sismos presentan sus epicentros principalmente el interior del continente
casi en su totalidad (Figura 3.3a, crculos verdes): el primer grupo se encuentra
ubicado paralelo a la lnea de costa por debajo de los 8 S y ha dado origen a sismos
que frecuentemente fueron sentidos en superficie (15 de enero de 1958, 6.5Ms; 28 de
septiembre de 1968, 6.0Ms; 1 de enero de 1974, 6.4Ms; 16 de febrero de 1979,
6.9Ms; 28 de marzo de 1982; 6.1Ms; 18 de abril de 1993, 6.0Ms). Como segundo
grupo se tiene los sismos con epicentros en el interior del continente y a lo largo de la
Cordillera Oriental y zona Subandina de las regiones norte y centro. Estos sismos
raramente producen daos en superficie. El tercer grupo de sismos se encuentra en la
regin sur del Per sobre la Cordillera Occidental y el Altiplano, siendo esta regin
la de mayor ndice de sismicidad a estos niveles de profundidad (Tavera y Buforn,
1998).
3.3.3 Sismicidad con foco Profundo

El mayor nmero de sismos con foco profundo ocurren en la regin central del
Per (Figura 3.3a, crculos azules), limite con Brasil entre las latitudes de 6 a 11 S
originando una distribucin de norte a sur sobre una longitud de 500 Km de
profundidad; mientras que, en la regin sur, frontera con Bolivia, el nmero de
sismos es menor. La mxima profundidad de estos sismos es del orden de 550 a 670
Km. Prcticamente no son sentidos en superficie. Aunque su origen todava tiene sus
controversias, se considera que ellos tendran su origen en transformaciones de fase
de olivino a espinel, y/o a un trozo de litosfera fra que flota en el interior del manto.
3.4 PERFILES SSMICOS

Con la finalidad de analizar la distribucin de la sismicidad en funcin de la
profundidad de sus focos, se ha elaborado tres perfiles ssmicos perpendiculares a la
fosa, para las regiones norte (perfil AA`), centro (perfil BB) y sur (perfil CC) del
Per.
28
3.4.1 Regin &orte: la Figura 3.3b, corresponde al perfil elaborado para la regin
Norte de Per, y en ella se observa que los focos aumentan su profundidad conforme
se distribuyen de Oeste a Este. Esta sismicidad se inicia en la fosa y se distribuye con
alineamiento con una inclinacin de 10-15 a 30 aproximadamente hasta alcanzar
una profundidad de 200 km. Para luego hacerse horizontal hasta una distancia de
870 km desde la fosa (F). En esta regin se observa una disminucin del nmero de
los sismos entre 320 y 480 km de distancia desde el punto A, sugiriendo que esta
rea puede corresponder a una zona assmica. La distribucin de los sismos en
funcin de la profundidad de sus focos, describe la subduccin de la placa ocenica
bajo la placa continental.
3.4.2 Regin Centro: el perfil para esta regin (Figura 3.3c), muestra que los sismos
tienen la misma distribucin en profundidad que los de la regin norte; sin embargo,
los sismos siguen una pendiente que va entre los 25 a 30 aproximadamente y a
partir de los 450 km de distancia desde el punto B, los focos de los sismos siguen una
distribucin horizontal hasta una distancia de 500 km aproximadamente. Asimismo,
en este perfil se observa la presencia de sismos con foco profundo, los mismos que se
distribuyen entre 550-670 km de profundidad y corresponden a los sismos que
ocurren en el lmite Per Brasil.
3.4.3 Regin Sur: para la regin sur de Per (Figura 3.3d), se observa la presencia
de un mayor nmero de sismos, que se distribuyen en profundidad siguiendo una
pendiente que va entre los 25 a 30 hasta alcanzar 300 km de profundidad. Entre
300 y 500 km de profundidad, existe ausencia de sismos. Estn ellas presentes a la
distancia de 650 metros de distancia del punto C (nivel de referencia) y a una
profundidad de 550 a 700 km respectivamente.
29
Figura 3.3: a) Mapa de sismicidad del Per para el periodo 1964-2008 (Mw4)
(Tavera y Bernal, 2009). Los crculos rojos indican los sismos con foco superficial (h 60
Km.), crculos verdes sismos con foco intermedio (61 < h 300 Km.), crculos azules sismos
con foco profundo (h >300 Km.). ). En ella se muestra los perfiles para las los regiones
3orte (AA), Centro (BB) y Sur (CC) del Per. a) Regin 3orte, b) Regin Centro, c)
Regin Sur. Donde F= Fosa. LC= Lnea de Costa. En todo el borde occidental del Per.
30
3.5 GEOMETRIA DE LA SUBDUCCI
La subduccin de placas es un proceso presente en los lmites convergentes de

una placa ocenica bajo otra continental, puesto que la primera es de mayor peso
especfico y ms fra y subduce bajo la segunda de menor peso especfico debido a su
mayor grosor cortical. El plano de friccin entre las dos placas que convergen se
llama plano de Benioff, y es en donde se concentran los focos o hipocentros de los
sismos. Debido a los diferentes niveles de velocidad y caractersticas dinmicas de la
tectnica de placas de cada regin, la geometra de subduccin es diferente en
funcin de la variacin en el ngulo de subduccin a lo largo del borde occidental del
Per.
Los estudios realizados para conocer el modelo la geometra del proceso de

subduccin son diversos, pero en conjunto han permitido tener una idea general
sobre la forma como se realiza este proceso. Asimismo, la distribucin espacial de
sismos con origen en la friccin de placas y la deformacin interna de la placa de
Nazca, han permitido a muchos autores configurar la geometra de la misma
utilizando datos telessmicos, regionales y locales. (Grange et al, 1984; Schneider y
Sacks, 1987; Rodrguez y Tavera, 1991; Cahill y Isacks, 1992; Tavera y Buforn,
1998; Bernal, 2002).
Los estudios realizados por Tavera y Buforn (2001) y Bernal (2002) han
mostrado, en el caso del Per, que la subduccin de la placa de Nazca por debajo de
la Sudamericana (Figura 3.4) se produce siguiendo dos tendencias. La primera
muestra que la placa de Nazca, en las regiones norte y centro, subduce con un ngulo
medio de 30 hasta una profundidad de 120 Km a partir de la cual se hace horizontal
hasta distancias de 700 Km desde la lnea de fosa. En la regin sur, la placa subduce
de manera contina hasta una profundidad de 250 Km. El cambio observado entre
ambos estilos de subduccin se debe a que la placa soporta una contorsin a altura de
la latitud de 16 Sur.
31
Figura 3.4: Esquema 3D que muestra la figura de la geometra de la placa de 3azca.

Bernal y Tavera, (2002).
Uno de los factores que influye en la geometra de la zona de subduccin es la

presencia de la dorsal de Nazca en el Sur del Per, la dorsal de Carnegie en el Sur de
Ecuador y posiblemente el Plato Inca propuesto por Gutscher et al., (1999) y
Hampel, (2002). La dorsal de Nazca en su impacto (11.2 Ma a 11 Latitud Sur.) y su
migracin, dio como resultado la subduccin subhorizontal hasta su posicin actual
entre 14 y 16 Latitud Sur (en el Norte-Centro del Per) y Normal (en el Sur). Al
sur de 17 adems, la dorsal es el lmite de la actividad volcnica en el Sur y la
deflexin de Abancay; as como la Dorsal de Carnegie lo es en el Sur de Ecuador
(Hampel, 2002).
Gutscher et al. (1999), mostraron la posible existencia de una meseta

actualmente subducida debajo del Nororiente Peruano, denominada Meseta Inca
(Figura 3.5a y 3.5b), basndose en el anlisis de la sismicidad y la reconstruccin de
la Placa Faralln. Los autores dedujeron que esta meseta Inca junto con la Dorsal de
Nazca y la Dorsal de Iquique son un espejo con respecto a la meseta Marquezas, la
Dorsal de Tuamoto y la Dorsal Austral, todas ellas pertenecientes a la Placa del
Pacfico respectivamente.
32
Los autores enuncian que estas estructuras probablemente produjeron una

subduccin tipo flat slab en la Placa de Nazca (entre 2.5S y 12S), ya que al ser
menos densas que la corteza ocenica producen un levantamiento del slab (Fig.
3.5b), y que adems coincide con la ausencia de vulcanismo. La Dorsal de Nazca y la
Meseta Inca comenzaron su subduccin en forma oblicua durante el Mioceno medio
(11 15 Ma) con una tasa de 7 a 7.5 cm/ao (Hampel, 2002).
Rosenbaum et al. (2005), mostraron una posible relacin entre la generacin

de un gran nmero de depsitos metalognicos de los Andes Peruanos durante el
Mioceno (15 - 5 Ma) y la subduccin de la Placa de Nazca bajo la Placa
Sudamericana.
La
distribucin
espacial
temporal
de
estos
depsitos
corresponderan con el arribo de la Dorsal de Nazca en el centro del Per y la

actualmente subducida meseta Inca en el norte del Per (Flat slab). Al desarrollar la
reconstruccin de la Placa de Nazca, los autores muestran la existencia de una rpida
reaccin metalognica originada por la subduccin de estas anomalas topogrficas,
reflejada en la generacin de depsitos minerales hidrotermales en las proximidades
de estas estructuras. Los autores evaluaron 382 localidades metalognicas formados
desde hace 35 Ma, finalmente excluidas todas las zonas cuyas edades de
mineralizacin no fueran Miocnicas. La distribucin temporal de estos depsitos se
muestra claramente que tiene un inicio Miocnico (15 Ma), Por ejemplo Yanacocha
cuya edad de mineralizacin es de 13 Ma; el hecho de que no hay una mineralizacin
ms joven en la zona de flat slab desde 10 Ma se podra explicar porque en ese
momento la Meseta Inca se encontraba totalmente en subduccin.
33
a)
b)
Figura 3.5: a) Mapa de las lneas de contorno del slab, modificado de los modelos de
Grange et al. (1984) para el sur y para el norte y centro del Per. Los tringulos rojos
indican las zonas de vulcanismo activo y los blancos representan volcanes extintos. b)
Modelo 3D de la zona de subduccin elaborado por Gutscher et al., (1999), mostrando las
influencias de la dorsal de 3azca y el Plato Inca perdido.
34
Captulo III: Caractersticas Ssmicas del Per
3.6 DISTRIBUCI DE REAS DE RUPTURA

La constante ocurrencia de sismos de gran magnitud frente al borde oeste de
Per, y el conocimiento de las dimensiones de sus reas de ruptura, ha permitido
identificar la presencia de zonas con diferentes caractersticas ssmicas. En la Figura
3.6 se presenta la distribucin de las reas de ruptura de los grandes sismos ocurridos
en todo el borde oeste de Amrica del Sur, para los siglos XIX, XX y XXI.
En el siglo XIX los ms grandes sismos ocurrieron en la regin sur del Per
(1868) y norte de Chile (1877), ambos con magnitudes Mw de 9.0. Hacia el norte y
sur de ambas reas de ruptura es notoria la ausencia total de sismos. Durante el siglo
XX, todos los sismos ocurrieron al norte y sur de las reas involucradas en los sismos
de 1868 y 1877, sobresaliendo el gran sismo de Chile de 1960, quizs el mayor para
el cual se tiene informacin ssmica instrumental (magnitud Mw de 9.5). Tambin es
notoria la presencia de una importante ausencia de sismos en la regin norte del Per.
En el siglo XXI se han producido 3 grandes sismos, as se tiene hacia el sur del Per
el terremoto de Arequipa del 2001 (7.9Mw), por sus caractersticas y longitud de
ruptura (Tavera y Bernal, 2005) no es proporcional al ocurrido en 1868.
Seguidamente y despus de seis aos se produce otro gran sismo en Pisco (8.0Mw),
y seguidamente frente a las costas de la regin sur de Chile ocurre el terremoto de
Biobio del 2010 (8.8Mw).
De acuerdo a lo escrito, las regiones centro y sur del Per son las de mayor
potencial ssmico en razn de que en ambas, los sismos parecen repetirse con
perodos de tiempo comprendidos mayores a 150 aos. Contrariamente, es posible
que en la regin norte sean necesarios perodos de acumulacin de energa mayores a
500 aos para que se produzca un gran sismo.
35
Figura 3.6: Distribucin de las reas de ruptura de grandes sismos ocurridos en el borde
Oeste de Sudamrica durante el siglo XIX, XX, XXI (Tavera y Bernal, 2005).
36
Una distribucin temporal ms detallada de los terremotos ocurridos en Per ha

sido representada en funcin de su longitud de ruptura y de su ubicacin en latitud vs
perodo de ocurrencia. En dicha figura puede observarse claramente el diferente
potencial ssmico de cada regin.
Regin &orte: en esta regin (Figura 3.7), solamente se tiene informacin
sobre la ocurrencia del sismo de 1619 (Trujillo) y despus de esto en la
regin norte no se ha vuelto a experimentar la ocurrencia de un nuevo
terremoto de gran magnitud, esta regin presenta un silencio ssmico hasta el
presente de 392 aos.
Regin Central: la actividad ssmica en la parte central del Per es bastante

compleja debido a la irregularidad de las longitudes de ruptura de los sismos
ocurridos en los ltimos 70 aos, En esta regin se identifica dos
agrupamientos de sismos que han generado longitudes de ruptura entre 100 y
500 Km (Figura 3.7). El primer grupo considera los sismos ocurridos entre
los aos 1650 y 1750, para luego aparecer un perodo de silencio ssmico del
orden de 190 aos que se rompe con la ocurrencia de 6 sismos hasta el ao
2007, a partir del cual volvera a iniciar el perodo de acumulacin de energa.
Regin Sur: los terremotos en esta regin definen dos agrupamientos, el

primero considera los terremotos ocurridos hasta el ao 1604 y a partir de esta
fecha se observa un periodo de silencio ssmico de 180 aos (Figura 3.7), para
luego ocurrir hasta dos grandes terremotos en los aos 1784 y 1868. Los
sismos estuvieran precedidos por otros de menor magnitud. As, antes del
sismo de 1604 (450 km. de longitud de ruptura) ocurri otro menor en 1582;
antes del sismo de 1868 (500 km. de longitud de ruptura) se produjo otro, en
1784, tambin de menor magnitud y finalmente, en el ao 2001 se tuvo lugar
otro sismo pequeo que podra ser el precursor de otro mucho mayor.
37
En esta figura tambin puede observarse que existen perodos en los cuales se
producen sismos en la regin sur y otros en los cuales la ocurrencia es nula, siendo
esta distribucin contraria a la que se observa en la regin central.
Figura 3.7: Grfico que muestra la conclusin de la longitud de ruptura de los grandes e
histricos sismos de subduccin registrados en el borde occidental del Per vs tiempo. Las
barras de color rojo indican la estimacin de la longitud de ruptura, parte superior de cada
barra se muestra el ao que ocurri el terremoto. Grafico Modificado de Dorbath et al,
(1990).
38
Captulo IV: Modelos de Procesos de Ruptura
CAPTULO IV
MODELOS DE PROCESOS DE RUPTURA

El proceso de ruptura de un sismo es por lo ms heterogneo y por lo tanto,
difcil de describir. Las caractersticas de la ruptura de los sismos han sido estudiadas
y discutidas por muchos investigadores, quienes han propuesto diferentes modelos
para explicar su gnesis.
El primer modelo mecnico para explicar la ocurrencia de terremotos, fue
presentado por Reid en 1911 para poder explicar la gran fractura producida por el
terremoto de San Francisco de 1906. Este modelo se conoce como la teora del rebote
elstico y propone que los terremotos ocurren en regiones sujetas a deformaciones
debida a causas externas, generalmente esfuerzos tectnicos regionales provocados
por el movimiento de placas litosfericas. Cuando los esfuerzos acumulados en un
lugar superan la resistencia de la roca se produce su fracturacin y/o un
desplazamiento en zonas de debilidad preexistentes. La zona de facturacin se
denomina falla ssmica. La deformacin producida por los esfuerzos tectnicos es
debida a la dinmica terrestre.
La teora del "rebote elstico" (Reid, 1911), que est ilustrada en la Figura 4.1,
establece que existen ciertas zonas preferenciales de la corteza terrestre (Figura 5.1a)
donde se van acumulando lentamente grandes esfuerzos que son soportados por los
materiales (rocas) que la constituyen. Estos esfuerzos ocasionan en las rocas
39
deformaciones elsticas cada vez mayores (Figura 5.1b) hasta que se supera la
resistencia de las mismas (Figura 5.1c), y se produce entonces una liberacin casi
instantnea de la energa acumulada a travs del tiempo.
El resultado de este mecanismo es la propagacin de la energa liberada, en

forma de ondas ssmicas y el retorno a un estado de equilibrio elstico (Figura 4.1d)
de la zona previamente sometida a esfuerzos, con la presencia de una fractura o falla
geolgica, muchas veces visible en la superficie de la tierra. Este modelo mecnico
para explicar el origen de los terremotos fue aceptado inmediatamente, pero qued
sin aclarar el por qu de la existencia de zonas preferenciales de concentracin de
esfuerzos.
Figura 4.1: Representacin esquemtica del modelo de rebote elstico de Red (1911). a)
Bloque en estado de reposo y acumulacin de esfuerzos, b) Etapa de deformacin debido a
los esfuerzos acumulados, c) Instante de la ruptura, d) Rebote elstico para una nueva
posicin de equilibrio.
40
4.1 MODELOS DE FUETE SISMICA
A la ocurrencia de un terremoto, la propagacin del frente de ruptura sobre el

plano de falla, depende principalmente del estado de esfuerzos a la que se encuentra
sometida la zona a fallar. Para ello se han desarrollado varios modelos para explicar
la complejidad de los procesos de ruptura: as tenemos el modelo de Barreras (Dass y
Aki 1977), y el de Asperezas (Kanamori, 1981). A continuacin se describe cada uno
de ellos.
4.1.1 MODELO DE BARRERAS
Las Barreras son zonas de las fallas que oponen una fuerte resistencia a la
ruptura y pueden as reducir la velocidad de propagacin de la ruptura o incluso
detenerla completamente, (Dass y Aki 1977). Las barreras constituyen un modelo
sencillo de introducir una cierta heterogeneidad a la ruptura y explicar la complejidad
de los sismogramas observados.
El modelo de Barreras puede ser mejor explicado a partir del estado de
esfuerzos antes y despus de ocurrir un terremoto, as como el deslizamiento que se
produce. De acuerdo a la Figura 4.2a, y suponiendo que las zonas que se han roto
durante el terremoto son dos, (representadas por las reas de color blanco) en los
perfiles adyacentes, se puede observar que los esfuerzos antes del terremoto son
uniformes y si hay variacin esta es pequea y relativamente importante. Durante el
terremoto, la ruptura se produce en las reas sin sombrear y en ellas los esfuerzos se
relajan y caen hasta llegar a un valor cero, producindose una transferencia de
esfuerzos hacia los bordes y uniones entre las fallas, en donde aumenta la tensin
para impedir que la falla se deslice. En el caso del deslizamiento, esto solo se
produce en las dos reas sin sombrear sobre la falla, esto debido a que las Barreras
reducen enormemente el movimiento de sus compartimentos. De manera general este
modelo, describe mejor los terremotos asociados a deformaciones corticales (fallas),
donde se produce eventos ssmicos con magnitud moderada (6.0).
41
4.1.2 MODELO DE ASPEREZAS

Se define como el lugar donde se inicia el proceso de ruptura en un medio muy
heterogneo, debido a que ya hubo terremotos anteriormente, entonces esta zona se
encuentra libre de esfuerzos (Tavera, 1992).
En este modelo, el estado inicial es una gran falla que ha sufrido una serie de
terremotos, pero que en su interior aun queda una serie de reas que aun no se han
roto. En la Figura 4.2b, se puede ver que toda la zona ha liberado energa (zonas
blancas), con los terremotos precedentes, quedando solo dos reas que no se han roto
(zonas oscuras). El nuevo terremoto se iniciar en estas dos reas o en cadena. En
trminos de esfuerzos, la situacin inicial es muy diferente al modelo de Barreras;
casi la totalidad de la falla ya liber energa y los esfuerzos se concentran en las
reas resistentes. Durante los terremotos posteriores, esta energa se liberar
totalmente y la situacin final es de esfuerzos casi uniformes. El deslizamiento
producido en este modelo es similar al de Barreras; en las reas de cada de esfuerzos
apenas si hay deslizamiento (Tavera, 1992).
Los grandes terremotos ocurridos en el periodo histrico e instrumental en el

borde occidental del Per, fueron debidos principalmente al proceso de subduccin y
segn los modelos expuestos obedeceran al modelo de asperezas.
42
Figura 4.2: a) Modelo de Barreras (Das y Aki, 1977), b) Modelo de asperezas (Kanamori, 1981), donde D= Deslizamiento,
= Esfuerzo. Modificado por Tavera, 1992.
43
El modelo de asperezas, es atribuido a la presencia de los gap ssmicos, los

mismos que se encuentran bien representados en las zonas de subduccin (Kanamori,
1981). En la Figura 4.3 se muestra un esquema en 3D para la zona de subduccin del
Per, en donde las Asperezas estn representadas por las reas de color rojo, las
mismas que indican reas con acumulacin de energa elstica. De manera detallada
en la Figura 4.3 se esquematiza el mecanismo de ruptura que se desarrolla en una
zona de subduccin; es decir, la zona de interaccin de la placa ocenica (Nazca) con
la placa continental (Sudamericana). En un tiempo (t0) en esta superficie se producen
terremotos de diversa magnitud (estrellas de color negro), todas distribuidas de
manera heterognea sobre la superficie (rea de color amarillo) restando reas que no
se han roto (reas de color rojo) y que actan o se comportan como Barreras. Pasado
un tiempo (t0 + t1), las Barreras llegan a comportarse como Asperezas, siendo estas,
las mismas fuentes de terremotos (estrellas de color verde), completndose de esta
manera, el ciclo de los terremotos.
Figura 4.3: Esquema en 3D donde se muestra la distribucin espacial de asperezas en

una zona de subduccin.
44
El mecanismo de ruptura de un terremoto, como se ha descrito, es muy

complejo y va depender principalmente del estado de esfuerzos y de las propiedades
fsicas del medio que est prximo a fracturarse. Los modelos de Asperezas y
Barreras explican mejor esta complejidad, el primero, asociado a los mecanismos de
ruptura de los terremotos que ocurren en las zonas de subduccin y el segundo, a los
terremotos asociados a fallas.
Los grandes terremotos ocurridos en el periodo histrico e instrumental en el

borde occidental del Per, fueron debidos principalmente al proceso de subduccin y
segn los modelos expuestos obedeceran al modelo de Asperezas. Los terremotos
ocurridos en el interior del continente, asociados a fallas geolgicas activas,
obedeceran al modelo de Barreras, sin embargo; estos son menos destructores a
comparacin de los que ocurren por el proceso de subduccin.
45
Captulo V: Metodologa para Clculo del Parmetro b
CAPTULO V
METODOLOGA PARA EL CLCULO DEL PARAMETRO b

La superficie de contacto entre dos placas no se encuentra en un estado
uniforme de esfuerzos, sino que existe una extensa liberacin de los mismos en
forma de calor y energa en algunas partes de la superficie de acoplamiento de
placas. Los sismos comprometen parte de estas superficies y dejan otras libres
llamadas asperezas. El terremoto debe originarse en una zona donde no se han
producido eventos por un tiempo, lo que describe la presencia de una aspereza o de
una zona con gran acumulacin de energa. A raz de estas observaciones Wiemer y
Wyss (1997) proponen una metodologa netamente probabilstica para identificar la
presencia
y la ubicacin geogrfica de zonas donde podra estar generando el
proceso de acumulacin de energa.
Este mtodo consiste en calcular los valores de a y b utilizando la

distribucin frecuencia-magnitud de los sismos Log"=a-bM, (Ishimoto e Iida, 1939;
Gutenberg y Richter, 1944), (Figura 5.1). El parmetro crtico de esta relacin viene
a ser
el valor de b que define la pendiente de la recta para la distribucin
frecuencia-magnitud de los sismos y que fsicamente caracteriza los niveles de

esfuerzos acumulados y liberados en una zona, y puede ser calculado utilizando una
base de datos ssmicos y el ajuste lineal de mnimos cuadrados o a partir del mtodo
de mxima verosimilitud (Aki, 1965; Utsu, 1965). Empleando los algoritmos
desarrollados en el programa
ZMAP (Wiemer y Zuiga, 1994), se calcula los
valores de los parmetros a y b para reas pequeas como parte de otra mayor
46
que previamente fue descretizada utilizando una malla, cuyos espaciamientos

dependern de la densidad de datos del catalogo ssmico y de las dimensiones del
rea de estudio.
Figura 5.1: Relacin de la distribucin Gutenberg- Richter que muestra el logaritmo del
nmero acumulado de sismos y el trazo de la pendiente, para identificar la magnitud de
completitud (Mc). Los cuadrados indican la frecuencia acumulativa; mientras que, los
tringulos la frecuencia no acumulativa.
5.1 DISTRIBUCIN FRECUENCIA-MAGNITUD
La distribucin de los sismos con respecto a sus magnitudes, muestra una auto
semejanza, esto parece ser el resultado de un proceso auto-similar y obedece a una
ley potencial o a una escala fractal. Esto implica la ausencia de caractersticas
asociadas al tamao del evento (limites tericos sobre el mximo o mnimo tamao
del evento). Por consiguiente, la frecuencia de los terremotos frente a su tamao (N),
tiende a tener una forma lineal fija.
47
Las siguientes frmulas empricas son conocidas en el Asia como la relacin

Ishimoto-Iida (Kulhanek, 2005, Ishimoto e Iida, 1939) y en, Occidente como la
relacin de Gutenberg-Richter (Kulhanek, 2005, Gutenberg y Richter, 1942) y est
definida como:
Log N= a bM
N= 10
a-bM
(5.1)
(5.2)
Donde:
- N es el nmero de sismos con una magnitud igual o mayor que M sobre un
periodo de observacin.
- a y b son constantes reales positivas.
La ecuacin (5.1) define la distribucin de los eventos con respecto a la

magnitud, para una cierta regin y un intervalo de tiempo. La ecuacin (5.2) tambin
usualmente es referida como la relacin magnitud-frecuencia (RMF) (Kulhanek,
2005), y representa una relacin lineal entre log N y M, tal como se muestra en la
Figura (5.1).
La ley de Gutenberg y Richter, es adoptada generalmente para definir la

recurrencia de sismos en una determinada regin y para su correcta estimacin, se
necesita disponer de informacin ssmica para un intervalo de tiempo bastante
grande. Kulhanek (2005) indica que las desviaciones lineales en la distribucin
frecuencia magnitud debe considerar lo siguiente:
-
Para magnitudes pequeas, los datos contenidos en los catlogos son

incompletos; sin embargo, el decremento del parmetro b (bajo el umbral
de magnitud), puede no ser slo un efecto artificial de un catlogo
incompleto, sino a que los sismos de magnitudes pequeas probablemente no
son tan nmeros.
Para magnitudes mayores, las desviaciones se deberan a la saturacin de las

escalas de magnitud. Este problema es asociado a la forma de medir la
48

magnitud de un sismo y en particular cuando se trata de magnitudes como mb

o Md (Stein y Wysession 2003).
-
Otra razn de las desviaciones lineales es el tamao de los catlogos

disponibles que por lo regular son muy cortos en tiempo y no llegan a incluir
la ocurrencia de los eventos mayores, con la consiguiente prdida de
informacin.
Por otro lado, Ziga y Wyss (2001) proponen que el valor absoluto del
parmetro b depende tambin de la escala de magnitud usada, ya que hay una clara
evidencia de que las escalas de magnitud cambian cuando se introduce una nueva
instrumentacin y/o nuevos procedimientos de anlisis, de la misma forma que se
modifican con el tiempo, y esto afecta las distribuciones de frecuencia-magnitud.
Finalmente, otra forma de representar esta relacin es de la siguiente manera:
= (M) = e-M
(5.3)
Donde:
-
N es el nmero promedio de terremotos con magnitud superiores a M

(funcin de densidad de probabilidad de M) en la unidad de volumen y
unidad de tiempo.
y b son constantes regionales.
La constante vara ampliamente de un lugar a otro del territorio, mientras que

se mantiene dentro de los lmites estrechos. La relacin que existe entre las
variables a, b y , son las siguientes:
a = Log
(5.4)
= (b/Log e) = b ln 10
(5.5)
Tanto los parmetros a, b como , , se obtienen de un ajuste de

mnimos cuadrados o por el mtodo de mxima verosimilitud.
49
5.1.1 Los Parmetros a y b

5.1.1.1 El Parmetro a: es una medida del nivel de sismicidad regional o
productividad ssmica (Kossobokov y Keilis-Borok, 2000). Esto est determinado
por la tasa de sismicidad, que para una cierta regin, depende del volumen y la
ventana de tiempo considerada. Matemticamente, expresa el logaritmo del nmero
de sismos con magnitudes mayor a cero; es decir, el total esperado de eventos en la
regin si se extrapola la relacin hasta las magnitudes ms pequeas, el parmetro
a puede ser calculado a partir de la siguiente relacin:
a = LogN log ( 10-bM)
(5.6)
5.1.1.2 El parmetro b: el valor de este parmetro se obtiene a partir del clculo

de la pendiente de la recta resultante del ajuste lineal entre el logaritmo de N y la
magnitud de los sismos (M). El valor del parmetro b se relaciona con las
caractersticas fsicas de cada regin ssmica de forma tal, que valores altos de b
implica el predominio de terremotos de magnitud pequea y por tanto, la regin
presenta poca resistencia a deformarse (Figura 5.2a). En cuanto a valores bajos de
b predominan los terremotos de magnitud mayor, indicando una mayor resistencia
del material a la ruptura (Figura 5.2b). Adems, esta constante vara de regin en
regin (Kossobokov y Keilis-Borok, 2000), y es considerada como una medida del
esfuerzo promedio en la misma (Wiemer y Wyss, 2002). Segn Scholz, (1968), el
valor del parmetro b tiene una fuerte relacin con el esfuerzo en el volumen a
considerar, en razn a que existe una relacin inversa entre el valor b y el nivel de
esfuerzo acumulado. Por ejemplo, para una determinada regin, menor valor de b
implica una menor pendiente en la relacin Gutenberg-Richter, que a su vez estara
relacionado con un nmero mayor de eventos de magnitudes elevados con relacin a
la cantidad de eventos menores, esto implica una mayor concentracin de esfuerzos
en la zona. De la misma forma, valores grandes de b implican una concentracin
de esfuerzos menores (Wiemer y Wyss 1997; Ziga y Wyss, 2001; Wiemer y Wyss,
2002). Por consiguiente, algunos autores como Bullen y Bolt (1987), afirman que
este parmetro puede ser utilizado como predictor de terremotos grandes ya que para
una determinada regin ssmica, el valor del parmetro b varia con el tiempo.
50
a)
b)
Figura 5.2: Representacin de los sismos en funcin de su frecuencia magnitud, a)

Valor alto del parmetro b implica predominio de sismos de magnitud pequea.
b) Valor bajo del parmetro b, predominio de sismos de magnitud grande.
(Udias y Mezcua, 1986).
51
Los parmetros fsicos que perturban los valores del parmetro b son: Los
niveles de esfuerzos y la heterogeneidad del medio. Estudios de laboratorio han
demostrado que el valor de b presenta una relacin inversamente proporcional al
esfuerzo (Scholz, 1968) situacin tambin demostrado para terremotos saciados en
Denver (EEUU), donde la variabilidad de la presin de poros fue conocida (Wyss,
1973). Esto significa que la magnitud (M) de un terremoto que acta a un
determinado volumen, depende de los niveles de esfuerzos necesarios para producir
la ruptura; esto porque la magnitud es inversamente proporcional al valor de b
(Aki, 1965; Utsu, 1965; Hamilton, 1967).
En la Figura 5.3 se muestra un ejemplo del experimento de rotura de rocas en

laboratorio, donde se observa que los valores elevados de tensin estn acompaados
de valores pequeos del parmetro b y los valores bajos de tensin, asociados a
valores elevados del parmetro b.
Figura. 5.3: Aplicacin de la presin de confinamiento, y esfuerzo cortante a una roca tipo
granito para ver la variacin del valor del parmetro b, (Amitrano, 2003)
52
5.2 ASPEREZAS Y VALORES M+IMOS DEL PARMETRO b

Parte de la hiptesis de que los valores anormalmente bajos del parmetro
b ayudan a estimar la presencia de zonas con alta concentracin de esfuerzos que
puedan considerarse como asperezas. El clculo de los valores del parmetro b
para el estudio de asperezas en una zona de subduccin como la del Per, se basa en
la distribucin frecuencia-magnitud de los sismos, y est asociado con el estado de
los esfuerzos presentes en las reas de acoplamiento entre dos placas Nazca y
Sudamericana. Por tanto, en las reas donde se han acumulado esfuerzos mnimos y
mximos (asperezas), presentan diferentes distribuciones para la relacin frecuenciamagnitud de los sismos (Wiemer, 1996), de tal manera que la distribucin espacial
del valor del parmetro b permitira identificar reas con altos esfuerzos (valores
bajos de b) y bajos esfuerzos (valores altos de b).
La Figura 5.4 permite ver el comportamiento de los valores del parmetro

b en el grfico de distribucin de frecuencia-magnitud, para dos regiones
ssmicamente diferentes. La curva de color azul identifica a una regin con baja
sismicidad y su curva presenta una pendiente muy alta, con sismos mayores entre
magnitudes de 6 a 6.5 Ms; por consiguiente, esta regin no est ligada o sujeta a
grandes esfuerzos.
Por el contrario, la curva de color rojo se caracteriza por
presentar una pendiente muy baja y a su vez presenta sismos con magnitudes
mximas del orden de 8 a 8.4 Ms, lo que indica que la zona est ligada a grandes
esfuerzos, y a la presencia de asperezas, de tal manera que en esta regin hay mayor
acumulacin de energa elstica y es ms propensa a generar sismos de gran
magnitud. El estudio en detalle del espacio-tiempo de la variacin del valor del
parmetro b permite inferir la presencia de futuras reas de ruptura; y de otras en
donde indica una alta heterogeneidad en el material o a valores
elevados de
gradiente trmico en el caso de volcanes.
53
Figura 5.4: Distribucin de frecuencia-magnitud de los sismos. Valores bajos del parmetro
b son asociados a reas o regiones con altos esfuerzos (asperezas) y valores altos, a
regiones con bajos esfuerzos. Mc es la magnitud mnima de completitud del catlogo
ssmico.
5.3 VARIACIO+ES DEL PARMETRO b
En general el valor del parmetro b oscila entre 0.5 y 1.5, siendo su valor
mayor ms cercano a la unidad. Los valores mayores a 1, caracterizan a los dorsales
meso ocenicas, y valores menores, a los rifts continentales y regiones de sismos
profundos. (Udias y Mezcua, 1986). El parmetro b puede utilizarse como una
forma de cuantificar el tamao promedio de las rupturas debido a que corresponde
a ley de potencias para las magnitudes y por lo tanto, para el rea de ruptura
promedio (Ziga y Wyss 2001). Utsu seal que el valor del parmetro b
comnmente es cercano a la unidad y vara aproximadamente en un intervalo entre
0.3 y 2.0. A pesar de que no existe acuerdo, sobre si el parmetro b aumenta o
disminuye antes de un terremoto fuerte, para muchos autores, a quienes les ha
llamado la atencin la variacin temporal de este parmetro, es claro que ste tiene
comportamientos particulares para cada regin y que acta como un premonitor de
terremotos (Suyehiro, 1966).
54
En la Tabla 5.1, se presenta los resultados del clculo del parmetro de b

para diferentes regiones tectnicas (zonas de subduccin, sistemas de fallas activas,
volcanes activos). Segn estos resultados los valores del parmetro b para el caso
de estudio de zonas de subduccin varan entre 0.4 a 1.6, los cuales indicaran la
concentracin de los niveles de esfuerzos altos y bajos en las zonas de contacto de
placas tectnicas. Asimismo, los valores de b, para el caso de los sistemas de
fallas activas, oscilan en el rango de 0.5 a 2, y en este caso, los valores mnimos
indicaran la presencia de anomalas estructurales asociados a valores con altos
niveles de esfuerzos (Asperezas). Investigaciones recientes revelan que los eventos
ssmicos con origen en fallas inversas estn asociados con valores bajos de b en
comparacin con las de fallas normales. Para el caso de los volcanes activos, los
valores oscilan en el rango de 0.8 a 2.2, estando los valores altos de b relacionados
a zonas con alto valor del gradiente trmico, y zonas con presencia de material
heterogneo.
Tabla 5.1: Fluctuacin de los valores del parmetro b para diferentes regiones
VOLCA+ES
FALLAS
ZO+AS DE
SUBDUCCI+
tectnicas.
Regin
Zona de Subduccin (Mxico)
Zona de Subduccin (Per)
Amrica Central
Zona de Subduccin (Per)
Zona de Subduccin (Nueva
Zelanda y Alaska)
Placa del Pacfico
falla Bocon (Venezuela)
Falla San Jacinto-Elsionore
Falla San Andrs
Falla Calaveras
Yake-Dake, Montaa Hida
Monte St. Helena (Washington)
Monte Redoubt (Alaska)
Montaa Mammoth (California)
Monte Etna (Italia)
Kilauea (Hawaii)
"b"
0.5-1.5
0.6-1.2
0.6-1.6
0.5-1.1
Referencia
Zuiga y Wyss (2001)
Heras y Tavera (2002)
Monterroso y Kulhanek (2003)
En este estudio (2011)
0.8-1.4
0.4-1.0
0.6-2.0
0.5-1.5
0.5-1.1
0.6-1.1
0.8-1.51.7
1-2.2
1-2.2
1-1.8
1.3-3
0.5-1.5
Wiemer y Benoit (1996)

Katsumata (2006)
Bendito et al., (2000)
Wyss et al., (2000)
Wiemer y Wyss (1997)
Wiemer y Wyss (1997)
Enescu y Ito (2002)
Wiemer y McNuutt (1997)
Wiemer y McNuutt (1997)
Wiemer et al (1998)
Murru et al., (1999)
Wyss y Nagamine (2001)
55
5.4 MTODOS PARA EL CLCULO DE a Y b

Para calcular los parmetros de a y b existen dos mtodos, y son:
- Mnimos cuadrados
- Mxima verosimilitud
5.4.1 Mtodo de Mnimos Cuadrados

El mtodo de ajustes de mnimos cuadrados emplea la suma de los cuadrados.
Para un valor x, la recta de ajuste proporciona un valor diferente de y, el mismo que
ser positivo para algunos puntos y negativo para otros, puesto que esto se dispone
alrededor de la recta. Por lo tanto, la suma de estas diferencias para todos los puntos
es poco significativa; es por ello que se emplea la suma de los cuadrados de las
diferencias con la finalidad de medir las discrepancias entre la recta y los puntos
(Ecuacin 5.7).
(5.7)
Para calcular el valor del parmetro b este mtodo realiza un ajuste de la

distribucin de frecuencia-magnitud de sismos, el mismo que se extiende a partir del
punto de mxima curvatura de esta distribucin, medido por la derivada y que en la
figura se representa por (Mc) hasta la mxima magnitud de los datos ssmicos.
Mediante el mtodo de mnimos cuadrados se calcul la recta de mejor ajuste a
partir de la dispersin de datos que resulta de graficar Log(N) vs. M, cuya pendiente
seria el valor de b y el intercepto con el eje M la constante a o magnitud mxima
esperada (Mmax). Este mtodo ajusta los valores a una recta en funcin de la
densidad de datos que existan en una determinada zona (Press et al, 1986). Para el
borde occidental del Per el ajuste de mnimos cuadrados para la expresin de la
relacin Gutenberg-Richter se muestra en la Figura 5.5.
56
5.4.2 Mtodo de Mximo Verosimilitud

El mtodo de mxima verosimilitud (Figura 5.6), estima el valor del parmetro
b basado en una magnitud media de sismos (Aki, 1965 y Utsu 1965). El mtodo
ajusta la recta al valor medio de los datos sobre la magnitud mnima de
homogeneidad, incluida la mxima magnitud observada y normalizada, a partir de la
siguiente expresin (ecuacin 5.8):
(5.8)
Donde:
-
<M> es el valor medio de todas las magnitudes dentro del rango

seleccionado.
Mmin es la magnitud mnima.
El mtodo de mxima verosimilitud es una tcnica de mayor solidez que el de

mnimos cuadrados, y se acerca ms al valor real de b. Ajusta la recta al valor
medio de los datos sobre la magnitud mnima de homogeneidad incluida la mxima
magnitud observada (Weichert, 1980; Press et al, 1986), lo que a su vez normaliza el
aporte que hacen los sismos de diferentes magnitudes cuando el nmero de sismos
grandes es variable (Sha y Bolt, 1980). Los resultados de los clculos del valor del
parmetro b para el borde occidental del Per, mediante los mtodos de Mxima
Verosimilitud (MMV) y el mtodo de Mnimos Cuadrados (MMC) son muy
similares, siendo estas de: b=0.675+/- 0.02 (MMV) y b=0.687+/- 0.03 (MMC).
El mtodo que presenta mejores lmites de confianza en los clculos de los

parmetros a y b es el de Mxima Verosimilitud (Aki, 1965); por lo tanto, para
este estudio se emplear dicho mtodo para el clculo de los parmetros a y b.
57
Figura 5.5: Clculo del valor del parmetro b por el mtodo de mnimos cuadrados para
el borde occidental del Per.
Figura 5.6: Clculo del valor del parmetro b por el mtodo de mxima verosimilitud
para el borde occidental del Per.
58
5.5 METODOLOGA PARA EL CLCULO DE a Y b
Para poder calcular los parmetros a y b de la relacin frecuenciamagnitud de los sismos, periodos de recurrencia y densidad de probabilidad
siguiendo el procedimiento de Wiemer y Wyss (1997), se debe de considerar lo
siguiente:
1. Los estudios realizados en diferentes regiones tectnicamente activas, han

demostrado que el buen clculo del valor de b depender principalmente de
la seleccin y calidad de los datos ssmicos, nmero de eventos vs cobertura,
clculo de la magnitud de completitud, seleccin del radio de muestreo, y del
buen ajuste de la distribucin frecuencia magnitud de los sismos. Por lo
tanto, se debe contar con un catlogo ssmico para un buen periodo de tiempo
seguidamente realizar un anlisis estadstico y desarrollar lo siguiente:
a. Homogeneizar las magnitudes
b. Eliminar las rplicas
c. Fijar la magnitud de completitud (Mc) del catlogo.
2. El algoritmo que permite calcular los valores de a y b en el programa

ZMAP requiere de una serie de parmetros como las dimensiones de la grilla
(x, y), magnitud mnima de completitud del catlogo, dimensin del radio
para el trazado de crculos, etc. Estos parmetros bsicamente dependern de
la densidad de datos que se cuente y de las caractersticas ssmicas de la
regin de estudio, tales como longitudes de ruptura y reas de ruptura de los
eventos principales, fuentes sismognicas, distribucin espacial de la
sismicidad.
3. El procedimiento a seguir se muestra en la Figura 5.7. Primero se cuenta con

un mapa ssmico homogneo y con una magnitud de completitud definida
(Figura 5.7a), seguidamente se procede a crear un reticulado de celdas
(Figura 5.7b), o mallas bordeando toda la regin de estudio, tal como se
ilustra en la Figura 5.7c.
59
4. A partir de cada nodo de las celdas se generan reas de circunferencias con un

radio de muestreo constante, el cual considera como origen el nodo de cada
celda (Figura 5.7d). El espaciamiento de la malla, y el volumen circular
dependern principalmente del nmero de sismos contenidos en el rea de
estudio, es por eso que los parmetros de entrada deben ser bien definidos, y
as reducir los errores en el clculo del valor de b. Para este estudio se fij
espaciamientos de 0.1 a 0.1 equivalente a 10 Km de largo y ancho. Se
consider volmenes esfricos de 150 Km de radio, en base a la superficie
disponible sobre las reas de ruptura (100-500 Km), que produjeron los
grandes terremotos en el periodo histrico e instrumental asociados a la zona
de friccin de placas (Heras, 2002). La curva de distribucin de la relacin
Gutenberg y Richter es calculada para cada volumen circular y a partir de esta
curva, los valores a y b son calculados con el mtodo de mxima
verosimilitud, descrito anteriormente.
5. Los resultados de los diferentes valores de b, Pr (periodos de retorno), y

densidad de probabilidad local de ocurrencia, son presentados en mapas
diferenciados por escala de colores (Figura 5.7e). Los colores sern
proporcionales al valor obtenido para cada parmetro a caracterizar. As, las
reas de color rojo representan valores anmalos de b y estos a su vez
definen la forma y geometra de las asperezas; mientras que, los colores
azules corresponden a valores altos de b y definen reas con poco o
ninguna posibilidad de generar grandes sismos.
60
Figura 5.7: Ejemplo de la secuencia seguida para el clculo del parmetro b. a) Mapa
ssmico homogneo, b) Creacin de la grilla, c) Delimitacin del rea de estudio, d) Radio
de muestreo, e) Resultados finales cuyos valores se muestran en una escala de colores.
.
61
5.6 PERIODOS DE RECURRE+CIA
Tal como se ha descrito anteriormente, los valores bajos del parmetro b

permiten identificar la existencia de posibles asperezas en una regin potencialmente
ssmica. Para estas asperezas, se puede extrapolar el periodo de recurrencia de
futuros terremotos con magnitudes Ms7.0 en base a los valores de a y b
obtenidos de la distribucin frecuencia-magnitud de los sismos.
En la Figura 5.8 se muestra una ilustracin esquemtica de la distribucin

frecuencia-magnitud de los sismos en forma acumulativa
para
tres modelos
propuestos. El primero por Gutenberg y Richter, (1954) que muestra la distribucin

asumida para el clculo estndar del periodo de recurrencia (Tr) para un valor de b=1
(Lnea de roja) y un intercepto en NMax=1. El segundo modelo
muestra la
distribucin observada para las fallas de tipo strike-slip usando el rea de ruptura
(Wesnousky, 1994; Stirling et al., 1996), con una pendiente de aproximadamente b=1
(Lnea discontinua de color azul) y un intercepto de Nma=1, donde Ma es la

magnitud de la rplica ms grande de un sismo ocurrido en un falla. El ltimo
modelo es el propuesto por Wiemer y Wyss, (1997) en donde consideran la
distribucin frecuencia-magnitud dentro de una aspereza; es decir, dentro de una rea
con un valor bajo de b0.5 (lnea de color negro) y esta se intercepta en Mmax = 1.
De los tres modelos definidos, en este estudio se aplica el tercero a fin de

calcular el periodo de recurrencia de un terremoto de magnitud elevada, esto debido
a que considera los valores de a y b calculados dentro del rea de una aspereza;
por lo tanto, la estimacin del periodo de recurrencia es ms real. Con respecto a los
dos primeros modelos estas consideran un valor promedio estndar de b=1, y no
estara en el rango de los valores que consideran a una aspereza.
62
Figura 5.8: Modelos de frecuencia magnitud de los sismos para la extrapolacin de

terremotos con magnitud elevada (M >7.2) y clculo del periodo recurrencia: (Gutenberg y
Richter, 1954); Wesnousky, 1994; Stirling et al.; 1996 y Wiemer y Wyss, 1997).
Segn el modelo de frecuencia magnitud de los sismos desarrollado por

Wiemer y Wyss (1997), el clculo del periodo de recurrencia de futuros terremotos
para un rango determinado de magnitud, es calculado en base a los valores de las
constantes a y b de la relacin Gutenberg y Richter (ecuacin 5.9), y esta se
expresa como :
(5.9)
(5.10)
Donde:
-
T es el periodo de datos considerado en la base de datos y
M, la magnitud mxima del sismo cuyo periodo de recurrencia se dara

conocer.
Aqu a vara en funcin del espacio, lo que es evidente en cualquier mapa de

sismicidad. Por otro lado b tambin vara fuertemente sobre cortas distancias
63
(aunque en grandes volmenes el valor promedio tiende a mostrarse al valor de 1),

entonces TL(Mmax) tambin es una funcin del espacio.
5.7 PROBABILIDAD DE RECURRE+CIA

La probabilidad local anual por kilmetro cuadrado PL(Mmax) para un sismo de
magnitud elevada, se define como la inversa de la ecuacin del periodo de
recurrencia TL( M).
(5.11)
Donde:
-
Mmax es la magnitud mxima.
A es el rea en Km2, entonces los mapas de Pr y PL son equivalentes.
Otra forma de cuantificar la probabilidad de recurrencia de uno o ms sismos

con magnitud mayores para un periodo de tiempo (T), es empleando la ecuacin
-"T
(5.12), que esta expresado en la forma exponencial e
donde N1 es la relacin
Gutenberg y Richter expresado de la forma: "1 = 10(a1 +bx ) para lo cual, la constante a1
es asumida a una correccin de tiempo tc, valor valor correspondiente al periodo de
tiempo del catalogo ssmico utilizado para el estudio, as a1 ser expresado de la
forma: a1 = a log tc .
Pr (M ) = 1 e "1T
(5.12)
Donde:
-
N1 es el nmero de sismos esperado por ao.
T el periodo de tiempo en aos.
64
5.8 APLICACIO+ES DEL PARMETRO b
La metodologa propuesta por Wiemer y Wyss (1997) ha sido aplicada para el

estudio de la sismicidad a escala local y regional en varias regiones tectnicas ligadas
directamente a la actividad ssmica y volcnica, Esta metodologa ha sido empleado
en numerosos estudios para el anlisis de asperezas y cambios en el ndice de
sismicidad de una determinada regin ssmica (Wiemer y Beniot 1996, Wyss et al.
2000, Monterroso y Kulhnnek, 2003). As mismo, se emple la metodologa para
estudios de evaluacin de catlogos ssmicos (Zuiga y Wyss, 1995), y para analizar
la distribucin frecuencia-magnitud de los sismos en funcin del tiempo y espacio.
Los resultados de la aplicacin del mtodo revelan una variacin

estadsticamente significativa del valor del parmetro b de al menos de 0.4 a 2.0 en
la distribucin de frecuencia-magnitud de los sismos. A travs del programa ZMAP
esos cambios pueden ser calculados en detalle en dos o en tres dimensiones, y los
cambios temporales del valor del parmetro b pueden dar clculos cuantitativos.
Las variaciones ms grandes de los valores de b se encontraron en nueve volcanes
como el Monte Santa Helena, Washington, Monte Redoubt, Alaska, Montaa
Mammoth, California, Monte Etna, Italia, El flanco Sur y Regin Kaoiki cerca al
volcn Kilauea, Hawai. Las anomalas con valores altos de b tambin se
encontraron en la superficie de la zona de subduccin en Alaska, New Zelanda y
Japn a profundidades de 100 150 Km, sugiriendo la posibilidad de ser regiones de
deshidratacin con fusin parcial.
Con respecto a estudios relacionados en fallas, las asperezas se han

caracterizado por presentar valores bajos de b (b<0.6) en volmenes con radios de
5 a 40 Km a largo de la falla San Andrs en la parte Sur y Central de California, en
Japn y Mxico. En zonas con rplicas se encontr valores de b ms bajos
simplemente fuera del rea de ruptura del sismo principal.
65
5.8.1 Estudios en las Zonas de Subduccin

5.8.1.1 Zona de Subduccin de Mxico
A lo largo de la zona de subduccin de Mxico (placa Pacfico y Caribe),
Ziga y Wyss (2001) analizaron la distribucin de los valores del parmetro b con
la finalidad de identificar asperezas y conocer los periodos de retorno para sismos
mayores de 7.2 Ms. Los resultados sugieren la presencia de asperezas con valores de
b<0.8 y periodos de recurrencia del orden de 20 a 60 aos para sismos ocurridos hasta
una profundidad de 60 km. Finalmente, se muestra la presencia de 7 asperezas, de las
cuales las ubicadas en la Costa de Guerrero y Sur de Chiapas representan las ms
importantes debido a su alta probabilidad de recurrencia: valores bajos de b (0.8) y
periodos de retorno de 20 aos (Figura 5.9).
Figura 5.9: a) Mapa de valores del parmetro b a lo largo de la Costa occidental de

Mxico, b) mapa de periodos de recurrencia local estimados probabilsticamente para
sismos con magnitud de 7.2 Ms, usando los valores de b, c) mapa de densidad
probabilstica PL para una magnitud M 7.2. (Zuiga y Wyss, 2001).
66
5.8.1.2 Zona de Subduccin Alaska y +ueva Zelanda
En la zona de subduccin de Alaska y Nueva Zelanda, se observa un

incremento del 45% en los valores del parmetro b para la sismicidad ocurrida a
una profundidad entre 90 y 100 Km, lo cual podra ser debido al incremento de la
deshidratacin de la placa que subduce, y que provocara un aumento de la presin
del poro para reducir los esfuerzos efectivos y junto con ello el incremento del valor
de b. Los valores altos de b, mayores a 1.4, estn relacionados con el volcanismo
presente en estas zonas de subduccin y los valores bajos de b indicaran la presencia
de asperezas (Figura 5.10). En los perfiles presentados (Alaska) se puede ver que los
valores altos de b se encuentran bajo el volcn, coincidiendo con el rea de
deshidratacin de la placa que subduce.
Figura 5.10: Aplicaciones del mtodo de Wiemer y Wyss a las zonas de subduccin de
Alaska y "ueva Zelanda. a) Distribucin de valores del parmetro b en profundidad para
la zona de subduccin de Alaska, b) Variacin de b en la zona de subduccin de Alaska y
"ueva Zelanda.
67
5.8.1.3: Zona de Subduccin de Per
Se realiz el anlisis de la distribucin espacial de la sismicidad asociado al

proceso de interaccin de las placas de Nazca y Sudamericana, con la finalidad de
identificar las reas mas probables para la ocurrencia de futuros terremotos de
magnitud Ms7.2. Los valores del parmetro b fueron calculados utilizando
catlogos ssmicos para el periodo de 1964-1999.
Los resultados obtenidos
corresponden a Heras y Tavera (2002) quienes identificaron la presencia de siete

asperezas de diferentes dimensiones distribuidas a lo largo del borde oeste del Per
(Figura 5.11a). Estas asperezas se ubican: regin Norte; en el lmite de Per y
Ecuador, frente a la costa de Sur de Piura, costa de Ancash y costa Norte de lima. En
la regin Central, en el lmite de Lima e Ica, a lo largo de la costa de Ica y finalmente,
en la regin Sur frente a la costa de Arequipa. Las asperezas presentan un periodo de
retorno de 50 aos aproximadamente (Figura 5.11b).
De estas asperezas, la ubicada frente al departamento de Arequipa dio origen al

terremoto del 2001; es decir, un ao despus de presentado esta investigacin;
asimismo la aspereza ubicada en la costa de Ica dio origen al terremoto de Pisco del
2007, siete aos despus de realizada la investigacin.
68
Figura 5.11: Distribucin de valores anmalos del parmetro b en todo el borde oeste del
Per, donde las reas de color azul indican las asperezas, a) mapa de valores del parmetro
b, b) mapa de periodos de retorno,(Heras y Tavera, 2002).
69
5.8.2 Sistemas de Fallas

La aplicacin del mtodo del valor del parmetro b se realiz en diferentes
sistemas de fallas activas en el mundo. Por ejemplo, en la falla de San Andrs en
California, Wiemer y Wyss (1997) identifican la presencia de importantes asperezas
ubicados en Parkfield (Middle Mountain) y Morga Hill (Figura 5.12a) con valores
anmalos de b=0.46 y b=0.5 respectivamente, y periodos de retorno para la aspereza
de Parkfield (Figura 5.12b) de 23 aos; mientras que, para Morgan Hill de 78 aos,
siendo estos valores similares a los periodos medios de
retorno de terremotos
histricos para esta regin (22 y 72 aos).
Estudios ms detallados fueron realizados por Gerstenber et al. (2001),

quienes sugieren que los valores de b en la zona de California varan en funcin de
la profundidad en las fuentes ssmicas. As, para un rango de profundidad entre 0 a 5
Km el valor de b es de b>1.1 y para 8 a 15 Km, b es del orden b<0.8, lo cual
sugiere la presencia de una aspereza a este nivel de profundidad. En la falla de San
Jacinto-Elsinore (Figura 5.12c) al sur de California, se ha llegado a definir periodos
de retorno para terremotos con magnitud M>6 en base a la historia ssmica. Por
ejemplo, de los 6 terremotos histricos ocurridos entre 1899-1968, 5 de ellos
correlacionan con cuatro asperezas identificadas en dicho estudio. Siendo la ms
importante la aspereza con valores mnimos de b=0.5 presente en el segmento de
Anza (Sur de California).
70
Figura 5.12: Aplicacin del mtodo de Wiemer y Wyss (1997) en sistema de fallas a) Falla
Parkfield, b) Falla Morgan Hill, c) Falla San Jasinto Elsinori.
5.8.3 Aplicacin en Volcanes

Se tena conocimiento, desde hace mucho tiempo, que el comportamiento de
los valores del parmetro b en el interior de los volcanes era elevado, pero fue
sorpresa cuando Wiemer y McNutt (1997) demostraron que este valor en el interior
de los volcanes Mount St. Helens y Mt. Supr, se concentra en pequeos volmenes (r
1.5 Km) con valores normales y bajos de b. Desde entonces, el mtodo fue
aplicado para conocer la variabilidad del valor de b en volcanes como Long Valley
71
y Mammoth Mountain (Wiemer et al., 1998), Montserrat (Power et al,.1998), Etna

(Murru et al.,1999), Katmai (Jolly y McNutt, 1999), Mt, Redoubt (Wiemer et al.,en
preparacin) (Ver Figura 5.13), y adyacente en la zona rift al este de Kilauea (Wyss
et al.,2001).
La conclusin bsica es que los valores del parmetro b en los volcanes no

son tan altos como se pensaba, aunque hay anomalas que se encuentran incrustadas
en el medio. Estas anomalas son altamente significativas y en general, no pueden ser
explicadas con certeza debido a problemas en los catlogos ssmicos. Sin embargo, al
momento de realizar comparaciones de los valores altos de b con otros datos
Geofsicos y Geodsicos, con respecto a la localizacin y extensin de la cmara
magmtica, se lleg a las siguientes conclusiones:
En los alrededores de la cmara magmtica del volcn, siempre se encontr

valores altos de b, particularmente siempre existan anomalas superficiales
con valores altos de b en estos casos no se sabe ciertamente si el magma
est asociado o no con las anomalas de valores de b; a pesar que la
ausencia de una anomala bastantemente fuerte sugiere la ausencia del
magma. Los factores que causen la ausencia de valores altos de b pueden
ser la presencia de fluidos calientes en los sistemas geotermales o
agrietamientos extensos en que un volumen pudo haber obtenido en
erupciones pasadas. Por supuesto que la cmara magmtica y los volmenes
contienen material en fusin parcial, por lo tanto estos no son capaces de
generar terremotos.
En muchos de los volcanes estudiados hasta el momento, las anomalas de los

valores de b no siempre se encuentran localizados directamente bajo el
cono volcnico, (Off-Izu, Etna, Mammoth, Kilauea east zone, Redoubt), lo
cual implica que algunos volcanes presentan conductos de magmas muy
complejos, o que la profundidad de las cmaras magmticas se encuentran en
un rango de 1 a 13 kilmetros de profundidad. Otros volcanes presentan
anomalas en mayor profundidad, lo cual sugiere la heterogeneidad de la
cmara magmtica.
72

Figura 5.13: Imgenes en 3D de la distribucin de los valores del parmetro b debajo y

adyacente a los volcanes. El origen de la actividad ssmica presente en cada volcn se
encuentra en las anomalas de color rojo. a) Monte Santa Helena, Washington, b) Monte
Redoubt, Alaska, c) Montaa Mammoth, California, d) Monte Etna, Italia, e) El flanco Sur y
Regin Kaoiki cerca al volcn Kilauea, Hawai. Wiemer y Mc"utt (1997) (Wiemer et al.,
1998), Montserrat (Power et al., 1998), Etna (Murru et al., 1999), Katmai (Jolly y Mc"utt,
1999), Mt, Redoubt (Wiemer et al., en preparacin).
73
Captulo VI: Anlisis y Seleccin de Datos (Aplicacin al borde occidental del Per)
CAPTULO VI
ALISIS Y SELECCI DE DATOS (APLICACI AL

BORDE OCCIDETAL DEL PER)
Para realizar algn estudio en el campo de la sismologa se debe contar con un
buen catalogo ssmico, ya que los resultados depender de la base de dato utilizado;
por lo tanto, se debe contar en su totalidad con todos los parmetros que caracterizan
a un sismo, tales como fecha, hora, latitud, longitud, magnitud, intensidad y
profundidad, todos calculados en las mismas condiciones a fin de constituir una base
homognea (Tavera y Agero, 2000). Cuanto ms precisos sean los datos, mayor es
la confiabilidad en los resultados obtenidos; caso contrario ellos pueden llevar a
errores para identificar, por ejemplo, patrones de sismicidad y/o estimar el nivel de
peligro ssmico en una determinada regin (Ziga y Wyss, 1995). Es generalmente
aceptado, que el clculo de los parmetros hipocentrales de los sismos contenidos en
un catalogo ssmico mejorar con el tiempo, ya que las tcnicas de registro y la
cobertura de las redes ssmicas son mejores. Adems los catlogos contienen
progresivamente mayor nmero de eventos pequeos y la magnitud de completitud
tiende a disminuir (Ziga y Wyss, 1995; Habermann, 1983). En este captulo se
procede a analizar la base de datos ssmicos proveniente del IGP para realizar el
estudio de la variacin del valor del parmetro b en el borde occidental del Per.
74
6.1 BASE DE DATOS SSMICO

6.1.1 Catalogo ssmico del Per
Para el Per, existe hasta tres Catlogos Ssmicos que contienen los sismos
ocurridos entre los aos 1471 y 1982. Estos catlogos son:
- Catlogo Ssmico Repblica del Per (1471 1982)

Por: Leonidas Ocola.
Proyecto SISAN.
Publicado en 1984.
- Catlogo Ssmico del Per (1500 1984)

Por: A. Espinoza, L. Casaverde, J. Michel, J. Alva, J. Vargas-Newman.
Instituto Geogrfico Nacional de Espaa, US Geological Survey, Universidad
Catlica del Per, Universidad Nacional de Ingeniera.
Publicado en 1985.
- Catlogo Ssmico del Per (1500 1982)

Por: Daniel Huaco.
Instituto Geofsico del Per.
Proyecto SISRA.
Publicado en 1986.
El anlisis de la informacin contenida en estos catlogos, permiti observar

que el nmero y los parmetros hipocentrales de algunos sismos eran diferentes en
cada uno de ellos; por lo tanto, surgi la necesidad de integrar esta informacin en un
catlogo nico tomando como base el publicado por el Instituto Geofsico del Per
(Huaco, 1986). En el ao 2000 el
departamento de Sismologa del Instituto
Geofsico del Per, realiza un trabajo de revisin y actualizacin del Catalogo

Ssmico para el periodo 1471-1982 (Tavera y Agero, 2000).
75
La base de datos utilizada en esta investigacin proviene del catlogo del

Instituto Geofsico del Per (IGP), y su distribucin espacial se muestra en la Figura
6.1. En general, este catalogo tiene las siguientes caractersticas:
Comprende el periodo instrumental y para el clculo de sus parmetros

hipocentrales donde el IGP cuenta con una densa red de estaciones ssmicas.
Mejora en el clculo de los parmetros hipocentrales y en la estimacin de la

magnitud.
Se considera los sismos con focos asociados a la convergencia de placas,

todos con profundidad de hasta 60 Km.
Figura 6.1: Mapa ssmico del Per para el periodo 1970-2010 con profundidad 60Km. Los
crculos de color rojo representan los epicentros para los sismos con Ms7.0 ocurridos
durante este periodo.
76
Una de las caractersticas tpicas de un catlogo, como la del IGP, es que

contienen eventos que no evolucionan segn la tasa uniforme en el espacio-tiempo,
debido a que contienen agrupaciones de sismos denominados clusters que pueden
ser observados en el grfico del Nmero Acumulado vs Magnitud de Sismos (Figura
6.2a). El origen de estos clusters pueden deberse al incremento o disminucin de
estaciones ssmicas, lo que causa el aumento o disminucin del nmero de sismos
registrados y por ende de la tasa de sismicidad. Otro factor que influye es ocurrencia
de rplicas causadas por un evento ssmico de gran magnitud en los alrededores del
hipocentro, llegando enmascarar la sismicidad intrnseca de la zona, y que
incremente o disminuye paulatinamente la sismicidad conforme avanza el tiempo.
Segn la Figura 6.2a la curva presenta una evolucin casi constante desde el ao
1970 hasta 1984, a partir del cual se nota un salto en la curva e incremento de la
sismicidad, en el caso del IGP, que posiblemente se deba al incremento de estaciones
ssmicas. La curva a partir de este ao presenta tres deflexiones muy notorias
originados por las replicas de los terremotos de 1996, 2001 y 2007. Luego la
evolucin de la tasa de sismicidad empieza a presentar un comportamiento contante
y normal hasta el ao 2010, fecha final de la base de datos utilizados en este estudio.
La Figura 6.2b muestra la distribucin temporal de la sismicidad en funcin de

su magnitud. Para el periodo de 1970 a 1984, la data no es homognea y ello se debe
bsicamente a que en este periodo el nmero de estaciones ssmicas en el Per es
mnima sugiriendo una magnitud de completitud del orden de 4.5mb. En general, a
este lapso de tiempo puede ser considerado como el periodo semi-instrumental en el
cual el Per no contaba con una buena cobertura de estaciones ssmicas. A partir de
1985, se produce el incremento de datos debido a que a partir de esa fecha inicia el
aumento de instrumental en el Per. y a la fecha se cuenta con bastante informacin
sobre los parmetros hipocentrales de los sismos, as como la mejora en el clculo
debido a la buena cobertura de estaciones ssmicas y por consiguiente, la magnitud
de completitud disminuye a valores entre 3.0 y 3.2 mb.
77
La Figura 6.2c muestra el histograma del nmero de eventos asociados dentro

del rea de estudio para el periodo 1970-2010. El anlisis del grfico permite
identificar y analizar los cambios en el incremento del nmero acumulado de sismos
y que corresponden a las rplicas de los principales terremotos ocurridos en 1970,
1974, 1996, 2001 y 2007.
Figura 6.2: Anlisis del catalogo ssmico del IGP para el periodo 1970-2010 a) Curva
acumulativa de sismos, los crculos de color rojo indican los cambios en el incremento de la
sismicidad por los terremotos ocurridos, b) Distribucin de sismos en funcin de la
magnitud y tiempo, c) Histograma de sismos.
78
6.2 HOMOGEEIZACI DE LA MAGITUD
Un factor de suma importancia es la homogeneidad del catlogo, los catlogos

pueden ser homogneos en funcin del tiempo, mediante la eliminacin de todos los
eventos que se encuentren por debajo de una magnitud mnima (magnitud mnima de
completitud o de corte, Mc), que es un parmetro importante. Sin embargo este
parmetro cambia con el tiempo en la mayora de los catlogos, por lo que se tiene
que tener especial cuidado al eliminar eventos bajo una magnitud nica
Debido a que las magnitudes estn definidas en funcin de la respuesta de un

sismgrafo o de tipo de ondas que corresponde a un rango de periodos determinados,
se produce el fenmeno de la saturacin de la escala. Cmo a medida que aumenta el
tamao de un terremoto, el mximo de amplitud se va desplazando hacia frecuencias
ms bajas, estas escapan de la deteccin del instrumento y tipo de ondas utilizadas,
produciendo la saturacin de escala que no responde a magnitudes mayores. Puesto
que, los sismos con magnitud elevada presentan superficies de ruptura extensas y por
lo tanto, una mayor irradiacin de energa con periodos ms largos.
En estas
condiciones las magnitudes de ondas de cuerpo (mb) estimaran incorrectamente el

tamao del terremoto, siendo la mxima magnitud posible de cuantificar de 6.0-6.6.
En cambio para las ondas superficiales Ms con periodos mayores, la magnitud
mxima capaz de estimar es alrededor de 8.5-8.7.
Segn los resultados de los estudios realizados sobre la variacin del parmetro
b, se ha mostrado de que el valor de b vara considerablemente de acuerdo al
tipo de escala utilizada para medir la magnitud de los sismos, ya que las diferentes
escalas al utilizar la amplitud del registro de los sismos, resultan ser frecuentemente
afectadas debido a las propiedades fsicas de cada regin que atenan en mayor o
menor grado a las ondas (Zuiga y Wyss, 1995; Zuiga y Wiemer, 1999).
79
En la Figura 6.3 se muestra un ejemplo de la saturacin de las escalas mb a Ms.

La magnitud mb ser constante en el rango de 6.0-6.6 y aunque la amplitud de la
seal y el tamao de sismos aumente (E+26, E+27, E+28) siempre se obtiene el
mismo valor de magnitud, contrariamente para la magnitud Ms, la saturacin de
escala se produce entre 7.0-8.0 Ms.
Figura 6.3: Relaciones tericas entre el espectro de la fuente y la frecuencia para diferentes
rangos de magnitud mb y Ms. (Geller, 1976).
80
A fin de resolver este inconveniente para la base de datos utilizada en este

estudio, se procedi unificar los datos ssmicos a la escala de magnitud Ms. Para tal
objetivo, se ha seleccionado de la base de datos del catalogo IGP un total de 145
sismos con magnitud entre 4.5 y 6.6 expresados en las escalas mb y Ms. Para sismos
de gran magnitud M7.0 se ha tomado los datos directamente del catalogo CMT,
para luego realizar una regresin lineal entre ambas escalas, obteniendo la siguiente
relacin (Figura 6.4). Esta ecuacin permitir convertir los valores de magnitud mb a
Ms para el clculo respectivo de los valore de b.
Ms = 1.744 x mb - 4.1448
(6.1)
Figura 6.4: Regresin lineal para las magnitudes mb y Ms.
81
6.3 MAGITUD MIMA DE COMPLETITUD (Mc)
Mc es aquella magnitud para la cual todos los eventos que ocurren por arriba
de esta son reportados. La determinacin de Mc de un conjunto de datos es
relativamente simple y depende de la distribucin emprica de la magnitud (log N= abM) y es determinada mediante la grfica de nmero acumulativo de eventos contra
magnitud. Estos datos son ajustados a una lnea recta y Mc es el nivel para el cual los
datos caen por debajo de esta lnea.
6.4 ELIMIACI DE RPLICAS
Las rplicas son eventos ssmicos de magnitudes menores y que ocurren

posteriores a un evento principal en una vecindad pequea en tiempo y espacio. Para
su anlisis se ha utilizado la ley de Omori (a la que posteriormente se ha validado
con algunas teoras), la misma que permite relacionar el nmero de rplicas y el
tiempo transcurrido. El resultado muestra el decrecimiento del nmero de rplicas
como el inverso del tiempo transcurrido despus del evento principal. Otra
observacin comn, es que la magnitud de la rplica mayor es generalmente un
grado o ms, menor a la magnitud del evento principal. El tiempo en el que se
pueden observar la ocurrencia de rplicas depende de la magnitud del evento
principal, pero en general va de las pocas horas hasta meses, e incluso aos para
sismos con magnitudes mayores a 7.0 grados. Sin embargo, la rplica mayor ocurre
por lo general unas horas a unos pocos das despus el evento principal.
6.4.1 Relacin Omori-Utsu

Omori (1894a, b), propuso una frmula para representar el decaimiento del
nmero de rplicas con el tiempo, a partir del decremento de las frecuencias de
medio da y mensuales de las rplicas que siguieron el sismo de 1891 de Nobi,
(Japn). Omori mostr que la frecuencia de las rplicas n(t) por intervalo de tiempo
(t) es bien representada por:
(6.2)
82
Donde:
-
n(t) es la frecuencia de rplicas.
K y c son constantes.
t es el tiempo de origen del evento principal.
La explicacin fsica de la frmula relaciona observaciones como el posdeslizamiento sobre la ruptura de una falla durante el evento principal y la emisin
acstica despus de la fractura de una muestra de roca, etc. (Utsu et al., 1995). Utsu
et al. (1995) y Utsu (1961) enfatizan que el decaimiento del nmero de rplicas de
varios eventos ssmicos, es ms rpido que aquel esperado en la frmula general de
Omori. La pendiente de la curva de nmero acumulativo de sismos n(t), graficado
contra log t, tiende a decrecer con el tiempo (Utsu et al., 1995):
(6.3)
Donde:
-
p es la tasa de decaimiento de la secuencia de rplicas.
Utsu et al. (1995) y Utsu (1961) llamaron a la ecuacin
(6.3) frmula
modificada de Omori (actualmente denominada relacin Omori-Utsu). El valor p

describe la tasa de decaimiento de la secuencia de rplicas. En trminos relativos, los
valores altos de p indican un rpido decaimiento de la secuencia de rplicas y un
valor bajo, el decaimiento lento de la secuencia.
6.4.2 Algoritmo de Reasenberg

El algoritmo de Reasenberg (1985), conocido como el mtodo de enjambres,
permite depurar las rplicas asociados a grandes eventos ocurridas en una regin.
Este mtodo asocia dos eventos, analizando las diferencias del tiempo de ocurrencia
en la zona de interaccin de los sismos. El algoritmo se basa en la ley de Omori,
magnitudes de los eventos y los principios del momento de segundo orden.
El momento de segundo orden (tambin llamado funcin de correlacin de

punto doble) es un mtodo relativamente poderoso para la identificacin de patrones
83
de sismicidad coherentes de pares precursores de pares de eventos, y es usado para

describir la interdependencia estadstica de paridad de sismos. Y la unidad elemental
de la interaccin de un sismo.
En el anlisis de un catlogo, el j-simo evento ej est representado por un
vector de 5 elementos:
xj= xjk, k=1,5
(6.4)
Donde:
-
(xj1, xj2) especifica el epicentro.
xj3 la profundidad.
xj4 la magnitud.
xj5 el tiempo de origen.
Entonces, el momento de segundo orden esta expresado como:

m2(x1,x2) = m1(x1) m(x2x1)
(6.5)
Donde:
- m(x2x1) es el momento condicional del proceso.
Por ejemplo, la esperanza del valor de un nmero de eventos con parmetro x2

dado que ha ocurrido otro evento en x1. La distribucin de pares de eventos en el
catlogo de acuerdo al momento de segundo orden es m2(x1,x2) siendo igual a la
esperanza del valor del nmero de pares de sismos (ei, ej). En la prctica, m2(x1,x2)
es calculado por la suma en el catlogo, cuando ei y ej son dibujados desde la misma
poblacin de eventos, m2 es la matriz de autocorrelacin y es simtrica con respecto
a i y j. Cuando ei y ej son dibujados de distintos subconjuntos del catlogo, m2 es
una matriz (no simtrica) de correlacin de cruz.
Reasenberg (1985), modela una secuencia de rplicas como un proceso

Poissoniano dependiente del tiempo, calculando la probabilidad P de observar uno o
ms eventos en un mismo intervalo de tiempo, determinando el intervalo necesario
para esperar el prximo evento en la secuencia, basado en la confianza de P. Para
P=0.95, es:
84
= 3t / 10 2(M 1) / 3
(6.6)
Donde:
- M = Mmx Mmn
- Mmx corresponden a las magnitudes mnimas del catalogo.
- Mmn corresponden a las magnitudes mximas del catlogo.
El nmero de enjambres y el nmero de eventos agrupados varan en un amplio

intervalo de P (correspondiente a la duracin de las zonas de interaccin); sin
embargo, para P < 0.8 (que corresponde a una breve interaccin de las zonas), el
nmero de eventos asociados con un enjambre decrece; mientras que, el nmero de
enjambres incrementa. En la actualidad la depuracin de rplicas de una base de
datos ssmicos se hace utilizando la tcnica de Reasenberg.
Debido a la gran importancia de tener una base de datos uniforme y

homogneo, para cumplir los requerimientos en un estudio probabilstico, se realiz
la depuracin de los eventos ssmicos dependientes, dado que en la mayora de los
eventos dependientes son rplicas, producto de los grandes terremotos ocurridos en
el borde occidental del Per para ello, se utiliz el algoritmo desarrollado por
Reasenberg (1985), basado en la relacin de Omori-Utsu y el momento de segundo
orden. Los parmetros utilizados fueron los propuestos en el trabajo original de
Reasenberg, y que incluyen errores en la localizacin horizontal de 10 km y 33 Km
en la vertical. Los valores de los tiempos mnimos y mximos de das del contenido
de eventos dependientes se consideraron de max=100 das y min=0.5 das y para una
magnitud de completitud de 3.5M.
El anlisis de la data ssmica utilizado, permite constatar que la dependencia de

los resultados en los parmetros no es muy significativa para los datos utilizados en
este estudio, con excepcin de la dependencia en el error de localizacin. Este ltimo
es un valor mnimo, pues se espera que las incertidumbres sean mayores al mismo,
pero al aumentar el error, los resultados no muestran enjambres diferentes ya que el
algoritmo tiende a agrupar eventos pertenecientes a otros enjambres cercanos.
85
En la figura 6.5 se muestra los resultados obtenidos de la aplicacin de este

algoritmo a la base de datos utilizados en este estudio, y se puede observar la
presencia de al menos seis enjambres importantes, dentro de los cuales se
encontraron los eventos principales de los sismos con magnitudes entre 7 y 8Mw
ocurridos en 1970, 1974, 1996, 1996, 2001 y 2007.
Figura 6.5: Mapa de sismicidad en donde se identifican los enjambres identificadas

mediante el algoritmo de Reasenberg (1985), para el borde occidental del Per. Las zonas
de color violeta indican los enjambres y los crculos de color rojo indican los sismos 7.0M.
86
6.5 BASE FIAL DE DATOS SSMICOS
La base de datos ssmicos utilizado en este estudio fue analizada y evaluada

para un intervalo de profundidad menor a 60 Km debido a que en este rango es
donde se producira la mayor friccin de placas y mayor liberacin de energa a lo
largo de todo la superficie de acoplamiento ssmico del proceso de subduccin
(Guardia y Tavera, 2010). Los valores dependientes o rplicas fueron depurados
aplicando el algoritmo de Reasenberg (1985) y como producto final, en las Figuras
(6.6, 6.7 y 6.8) se muestran los resultados obtenidos.
En la Figura 6.6a se muestra el grfico de nmero acumulativo vs magnitud,

calculado a partir de la relacin Gutemberg y Richter (1954), determinndose el
valor de b=0.661+/- 0.02 con una magnitud mnima de completitud de 3.8Ms, para la
base de datos a utilizarse en ente estudio. En la Figura 6.6b se muestra el grfico de
la distribucin espacial de la sismicidad, donde los valores de la magnitud son
considerados desde la magnitud mnima de completitud (Mc) 3.8Ms. Y en este caso
se observa una distribucin homognea de la sismicidad para magnitudes mayores a
Mc.
a)
b)
Figura 6.6: a) Distribucin frecuencia-magnitud b). Distribucin espacial de la sismicidad

para la base de datos del catalogo ssmico del IGP para el periodo1970-2010.
87
En la Figura 6.7a se muestra el grfico de histograma de los eventos ssmicos

seleccionados y utilizados, observe que el mayor nmero de sismos corresponde a los
aos en que ocurrieron grandes sismos (1996, 2001, 2007). El nmero de eventos
independientes acumulados para la regin de estudio (periodo de 1970-2010), es
mostrado en la Figura 6.7b, y en ella se observa que la pendiente es relativamente
constante en el tiempo, lo cual sugiere que puede asumirse que la frecuencia de
eventos es homognea en el tiempo en toda la regin de estudio. Esta distribucin
valida la calidad del catalogo utilizado.
a)
b)
Figura 6.7: a) Histograma de sismos de la base final de datos b) Distribucin acumulada

de sismos en funcin del tiempo para el periodo1970-2010.
En conclusin, la base de datos utilizada en este estudio corresponde al

catlogo ssmico del IGP para el periodo 1970-2010 y considera un total de 1367
sismos con epicentros distribuidos a lo largo del borde occidental del Per. La
magnitud mnima de completitud del catalogo es de 3.8 Ms. (Figura 6.8).
88
Figura 6.8: Distribucin epicentral de los eventos ssmicos que constituyen el catlogo
utilizado en este estudio: periodo1970-2010.
89
Captulo VII: Distribucin del Parmetro b (Aplicacin al borde occidental del Per)
CAPTULO VII
DISTRIBUCI DEL PARMETRO b (APLICACI AL BORDE

OCCIDETAL DEL PER)
Para el clculo de los valores del parmetro b y periodos de retorno local

aplicando la metodologa desarrollada por Wiemer y Wyss (1997) al borde
occidental del Per, se ha utilizado el catlogo ssmico del Instituto Geofsico del
Per para el periodo de 1970-2010 (Tavera y Agero 2010). Este catlogo contiene
1367 eventos ssmicos con magnitud de completitud de 3.8Ms y epicentros ubicados
entre la fosa Per-Chile y la lnea de costa; es decir, en el borde occidental del Per.
En este estudio, la distribucin espacial de los valores del parmetro b en el

borde occidental del Per va a permitir delimitar y ubicar las reas con valores
anmalos (valor bajo de b), y a partir de ellas identificar la presencia de asperezas
que en el futuro podran generar sismos de gran magnitud. Posteriormente, una vez
obtenidos estos valores se proceder a calcular los periodos de retorno y la densidad
de probabilidad de ocurrencia para sismos con magnitudes mayores a 7.0Ms.
Adicionalmente, se estimar las dimensiones de las asperezas, a fin de calcular la
magnitud de los sismos que podran producir, informacin importante para establecer
medidas preventivas que podran ayudar a disminuir el riesgo ssmico en estas
regiones.

90
7.1 DISTRIBUCI DE VALORES DEL PARMETRO b
Para el borde occidental del Per, los valores del parmetro b varan entre
0.5 a 1.1, tal como se muestra en el histograma de la Figura 7.1. La distribucin
espacial de estos valores, ha permitido caracterizar y localizar geogrficamente la
presencia de hasta cinco reas con alta concentracin de energa elstica (reas de
color rojo)
debido a la friccin y acoplamiento de las placas de Nazca y
Sudamericana (Figura 7.2). Estas reas denominadas asperezas, serian fuente de

futuros terremotos a ocurrir y por sus caractersticas, tendran magnitudes Ms
superiores a 7.0.
Figura 7.1: Histograma de los valores del parmetro b observados en cada nodo de la
malla para toda el rea de estudio.

91
De acuerdo a la Figura 7.2 las reas significativamente anmalas con valores
bajos del parmetro b se ubican en:

Regin Sur del Per, cubriendo las costas de los departamentos de Tacna,
Moquegua y parte sur del departamento de Arequipa.
Al norte del departamento de Arequipa, frente a la localidad de Yauca.
Exactamente frente a la costa del departamento de Lima.
Al norte del departamento de Lima y frente al departamento de Ancash.
Frente a las costas de los departamentos de Lambayeque y Piura.
Sobre el mapa de valores del parmetro b se ha ploteado los epicentros de los
terremotos con magnitudes 7.0 (estrellas), ocurridos en los ltimos 40 aos, estando
los terremotos de 1970 y 1996 alrededor de las reas de valores anmalos del
parmetro b, solo el terremoto de 1974 se encuentra dentro de la aspereza ubicado
frente al departamento de Lima, pero a diferencia de los anteriores, el proceso de
ruptura de este terremoto fue asociado a procesos de deformacin por extensin que
se habran producido en el interior de la placa de Nazca (Stauder, 1975).
Los epicentros de los sismos de Nazca, Arequipa y Pisco (1996, 2001, 2007),
se distribuyen dentro de las reas de color azul, es decir, zonas de bajos esfuerzos
(altos valores de b), lo que indica que los sismos ya ocurridos liberaron la energa
en su totalidad y se encuentran en un estado post-ssmico o etapa de acumulacin de
energa.
De acuerdo a estos resultados, las reas con valores anmalas del parmetro
b o asperezas (reas de color rojo), quedan bien determinadas y sugieren
longitudes de ruptura equivalentes a los terremotos histricos ocurridos en 1746 en el
caso de Lima y 1868, y 1913 para la zona sur de Per, las dos primeras con
magnitudes entre 8.0 y 9.0Mw aproximadamente. Para la aspereza de la zona norte
de Per, no es posible asociarlo con algn sismo histrico importante, aunque el
sismo 1619, segn Silgado (1978), habra tenido una longitud de ruptura de 100 a
150 Km, y no se repite desde hace 392 aos.

92
Figura 7.2: Distribucin espacial de los valores del parmetro b en el borde occidental
del Per. Las reas de color rojo indican zonas de acumulacin de energa elstica. Las
estrellas amarillas los epicentros de los principales terremotos ocurridos en el periodo
1970-2010, y las barras indican las longitudes de las reas de ruptura estimada para los
terremotos histricos.

93
7.2 RESULTADOS DEL PERIODO DE RECURRECIA

El mapa de periodos de recurrencia de la Figura 7.3 fue calculado a partir de la
ecuacin 5.10 (Capitulo V) para una magnitud
Mmax=7.0
considerando la
distribucin de los valores del parmetro b (Figura 7.2). El resultado proporciona

un mapa con reas o asperezas definidas con mayor resolucin, siendo las mismas
coherentes con el mapa de valores del parmetro b. Las reas con valores mnimos
de tiempo de recurrencia local oscilan entre 40 y 70 aos y presentan las siguientes
caractersticas.
Aspereza (A1).- Aspereza localizada en la regin sur del Per, frente a las costas de
los departamentos de Tacna, Moquegua y Parte Sur del departamento de Arequipa,
(16.436 a 19.129 S): periodo de retorno de 50 aos.
Aspereza (A2).- Aspereza localizada al norte del departamento de Arequipa, frente a

la localidad de Yauca, (15.564 a 16.147 S): periodo de retorno entre 60 y 70 aos.
Asperezas (A3).- Aspereza ubicada en la regin centro del Per, frente a la costa del
departamento de Lima (12.130 a 13.367 S): periodo de retorno de 50 y 70 aos.
Asperezas (A4).- Aspereza ubicada frente a las costas de los departamentos de

Ancash y parte norte del departamento de Lima (9.417 a 11.914 S): periodo de
retorno de 50 a 70 aos.
Asperezas (A5).- Aspereza ubicada en la parte norte del Per, exactamente frente a
las costas de los departamentos de Piura y Lambayeque (5.397 a 7.132S): periodo
de retorno de 60 aos.

94
Figura 7.3: Mapa de periodos de retorno local para asperezas en el borde occidental del
Per. Las reas de color rojo indican las asperezas y estn enumeradas, de Sur a 3orte,
como A1, A2, A3, A4, A5.

95
7.3 RESULTADOS DE LA PROBABILIDAD DE RECURRECIA

Adicionalmente, se procedi a calcular la probabilidad de ocurrencia a partir de
las ecuaciones (5.11 y 5.12) (Captulo V) para un sismo con magnitud Mmax=7.0, y
los resultados obtenidos se muestran en la Figura 7.4a (Mapa de la distribucin
espacial de densidad de probabilidad anual por km2) y 7.4b (Mapa de probabilidad
recurrencia anual para un periodo de retorno promedio de 50 aos). De acuerdo a la
escala de colores, las reas de color rojo corresponden a las asperezas que presentan
la ms alta probabilidad de producir un sismo de magnitud Ms7.0. Y ellas se
encuentran ubicadas en:
Frente a las costas de los departamentos de Tacna, Moquegua y parte sur del
departamento de Arequipa (Aspereza 1).
Ubicada frente a los departamentos de Lima y Ancash que en conjunto forman

una gran aspereza (Aspereza 3 y 4).
Finalmente, la aspereza A5 ubicada en la zona norte del Per, presenta una baja
probabilidad de ocurrencia; por lo tanto, la regin no presenta mayor peligro a la
ocurrencia de sismos de grandes magnitudes en los prximos 50 aos.

96
Figura 7.4: Mapa de Probabilidades de recurrencia para el borde occidental del Per. Las reas de color rojo indican zonas ms
probables a la ocurrencia de sismos de gran magnitud (Ms 7.0Ms). a) Mapa de densidad de probabilidad local por unidad de
rea. b) Mapa de probabilidad de recurrencia anual para un periodo de retorno promedio de 50 aos.

97
7.4 ESTIMACI DE LA MAGITUD A PARTIR DEL REA DE

ASPEREZAS
La distribucin espacial de los valores del parmetro b y de periodos de
retorno, permiti determinar la ubicacin geogrfica y la geomtrica de las asperezas
(reas con alta concentracin de esfuerzos), presentes en el borde occidental del
Per. De manera complementaria, es posible determinar la magnitud de los sismos
que se produciran. A partir de las dimensiones del rea de las estas asperezas,
asumiendo que ellas en conjunto forman parte de toda el rea involucrada en la
ruptura. La magnitud momento Mw puede ser estimada utilizando la ecuacin de
Kanamori y Anderson (1975).
(7.1)
(7.2)
En la Figura 7.5 se muestra un ejemplo del clculo del momento ssmico Mo, a
partir del rea de rplicas correspondiente al rea de ruptura asociada al sismo del 23
de junio del 2001. El rea de ruptura fue de 370x150 Km2 y segn la ecuacin (7.1),
se obtuvo un momento ssmico de 1.1x10^28, equivalente a una magnitud momento
de Mw=7.9 (Tavera et al, 2001).
Figura 7.5: Distribucin espacial de las rplicas del sismo del 2001, ocurridos entre el 25
de junio y el 19 de Julio. Las estrellas representan a los epicentros del sismo principal y las
3 replicas de mayor magnitud. Las lneas discontinuas indican el rea de replicas (Tavera et
al, 2001).
98
En este estudio se ha procedido a ajustar el rea de las asperezas a un

rectngulo a fin de facilitar los clculos (Figura 7.6), obteniendo los siguientes
resultados.
Aspereza (A1).- El rea estimada para esta aspereza es de 89797 Km2 y la
magnitud que podra alcanzar el sismo sera del orden de 8.8Mw.
Aspereza (A2).- Esta aspereza presenta un rea aproximada de 4912 Km2
equivalente a una magnitud estimada de 7.5Mw.
Aspereza (A3).- Esta aspereza presenta un rea de 23943 Km2 y su magnitud
estimada es de 8.2Mw.
Tabla 7.1. Resultados del clculo de la magnitud estimada a partir del rea de la
asperezas en todo el borde Occidental del Per.
N
Asperezas
rea
(Km2)
Mo
(dina-cm)
Mw (estimada)
89797
1,88361E+29
8.8
4912
2,40983E+27
7.5
23943
2,59338E+28
8.2
46650
7,05303E+28
8.5
8207
5,20445E+27
7.7

99
Figura 7.6: Estimacin de la magnitud en funcin del rea de las asperezas identificadas en
este estudio y valores obtenidos utilizando la ecuacin de Kanamori y Anderson (1975).

100
Captulo VIII: Discusin e Interpretacin
CAPTULO VIII
DISCUSI E ITERPRETACI
En la actualidad, los estudios sobre la prediccin de terremotos sigue siendo
una gran controversia; sin embargo es conocido que para predecir un terremoto se
debe disponer de la siguiente informacin
Prediccin = Lugar + Tamao + Fecha
De estos parmetros, solo los dos primeros han sido posibles de identificar y
cuantificar, restando aun poder definir la fecha de ocurrencia de los sismos. En tal
sentido, es factible afirmar que la ciencia solo ha llegado a pronosticar la ocurrencia
de los sismos y para ello, se han aplicado diversas metodologas puramente
estadsticas haciendo uso de catlogos ssmicos y recientemente, cuantitativas con el
uso de informacin ssmica, GPS e interferometra.
Para el Per, este tipo de estudios y sobre todo los que consideran el
seguimiento de los grandes sismos y la ubicacin geomtrica de sus reas de ruptura,
se iniciaron en los aos 72 con los primeros mapas de Lagunas Ssmicas
propuestos por Kelleher (1972) para todo el borde de la placa del Pacfico. Desde
estos aos, se identificaron como futuras reas a ser afectadas por sismos las que se
encuentran en la regin sur del Per y norte de Chile, as como, la que se encontraba
frente al departamento de Ica. Estudios posteriores fueron realizados por Heras y
Tavera (2000), (Publicado en el 2002), haciendo uso de la metodologa de Wiemer y
101
Wyss (1997) para llegar a identificar la presencia de asperezas que daran origen a
futuros terremotos en la costa de Tacna, Moquegua, Arequipa (7), Lima (3,4) e Ica
(5) (Figura 8.1a). Posteriormente, revisiones ms detalladas sobre la distribucin
espacial de las reas de ruptura de grandes sismos, permitieron a Tavera y Bernal
(2005) definir con mayor precisin la presencia de reas de futuros terremotos en la
costa de Arequipa, Moquegua y Tacna, y otra en la costa norte del departamento de
Ica (Figura 8.1b).
En los aos 2001 y 2007 ocurrieron los terremotos de Arequipa y Pisco (Figura
8.1a, b), con magnitudes del orden de 8.0Mw involucrando en su proceso de ruptura,
las reas previamente definidas en los estudios ya indicados.

102
ATES
a.
b.
DESPUES
Figura 8.1: Lagunas ssmicas antes y despus de los terremotos del 2001 y 2007. a) Distribucin de las reas de ruptura y asperezas
asociadas al proceso de subduccin definidas por Heras y Tavera (2002). Y b) Distribucin de lagunas ssmicas segn Tavera y Bernal
(2005).
103
En este estudio, se continua con el trabajo iniciado por Heras y Tavera (2000)
llegndose a re-evaluar la presencia de las asperezas ya identificadas y proponer
otras que podran presentarse despus de ocurrido los sismos del 2001 y 2007. Los
resultados obtenidos en este estudio, permiten definir con mayor precisin la
presencia de asperezas ubicadas exactamente frente a los departamentos de
Arequipa-Moquegua-Tacna, Lima- Ancash y Lambayeque. De acuerdo a las
dimensiones de las reas de las asperezas, se estima el tamao de los terremotos que
podra ocurrir y afectar a estas ciudades y departamentos.
Por otro lado, estudios realizados con mediciones de GPS en las regiones norte
y centro de Per (Figura 8.2a), han permitido a Villegas (2009), identificar la
presencia de dos reas de mximo acoplamiento ssmico frente a la costa de Lima,
estas reas son similares, en cuanto a su ubicacin a las propuestas en este estudio
(Figura 8.2b).

104
a.
b.
Figura 8.2: Comparacin de resultados obtenidos con datos de GPS y anlisis del valor del parmetro b a) Zonas de mayor
acoplamiento de placas segn Villegas (2009), y b) Asperezas identificados en este estudio. Las reas de color rojo indican la
ubicacin de las asperezas.
105
Del mismo modo, estudios realizados en la regin centro y sur de Per con
datos de GPS, ha permitido a Perfettini, (2010), identificar nuevas reas de
acoplamiento ssmico frente a las costa de los departamentos de Lima, Arequipa,
Moquegua, Tacna, Arica e Iquique (Figura 8.3), as como frente al departamento de
Lima, siendo la ubicacin de las reas propuestas similares a las identificadas en este
estudio (Figura 8.3).
Perfettini, (2010)
Este estudio
Figura 8.3: Comparacin entre los resultados obtenidos en este estudio (valores de b) y
estudios de GPS, (Perfettini, 2010).

106
En conclusin, es notorio que las reas anmalas para valores del parmetro
b asociados a asperezas, en cuanto a ubicacin y geometra, son similares a las
propuestas por otros autores utilizando metodologas estadsticas y tcnicas de GPS.
En todos los casos, las magnitudes estimadas a partir del rea de la aspereza o de
acoplamiento mximo, permiten considerar que los sismos a ocurrir en el futuro
tendran valores mayores a 8.0Mw. En el caso del sur del Per, el sismo a ocurrir
tendra una magnitud del orden de 8.8Mw y afectara a la regin norte de Chile, y en
el caso de la regin central, sera de 8.5Mw. Para la aspereza pequea identificada
frente a la costa de la localidad de Yauca, el sismo tendra una magnitud cercana a
7.5Mw. De estas asperezas, las dos primeras tienen mayor probabilidad de
ocurrencia (Tabla 8.1).
Tabla 8.1: Asperezas con alta probabilidad de generar un sismo de gran magnitud en el
borde occidental del Per.

Asperezas
Localizacin
Coordenadas
Latitud Sur ()
Arequipa,
Moquegua y Tacna
16.436 - 19.129
Costas de Lima y
Ancash
9.417 - 13.367

107
Conclusiones
COCLUSIOES
Los resultados obtenidos en el presente estudio han permitido llegar a las
siguientes conclusiones:
1. Los valores del parmetro b a lo largo del borde occidental del Per varan
entre 0.5 a 1.1 y los periodos de recurrencia para terremotos con magnitud mayor
a 7.0Ms, oscilan entre 50 y 70 aos. Estos valores sugieren que la regin de
estudio es, por dems, compleja y heterognea en su proceso de deformacin y
acumulacin de energa.
2. La distribucin espacial de los valores del parmetro
b, han permitido
identificar la ubicacin geogrfica de 5 asperezas que representan ser fuentes

probables de producir sismos con magnitudes mayores a 7.0Ms en el futuro.
Estas asperezas son:
Aspereza ubicada en la regin sur del Per, frente a las costas de los
departamentos de Tacna, Moquegua y parte Sur del departamento de
Arequipa, (16.436 a 19.129 S).
Aspereza ubicada al norte del departamento de Arequipa, frente a la
localidad de Yauca, (15.564 a 16.147 S).
Aspereza ubicada en la regin centro del Per, frente a la costa del
departamento de Lima (12.130 a 13.367 S).
Aspereza ubicada frente a las costas de los departamentos de Ancash y
parte norte del departamento de Lima (9.417 a 11.914 S).

108
Conclusiones
Aspereza ubicada en la parte norte del Per, exactamente frente a las

costas de los departamentos de Piura y Lambayeque (5.397 a 7.132S).
3. Las probables reas de ruptura asociadas a los grandes sismos ocurridos en el

periodo histrico (1619, 1746, 1868, 1913), coinciden con la distribucin de las
reas de valores anmalos del parmetro
b; es decir, con las asperezas
identificadas en este estudio.
4. Las asperezas identificadas presentan variados periodos de retorno y magnitudes,

segn se detalla a continuacin.
Localizacin
Ubicacin
Latitud Sur ()
Arequipa, Moquegua y
Tacna
16.436 19.129
50
Norte de Arequipa
15.564 16.147
70 - 80
Costas de Lima
12.130 13.367
Ancash y Norte de Lima

Piura y Lambayeque
(Mw)
Estimado
Sismos
Histricos
8.8
1868
49
7.5
1913
70 - 80
75
8.2
1746
9.417 11.914
50 - 80
75
8.5
1746
5.397 7.132
60
25
7.7
1619
Pr (M7.0)
(aos)
PL (M7.0)
(%)
75
Pr: Periodo de Retorno; PL: Probabilidad de Recurrencia; Mw: Magnitud Momento.
5. Las asperezas identificadas en este estudio presentan geometras y ubicaciones

geogrficas similares a los propuestos por otros autores basados en datos de GPS.
6. Las asperezas con mayor probabilidad a generar un sismo de gran magnitud se

encuentran
frente a los departamentos de Lima, Ancash y frente a los
departamentos Arequipa, Moquegua y Tacna.

109
Recomendaciones
RECOMEDACIOES
De acuerdo a los resultados obtenidos, es importante que las autoridades de las
regiones o departamentos involucrados en este estudio, tomen las medidas necesarias
de prevencin y mitigacin.
reas probables de ruptura ssmica en el borde occidental del Per a partir de la

110
Bibliografa
BIBLIOGRAFA
Amitrano, D. (2003). Brittle-ductile transition and associated seismicity:
Experimental and numerical studies and relationship with the b-value. J.
Geophys. Res., 108(B1), 2044, doi: 10.1029/2001JB000680.
Aki, K. (1965). Maximum likelihood estimate of b in the formula log . = a bM
and this confidence limits. Bull. of the Earthquake Research Institute. Vol. 43
pp.237- 239.
Barriga, P. Victor, M. (1951). Los Terremotos en Arequipa, 1582-1868: Documentos
de los Archivos de Arequipa y de Sevilla. Arequipa: La Colmena, S.A., 1951.
Barriga, V. (1939). Libro de los cabildos de la ciudad de Arequipa. Sesin del 10 de
Agosto de 1555, reproducida en documentos para la Historia de Arequipa.
Barazangi, M., Isacks, B. (1976). Spatial distribution of earthquakes and subduction

of the Nazca plate beneath South America. Geology, 4: 686-692.
Bernal, I. (2002). Aproximacin de un modelo detallado de la sismicidad en el
Per: Caractersticas y evaluacin de la energa ssmica liberada. Tesis de
Ingeniera U...S.A. 170 pg.
Bernal, I. (2000). Caractersticas de la sismicidad en la regin sur de Per. Revista
de Trabajos de Investigacin, C.DG-Biblioteca, Instituto Geofsico del Per,
p. 69-80.
Bernal, I., Tavera H. (2002). Geodinamica, sismicidad y energa ssmica en Per,
Monografia CDNG-Sismologia, Instituto Geofisico del Per, Lima, Per, 70
pp.
Bevis, M., and Isacks, B.L., (1984). Hypocentral trend surface analysis : probing the
geometry of Benioff zones, J. Geophys. Res., 89, 6153-6170.
Bullen, K. A., Bolt, B. A., (1987). An Introduction to the Theory of Seismology.
Cambridge, 4 ed., pp. 371395.
Cahill, T. y Isacks, B. L. (1992). Sismicidad y la forma de la subduccin de la placa
de Nazca. Journal of Geophysical Research, vol. 97, no. B12, p. 17 ,503-17,
529.
Dass, S., Aki K., (1977). Fault planes with Barriers a versatile earthquakes model. J.
Geophys. Res. 82.
Dalmayrac, B., Molnar, P. (1981). Parallel thrust and normal faulting in Peru and
constraints the state on stress: earth and Plan. Sc. Let, 55, p.473-481.
DeMets, C., Gordon, R., Aarhus, A., Stein, S. (1990). Current plate motions.
Geophys. J. Int., 101, 425-478.
reas probables de ruptura por sismos en el borde occidental de Per, a partir de la
variacin del valor de b.
111
Bibliografa
Dorbath, L., Cisternas, A. and Dorbath, C. (1990). Assessment of the size of

large historical earthquakes in Peru. Bull. Seismol. Soc. Am., 88, 551-576.
Ericksen, G. E., Plafker, G., Fernandez, J. (1970). Preliminary report on the
geologic events associated with the May 31, 1970, Peru earthquake. U. S.
Geological Survey, Washington, DC, Geological Survey Circular, 639, 25 pp.
Gutscher, M. A., Olivet, J., Aslanian, D., Maury, R., y Eissen, J. P. (1999). The
"lost Inca Plateau": Cause of the flat subduction beneath Peru?, Earth
Planetary Sciences Letters, 171: 335-341.
Geller, R. J. (1976). Scaling relations for earthquake source parameters and
magnitudes. Bull. Seismol. Soc. Am., 66.
Guardia P., Tavera H. (2010). Inferencias sobre la superficie de contacto entre
placas en Per: Nazca y Sudamericana. XV Congreso peruano de Geologa
Cusco-2010.
Gutenberg, R., Richter, C.F. (1944). Frequency of earthquakes in California. Bull.
Seism. Soc. Am. 34, 185-188.
Gutenberg, B., Richter, C. F., (1942). Earthquake magnitude, intensity, energy
and acceleration. Bulletin of the Seismological Society of America. Vol. 32,
.o. 3, pp 163-191.
Gutenberg, B., Richter, C. F., (1954). Seismicity of the Earth (and associated
phenomena). Princeton University Press, Princeton, NJ, EUA.
Gerstenber, M., Wiemer, S., Giardini, D. (2001). A systematic test of the

hypothesis that the b value varies with depth in California. Geoph. Res. Letts.
, 28, 57-60.
Gerstenberger, M.C., Wiemer, S., Giardini, D., Hauksson, E., Jones, L.M.
(2001b). Timedependent hazard assessment for California in near real-time.
Seism. Res. Letts., 72, 273.
Grange, F., Hatzfeld, D., Cunningham, P., Molner, P., Roecker, S.W., Suarez,
G., Rodriguez, A.B., and Ocola, L., (1984). Tectonic implications of the
microearthquake seismicity and f a u l t plane solutions in southern Peru, J.
Geophys. Res., 89, 6139-6152.
Hasegawa, A., Sacks, I.S. (1981). Subduction of the Nazca plate beneath Peru as
determined from seismic observations, J. Geophys. Res., 86, 4971-4980.
Hampel, A. (2002). The migration history of the Nazca Ridge along the Peruvian active
margin: A re-evaluation, Earth and Planetary Science Letters, 203: 665679.

112
Bibliografa
Hamilton, R. M. (1967). Mean magnitude of an earthquake sequence. Bull. Seism.

Soc. Am. 57, 1115-1116.
Habermann, R. E. (1983). Teleseismic Detection in the Aleutian Island Arc.
Journal of Geophysical Research, Vol. 88, .o. B6, pages 5056-5064.
Heras, H. (2002). Anlisis de la distribucin del valor de b en la zona de subduccin
de Per. Tesis de grado U...S.A 41-106 Pgs.
Heras, H., Tavera, H. (2002). Localizacin de reas probables a ser afectadas por
grandes sismos en el borde oeste de Per: Estimacin a partir de periodos de
retorno local basado en distribucin de valores de b. Centro .acional de
Datos Geofsicos, boletn de la S.G.P. p. 7-16.
Huaco, D. (1986). Catalogo Ssmico de Per (1500-1982). Instituto Geofsico del
Per (Proyecto SISRA).
Ishimoto, M., Iida, K. (1939). Observation of earthquakes registered with the micro
seismograph constructed recently. Bull. Earthquake Res. Inst. Univ Tokyo,
443-478.
James, D. (1978). Subduction of the Nazca plate beneath central Per Geology 6, 174178.
Jolly, A. D., Mc<utt, S. R. (1999). Seismicity at the volcanoes of Katmai National
Park, Alaska; July 1995- December 1997. Journal of Volcanologist and
Geothermal Research, 93, 173-190.
Kanamori, H., Anderson, D. L. (1975). Theoretical basis of some empirical
relations in seismology. Bull. Seismology. Soc. Amr., 65, 1073-1096.
Kanamori, H. (1981). The nature of seismicity patters before large earthquake In:
Earthquake prediction International Review. Am. Gephys. Union.
Katsumata, K. (2006). Imaging the high b-value anomalies within the sub ducting
plate in the Hokkaido corner. Earth Planets Space, 58,-e49-e52.
Kendrick, E., Bevis R., Smalley, B., Brooks, B., Vargas, E., Lauria, L., Souto, F.
(2003).The Nazca-South America Euler vector and its rate of change. Journal
of South American Earth Sciences Volume 16, Issue 2, Pages 125-131.
doi:10.1016/S0895- 9811(03)00028-2.
Kelleher, J. (1972). Rupture Zones of large South American earthquakes and some
predictions. J. Geophysics. Res., 77, 2087-2103.
Kulhanek, O. (2005). Seminar on b-value. Department of Geophysics, Charles
University, Prague.

113
Bibliografa
Kossobokov, V. G., Keilis-Borok, V. I. (2000). Implications of a statistical physics

approach for earthquake hazard assessment and forecasting. Pure Appl.
Geophys., 157: 2323-2349.
Mgard, F. (1978). tude Gologique des Andes du Prou Central: Mmoires
ORSTOM, Paris, v. 86, 310 p. Mgard, F., R. Marocco, JC Vicente, C.
Muoz, R.
Marocco, R. (1980). Geologie des Andes Peruviennes, un segment E-W de la chane
des Andes peruviennes : La dflexin dAbancay. Etude gologique de la
Cordillera orientale et des Hauts Plateaux entre Cuzco et San Miguel Sed du
Prou. Doc. O.R.S.T.O.M. 94 pp.
Monterroso, D. A., Kulhnnek, O. (2003). Spatial Variations of b-values in the
Subduction zone of Central America. Geophy. Int. Universidad .acional
Autonoma de Mexico. pp. 575-587.
Murru, M., Montuori, C., Wyss, M. Privitera, E. (1999). The location of magma
chambers at Mt. Etna, Italy, mapped by b-values. Geoph. Res. Letts. , 26,
2553-2556.
<orabuena, E., Timothy, H., Dixon., Seth Stein., Christopher, G., Harrison, A.
(1999). Decelerating Nazca-South America and Nazca-Pacific Plate Motions.
Geophysical Research Letters, Vol. 26 pp. 3405-3408.
Ocola, L. (1984). Catlogo ssmico del Per. Proyecto SISA., Instituto Geofsico del
Per, Lima, pp. 464.
Omori, F. (1894a). On aftershocks. Rep. Imp. Earthq. Inv. Com., 2, 103-138.
Omori, F. (1894b). "On the aftershocks of earthquakes". Journal of the College of
Science, Imperial University of Tokyo, vol. 7, pages 111200.
Perfettini, U. (2010). Desplazamiento Ssmico y Assmico en la zona del terremoto
de Pisco: Avances en la Comprensin del Ciclo Ssmico. XV Congreso
peruano de Geologa Cusco-2010..
Press W. H., Flannery B. P., Teukolsky, S. A., Vetterling W. T. (1986.). The
standard error of the magnitude-frequency b value. Bull. Seism. Soc. Am. Vol.
72, p. 1677-1687.
Power, J.A., Wyss M., Latchman, J. L. (1998). Spatial variations in frequencymagnitude distribution of earthquakes at Soufriere Hills volcano, Montserrat,
West Indies, Geoph. Res. Letts. , 25, 3653-3656.
Polo, T. (1904). Sinopsis de los temblores y volcanes de Per. Boletin de la Sociedad
Geogrfica de Lima. Vol. VII-XVI, 323pags.

114
Bibliografa
Reasenberg, P. (1985). Second-Order Moment of Central California Seismicity,

1969-1982. Journal Geophysical Research, vol. 90, .o. B7, Pages 54795495.
Reid, H. F. (1911). The elastic rebound theory of earthquake, Bulletin of the
Department of Geology, 6, 19.
Rodriguez, L., Tavera, H. (1991). Determinacin con alta resolucin de la
geometra de la zona Wadati-Benioff en el Per Central. Revista Brasileira de
geofsico. 9, 141-159.
Rosenbaum, G., Giles, D., Saxon, M., Betts, G., Weinberg, R., Duboz, C. (2005).
Subduction of the Nazca Ridge and the Inca Plateau: Insights into the
formation of ore deposits in Peru. Earth and Planetary Science Letters 239
(2005) 18 32.
Silgado, E. (1978). Historia de los sismos ms notables ocurridos en el Per, 15131974. Instituto Geolgico Minero del Per, Lima, 131 pp.
Sha, Y., y Bolt, B.A. (1980). The standard error of the magnitude-frequency b value.
Bull. Seism. Soc. Am.Vol. 72, p. 1677-1687.
Sebrier M., Mercier J., Mergad F., Labaucher G. (1985). Quaternary fault and
reverse faulting and the state or stress in the central Andes of Peru. Tectonics,
4, 739-780.
Soler, P., Sebrier M. (1990). At present the crustal thickness reaches 55-60 km
below the highest part of the Western Cordillera decreasing both eastwards.
Stauder, W. (1975). Subduction of the Nazca plate under Peru as evidenced by focal
mechanism and by seismicity. Journal of Geophysical Research, 80: 10531064.
Shcherbakov, R., Turcotte, D., Rundle, L., (2004). A generalized Omoris law for
earthquake aftershock decay. Geophysical Research Letters, Vol. 31,
L11613.
Scholz, C. H. (1968). The Frequency-Magnitude Relation of Micro fracturing in
Rock and its Relation to Earthquakes. Bulletin of Seismological Society of
America, Vol. 58, .o. 1, pp 399-415.
Schneider, J. y Sacks, S. (1987). El estrs en la placa de Nazca debajo de retorcido
del sur del Per a causa de terremotos locales, Geophys J.. Res. 92, 1388713902.
Suyehiro, S. (1966). Difference between aftershocks and foreshocks in the
relationship of magnitude to frequency of occurrence for the great Chilean
earthquake of 1960. Bull. Seismol. Soc. Am., 56, 185-200.
115
Bibliografa
Stirling, M. W., Wesnousky S. G., Shimazaki K. (1996). Fault trace complexity,

cumulative slip, and the shape of the magnitude-frequency distribution for
strike slip faults. A global survey, Geophy. J. Int., 124, 833-868.
Stein, S., Wysession, M. (2003). An introduction to Seismology, Earthquakes and
Earth Structure, Blackwell Publishing, ISBN: 9780865420786.
Tavera, H., Agero, C. (2000). Sismos Sensibles Ocurridos en Per. Boletn del
Centro .acional de Datos Geofsicos-IGP, Lima, 50 p.
Tavera, H., Agero C. (2010). Catalogo ssmico del Per 1961-2010. Boletn del
Centro .acional de Datos Geofsicos-IGP, Lima.
Tavera, H., Salas, H., Jimenez, C., Antayhua, Y., Vilcapoma, L., Millones, J.,
Bernal, I., Zamudio, Y., Carpio, J., Agero, C., Perez-Pacheco, I. and
Rodriguez, S. (2001). El terremoto de Arequipa del Junio de 2001.
Instituto Geofsico del Per, 16 p.
Tavera, H. (1992). El Proceso de Ruptura Ssmica: Barrera o Aspereza?. Boletn
Sociedad Geolgica del Per. Vol. 83, p 69_73.
Tavera, H., Bernal, I. (2005). Distribucin Espacial de reas de Ruptura y lagunas
ssmicas en el Borde Oeste del Per. Volumen Especial . 6 Alberto Giesecke
Matto p. 89-102.
Tavera, H., Bernal, I. (2009). Mapa de sismicidad del Per para el periodo 1964-2008
(Mw4). Instituto Geofsico del Per.
Tavera, H. (2003). La Prediccin Tectnica. Centro .acional de Datos geofsicos.
(C.DG) p. 23 - 25.
Tavera, H., Buforn, E. (1998). Sismicidad y Sismotectnica del Per. Fsica de la
Tierra UCM, . 10, 187-219.
Tavera, H., y Buforn, E. (2001). Fuente mecanismo de los terremotos en Per.
Diario de Sismologia, 5, 519-540.
Udias, A., Mezcua, J. (1986). Fundamentos de Sismologa, San Salvador. El
Salvador,127 _ 128 Pgs.
Utsu, T. (1961). A statically study on the occurrence of aftershocks. Geophysics
Mag., 30, 521-605.
Utsu, T., Ogata Y., Matsura, R. S. (1995). The Centenary of the Omori Formula
for Decay Law of Aftershock Activity. Journal Physics of The Earth, vol.
43, .o. 1, p. 1-33.

116
Bibliografa
Utsu, T. (1965). A method for determining the value of b in a formula log . = a-bM
showing the magnitude frequency for earthquakes. Geophy Bull. Hokkaido
Univ., 99-103.
Villegas, J. C. (2009). Processing and Analysis of the LISN permanent GPS
network, Masters theses Universit de .ice - Sophia Antipolis-France.
Wesnousky, S. G. (1994). The Gutenberg-Richter or characteristic earthquake
distribution, which is it?, Bull. Seismol. Soc. Am., 84.
Wiemer, S., Wyss, M. (1994). Seismic quiescence before the 1993 M=7.5 Landers
and M=6.5 Big Bear (California) earthquake. Bull. Seism. Soc. Am. 84, 3,
900-916.
Wiemer, S., R. Zuniga, Z. (1994). MAP-a softwarep ackaget o analysiss eismicity
(abstract) E, OS, Trans.,A GU, 75, 456.
Wiemer, S. and Benoit, J. (1996). Mapping the b-value anomaly at 100 km depth in
the Alaska and New Zealand subduction zones. Geophys. Res. Let., 23: 15571560.
Wiemer, S., Mc<utt, S. (1997).Variations in frequency-magnitude distribution with
depth in two volcanic areas: Mount St. Helens, Washington, and Mt. Spurr,
Alaska. Geoph. Res. Letts. , 24, 189-192.
Wiemer, S., Katsumata, K. (1999). Spatial variability of seismicity parameters in
aftershocks zones. Journal Geophysical Research, Vol. 104, .o. B6, pages
13, 135- 13, 151.
Wiemer, S. (1996). Analysis of seismicity: New techniques and case studies.
Dissertation thesis, 151 pp, University of Alaska, Fairbanks, Alaska.
Wiemer, S., Wyss, M. (2002). Mapping spatial variability of the frequencymagnitude distribution of earthquakes. Advances in Geophysics, 45, 259-302
.
Wiemer, S., Wyss, M. (1997). Mapping the frequency-magnitude distribution in
asperities: an improved technique to calculate recurrence times?. Jour.
Geophys. Res. 102, 15115-15128.
Wiemer, S., Mc<utt, S.R., Wyss, M. (1998). Temporal and three-dimensional
spatial analysis of the frequency-magnitude distribution near Long Valley
caldera, California. Geophys. J. Int., 134, 409 421.
Weichert, D. (1980). Estimation of earthquake recurrence parameters for unequal
observations for different magnitudes. Bull. Seism. Soc. Am. 70, p. 13371346.

117
Bibliografa
Wyss, M., Wiemer, S., Zuiga, R. (2001). A tool for Analyses of Seismicity
Patterns. Tutorial de Zmap. P 3 63.
Wyss, M. (1973). Towards a physical understanding of the earthquake frequency
distribution. Geophys. J. R. Astron. Soc., 31, 341-359.
Wyss, M., Wiemer, S., Zuiga, R. (2000). Minimum magnitude of completeness in
earthquake catalog: Examples from Alaska, the western US and Japan. Bull
Seism. Soc. Am., 90, 859-869.
Wyss, M., <agamine, K., Klein, F.W., Wiemer, S. (2001). Evidence for magma at
intermediate crustal depth below Kilauea's East Rift, Hawaii, based on
anomalously high b-values. J. Volc. Geoth.Res., 106, 23-37.
Zuiga, R., Wyss, M. (2002). Mapping spatial variability of the frequencymagnitude distribution of earthquakes. Adv. 45, 259-302.
Zuiga, R., Wyss, M. (2001). Most and least likely locations of large to great
earthquakes along the pacific coast of Mexico, estimate from local recurrence
times based on b-values. Geophys. Jour. Intern. 90. in press, 2000.
Zuiga, F. R; Wyss, M. (1995). Inadvertent Changes in Magnitude Reported in
Earthquake Catalogs: Their Evaluation through b-Value Estimates. Bulletin
of the Seismological Society of America, Vol. 85, .o. 6, pp. 1858-1866.
Zuiga, F. R., Wiemer, S. (1999). Seismicity patterns: are they always related to
natural causes?, Pageoph, 155, 713-726.

118
Anexo A
AEXO A
Se adjuntan el resumen extendido del resultado preliminar de Investigacin

desarrollado por el autor, durante la etapa de elaboracin del presente estudio. El
mismo fue presentado en forma oral en el Congreso peruano de Geologa (2010).
Condori, C. y Tavera, H. (2010). Anlisis y Localizacin de Probables reas de
ruptura por sismos, en el borde occidental de Per: distribucin de valores de b. XV
Congreso Peruano de Geologa. Resmenes Extendidos. Sociedad Geolgica del
Per, Pub. Esp. & 9 (2010), Cusco p. 1193-1195

119
Anexo A

120
Anexo A

121
Anexo A

122
Anexo A

123

Tavera Condori 2011

Cargado por

Copyright:

Formatos disponibles

Tavera Condori 2011

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Tavera Condori 2011

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD NACIONAL DE SAN AGUSTN DE AREQUIPA

FACULTAD DE GEOLOGA, GEOFSICA Y MINAS

REAS PROBABLES DE RUPTURA SSMICA E EL BORDE

1.3 TRABAJOS ANTERIORES

CAPTULO II: CONTEXTO GEODINMICO DEL PER

2.1.1 Dorsal de Nazca

2.1.2 La Fractura de Mendaa

2.1.3 La Fosa Marina

2.1.4 La Cordillera Andina

2.1.5 La Cadena Volcnica

2.1.6 Los Sistemas de Fallas

2.2 ZONIFICACIN MORFOLGICA

2.2.1 La Cordillera Costanera

2.2.2 La Zona Costera

2.2.3 La Cordillera Occidental

2.2.4 Zona Intercordillerana

2.2.5 La Cordillera oriental

2.2.6 Zona Subandina

2.2.7 Llanura Amaznica

CAPTULO III: CARACTERSTICAS GENERALES DE LA SISMICIDAD EN

3.2 PERIODO INSTRUMENTAL

3.3 DISTRIBUCIN ESPACIAL

3.3.1 Sismicidad con foco Superficial

3.3.2 Sismicidad con foco Intermedio

3.3.3 Sismicidad con foco Profundo

3.4 PERFILES SSMICOS

3.4.1 Regin Norte

3.4.2 Regin Centro

3.4.3 Regin Sur

3.5 GEOMETRA DE LA SUBDUCCIN

3.6 DISTRIBUCIN DE REAS DE RUPTURA

CAPTULO IV: MODELOS DE PROCESOS DE RUPTURA

4.1.1 MODELO DE BARRERAS

4.1.2 MODELO DE ASPEREZAS

CAPTULO V: METODOLOGA PARA EL CLCULO DEL PARAMETRO b

5.1.1 Los Parmetros a y b

5.2 ASPEREZAS Y VALORES MNIMOS DEL PARAMETRO b

5.3 VARIACIONES DEL PARMETRO b

5.4 MTODOS PARA EL CALCULO DE a y b

5.4.1 Mtodo de Mnimos Cuadrados

5.4.2 Mtodo de Mximo Verosimilitud

5.5 METODOLOGIA PARA EL CALCULO DE a y b

5.6 PERIODOS DE RECURRENCIA

5.7 PROBABILIDAD DE RECURRENCIA

5.8 APLICACIONES DEL PARMETRO b

5.8.1 Estudios en las Zonas de Subduccin

5.8.1.1 Zona de Subduccin de Mxico

5.8.1.2 Zona de Subduccin Alaska y Nueva

5.8.1.3 Zona de Subduccin de Per

5.8.2 Sistemas de Fallas

5.8.3 Aplicacin en Volcanes

CAPTULO VI: ANLISIS Y SELECCIN DE DATOS (APLICACIN AL

6.2 HOMOGENEIZACIN DE LA MAGNITUD

6.3 MAGNITUD DE COMPLETITUD (Mc)

6.4 ELIMINACIN DE RPLICAS

6.4.1 Ley de Omori-Utsu

6.4.2 Algoritmo de Reasenberg

6.5 BASE FINAL DE DATOS SSMICOS

CAPTULO VII: DISTRIBUCIN DEL PARMETRO b (APLICACIN AL

REAS PROBABLES DE RUPTURA SSMICA E EL BORDE

1.3 TRABAJOS ATERIORES

COTEXTO GEODIMICO DEL PER

2.1 PRICIPALES RASGOS TECTICOS

2.2 ZOIFICACI MORFOLGICA

CARACTERISTICAS GEERALES DE LA SISMICIDAD E