Universidad de Chile
Facultad de Ciencias Físicas y Matemáticas
Departamento de Geofísica
GF4001 - Sismología
Tarea N°3
Cálculo de Magnitudes Local
y de Momento.
Estudiante:
Martín Fuentealba
Profesor:
Mario Pardo
Auxiliares:
Daniela Benavente
Bruno Quezada
28 de junio de 2019
Santiago, Chile.
GF4001 Sismología Tarea N°1
Tabla de contenido
1. Introducción ................................................................................................................................ 2
1.1 Antecedentes ....................................................................................................................... 2
1.2 Problema a resolver ............................................................................................................. 3
1.3 Memoria de cálculo ............................................................................................................. 3
1.4 Alcances .............................................................................................................................. 4
1.5 Objetivos ............................................................................................................................. 4
2. Metodología ................................................................................................................................ 5
3. Resultados ................................................................................................................................... 7
4. Conclusión................................................................................................................................... 9
5. Referencias ................................................................................................................................ 10
6. Anexos....................................................................................................................................... 10
GF4001 Sismología Tarea N°1
1. Introducción
1.1Antecedentes
Los terremotos se pueden producir dentro de un amplio espectro de tamaños. Sismicidad se induce
en la minería subterránea1 que involucra fracturas del orden de 100 metros mientras que la zona de
fractura producida durante el evento de Valdivia 1960 alcanzó ordenes de magnitud de hasta 106
metros.
Dado lo anterior, la magnitud es un parámetro que mide el tamaño de un terremoto,
relacionada con la energía liberada en forma de ondas elásticas. De esta manera es posible la
comparación entre eventos. Se tiene que la magnitud es una variable cuantificable y, por lo tanto,
única y determinada de manera objetiva. La magnitud se determina mediante registros sismográficos
del evento.
Variadas son las escalas expuestas en la literatura, debido a que dependen del tipo de onda
considerado, del periodo de esta, del instrumento con el que se registra la onda, la distancia al foco,
entre otros. Dentro de las más usadas se encuentran la magnitud body wave 𝑚𝑏 , la magnitud de ondas
superficiales 𝑀𝑆 , la magnitud local (o de Richter) 𝑀𝑙 y la magnitud de momento sísmico 𝑀𝑊 . Este
reporte se centra en las últimas dos únicamente.
La magnitud local fue propuesta por (Richter, 1935), quien la define calibrada a partir de
terremotos de California, Estados Unidos, ocurridos en un radio de 600 kilómetros, registrados en
sismógrafo de torsión de Wood-Anderson.
Esta escala de magnitud se basa sobre el concepto de que, conociendo la distancia del
sismógrafo al terremoto y determinando la máxima amplitud de la señal registrada por el sismógrafo,
es posible estimar el tamaño relativo del evento según una ecuación empírica.
Las limitaciones de la magnitud local de momento se asocian a que en California los sismos
son, en general, superficiales. Lo anterior, sumado a que la expresión no considera correcciones por
profundidad, significa que no se estima de manera precisa el tamaño de los sismos de profundidad
considerable.
Por otra parte, las escalas 𝑚𝑏 y 𝑀𝑆 tienen la limitación de que, ante eventos con extensas
zonas de ruptura que liberan mucha energía en un largo periodo, subestiman el tamaño real del sismo.
Para solucionar el problema anterior, (Hanks & Kanamori, 1979) proponen una escala que se
determina a partir del momento sísmico. El momento sísmico es un parámetro que da cuenta de la
geometría de la falla, su rigidez y el desplazamiento de los bloques. Es una medida más consistente
para estimar el tamaño del terremoto. Finalmente, la magnitud se obtiene a partir una expresión
empírica.
1
Principalmente por tronadura. Estallidos de roca producidos por un campo de esfuerzos alto (alto SRF) puede
generar zonas de fracturamiento mayores.
GF4001 Sismología Tarea N°1
La ventaja de 𝑀𝑊 es que no se satura. Por lo tanto, es idónea para la comparación entre
grandes eventos. (Pardo, 2019)
1.2 Problema a resolver
Se pide calcular la magnitud local y de momento para un sismo ocurrido en el norte de Chile, cuyas
coordenadas para su hipocentro se presentan en la Tabla 1.
Tabla 1: Hipocentro de sismo a estudiar
Latitud [°] -20.12
Longitud [°] -69.33
Profundidad [°] 90
Se tienen los registros sismográficos de cinco estaciones, cuyas coordenadas y distancia
epicentral al evento se explicita en la Tabla 2. Para la magnitud local se tienen los sismogramas sin
la respuesta instrumental y deconvolucionados en velocidad y los sismogramas sintéticos de Wood-
Anderson. Para la magnitud de momento se tienen los sismogramas sin la respuesta instrumental en
desplazamiento, los sismogramas en desplazamiento de la fase S y los registros de amplitud de la
transformada de Fourier en frecuencia.
Se considera que la respuesta de la superficie libre (F) y el patrón de radiación promedio (R)
se asumen como constantes y tienen valores de F = 2 y R = 0.67.
Tabla 2: Ubicación estaciones sismográficas
Estación Latitud [°] Longitud [°] 𝚫 [Km]
MNMCX -19.13 -69.59 113.7
PB 07 -21.72 -69.88 185.9
PSGCX -19.60 -70.12 100.7
PATA -20.69 -70.00 93.0
UNAP -20.24 -70.14 84.8
1.3 Memoria de cálculo
• Magnitud local 𝑀𝑙
𝑀𝑙 = log 𝐴 − log 𝐴0 (1)
Donde:
𝐴 [𝑚𝑚]: Amplitud máxima de un sismo medida por un
sismómetro Wood-Anderson.
𝐴0 [𝑚𝑚]: Amplitud máxima de un sismo medido por un
W-A que cumple que a 100 km de distancia
epicentral esta es 0.001 mm.
GF4001 Sismología Tarea N°1
• Curva de atenuación log(𝐴0 )
Δ
log 𝐴0 = − log ( )−3 (2)
100
Donde:
Δ [𝐾𝑚]: Distancia epicentral.
• Magnitud de momento 𝑀𝑊
2 (3)
𝑀𝑊 = log 𝑀0 − 6,07
3
Donde:
𝑀0 [𝑁𝑚]: Momento sísmico.
• Momento sísmico (en Nm)
4𝜋𝑣 3 𝜌 (4)
𝑀0 = ⋅ Δ ⋅ √Ω2E + Ω2𝑁 + Ω2𝑍
𝐹⋅𝑅
Donde:
𝑣 [𝐾𝑚/𝑠]: Velocidad de la onda.
𝜌 [𝑔/𝑐𝑚3 ]: Densidad del medio a la profundidad de la
fuente.
F [-] Respuesta de la superficie libre.
R [-] Patrón de radiación promedio.
Ω [𝑚]: Valor del plateau observado en el espectro de
frecuencias para las componentes (E, N, Z).
1.4Alcances
En el presente reporte se presenta la metodología empleada, los resultados y gráficos obtenidos y las
conclusiones que surgen a partir del trabajo realizado.
1.5 Objetivos
El objetivo de este trabajo es el de calcular la magnitud local y de momento para un sismo ocurrido
en el norte de Chile, a partir de los registros de cinco estaciones.
GF4001 Sismología Tarea N°1
2. Metodología
En primer lugar, se ubica el epicentro del evento, así como también, las cinco estaciones sismográficas
con las que se trabaja. En la Figura 1 se muestran sus ubicaciones.
Figura 1: Ubicación de epicentro y estaciones sismográficas
Luego, se procede al cálculo de la magnitud del evento según la escala local y la escala de
momento sísmico.
Cálculo de magnitud local
Para el cálculo de la magnitud local, primero se grafícan los sismogramas sin la respuesta
instrumental, deconvolucionados a velocidad, presentados anexos en este informe en el Gráfico 1. A
partir de estos, se obtienen los sismogramas sintéticos de Wood-Anderson, presentados en el Gráfico
2 de la sección Anexos.
GF4001 Sismología Tarea N°1
Con los datos de W-A, se determina la mitad de la máxima amplitud peak-to-peak para cada
componente (E, N, Z).
Una vez obtenidas las amplitudes A, y calculada la curva de atenuación para cada estación,
se obtiene la magnitud local a partir de la ecuación (1). La magnitud final será el promedio de las
cinco estaciones.
Cálculo de magnitud de momento
Al igual que para el caso anterior, se grafican en primer lugar los registro sin la respuesta instrumental
(ver Gráfico 3 en Anexos). A partir de estos, se obtienen los sismogramas que muestran sólo la fase
S, presentados en el Gráfico 4. Se aplica la transformada de Fourier en el dominio de las frecuencias
y se obtiene el registro de amplitud.
Se grafica el último registro en escala log-log, el cual se muestra en el Gráfico 5. Se determina
el valor del plateau de la amplitud para el espectro de bajas frecuencias de cada estación en las tres
componentes (Gráfico 6, Gráfico 7, Gráfico 8, Gráfico 9 y Gráfico 10).
Con el valor anterior, y con la distribución de velocidad de la onda S y densidad en
profundidad, se obtiene el momento sísmico del evento. Aplicando la ecuación (3) se obtiene la
magnitud de momento para el terremoto estudiado. La magnitud final será el promedio de las cinco
estaciones.
GF4001 Sismología Tarea N°1
3. Resultados
Cálculo de magnitud local
Se obtiene la curva de atenuación y las magnitudes locales por estación, las cuales se presentan en la
siguiente tabla.
Tabla 3: Magnitud local en las cinco estaciones sismográficas
Estación 𝐥𝐨𝐠 𝑨𝟎 𝑴𝒍
MNMCX -3.06 5.68
PB 07 -2.97 5.42
PSGCX -3.27 5.10
PATA -3.00 5.72
UNAP -2.93 5.37
Lo anterior resulta en una magnitud 𝑀𝑙 de 5.46 promedio.
Cálculo de magnitud de momento
Luego de determinar los valores para el plateau en el gráfico loglog de la transformada de Fourier en
frecuencia, se obtiene que las amplitudes varían en el rango de 1 ⋅ 10−4 y 3 ⋅ 10−4 metros.
Considerando la distribución de velocidad de la onda S y la densidad en profundidad
mostradas en la Figura 2 y Figura 3, se calcula el momento sísmico y, a partir de este, la magnitud de
momento para cada estación. Se tiene que la información de las figuras anteriores llega hasta los 70
km únicamente. Debido a esto, se determina la magnitud de momento para dos casos: el primero
considerando las propiedades a la profundidad 70 km y el segundo considerando que las propiedades
varían linealmente según el ajuste presentado en las figuras y obteniéndolas para 90 km de
profundidad. Los valores de magnitud para ambos casos se presentan en la tabla Tabla 4.
Velocidad de onda de corte [km/s]
0 1 2 3 4 5 6
0
10
20
Profundidad [km]
30
40
50
60 y = 37,759x - 119,19
R² = 0,9452
70
80
Figura 2: Distribución de velocidad de onda de corte en profundidad
GF4001 Sismología Tarea N°1
Densidad [g/cm3]
0 0,5 1 1,5 2 2,5 3 3,5 4
0
10
20
Profundidad [km]
30
40
50
60
y = 90,668x - 241,63
70 R² = 0,6477
80
Figura 3: Distribución de densidad en profundidad
Tabla 4: Magnitud de magnitud en las cinco estaciones sismográficas
Estación 𝑴𝑾 (𝒛 = 𝟕𝟎 𝒌𝒎) 𝑴𝑾 (𝒛 = 𝟗𝟎 𝒌𝒎)
MNMCX 5.34 5.49
PB 07 5.29 5.45
PSGCX 5.30 5.46
PATA 5.28 5.44
UNAP 5.23 5.38
A partir de lo presentado en la tabla anterior, resulta que para el primer caso se tiene una
magnitud de momento de 5.29 promedio y para el segundo caso de 5.44 promedio.
GF4001 Sismología Tarea N°1
4. Conclusión
A partir del trabajo realizado según la metodología presentada y los resultados reportados en este
informe, es posible señalar que:
• El cálculo de la magnitud local resulta muy conveniente para la estimación rápida de la
magnitud de eventos sísmicos. A partir del filtrado de los registros y su procesamiento para
generar los sismogramas sintéticos de W-A, es posible la determinación de la magnitud de
un sismo en un corto periodo de tiempo y con pocas variables que analizar.
• Por otro parte, y luego del filtrado y procesamiento de los registros sin la respuesta
instrumental, es necesaria la determinación de la distancia epicentral de la estación al foco y,
a partir de esta, la obtención de modelos de distribución de velocidades de onda y densidad
del material para poder determinar la magnitud de momento. Si bien la metodología para esta
escala considera pasos adicionales, se cuantifican los efectos propios de la profundidad del
evento.
• Lo anterior se confirma al notar el leve aumento en la magnitud de momento del caso 1 al
caso 2, producido por el aumento en la velocidad de onda de corte y la densidad debido al
aumento en la profundidad considerada.
• Si se compara el resultado obtenido para la magnitud local y para la magnitud de momento
del segundo caso, es posible observar son bastante similares, con una variación de sólo 0.02
puntos de magnitud. Lo anterior sugiere que los efectos de la profundidad se cuantifican bien
en la escala de magnitud local, a partir de la curva de atenuación, para la magnitud del sismo
estudiado. Para realizar una comparación más completa, es necesario considerar eventos de
mayor magnitud y con variabilidad en la profundidad del foco para notar la influencia de esta
sobre la magnitud.
GF4001 Sismología Tarea N°1
5. Referencias
Hanks, T., & Kanamori, H. (1979). Hanks, T. C., & Kanamori, H. (1979). A moment magnitude
scale. Journal of Geophysical Research: Solid Earth, 84(B5), 2348-2350.
Pardo, M. (2019). Clase X - Magnitud. Santiago, Chile.
Richter, C. (1935). An instrumental earthquake magnitude scale. Bulletin of the Seismological Society
of America, 25(1), 1-32.
6. Anexos
Gráfico 1: Sismogramas de las cinco estaciones en componentes E, N, Z sin la respuesta instrumental deconvolucionados en velocidad
GF4001 Sismología Tarea N°1
Gráfico 2: Sismogramas sintéticos de W-A de las cinco estaciones en componentes horizontales
GF4001 Sismología Tarea N°1
Gráfico 3: Sismogramas de las cinco estaciones en coordenadas E, N, Z sin la respuesta instrumental en desplazamiento
GF4001 Sismología Tarea N°1
Gráfico 4: Sismograma de las cinco estaciones en coordenadas E, N, Z en desplazamiento para la fase S
GF4001 Sismología Tarea N°1
Gráfico 5: Transformada de Fourier en frecuencia de las cinco estaciones en componentes E, N, Z
GF4001 Sismología Tarea N°1
Gráfico 6: Transformada de Fourier en frecuencia de la estación MNMCX en componentes E, N, Z
GF4001 Sismología Tarea N°1
Gráfico 7: Transformada de Fourier en frecuencia de la estación PB07 en componentes E, N, Z
GF4001 Sismología Tarea N°1
Gráfico 8: Transformada de Fourier en frecuencia de la estación PSGCX en componentes E, N, Z
GF4001 Sismología Tarea N°1
Gráfico 9: Transformada de Fourier en frecuencia de la estación PATA en componentes E, N, Z
GF4001 Sismología Tarea N°1
Gráfico 10: Transformada de Fourier en frecuencia de la estación UNAP en componentes E, N, Z
