HECTOR WILLIAM ARBOLEDA
INGENIERIA TOPOGRAFICA
UNIVERSIDAD DEL VALLE
CALCULO Y ZONIFICACION DE LOS TENSORES DE ESFUERZOS GENERADOS
POR LA PLACA NAZCA EN EL DEPARTAMENTO DEL VALLE DEL CAUCA,
COLOMBIA
1
2
Dedicado a
Mis padres, Héctor David Arboleda Millán y Graciela Orejuela Valero, por haberme
apoyado siempre, por los valores inculcados que me permiten ser una persona de
bien y sobre todo por su amor…
A mi esposa Carolina Perdomo Vargas y a mis hijos Fabiana y Sergio, por ser la
energía que impulsa mi vida y por llenarla de felicidad, a Carolina por su entrega y
compañía, a Fabiana por enseñarme que existen diferentes formas de querer y a
Sergio por llenar mi corazón de ilusión, los AMO profundamente.
Agradezco
A mi hermano José Danilo Arboleda Orejuela por su apoyo y ejemplo, a Trinidad
Vargas Becerra y Etelvina Becerra por su cariño, por recibirme con los brazos
abiertos en su hogar, y brindarme cobijo como el de una madre.
A mis abuelos Julio, Graciela, Claudio y Julia por su compañía, cariño e historias.
A la profesora Olga Baquero su apoyo, disposición de enseñar y su dirección en la
ejecución de este trabajo de grado, a Gustavo Salazar por resolver mis
inquietudes, a mis familiares, amigos y compañeros por sus palabras de aliento.
A todos Gracias.
3
Contenido
Índice de ilustraciones _________________________________________________________ 4
Índice de tablas _______________________________________________________________ 5
1. Introducción ______________________________________________________________ 7
1.1. Planteamiento del problema ____________________________________________ 8
1.2. Antecedentes _________________________________________________________ 8
1.3. Objetivos ____________________________________________________________ 10
1.1.1. Objetivo General _________________________________________________ 10
1.1.2. Objetivos Específicos ____________________________________________ 11
1.4. Justificación _________________________________________________________ 11
2. Marco teórico ____________________________________________________________ 12
2.1. Geología estructural ____________________________________________________ 12
2.2. Mecanismos focales ________________________________________________ 13
2.3. Esfuerzos Tectónicos _______________________________________________ 17
2.4. Régimen de esfuerzos ______________________________________________ 19
2.5. Modelos de inversión _______________________________________________ 20
2.6. Sistemas de Información Geográfica _________________________________ 21
2.7. Software Libre _____________________________________________________ 22
3. Inversión amortiguada de esfuerzos tectónicos ____________________________ 22
3.1. Análisis de mecanismos focales de terremotos como modelo de
inversión. ________________________________________________________________ 22
4. Desarrollo del proyecto ___________________________________________________ 29
4.1. Descripción de la zona de estudio _____________________________________ 29
4.2. Distribución de mecanismos focales y la geología estructural. ___________ 31
4.3. Geología estructural de la zona de estudio. ____________________________ 31
4.3.1. Zona Subducción ________________________________________________ 31
4.3.2. Naya-Micay ______________________________________________________ 32
4.3.3. Falla El Tambor __________________________________________________ 33
4.3.4. Sistema de Falla Dagua Calima ____________________________________ 33
4.3.5. Falla Garrapatas __________________________________________________ 34
4.3.6. Argelia __________________________________________________________ 34
4.3.7. Falla Toro ________________________________________________________ 35
4.3.8. Falla Ibagué ______________________________________________________ 35
4
4.3.9. Falla Palestina ___________________________________________________ 35
4.3.10. Falla Cucuana __________________________________________________ 36
4.3.11. Sistema de Falla Romeral _______________________________________ 36
4.3.12. Falla Montenegro _______________________________________________ 37
4.3.13. Falla Armenia __________________________________________________ 37
4.3.14. Falla Paraíso ___________________________________________________ 38
4.3.15. Falla Piendamo _________________________________________________ 38
4.3.16. Falla Córdoba – Navarco ________________________________________ 39
4.4. Mecanismos focales. _________________________________________________ 39
4.5. Mapa mecanismos focales vs fallas. ___________________________________ 49
4.6. Aplicación del modelo de inversión amortiguada. _______________________ 54
4.7. Mapa tensor de esfuerzos _____________________________________________ 55
4.8. Clasificación de los tensores de esfuerzos _____________________________ 67
5. Conclusiones ____________________________________________________________ 70
6. Bibliografía ______________________________________________________________ 71
ANEXOS _____________________________________________________________________ 73
Índice de ilustraciones
Ilustración 1 Relación entre el esfuerzo y la deformación a la que se somete un medio elástico, alcanzando
una deformación elástica, deformación plástica llegando al límite de ruptura. _______________________ 12 Ilustración 2 Sistema de fallas conjugadas, en relación con los tipos de esfuerzo del Modelo de fracturación
de Anderson: A) fallas normales, B) fallas en dirección (desgarres) y C) fallas inversas (cabalgamientos).__ 13 Ilustración 3 Diagrama de una Mecanismo Focal _______________________________________________ 14 Ilustración 4 Registro del primer movimiento de la onda P, Fuente Matt Miller, 2014 _________________ 14 Ilustración 5 Representación Estereográfica del Mecanismo Focal mediante la esfera de Schmidt, a)
Representación estereográfica de los datos sísmicos, b) Trazados del gran circulo en una red estereográfica,
c) Identificación de los cuadrantes. Fuente Matt Miller, 2014 _____________________________________ 15 Ilustración 6 Componentes del Tensor de Esfuerzos, fuente REDUCA 2010 __________________________ 18 Ilustración 7 Representación de Superficie de Falla, fuente D. Legrand, U. Chile ______________________ 18 Ilustración 8 Régimen de Esfuerzos, Tomado de World Stress Map. ________________________________ 20 Ilustración 9 Diagrama de Variables en el Bloque de Masa. Héctor William Arboleda __________________ 24 Ilustración 10 Zona de Estudio, Fuente: OSSO 2013 _____________________________________________ 30 Ilustración 11 Mapa tectónico y fallas de la zona de estudio. Trenkam et.al. _________________________ 32 Ilustración 12 Consulta de Datos Catalogo CMT ________________________________________________ 40 Ilustración 13 Tabla de Atributos de Mecanismos Focales utilizados en el desarrollo del estudio. ________ 50 Ilustración 14 Tabla de Atributos de Fallas ____________________________________________________ 51 Ilustración 15 Mapa Mecanismos Focales VS Fallas _____________________________________________ 52 Ilustración 16 Formato entrada datos Mecanismo a GMT _______________________________________ 53
5
Ilustración 17 Mapa Mecanismos Focales, Software GTM________________________________________ 54 Ilustración 18 Interfaz gráfica catalogo CMT __________________________________________________ 55 Ilustración 19 Curva de Compensación para Grilla 0,1°. Grafica generada en Software GMT ____________ 59 Ilustración 20 Parte 1 archivo de salida SATSI, parámetro de amortiguamiento 0,8 ___________________ 59 Ilustración 21 Parte 2 archivo de salida SATSI, parámetro de amortiguamiento 0,8 ___________________ 60 Ilustración 22 Parte 3 archivo de salida SATSI, parámetro de amortiguamiento 0,8 ___________________ 60 Ilustración 23 Componentes principales para Tensores de Esfuerzos Calculados sin parámetro de
amortiguamiento, para el plano de falla 1. ___________________________________________________ 64 Ilustración 24 Componentes principales para Tensores de Esfuerzos Calculados sin parámetro de
amortiguamiento, para el plano de falla 2. ___________________________________________________ 64 Ilustración 25 Componentes principales para Tensores de Esfuerzos Calculados con un parámetro de
amortiguamiento de 0,8, para el plano de falla 1. ______________________________________________ 65 Ilustración 26 Componentes principales para Tensores de Esfuerzos Calculados con un parámetro de
amortiguamiento de 0,8, para el plano de falla 2. ______________________________________________ 65 Ilustración 27 Componentes principales para Tensores de Esfuerzos Calculados con un parámetro de
amortiguamiento de 3, para el plano de falla 1. _______________________________________________ 66 Ilustración 28 Componentes principales para Tensores de Esfuerzos Calculados con un parámetro de
amortiguamiento de 3, para el plano de falla 2. _______________________________________________ 66 Ilustración 29 Herramienta Grouping Analysis _________________________________________________ 67 Ilustración 30 Mapa de Zonificacion de Tensores de Esfuerzos. ___________________________________ 68 Ilustración 31 Tabla Clasificación de Tensores _________________________________________________ 69
Índice de tablas
Tabla 2.1.1 Ejemplo Solución de un Mecanismo Focal ..................................................................................... 16 Tabla 4.1 Descripción de la Zona de Subducción .............................................................................................. 32 Tabla 4.2 Descripción Falla Naya-Micay ........................................................................................................... 33 Tabla 4.3 Descripción Falla El Tambor .............................................................................................................. 33 Tabla 4.4 Descripción Sistema de Falla Dagua Calima ..................................................................................... 34 Tabla 4.5 Descripción Falla Garrapatas ............................................................................................................ 34 Tabla 4.6 Descripción Falla Argelia................................................................................................................... 34 Tabla 4.7 Descripción Falla Toro ....................................................................................................................... 35 Tabla 4.8 Descripción Falla Ibagué ................................................................................................................... 35 Tabla 4.9 Descripción Falla Palestina ............................................................................................................... 36 Tabla 4.10 Descripción Falla Cucuana .............................................................................................................. 36 Tabla 4.11 Descripción Sistema de Falla Romeral ............................................................................................ 37 Tabla 4.12 Descripción Falla Montenegro ........................................................................................................ 37 Tabla 4.13 Descripción Falla Armenia .............................................................................................................. 37 Tabla 4.14 Descripción Falla Paraíso ................................................................................................................ 38 Tabla 4.15 Descripción Falla Piendamo ............................................................................................................ 38 Tabla 4.16 Descripción Falla Córdova - Navarco .............................................................................................. 39 Tabla 4.17 Mecnismos Focales desde 1976 hasta 2018 .................................................................................. 41 Tabla 4.18 Mecanismos Focales ....................................................................................................................... 56 Tabla 4.19 Tabla de Datos para Parámetro de Amortiguación ........................................................................ 57 Tabla 4.20 Tensor de Esfuerzos ........................................................................................................................ 60 Tabla 4.21 Componentes para graficar tensores de esfuerzos mediante la variable –Sx del Software GMT. . 62
6
7
1. Introducción
El departamento del Valle del Cauca se localiza sobre una zona sísmica activa, debido a
la interacción fundamentalmente entre las placas tectónicas Nazca y Sur América. La
placa Nazca subduce bajo la placa Sur América, la cual se mueve a una velocidad
aproximada de 58 mm/año en una dirección W-E; dicho movimiento hace que se
produzca una zona de alta complejidad tectónica, representada en grandes
deformaciones asociadas a fallas activas, que finalmente son las responsables de los
sismos ocurridos en la región.
La interacción entre las placas genera esfuerzos tectónicos, podemos considerar la
corteza terrestre como un medio elástico que acumula energía, hasta que la fuerza de
resistencia de la roca es superada y se rompe liberándola, esta energía que pude ser
distribuida en formas de onda, la cual cuando se librera y llega a la superficie es
monitoreada y analizada mediante estaciones de sensores sísmicos, un producto de este
análisis son los mecanismos focales que describen como, posiblemente, se generó la
ruptura y liberación de los esfuerzos generados sobre la masa de roca.
Los esfuerzos en la corteza terrestre se afectan por procesos tectónicos de acumulación,
así como por los cambios de esfuerzos causados por grandes terremotos. Los estudios
regionales de estrés cortical y su variación espacio-temporal pueden proporcionar una
idea de la mecánica sísmica, la carga tectónica y las interacciones de estrés sísmico Se
realizó un mapa de esfuerzos del suroccidente en el Valle del Cauca que representa la
variación en la orientación de los esfuerzos principales y/o su factor de forma que es
activo en un determinado intervalo de tiempo en una zona determinada. (Olaiz 2006).
El estudio de los esfuerzos, de la deformación acumulada y de la actividad sísmica son de
gran importancia para entender tanto los mecanismos de la tectónica de placas en un
área determinada como los sismos generados, permitiendo determinar fuentes,
comportamientos, magnitudes probables y así establecer planteamientos hipotéticos de
generación de sismos y es muy importante para la evaluación de la peligrosidad sísmica
Para este estudio se obtuvieron del Catálogo mundial CMT, Centroid-Moment-Tensor por
sus siglas en inglés 43 mecanismos focales, calculados a partir sismos generados en el
Valle del Cauca y zonas aledañas, los mecanismos focales tienen información sobre la
deformación y mediante solución del problema inverso se obtuvieron los tensores de
8
esfuerzo. El cálculo se hizo en una grilla y para evitar discontinuidades espaciales de
estos tensores se hace una inversión amortiguada que minimiza las diferencias de los
resultados entre sub-áreas adyacentes y de esta manera se realizó un mapa de
esfuerzos. De esta manera se podrá obtener un mapa de esfuerzos que muestra las
direcciones en el plano horizontal de los esfuerzos principales.
1.1. Planteamiento del problema
El Valle del Cauca localizado al Sur Occidente Colombiano se encuentra sobre la zona de
subducción de la placa Nazca en la placa Sur Americana; la interacción de estas placas
generan esfuerzos tectónicos al interior de la corteza y una sismicidad alta que a su vez
puede modificar el régimen de esfuerzos tectónicos al interior de la corteza terrestre. Los
sismos generados son producto de la fracturación y movimiento de las rocas cuando los
esfuerzos superan el régimen elástico y plástico de estas. Esta deformación de las rocas
puede modelarse como mecanismos focales o a partir del cálculo del tensor momento a
partir de los registros sísmicos. Esta deformación está directamente relacionada con los
esfuerzos, es posible entonces, a partir de los mecanismos focales de los sismos de una
zona calcular el régimen de esfuerzos Para poder establecer las condiciones actuales del
régimen de esfuerzos se propone aplicar el modelo de inversión amortiguada (Hardebeck,
2006), para obtener la distribución espacial de los esfuerzos tectónicos.
1.2. Antecedentes
Tabares et al. (1999). demostró que mediante el cálculo del Tensor de Momento Sísmico
(12 mecanismos focales) asociados a igual número de sismos, se pueden calcular las
velocidades de desplazamiento en cada uno de los ejes, permitiendo asociar la actividad
sísmica de la región al movimiento de placas. Lo aplicó específicamente al movimiento de
subducción de la placa de Nazca bajo la placa de Sur América, concluyendo que esta
placa se mueve a 41mm/año, con un azimut de 133
Vargas et.al. en 2002 plantean que en Colombia se pueden observar tres dominios
sísmicos y fisiográficos diferentes: Llanos Oriéntales, Región Andina y Región Caribe y
que, la principal actividad sísmica se presenta en la Región Andina que la comprenden un
9
cinturón montañoso que se divide en tres cordilleras, Oriental, Central y Occidental.
Estimaron el estado de deformación y esfuerzos en el territorio colombiano realizando un
análisis de la sismicidad con datos obtenidos por INGEOMINAS entre 1992 y 1999,
correlacionada con datos de GPS y 96 soluciones de Mecanismos Focales de la
Universidad de Harvard para el periodo 1976 a 2000 con magnitudes Mw≥ 5 realizando la
inversión de los datos de estos mecanismos con el método de Reches que permite
seleccionar el tensor de esfuerzos adecuado entre lo que se calculan para diferentes
posibles coeficientes de fricción de la roca. Finalmente concluyeron que la placa Nazca y
Caribe ejercen esfuerzos compresivos sobre la placa Sur Americana, además la
presencia de tres bloques o micro placas: Panamá-Costa Rica, Andes y Choco y que el
sistema de falla del bloque andino y el proceso de subducción de las placas Nazca y
Caribe son generadores de la actividad sísmica de la región.
En el 2005 INGEOMINAS presenta los resultados del estudio de microzonificación
sísmica para la ciudad de Santiago de Cali, en dicho estudio realizaron un análisis
probabilístico de la amenaza sísmica, para la ejecución de este análisis fue necesario
entre otros hacer una caracterización sismo-tectónica se requiere un modelo sismo-
tectónico determinado por la geometría de las fuentes y sus características como lo son:
la localización, magnitud, profundidad. Esto lo lograron realizando un levantamiento neo-
tectónico detallado en el que basados en cartografía geológica base, aerofotografías en
las que se identificaron indicadores cinemáticos y que posteriormente se realizaron visitas
de campo, estableciendo la geología, el marco tectónico regional y la tectónica activa,
para con este último determinar los campos de esfuerzos. Se analizaron de 277
mecanismos focales suministrados por el catalogo CMT (Centroid-Moment Tensor) que
permitieron calcular las trayectorias de los esfuerzos con radios de interpolación altos ya
que la distribución de los datos no es homogénea, también se realizó el cálculo del tensor
de esfuerzos mediante el método de diedros rectos e inversión directa para poblaciones
de estrías de falla.
Olaiz (2006), realiza sobre la península ibérica una inversión de los esfuerzos activos
mediante mecanismos focales calculados a partir del momento tensor sísmico. En este
trabajo utilizo la metodología para mecanismos focales calculados a partir de las primeras
llegadas, de este modo la solución en la inversión es mejor y explica un mayor número de
fallas, aunque los resultados tensoriales son similares tanto si se introducen solo los
planos reactivados o los dos (reactivado y neo-formado). Posteriormente se procede a
10
calcular mediante el método de deslizamiento los valores de máximo acortamiento
horizontal y el factor de forma (Dey y K`), se analizó la distribución espacial de los
mecanismos focales en función del factor de forma realizando una interpolación de los
valores de K` para obtener mapas de trayectoria y régimen de esfuerzos. Luego se
zonifica la zona de estudio en 6 áreas con interés para la inversión del tensor de
esfuerzos con el método de inversión de esfuerzos de Reches et al. (1992) para realizar
un posterior análisis de los resultados obtenidos.
Entre los diferentes métodos de cálculo para el tensor de esfuerzo, Hardebeck et al.
(2006), propone un nuevo modelo de inversión de los esfuerzos derivados de
mecanismos focales para usarlo en regiones, basada en las orientaciones de los
esfuerzos y utilizando un parámetro de amortiguamiento; divide la zona de estudio en sub-
áreas más pequeñas y realiza el cálculo independiente de cada sub-área de los
esfuerzos, lo que puede producir una variabilidad aparente entre las sub-áreas que es
más un error de adaptación de los datos ruidosos y para solucionar este problema se ideo
un método de inversión amortiguado que le permitió simultáneamente invertir los tensores
de esfuerzos mientras que se mantienen la variabilidad de cada sub-área. La
amortiguación remueve los errores por rotaciones, mientras resuelve estrés de rotación
mejor que un modelo suavizado o ventana en movimiento. Este método, se fundamenta
en uno de los métodos de cálculo del tensor de esfuerzos único que mejor se ajusta a una
población de mecanismos focales (Michael 1984), método elegido por que utiliza una
inversión lineal.
1.3. Objetivos
1.1.1. Objetivo General
Realizar el cálculo y zonificación de los esfuerzos tectónicos el Valle del Cauca y zonas
aledañas debido a la Plazca Nazca.
11
1.1.2. Objetivos Específicos
1. Hacer un mapa con distribución espacial de los Mecanismo Focales y la geología
estructural de la Zona de Interés.
2. Realizar la inversión de los mecanismos focales aplicando el método de Inversión
Amortiguada de Esfuerzos.
3. Realizar la zonificación de los tensores de esfuerzos tectónicos.
1.4. Justificación
En las últimas décadas el estudio del riesgo sísmico de región ha sido abordado desde
muchas disciplinas y enfoques, especialmente las vinculadas con las Ciencias de la
Tierra, la Ingeniería y las ciencias sociales con la vulnerabilidad social. Debido a que no
es posible predecir la ocurrencia de un terremoto, todos los estudios se encaminan en la
mitigación del riesgo asociado a un terremoto para lo cual es importante tener el régimen
de esfuerzos para ver la acumulación y liberación de esfuerzos dentro del límite de las
placas tectónicas, y evaluar la peligrosidad sísmica.
En esta línea, una de las variables fundamentales del peligro es el estado de esfuerzos
sobre una región. La distribución de esfuerzos tectónicos demuestra de manera general
los procesos tectónicos al interior de la corteza, que se ven afectados por la ocurrencia de
sismos y/o procesos tectónicos de acumulación. Los estudios regionales de esfuerzos en
la corteza, tanto espaciales como temporales pueden dar una idea de la mecánica de los
terremotos, carga tectónica y las interacciones de los esfuerzos del terremoto (Hardebeck
y Michael 2006). Mediante análisis posterior de la distribución de los esfuerzos tectónicos
se podrá establecer de manera detallada la Amenaza sísmica de la región, que
posteriormente llevara a establecer el riesgo sísmico.
12
2. Marco teórico
2.1. Geología estructural
La Geología estructural estudia, fundamentalmente, las rocas deformadas. La
deformación es la respuesta de las rocas cuando se las somete a un determinado estado
de esfuerzos.
Ilustración 1 Relación entre el esfuerzo y la deformación a la que se somete un medio elástico, alcanzando una deformación elástica, deformación plástica llegando al límite de ruptura.
Una falla es una discontinuidad que se forma debido a la fractura de grandes bloques de
rocas en la Tierra cuando las fuerzas tectónicas superan la resistencia de las rocas
(INPRES 2014). Esta discontinuidad puede generar movimientos verticales, horizontales
o combinación de estos, el Valle del Cauca presenta una compleja condición estructural
en donde “la mayoría de las rocas se encuentran falladas y diaclasadas” (IGAC 2009) en
diversas direcciones.
Generalmente estas fallas se presentan en los límites de las placas tectónicas, pero
también en los lugares donde la fuerza haya superado la resistencia del material rocoso y
se rompa. Existen tres tipos de falla a los que se asocian los movimientos de ruptura,
falla normal, falla inversa y falla de transformación.
13
Ilustración 2 Sistema de fallas conjugadas, en relación con los tipos de esfuerzo del Modelo de fracturación de Anderson: A) fallas normales, B) fallas en dirección (desgarres) y C) fallas inversas (cabalgamientos).
2.2. Mecanismos focales
El Mecanismo focal derivado de la ocurrencia de un sismo proporciona información
importante entre los que se incluye el tiempo de origen, la ubicación del epicentro, la
profundidad local, el momento sísmico, magnitud y orientación espacial, descrita por la
14
orientación de la falla y la dirección del movimiento sobre la falla, estas últimas
componentes se pueden observar gráficamente mediante un balón de playa, que muestra
dos soluciones del mecanismo focal, estas dos soluciones están dadas por dos planos, un
plano de falla y uno auxiliar, por lo que se puede relacionar con el tipo de falla de la región
donde se genera el sismo.
Los diagramas de balón de playa son representación estereográfica que muestran dos
cuadrantes negros y dos blancos, separados por dos círculos perpendiculares entre sí,
que serán los planos nodales (plano de falla y plano auxiliar), según sea la orientación de
estos planos respecto al norte (orientación de la falla), el ángulo de la falla con respecto a
la horizontal (superficie) y el ángulo del movimiento sobre el plano de falla, se pueden
asociar estos mecanismos a un tipo de falla.
Ilustración 3 Diagrama de una Mecanismo Focal
La solución del mecanismo focal se puede obtener mediante el análisis de la llega de
ondas. Cuando se genera un sismo, se libera energía en forma de onda, que recorre la
corteza terrestre, desde el foco hasta que arriba a los sismogramas dispuestos en la
superficie, a estos llegan dos tipos de onda, denominadas P y T, las que permiten
determinar si al sismograma llegan ondas producto de fuerzas de compresión o distensión
respectivamente, para establecer el mecanismo focal se analizan las llegadas de la onda
P, observando si el registro de la componente vertical en cada sismograma, si el primer
movimiento detectado fue arriba, si fue abajo o sin señal aparente (ilustración 4).
SIN SEÑAL ARRIBA ABAJO
Ilustración 4 Registro del primer movimiento de la onda P, Fuente Matt Miller, 2014
Lateral Normal Inversa Inversa-Oblicua
15
Para representar los datos de las llegadas de onda a los sismogramas se utiliza la red
estereográfica de Schmidt. Cada señal de llegada a cada sismograma será representado
por uno de los tres símbolos establecidos de la siguiente manera, un circulo en blanco
para las llegadas abajo, un circulo negro para las llegadas arriba y una x si el primer
movimiento de llegada no se logra diferenciar, el símbolo es localizado sobre una línea
que se extiende desde el centro hacia el azimut de la estación relativa al foco del sismo, y
el ángulo de despegue define la distancia angular desde el centro al símbolo. Para
observar gráficamente el proceso se presenta el siguiente ejemplo tomado de documento
notas de clase Sismología Aplicada de Matt Miller (2014). En el que el símbolo asociado
con la estación A está a 60° del centro, a lo largo de una línea dirigida hacia el azimut 50°
(Ilustración 4. Ejemplo de Representación Estereográfica).
Ilustración 5 Representación Estereográfica del Mecanismo Focal mediante la esfera de Schmidt, a) Representación estereográfica de los datos sísmicos, b) Trazados del gran circulo en una red estereográfica, c) Identificación de los
cuadrantes. Fuente Matt Miller, 2014
Luego de representar los datos sísmicos en la esfera estereográfica, se puede identificar
dos arcos que conforman dos círculos perpendiculares entre sí, separando los círculos
negros de los blancos, pasando cerca o sobre los símbolos x, estos serán los planos
nodales, uno será denominado plano de falla y el otro plano auxiliar, dependiendo de la
geología estructural del sector.
La solución geométrica del mecanismo focal está dada por tres elementos que permiten
especificar: la orientación del plano de falla, ángulo del plano de falla con respecto al
norte, el ángulo de buzamiento, ángulo del plano de falla con respecto a la horizontal y el
16
ángulo del deslizamiento del bloque sobre el plano de falla, estos ángulos son
normalmente denominados Strike, Dip y Slip.
Tabla 2.1.1 Ejemplo Solución de un Mecanismo Focal
Mecanismo 1 Región SUR PANAMA Fecha 1985/06/10
060185A
Latitud Longitud Profundidad Magnitud
3,24 -78,99 26 5,5
Solución Strike Dip Slip
Plano de Falla 32 19 125
Plano Auxiliar 176 74 79
Estos mecanismos focales son obtenidos de catálogos de datos, entre las que se
encuentran: Centro Internacional de Sismología (ISC por sus siglas en ingles), el Boletín
del Centro Internacional de Sismología es la primera salida del ISC y es considerado
como un registro definitivo de la sismología de la tierra, Este boletín se contiene datos
desde 1990 hasta el día presente, el boletín es revisado manualmente por los analistas
del ISC y usualmente tiene un retraso de 24 meses, este boletín es la contribución de
varias agencias de sismología de todo el mundo, hasta la fecha 565 agencias han
contribuido con este boletín a lo largo de su historia.
Otro catalogo el del Centro de Datos de Terremotos del Sur de California (SCEDC por sus
siglas en ingles), este centro opera en el Laboratorio de Sismología de Caltech y es el
principal archivo de datos sismológicos para el sur de California, cuenta con datos desde
1932 hasta la fecha, los cuales inicialmente son del Caltech y el USGS (Servicio
Geológico de los Estados Unidos) y es mantenido por el SCEDC, su misión es garantizar
la accesibilidad a los datos sísmicos para cualquier región de los Estados Unidos.
El proyecto Global Centroid-Moment-Tensor (CMT) supervisado por el investigador
principal Göran Ekström en el Observatorio Terrestre Laomnt-Doherty (LDEO) de la
Universidad de Columbia, este proyecto fue fundado por Adam Dziewonski en la
Universidad de Harvard y funciono allí como el proyecto Harvard CMT desde 1982 hasta
el 2006, este proyecto de investigación avanza en el mundo con el nombre de “The Global
CMT Project”, este catálogo cuenta con información sísmica que data desde 1976 hasta la
actualidad, con más de 25.000 tensores de momento para terremotos de todo el mundo.
17
2.3. Esfuerzos Tectónicos
Los esfuerzos en la corteza terrestre son afectados por procesos de carga tectónica, así
como los cambios de esfuerzos por grandes terremotos (Hardebeck 2006). Estos
procesos de carga y cambios de esfuerzos generan una serie de fuerzas primarias que se
conducen a través de varios kilómetros al interior de la corteza terrestres y son
evidenciados mediante movimientos y deformaciones. El resultado del balance de estos
esfuerzos en una zona determinada, controla la dinámica de deformación de dicha área
(Olaiz 2006).
El esfuerzo (σ) se define como la fuerza (F) por unidad de área (A).
𝜎 = 𝐹→
𝐴→
(2,3,1)
En el caso de la corteza terrestre que es un medio elástico, se habla de tensor de
esfuerzo, al aplicar una fuerza sobre el medio rocoso y continuar aplicando cada vez más
y más fuerza este llegara al punto de ruptura lo que daría origen a una falla que se
rompería en una dirección según la dirección de la fuerza, es por eso que el tensor de
esfuerzos se define con direcciones. El tensor de esfuerzos puede ser estudiado en 1D,
2D o 3D en el que los esfuerzos se presentan como un tensor. Un tensor de esfuerzos es
un operador vectorial lineal de N componentes y el número de componentes viene
definido por la expresión:
𝑁 = 𝑛𝑚 (2,3,2)
Donde n es la numero de dimensiones en el espacio y m es el orden del tensor.
Así mismo, los esfuerzos de la corteza son también fuerzas impulsoras de terremotos. De
manera que el estudio de estos esfuerzos de la corteza y sus variaciones espaciales y
temporales pueden proveer información acerca de la mecánica del terremoto, la carga
tectónica y las interacciones del esfuerzo de los terremotos.
Para determinar el estado de los esfuerzos se deben definir orientados en tres
dimensiones ortogonales entre sí, como está definido por la fuerza que tiene tres
componentes y la superficie que también tiene 3 componentes, resultan nueve
componentes, de las cuales tres son perpendiculares a la superficie que serán los
18
esfuerzos normales y seis actúan paralelos a la superficie que serán los esfuerzos de
cizalla.
Ilustración 6 Componentes del Tensor de Esfuerzos, fuente REDUCA 2010
De acuerdo con el sistema de referencia anterior (Ilustración 6), la fuerza se puede
expresar como:
𝐹𝑖 = 𝜎𝑖𝑗 𝑛𝑗 𝑆 (2,3,3)
Donde σ = σij es el tensor de esfuerzos que tiene dos índices, i relacionado con la
dirección de la fuerza y j relacionado con la orientación de la superficie.
Para una superficie con orientación cualquiera como se muestra en la ilustración 7.
Ilustración 7 Representación de Superficie de Falla, fuente D. Legrand, U. Chile
Se tiene una representación matricial de la fuerza de la siguiente manera:
(
𝐹𝑥𝐹𝑦𝐹𝑧
) = (
𝜎𝑥𝑥 𝜎𝑦𝑥 𝜎𝑧𝑥𝜎𝑥𝑦 𝜎𝑦𝑦 𝜎𝑧𝑦𝜎𝑥𝑧 𝜎𝑦𝑧 𝜎𝑧𝑧
)(
𝑛𝑥𝑛𝑦𝑛𝑧)𝑆 (2,3,4)
19
2.4. Régimen de esfuerzos
Los esfuerzos se pueden clasificar en tres tipos principales y un par de casos intermedios,
esta clasificación está en función de los ejes principales y la relación existente entre sus
magnitudes relativas. Partiendo de que las magnitudes son definidas usando una
notación estándar con esfuerzos compresivos positivos y S1>S2>S3, de manera que S1
es el máximo y S3 la componente principal mínima de compresión del esfuerzo. A
continuación se describen los tipos de régimen de esfuerzos.
FN: Falla Normal
NS: Predomina: Falla Normal,
con component strike
LS: Falla Lateral (incluye
componentes menores normales
o inversas)
IS: Falla Inversa con
componente Strike
20
TF: Falla Inversa
Ilustración 8 Régimen de Esfuerzos, Tomado de World Stress Map.
2.5. Modelos de inversión
El tensor de esfuerzos puede determinarse a partir del análisis de los siguientes diferentes
tipos de datos:
Análisis de mecanismos focales de terremotos: (Michael 1987, Rivera y
Cisternas 1990, Hardebeck 2006, Olaiz 2006, etc.). Es el único método que
permite la estimación del estado de esfuerzos en todo el rango de profundidades
de la litosfera.
Roturas de sondeos: utiliza las concentraciones naturales de esfuerzos alrededor
de los sondeos.
Medidas de esfuerzos por sobre-perforación: medidas de esfuerzos o
deformación de un volumen de roca cuando se aísla del resto del macizo.
Fracturación hidráulica: se somete a un sondeo a presión hidráulica, hasta que
se desarrolla una fractura de tensión paralela a la dirección máxima de
compresión. Solo informa de la orientación de σHmax.
Sin embargo, son los métodos de inversión de esfuerzos o análisis de población de fallas
los que mejor caracterizan el estado de esfuerzos activos en la corteza superior. Este tipo
de análisis fue inicialmente desarrollado para fallas, aunque es inmediata su aplicación a
mecanismos focales de terremotos.
Para un tensor de esfuerzos dado, el problema consiste en determinar la orientación el
sentido del movimiento, conocida la orientación de un plano de falla, por el contrario, el
problema inverso consiste en determinar el tensor de esfuerzos, a partir de las orientación
y sentidos de una población de fallas o de mecanismos focales de terremotos (Olaiz
2006). Entre las técnicas de inversión que existen, el modelo de inversión amortiguada de
Hardebeck (2006) se destaca sobre otras técnicas, a pesar que el método es similar a
21
otras, en la prueba del tablero de ajedrez, la inversión amortiguada elimina los errores de
rotación mostrados por una inversión no amortiguada, al tiempo mejora las tensiones
mucho mejor que un modelo de suavizado o una inversión de ventana móvil.
2.6. Sistemas de Información Geográfica
Los sistemas de información geográfica permiten la integración de información
cuantitativa y cualitativa asociada a una variables espaciales que se pueden distinguir en
dos tipos discretas y continuas, las primeras serán objetos puntuales, una casa por
ejemplo y las segundas serán datos como los de temperatura o precipitación, variables
continuas que se distribuyen sobre todo el espacio. Estas variables se pueden almacenar
en dos tipos de datos, uno raster (pixel) y otro vectorial (punto, línea o polígono), a estas
entidades vectoriales se les pueden asociar n cantidad de datos que describen la entidad
y serán denominados atributos.
Un sismo trae asociada una cantidad de información como es su magnitud (diferentes
tipos), profundidad, año, mes, día, hora, fuente y coordenadas, estas últimas son las que
permiten geo-referenciar espacialmente en un sistema acorde con una región
determinada, esta información es publicada por las autoridades sismológicas en catálogos
sísmicos, estos catálogos son una base de datos con información de diferentes sismos a
lo largo de un periodo de tiempo, estas características permiten distribuir espacialmente
los sismos. Siguiendo con el proceso de la información sísmica se pude seguir
adquiriendo atributos asociados a los sismos o también generar nuevas entidades con
datos asociados que pueden ser asociados al SIG y ser manipulados en diferentes capas.
Esto permite hacer un análisis posterior de la información, hacer evaluación de variables
específicas, análisis de correspondencia o que demuestren una característica
predominante, zonificaciones, muestra de únicas entidades etc. Como producto final se
pude generar mapas con información a diferentes niveles, generales y específicos, de una
o varias entidades, de una selección de datos y a diferentes escalas, según la necesidad
del usuario final de la información.
22
2.7. Software Libre
El software que se utiliza para aplicar el Modelo de inversión amortiguada es el SATSI
(Spatial and Temporal Stress Inversion) es una versión modificada del código de Michael
(JGR 1984, 1987) que permite invertir datos espaciales y/o temporales del campo de
esfuerzos de mecanismos focales. Esta inversión encuentra la solución menos compleja
de un campo de esfuerzos que es consistente con los datos, el software contiene
directorios en código C, que implementa el modelo de inversión para 2 y 4 dimensiones
de los campos de esfuerzos, una vez procesados los datos estos podrán ser visualizados
en el Software Libre GMT. Desarrollado por Jeanne Hardebeck y Andy Michael (2006).
El software GMT (Generic Mapping Tools), es una librería de herramientas desarrollada
para sistema operativo UNIX, ejecutada mediante una terminal de comandos, que permite
generar mapas a partir de ficheros con datos. El GMT permite manipular datos, hacer
operaciones matemáticas, generar mapas de contornos, sombreados de relieve, mapas
vectoriales, utilizados por ejemplo para visualizar los tensores de esfuerzos, es un
software libre potente, aunque no tiene una interfaz grafica.
3. Inversión amortiguada de esfuerzos tectónicos
3.1. Análisis de mecanismos focales de terremotos como modelo de
inversión.
Existen diferentes métodos de inversión que permiten obtener los esfuerzos tectónicos
analizando los mecanismos focales o una población de fallas, muchos de estos métodos
parten de asumir que posterior al proceso de facturación aparecen una serie de
propiedades geométricas independiente del cambio de escala, lo que permite realizar
inversiones de datos de fallas o mecanismos focales sin escalarlos según su magnitud,
Según Olaiz (2006), los primeros métodos de inversión para poblaciones de fallas fueron
propuestos por Carey y Brunier en 1974, en donde se asumía que la estría coincidía con
la dirección del máximo esfuerzo de cizalla hipótesis discutida ya que en la práctica los
resultados obtenidos a partir de mecanismos focales son altamente congruentes con una
gran variedad de situaciones tectónicas, estas metodologías son basadas en la relaciones
esfuerzo-deslizamiento descritas por Wallace (1951) y Bott (1959). Otros autores como
23
Anderson (1951), estableció una relación entre el estado esfuerzos, la simetría y la
orientación de las fracturas de las rocas, denominado Modelo de Fracturación de
Anderson. Existe otra formulación en la que se estableció una relación entre la orientación
y la forma del elipsoide de esfuerzos que actúa sobre un plano preexistente y se expresa
matemáticamente por la Ecuación de Bott. Otro método es el de los Diedros Rectos,
método grafico más utilizado y aplicable tanto a mecanismos focales como a fallas.
También está el modelo de deslizamiento que busca explicar las fracturas se disponen
según una simetría ortorrómbica con respecto a los ejes fundamentales del elipsoide de
deformación, en condiciones tri-axiales de deformación frágil natural y experimental (Olaiz
2006). Este modelo es deducido a partir del Modelo de Fracturación de Anderson y del
criterio de fractura de Navier-Coulomb.
Para el presente trabajo de Cálculo y zonificación de los tensores de esfuerzos generados
por la placa nazca en el departamento del Valle del Cauca, Colombia, se utiliza el método
de inversión de Hardebeck (2006), quien plantea un modelo de inversión amortiguado a
escalas regionales, este modelo es derivado de Michael (1987), quien plantea un modelo
de inversión para los datos de deslizamiento de una falla, el ángulo de dirección de
deslizamiento (θ), ángulo deslizamiento del bloque con respecto a la superficie (δ) y el
ángulo de deslizamiento sobre el plano de falla (λ) (Ilustración 9), primero se plantea la
suposición de que la magnitud de la tracción tangencial es similar en los diferentes planos
de falla en el momento de ruptura, segundo las magnitudes relativas del tensor de
esfuerzos isotrópico y desviador se asume que son por fractura simple, este es el primer
intento de aplicar este tipo de inversión a datos de deslizamientos obtenidos de procesos
dúctiles como lo son las fallas. Esta nueva inversión también modifica los ejes del tensor
de esfuerzos, limitando a que un eje principal sea vertical y los otros dos ejes principales
sean horizontales, diferentes estudios apoyan la tesis de que uno de los ejes principales
del tensor está cerca de la vertical, esta restricción podrá ser útil cuando se analicen datos
dispersos debido a que esto reduce el número de variables en la inversión.
24
Ilustración 9 Diagrama de Variables en el Bloque de Masa. Héctor William Arboleda
El objetivo fue encontrar un solo tensor de esfuerzos uniforme que tenga mayor
probabilidad de generar la falla. El tensor de esfuerzos se puede separar en dos
componentes, una componente isotrópica y una componente desviadora, restringir el
tensor de esfuerzos a la componente isotrópica no puede ser determinado solo con los
datos de deslizamiento de la falla, el tensor de esfuerzos isotrópico está compuesto por
las componentes principales medias del tensor de esfuerzos principales, para este modelo
se realiza la inversión para el tensor de esfuerzo desviador (Deviatoric), el tensor de
esfuerzo desviador es el que causa la rotación o distorsión de la masa, siendo siempre
máximo positivo el que representa la compresión y siempre mínimo negativo el que
representa la extensión, así mismo solo se pueden determinar las magnitudes relativas
del tensor de esfuerzos y no las absolutas. Para hallar la componente desviadora del
tensor de esfuerzos solo se considera la orientación de deslizamiento respecto a la
orientación de la falla.
La base de esta inversión es que la dirección de la tracción tangencial sobre el plano de
falla tiende a ser paralela a la dirección del movimiento del deslizamiento (�̂�).
�̂� = �⃗⃗� (�̂�,𝜎)
|�⃗� (�̂�,𝜎)|= �̂� (3,1);
X3
X2
X1
Norte
λ
β
θ
δ
�̂�
�̂�
�̂�
25
Donde, �⃗⃗� (�̂�, 𝜎) es la tracción tangencial sobre el plano de falla con normal unitario �̂�,
debido al tensor de esfuerzos desviador σ; y �̂� es el vector unitario asociado con �⃗⃗� .
Con la información de una sola falla no es posible restringir el tensor de esfuerzos al
tensor desviador σ, es necesario resolver (3,1) para varias fallas al mismo tiempo y así
poder determinar un solo tensor de esfuerzos σ que mejor satisfaga todas las fallas. Esto
es lo mismo que asumir que el tensor de esfuerzos es el mismo para toda una región.
Adicionalmente es necesario asumir la naturaleza de �⃗⃗� (�̂�, 𝜎) como la tracción total menos
la tracción normal.
�⃗� = 𝜎�̂� − [(𝜎�̂�) ∙ �̂�]�̂� (3,2);
Nótese que �⃗� es lineal con respecto a σ, implicando que |�⃗� | no es lineal respecto a σ, de
esta manera |�⃗� | si puede ser removida de (3,1) ya que una inversión lineal es posible para
σ. Las inversiones lineales son preferibles por tres razones: computacionalmente son más
eficientes, no requieren una estimación inicial y produce fácilmente los límites de
confianza (Michel 1984). Para remover |�⃗� | debe ser determinada. Además de la
dirección de deslizamiento, otro dato disponible para cada falla es el deslizamiento, esto
implica que la magnitud de la tracción tangencial para cada plano de falla es significante.
Si fuese pequeña sobre un plano, se podría determinar que la deformación tomara
preferiblemente un lugar arriba, lo más fácil es asumir que |�⃗� | es constante a lo largo de
todo el plano de falla, como solo se pueden hallar las magnitudes relativas del tensor de
esfuerzo, se puede suponer que |�⃗� | = 1. Para este proceso de inversión se requiere que
las magnitudes de tracción tangenciales sean similares no precisamente 1, así que |�⃗� |
para cada plano de falla se ajusta minimizando la diferencia |�⃗� | − 1.
Mientras |�⃗� | es la función de orientación del plano de falla, se puede creer que un
conjunto de fallas que se deslizan en eventos relacionados tendrán tracciones similares.
Tomando |�⃗� | = 1 en combinación con (3,1) resultado de �⃗� = �̂�. Se toma en combinación
con la restricción de que el tensor de esfuerzos isotrópico es cero (𝜎33 = −(𝜎11 + 𝜎22)), y
la ecuación (4,2), se plantea el siguiente ecuación de inversión lineal.
26
(
𝑛1 − 𝑛13 + 𝑛1𝑛3
2 𝑛2 − 2𝑛2𝑛12 𝑛3 − 2𝑛3𝑛1
2 −𝑛1𝑛22 + 𝑛1𝑛3
2 −2𝑛1𝑛2𝑛3−𝑛2𝑛1
2 + 𝑛2𝑛32 𝑛1 − 2𝑛1𝑛3
2 −2𝑛1𝑛2𝑛3 𝑛2 − 𝑛23 + 𝑛2𝑛3
2 𝑛3 − 2𝑛3𝑛22
𝑛1 − 𝑛13 + 𝑛1𝑛3
2 −2𝑛1𝑛2𝑛3 𝑛1 − 2𝑛1𝑛32 −𝑛2
2𝑛3 − 𝑛3 + 𝑛33 𝑛2 − 2𝑛2𝑛3
2
)
(
𝜎11𝜎12𝜎13𝜎22𝜎23)
=
(𝑆1𝑆2𝑆3
) (3,3)
La ecuación (3,3) contiene a la ecuación (3,2), así que obliga a que la tracción tangencial
prevista sea en el plano de falla. De forma simplificada la ecuación (3,3) se puede ver
como:
𝐺𝑚 = 𝑑 (3,4)
Donde, G es la matriz derivada del vector normal al plano de falla, m es el vector de las
componentes del tensor de esfuerzos y d es el vector de dirección de deslizamiento de la
falla. La Ecuación (3,4) se expande para cada uno de los datos de fallas de la siguiente
manera:
(
𝐺1𝐺2𝐺3⋮𝐺𝑘)
𝑚 =
(
𝑆1𝑆2𝑆3⋮𝑆𝑘)
(3,5)
Finalmente la ecuación (3,5) puede ser resuelta para los elementos del tensor de
esfuerzos y los errores estimados por mínimos cuadrados. Una vez conocido el tensor de
esfuerzos, se pueden calcular los esfuerzos principales y sus direcciones.
Michel estableció un modelo de inversión determinado a partir de la geometría de la falla y
su información asociada al deslizamiento, asumiendo que la tracción tangencial al plano
de falla es similar en los diferentes planos de falla, y en el que se restringe la inversión al
tensor de esfuerzos desviador. Esto permite hacer una relación entre la tracción
tangencial, el vector del movimiento del deslizamiento y la normal del plano de falla,
hallando el modelo de inversión simplificado que más se acerca a un tensor de esfuerzos
que sea capaz de general la falla. Logrando el objetivo de tener un modelo de inversión
27
que minimiza la diferencia entre el vector tangencial al plano de falla y el vector de
deslizamiento.
Hardebeck (2006) partiendo del modelo de inversión de Michel (1984), plantea un modelo
de inversión amortiguado para un conjunto de datos de tensores de esfuerzos aplicado a
una escala regional, esta metodología permite obtener la orientación de los tensores de
esfuerzos con la mínima complejidad necesaria para ajustar los datos. El modelo menos
complejo es a veces encontrado mediante una inversión amortiguada es el que minimiza
la suma ponderada de dos valores, el primero es el desajuste de los datos y el segundo la
longitud del modelo. Esto permite establecer el grado de heterogeneidad en la solución,
la cual depende directamente de los datos y la heterogeneidad en las direcciones de los
esfuerzos. En áreas donde la heterogeneidad es condicional para los datos, las
inconsistencias de datos inadaptados para un modelo suavizado son mayores que las
inconsistencias de longitud del modelo para un modelo heterogéneo. En áreas donde la
heterogeneidad no es condicional para los datos, las inconsistencias de longitud del
modelo para un modelo heterogéneo son mayor que las inconsistencias de datos
inadaptados para un modelo suavizado, mientras que un modelo de inversión
amortiguado, es un modelo de suavizado adaptativo que es más sofisticado que hacer un
simple modelo de suavizado, ya que se discriminan entre las variaciones condicionadas
por los datos y que deben mantenerse en el modelos y las variaciones que no son
indispensables y pueden ser removidas. El parámetro de amortiguamiento es
seleccionado de la curva de compensación generada entre la longitud del modelo
innecesario y el desajuste excesivo de los datos, seleccionando el parámetro que genera
un equilibrio entre ellos.
Así como Michel (1984), se planteó establecer un modelo de inversión que determine un
solo tensor de esfuerzos uniforme que tenga mayor probabilidad de generar la falla,
Hardebeck (2006), tiene como objetivo producir un modelo de campo de esfuerzos menos
complejo para ajustar los de datos de las observaciones, para lograrlo se introduce un
método de inversión amortiguado. En este proceso se divide la región de la zona de
estudio en pequeñas áreas y se hace la inversión simultáneamente lo cual minimiza la
diferencia de los esfuerzos entre las áreas contiguas. Con un adecuada amortiguación,
se obtiene un modelo del campo de esfuerzos a nivel regional que solo contiene
variaciones en el campo de esfuerzos que son requeridas por los datos.
28
Este nuevo modelo inversión amortiguada se amplió a realizar el proceso para varios
datos de mecanismos focales simultáneamente, partiendo de la ecuación (3.4). Este
modelo de inversión fue expandido a una escala regional, acá se halla el tensor de
esfuerzos sobre un punto I x J de una grilla de puntos de dos dimensiones, para este
proyecto se hallaron los esfuerzos en dos dimensiones, pero este proceso puede ser
expandido fácilmente a tres o cuatro dimensiones incluyendo la profundidad y/o el tiempo.
El modelo inicial es modificado al obtener un tensor de esfuerzos para el punto en la fila i-
esima y la j-esima columna, esto es representado por un vector, en donde cada
componente del vector corresponde a la solución en un punto de la grilla, esta expansión
se ve representada de la siguiente manera:
(
𝐺11𝐺12 0
⋯𝐺21
𝐺220 ⋯
𝐺𝐼𝐽)
(
𝑚11𝑚12⋯𝑚21𝑚22…𝑚𝐼𝐽)
=
(
𝑑11𝑑12…𝑑21𝑑22…𝑑𝐼𝐽)
A la solución de este modelo por mínimos cuadrados se añade una restricción de planitud,
con lo que se minimizara la diferencia de los tensores de esfuerzos entre los puntos
adyacentes, sin incluir los puntos adyacentes en diagonal, solo los cuatro puntos laterales,
esto se hace mediante una matriz de amortiguación D, dada por la siguiente expresión:
𝐷 =
(
𝐼 −𝐼 ⋯ 0 0 ⋯𝐼 0 ⋯ −𝐼 0 ⋯0 𝑖 ⋯ 0 −𝐼 ⋯
⋯0 0 ⋯ 𝐼 −𝐼 ⋯
⋯ )
, (3,7)
Donde I es la matriz identidad de 5X5, la matriz D tiene M = 5 * I * J columnas y 5 filas por
cada par de puntos adyacentes en la grilla de puntos. Cada fila de la matriz tiene dos
entradas diferentes de cero, una +I y una -I para un mismo tensor de esfuerzos para un
par de puntos adyacentes. Esta inversión amortiguada disminuye la suma de los valores
del desajuste de los datos y la longitud del modelo. La longitud del modelo es
representada por la norma L2 del vector que contiene las diferencias entre cada
29
componente del tensor de esfuerzos para cada par adyacente de puntos de la grilla, dada
por:
‖𝐷𝑚𝑡𝑜𝑡𝑎𝑙‖2 (3.8)
Y la predicción del desajuste de los datos del modelo:
‖𝐺𝑡𝑜𝑡𝑎𝑙𝑚𝑡𝑜𝑡𝑎𝑙 − 𝑑𝑡𝑜𝑡𝑎𝑙‖2 (3.9)
La solución de mínimos cuadrados para el modelo amortiguado será la siguiente:
{𝐺𝑡𝑜𝑡𝑎𝑙𝑇 𝐺𝑡𝑜𝑡𝑎𝑙 + 𝑒
2𝐷𝑇𝐷}𝑚𝑡𝑜𝑡𝑎𝑙 = 𝐺𝑡𝑜𝑡𝑎𝑙𝑇 𝑑𝑡𝑜𝑡𝑎𝑙 (3.10)
Donde e es un parámetro de amortiguación escalar que controla el peso relativo
minimizando el desajuste de los datos y la longitud del modelo. Este parámetro se
obtiene de la curva de compensación entre el desajuste de los datos y la longitud del
modelo, seleccionando el valor que minimiza conjuntamente el desajuste y la longitud del
modelo. Como Gtotal y D son matrices dispersas, la anterior ecuación (3.10) se resuelve
usando una implementación del método del gradiente conjugado.
4. Desarrollo del proyecto
4.1. Descripción de la zona de estudio
La zona de estudio comprende el departamento del Valle del Cauca y zonas aledañas. El
Departamento del Valle del Cauca está situado al suroccidente del país, formando parte
de las regiones andina y pacífica. Cuenta con una superficie de 22.140 km2. Las zonas
aledañas son porciones de departamento con los que limita de la siguiente manera: por el
Norte con los departamentos de Chocó, Caldas y Quindío; por el Este con los
departamentos del Quindío y Tolima, por el Sur con el departamento del Cauca y por el
Oeste con el océano Pacífico y el departamento del Chocó.
30
Ilustración 10 Zona de Estudio, Fuente: OSSO 2013
Esta zona está bajo la influencia directa de la zona de subducción de la placa de Nazca al
limitar con la placa de Sur América, a esta falla inversa están asociadas las geo-formas
que se presentan en la zona de estudio y hace que tenga una actividad sísmica
permanente. La zona de estudio (Valle del Cauca y Zonas aledañas) se limita por los
paralelos 2.9º y 5.3º y los meridianos -75.6º a -79.5º.
Las principales direcciones que siguen las fallas sobre el Valle del Cauca son N-S y E-W,
entre las que se encuentran las siguientes: falla Garrapatas, falla Dagua-Calima, falla El
Tambor, Falla Naya-Micay, Falla Toro, Falla Argelia, Falla Paraíso y el Sistema de Fallas
Romeral, adicionalmente la zona de estudio se encuentra bajo la influencia de las Fallas
Armenia, Montenegro, Ibagué, Cucuana, Piendamo y La Dina, y sobre el pacifico se
encuentra la Zona de Subducción de la placa de Nazca.
31
4.2. Distribución de mecanismos focales y la geología estructural.
La zona de estudio se encuentra bajo constante influencia de la zona de subducción de la
placa nazca, adicionalmente cuenta con una serie de fallas en diferentes direcciones, los
movimientos de una falla genera un sismo, esta ruptura se puede representar mediante
un mecanismo focal, de manera que habrá una estrecha relación entre el mecanismo
focal y dirección y tipo de movimiento de la falla que se pueda encontrar cerca de él.
Hacer una relación entre los mecanismos focales y las fallas de la zona de estudio
permitirá realizar un análisis de la dinámica de la geología estructural de la zona de
estudio.
4.3. Geología estructural de la zona de estudio.
La zona de estudio comprendida por el departamento del Valle del Cauca y sus
alrededores se encuentra principalmente bajo la influencia de la zona de subducción de la
placa de Nazca con la placa continental, situación a la que se debe que sea una zona con
una alta actividad sísmica y generación de las principales fallas que se describen a
continuación.
4.3.1. Zona Subducción
Ubicada aproximadamente a 100 Km al oeste de la costa pacífica de Colombia, la falla se
extiende de norte a sur sobre la costa pacífica Colombiana. La geometría de la zona de
subducción se define de norte a sur por un conjunto de hipocentros de terremotos más
densos entre los 5° Latitud Norte y el límite entre Colombia – Panamá. El límite Norte se
encuentra a 7° Latitud Norte a lo largo de una línea imaginaria este-oeste que se proyecta
a través de Darién y el nido de Bucaramanga. Las fallas de tipo normal y de empuje se
presentan de manera perpendicular a los planos nodales a lo largo de la fosa de
subducción. La zona de subducción se divide en tres secciones, la primera costa afuera
de Colombia, segunda la base de la inmersión y tercera la dirección de convergencia de
las placas.
32
Tabla 4.1 Descripción de la Zona de Subducción
Longitud Km
Rumbo Inclinación Sentido Movimiento
Tasa de Movimiento
Ultimo Movimiento
Compilador Afiliación Fecha
670,2 N12,4°E+- 31° Variable, bajo al Este
Empuje bajo al Este
70 mm/año Histórico 1906, 1942, 1958, 1979
Gabriel Paris
Univalle 06/2000
Ilustración 11 Mapa tectónico y fallas de la zona de estudio. Trenkam et.al.
4.3.2. Naya-Micay
Ubicada al suroeste y paralela a la zona costera del pacifico colombiano sobre el Valle del
Cauca. Esta falla desplaza rocas sedimentarias marinas y no marinas del Plioceno.
Localmente se compensa con depósitos aluviales Cuaternarios. En general no se elevan
sedimentos terciarios al este de la falla y de cuaternario al oeste de la falla. La falla
parece ser una continuación al norte de la falla Remolino-El Charco.
33
Tabla 4.2 Descripción Falla Naya-Micay
Longitud Km
Rumbo Inclinación Sentido Movimiento
Tasa de Movimiento
Ultimo Movimiento
Compilador Afiliación Fecha
670,2 N34,1°E±12° Este
Dextral (Lateral Derecha) Oblicua
0,2-1,0 mm/año
Cuaternario <1,6
Gabriel Paris
Ingeominas 09/1993
4.3.3. Falla El Tambor
La Falla bordea la vertiente occidental de la Cordillera Occidental. Gran parte de su
recorrido contiene rocas volcánicas y sedimentarias del Cretácico sobre el este y al
costado oeste del Neógeno.
Tabla 4.3 Descripción Falla El Tambor
Longitud Km
Rumbo Inclinación Sentido Movimiento
Tasa de Movimiento
Ultimo Movimiento
Compilador Afiliación Fecha
117,5 N26,1°E±9° Este Alto
Inversa, Dextral? (Lateral Derecha
0,2 mm/año
Holoceno / post glacial <1,5 Ka
Gabriel Paris
Ingeominas 09/1993
4.3.4. Sistema de Falla Dagua Calima
Esta falla se forma sobre roca oceánica del Cretáceo en la parte media de la vertiente
oriental de la Cordillera Occidental, se localiza al suroccidente de Colombia y al nororiente
de Cali. La falla coloca roca volcánica oceánica contra roca sedimentaria, causando
alrededor de 200m de desplazamiento de una superficie terciaria de erosión o penillanura
a través de la falla. Gran parte del movimiento se originó durante la Orogenia Andina. A
nivel local, hacia la mitad sur de la falla, el escarpe orientado al este se asocia con una
graben de 3 – 5 Km. De ancho por unos 20 Km. De Lardo que está parcialmente lleno de
depósitos fluviales del Cuaternario.
34
Tabla 4.4 Descripción Sistema de Falla Dagua Calima
Longitud Km
Rumbo Inclinación Sentido Movimiento
Tasa de Movimiento
Ultimo Movimiento
Compilador Afiliación Fecha
214,9 N18,8°E±20° 70°E–SE Normal < 0,2 mm/año
Cuaternario <1,6Ma
Gabriel Paris
Ingeominas 09/1993
4.3.5. Falla Garrapatas
La falla tiene un valle muy bien desarrollado en forma de V, a lo largo de las partes
superiores de los ríos Garrapatas y Las Vueltas. La falla se encuentra al occidente de la
Ciudad de Buga y se extiende entre el eje de la Cordillera Occidental y el Cerro Paraguas
y se desplaza en rocas volcánicas y sedimentarias oceánicas. Probablemente se conecta
con la Falla Argelia.
Tabla 4.5 Descripción Falla Garrapatas
Longitud Km
Rumbo Inclinación Sentido Movimiento
Tasa de Movimiento
Ultimo Movimiento
Compilador Afiliación Fecha
140,7 N60,8°E±14° 50°NW Inverso / Probable Dextral
< 0,2 mm/año
Cuaternario <1,6Ma
Gabriel Paris
Ingeominas 09/1993
4.3.6. Argelia
Ubicada en la Cordillera Occidental, al oriente de la ciudad de Pereira. La falla corta
rocas oceánicas del Mesozoico que fueron unidas al continente durante el Cretácico
Tardío. La falla Argelia muestra algunas características neotectónicas y probablemente
se conecta con la falla de Garrapatas.
Tabla 4.6 Descripción Falla Argelia
Longitud Km
Rumbo Inclinación Sentido Movimiento
Tasa de Movimiento
Ultimo Movimiento
Compilador Afiliación Fecha
59,5 N14,5°E±18° Moderada alta al W
Inversa / Probable Dextral
< 0,2 mm/año
Cuaternario <1,6Ma
Gabriel Paris
Estudios Geológico Mineros
06/1998
35
4.3.7. Falla Toro
La falla Toro corta rocas oceánicas unidas a la Cordillera Occidental, cerca del Valle del
rio Cauca. Es una de las fallas del sistema de Falla Cauca-Patía que limita con el lado
oriental de la cordillera Occidental a lo largo de la mayor parte de la longitud.
Tabla 4.7 Descripción Falla Toro
Longitud Km
Rumbo Inclinación Sentido Movimiento
Tasa de Movimiento
Ultimo Movimiento
Compilador Afiliación Fecha
61,8 N6,6°E±8° Desconocida Sinestral lateral Izquierdo
< 0,2 mm/año
Cuaternario <1,6Ma
Gabriel Paris
Estudios Geológico Mineros
06/1998
4.3.8. Falla Ibagué
La falla Ibagué atraviesa la parte central y vertiente oriental de la Cordillera Central de
Colombia, cerca de la ciudad de Ibagué. La falla corta rocas Precámbricas y la edad
Jurásica Ibagué Batolito. La falla golpea WSW-ENE controlando el rio Cócora,
desplazando y deformando el Neógeno abanico de Ibagué, que es de origen sedimentario
volcánico.
Tabla 4.8 Descripción Falla Ibagué
Longitud Km
Rumbo Inclinación Sentido Movimiento
Tasa de Movimiento
Ultimo Movimiento
Compilador Afiliación Fecha
123,9 N67,9°E±11° Vertical Común
Dextral Oblicua ligeramente
1 – 5 mm/año
Holoceno / post glacial <15 Ka
Gabriel Paris
Ingeominas 09/1993
4.3.9. Falla Palestina
La falla Palestina se extiende desde el norte del Departamento de Antioquia a la zona sur
del Volcán Nevado del Ruiz. Se extiende principalmente a lo largo de la vertiente oriental
de la Cordillera Central de Colombia, desplazando rocas cristalinas metamórficas del
Paleozoico y en menor cantidad rocas plutónicas del Mesozoico. Las rocas más antiguas
se encuentran principalmente en el bloque occidental, que eleva una superficie de erosión
36
probablemente Mioceno, cuyos erosionados remanentes son caracterizados por planos
alineados de bancos y espuelas estrechos.
Tabla 4.9 Descripción Falla Palestina
Longitud Km
Rumbo Inclinación Sentido
Movimiento Tasa de
Movimiento Ultimo
Movimiento Compilador Afiliación Fecha
369,6 N17,8°E±11° W Alto Lateral
Izquierda Inversa
0,2 – 1 mm/año
Cuaternario <1,6 Ma
Gabriel Paris
Ingeominas 09/1993
4.3.10. Falla Cucuana
La falla Cucuana es paralela y al sur de la falla de Ibagué en la Cordillera Central,
atravesando rocas metamórficas del Paleozoico, rocas ígneas del Mesozoico y camas
terciarias. La falla se extiende en la Cordillera Oriental a través del valle del Magdalena
cortando camas del Cretácico. El relleno sedimentario y volcánico del Cuaternario del
valle del Rio Magdalena no se ha reportado deformaciones por la falla.
Tabla 4.10 Descripción Falla Cucuana
Longitud Km
Rumbo Inclinación Sentido
Movimiento Tasa de
Movimiento Ultimo
Movimiento Compilador Afiliación Fecha
369,6 N17,8°E±11° W Alto Lateral
Izquierda Inversa
0,2 – 1 mm/año
Cuaternario <1,6 Ma
Gabriel Paris
Estudios Geológico Mineros
06/1998
4.3.11. Sistema de Falla Romeral
Este sistema de fallas de casi 700 Km. de longitud, se compone de tres o cuatro fracturas
regionales paralelas que forman una zona de transición entre las rocas oceánicas hacia el
occidente y rocas continentales hacia el oriente. La geología del dominio occidental
consiste en un cinturón ofiolita con gabbroic oceánico, rocas basálticas y sedimentarias
del Cretáceo. El dominio oriental se compone principalmente de pizarras metamórficas
continentales, oceánicas y roca continental, principalmente del Paleozoico. El sistema de
Falla Romeral forma un cinturón de deformación de 20 km a 40 km de ancho paralelo a la
vertiente occidental de la Cordillera Central de Colombia; se extiende desde el sur del
golfo de Guayaquil en Ecuador hasta el norte del mar caribe en Colombia.
37
Tabla 4.11 Descripción Sistema de Falla Romeral
Longitud Km
Rumbo Inclinación Sentido
Movimiento Tasa de
Movimiento Ultimo
Movimiento Compilador Afiliación Fecha
697,4 N17,6°E±16° 65°E Lateral
Izquierda
Gabriel Paris
Estudios Geológico Mineros
06/1998
4.3.12. Falla Montenegro
La falla Montenegro hace parte del sistema de Falla Romeral, corriendo a través de la
vertiente occidental de la Cordillera Central. La falla Montenegro se encuentra al
accidente de la ciudad de Armenia. La falla corta transversalmente y deforma los
depósitos volcánicos del Pleistoceno y los depósitos volcánicos sedimentarios del abanico
del Quindío que cubre cerca de 400 Km2.
Tabla 4.12 Descripción Falla Montenegro
Longitud Km
Rumbo Inclinación Sentido
Movimiento Tasa de
Movimiento Ultimo
Movimiento Compilador Afiliación Fecha
21,7 N25,1°E±9° W alto Lateral
Izquierda < 0,2
mm/año Cuaternario?
>1,6 Ma Gabriel
Paris
Estudios Geológico Mineros
06/1998
4.3.13. Falla Armenia
La falla Armenia es parte del sistema de fallas de Romeral en la vertiente occidental de la
Cordillera Central de Colombia. La falla cruza la ciudad de Armenia y desplaza depósitos
volcánicos del Plioceno-Pleistoceno y depósitos sedimentarios volcánicos del abanico del
Quindío, que abarca aproximadamente 400 Km2. Las características geométricas y
neotectónicas de las fallas Montenegro y Armenia son muy similares.
Tabla 4.13 Descripción Falla Armenia
Longitud Km
Rumbo Inclinación Sentido
Movimiento Tasa de
Movimiento Ultimo
Movimiento Compilador Afiliación Fecha
33 N23,2°E±11° W Muy
alto Lateral
Izquierda < 0,2
mm/año Cuaternario?
>1,6 Ma Gabriel
Paris
Estudios Geológico Mineros
06/1998
38
4.3.14. Falla Paraíso
La falla Paraíso hace parte del sistema de fallas de Romeral en el suroccidente de
Colombia. Se encuentra en la vertiente occidental de la Cordillera Central, al oriente de la
ciudad de Palmira. La falla desplaza abanicos aluviales y flujos de escombros en la
frontera oriental del Valle del Cauca. Al norte del rio Amaime, la falla parece ser más
activo en Cuaternario que en la porción sur del rio. Dos trincheras fueron abiertas en la
parte norte: la Trinchera Venecia y la Trinchera Piedechinche, cada una muestra
deformación tectónica del Holoceno.
Tabla 4.14 Descripción Falla Paraíso
Longitud Km
Rumbo Inclinación Sentido
Movimiento Tasa de
Movimiento Ultimo
Movimiento Compilador Afiliación Fecha
33 N12,5°E±3° 12°NE –
30°E Lateral
Derecho 0,2 -1
mm/año
Holoceno post glacial
<15 Ka
Gabriel Paris
Univalle 08/1999
4.3.15. Falla Piendamo
La falla Piendamo es parte del sistema de falla Romeral en el suroccidente de Colombia;
que se encuentra en la base de la parte delantera de la montaña de la vertiente occidental
de la Cordillera Central, al norte de la ciudad de Popayán. La falla desplaza depósitos
volcánicos piroclásticos y flujos de lodo de la formación Terciaria – Cuaternaria de
Popayán. Forma un bloque topográfico sobresaliente y con una tectónica limitada por dos
fallas que se cruzan.
Tabla 4.15 Descripción Falla Piendamo
Longitud Km
Rumbo Inclinación Sentido
Movimiento Tasa de
Movimiento Ultimo
Movimiento Compilador Afiliación Fecha
28,3 N8,4°W±18° Medio a
Alto al E y NO
Inversa Lateral
Derecha
0,2 - 1 mm/año
Cuaternario? <1,6 Ma
Gabriel Paris
Estudios Geológico Mineros
06/1998
39
4.3.16. Falla Córdoba – Navarco
Esta sección del sistema de fallas de Romeral está formado por las fallas de Córdova y
Navarco, que son hebras orientales del sistema de fallas Romeral, al sur de la ciudad de
Armenia. Estas fallas se encuentran dentro del epicentro del terremoto de Armenia del 25
de enero de 1999. Las fallas se extienden a través de rocas esquilada cataclástica,
basáltica no deformada y oceánicas sedimentarias del Cretácico, que afloran en la
vertiente oriental de la Cordillera Central de Colombia.
Tabla 4.16 Descripción Falla Córdova - Navarco
Longitud Km
Rumbo Inclinación Sentido
Movimiento Tasa de
Movimiento Ultimo
Movimiento Compilador Afiliación Fecha
21,2 N18,5°W±4° Vertical Inversa Lateral
Derecha
< 0,2 mm/año
Histórico 1999
Gabriel Paris
Univalle 06/1999
4.4. Mecanismos focales.
Para el desarrollo de este estudio se utilizaron 43 , del Catálogo Global Centroid Moment
Tensor (CMT) de la Universidad de Harvard es una recopilación de los mecanismos
focales de terremotos desde 1976 hasta la actualidad y permite realizar una búsqueda
rápida en la página web: http://www.globalcmt.org/CMTsearch.html.
Se tomaron los mecanismos focales con magnitudes entre 0 y 10 desde el 1 de enero
1976 hasta el 28 de febrero de 2018 entre los paralelos 2.9º y 5.3º y los meridianos -75.6º
a -79.5º. (Ilustración 11), que comprende el Departamento del Valle del Cauca. Este
catálogo
Los campos de la base de datos son:
Región: Zona geográfica
Profundidad (km)
Fecha
Coordenadas geográficas: Latitud Longitud
40
Magnitud: Kanamori (energía)
Se dan las soluciones con el plano de falla y el plano auxiliar
Strike: Azimut de una línea creada por la intersección de un plano de falla y una
superficie horizontal,
Dip: Buzamiento (grados) Angulo entre la falla y un plano horizontal
Slip o Rake: Dirección (grados)del movimiento del bloque colgante de la falla durante la
ruptura y es medido sobre el plano de la falla
Ilustración 12 Consulta de Datos Catalogo CMT
41
Tabla 4.17 Mecnismos Focales desde 1976 hasta 2018
Mecanismo 1 Región SUR PANAMA Fecha 1985/06/10
060185A
Latitud Longitud Profundidad Magnitud
3,24 -78,99 26 5,5
Solución Strike Dip Slip
Plano de Falla 32 19 125
Plano Auxiliar 176 74 79
Mecanismo 2 Región SUR PANAMA Fecha 1987/01/25
012587A
Latitud Longitud Profundidad Magnitud
3,30 -79,31 15,4 5,7
Solución Strike Dip Slip
Plano de Falla 19 32 -109
Plano Auxiliar 221 60 -79
Mecanismo 3 Región COSTA OESTE COLOMBIA
Fecha 1988/09/20
092088A
Latitud Longitud Profundidad Magnitud
4,99 -77,68 24,1 5,8
Solución Strike Dip Slip
Plano de Falla 14 18 104
Plano Auxiliar 179 73 85
Mecanismo 4 Región COLOMBIA Fecha 1988/11/29
112989B
Latitud Longitud Profundidad Magnitud
4,72 -77,07 79.8 5,7
Solución Strike Dip Slip
Plano de Falla 220 51 3
Plano Auxiliar 128 88 141
42
Mecanismo 5 Región COSTA OESTE COLOMBIA
Fecha 1991/11/19
111991G
Latitud Longitud Profundidad Magnitud
4,80 -77,18 19,1 7,2
Solución Strike Dip Slip
Plano de Falla 13 13 95
Plano Auxiliar 188 77 89
Mecanismo 6 Región COSTA OESTE COLOMBIA
Fecha 1991/12/10
121091D
Latitud Longitud Profundidad Magnitud
4,74 -77,48 20,7 5,2
Solución Strike Dip Slip
Plano de Falla 45 22 127
Plano Auxiliar 186 72 76
Mecanismo 7 Región COLOMBIA Fecha 1992/08/15
081592C
Latitud Longitud Profundidad Magnitud
5,08 -75,88 125,9 5,9
Solución Strike Dip Slip
Plano de Falla 228 22 -71
Plano Auxiliar 28 69 -97
Mecanismo 8 Región COLOMBIA Fecha 193/10/12
101293ª
Latitud Longitud Profundidad Magnitud
4,15 -76,98 90,4 5,2
Solución Strike Dip Slip
Plano de Falla 234 44 77
Plano Auxiliar 73 47 103
Mecanismo 9 Región SUR DE PANAMA Fecha 1994/06/03
060394C
Latitud Longitud Profundidad Magnitud
3,21 -78,78 15 5,8
Solución Strike Dip Slip
Plano de Falla 10 42 -123
Plano Auxiliar 231 56 -63
43
Mecanismo 10 Región COLOMBIA Fecha 1994/06/06
060694J
Latitud Longitud Profundidad Magnitud
2,93 -75,94 15 6,8
Solución Strike Dip Slip
Plano de Falla 206 76 170
Plano Auxiliar 299 80 14
Mecanismo 11 Región COLOMBIA Fecha 1995/02/08
020895E
Latitud Longitud Profundidad Magnitud
4,09 -76,36 68,2 6,3
Solución Strike Dip Slip
Plano de Falla 210 21 -90
Plano Auxiliar 30 69 -90
Mecanismo 12 Región COSTA OESTE COLOMBIA
Fecha 1995/04/12
041295B
Latitud Longitud Profundidad Magnitud
3,02 -78,11 54,1 5,0
Solución Strike Dip Slip
Plano de Falla 193 50 51
Plano Auxiliar 65 53 127
Mecanismo 13 Región COLOMBIA Fecha 1995/08/19
081995C
Latitud Longitud Profundidad Magnitud
5,22 -75,69 128,7 6,5
Solución Strike Dip Slip
Plano de Falla 204 30 -99
Plano Auxiliar 34 61 -85
Mecanismo 14 Región COLOMBIA Fecha 1996/09/11
091196B
Latitud Longitud Profundidad Magnitud
4,59 -76,90 118,2 5,3
Solución Strike Dip Slip
Plano de Falla 306 51 144
Plano Auxiliar 61 63 45
44
Mecanismo 15 Región COLOMBIA Fecha 1997/02/19
021997B
Latitud Longitud Profundidad Magnitud
4,78 -76,50 118,1 5,8
Solución Strike Dip Slip
Plano de Falla 138 42 105
Plano Auxiliar 299 50 77
Mecanismo 16 Región COLOMBIA Fecha 1997/12/11
121197C
Latitud Longitud Profundidad Magnitud
4,11 -75,84 189,5 6,3
Solución Strike Dip Slip
Plano de Falla 249 53 31
Plano Auxiliar 140 66 139
Mecanismo 17 Región COLOMBIA Fecha 1999/01/25
012599F
Latitud Longitud Profundidad Magnitud
4,58 -75,75 27,7 6,1
Solución Strike Dip Slip
Plano de Falla 8 65 -21
Plano Auxiliar 107 71 -153
Mecanismo 18 Región COLOMBIA Fecha 1999/01/25
012599G
Latitud Longitud Profundidad Magnitud
4,31 -75,74 33,8 5,5
Solución Strike Dip Slip
Plano de Falla 17 67 -23
Plano Auxiliar 116 69 -155
Mecanismo 19 Región COLOMBIA Fecha 2001/09/22
092201A
Latitud Longitud Profundidad Magnitud
4,03 -76,17 178,0 5,9
Solución Strike Dip Slip
Plano de Falla 247 53 31
Plano Auxiliar 137 65 139
45
Mecanismo 20 Región COSTA OESTE COLOMBIA
Fecha 2002/08/08
080802D
Latitud Longitud Profundidad Magnitud
5,17 -77,62 22,3 5,7
Solución Strike Dip Slip
Plano de Falla 12 45 -140
Plano Auxiliar 251 62 -52
Mecanismo 21 Región SUR PANAMA Fecha 2002/12/21
122102A
Latitud Longitud Profundidad Magnitud
3,76 -78,90 15 5,2
Solución Strike Dip Slip
Plano de Falla 36 40 -70
Plano Auxiliar 191 53 -106
Mecanismo 22 Región COSTA OESTE COLOMBIA
Fecha 2003/08/22
082203C
Latitud Longitud Profundidad Magnitud
3,07 -77,73 31,4 5,0
Solución Strike Dip Slip
Plano de Falla 253 45 -73
Plano Auxiliar 50 47 -106
Mecanismo 23 Región COSTA OESTE COLOMBIA
Fecha 2003/11/05
110503B
Latitud Longitud Profundidad Magnitud
5,14 -77,81 27,6 5,9
Solución Strike Dip Slip
Plano de Falla 15 38 -137
Plano Auxiliar 248 65 -60
Mecanismo 24 Región COSTA OESTE COLOMBIA
Fecha 2004/11/15
111504C
Latitud Longitud Profundidad Magnitud
4,72 -77,57 16 7,2
Solución Strike Dip Slip
Plano de Falla 21 11 114
Plano Auxiliar 177 79 85
46
Mecanismo 25 Región COLOMBIA Fecha 2005/03/08
200503081559A
Latitud Longitud Profundidad Magnitud
4,50 -75,92 60,0 5,0
Solución Strike Dip Slip
Plano de Falla 265 70 12
Plano Auxiliar 171 79 159
Mecanismo 26 Región COLOMBIA Fecha 2005/04/21
200504210339A
Latitud Longitud Profundidad Magnitud
5,18 -76,32 116,2 5,2
Solución Strike Dip Slip
Plano de Falla 232 56 1
Plano Auxiliar 141 90 146
Mecanismo 27 Región SUR DE PANAMA Fecha 2007/03/17
200703172243A
Latitud Longitud Profundidad Magnitud
4,61 -78,53 13,1 6,0
Solución Strike Dip Slip
Plano de Falla 55 37 -67
Plano Auxiliar 206 56 -107
Mecanismo 28 Región SUR DE PANAMA Fecha 2007/03/18
200703180211A
Latitud Longitud Profundidad Magnitud
4,69 -78,53 12 6,2
Solución Strike Dip Slip
Plano de Falla 38 36 -69
Plano Auxiliar 193 57 -104
Mecanismo 29 Región COSTA OESTE COLOMBIA
Fecha 2007/09/10
200709100149A
Latitud Longitud Profundidad Magnitud
3,08 -78,12 18,9 6,7
Solución Strike Dip Slip
Plano de Falla 54 23 -95
Plano Auxiliar 239 67 -88
47
Mecanismo 30 Región SUR DE PANAMA Fecha 2008/02/23
200802230323A
Latitud Longitud Profundidad Magnitud
4,06 -78,75 12 5,1
Solución Strike Dip Slip
Plano de Falla 31 38 -98
Plano Auxiliar 221 52 -84
Mecanismo 31 Región COLOMBIA Fecha 2009/09/08
200909080826A
Latitud Longitud Profundidad Magnitud
4,77 -76,71 82,7 4,8
Solución Strike Dip Slip
Plano de Falla 273 46 115
Plano Auxiliar 60 49 66
Mecanismo 32 Región COLOMBIA Fecha 2009/12/18
200912181430A
Latitud Longitud Profundidad Magnitud
3,08 -76,30 169,2 5,1
Solución Strike Dip Slip
Plano de Falla 110 46 -172
Plano Auxiliar 14 84 -44
Mecanismo 33 Región COLOMBIA Fecha 2010/01/29
201001291752A
Latitud Longitud Profundidad Magnitud
4,83 -76,07 128 4,9
Solución Strike Dip Slip
Plano de Falla 195 41 -72
Plano Auxiliar 352 51 -105
Mecanismo 34 Región COLOMBIA Fecha 2013/08/27
201308272322A
Latitud Longitud Profundidad Magnitud
4,08 -76,73 68,1 4,8
Solución Strike Dip Slip
Plano de Falla 63 42 136
Plano Auxiliar 189 63 57
48
Mecanismo 35 Región SUR DE PANAMA Fecha 2013/10/24
201310241444A
Latitud Longitud Profundidad Magnitud
5,21 -78,05 22,5 5,0
Solución Strike Dip Slip
Plano de Falla 82 71 -9
Plano Auxiliar 175 82 -160
Mecanismo 36 Región COLOMBIA Fecha 2014/01/05
201401050336A
Latitud Longitud Profundidad Magnitud
4,66 -76,72 25,5 5,5
Solución Strike Dip Slip
Plano de Falla 288 27 -13
Plano Auxiliar 30 84 -117
Mecanismo 37 Región COLOMBIA Fecha 2014/01/09
201401091250A
Latitud Longitud Profundidad Magnitud
3,02 -75,92 24,7 5,0
Solución Strike Dip Slip
Plano de Falla 43 82 173
Plano Auxiliar 134 83 8
Mecanismo 38 Región COLOMBIA Fecha 2014/09/20
201409200649A
Latitud Longitud Profundidad Magnitud
3,08 -79,41 16,8 4,7
Solución Strike Dip Slip
Plano de Falla 8 35 -106
Plano Auxiliar 206 57 -80
Mecanismo 39 Región COLOMBIA Fecha 2015/02/21
201502212226A
Latitud Longitud Profundidad Magnitud
4,80 -76,06 102,7 5,2
Solución Strike Dip Slip
Plano de Falla 93 51 -179
Plano Auxiliar 2 89 -38
49
Mecanismo 40 Región COSTA OESTE COLOMBIA
Fecha 2015/07/19
201507191116A
Latitud Longitud Profundidad Magnitud
2,96 -79,26 12 5,0
Solución Strike Dip Slip
Plano de Falla 17 10 101
Plano Auxiliar 186 81 88
Mecanismo 41 Región COLOMBIA Fecha 2016/01/22
201601220241A
Latitud Longitud Profundidad Magnitud
5,14 -78,26 18,9 5,0
Solución Strike Dip Slip
Plano de Falla 15 42 -121
Plano Auxiliar 234 55 -65
Mecanismo 42 Región COSTA OESTE COLOMBIA
Fecha 2016/12/06
201612060451A
Latitud Longitud Profundidad Magnitud
4,77 -77,34 31 4,9
Solución Strike Dip Slip
Plano de Falla 78 41 -120
Plano Auxiliar 296 56 -67
Mecanismo 43 Región COLOMBIA Fecha 2018/02/19
201802190150A
Latitud Longitud Profundidad Magnitud
4,24 -76,53 109.0 4,8
Solución Strike Dip Slip
Plano de Falla 60 79 -7
Plano Auxiliar 151 83 -169
4.5. Mapa mecanismos focales vs fallas.
El mapa de Mecanismos Focales y Fallas geológicas se generó en el software ArcGis
10.3.1, se tomó un archivo en formato .xls con los Mecanismos Focales de la zona de
estudio, cada entidad espacial tiene asociado lo información característica del
mecanismos focal como profundidad, magnitud y ángulos para los planos de falla
(Ilustración 13), esta tabla de datos fue exportada a formato shp y layer, este ultimo con el
50
objetivo de generar los balones de playa, ArcGis tiene disponible una herramienta para
generar los balones de playa, desarrollada por Tiesheng Wu, Hainan Research Center of
Engineering Earthquake, la cual se puede descargar libremente, esta herramienta se
encuentra en lenguaje phyton, al momento de generar los balones de playa, debe
seleccionarse los mecanismos en el formato layer, luego se seleccionan las columnas del
archivo correspondientes los ángulos (Strike, Dip y Slip/Rake), esta herramienta genera 5
archivos en formato .shp, uno de los cuales contiene las componentes dilatación, T y
compresión P, que permite visualizar un símbolo como un balón de playa. (Ilustración 15),
otro de los archivos contiene las componentes P y T en formato vectorial, indicando la
dirección y magnitud de las mismas.
Ilustración 13 Tabla de Atributos de Mecanismos Focales utilizados en el desarrollo del estudio.
51
El mapa de fallas se generó a partir de la digitalización tomando como referencia el Mapa
de Fallas y Pliegues Cuaternarias de Colombia y Regiones Oceánicas Adyacentes, en el
marco del Proyecto Internacional de la Litosfera, Grupo de Trabajo II-2, Principales Fallas
Activas en el mundo, un proyecto de Cooperación entre el Servicio Geológico de los
Estados Unidos (USGS) y Gabriel Paris, adicionalmente se compilo junto con la base de
Datos de Fallas Activas de Ingeominas (2001) y la Base cartográfica de Fallas 2007 del
IGAC. A cada falla que se identificó en el mapa se asoció información característica de la
misma (Ilustración 14). Identificando un total de 21 fallas o sistemas de falla, como es el
caso de la Falla Romeral,
Ilustración 14 Tabla de Atributos de Fallas
En el Mapa de Mecanismos Focales VS Fallas, se pueden observar la distribución
espacial de los Mecanismos Focales en naranja y de las Fallas en rojo, se puede
determinar que existen dos mecanismos focales en el territorio del Valle del Cauca el
resto se ubican en los alrededores principalmente en la zona norte y en el mar pacifico en
cercanías a la línea límite de las placas Nazca y Continental.
52
Ilustración 15 Mapa Mecanismos Focales VS Fallas
También mediante el software libre Generic Map Tools GMT, se graficaron los
mecanismos focales de la zona de estudio, usando la función PSMECA se graficaron los
mecanismos focales, teniendo en cuenta que los datos solo tiene una solución del plano,
(Azimut, Deslizamiento y Buzamiento), y el orden de entrada mediante un archivo *.txt es:
coordenadas, profundidad, solución del Mecanismos y magnitud, (Ilustración 11),
adicionalmente se categorizaron por colores, según su profundidad azul= 100 – 300 km,
verde= 30 – 100 km y rojo= 0 – 30 km (Ilustración 12). Y el tamaño del Balón de playa se
establece automáticamente de acuerdo con la magnitud del mecanismo.
53
Ilustración 16 Formato entrada datos Mecanismo a GMT
Finalmente se obtiene como resultado el siguiente mapa de Mecanismos focales en
coordenadas geográficas, en el que se observan los mecanismos de diferentes tamaños,
según su magnitud y con un código de colores según su profundidad.
Este software libre tiene una base cartográfica de las costas y de la división política de los
países, en este caso departamental, la cual se encuentra en tres resoluciones, liviana,
media y completa, para esta salida grafica de los mecanismos focales se seleccionó la
base completa, que fue usada como base para la generación de la cartografía de los
mecanismos focales y de los tensores de esfuerzos. Así mismo se asignó un color azul a
los mecanismos, sin que este signifique algo, solo se busca resaltar el balón de playa que
genera cada uno de los mecanismos, el software admite diferentes formatos para
generación grafica del mecanismo focal, el cual debe seleccionarse de acuerdo con el
archivo de datos seleccionado, para tal caso se utilizó el –Sd convención para los
mecanismos suministrados por Harvar CMT, catalogo utilizado para el desarrollo de este
proyecto, el cual no esta representado por los ángulos ya mencionados anteriormente, si
no por la mejor doble dupla definida del momento tensor con el siguiente encabezado X Y
depth mrr mtt mff mrt mrf mtf exp newX newY [event_title], en donde X y Y son las
coordenadas del mecanismo, la siguiente variable la profundidad, seguida por las seis
componentes del momento tensor, luego la exponencial, las componentes newX y newY
no fueron utilizadas en esta salida grafica, pero sirven para asignar nuevas coordenadas
al balón de playa extendiendo una línea desde su coordenada original y finalmente el
título del evento el cual no fue incluido. Como resultado se obtiene la siguiente salida
gráfica.
54
Ilustración 17 Mapa Mecanismos Focales, Software GTM
4.6. Aplicación del modelo de inversión amortiguada.
El método presentado en el numeral aplicó a datos de mecanismos focales del Valle del
Cuaca y sus alrededores, en total 43 Mecanismos focales fueron usados. Estos
Mecanismos Focales fueron tomados del catálogo Global Centroid Moment Tensor (CMT)
de la Universidad de Harvard y fueron derivados de sismos ocurridos desde 1978 hasta
2018. Este catálogo es una recopilación de los mecanismos focales de terremotos desde
1976 hasta la actualidad y permite realizar una búsqueda rápida mediante la siguiente
página web: http://www.globalcmt.org/CMTsearch.html.
55
Ilustración 18 Interfaz gráfica catalogo CMT
4.7. Mapa tensor de esfuerzos
Los cálculos desarrollados para obtener los tensores de esfuerzos mediante la
metodología anteriormente enunciada son realizados en el software SATI. Este software
fue desarrollado por Hardebeck (2006), y está disponible en la página web del Servicio
Geológico de los Estados Unidos en el siguiente link:
http://earthquake.usgs.gov/research/software/#SATSI. Este software permite realizar una
inversión de los mecanismos focales con datos espaciales y/o temporales para un campo
de esfuerzos, el modelo de inversión busca el campo de esfuerzos menos complejo
consistente con los datos y obtener datos que permitirán calcular los tensores de
esfuerzos de terremotos en 2, 3 y 4 dimensiones, para el caso de este estudio se realizó
en 2 dimensiones, donde se tiene en cuenta las coordenadas del punto (NORTE, ESTE) y
los tres ángulos del mecanismo focal Azimut, Buzamiento y Deslizamiento (Strike, Dip y
Rake) y la magnitud de los mecanismos focales. Las variaciones a 3 y 4 dimensiones
comprenden datos de profundidades y datos temporales de los mecanismos focales.
En la metodología planteada el cálculo de los tensores de esfuerzos se realiza dividiendo
el área de estudio en una grilla, para el presente trabajo se realizó el cálculo de los
tensores de esfuerzos con una grilla de ancho y alto 0,1° Cada Mecanismo Focal es
56
asignado al punto de grilla más cercano, en la tabla 4.7-1 se observa los mecanismos
focales utilizados en el desarrollo.
Mediante la generación de la grilla utilizando el Software ArcMap, se asignó un par de
coordenadas X,Y que se encuentran entre 1 y 39 en la ordenada X y entre 1 y 28 en la
ordenada Y.
Tabla 4.18 Mecanismos Focales
Num Lat Long Depth Strike Dip rake - slip Strike Dip rake - slip Mw
1 3.24 -78.99 26 32 19 125 176 74 79 5.5
2 3.3 -79.31 15.4 19 32 -109 221 60 -79 5.7
3 4.99 -77.68 24.1 14 18 104 179 73 85 5.8
4 4.72 -77.07 79.8 220 51 3 128 88 141 5.7
5 4.8 -77.18 19.1 13 13 95 188 77 89 7.2
6 4.74 -77.48 20.7 45 22 127 186 72 76 5.2
7 5.08 -75.88 125.9 228 22 -71 28 69 -97 5.9
8 4.15 -76.98 90.4 234 44 77 73 47 103 5.2
9 3.21 -78.78 15 10 42 -123 231 56 -63 5.8
10 2.93 -75.94 15 206 76 170 299 80 14 6.8
11 4.09 -76.36 68.2 210 21 -90 30 69 -90 6.3
12 3.02 -78.11 54.1 193 50 51 65 53 127 5
13 5.22 -75.69 128.7 204 30 -99 34 61 -85 6.5
14 4.59 -76.9 118.2 306 51 144 61 63 45 5.3
15 4.78 -76.5 118.1 138 42 105 299 50 77 5.8
16 4.11 -75.84 189.5 249 53 31 140 66 139 6.3
17 4.58 -75.75 27.7 8 65 -21 107 71 -153 6.1
18 4.31 -75.74 33.8 17 67 -23 116 69 -155 5.5
19 4.03 -76.17 178 247 53 31 137 65 139 5.9
20 5.17 -77.62 22.3 12 45 -140 251 62 -52 5.7
21 3.76 -78.9 15 36 40 -70 191 53 -106 5.2
22 3.07 -77.81 31.4 253 45 -73 50 47 -106 5
23 5.14 -77.81 27.6 15 38 -137 248 65 -60 5.9
24 4.72 -77.57 16 21 11 114 177 79 85 7.2
25 4.5 -75.92 60 265 70 12 171 79 159 5
26 5.18 -76.32 116.2 232 56 1 141 90 146 5.2
27 4.61 -78.53 13.1 55 37 -67 206 56 -107 6
28 4.69 -78.53 12 38 36 -69 193 57 -104 6.2
29 3.08 -78.12 18.9 54 23 -95 239 67 -88 6.7
30 4.06 -78.75 12 31 38 -98 221 52 -84 5.1
57
31 4.77 -76.71 82.7 273 46 115 60 49 66 4.8
32 3.08 -76.2 169.2 110 46 -172 14 84 -44 5.1
33 4.83 -76.07 128 195 41 -72 352 51 -105 4.9
34 4.08 -76.73 68.1 63 42 136 189 63 57 4.8
35 5.21 -78.05 22.5 82 71 -9 175 82 -160 5
36 4.66 -76.72 25.5 288 27 -13 30 84 -117 5.5
37 3.02 -75.92 24.7 43 82 173 134 83 8 5
38 3.08 -79.41 16.8 8 35 -106 206 57 -80 4.7
39 4.8 -76.06 102.7 93 51 -179 2 89 -38 5.2
40 2.96 -76.26 12 17 10 101 186 81 88 5
41 5.14 -78.26 18.9 15 42 -121 234 55 -65 5
42 4.77 -77.34 31 78 41 -120 296 56 -67 4.9
43 4,24 -76,53 109.0 60 79 -7 151 83 -169 4.8
Una vez se cuenta con los datos de los mecanismos focales con el par de coordenadas
nuevas de acuerdo con la grilla, se realiza la instalación del Software SATSI el cual tiene
en su última versión la generada en SATSI_140818. El software contiene una serie de
ficheros que deben se iniciados en el sistema operativo Linux, y mediante un archivo de
programación en PERL, ejecutara el cálculo de los tensores de esfuerzos utilizando un
modelo amortiguado, por esta razón antes de procesar los datos, se debe definir el
parámetro de amortiguación, para lo cual se generó la gráfica de la curva de
compensación, la que permite establecer un punto de equilibrio entre la variabilidad de los
datos y la longitud del modelo. Los resultados fueron graficados median el software Excel
de oficce y el software libre Generic Map Tools GMT, y los resultados se muestran a
continuación.
Tabla 4.19 Tabla de Datos para Parámetro de Amortiguación
Parametro Amortiguamiento
Varianza de los Datos
Longitud del Modelo
0 1.16E-07 0.322968
0.02 1.22E-07 0.208743
0.05 3.11E-06 0.162239
0.1 1.73E-05 0.0844702
0.2 0.00027498 0.0841785
0.3 0.00136488 0.0829643
58
0.4 0.00410305 0.0799413
0.6 0.0163855 0.0668515
0.8 0.03447 0.0489844
1 0.0516545 0.0336774
1.2 0.0651129 0.022985
1.4 0.0750277 0.0159513
1.6 0.0822923 0.011309
1.8 0.0876716 0.00817053
2 0.0916923 0.00599492
2.25 0.0953385 0.004147
2.5 0.0978878 0.00292385
2.75 0.099673 0.0021001
3 0.100929 0.00153614
3.5 0.102461 0.00086557
4 0.10327 0.00051899
5 0.103986 0.00021686
6 0.104252 0.00010536
8 0.104424 3.35E-05
10 0.104472 1.37E-05
20 0.104503 8.59E-07
50 0.104505 2.20E-08
Cada punto sobre la gráfica representa un parámetro de amortiguación que va desde 0
hasta 50, para calcular el parámetro de amortiguamiento se utiliza el archivo
invert_damping_parameter.pl, script perl que sirve para automatizar el proceso,
internamente a este archivo hay que indicar el nombre del archivo con los datos de
entrada y el archivo con los datos de salida, este archivo ejecuta el programa
satsi_2D_tradeoff. Con el objetivo de conocer los ajustes que realiza el programa a los
datos mediante el parámetro de amortiguamiento, se realizó el procesamiento de los
datos con parámetro 0, 0,8 y 3, de acuerdo a la gráfica el parámetro de 0,8, es el mejor
parámetro de amortiguamiento ya que presenta la mejor opción para la longitud del
modelo y la varianza, el proceso de clasificación se realizara con los resultados de este
parámetro.
59
Ilustración 19 Curva de Compensación para Grilla 0,1°. Grafica generada en Software GMT
El procesamiento se realiza con los datos de los dos planos de falla de los mecanismos
focales para cada parámetro de amortiguamiento seleccionado, el archivo de salida tiene
tres partes, en la primera parte aparece en Tensor de esfuerzos con las coordenadas de
grilla y las 6 componentes respectivas, en la segunda parte están los datos del
mecanismo focal (strike, dip y rake), el ajuste angular realizado, la magnitud tau y las
coordenadas.
Ilustración 20 Parte 1 archivo de salida SATSI, parámetro de amortiguamiento 0,8
00.02
0.05
0.1
0.2
0.3
0.4
0.60.8
11.21.41.61.822.252.52.753 3.5 4 5 6 8 10 20 50
3,73E-09
1,86E-08
9,31E-08
4,66E-07
2,33E-06
1,16E-05
5,82E-05
2,91E-04
1,46E-03
7,28E-03
3,64E-02
1,82E-01
9,09E-011E-080,00000010,0000010,000010,00010,0010,010,11
VA
RIA
NZA
LONGITUD DEL MODELO
Parametro de AmotiguamientoVarianza VS Longitud
COORDINATES ARE EAST,NORTH,UP.
stress tensors are:
X Y See Sen Seu Snn Snu Suu
2 7 0.523316 0.279767 0.055668 1.031642 -1.004499 -1.554959
3 9 0.542225 0.331735 0.032236 1.041870 -0.914127 -1.584095
60
La tercera parte del archivo indica las estadísticas de los ajustes realizados sobre los
datos de entrada, se indican cuatro parámetros, el ajuste medio del angulo, la desviación
estándar, el promedio de la magnitud Tau y la desviación estándar.
Ilustración 21 Parte 2 archivo de salida SATSI, parámetro de amortiguamiento 0,8
Ilustración 22 Parte 3 archivo de salida SATSI, parámetro de amortiguamiento 0,8
La tabla 4.20, muestra 42 tensores de esfuerzos para los 43 mecanismos focales con los
que se cuenta, cuando en un punto de grilla se adjudica más de un mecanismo focal, el
resultado se ajusta a un solo tensor de esfuerzos en ese punto
Tabla 4.20 Tensor de Esfuerzos
X Y See Sen Seu Snn Snu Suu
2 7 0.523316 0.279767 0.055668 1.031642 -1.004499 -1.554959
3 9 0.542225 0.331735 0.032236 1.04187 -0.914127 -1.584095
6 8 0.196417 0.153742 -0.065032 0.44418 -1.29655 -0.640597
7 14 0.031574 0.246986 -0.319578 1.053137 -0.314345 -1.084711
8 8 0.361469 0.308176 0.205698 0.614738 -0.99595 -0.976207
9 17 0.131696 0.399232 0.01453 0.991425 -0.6156 -1.123121
11 22 -0.020149 0.614597 -0.380069 1.175844 -0.436127 -1.155695
11 23 -0.065079 0.578712 -0.4071 1.149787 -0.440922 -1.084708
13 27 -0.24883 0.946985 -0.010605 0.87572 -0.959847 -0.626889
dip direction, dip, rake, fit angle, mag tau, X, Y
32.0 19.0 125.0 12.3 0.79 6 8
19.0 32.0 -109.0 2.0 0.94 3 9
14.0 18.0 104.0 9.0 0.81 19 26
fit angle mean= 8.768213 standard deviation= 8.497653
avg tau= 0.818444 , std. dev.= 0.118988
61
15 6 -0.063447 -0.168689 0.566362 -0.276396 -0.294827 0.339844
15 7 0.306786 0.147836 0.501728 0.292292 -0.074342 -0.599078
16 28 -0.285297 1.140472 0.126821 0.779495 -0.96301 -0.494198
18 7 0.451828 0.109674 0.247737 0.415705 -0.363124 -0.867533
18 27 -0.189655 0.815923 0.141055 0.737551 -1.108627 -0.547896
19 26 -0.325629 0.514736 -0.057383 0.488411 -1.256114 -0.162782
20 23 -0.345596 0.268487 -0.041812 0.446724 -1.129153 -0.101128
20 28 -0.22595 0.695047 0.175467 0.602945 -1.139525 -0.376995
21 23 -0.394143 0.191195 -0.104246 0.388887 -1.171011 0.005255
23 24 0.014592 0.280089 -0.012485 0.698027 -1.034551 -0.712618
24 24 -0.337027 0.160056 -0.093319 0.518299 -1.077757 -0.181271
25 23 -0.572281 0.069473 0.021909 0.552286 -0.914268 0.019995
26 18 -0.720156 -0.21117 -0.178488 -0.006305 -0.600279 0.726461
27 22 -0.95427 -0.032239 0.020842 0.458443 -0.466481 0.495827
29 17 -0.631045 -0.130154 -0.055917 0.143916 -0.927166 0.487129
29 23 -0.51667 -0.045171 -0.396702 0.420841 -0.748981 0.09583
29 24 -0.828428 -0.261367 -0.217737 0.2965 -0.408733 0.531927
31 24 -0.819658 0.23325 -0.128491 0.492784 -0.652537 0.326874
32 17 -0.819123 0.142864 -0.067013 0.053259 -0.48284 0.765864
33 6 -0.523986 0.201456 -0.385359 0.469206 -0.819714 0.05478
33 28 -0.068138 -0.082458 -0.068488 0.00004 -0.969703 0.068098
34 7 -0.560118 0.149972 -0.137719 0.747037 -0.802586 -0.186919
34 16 -0.004104 -0.159416 -0.057618 0.037856 -0.986426 -0.033752
35 24 -0.600974 0.120013 -0.210708 0.344476 -0.870382 0.256498
37 5 -0.460436 -0.117468 -0.27851 0.862237 -0.96194 -0.4018
37 6 0.680248 0.538619 0.27908 -0.820067 -0.329083 0.139819
37 21 0.481901 0.251786 0.341885 -0.68619 -0.449651 0.204289
37 27 -0.57357 0.246382 -0.20951 0.642667 -0.855377 -0.069097
38 17 -0.403035 0.071442 -0.546691 0.760413 -0.768095 -0.357378
39 19 -0.674102 0.123404 -0.210844 0.396597 -0.900092 0.277506
39 22 -0.816471 -0.38581 -0.437123 0.542066 -0.597462 0.274405
39 28 -0.734786 -0.388235 -0.489863 0.522419 -0.61207 0.212366
62
Los tensores de esfuerzos los vemos en una matriz de la siguiente manera:
(𝑆𝑒𝑒 𝑆𝑒𝑛 𝑆𝑒𝑢𝑆𝑛𝑒 𝑆𝑛𝑛 𝑆𝑛𝑢𝑆𝑢𝑒 𝑆𝑢𝑛 𝑆𝑢𝑢
) (4.1)
A partir de esta matriz 3X3, se calculan las magnitudes de los principales esfuerzos
horizontales, este resultado ya podrá ser graficado en GMT mediante las cruces de
tensión –Sx, el cual requiere el par de coordenadas, el Eps1, componente de tensión
máxima, Eps2 componente horizontal máxima de compresión. y el angulo theta el angulo
de la orientación medido a partir el Norte
Tabla 4.21 Componentes para graficar tensores de esfuerzos mediante la variable –Sx del Software GMT.
ESTE OESTE COMPONENTES
-79.41 3.08 1.000000 0.345750 -66.127275
-79.31 3.3 1.000000 0.312065 -63.491309
-78.99 3.24 1.000000 0.237294 -64.430504
-78.9 3.76 1.000000 -0.022534 -77.097064
-78.78 3.21 1.000000 0.188636 -56.169276
-78.75 4.06 1.000000 -0.021859 -68.557968
-78.53 4.61 1.000000 -0.194840 -67.107833
-78.53 4.69 1.000000 -0.214737 -68.193537
-78.26 5.14 1.000000 -0.556899 -60.349921
-78.11 3.02 0.080021 -1.000000 28.870182
-78.12 3.08 1.000000 0.338565 -43.596788
-78.05 5.21 1.000000 -0.671786 -57.512071
-77.81 3.07 1.000000 0.592044 -40.324134
-77.81 5.14 1.000000 -0.548083 -59.802501
-77.68 4.99 1.000000 -0.779310 -64.167336
-77.57 4.72 1.000000 -0.808880 -72.936789
-77.62 5.17 1.000000 -0.622147 -60.403537
-77.48 4.74 0.988009 -1.000000 -76.985795
-77.34 4.77 1.000000 -0.107159 -70.330108
-77.18 4.8 1.000000 -0.668770 -79.740636
63
-77.07 4.72 0.965320 -1.000000 -86.478255
-76.98 4.15 0.066180 -1.000000 74.694955
-76.9 4.59 0.480812 -1.000000 88.693381
-76.73 4.08 0.253236 -1.000000 80.716447
-76.72 4.66 0.815302 -1.000000 87.247882
-76.71 4.77 0.399758 -1.000000 77.538273
-76.53 4.24 0.619861 -1.000000 -80.216267
-76.5 4.78 0.090340 -1.000000 -80.932510
-76.36 4.09 0.902750 -1.000000 -78.959447
-76.26 2.96 0.447595 -1.000000 56.230357
-76.32 5.18 1.000000 -0.755349 -83.538224
-76.2 3.08 1.000000 -0.810026 48.748673
-76.17 4.03 0.583587 -1.000000 -82.877497
-76.07 4.83 1.000000 -0.539530 84.964023
-75.94 2.93 0.859253 -1.000000 -17.839323
-75.92 3.02 0.723243 -1.000000 -11.660602
-75.92 4.5 1.000000 -0.899958 -78.972222
-75.88 5.08 1.000000 -0.532707 -86.499259
En la tabla de los tensores de esfuerzos podemos observar que la orientación se
encuentra predominando hacia el oriente, ángulos entre -60 y – 80, precisamente estas
variaciones son las que ajusta mediante la amortiguación el proceso desarrollado, a
continuación se presentan los resultados obtenidos para el cálculo de los tensores de
esfuerzos de la zona de estudio mediante el método de inversión amortiguado. Las
componentes de esfuerzo principal horizontal fueron calculados en Octave, que es un
software libre simil de Matlab, mediante un script que calcula las magnitudes de los
vectores principales en las ordenadas X y Y, y la orientación.
64
Ilustración 23 Componentes principales para Tensores de Esfuerzos Calculados sin parámetro de amortiguamiento, para el plano de falla 1.
Ilustración 24 Componentes principales para Tensores de Esfuerzos Calculados sin parámetro de amortiguamiento, para el plano de falla 2.
65
Ilustración 25 Componentes principales para Tensores de Esfuerzos Calculados con un parámetro de amortiguamiento de 0,8, para el plano de falla 1.
Ilustración 26 Componentes principales para Tensores de Esfuerzos Calculados con un parámetro de amortiguamiento de 0,8, para el plano de falla 2.
66
Ilustración 27 Componentes principales para Tensores de Esfuerzos Calculados con un parámetro de amortiguamiento de 3, para el plano de falla 1.
Ilustración 28 Componentes principales para Tensores de Esfuerzos Calculados con un parámetro de amortiguamiento de 3, para el plano de falla 2.
67
Los resultados de una inversión sin amortiguamiento (Ilustración 23 y 24) demuestran la
constricción del esfuerzo a las áreas de las grillas, limitando los esfuerzos a un área
específica sin considerar que los esfuerzos que se generan en su alrededor, mientas que
los resultados con un amortiguamiento de 0,8, se presenta un mejor ajuste a los tensores
de esfuerzos, en el que los ajustes angulares tienen un promedio de 8,7°, una desviación
estándar de 8,49, estos resultados también se evidencian para el segundo plano de falla,
finalmente los resultados para el parámetro de amortiguamiento 3, demuestran un calculo
del tensor de esfuerzos sobre-amortiguado que limita las orientaciones del tensor y sus
magnitudes.
4.8. Clasificación de los tensores de esfuerzos
Tomando como resultado el tensor de esfuerzos calculado con el parámetro de amortiguamiento
0,8, se realiza clasificación de los componentes principales y la orientación para cada tensor,
mediante la herramienta de ArcGis Grouping Analysis, en la que se ingresa el archivo de los
tensores de esfuerzos, genrados a partir de datos X, Y.
Ilustración 29 Herramienta Grouping Analysis
Se establece la cantidad de clases en este caso 4, que tuvo en cuenta las profundidades como
campos de análisis y no se constriño espacialmente, generando 4 clases de acuerdo a la
variabilidad y vecindad de los datos que las componen, posterior a esta clasificación se generaron
68
las regiones con la herramienta Desviación Direccional (Elipse de Desviación Estándar), el cual
calculo una regiones elipsoidales en las que los datos tienen la misma clase, debido a que estas
regiones están limitadas a la geometría del elipsoide se quedan algunos puntos por fuera de la
región a pesar que hacen parte de está.
Ilustración 30 Mapa de Zonificacion de Tensores de Esfuerzos.
69
La clasificación de los tensores teniendo en cuenta su profundidad realiza en 4 rangos de
profundidades, clase 1 de 54,1 a 90,4, clase 2 de 109 a 128,7, clase 3 de 12 a 33,8 y clase 4 de
169,2 a 189,5, esto permite establecer una relación entre la dinámica de la corteza terrestre con
respecto a las profundidades.
Ilustración 31 Tabla Clasificación de Tensores
70
5. Conclusiones
El mapa de esfuerzos activos obtenido, a partir de la inversión de mecanismos
focales, utilizando el método de inversión amortiguada, permite conocer las
variaciones en la dirección del tensor de compresión principal horizontal.
El método de inversión amortiguado ajusta el error que se presenta entre las sub-
áreas adyacentes, representando el régimen de esfuerzos a una escala regional,
la cual es considerada como un medio elástico que es continuo y no se limita a las
sub-áreas establecidas para ejercer la inversión del mecanismo focal.
El mapa de manera general muestra una dirección predominante noreste, con una
pequeña variación en la dirección en la región sureste de la zona de estudio,
manteniéndose en la dirección predominante.
Se evidencia la influencia de la placa de nazca que subduce bajo la placa
suramericana, generándose esfuerzos a lo largo de esta con similitudes en la
magnitud y dirección que se evidencia visualmente.
Es importante tener en cuenta las condiciones en las que se generan los tensores
de esfuerzos y los planos de los mecanismos focales de los que calculamos los
tensores, para así poder entender e interpretar los resultados.
La información cartográfica generada de tensores de esfuerzos mediante el
método de inversión amortiguada de mecanismos focales es fundamental para la
interpretación y análisis del régimen de esfuerzos de la zona de estudio.
La cartografía generada permite que un experto interprete la dinámica de la
corteza terrestre, mediante los tensores de esfuerzos y la correlación con las fallas
y los mecanismos focales.
71
6. Bibliografía
Heidbach, O., Tingay, M., Barth, A., Reinecker, J., Kurfeß, D., Müller, B., The World Stress Map based on the database release 2008, equatorial scale 1:46,000,000, Commission for the Geological Map of the World, Paris, doi:10.1594/GFZ.WSM.Map2009, 2009.
Heidbach, O., Tingay, M., Barth, A., Reinecker, J., Kurfeß, D. and Müller, B., 2008The World Stress Map database release 2008 doi:10.1594/GFZ.WSM.Rel2008, WSM 2008 Data Base
Heidbach, O., Tingay, M., Barth, A., Reinecker, J., Kurfeß, D. & Müller, B., 2010. Global crustal stress pattern based on the World Stress Map database release 2008, Tectonophysics, 462, doi:10.1016/j.tecto.2009.1007.1023.
International Seismological Centre 2015., On-line Bulletin, http://www.isc.ac.uk, Internatl. Seismol. Cent., Thatcham, United Kingdom,
SCEDC. 2013 Southern California Earthquake Center. Caltech.Dataset. doi:10.7909/C3WD3xH1
Michael, A. J. (1984), Determination of stress from slip data: Faults and folds, J. Geophys. Res., 89, 11,517– 11,526.
Tabares, L. M., Mora H., y Salcedo E. 1999. Actividad Sísmica y tasa de deformación sismotectónica en la zona del Viejo Caldas, Colombia. Rev. Acad. Colomb.Cienc. CIENCIAS DE LA TIERRA. 23(88): 359-373. ISSN 0370-3908.
Rodríguez-Pascua, M. A. y De Vicente, G. 2001. Estado de esfuerzos actual en la Cordillera Bética Oriental. Boletín Geológico y Minero, 112 (4): 79-95. ISSN: 0366-0176.
Vargas-Jiménez C. A., Pujades L. G., A. Ugalde y J. A. Canas (2002), Estado de deformación y esfuerzos en el territorio colombiano. Rev. Acad. Colomb. Cienc. 26(100): 373-391; 2002. ISSN 0370-3908.
INGEOMINAS 2005. Estudio de microzonificación símica Santiago de Cali
Muñoz A. y Vicente G. (2006), Notas de clase asignatura Análisis del Estado de esfuerzos activos en la corteza, Programa de Doctorado curso 2005 – 2006
Hardebeck, J. L., Michael A. J. 2006., Damped regional-scale stress inversions: Methodology and examples for southern California and the Coalinga aftershock sequence, J. Geophys. Res., 111, B11310, doi:10.1029/2005JB004144.
C. López Fernández, et al, (2008). Zonación sismo-tectónica del NO de la Península Ibérica. Geo-Temas, Vol. 10.
Reduca (Geología). Serie Tectónica 2010. 2 (6): 1-91, 2010. ISSN: 1989-6557.
Mata S. M. 2011. Análisis sismológico del sismo del 12 de septiembre de 2009, en la costa centro occidental venezolana. Trabajo de grado Magister Scientiae en Física Fundamental, Universidad de los Andes,
72
Omaña I. D. 2011. Modelado cortical con datos magnéticos y sismológicos entre plataforma Deltana y Prisma de Barbados, proyecto de grado Ingeniero Geofísico, Universidad Simón Bolívar, Sartenejas,
Martínez M. R. 2011. Estudio sísmico de la estructura cortical en el Bloque Jalisco a partir de registros locales del proyecto Mars, Tesis Maestro en Ciencias de la Tierra, Universidad Nacional Autónoma de México,
Legrand.D., Tensor de Esfuerzos, notas de clase, Departamento de Geofísica, Universidad de Chile.
J. A. Rodríguez Losada, Mecanismo Focal de un Terremoto, Notas de Clase, Departamento de Edafología y Geología, Universidad de la Laguna, La Laguna, S/C de Tenerife
Stein & Wyssession, Tensor de Momento Sísmico, pp. 239-246. Notas de clase.
Miller M., (2014) Notas de clase Sismología Aplicada y de Exploración, Departamento de Geofísica, Facultad de Ciencias Físicas y Matemáticas. Tomado de <http://mttmllr.com/sismologia.htm> 18072018
Trenkamp R., Kellogg J. N.,Freymueller J.T., Mora H.P. 2002. Wide Plate margin defrmation southern Central America and north western South America, CASA GPS observations. J.South America Earth Sciences 15 (2002) 157-171
Heidbach, O., Tingay, M., Barth, A., Reinecker, J., Kurfeß, D., Müller, B., The World tress Map based on the database release 2008, equatorial scale 1:46,000,000, Commission for the Geological Map of the World, Paris, doi:10.1594/GFZ.WSM.Map2009, 2009.
INPRES, Instituto Nacional de Prevención Sísmica Argentina, Fallas Geológicas, 2014
Estudio general de suelos y zonificación de tierras del departamento del Cauca Escala 1:100.000 IGAC, 2009, Poster WSM
73
ANEXOS
74
ANEXO I – SATSI
El software SATSI es un desarrollo realizado por Jane Hardebeck (2006), la última
versión de este software de libre distribución es del 2014, se puede localizar en la
página del servicio geológico de los Estados Unidos, este software permite realizar
la inversión amortiguada de un Mecanismo Focal, esta inversión puede realizarse
en 2D, con una solución del mecanismo de falla, dado por los tres ángulos Dip,
Strike y Rake, 3D o 4D con información de altura y/o tiempo.
El Software SATSI contiene unos de scripts que deben ser compilados
previamente al uso del software, estos son los siguientes archivos: satsi_2D.c,
satsi_2D_tradeoff.c, satsi_4D.c, satsi_4D_tradeoff.c, satsifast_2D.c,
satsifast_4D.c, bootmech_2D.c, bootmech_4D.c, boot_uncert.c,, esta compilación
se ejecuta mediante el archivo Makefile.
También cuenta con otros archivos de subrutinas que utiliza el software para su
ejecución, estos archivos son: dirplg.c, eigen.c, leasq_sparse.c, myrand.c,
slfast_2D.c, slfast_4D.c, sort.c, stridip.c, switchsub.c.
Finalmente a la hora de procesar se utilizan los archivos:
invert_damping_parameter.pl
invert_smooth.pl
invert_smooth_uncert.pl
Estos archivos funcionan como un icono de otro programa, al compilarlos ejecutan
los archivos .c, automatizando el proceso de ingresar los datos e indicar el archivo
de salida, debe tener en cuenta a los archivos .pl, se deben modificar agregando ./
al inicio del llamado de los archivos .c, esto con el objetivo de que puedan ser
compilados.
Adicionalmente se cuenta con un archivo de test_data.txt, con datos de
mecanismos focales de prueba, así mismo se cuenta con los archivos con los
resultados que generan los archivos de perl.
75
Finalmente mediante el Generic Map Tools, se pueden graficar los resultados, se
debe tener en cuenta que previamente al dibujo de los tensores de esfuerzos se
debe calcular las componentes horizontales principales, para la ejecución de este
proyecto se utilizó el siguiente script compilado en Octave.
sten=[... ...];
for i=1:Cantidad de Tensores
a=sqrt(i);
M=[sten(i,1) sten(i,2) sten(i,3) ; sten(i,2) sten(i,4) sten(i,5) ;
sten(i,3) sten(i,5) sten(i,6)];
[V,D]=eig(M);
[Y,I]=sort([D(1,1) D(2,2) D(3,3)]);
slen=D(I(3),I(3))-D(I(1),I(1));
M2=[sten(i,1) sten(i,2); sten(i,2) sten(i,4)]/slen;
[V,D]=eig(M2);
[Y,I]=sort([D(1,1) D(2,2)]);
S3=V(:,I(1));
trend=(180/pi)*atan2(S3(1),S3(2));
while (trend<-90)
trend=trend+180;
endwhile
while (trend>90)
trend=trend-180;
endwhile
s3mag=(D(I(2),I(2)));
s1mag=(D(I(1),I(1)));
len=max([abs(s3mag) abs(s1mag)]);
s3mag=s3mag/len;
s1mag=s1mag/len;
printf ("%12.6f %12.6f %12.6f %12.6f %12.6f\n " ,
lon(i),lat(i),s3mag,s1mag,trend);
endfor
Los datos de los tensores ingresados en sten, deben ser solo las 6 componentes
del tensor.
76
ANEXO II – GENERIC MAP TOOLS
El software Generic Map Tools es una librería de herramientas que permite
generar salidas gráficas, a continuación se presentan las librerías utilizadas y el
código implementado para la generación de las salidas
Todas las saldas graficas tienen como cartografía base, la división política
departamental de Colombia, para esto se crea un mapa base con la librería
PSCOAST, se utilizaron los datos de la base en formato full. Este mapa base será
llamado en todas las salidas.
Mecanismos Focales
Se utiliza la librería PSMECA.
gmt meca [<table>] -J<args> -R<west>/<east>/<south>/<north>[+r] -S<format><scale>[/<fontsize>[/<justify>/<offset>/<angle>/<form>]] [-B<args>] [-C[<pen>][P<pointsize>]] [-D<depmin>/<depmax>] [-E<fill>] [-G<fill>] [-L<pen>] [-M] [-Fa[<size>[/<Psymbol>[<Tsymbol>]]] [-Fe<fill>] [-Fg<fill>] [-Fo] [-Fr<fill>] [-Fp[<pen>]] [-Ft[<pen>]] [-Fz[<pen>]] [-N] [-T<nplane>[/<pen>]] [-U[[<just>]/<dx>/<dy>/][c|<label>]] [-V[<level>]] [-W<pen>] [-X[a|c|r]<xshift>[<unit>]] [-Y[a|c|r]<yshift>[<unit>]] [-Z<cpt>] [-di<nodata>] [-e[~]<pattern>] [-h[i|o][<nrecs>][+c][+d][+r<remark>][+t<title>]] [-i<cols>[+l][+s<scale>][+o<offset>][,...]] [-t<transp>] [-:[i|o]] [--PAR=<value>] -J = Sistema de Proyeccion del Mapa -Jm|M[<lon0>/[<lat0>/]]<scl>|<width> (Mercator) -R = Coordenadas Minimas y Maximas de la Region de los Datos -R<xmin/xmax/ymin/ymax>+<unit>
77
-S = Tipo de Formato y Tamaño del Simbolo c Focal mechanisms in Harvard CMT convention X Y depth strike1 dip1 rake1 strike2 dip2 rake2 moment newX newY [event_title] with moment in 2 columns : mantissa and exponent corresponding to seismic moment in dynes-cm d Best double couple defined from seismic moment tensor (Harvard CMT, with zero trace): X Y depth mrr mtt mff mrt mrf mtf exp newX newY [event_title] p Focal mechanism defined with: X Y depth strike1 dip1 strike2 fault mag newX newY [event_title] fault = -1/+1 for a normal/inverse fault m Seismic moment tensor (Harvard CMT, with zero trace): X Y depth mrr mtt mff mrt mrf mtf exp newX newY [event_title] Optionally add /fontsize[/offset][u] [Default values are /9/3.000000p] fontsize < 0 : no label written; offset is from the limit of the beach ball. By default label is above the beach ball. Add u to plot it under. -B = Especifica la configuracion del marco del mapa base y los parametros de los ejes. (1) Frame settings are modified via an optional single invocation of -B[<axes>][+g<fill>][+n][+o<lon>/<lat>][+t<title>] (2) Axes parameters are specified via one or more invocations of -B[p|s][x|y|z]<intervals>[+l<label>][+p<prefix>][+u<unit>] <intervals> is composed of concatenated [<type>]<stride>[<unit>][l|p] sub-strings See psbasemap man page for more details and examples of all settings. -C = Offset a la coordenada del Mecanismo Focal especificada en las ultimas dos columnas del del archivo, antes del label, los atributos por default son seleccionados por - W -D = Imprime eventos entre una profundidad y minima y una profundidad maxima -E = Selecciona el color usado para las partes extensivas, por defecto el Blanco -G = Selecciona el color para las partes compresivas por defecto es Negro <r/g/b> (each 0-255) for color or <gray> (0-255) for gray-shade [0]. -L = Selecciona el atributo lapiz para linea extioror otra mas que por defecto seleccionada por -W -M = El mismo tamaño para cualquier magnitud, tamaño es dado con -S.
78
-F = Selecciona varios atributos de los símbolos dependiendo del Modo a Plot axis. Default symbols are circles; otherwise append <size>[/<Psymbol>[<Tsymbol>]. g Append color used for P_symbol [default as set by -G]. e Append color used for T_symbol [default as set by -E]. p Draw P_symbol outline using the default pen (see -W) or append pen attribute for outline. t Draw T_symbol outline using the default pen (see -W) or append pen attribute for outline. o Use psvelomeca format (Without depth in third column). r Draw box behind labels. z Overlay zero trace moment tensor using default pen (see -W) or append outline pen. -N = No salta o corta simbolos que caigan por fuera del borde del mapa (por defecto se igoran los que esten por fuera). -Tn[/<pen>] Draw nodal planes and circumference only to provide a transparent beach ball using the default pen (see -W) or sets pen attribute. n = 1 the only first nodal plane is plotted. n = 2 the only second nodal plane is plotted. n = 0 both nodal planes are plotted. If moment tensor is required, nodal planes overlay moment tensor. -U = Plot Unix System Time stamp [and optionally appended text]. You may also set the reference points and position of stamp [BL/-1.905c/-1.905c]. Give -Uc to have the command line plotted [OFF]. -V = Change the verbosity level (currently v). Choose among 6 levels; each level adds more messages: q - Quiet, not even fatal error messages. n - Normal verbosity: only error messages. c - Also produce compatibility warnings [Default when no -V is used]. v - Verbose progress messages [Default when -V is used]. l - Long verbose progress messages. d - Debugging messages. -W = Set pen attributes [default,black]. -X -Y = Shift origin of plot to (<xshift>, <yshift>). -Z = Use CPT to assign colors based on depth-value in 3rd column.
79
gmt meca [DATOS] -Jm -R-80/-70/4/2 -S<format><scale>[/<fontsize>[/<justify>/<offset>/<angle>/<form>]] [-B<args>] [-C[<pen>][P<pointsize>]] [-D<depmin>/<depmax>] [-E<fill>] [-G<fill>] [-L<pen>] [-M] [-Fa[<size>[/<Psymbol>[<Tsymbol>]]] [-Fe<fill>] [-Fg<fill>] [-Fo] [-Fr<fill>] [-Fp[<pen>]] [-Ft[<pen>]] [-Fz[<pen>]] [-N] [-T<nplane>[/<pen>]] [-U[[<just>]/<dx>/<dy>/][c|<label>]] [-V[<level>]] [-W<pen>] [-X[a|c|r]<xshift>[<unit>]] [-Y[a|c|r]<yshift>[<unit>]] [-Z<cpt>] [-di<nodata>] [-e[~]<pattern>] [-h[i|o][<nrecs>][+c][+d][+r<remark>][+t<title>]] [-i<cols>[+l][+s<scale>][+o<offset>][,...]] [-t<transp>] [-:[i|o]] [--PAR=<value>] Código utilizado para generar la salida grafica pscoast -R-80/-75/2/6 -JM18cm -Ba1f0.5/a1f0.5 -Df -G240 -W1/0 -P -Na -K > meca14.ps psmeca mecaformat -JM -R -Sd0.15 -W1/234/1/0 -Gblue -O >> meca14.ps
Tensor de esfuerzos
Se utiliza la librería PSVELO
psvelo [meca] 6.0.0_r20254 [64-bit] - Plot velocity vectors, crosses, and wedges on maps usage: psvelo [<table>] -J<args> -R<west>/<east>/<south>/<north>[+r] [-A<vecpar>] [-B<args>] [-D<sigscale>] [-G<fill>] [-K] [-L] [-N] [-O] [-P] [-S<symbol><scale><fontsize>] [-U[[<just>]/<dx>/<dy>/][c|<label>]] [-V] [-W<pen>] [-X[a|c|r]<xshift>[<unit>]] [-Y[a|c|r]<yshift>[<unit>]] [-di<nodata>] [-e[~]<pattern>] [-h[i|o][<nrecs>][+c][+d][+r<remark>][+t<title>]] [-i<cols>[+l][+s<scale>][+o<offset>][,...]] [-t<transp>] [-:[i|o]] [--PAR=<value>] -J Select map proJection. (<scale> in cm/degree, <width> in cm) Append h for map height, or +|- for max|min map dimension. Azimuthal projections set -Rg unless polar aspect or -R<...>r is set.
80
-Ja|A<lon0>/<lat0>[/<hor>]/<scl (or <radius>/<lat>)|<width> (Lambert Azimuthal EA) -Jb|B<lon0>/<lat0>/<lat1>/<lat2>/<scl>|<width> (Albers Conic EA) -Jcyl_stere|Cyl_stere/[<lon0>/[<lat0>/]]<lat1>/<lat2>/<scl>|<width> (Cylindrical Stereographic) -Jc|C<lon0>/<lat0><scl>|<width> (Cassini) -Jd|D<lon0>/<lat0>/<lat1>/<lat2>/<scl>|<width> (Equidistant Conic) -Je|E<lon0>/<lat0>[/<horizon>]/<scl (or <radius>/<lat>)|<width> (Azimuthal Equidistant) -Jf|F<lon0>/<lat0>[/<horizon>]/<scl (or <radius>/<lat>)|<width> (Gnomonic) -Jg|G<lon0>/<lat0>/<scl (or <radius>/<lat>)|<width> (Orthographic) -Jg|G[<lon0>/]<lat0>[/<horizon>|/<altitude>/<azimuth>/<tilt>/<twist>/<Width>/<Height>]/<scl>|<width> (General Perspective) -Jh|H[<lon0>/]<scl>|<width> (Hammer-Aitoff) -Ji|I[<lon0>/]<scl>|<width> (Sinusoidal) -Jj|J[<lon0>/]<scl>|<width> (Miller) -Jkf|Kf[<lon0>/]<scl>|<width> (Eckert IV) -Jks|Ks[<lon0>/]<scl>|<width> (Eckert VI) -Jl|L<lon0>/<lat0>/<lat1>/<lat2>/<scl>|<width> (Lambert Conformal Conic) -Jm|M[<lon0>/[<lat0>/]]<scl>|<width> (Mercator) -Jn|N[<lon0>/]<scl>|<width> (Robinson projection) -Jo|O (Oblique Mercator). Specify one of three definitions: -Jo|O[a|A]<lon0>/<lat0>/<azimuth>/<scl>|<width> -Jo|O[b|B]<lon0>/<lat0>/<lon1>/<lat1>/<scl>|<width> -Jo|Oc|C<lon0>/<lat0>/<lonp>/<latp>/<scl>|<width> -Jpoly|Poly/[<lon0>/[<lat0>/]]<scl>|<width> ((American) Polyconic) -Jq|Q[<lon0>/[<lat0>/]]<scl>|<width> (Equidistant Cylindrical) -Jr|R[<lon0>/]<scl>|<width> (Winkel Tripel) -Js|S<lon0>/<lat0>/[<horizon>/]<scl> (or <slat>/<scl> or <radius>/<lat>)|<width> (Stereographic) -Jt|T<lon0>/[<lat0>/]<scl>|<width> (Transverse Mercator) -Ju|U[<zone>/]<scl>|<width> (UTM) -Jv|V<lon0>/<scl>|<width> (van der Grinten) -Jw|W<lon0>/<scl>|<width> (Mollweide) -Jy|Y[<lon0>/[<lat0>/]]<scl>|<width> (Cylindrical Equal-area) -Jp|P[a]<scl>|<width>[/<origin>][r|z] (Polar [azimuth] (theta,radius)) -Jx|X<x-scl>|<width>[d|l|p<power>|t|T][/<y-scl>|<height>[d|l|p<power>|t|T]] (Linear, log, and power projections) (See psbasemap for more details on projection syntax) -R Specify the min/max coordinates of your data region in user units. Use dd:mm[:ss] for regions given in arc degrees, minutes [and seconds]. Use -R<xmin/xmax/ymin/ymax>+<unit> for regions given in projected coordinates. with <unit> selected from e|f|k|M|n|u. Use [yyy[-mm[-dd]]]T[hh[:mm[:ss[.xxx]]]] format for time axes.
81
Append +r if -R specifies the coordinates of the lower left and upper right corners of a rectangular area. Use -R<gridfile> to use its limits (and increments if applicable). Use -Rg and -Rd as shorthands for -R0/360/-90/90 and -R-180/180/-90/90. Derive region from closed polygons from the Digital Chart of the World (DCW): Append a comma-separated list of ISO 3166 codes for countries to set region, i.e., <code1>,<code2>,... etc., using the 2-character ISO country codes (see pscoast -E+l for list). To select a state of a country (if available), append .state, e.g, US.TX for Texas. To select a whole continent, give =AF|AN|AS|EU|OC|NA|SA as <code>. Use +r to modify the region from polygon(s): Append <inc>, <xinc>/<yinc>, or <winc>/<einc>/<sinc>/<ninc> to round region to these multiples; use +R to extend region by those increments instead [0]. Or use -R<code><x0>/<y0>/<n_columns>/<n_rows> for origin and grid dimensions, where <code> is a 2-char combo from [T|M|B][L|C|R] (top/middle/bottom/left/center/right) and grid spacing must be specified via -I<dx>[/<dy>] (also see -r). OPTIONS: <table> is one or more data files (in ASCII, binary, netCDF). If no files are given, standard input is read. -B Specify both (1) basemap frame settings and (2) axes parameters. (1) Frame settings are modified via an optional single invocation of -B[<axes>][+g<fill>][+n][+o<lon>/<lat>][+t<title>] (2) Axes parameters are specified via one or more invocations of -B[p|s][x|y|z]<intervals>[+l<label>][+p<prefix>][+u<unit>] <intervals> is composed of concatenated [<type>]<stride>[<unit>][l|p] sub-strings See psbasemap man page for more details and examples of all settings. -A Specify arrow head attributes: Append length of vector head, with optional modifiers: [Left and right are defined by looking from start to end of vector] +a<angle> to set angle of the vector head apex [30] +b to place a vector head at the beginning of the vector [none]. Append t for terminal, c for circle, s for square, a for arrow [Default], i for tail, A for plain arrow, and I for plain tail. Append l|r to only draw left or right side of this head [both sides]. +e to place a vector head at the end of the vector [none]. Append t for terminal, c for circle, s for square, a for arrow [Default], i for tail, A for plain arrow, and I for plain tail.
82
Append l|r to only draw left or right side of this head [both sides]. +g<fill> to set head fill or use - to turn off fill [default fill]. +h sets the vector head shape in -2/2 range [0]. +j<just> to justify vector at (b)eginning [default], (e)nd, or (c)enter. +l to only draw left side of all specified vector heads [both sides]. +m[f|r] to place vector head at mid-point of segment [Default expects +b|+e]. Specify f or r for forward|reverse direction [forward]. Append t for terminal, c for circle, s for square, or a for arrow [Default]. Append l|r to only draw left or right side of this head [both sides]. +n<norm> to shrink attributes if vector length < <norm> [none]. +o[<plon/plat>] sets pole [north pole] for great or small circles; only give length via input. +p[-][<pen>] to set pen attributes, prepend - to turn off head outlines [default pen and outline]. +q if start and stop opening angle is given instead of (azimuth,length) on input. +r to only draw right side of all specified vector heads [both sides]. +s if (x,y) coordinates of tip is given instead of (azimuth,length) on input. +t[b|e]<trim(s)>[unit] to shift begin or end position along vector by given amount [no shifting]. Default is 9p+gblack+p1p -D Multiply uncertainties by <sigscale>. (Se and Sw only)i -E Set color used for uncertainty wedges in -Sw option. -G Specify color (for symbols/polygons) or pattern (for polygons). fill can be either 1) <r/g/b> (each 0-255) for color or <gray> (0-255) for gray-shade [0]. 2) p[or P]<iconsize>/<pattern> for predefined patterns (0-90). -K Allow for more plot code to be appended later [CLASSIC MODE ONLY]. -L Draw line or symbol outline using the current pen (see -W). -N Do Not skip/clip symbols that fall outside map border [Default will ignore those outside]. -O Set Overlay plot mode, i.e., append to an existing plot [CLASSIC MODE ONLY]. -P Set Portrait page orientation [OFF]; [CLASSIC MODE ONLY]. -S Select symbol type and scale. Choose between: e Velocity ellipses: in X,Y,Vx,Vy,SigX,SigY,CorXY,name format. r Velocity ellipses: in X,Y,Vx,Vy,a,b,theta,name format. n Anisotropy : in X,Y,Vx,Vy. w Rotational wedges: in X,Y,Spin,Spinsig. x Strain crosses : in X,Y,Eps1,Eps2,Theta. -U Plot Unix System Time stamp [and optionally appended text]. You may also set the reference points and position of stamp [BL/-1.905c/-1.905c]. Give -Uc to have the command line plotted [OFF]. -V Change the verbosity level (currently v). Choose among 6 levels; each level adds more messages: q - Quiet, not even fatal error messages.
83
n - Normal verbosity: only error messages. c - Also produce compatibility warnings [Default when no -V is used]. v - Verbose progress messages [Default when -V is used]. l - Long verbose progress messages. d - Debugging messages. -W Set pen attributes [default,black]. -X -Y Shift origin of plot to (<xshift>, <yshift>). Prepend r for shift relative to current point (default), prepend a for temporary adjustment of origin, prepend f to position relative to lower left corner of page, prepend c for offset of center of plot to center of page. For overlays (-O), the default setting is [r0], otherwise [f2.54c]. -di Replace any <nodata> in input data with NaN. -e Only accept input data records that contain the string "pattern". Use -e~"pattern" to only accept data records that DO NOT contain this pattern. If your pattern begins with ~, escape it with \~. To match against extended regular expressions use -e[~]/regexp/[i] (i for case-insensitive). Give +f<file> for a file list with such patterns, one per line. To give a single pattern starting with +f, escape it with \+f. -h[i][<n>][+c][+d][+r<remark>][+t<title>] Input/output file has [0] Header record(s) [OFF] Optionally, append i for input only and/or number of header records [0]. -hi turns off the writing of all headers on output. Append +c to add header record with column information [none]. Append +d to delete headers before adding new ones [Default will append headers]. Append +r to add a <remark> comment to the output [none]. Append +t to add a <title> comment to the output [none]. (these strings may contain \n to indicate line-breaks) For binary files, <n> is considered to mean number of bytes. -i Sets alternate input column order and optional transformations [Default reads all columns in order]. Append list of columns; t = trailing text. Use -in for just numerical input. -t Set the layer PDF transparency from 0-100 [Default is 0; opaque]. -: Swap 1st and 2nd column on input and/or output [OFF/OFF]. -^ (or -) Print short synopsis message. -+ (or +) Print longer synopsis message. -? (or no arguments) Print this usage message. --PAR=<value> Temporarily override GMT default setting(s) (repeatable). (See gmt.conf man page for GMT default parameters). psvelo [<table>] -J<args> -R<west>/<east>/<south>/<north>[+r] [-A<vecpar>] [-B<args>] [-D<sigscale>] [-G<fill>] [-K] [-L] [-N] [-O] [-P] [-S<symbol><scale><fontsize>]
84
[-U[[<just>]/<dx>/<dy>/][c|<label>]] [-V] [-W<pen>] [-X[a|c|r]<xshift>[<unit>]] [-Y[a|c|r]<yshift>[<unit>]] [-di<nodata>] [-e[~]<pattern>] [-h[i|o][<nrecs>][+c][+d][+r<remark>][+t<title>]] [-i<cols>[+l][+s<scale>][+o<offset>][,...]] [-t<transp>] [-:[i|o]] [--PAR=<value>] Código utilizado para generar la salida grafica psvelo tensor2018 -JM17cm -R0/40/0/30 -A -Ba5 -D -Gred -K -L -O -P -Sx -W