SECRETARÍA DE GOBERNACIÓN CORDINACIÓN NACIONAL …...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

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