1
SECRETARÍA DE GOBERNACIÓN
CORDINACIÓN NACIONAL DE PROTECCIÓN CIVIL MÉXICO
CENTRO NACIONAL DE PREVENCIÓN DE DESASTRES
ANÁLISIS GEOFÍSICO DEL VOLCÁN POPOCATÉPETL
Dirección de Investigación
Subdirección de Riesgos Sísmicos
Jefatura de Investigación Aplicada a Riesgos Sísmicos
DICIEMBRE 2015
Elaboró: Sergio Alberto Galaviz Alonso
2
Se agradece la asesoría y colaboración por la adquisición, procesamiento
e interpretación de los datos obtenidos en este trabajo a los
investigadores: el M.C. Leobardo Salazar Peña, profesor del Instituto
Politécnico Nacional, al Dr. Carlos Valdés González, Director General
del CENAPRED, a la Ingeniera Geóloga Ariadna Hernández Oscoy de la
Subdirección de Riesgos Volcánicos de la Dirección de Investigación y
al Ing. Miguel Ángel Cruz Pliego, de la Subdirección de Sistemas de
Información de Riesgos de la Dirección de Análisis y Gestión de
Riesgos.
3
Créditos
Procesamiento de Imágenes
Ing. Sergio Alberto Galaviz Alonso
Tec. Lucio Cárdenas González
Simulación, Modelación y Animación Computacional
M.C. Leobardo Salazar Peña
Ing. Sergio Alberto Galaviz Alonso
Ing. Ariadna Hernández Oscoy
Adquisición de Datos Gravimétricos
Ing. Sergio Alberto Galaviz Alonso
M.C. Leobardo Salazar Peña
Ing. Ariadna Hernández Oscoy
Ing. Miguel Ángel Cruz Pliego
Procesamiento de Datos Gravimétricos
M.C. Leobardo Salazar Peña
Ing. Sergio Alberto Galaviz Alonso
M.C. Juan Carlos Jiménez Velázquez
Procesamiento de Datos Geodésicos
M.C. Juan Carlos Jiménez Velázquez
Ing. Bernabé Juárez García
Actualización de la Localización de Eventos Sísmicos
Ing. Gema Victoria Caballero Jiménez
Análisis de Datos Geofísicos
Ing. Sergio Alberto Galaviz Alonso
M.C. Leobardo Salazar Peña
4
ÍNDICE
1.- INTRODUCCIÓN……………………………………………………………………………..5
2.-IMÁGENES SATELITALES LANDSAT 7 ETM+ Y EL MODELO DIGITAL DE
ELEVACIONES…………………………………………………………………………….…….6
3.-SISMICIDAD REGIONAL.………………………………………......................................…9
4.-ANÁLISI DE LOS EVENTOS DE LA RED DE MONITOREO DDEL VOLCÁN
POPOCATÉPETL……….………………………………………………………………………11
5.- CAMPAÑAS DE ADQUISISCIÓN GRAVIMÉTICA……………………………………….19
ANEXO DE MAPAS ……………………………………………………………………………25
6.- PREPARACIÓN DE DATOS PARA LA SIMULACIÓN COMPUTACIONAL Y EL
MODELO GEOFÍSICO…………………………………………………………………………28
7.- ACCIONES A SEGUIR……………………………………………………………………..30
BIBLIOGRAFÍA…...…………….……………………………………………………………..31
5
1.- Introducción
El volcán Popocatépetl, ubicado en el Eje Neovolcánico Transversal Mexicano, en la parte austral
de la Sierra Nevada, con localización, 19.02° de latitud norte, 98.62° de longitud oeste y una
altura de 5,452 metros sobre el nivel del mar (m.s.n.m.), es uno de los volcanes más
monitoreados del país, principalmente por la etapa de reactivación que dio inicio el 21 de
diciembre de 1994, y por los más de 25 millones de personas que habitan en las inmediaciones de
su cráter, convirtiéndolo en uno de los volcanes más peligrosos del mundo.
La historia geológica del Popocatépetl, indica que su actividad se ha desarrollado en un lapso
mayor a medio millón de años y ha presentado etapas de crecimiento y destrucción de domos y,
al menos, la formación de tres volcanes previos, Nexpayantla, de aproximadamente 400,000
años, Ventorillo alrededor de 23,000 años y El Fraile de 14,500 años, así como de eventos
efusivos y explosivos, que han afectado severamente a los asentamientos humanos del pasado, así
como actualmente a las comunidades de los estados aledaños de México, Puebla y Morelos.
En la actualidad, para el monitoreo de la actividad del volcán Popocatépetl, se hace uso de
estaciones geodésicas de dos y tres componentes, estaciones sísmicas de periodo corto y de banda
ancha de tres componentes, estudios gravimétricos y magnetométricos, imágenes satelitales,
cámaras visuales y térmicas, detectores de gases y el análisis químico de muestras de
manantiales.
El presente trabajo hace uso de la información de las señales que nos proporcionan las imágenes
satelitales y visuales, los registros de las estaciones geodésicas y sísmicas, así como de cuatro
campañas de adquisición gravimétrica, tres por el Instituto Politécnico Nacional, (IPN) en
conjunto con el Centro Nacional de Prevención de Desastres (CENAPRED) y una por la
Universidad Autónoma de México (UNAM). El objetivo estriba en la identificación de las
estructuras sismogenéticas que la sismicidad histórica permite observar, así como la generación
de un modelo geológico-geofísico, el cual se use como insumo para la simulación computacional
el modelado y la animación 2D y 3D del fenómeno sísmico que permita comprender el estado
sismotectónico que rige la deformación y la configuración regional estructural .
6
2.- Imágenes Satelitales LANDSAT 7 ETM+ y el Modelo Digital de Elevaciones
En 2014, se hizo la instalación de la estación Encinos*, en el flaco NE del volcán Popocatépetl,
para ello se analizó la imagen satelital LANDSAT 7 ETM+ (Figura 1), mediante la técnica de
falso color, con el software Geomatica 2013. Con dicha técnica es posible identificar rasgos y
objetos que no se observan a simple vista (RGB 321), para ello se hicieron combinaciones de las
ocho bandas espectrales, generando imágenes de falso color, utilizando las combinaciones que
resaltan las estructuras geológicas (RGB 531, 135 y 652), la vegetación (RGB 432 y 741) y el
uso de suelo (RGB 652 y 752). En dicho análisis se identificó un rasgo estructural
correspondiente al plano de una falla que cruza al cráter con dirección NE45SW, siendo la falla el
objetivo de este análisis.
Figura 1.- Análisis de Imágenes Satelitales LANDSAT 7 ETM+, y MDE del INEGI
* PAT 2014 oficio 8_140128_RS_Istalación_de_la_estación_Encinos.pdf
7
A partir del 2015 se ha hecho uso del Modelo Digital de Elevaciones (MDE) de resolución
espacial de tres metros (Figura 2), del Instituto Nacional de Estadística, Geografía e Informática
(INEGI) para el análisis de las geoformas de las inmediaciones del volcán, para ello se usó el
software Generic Mapping Tools (GMT) y de tres paletas de colores para la identificación de los
principales rasgos de interés.
Figura 2.- Modelos 2D volcanes Popocatépetl e Iztaccíhuatl
Posteriormente se generó un modelo 3D con exageración en la vertical de 1.5, y se hizo un
recorte del volcán Popocatépetl generando vistas 3D cada grado para tener un giro de 360 sobre
la horizontal y de 90 sobre la vertical, pudiendo observar la sinuosidad de sus flancos, los
derrumbes del sur, la colada del noreste, las escorrentías y los estructuras que se encuentran
alrededor de los flacos (Figura 3).
Ambos insumos son utilizados para el análisis geomorfológico del terreno y la selección de áreas
de interés para el modelado computacional.
8
Figura 3.- Modelo 3D del Volcán Popocatépetl
9
Usando el MDE y la información de carácter abierta por el INEGI sobre los Recurso Naturales,
se generaron los mapas con las fallas y fracturas (Figura 4),
Figura 4.- MDE del Popocatépetl con fallas y fracturas del catálogo del INEGI
3.-Sismicidad Regional en la Zona Reportada por el SSN
De igual manera que con la fallas y fracturas, se combinó el MDE y los sismos reportados por el
Servicio Sismológico Nacional (SSN) desde 1998 a la fecha, para generar un modelo 2D (Figura
5) y 3D (Figura 6 y 7) de la sismicidad en el área de interés. Los circunferencias, representan la
sismicidad reportado de 1998 al 2005 y los círculos con relleno la sismicidad de 2006 a la fecha.
Los sismos de M0≥M<1M son de color azul, los de M1≥M<2M de color verde y los de
M2≥M<3M de color rojo. En la Figura 6 y 7 se puede apreciar la profundidad de ocurrencia; en
10
este caso en cubos con las mismas características que los círculos, que ha calculado el SSN. La
Figura 6 muestra la cara este de los volcanes, con un azimuth de 70 grados y la Figura 7 muestra
la cara oeste con azimuth de 250 grados.
Figura 5.- Mapa 2D de la sismicidad reportada por el SSN de 1998 al 2015
Figura 6.- Mapa 3D, Azimut de70°, con la sismicidad reportada por el SSN de 1998 al 2015
11
Figura 7.- Mapa 3D, Azimut de 250°, con la sismicidad reportada por el SSN de 1998 al 2015
4.-Análisi de los Eventos de la Red de Monitoreo Sísmico del Volcán Popocatépetl
En la actualidad el volcán Popocatépetl cuenta con una red de observación sísmica de registro
continuo, las estaciones se encuentran estratégicamente ubicadas en las inmediaciones del
edificio volcánico, en la tabla 2, se muestran las claves y coordenadas de las estaciones, en la
figura 8 se muestra la ubicación de las estaciones sísmicas de la Red de Monitoreo Sísmico del
Volcán Popocatépetl; los triángulos negros con contorno azul corresponden a sensores de periodo
corto y con contorno verde los sensores de banda ancha. Se ha utilizado una paleta de colores,
con el objetivo de resaltar la altura a las que se encuentran las estaciones de la red.
Tabla 2 Características de la red de observación sísmica de volcán Popocatépetl
Estación Clave Latitud Longitud Altura
(m.s.n.m.)
Banda Ancha Periodo Corto
Calo IIB 18.96 -98.34 2500 sensor triaxial
Canario PPP 19.04 -98.62 4313 sensor triaxial sensor triaxial
Chipiquixtle PPX 19.00 -98.65 3990 sensor triaxial sensor triaxial
Colibri PCC 18.98 -98.55 2670 sensor triaxial sensor triaxial
Encinos PPEN 19.04 -98.58 3655 sensor triaxial
12
Tetexcalo PPT 18.97 -98.62 3138 sensor triaxial
Tlamacas PPM 19.06 -98.62 4009 sensor triaxial sensor triaxial
Figura 8: Ubicación de la red de observación sísmica de volcán Popocatépetl
El registro sísmico volcánico, utilizado en este documento, data del dos de enero de 1995 hasta el
26 de diciembre de 2015, con un total de 2979 eventos volcanotectónicos localizados y referidos
a la estación Canario (PPP); cabe mencionar que 24 eventos no fueron incluidos debido a que
carecían de profundidad o magnitud. La tabla 2 muestra el número de eventos con magnitudes
por año y en el Anexo A se incluyen los datos de los eventos sísmicos registrados y localizados
por el personal del CENAPRED.
13
Tabla 2 Eventos Registrados por el CENAPRED de 1995 a 2015
Año No de VT´s M0≥M<M1 M1≥M<M2 M2≥M<M3 M3≥M<M4 M4≥M<M5
1995 107 NaN 15 92 NaN NaN
1996 166 NaN 22 139 5 NaN
1997 262 NaN 58 196 8 NaN
1998 152 NaN 44 104 4 NaN
1999 177 NaN 56 112 9 NaN
2000 144 NaN 36 107 1 NaN
2001 191 NaN 20 164 7 NaN
2002 275 NaN 33 239 3 NaN
2003 80 NaN 3 74 3 NaN
2004 60 NaN 7 53 NaN NaN
2005 95 NaN 6 89 NaN NaN
2006 138 NaN 6 130 2 NaN
2007 124 NaN 9 113 2 NaN
2008 78 NaN 5 72 1 NaN
2009 85 NaN 5 79 1 NaN
2010 46 NaN 3 43 0 NaN
2011 27 NaN 2 24 1 NaN
2012 95 NaN 75 18 2 NaN
2013 242 10 64 138 29 1
2014 241 NaN 60 159 22 NaN
2015 186 NaN 126 60 NaN NaN
Total 2971 10 655 2205 100 1
A continuación se muestra la distribución epicentral de los eventos sísmicos localizados por el
CENAPRED (Figura 9), en la imagen se observan dos zonas con cúmulos de eventos, la primera
de ellas ubicada debajo del cráter con mayor densidad de eventos y una segunda al sureste del
cráter el cual podría estar asociado a la actividad de la falla con rumbo NE45SW observada con
percepción remota y análisis geomorfológico. Los sismos de M0≥M<M1 (Figura 10) son de color
cian, de M1≥M<M2 (Figura 11) de color verde, de M2≥M<M3 (Figura 12) de color azul, de
14
M3≥M<M4 (Figura 13) de color rojo y de M4≥M<M5 (Figura 14) de color magenta.
Figura 9: Distribución espacial de los eventos volcanotectónicos del volcán Popocatépetl
Figura 10: Distribución espacial de los eventos volcanotectónicos M0≥M<M1 del volcán Popocatépetl
15
Figura 11: Distribución espacial de los eventos volcanotectónicos M1≥M<M2 del volcán Popocatépetl
Figura 12: Distribución espacial de los eventos volcanotectónicos M2≥M<M3 del volcán Popocatépetl
16
Figura 13: Distribución espacial de los eventos volcanotectónicos M3≥M<M4 del volcán Popocatépetll
Figura 14: Distribución espacial de los eventos volcanotectónicos M4≥M<M5 del volcán Popocatépetl
17
De las últimas cinco figuras se pueden observar cúmulos y lineamientos asociados a cada
magnitud; en la Figura10, correspondiente a eventos de M0≥M<M1, se observan dos
distribuciones, un cúmulo de eventos que se encuentran por debajo del cráter y dos eventos
aislados ubicados en un lineamiento estructural NW45SE; en la Figura 11, correspondiente a
eventos de M1≥M<M2, se observan una mayor cantidad de eventos con las mismas
características que los eventos de la Figura 10, aunado a esto se observa eventos con distribución
irregular en la porción sureste y norte del cráter. En la Figura 12 se observan los eventos de
M2≥M<M3, los cuales son los más abundantes, cuentan con dos cúmulos saturados ubicados
debajo y en la porción sureste del cráter, de igual manera se observan distribuciones lineales
norte-sur del cráter así como del lineamiento NW45SE y en menor cantidad en el lineamiento
NE45SW. La Figura 13, correspondiente a eventos de M3≥M<M4, muestra una clara
distribución espacial en dos cúmulos de los eventos calculados, el primero de ellos de bajo del
cráter y el segundo de ellos correspondiente al lineamiento estructural NE45SW, el cual cuenta
con una cantidad mayor de eventos en la porción sureste del cráter. La figura 14 muestra un
evento de magnitud 4.1 en la porción norte del cráter.
Después de analizar la distribución espacial en 2D se continuó con el análisis 3D de los
hipocentros en tiempo y su relación geomorfológica. Para ello se roto 360 grados los eventos de
la base del CENAPRED con el MDE del INEGI, clasificados por magnitud; cabe mencionar que
solo se tomaron en cuenta los eventos ubicados entre las latitudes 18.8° a 19.2°, longitudes -98.9
a -98.3 y altura -15000 a 4313, esta última altura, referida a la estación Canario, la cual es la
estación de referencia para la localización de la sismicidad registrada de la red.
Después de rotar las estructuras volcánicas con la sismicidad registrada, se observó que el
azimuth de 70° (Figura 15) y de 250° (Figura 16), correspondiente a su antípoda azimuthal,
muestran la distribución de cúmulos y lineamientos interrelacionados, se observa un cúmulo
asociado al cuerpo magmático debajo del cráter y el segundo cúmulo de sismicidad en la zona SE
del Popocatépetl y alineado a las estructuras geológicas denominadas Iztaccíhuatl. El Anexo B es
un video del modelo 3D y la sismicidad en tiempo el cual muestra la interdependencia sísmica
entre las estructuras geológicas ya mencionadas.
18
Figura 15.- Mapa 3D, Azimut de 70°, Eventos Sísmicos Localizados por el CENAPRED de 1995
al 2015.
Figura 16.- Mapa 3D, Azimut de 250°, Eventos Sísmicos Localizados por el CENAPRED de
1995 al 2015
19
5.- Campañas de Adquisición Gravimétrica
En 2014 se realizó el 1° levantamiento gravimétrico de la Subdirección de Riesgos Sísmicos de la
comunidad de San Nicolás de los Ranchos, Puebla, a la estación Encinos (Figura 17), con un
Gravímetro GC5 del Instituto Politécnico Nacional (IPN) y con la asesoría del M.C. Leobardo
Salazar Peña (Figura 18), con el objetivo de comenzar a obtener datos gravimétricos de las
inmediaciones del volcán, este procesamiento* fue realizado por el M.C. Juan Carlos Jiménez
Velázquez, Jefe de Departamento de esta subdirección, de esta subdirección.
Figura 17.- Primer Campaña de Adquisición de Datos Gravimétrico, S. Nicolas de los Ranchos-
Encinos.
*PAT 2014 oficio 6_140701_RS_Informe Gravimetria 2014.doc
20
Figura 18.- 1°Campaña de Adquisición de datos Gravimétricos por M.C. Leobardo Salazar Peña
(IPN) y Alberto Galaviz (CENAPRED).
Posteriormente el 15 de febrero de este año se realizó la adquisición de un segundo perfil
Amecameca, estado de México, al refugio de Tlamacas (Figura 19), también en conjunto con el
IPN (Figura 20)
Figura 19.- Segunda Campaña de Adquisición de Datos Gravimétrico, Amecameca-Tlamacas.
21
Figura 20.- 2° Campaña de Adquisición de datos Gravimétricos por M.C. Leobardo Salazar
Peña (IPN) y Alberto Galaviz (CENAPRED).
Finalmente el 29 y 30 de diciembre de 2015, con el objetivo de identificar la posible estructura
asociada a una falla con rumbo NE45SW se realizaron dos perfiles gravimétricos, el primero de
ellos denominad A-A’ de Amecameca, México, a Tétela del Volcán, Morelos y el segundo de
ellos de San Buena Aventura a Huejotzingo, Puebla denominado perfil B-B’ (Figura 21 y 22)
Figura 20.- Tercera Campaña de Adquisición de Datos Gravimétrico, Perfiles A-A’ y B-B’
22
Figura 21.- 3°Campaña de Adquisición de datos Gravimétricos por M.C. Leobardo Salazar
Peña (IPN), Alberto Galaviz, Ariadna Hernández Oscoy y Miguel Ángel Cruz Pliego
(CENAPRED)
Con los datos gravimétricos obtenidos, se procedió a generar un script para el procesamiento de
los datos observados; para ello se hizo la conversión de horas a horas fracción de los datos
adquiridos así como las horas de las bases, la corrección por deriva del instrumento, la corrección
de base, la gravedad teórica, la corrección por aire libre y la corrección de la loza de Bouger, para
poder sacar la anomalía de Bouger Simple. Los datos observados, y los pasos en el
procesamiento, (realizado con GAWK, Fortran y GMT) se incluyen en el ANEXO C.
En la Figura 22 correspondiente al perfil A-A’ se aprecian los datos observados, los datos con la
Anomalía de Bouger Simple y las elevaciones de las estaciones adquiridas; en dicha gráfica se
aprecia que la anomalía está asociada a una falla de escalón, entre los puntos 10 y 15 se
observaría el plano de la falla; mientras que en la Figura 23, correspondiente al perfil B-B’,
mismo proceso que en la figura anterior, se aprecia un salto en los datos observados sin embargo
al procesarlo se observa un aplanamiento entre los puntos 15 y 20; por lo que el levantamiento de
Ambos perfiles deberían de extenderse hacia el norte, así como el cierre entre A’-B y el cierre
entre B-A’, así como la adquisición en las rutas que permitan subir los flancos del volcán para
tener anomalías asociadas en el cuerpo del volcán; así como su adquisición anual para observar la
migración del magma en el tiempo, esto último con el objetivo de tener un modelo realista de la
geodinámica de los reservorios magmáticos, e incursionar en los modelos hidrodinámicos del
23
mismo.
Figura 22 Datos Observaos, Anomalía de Bouger Simple y Altura del Perfil A-A’.
24
Los datos gravimétricos fueron graficados en Surfer10 (Figura 23), y para su correcta
interpretación se consultó el fragmento de la carta de Anomalía de Bouger Simple y Aire Libre en
Mar de la República Mexicana (Figura 24), Así como del perfil N-S regional de Sandoval (2003)
et.al. (Figura 25). Finalmente se muestran los datos procesado incluyendo las tres campañas antes
mencionadas y una campaña más hecha por los estudiantes de la Facultad de Ingeniería de
Universidad Nacional Autónoma de México** (UNAM) (Figura 26).
Figura 23 Datos Observados, Anomalía de Bouger Simple y Altura del Perfil B-B’.
*PAT 2015
** PAT 2014 oficio 6_140701_RS_Informe Gravimetria 2014.doc
25
Figura 23.- Anomalía de Bouger Simple en 2D y 3D, de la Campaña Gravimétrica Perfiles A-A’
y B-B’
26
Figura 24.- Fragmentos del mapa de Anomalía de Bouger Simple en Continente y Aire Libre en
Mar de la República Mexicana, foto tomada del Instituto de Geofísica, UNAM, por M.C. Juan
Carlos Jiménez
Figura 25.- Perfil Gravimétrico del S-N de las inmediaciones del Volcán Popocatépetl Sandoval
2003et.al
27
Figura 26.- Anomalía de Bouger Simple de las campañas de adquisición gravimétrica
CENAPRED-IPN y de la campaña de adquisición de los estudiantes de la Facultad de
Ingeniería UNAM. Elavorado por el M.C. Juan Carlos Jiménez.
28
De la Figura 23 se puede observar la Anomalía de Bouger Simple asociada a la estructura
geológica identificada en 2014, y analizada en este 2015; derivado de este estudio geofísico, es
posible hacer una simulación y modelado computacional de la interrelación entre ambos cuerpos
sismogenéticos, la falla y el volcán Popocatéptl.
6.-Preparación de Datos para la Simulación Computacional y el Modelo Geofísico
Derivado del análisis 4D de los eventos sísmicos y la gravimetría procesada se ha podido
comenzar con la simulación computacional y la modelación sismotectónica. Para ello se utiliza el
programa maestro del M.C. Leobardo Salazar Peña para la generación del modelo de velocidades
que actualmente se usa para la localización de los eventos sísmicos (Figura 27). El cono tiene una
velocidad promedio de 3.5 Km/s y con un espesor de 5 km., adyacentemente se puesto una capa
de velocidad de 4.5 Km/s y espesor de 500 mts. y por último un espacio semi-infinito de velocida
de 5 Km/s.
Figura 27.- Modelo de Velocidades del Volcán Popocatéptl
29
Posteriormente se generaran las instantáneas alusivas a una explosión con pulso de Lamb, una
explosión interna u un evento tipo A. (Figura 28)
Figura 28.- Instantánea alusiva a un pulso de LAMB para la generación de sismogramas
sintéticos.
30
7.- Acciones a Seguir
En este 2015 se han obtenido los modelos simotectónicos de las inmediaciones del volcán
Popocatépetl, se ha generado conocimiento del proceso geodinámico que rige la zona del volcán,
así como la generación de software que permite el procesamiento de datos gravimétricos y la
generación de señales sísmicas sintéticas del volcán, lo que permitirá en un futuro la antelación
ante la ocurrencia de los eventos sísmicos y la relación entre la falla normal escalonada con
componente lateral izquierda y la evolución de los volcanes situados en la Sierra Nevada.
Las acciones a seguir son las siguientes:
Extensión de los perfiles gravimétricos hacia el norte y sur, para comprobar el corrimiento
de la Sierra Nevada.
Adquisición anual de los levantamientos gravimétricos en el cono del volcán, como ya se
ha hecho con Encinos y Tlamacas, para el mapeo de la migración de magmas.
Generación del modelo de la falla normal con componente lateral izquierdo.
Generación de sismogramas Sintéticos con el objetivo de identificar dicha estructura de
carácter cortical y su magnitud asociada.
Mapa de peligro sísmico por sismo de tipo normal.
Generación de estudios de frontera en materia de adquisición de datos sismotectónicos,
monitoreo y modelación de estratovolcanos.
Comparación de las señales sísmicas con las sintéticas para la clasificación de las mismas.
Simulación, modelación y animación en tiempo cuasi real de la actividad volcánica.
Artículo Internacional: Estructura Sismotectónica de las Inmediaciones del Volcán
Popocatépetl, Implicaciones Sobre Peligro Sísmico.
Paleomapas y paleosismología de la Sierra Nevada
31
BIBLIOGRAFÍA CONSULTADA
Recursos Naturales Fallas y Fracturas
http://www.inegi.org.mx/geo/contenidos/recnat/geologia/infoescala.aspx
Catálogo de Sismos
http://www2.ssn.unam.mx:8080/website/jsp/catalogo1.jsp
Alfonso Muños-Martín, Gerardo de Vicente, (2010) Análisis de Esfuerzos tectónicos
Fallas y Sismos RE Simotectónica ISSN:1986-6557
Viniegra Heberlein Germín, Mecánica de los Cuerpos Deformables, 1ª edición 2011
isbs:978-607-02-2102-6
Mercier-Vergely Tectónica, Limusa edición 2010
Salazar L. (2004) Simulación Numérica de la Propagación de Ondas en Medios
Irregulares: Aplicaciones a Estructuras Volcánicas y Valles Aluviales