UNIVERSIDAD COMPLUTENSE DE MADRID

FACULTAD DE CIENCIAS FSICAS

DEPARTAMENTO DE FSICA DE LA TIERRA,
ASTRONOMA Y ASTROFSICA I
(GEOFSICA Y METEOROLOGA)

PROGRAMA OFICIAL DE POSTGRADO DE FSICA:

MSTER DE GEOFSICA Y METEOROLOGA
TRABAJO DE INVESTIGACIN

LOCALIZACIN DE TERREMOTOS A
PARTIR DE FORMAS DE ONDA
Almudena Gomis Moreno
Directora: Dra. Elisa Buforn Peir
Junio 2007

AGRADECIMIENTOS
En primer lugar quisiera agradecer al Departamento de Fsica de la
Tierra, Astronoma y Astrofsica I de la Facultad de Fsicas de la Universidad
Complutense de Madrid el haberme facilitado los medios para poder llevar a
cabo este trabajo de investigacin del Mster de Geofsica y Meteorologa. En
particular, querra mostrar mi ms profundo agradecimiento a la directora de
este trabajo, la Dra. Elisa Buforn Peir, sin cuya ayuda y consejos no se
podra haber realizado.
Este trabajo ha sido llevado a cabo gracias al proyecto ERSE
(Escenarios Realistas de Riesgo Ssmico en Espaa) (REN2003-5178-C03-01),
financiado por el Ministerio de Ciencia y Tecnologa (en la actualidad
Ministerio de Educacin y Ciencia) y al proyecto RISTE (Riesgo de Terremotos
y Tsunamis en Espaa) (CGL2006-10311-C03-01/BTE), financiado por el
Ministerio de Educacin y Ciencia.
Gracias a Hlne Lyon-Caen, Francesco Pacchiani y Aurore Franco, de
la Ecole Normal Suprieur de Pars, por dedicar su tiempo a ensearme los
conceptos fundamentales y los algoritmos, y facilitarme los programas para
llevar a cabo este trabajo.
Quisiera agradecer tambin al Instituto Geogrfico Nacional su
colaboracin con la aportacin de los registros e informacin instrumental
necesarios para el desarrollo de este proyecto.
Mis compaeros de estos dos aos han sido una parte importantsima
de este proceso. Gracias por todos los consejos, los nimos, los fantsticos
momentos que hacen que todo sea ms fcil, y por todas las cosas que he
aprendido de vosotros: Tatiana, Carmen, Juan Luis, Paco, Dani, Jacques,
Esther, Ana, Javi. Todos habis puesto vuestro granito de arena. Tatiana,
muchsimas gracias por el apoyo en todos los sentidos, incluso logstico, hasta
el ltimo minuto.

Gracias a Salva y Luca, que, en estos dos aos y hasta el ltimo

momento, me han sacado de ms de un apuro. Sin todas esas pequeas
ayudas del da a da no se podra acabar ningn trabajo de este tipo.
De manera especial quiero expresarle mi agradecimiento a Joaqun por
sus ratos de escucha, consejos e inters, cuando la geofsica no le toca ni de
cerca, y por sus intentos de atraerme al lado oscuro. Tentaciones de ese tipo
siempre son bienvenidas. Qu bonita es la Fsica!, verdad?
Carmen, gracias, porque has hecho que sienta que esta etapa de mi
vida no est aislada de todo lo dems, sino que sea continuidad de todo lo que
pasamos juntas en el colegio.
Llevo dos semanas intentando encontrar el agradecimiento adecuado
para Juan, y todo resulta insuficiente. Solo s que me siento muy afortunada
de que seas t el que tenga que aguantar mis enfados (bendita paciencia), y
con el que pueda compartir todas mis alegras.
En ltimo lugar quiero darles las gracias a mis padres, del los que he
aprendido a querer aprender. Me habis dado todo vuestro apoyo en los
momentos ms difciles y me habis ayudado a elegir el camino correcto en los
momentos claves.

Y puede que este sea un buen principio. Principio de incertidumbre*

Principio de incertidumbre. Ismael Serrano.

ii

LOCALIZACIN DE TERREMOTOS A
PARTIR DE FORMAS DE ONDA

Trabajo de Investigacin elaborado por Almudena Gomis Moreno, encuadrado

en el Programa Oficial de Postgrado de Fsica: Mster en Geofsica y Meteorologa,
impartido por los Departamentos de Fsica de la Tierra, Astronoma y Astrofsica I
y II de la Facultad de Ciencias Fsicas de la Universidad Complutense de Madrid.

Madrid, 15 de junio de2007

Almudena Gomis Moreno

iii

NDICE

1. Introduccin ....................................................................................... 1
2. Metodologa ........................................................................................ 3
2.1. Identificacin de familias o multipletes....................................... 6
2.2. Correlacin cruzada ................................................................... 8
2.3. Relocalizacin de los eventos.................................................... 10
3. Implementacin del algoritmo........................................................... 12
3.1. Determinacin de las familias de sismos .................................. 12
3.1.1. Clculo de la transformada wavelet de las ondas P .......... 12
3.1.2. Normalizacin de la wavelet............................................. 13
3.1.3. Clculo de las distancias de Levenstein ........................... 14
3.1.4. Creacin de un diagrama de rbol ................................... 14
3.1.5. Extraccin de los multipletes de sismos .......................... 15
3.2. Clculo de la correlacin cruzada............................................. 16
3.2.1. mkhdr ............................................................................. 16
3.2.2. mkdat ............................................................................. 16
3.2.3. mklist.............................................................................. 16
3.2.4. mktbl .............................................................................. 17
3.2.5. dbh ................................................................................. 17
3.3. Relocalizacin .......................................................................... 17
4. Aplicacin a la serie ssmica de 2002 de Grgal (Almera).................. 18
4.1. Introduccin ............................................................................ 18
4.2. Procesamiento de datos............................................................ 20
4.3. Resultados ............................................................................... 23
5. Conclusiones .................................................................................... 34
6. Bibliografa ....................................................................................... 36
7. Anexo 1: sismograma y sonograma ................................................... 39

iv

NDICE DE APNDICES

Se adjunta un CD-ROM en el que se tiene acceso a los siguientes apndices:

A. Programa en .bash para la automatizacin de los programas param, scale,
dist-spares y cluster-sparse
B. Programa para el clculo de coordenadas UTM
C. Tiempos de llegada manuales de la fase P
D. Resultados de la relocalizacin con Hypo71
E. Diagramas de rbol de la distancia de Levenstein
F. Extraccin de multipletes
G. Archivos con extensin .hdr
H. Archivos con extensin .dat
I. Otros
I.1. Lista de sismos de Grgal :gergal.lis
I.2. Lista de sismos de la familia A: mult0001.lis
I.3. Lista de sismos de la familia B: mult0002.lis
I.4. Informacin de la familia A: mult0001.table
I.5. Informacin de la familia B: mult0002.table
I.6. Retardos calculados para la familia A
I.7. Retardos calculados para la familia B

1. INTRODUCCIN
La ocurrencia de terremotos es uno de los fenmenos naturales que mayores
daos econmicos y humanos puede causar. El estudio detallado de la ocurrencia de
sismos es fundamental para poder prevenir y mitigar los efectos destructores de los
mismos. Uno de los problemas fundamentales en el estudio de terremotos es su
localizacin espaciotemporal, que permitir identificar las fuentes sismognicas de la
regin y asociarlas a las fallas activas.
Los primeros mtodos y ms sencillos para la localizacin de un sismo
(determinacin hipocentral) fueron los mtodos grficos, que requeran un mnimo de
tres estaciones para obtener la localizacin. Posteriormente se desarrollaron los
mtodos numricos.
Hoy en da, la mayora de los mtodos numricos que se utilizan se basan en la
minimizacin de los residuos entre los tiempos tericos y los observados de llegadas
de las diversas fases del sismograma: las ondas P, S, Lg, etc. Esta minimizacin
realizada para diferentes modelos de Tierra permite obtener las coordenadas del
epicentro, la profundidad y la hora origen del terremoto. Este es el fundamento de los
mtodos ms utilizados, como el Hypo71 (Lee y Lahr, 1972), HYPOINVERSE (Klein,
1978), etc.
De todos los parmetros a determinar, la profundidad es el que, con mayor
frecuencia, queda peor determinado, con grandes errores, al no converger el sistema
de ecuaciones y debiendo fijarse este parmetro. Este problema se presenta
normalmente cuando existe una mala geometra de la red y se ve muy afectado por las
discontinuidades en los modelos de Tierra, estructurados en capas plano-paralelas de
velocidad, lo que le resta exactitud a la relocalizacin hipocentral. Recientemente se
han desarrollado nuevas metodologas basadas en la utilizacin de formas de onda,
que permiten determinar con mayor precisin la profundidad.
Los mtodos que utilizan formas de onda se basan en la comparacin de la
forma de onda de distintos sismos en cada una de las estaciones en las que se han
registrado, agrupndolos en familias segn la similitud de sus formas de onda
(Poupinet et al., 1984; Ito, 1990). Una vez identificadas las diferentes familias, se
procede a computar el retardo entre las diferentes seales, usando mtodos de
correlacin cruzada (Waldhauser y Ellsworth, 2002; Hauksson y Shearer, 2005;
Shearer et al., 2005). La hiptesis en que se basa el mtodo consiste en que la
diferencia de tiempos de llegada de una misma fase de dos sismos, con un mecanismo
focal semejante y prximos en el espacio, es debida a la diferencia hipocentral de los

mismos. Conocido el retardo ser posible calcular la posicin relativa de los

hipocentros de estos terremotos mediante distintos procesos, como por ejemplo la
doble diferencia de Waldhauser y Ellsworth (2000). La aplicacin de estos mtodos
permite una relocalizacin ms exacta que la obtenida con los mtodos tradicionales,
minimizando la influencia del modelo de Tierra en la relocalizacin.
Hay que tener en cuenta que, debido a que las zonas ssmicamente ms activas
de Espaa son el sur de la Pennsula Ibrica, el mar de Alborn, el Golfo de Cdiz y el
norte de Marruecos, la cobertura azimutal de estaciones no es ptima, lo que produce
en algunos casos grandes errores en la profundidad. Una metodologa como la que
aqu se plantea permitir mejorar la relocalizacin hipocentral, lo que implica una
mejor definicin de las zonas sismognicas y su posterior aplicacin a estudios de
sismotectnica o de riesgo ssmico.
El objetivo de este trabajo es la relocalizacin por medio de formas de onda de
terremotos ocurridos en la regin Ibero-Magreb (sur de la Pennsula Ibrica, Golfo de
Cdiz, mar de Alborn y norte de frica). En concreto se centrar en la denominada
serie de Grgal (Almera), que comenz el 4 de febrero de 2002 con la ocurrencia de un
terremoto principal de magnitud mbLg = 5.1 que fue seguido por un total de 39 sismos
que se extendieron hasta el da 26 de dicho mes. Para llevar a cabo este proceso se
van a emplear registros digitales de banda ancha y de perodo corto.
La serie de Grgal ha sido determinada previamente por el Instituto Geogrfico
Nacional (IGN) utilizando mtodos tradicionales de lectura de fases P y S. Los
resultados obtenidos en este trabajo utilizan formas de onda, y sern comparados con
los resultados obtenidos con la metodologa tradicional utilizada por el IGN. El trabajo
se va a desarrollar segn el siguiente esquema:
1. identificacin de posibles familias,
2. determinacin hipocentral de parmetros focales por medio de formas de onda,
3. comparacin de los resultados con los obtenidos por el IGN.

2. METODOLOGA
La localizacin de un terremoto consiste en determinar las coordenadas
epicentrales, profundidad y hora origen del sismo (xo, yo, zo, to), minimizando en la
medida de lo posible sus errores. Para un grupo de terremotos, la localizacin puede
ser: a) absoluta, en la que se localiza de manera absoluta el hipocentro de los sismos,
o b) relativa, en la que los sismos se relocalizan en relacin a otro sismo o al
baricentro de un grupo de sismos. La precisin de la localizacin absoluta depende de
muchos factores, tales como la geometra de la red, la precisin de la toma del tiempo
de llegada de la fase y el conocimiento del medio (Pavlis, 1986).
El mtodo que se va a emplear en este trabajo se basa en la determinacin del
tiempo de retardo de una misma fase de dos sismos diferentes registrados en la misma
estacin, lo que se llama doble diferencia (Waldhauser y Ellsworth, 2000). Un esquema
de esta situacin se muestra en la figura 1. Puede decirse que la doble diferencia
proporciona una medida de la distancia entre los dos hipocentros (Menke, 1999). El
resultado es una relocalizacin relativa de los sismos. La ventaja de utilizar
relocalizaciones relativas es que minimizan los errores introducidos por el medio
(Poupinet et al., 1984; Frechet, 1985; Got et al., 1994).

Figura 1. Disposicin de los hipocentros y la estacin. DD indica doble diferencia.

Se define el retardo entre dos seales, o doble diferencia, como la diferencia del
tiempo de recorrido desde el foco del sismo a una misma estacin para una fase
determinada en dos sismos diferentes. Este retardo se debe nicamente a la diferencia
en las localizaciones hipocentrales de los sismos que las han generado, si se cumplen
las siguientes condiciones:
1. que sus mecanismos focales sean similares,
2. que sus focos estn prximos, en comparacin con la distancia foco-estacin,
de modo que las ondas hayan seguido aproximadamente la misma trayectoria
y hayan atravesado el mismo medio. Un criterio a seguir es que dos sismos
tendrn sismogramas incoherentes cuando sus hipocentros estn separados
ms de la cuarta parte de la longitud de la frecuencia predominante en el
sismograma (Geller y Mueller, 1980).
El conocimiento del medio siempre es necesario, lo que constituye una
dificultad aadida a los procesos de relocalizacin, ya que el medio atravesado por las
ondas rara vez es bien conocido. Pero en el caso de dos sismos con hipocentros muy
prximos la influencia del medio es mucho menor, ya que la influencia de las
anomalas de velocidad es igual para ambos y, por tanto, no tienen influencia directa
en la relocalizacin. No obstante, en ausencia de errores proporcionados por el medio,
en la relocalizacin siempre estara presente el error de la determinacin manual de
fases mediante un programa de procesado ssmico, que en los mejores casos es de tan
solo 0.1s (Lyon-Caen et al., 2004), lo que conlleva un error epicentral de 500 m y de
alrededor de 1 Km en profundidad. Por tanto, la gran ventaja del empleo de la
correlacin cruzada para el clculo de la doble diferencia es que las incertidumbres de
las coordenadas hipocentrales disminuyen en ms de un orden de magnitud en
comparacin con otras tcnicas que no emplean la doble diferencia (Waldhauser y
Ellsworth, 2000), ya que se elimina la incertidumbre que acarrea determinar de
manera manual el tiempo de llegada de las distintas fases que van a intervenir en la
relocalizacin. La precisin temporal de este mtodo puede llegar a ser de 1 ms frente
al rango de 10-30 ms en el caso de la determinacin manual del tiempo de llegada, lo
que proporciona una relocalizacin relativa entre terremotos con un error que oscila
entre unos metros a unas decenas de metros. El error del que no nos vamos a poder
desprender es el introducido por la geometra de la red. Una mayor cobertura azimutal
proporciona unos mejores resultados.

En el proceso de relocalizacin, el clculo del retardo se puede llevar a cabo

bsicamente de dos maneras distintas:
1. tomando el tiempo de llegada de la fase directamente sobre el sismograma y
calculando la diferencia,
2. empleando el mtodo de correlacin cruzada entre dos seales que son
similares.
La ventaja de la correlacin cruzada es que proporciona menor error que la
lectura de fases de un sismograma. La correlacin cruzada determina directamente el
retardo entre dos sismos, lo que conlleva asociado un error ms pequeo. Algunos
autores (Pacchiani, 2006) denominan relocalizacin relativa a la relocalizacin llevada
a cabo con la correlacin cruzada y relocalizacin absoluta a la que emplea el tiempo
de llegada de cada fase, si bien no hace referencia al proceso de relocalizacin en s.
Es decir, a si los sismos son relocalizados en funcin de otro sismo o de un baricentro,
o si lo son de manera absoluta.
El procedimiento a seguir en este trabajo es el que se muestra a continuacin:
1. Identificacin de familias o multipletes. Un multiplete, o familia, se define
como un conjunto de terremotos cuyos sismogramas tienen un alto grado de
similitud debido a que han tenido lugar muy prximos en el espacio y con un
mecanismo focal muy parecido. Gracias a esta caracterstica, se puede calcular
la doble diferencia para una misma fase mediante la correlacin cruzada. A dos
sismos pertenecientes a una misma familia se les denomina doblete.
2. Correlacin cruzada. Es la comparacin de la forma de onda de los terremotos
de un mismo multiplete. Se procede al clculo del retardo de la misma fase en
cada par de sismogramas. Esta diferencia temporal es fundamental para la
relocalizacin hipocentral. Es la etapa ms importante del proceso de
relocalizacin, pues de su precisin depender el xito de la relocalizacin.
3. Relocalizacin de los eventos. Mediante cualquier mtodo que emplee el
retardo en la llegada de una fase para cada par de estaciones. Si se relocalizan
los sismos en funcin de un sismo principal el mtodo se denomina master
event (Ito, 1985; Deichmann y Garca-Fernndez, 1992). Si por el contrario la
relocalizacin se lleva a cabo empleando todos y cada uno de los sismos de la
familia como evento maestro, el mtodo se

denomina master event

generalizado. Asimismo se puede emplear el programa de relocalizacin

HypoDD (Waldhauser y Ellsworth, 2000), en el que se puede llevar a cabo tanto

una relocalizacin a partir de retardos de correlacin cruzada, como con

retardos obtenidos a partir de medidas absolutas de los tiempos de llegada,
como con una combinacin de ambos tipos. La solucin se halla mediante
mnimos cuadrados, de manera iterativa.
A continuacin se describe en detalle cada una de estas etapas del proceso de
relocalizacin.
2.1. Identificacin de familias o multipletes
Para la agrupacin de los sismos en distintas familias es necesaria la obtencin
previa de sus sonogramas. Se define un sonograma como la representacin del
contenido en frecuencias de un sismograma frente al tiempo. Esta representacin se
obtiene al aplicar transformadas wavelet a los sismogramas correspondientes a los
eventos ssmicos de inters.
La transformada wavelet consiste en ir analizando la seal temporal
(sismograma) usando unas determinadas funciones (wavelet madre) de distintos
tamaos que se trasladan en el tiempo a lo largo de la seal original. As se conocer
el contenido de las distintas frecuencias para cada instante de tiempo.
La definicin matemtica de la transformada wavelet es la que sigue:

Ws ( , s ) =

1
s

t
x(t ) *
dt
s

[1]

donde es el parmetro de traslacin, s es el parmetro de escala, x(t) es la seal

procesada y (t) la wavelet madre. Hay muchos tipos diferentes de wavelets madre.
En este trabajo se va a usar la wavelet de Morlet o gaussiana (Goupillaud, Grossman y
Morlet,1984), definida como:

(t ) = e

t 2

e i 2f ot

[2]

que se encuentra centrada en to en el dominio del tiempo y en fo en el de frecuencias,

con una duracin t=s to en tiempo y f= fo/s en frecuencias.

A partir de la wavelet madre van a ir surgiendo otras funciones con distintos

valores de y de s para llevar a cabo la transformada. El parmetro de traslacin, ,
est relacionado con la posicin de la funcin wavelet sobre la seal que queremos
transformar: se refiere a la localizacin de la ventana de anlisis. Corresponde a la
informacin en el tiempo que va a proporcionar la transformada. Por su parte el
parmetro de escala, s, est relacionado con el contenido en frecuencias: el
escalamiento de una seal corresponde a una compresin o a una dilatacin de la
wavelet madre, por lo que, si la wavelet madre es estrecha, se estar obteniendo
informacin de las altas frecuencias, y por el contrario, si es ancha, se obtiene
informacin de las bajas frecuencias. En la figura 2 se muestra un ejemplo de una
seal y de una wavelet madre para valores de s de 1 y 5 , y valores de to de 2, 40, 90 y
140 s, donde la zona sombreada en naranja es la wavelet que se desplaza a lo largo de
la seal.

Figura 2. Aplicacin de una wavelet (zona sombreada) sobre una seal, para
diferentes valores de s y to. Las unidades del eje de ordenadas son cuentas, y en
abscisas segundos. (Modificado de Polikar, 1994).
La transformada wavelet se aplicada a cada uno de los sismogramas con los
que vamos a trabajar. El resultado de la transformacin wavelet es un sonograma.
El siguiente paso es determinar qu sismos forman parte de un mismo
multiplete, y para ello es necesario estudiar la similitud que hay entre los sonogramas
que se han obtenido al aplicar la transformada wavelet a cada uno de los
sismogramas. La similitud viene determinada por la distancia de Levenstein, que se
define como el nmero de cambios, supresiones o inserciones que se han de llevar a
cabo para transformar un patrn (en este caso los sonogramas) en otro (Levenstein,
1965). Para comprender mejor este proceso supongamos el siguiente ejemplo (figura
3): dos patrones formados por cuatro cuadrados de distinto color cada uno. Para

transformar uno en el otro bastar con intercambiar la posicin del cuadro blanco y la
del naranja. En este caso la distancia de Levenstein es 1.

Figura 3. Ejemplo del clculo de la distancia de Levenstein.

Por tanto, el siguiente paso a seguir es calcular la distancia de Levenstein para
cada par de sonogramas de una misma estacin.
Pero hay que tener en cuenta lo siguiente: cada sismo ha sido registrado en
varias estaciones, por lo que no es inmediato que dos sonogramas que sean similares
en una estacin impliquen que los sismos correspondientes sean similares. Para
concluir que dos terremotos forman parte de un mismo multiplete ser necesario un
alto grado de similitud de los sonogramas en el mayor nmero de estaciones posible.
Este nmero de estaciones ha de ser determinado dependiendo del problema que se
est tratando.
2.2. Correlacin cruzada
Una vez se han obtenido los multipletes o familias de sismos se procede a
calcular el retardo de las parejas de seales pertenecientes a un multiplete mediante
una correlacin en el dominio de las frecuencias (Poupinet et al, 1984; Got et al.,
1994), si bien tambin se puede llevar a cabo en el dominio temporal. La similitud
entre dos seales viene determinada por el ndice de correlacin, que vara entre 1 y
1. Un valor de 1 indica la mxima similitud, 0 la mnima. Un valor negativo significa
que la seal est invertida.
Para comprender el proceso partamos de dos seales dependientes del tiempo,
s1(t) y s2(t), entre las que se cumple la siguiente relacin:

s2 (t ) = C te s1 (t + )

[3]

Para ver el grado de correlacin entre dos seales, se comparan en el dominio

de las frecuencias. Por tanto, el siguiente paso es expresar las seales en el dominio
de la frecuencia, para lo que se aplica la transformada de Fourier:
+

S1 ( f ) = s1 (t )e 2ift dt

S 2 ( f ) = s 2 (t )e

2ift

dt = C

te

s1 (t )e

2if (t + )

dt = C S1 ( f ) e
te

2if

[4]

El grado de correlacin entre dos seales en el dominio de la frecuencia viene

dado por el ndice de coherencia:

(f )=

S12

[5]

A1 A2

donde S12 es el interespectro, y A1 y A2 los autoespectros, que se definen como

S12 ( f ) = S1 ( f ) S 2* ( f ) = C te S1 ( f ) e 2if

[6]

A1 = S1 S1*

[7]

A2 = S 2 S 2*

[8]

Para el clculo del retardo entre dos seales que son suficientemente similares
se emplea la fase del interespectro, que se obtiene de la expresin [6]:

12 ( f ) = 2f =

[9]

El retardo entre las dos seales es la pendiente de la recta de la fase del

interespectro frente a la frecuencia

. La pendiente est ponderada mediante unos

pesos w(f) que dependen del ndice de coherencia, dando as mayor peso all donde
existe mayor semejanza entre las dos seales:

2( f )
w( f ) =
(1 2 ( f ))

[10]

Por tanto, la regresin en cada estacin est limitada a un nmero

determinado de frecuencias alrededor del mximo del interespectro. Para que una
frecuencia se pueda considerar apta para tomar parte en la regresin, es necesario
que cumpla dos condiciones:
1. que la energa asociada a esa frecuencia sea por lo menos el 75 % de la energa
mxima del sismograma,
2. que la coherencia sea al menos de 0.8.
2.3. Relocalizacin de los eventos
Una vez que se han obtenido los retardos para cada par de sismos ya se puede
proceder a la relocalizacin. Para ello se puede usar cualquier mtodo que emplee la
doble diferencia.
En general, el tiempo de recorrido de una onda a una estacin se puede
expresar en funcin de las coordenadas del sismo (xo, yo, zo), las de la estacin i-sima
(xi, yi, zi), la hora origen (to) y la velocidad, v, como

t obs = t o +

( x i x o )2 + ( y i y o )2 + ( z i z o )2
v

= to +

i
= F (xo , y o , z o , t o )
v

[11]

En esta ecuacin son cuatro las incgnitas que se han de determinar

(xo, yo, zo, to). Se supone conocida la distribucin de velocidades y que la lectura del
tiempo de llegada de las fases est libre de error, pero esto no es cierto, por lo que es
preferible tener un sistema sobredeterminado. Incluso si esto no pudiera ser, se puede
llegar a encontrar la mejor solucin posible por medio de una inversin. Lo ms
habitual es tratar de minimizar iterativamente el residuo mediante mnimos
cuadrados. Se define el residuo como

Ri = t i

obs

ti

cal

[12]

donde i hace referencia al sismo. El tiempo calculado de llegada de una determinada

fase se puede expresar como

t cal = F ( xi , y i , z i , t i )

10

[13]

Como el problema directo no es lineal, para resolver el inverso se procede a

desarrollar F(xi, yi, zi, ti) en serie de Taylor (Geiger, 1910) en torno a unos parmetros
iniciales (xo+dxo, yo+dyo, zo+dzo, to+dto). Para facilitar la notacin llamemos

r
m =(m1, m2,

m3, m4)=(x0, y0, z0, t0); entonces el tiempo calculado, si nos quedamos en primer orden,
se puede expresar como

(x x )
(y y )
(z z )
r
Rki = t obs Fi m 0 = i o o x0 + i o o y 0 + i o o z 0 t 0
v i
v i
v i

( )

[14]

donde el subndice k hace referencia a la estacin.

La expresin anterior se puede abreviar como

Rki =

t ki r
r mi
m

[15]

Esta ecuacin no es apropiada para emplear la doble diferencia, pero se puede

definir una diferencia de residuos para dos sismos, i y j, en la misma estacin k-sima
(Frechet, 1985):

dR = t t
ij
k

i
k

j obs
k

t t
i
k

j cal
k

t ki r t kj r
= r mi r mi
m
m

[16]

En esta ocasin se pueden emplear tanto las dobles diferencias, como los
tiempos de llegada absolutos obtenidos de la observacin de sismogramas.
Este

problema

se

resuelve

mediante

inversin,

bien

mediante

la

descomposicin en valores singulares (SVD) para los grupos de sismos pequeos, bien
por el mtodo de los gradientes conjugados (LSQR), (Paige y Saunders, 1982) para
sistemas con ms ecuaciones.

11

3. IMPLEMENTACIN DEL ALGORITMO

La metodologa descrita en el apartado anterior se ha aplicado mediante
diversos programas implementados por varios autores. Estos programas han sido
adaptados para poder usarlos con los datos utilizados en este estudio, siendo
necesario modificar la razn de muestreo de nuestras seales (t), el formato de los
nombres de nuestros archivos de entrada, el rango de frecuencias de inters en
nuestro trabajo, y la estructura de carpetas, es decir, la disposicin de los directorios
y archivos en el ordenador. A continuacin se describen brevemente las principales
caractersticas de los algoritmos utilizados de acuerdo con la metodologa discutida en
el apartado anterior.
3.1. Determinacin de las familias de sismos.
Se ha utilizado la metodologa SPARS (Syntactic Pattern Recognition Scheme),
(Zhizhin, 1992), de reconocimiento sintctico de patrones. Para ello se ha
implementado un programa en .bash (Apndice A del CD-ROM que se adjunta) para la
ejecucin automtica de los cuatro primeros programas que obtienen las familias de
sismos. A continuacin se describen brevemente los algoritmos utilizados y los
parmetros de entrada, as como los resultados que proporciona cada uno de ellos.
3.1.1. Clculo de la transformada wavelet de las ondas P.
Se ha realizado mediante el algoritmo param (programa lenguaje en C; Zhizhin,
1992). Los parmetros de entrada de este programa son los siguientes:
1. n, el nmero de intervalos de frecuencia en los que se va a estudiar el
contenido frecuencial (oscila entre 1 y 20). Se han tomado 20 intervalos
para hacer un estudio ms exhaustivo de los sismogramas,
2. u, lmite superior del intervalo de frecuencias, expresado en Hz. No
debe ser mayor que la frecuencia de Nyquist, que se define como fN=1/2t,
donde t es el intervalo de muestreo, que en este caso es t = 0.02, y por
tanto fN = 25 Hz. Se ha tomado u=10 Hz, ya que se va a proceder a filtrar
los sismogramas entre 1 y 10 Hz, debido a que se ha observado que de este
modo se obtiene una correlacin cruzada ptima (Haukkson y Shearer,
2005),

12

3. l, frecuencia inferior del intervalo de frecuencias, tambin en Hz. No

puede ser menor que la frecuencia de resolucin, que se define como la
inversa de la duracin de la seal. Tomamos l=1 Hz, debido al filtrado
paso-banda al que sometimos a los sismogramas, tal y como se ha
explicado en el prrafo anterior,
4. d, duracin de las ventanas temporales. Debe ser mayor que el ritmo
de muestreo, pero menor que la duracin total del intervalo. Tomamos un
valor de 0.032 s, que es el ms comnmente empleado, ya que en nuestros
sismogramas t=0.02 s.
El resultado final de este programa es la transformada wavelet de la seal
cortada. En nuestro caso, la fase utilizada es la onda P. Se definen los parmetros
IOFFSET (seal antes de la llegada de la onda P) y WFLENGTH (longitud total de la
seal). Los valores utilizados son 20 y 256 puntos, respectivamente, lo que asegura
que la ventana utilizada contenga la onda P.
3.1.2. Normalizacin de la wavelet.
El algoritmo empleado en esta etapa es scale (programa en lenguaje C; Zhizhin,
1992). Con l se normaliza la transformada wavelet y se convierten las unidades a dB.
El resultado es el sonograma de la seal ssmica correspondiente. En la figura 4 se
muestra un ejemplo del sonograma para la componente z del terremoto del 5 de
febrero de 2002 en la estacin ACLR. En el sismograma, el eje de ordenadas son
cuentas, y en el sonograma, el intervalo de frecuencia empleado. En ambos casos el

Cuentas digitales

eje de abscisas son segundos.

t (s)

Figura 4. En la parte superior se muestra la ventana temporal utilizada para la

obtencin del sonograma parte inferior-del terremoto del 5 de febrero de 2002 en
la estacin ACLR. Se puede consultar esta figura en mayores dimensiones en el
Anexo 1 al final de este trabajo.

13

3.1.3. Clculo de las distancias de Levenstein.

Esta etapa se llev a cabo mediante el algoritmo distpar (programa en lenguaje
C; Zhizhin, 1992 modificado por Gaucher, 1996), cuya entrada son los archivos
KnowBaseName y TestBaseName, que contienen la lista de eventos que se van a
procesar, y LocBaseName, que contiene la lista de eventos con su localizacin
preliminar.
Asimismo, se implement un programa en MATLAB (Apndice B del CD-ROM)
que transforma coordenadas geogrficas en coordenadas UTM empleando el elipsoide
de referencia del WGS84 (sistema de coordenadas mundiales que data de 1984),
necesarias para el clculo de la distancia euclidiana, que es la que es calculada por el
programa.
El resultado de este programa son las distancias de Levenstein entre todos los
sismogramas, comparndolos de 2 en 2.
3.1.4. Creacin de un diagrama de rbol.
Este paso se lleva a cabo mediante el programa cluster-sparse (programa en
lenguaje C; Zhizhin, 1992).
Ordena las distancias de Levenstein obtenidas por el mtodo anterior en un
diagrama de rbol que muestra la similitud entre los sismos. Los sismos relacionados
con las ramas ms cortas del rbol son ms semejantes que los relacionados mediante
las ramas largas, entendindose que pueden pertenecer a una misma familia. En la
figura 5 se muestra un ejemplo de diagrama de rbol para la estacin EBAN. Las
distancias de Levenstein aparecen a la derecha, oscilando entre 3.396105 y
1.577106.

14

Tree[root].Dist = 1577359.000000
2002.02.04-20.13.32-------------------------------->
8.682e+05+-------->
2002.02.04-21.37.53----------------------->
|
|
6.240e+05+--------<
|
2002.02.04-20.44.31------------------->
|
|
5.010e+05+---<
|
2002.02.05-05.28.20-------------------<
|
1.106e+06+---------------->
2002.02.05-00.19.41----------------------->
|
|
6.242e+05+--->
|
|
2002.02.14-03.16.10-----------------------<
|
|
|
7.185e+05+>
|
|
2002.02.04-20.31.02----------------------->
||
|
|
6.184e+05+---<|
|
|
2002.02.10-02.42.39-----------------------<
|
|
|
7.585e+05+------------<
|
2002.02.10-19.28.02---------------------->
|
|
5.963e+05+-----<
|
2002.02.04-20.17.44---------------------->
|
5.740e+05+
|
2002.02.05-13.51.22----------------------<
|
1.577e+06+
2002.02.05-12.41.16--------------------------------->
|
8.840e+05+------------------------<
2002.02.04-21.43.21----------------------->
|
6.009e+05+---->
|
2002.02.04-22.11.25-----------------------<
|
|
7.586e+05+----<
2002.02.04-22.00.03------------->
|
3.396e+05+--------------<
2002.02.05-21.42.00-------------<

Figura 5. Ejemplo de un diagrama de rbol para la estacin EBAN. Los nmeros a

la izquierda identifican al sismo. Los restantes corresponden a las distancias de
Levenstein entre unos sismos y otros. Los crculos rojos indican la mayor y la
menor distancia.
3.1.5. Extraccin de los multipletes de sismos.
Esta fase se ha realizado mediante el programa mkclus (programa en lenguaje
C; adaptado de Pacchiani, 2006).
Obtiene las familias de sismos mediante la comparacin de los diagramas de
rbol de las distintas estaciones obtenidos en el apartado anterior. Determina qu
sismos se parecen y en qu estaciones, a fin de, imponiendo un nmero mnimo de
estaciones en las que los sismogramas deben parecerse, extraer los multipletes de
manera objetiva. A continuacin se describen los parmetros que controlan el proceso,
as como sus valores ptimos:
1. nmero de estaciones en que se busca coincidencia: dos sesmos similares
deben tener registros parecidos en el mximo nmero de estaciones posible.
Restringir mucho este parmetro provoca familias poco pobladas, pero se
asegura la similitud de los sismos que la constituyen. El inconveniente es que

15

se pueden discriminar sismos que pertenecen a esa misma familia, pero de los
que no se tiene un registro suficientemente bueno. Por el contrario, un valor
pequeo produce familias muy numerosas formadas por sismos que no tienen
por qu ser similares realmente. Por tanto, es necesario llegar a un
compromiso. En este trabajo se ha estudiado el rango comprendido entre 1 y 5
estaciones,
2.

distancia de Levenstein: indica el grado de similitud que se requiere para

considerar que en una estacin dos sismos son coincidentes. A menor umbral,
ms parecidos sern los sismos. Se ha trabajado con valores que oscilaban
entre un valor de 4105 y 10105.
Las familias propuestas por mkclus han sido confirmadas visualizando los

sismogramas correspondientes.
3.2. Clculo de la correlacin cruzada.
Como paso previo a la correlacin cruzada es necesario obtener una primera
localizacin. Se ha utilizado el programa Hypo71, que proporciona un fichero de salida
con extensin .vis, del que los programas utilizados en la etapa de correlacin leern
la informacin de cada sismo, como es el nombre de las estaciones en las que se tiene
registro, el azimut o el ngulo de incidencia. Para llevar a cabo esta relocalizacin es
necesario tomar en los sismogramas la hora de llegada de las fases implicadas en el
proceso de relocalizacin. En nuestro caso se trabajar con ondas P.
Los algoritmos empleados para llevar cabo la correlacin cruzada son:
3.2.1. mkhdr (programa en lenguaje Fortran; autor desconocido; modificado
de Pacchiani, 2006). Escribe las cabeceras de los sismos, que contienen las
coordenadas del hipocentro, hora origen, azimut, ngulo de incidencia y tiempo de
llegada de la P,
3.2.2. mkdat (programa en lenguaje Fortran; autor desconocido; adaptado de
Pacchiani, 2006). Corta los sismogramas con una ventana de una longitud
determinada. En este trabajo se ha usado la zona de la onda P, empleando una
ventana de 512 puntos de longitud, 100 puntos por delante de la llegada de la onda P.
Se modific el programa para poder visualizar los sismogramas cortados con cualquier
programa de hoja de clculo o similar,
3.2.3. mklist (programa en lenguaje Fortran; modificado de Pacchiani, 2006).
Genera una lista con los sismos pertenecientes a un mismo multiplete,

16

3.2.4. mktbl (programa en lenguaje Fortran; Got, 1999; modificado por

Pacchiani, 2006. Adaptado de Pacchiani). Crea dos archivos con informacin de las
estaciones y de los rayos: distancia, azimut, ngulo de incidencia y tiempo de viaje de
la onda P. Estos archivos son necesarios para el clculo del retardo,
3.2.5. dbh (programa en lenguaje Fortran; modificado de Got, 1999; basado en
el algoritmo de Frechet, 1985). Calcula el retardo de los sismos dos a dos, indicando la
coherencia entre ellos. Fue modificado para que el formato de salida coincidiera con el
de entrada del programa HypoDD, usado para la relocalizacin.
3.3. Relocalizacin.
Este paso se llev a cabo mediante el programa HypoDD (Waldhauser y
Ellsworth, 2000), al que se puede acceder de manera libre en la red. Se trata de un
algoritmo que permite realizar la relocalizacin relativa de sismos, tanto a partir de
medidas de correlacin cruzada, como a partir de tiempos absolutos del tiempo de
llegada de la onda P, o mediante un procedimiento mixto. Tan slo hubo que
compilarlo y ponerlo al da a partir de la fe de errores proporcionada por los autores.
Al ser un algoritmo suficientemente conocido, no se describe en detalle en esta
memoria.

17

4. APLICACIN A LA SERIE DE 2002 DE GRGAL (ALMERA).

El mtodo descrito se va a aplicar a la serie de Grgal con el objetivo de obtener
una localizacin ms exacta de los hipocentros, mejorando las profundidades con
respecto a la solucin calculada por el Instituto Geogrfico Nacional (IGN)
4.1. Introduccin
El terremoto principal de la serie de Grgal tuvo lugar el 4 de febrero de 2002 a
las 20:09:30.48 (localizacin del IGN) (figura 6) con magnitud mbLg=5.1, y fue sentido
con intensidad V en las poblaciones de Grgal, Ohanes y Tabernas. Fue precedido por
un terremoto precursor el mismo da a las 19:55:31.82 de magnitud mbLg=2.3 y
seguido por un total de 39 rplicas (tabla 1), cuyos epicentros (IGN) se muestran en la
figura 7. En esta localizacin se observa que los epicentros se distribuyen siguiendo
una lnea en direccin N-S a lo largo de una extensin de unos 160 Km de largo, con
una anchura de unos 50 Km aproximadamente. Las profundidades determinadas por
el IGN son superficiales (menores de 40 Km), si bien de 14 sismos no se ha podido
calcular la profundidad.

Figura 6. Sismicidad de la Pennsula Ibrica y norte de frica. Se representan

sismos con M3.0, en el perodo 1980-2004. En rojo los sismos superficiales (h
40 Km), en verde los de profundidad intermedia (40 Km h 150 Km), y en azul
los profundos (600 Km h 700 Km). La estrella negra determina la localizacin
del sismo principal de la serie de Grgal (modificada de Goded, 2006)

18

Figura 7. Epicentros de los 39 sismos considerados de la serie de Grgal segn la

localizacin inicial del IGN. El crculo azul corresponde al sismo principal. El
tamao de los puntos indica la magnitud del terremoto correspondiente. El
tringulo verde indica la posicin de la estacin ACLR, que es la ms cercana a la
zona de estudio.

TABLA 1
SISMOS DE LA SERIE DE GRGAL (IGN)
FECHA
04/02/2002
04/02/2002
04/02/2002
04/02/2002
04/02/2002
04/02/2002
04/02/2002
04/02/2002
04/02/2002
04/02/2002
04/02/2002
04/02/2002
04/02/2002
04/02/2002
04/02/2002
04/02/2002
04/02/2002
04/02/2002

HORA
19:55:31.82
20:09:30.48
20:13:33.20
20:15:16.51
20:17:45.12
20:18:23.92
20:20:09.25
20:24:59.86
20:31:02.64
20:34:08.35
20:35:24.40
20:40:14.88
20:44:32.54
21:09:35.20
21:37:54.77
21:43:21.71
22:00:04.88
22:11:26.89

LATITUD
37,0972
37,0931
37,0489
37,0575
37,1093
37,0358
37,0994
37,0911
37,0916
37,0922
37,0739
37,0144
37,0463
37,0439
36,9773
37,0932
37,1016
37,0778

19

LONGITUD
-2,5347
-2,5379
-2,5001
-2,4993
-2,5366
-2,5318
-2,5251
-2,5211
-2,5852
-2,5355
-2,5230
-2,5282
-2,5326
-2,5201
-2,5444
-2,5462
-2,5457
-2,5470

H (KM)
0,0f
0,8
11,3
0,0f
4,9
0,5
10,8
0,0f
1,3
0,0f
5,9
0,0f
2,7
0,0f
6
1,1
0
0,0f

mbLg
2,3
5,1
2,7
2,2
2,4
2,3
2,1
2,5
2,6
3,0
2,5
1,8
2,4
2,0
2,3
2,9
2,8
2,8

TABLA 1
CONTINUACIN
FECHA
04/02/2002
05/02/2002
05/02/2002
05/02/2002
05/02/2002
05/02/2002
05/02/2002
05/02/2002
05/02/2002
05/02/2002
05/02/2002
05/02/2002
05/02/2002
05/02/2002
06/02/2002
06/02/2002
10/02/2002
10/02/2002
14/02/2002
15/02/2002
26/02/2002

HORA
23:56:30.18
00:03:57.83
00:19:42.46
02:21:02.69
05:28:20.50
06:08:10.96
08:40:12.57
12:41:17.88
12:51:07.12
13:51:23.39
14:21:22.77
18:06:03.43
20:00:25.99
21:42:01.22
07:06:18.40
21:21:48.08
02:42:39.67
19:28:02.48
03:16:11.43
13:17:48.00
17:19:44.43

LATITUD
37,0955
37,0970
37,0749
37,0194
37,0622
37,0386
37,1054
37,0622
37,0703
37,0395
37,0494
37,0155
37,0485
37,0723
37,0587
37,0554
37,1118
37,0940
37,1079
37,0967
37,0924

LONGITUD
-2,5447
-2,5407
-2,5021
-2,5431
-2,5417
-2,5336
-2,5471
-2,5556
-2,5446
-2,5268
-2,5331
-2,5331
-2,5351
-2,5582
-2,5458
-2,5573
-2,5403
-2,5525
-2,5558
-2,5441
-2,5448

H (KM)
2,7
1,4
10,6
11,4
10,6
16,1
6,4
10,5
3,8
4,6
0,0f
7,8
0,0f
1,2
11,2
0,0f
0,0f
7,3
0,0f
0,0f
0,0f

MBLG
1,6
2,2
2,2
1,7
2,4
2,2
2,4
3,1
2,3
2,2
2,1
2,1
1,9
3,0
2,6
2,0
1,7
2,1
1,8
2,0
2,2

4.2. Procesamiento de datos

Para la relocalizacin de esta serie por medio de formas de onda, se han
utilizado los datos correspondientes a las estaciones digitales de banda ancha y de
perodo corto del Instituto Geogrfico Nacional (IGN) y de la red conjunta del Real
Observatorio de la Armada de San Fernando y de la Universidad Complutense de
Madrid (ROA-UCM). En las tablas 2 y 3 se muestran las coordenadas y caractersticas
de las estaciones empleadas en este trabajo. Su localizacin se muestra en la figura 8.
TABLA 2
ESTACIONES IGN
Estacin
EADA
EALH
EBAN
EBER
ECOG
EALB
EMAL
EGUA
ELOJ
ELUQ
ENIJ
ERON
EQES
ETOB
EVIA

Latitud
38.1673
37.8582
38.1710
36.8978
37.2772
35.9398
36.7620
36.8337
37.1480
37.5592
36.9715
37.0180
37.8028
38.6447
38.6387

Longitud
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N

4.5772
1.4197
3.7900
2.8897
3.5663
3.0343
4.4280
3.5653
4.1530
4.2678
2.2070
3.8050
3.0712
1.5478
2.5025

W
W
W
W
W
W
W
W
W
W
W
W
W
W
W

20

Muestreo
(muestras/s)
50
50
50
50
50
50
50
50
50
50
50
50
50
50
50

Tipo
banda ancha
perodo corto
perodo corto
banda ancha
perodo corto
banda ancha
banda ancha
perodo corto
perodo corto
perodo corto
perodo corto
perodo corto
banda ancha
banda ancha
perodo corto

TABLA 3
ESTACIONES ROA-UCM

-8

Estacin

Latitud

Longitud

ACBG
ACLR
ALB
CART
CEU
CFS
ELUQ
ORGV

36.7672
37.1885
35.9300
37.5868
35.9160
35.2300
37.5605
36.8605

N
N
N
N
N
N
N
N

-6

2.1950
2.5832
3.0300
1.0120
5.3300
2.4870
4.2668
3.4293

Muestreo
(muestras/s)
50
50
50
100 / 80
100 / 80
50
50
50

W
W
W
W
W
W
W
W

-4

EADA

0
ETOB

38
EALH

358

EQES

CART

ELU y LUQ
ACLR

ECOG
ELOJ
ERON
EMAL

Gergal
EBER

ENIJ

ORGV
EGUA

ACBG

MA

36

CEU

-6

-4

ITE
ED

N
RA

EO

36

ALB y EALB

MELI

-8

ancha
ancha
ancha
ancha
ancha
ancha
ancha
ancha

EBAN

PENINSULA IBERICA

LIJA

Banda
Banda
Banda
Banda
Banda
Banda
Banda
Banda

-2
EVIA

38

Tipo

AFRICA
CFS

-2

Figura 8. Distribucin de las estaciones ssmicas de la red IGN y ROA-UCM

empleadas en este trabajo. El punto rojo indica el epicentro del sismo principal de
la serie de Grgal (IGN).
De acuerdo con la metodologa descrita, la primera etapa es la identificacin de
la fase P en los sismogramas, para lo que se han utilizado los registros en formato SAC
(Seismic Analysis Code; Goldstein, 1987) (Apndice C del CD-ROM). En algunos de los
sismogramas se poda intuir la onda S, sin embargo, no se ha tenido en cuenta en este
trabajo por dos motivos fundamentales:

al no observar la fase S claramente, el error introducido al tenerla en cuenta

poda ser ms grande que las ventajas derivadas de considerar una nueva fase,

para poder llevar a cabo la correlacin cruzada es necesario observar la misma

fase en bastantes sismogramas y, a ser posible, en las mismas estaciones para
poder llevar a cabo una comparacin satisfactoria. Con el nmero de fases S
que se poda llegar a identificar, no se cumpla este requisito.

21

Para la identificacin de la fase P se han seleccionado las estaciones ptimas

considerando como tales las estaciones ms cercanas, hasta un mximo de 214 Km,
que son las que proporcionan un buen registro de la seal, con una relacin
seal/ruido baja. Especial inters tienen las estaciones ACLR, de la red ROA-UCM,
situada en Calar Alto (a 10 Km del epicentro del sismo principal) y las estaciones
EBER y ENIJ del IGN, que se encuentran a menos de 40 Km (figura 8).
A partir del tiempo de llegada de las ondas P se ha realizado una primera
localizacin con el programa Hypo71 (Apndice D del CD-ROM).
El siguiente paso consisti en remuestrear las seales ssmicas con el mismo
t, ya que la razn de muestreo oscila entre 50 y 100 muestras/s segn la estacin
(tablas 2 y 3). Es un requisito necesario para comparar unos sismogramas con otros,
tanto en el paso para obtener las familias de sismos, como para calcular la doble
diferencia. El muestreo elegido en este trabajo ha sido 50 muestras/s, es decir
t=0.02s. Tanto la identificacin de la onda P como el remuestreo se han realizado con
SAC.
Es importante comprobar que los espectros de los sismogramas remuestreados
no varan respecto a los sismogramas originales, para asegurarse de que el contenido
en frecuencias sigue siendo el mismo despus del remuestreo. En las figuras 9a y 9b
se puede ver un ejemplo, donde se puede apreciar cmo la resolucin ha disminuido
ligeramente al reducir el ritmo de muestreo, observndose menos detalles en la nueva
seal ssmica. Sin embargo, los espectros de la seal original y remuestreada son
idnticos, permaneciendo inalterados los valores de amplitud en su parte plana (2101
cuentas) y la frecuencia esquina (7 Hz).

Figura 9 a. Sismograma original y remuestreado a 50 muestras/s del sismo del 5 de

febrero de 2002 a las 05:28:20.50. Los crculos rojos indican dos instantes de
tiempo para los cuales se aprecia, a modo de ejemplo, una significativa prdida de
resolucin.

22

AMPLITUD

AMPLITUD

FRECUENCIA (Hz)

FRECUENCIA (Hz)

Figura 9 b. Espectros de amplitud del sismograma original y del remuestreado.

En esta fase del tratamiento de datos es necesario llevar a cabo un filtrado de
los registros, pues muchos de ellos llevan un ruido aleatorio asociado que puede
interferir a la hora de proceder a la correlacin cruzada. Asimismo, la correlacin
cruzada es ms fiable si se lleva a cabo un filtrado entre 1 y 10 Hz (Hauksson y
Shearer, 2005). De modo que se procedi a aplicar un filtro Butterworth paso-banda
con 3 polos, entre 1 y 10 Hz. Con este paso queda finalizado el tratamiento de datos y
se puede comenzar con el proceso de extraccin de familias.
4.3. Resultados
El siguiente paso consiste en obtener los sonogramas y normalizarlos, para
finalmente obtener las familias de sismos. Empleando un valor umbral de 6105 en la
distancia de Levenstein se han obtenido dos familias, a partir de los 39 sismos
iniciales de partida. La primera de ellas, la familia A, consta de 7 sismos, y la familia B
de 6 (tabla 4). Con los sismos restantes no se ha podido llevar a cabo la relocalizacin,
ya que han sido descartados debido a que, con la distancia de Levenstein elegida como
umbral, no formaban parte de ninguna familia. Se pueden consultar los diagramas de
rbol obtenidos en el Apndice E del CD-ROM, y el estudio de la influencia del nmero
de estaciones y de la distancia de Levenstein1 en la extraccin de familias en el
Apndice F.

ID
sismo
1
2
3
4
5
6
7

04
04
04
04
04
05
10

TABLA 4
FAMILIAS DE SISMOS
Familia A
Familia
Fecha
Hora
Fecha
/ 02 /2002
20:13:32
04 / 02 /2002
/ 02 /2002
20:18:22
04 / 02 /2002
/ 02 /2002
20:34:07
05 / 02 /2002
/ 02 /2002
21:43:21
05 / 02 /2002
/ 02 /2002
22:00:03
05 / 02 /2002
/ 02 /2002
21:42:00
05 / 02 /2002
/ 02 /2002
02:42:39

Apartado 3.1.5., p. 15.

23

B
Hora
20:44:31
21:37:53
05:28:20
06:08:10
08:40:11
12:51:06

En la figura 10 se muestran dos ejemplos de sismos de una misma familia, uno

para cada caso, y en dos estaciones distintas, ACLR y EBER. Se puede apreciar que el
contenido en frecuencias es prcticamente igual. La diferencia reside en la amplitud
de las fases, debido a que se trata de sismos de diferente magnitud. Para la familia A
se han seleccionado el sismo del 4 de febrero a las 20:18, que tuvo una magnitud
mbLg=2.3, y el del 10 de febrero a las 02:42, con mbLg=1.7. La diferencia en magnitud
es lo suficientemente grande como para apreciar diferencias en amplitud en los
sismogramas. En la familia B se han seleccionado sismos del 5 de febrero, a las 05:28
y 06:08, con magnitud mbLg de 2.4 y 2.2 respectivamente. En este caso la diferencia
de magnitud es menor, por lo que las amplitudes son ms similares.

Figura 10. Ejemplos de dos sismos pertenecientes a la familia A (en la parte

superior) y B (en la parte inferior).
Una vez obtenidas las diferentes familias de terremotos, se procede al clculo
de la doble diferencia para todos los pares de sismos de una misma familia (resultados
intermedios en los Apndices G, H e I del CD-ROM). Los resultados obtenidos de la
doble diferencia se muestran en las tablas 5 y 6, en las que se indica la identificacin
de los sismos que se van a correlacionar dentro del multiplete (ID1 e ID2), y el nmero
de estaciones en las que se ha llevado a cabo el clculo del retardo. A continuacin se
indica la estacin, el retardo (10-4 s) para esa estacin y el peso, definido segn la
ecuacin [10] que se le asigna a ese valor. Por ejemplo, en la primera lnea de la tabla,
ID1=1 e ID2=2 indican que se van a correlacionar los sismos 1 y 2 en seis estaciones
diferentes. En el segundo caso, ID1=1 e ID2=3, sealan una correlacin entre los
sismos 1 y 3 en 10 estaciones distintas. Y as sucesivamente.

24

TABLA 5
RETARDOS DE LA FAMILIA A
ID1
Sta
1
ACBG
EGUA
1
ACBG
EBER
ENIJ
EVIA
1
ACBG
EBER
ELUQ
EVIA
1
ACBG
EBER
ENIJ
EVIA
1
ACBG
EBER
ENIJ
1
ACLR
EGUA
2
ACBG
EGUA
2
ACBG
EGUA
2
ACBG
EGUA
2
ACBG
EGUA
2
ACLR
ENIJ
3
ACBG
EBAN
ELUQ
EVIA
3
ACBG
EBAN
ENIJ
EVIA
3
ACBG
EBAN
ENIJ

ID2
retardo
10-4 s
2
1945.7
2112.6
3
657.2
2203.1
2003.6
4381.9
4
2177.3
3967.3
1980.9
975.0
5
5134.3
5482.7
5189.5
1581.8
6
-368.6
5988.3
4591.6
7
2478.3
2427.7
3
-1588.6
4198.5
4
184.0
3728.3
5
3201.1
3585.1
6
1822.7
2071.8
7
86.9
-83.6
4
1316.1
-3522.3
-714.3
-358.8
5
2438.0
2204.3
2232.7
-3343.4
6
2898.4
3853.4
2069.9

n de estaciones en
Peso
Sta
w
6
0.8291
ACLR
0.8689
ENIJ
10
0.7723
ACLR
0.7833
EGUA
0.5865
ORGV
0.647
10
0.6850
ACLR
0.8887
ECOG
0.6485
ENIJ
0.5335
10
0.8390
ACLR
0.9357
ECOG
0.9280
ORGV
0.7438
9
0.6382
ACLR
0.8693
ECOG
0.7980
ORGV
5
0.9105
EBER
0.9238
ENIJ
6
0.8618
ACLR
0.7796
ENIJ
6
0.8280
ACLR
0.8216
ENIJ
6
0.9206
ACLR
0.7770
ENIJ
6
0.8385
ACLR
0.7418
ENIJ
4
0.9687
EBER
0.9492
10
0.7965
ACLR
0.8584
EBER
0.6672
ENIJ
0.8081
10
0.8437
ACLR
0.9136
EBER
0.7821
ORGV
0.7257
8
0.9278
ACLR
0.7316
EBER
0.8156
ORGV

que se lleva a cabo la correlacin cruzada.

retardo
Peso
Sta
retardo
10-4 s
10-4 s

Peso

2389.1
4325.1

0.9678
0.8847

EBER
ORGV

1989.4
429.3

0.9599
0.7085

3110.4
2478.7
-3288.7

0.7195
0.8235
0.8138

EBAN
ELUQ
EQES

3341.9
-463.7
-289.2

0.8138
0.6420
0.7045

0
9344.8
4528.5

0
0.6602
0.6831

EBAN
EGUA
ORGV

9036.9
2511.9
-2356.9

0.8312
0.7970
0.8623

5496.1
7328.8
-1316.8

0.9038
0.7503
0.8189

EBAN
EGUA
EQES

5314.8
5694.0
3524.6

0.9389
0.9001
0.7068

3621.3
-291.6
349.2

0.8164
0.6461
0.7025

EBAN
EGUA
ERON

9703.8
412.2
8366.7

0.6126
0.7969
0.7988

255.7
2154.8

0.8800
0.9335

ECOG

835.4

0.6620

816.6
808.1

0.8239
0.6498

EBER
ORGV

231.2
131.7

0.8520
0.7487

-729.6
709.9

0.8270
0.8318

EBER
ORGV

1974.0
869.0

0.9328
0.7303

3146.2
2955.7

0.9517
0.9853

EBER
ORGV

3521.8
3145.7

0.9808
0.8362

1216.6
2798.2

0.8183
0.8997

EBER
ORGV

3970.7
3730.4

0.9364
0.8030

155.3

0.8883

EGUA

1863.9

0.7654

513.2
1914.7
2744.5

0.9017
0.8834
0.8665

EALH
EGUA
ORGV

-951.0
1354.2
1557.3

0.9720
0.9595
0.8107

0
3322.0
3046.8

0
0.8853
0.9252

EALH
EGUA
EQES

3205.2
3256.2
4139.9

0.9590
0.9451
0.8153

3638.8
3649.7
-152.4

0.9382
0.9638
0.7778

EALH
EGUA

3701.3
3597.7

0.9859
0.9464

25

TABLA 5
(CONTINUACIN)
3
ACBG
ENIJ
4
ACBG
EALH
CART
ENIJ
4
ACBG
EBAN
ECOG
ORGV
4
ACLR
EGUA
5
ACBG
EBAN
ECOG
ORGV
5
ACLR
EGUA
6
ACLR
EGUA

7
2356.5
-818.9
5
879.3
4380.4
985.4
958.9
6
3411.6
-199.1
76.5
360.4
7
-864.4
-194.0
6
485.3
2180.0
612.1
657.3
7
-3064.0
0
7
-1435.6
-557.5

4
0.8319
0.7735
12
0.9010
0.9899
0.7356
0.8624
10
0.8074
0.9053
0.8445
0.7865
5
0.7865
0.8244
10
0.8816
0.8308
0.9825
0.9432
5
0.9804
0
5
0.7699
0.7496

EBER

-10.2

0.7255

EGUA

-23.6

0.8832

ACLR
EBAN
ECOG
ORGV

543.4
-4451.4
1498.0
46.0

0.7351
0.9330
0.8856
0.8449

EADA
EBER
EGUA
EVIA

2897.6
1519.2
1909.3
539.3

0.7352
0.9615
0.9450
0.6681

ACLR
EBER
EGUA

-321.9
1928.1
0

0.7660
0.9688
0

EALH
CART
ENIJ

4622.8
4404.5
2124.7

0.9696
0.7306
0.8543

EBER
ENIJ

-1826.7
-793.8

0.8356
0.8093

ECOG

-1662.9

0.6702

ACLR
EBER
EGUA

-1817.3
411.3
411.6

0.9048
0.9955
0.9639

EALH
CART
ENIJ

444.5
-1701.8
-176.7

0.9289
0.8780
0.9457

EBER
ENIJ

-3365.5
-3045.0

0.8944
0.9557

ECOG

-3137.8

0.7649

EBER
ENIJ

-6928.6
-2861.2

0.8126
0.9075

ECOG

-3771.5

0.7257

TABLA 6
RETARDOS DE LA FAMILIA B
ID1
Sta
1
ACBG
EBER
ENIJ
1
ACBG
EBER
ENIJ
1
ACBG
EBER
ORGV
1
ACBG
EGUA
1
ACLR
EGUA
2
ACBG
EBER
ENIJ
2
ACBG
EBER
ORGV

ID2
retardo
10-4 s
2
-3841.3
-2154.8
1179.9
3
7044.1
-1263.6
1610.7
4
-325.3
-1123.7
-1113.6
5
-1361.4
-327.0
6
-955.4
-281.3
3
-1000.8
886.6
298.8
4
907.2
1046.1
831.1

n de estaciones en que se lleva

Peso
Sta
retardo
10-4 s
9
0.8192
ACLR
-2927.4
0.9879
ECOG
-1603.3
0.9079
ORGV
1015.5
8
0.7334
ACLR
-4065.7
0.9813
ECOG
78.7
0.9173
ORGV
-1351.7
7
0.6993
ACLR
-3725.7
0.9568
EGUA
-1208.1
0.8492
6
0.8050
ACLR
1124.9
0.7436
ENIJ
5017.9
6
0.8514
EBAN
619.7
0.8128
ENIJ
1628.0
8
0.7980
ACLR
-1109.0
0.9854
ECOG
-33.4
0.9694
ORGV
662.8
7
0.7694
ACLR
-778.2
0.9660
EGUA
775.6
0.8489

26

a cabo la correlacin cruzada.

Peso
Sta
retardo
10-4 s

Peso

0.8560
0.7955
0.7991

EBAN
EGUA
EQES

-165.3
-2352.0
1333.7

0.7934
0.7838
0.6909

0.8653
0.8859
0.8197

EBAN
EGUA

-1038.3
-1511.7

0.7939
0.9146

0.8742
0.8773

EBAN
ENIJ

0
1829.9

0
0.9153

0.8134
0.9231

EBER
ORGV

1863.6
1899.5

0.9338
0.7911

0.7104
0.6862

EBER
ORGV

1931.5
1903.8

0.9864
0.6903

0.9626
0.8521
0.7733

EBAN
EGUA

-961.9
565.1

0.8985
0.8957

0.9712
0.8446

EBAN
ENIJ

-664.1
568.0

0.8399
0.9856

TABLA 6
(CONTINUACIN)
2
ACBG
EGUA
2
ACLR
EGUA
3
ACBG
EBER
ORGV
3
ACBG
EGUA
3
ACLR
EGUA
ERON
4
ACBG
EGUA
4
ACLR
EGUA
EVIA
5
ACLR
ENIJ

5
3875.8
3708.4
6
1934.1
1817.1
4
2026.4
160.9
1884.1
5
3056.9
3169.5
6
2928.7
0
9387.2
5
1255.7
3053.7
6
2604.2
2844.5
-1978.8
6
0
681.8

6
0.8729
0.9179
6
0.9884
0.7906
8
0.8771
0.9695
0.7802
6
0.7356
0.9637
8
0.9598
0
0.6020
6
0.6615
0.8441
7
0.9711
0.7854
0.7198
5
0
0.7405

ACLR
ENIJ

4081.1
3799.7

0.8236
0.9643

EBER
ORGV

4026.7
3868.4

0.9317
0.8631

EBAN
ENIJ

5066.8
482.6

0.7274
0.7461

EBER
ORGV

3945.9
-1116.9

0.9624
0.7839

ACLR
EGUA
EVIA

336.2
159.7
-4123.0

0.9964
0.8908
0.6388

EBAN
ENIJ

268.8
270.3

0.9268
0.9729

ACLR
ENIJ

3852.8
3457.6

0.9211
0.9652

EBER
ORGV

3160.6
-372.1

0.9506
0.6873

EBAN
ENIJ
EVIA

2993.2
4024.5
-1194.1

0.7105
0.7062
0.6053

EBER
ORGV

3116.0
-541.2

0.9713
0.7814

ACLR
ENIJ

3500.5
3209.0

0.8994
0.9688

EBER
ORGV

2979.3
3039.9

0.9472
0.8624

EBAN
ENIJ

2608.8
-80.7

0.7320
0.7392

EBER
ORGV

2974.1
-601.9

0.9553
0.7202

EBER
ORGV

-13.0
-1846.6

0.9779
0.8392

EGUA

-34.3

0.7563

Un valor nulo del retardo en una estacin indica que en esa estacin los
registros han sido lo suficientemente distintos como para no poder comparar la seal,
y por tanto se les asocia un peso nulo.
Se puede apreciar que el peso no depende nicamente de la distancia entre la
estacin y el hipocentro. Si as fuera, las estaciones ms alejadas tendran asociado
un peso menor. Pero, tal y como se explic en la metodologa, el peso depende del
ndice de coherencia entre las dos seales que se correlacionan, dando mayor peso a
seales ms coherentes entre s. Por tanto se observa que, en algunos casos, a
retardos calculados en estaciones ms lejanas se les ha asociado un peso mayor.
No obstante, en trminos generales, s se observa que en las estaciones ms
cercanas se obtienen retardos con un mayor peso asociado, como es el caso de ACLR,
ENIJ o EBER, y retardos con pesos menores en las ms lejanas, particularmente en
EVIA. La explicacin reside en que, para estaciones ms cercanas las seales que se
registran generalmente son mejores, con una mejor relacin seal/ruido, lo que
conlleva que la correlacin cruzada sea mejor, y por tanto, el peso tambin mayor.
La excepcin es la estacin EALH, situada a unos 130 Km del epicentro del
baricentro de la familia A, que proporciona unos pesos entre 0.9 y 1.0 (frente a pesos

27

entre 0.6 y 0.7 en estaciones ms cercanas), si bien la correlacin solo ha sido posible
para 5 pares de sismos, por lo que no se considera significativo este resultado. La
estacin que proporciona en general los peores retardos es ELUQ, con pesos muy
bajos, que varan entre 0.6 0.7. La estacin que proporciona el peso menor de todos
es EVIA, que es la estacin ms alejada del epicentro del sismo principal.
Con estos resultados se procede a llevar a cabo la relocalizacin mediante el
programa HypoDD para cada una de las familias por separado. El modelo de Tierra
empleado fue obtenido por Daobeitia et al. (1998) (figura 11). Es un modelo especfico
para el rea de estudio, frente al modelo empleado por el IGN, de carcter ms
general. Se ha elegido este modelo de Tierra a partir de los resultados obtenidos por
Rodrguez Abad (2004), que concluye que este modelo es ms adecuado para la
relocalizacin del sismo principal de esta serie.

Velocidad
(Km/s)
5.70
5.75
5.90
5.95
6.70
6.75
7.90

Profundidad
(Km)
0.00
8.00
12.00
18.00
23.00
29.00
35.00

Figura 11. Modelo de Tierra empleado en este trabajo. (Daobeitia, 1998).

Los resultados de la relocalizacin se muestran en las tablas 7, 8, 9 y 10,
comparados con la localizacin inicial dada por el IGN.

TABLA 7
RELOCALIZACIN FINAL FAMILIA A.
ID
1
2
3
4
5
6
7

Fecha
04
04
04
04
04
05
10

/
/
/
/
/
/
/

02
02
02
02
02
02
02

/2002
/2002
/2002
/2002
/2002
/2002
/2002

Hora origen
(hh:mm:ss.ss)
20:13:32.88
20:18:23.00
20:34:07.84
21:43:21.00
22:00:03.60
21:42:00.08
02:42:39.18

RMS

Latitud

Longitud

0.70
0.24
0.47
0.56
0.49
0.50
0.23

37.0903
37.0944
37.0951
37.0941
37.0944
37.0929
37.1009

-2.5656
-2.5580
-2.5622
-2.5661
-2.5611
-2.5619
-2.5600

28

Profundidad
(Km)
9.1 1.4
9.0 4.5
10.0 1.3
9.3 1.3
8.8 1.3
10.0 1.3
10.4 2.2

TABLA 8
LOCALIZACIN IGN FAMILIA A.
ID
1
2
3
4
5
6
7

Fecha
04
04
04
04
04
05
10

/
/
/
/
/
/
/

02
02
02
02
02
02
02

/2002
/2002
/2002
/2002
/2002
/2002
/2002

Hora origen
(hh:mm:ss.ss)
20:13:33.20
20:18:23.92
20:34:08.35
21:42:21.71
22:00:04.88
21:42:01.22
02:42:39.67

RMS

Latitud

Longitud

0.65
0.73
0.75
0.97
0.68
0.70
0.40

37.0489
37.0358
37.0922
37.0932
37.1016
37.0723
37.1118

-2.5001
-2.5318
-2.5355
-2.5462
-2.5457
-2.5582
-2.5403

Profundidad
(Km)
11.3 0.5
0.5 2.5
0.0f 0.0
1.1 1.1
0.0 1.3
1.2 1.0
0.0f 3.7

TABLA 9
RELOCALIZACIN FINAL FAMILIA B.
ID
1
2
3
4
5
6

Fecha
04
04
05
05
05
05

/
/
/
/
/
/

02
02
02
02
02
02

/2002
/2002
/2002
/2002
/2002
/2002

Hora origen
(hh:mm:ss.ss)
20:44:32.02
21:37:54.00
05:28:20.33
06:08:10.76
08:40:11.86
12:51:06.21

RMS

Latitud

0.44
0.36
0.43
0.36
0.38
0.39

37.0742
37.0722
37.0659
37.0646
37.0660
37.0670

Longitud
-2.5723
-2.5687
-2.5682
-2.5680
-2.5682
-2.5690

Profundidad
(Km)
6.7 2.1
7.6 1.8
8.1 1.8
7.7 1.8
5.6 7.0
7.3 1.9

TABLA 10
LOCALIZACIN IGN FAMILIA B.
ID
1
2
3
4
5
6

Fecha
04
04
05
05
05
05

/
/
/
/
/
/

02
02
02
02
02
02

/2002
/2002
/2002
/2002
/2002
/2002

Hora origen
(hh:mm:ss.ss)
20:44:32.54
21:37:54.77
05:28:20.50
06:08:10.96
08:40:12.57
12:51:07.12

RMS

Latitud

Longitud

0.97
0.76
0.75
0.52
0.60
0.53

37.0463
36.9773
37.0622
37.0386
37.1054
37.0703

-2.5326
-2.5444
-2.5417
-2.5336
-2.5471
-2.5446

Profundidad
(Km)
2.7 2.5
6.0 2.2
10.6 0.3
16.1 3.5
6.4 2.7
3.8 2.3

El nmero de correlaciones cruzadas de las que se dispone para cada sismo

vara en cada caso (tabla 11).

TABLA 11
NMERO DE CORRELACIONES DE CADA SISMO
Familia A
Sismo
N de
correlaciones
1
49
2
34
3
48
4
53
5
53
6
48
7
28

Sismo
1
2
3
4
5
6

29

Familia B
N de
correlaciones
36
36
38
35
29
32
36

A partir de los resultados mostrados en las tablas 7, 8, 9, 10 y 11 se concluye

que:
1. los valores en latitud y longitud son similares a los obtenidos por el IGN. No
ocurre lo mismo con el valor de la profundidad. Para la familia A la
profundidad de todos los terremotos est comprendida entre 8.8 y 10.4 Km,
mientras que las profundidades del IGN tienen un mayor rango de variacin
(entre 1.1 y 11.3 Km) o no se ha podido obtener la profundidad (valores de 0.0f
Km). En la familia B tambin se observa una clara mejora en la profundidad
frente a las soluciones del IGN: todos los sismos de nuestra solucin tienen
profundidades comprendidas entre 5.6 y 8.1 Km, frente al rango de 2.7 a 16.1
Km del IGN,
2. los errores obtenidos en profundidad en la familia A son ligeramente mayores
para el sismo 1 (1.4 Km frente a 0.5 Km del IGN) y similares en los sismos 4 y
6 (1.3 Km frente a 1.1 Km del IGN). En el resto de los sismos no es posible
hacer esta comparacin ya que para esos terremotos el IGN no determin la
profundidad. En la familia B los errores son menores salvo en los sismos 3 y 5,
en que nuestros errores son sensiblemente mayores,
3. el valor cuadrtico medio (RMS) de la hora origen en la familia A oscila entre
(0.70 0.23) frente a (0.97 0.40) del IGN, lo que indica una clara disminucin
salvo para el sismo 1. En el caso de la familia B, se encuentra que el RMS
oscila entre (0.36 - 0.44) frente a (0.52 - 0.97) de la solucin del IGN,
4. los sismos para los que se dispone de un mayor nmero de retardos llevan
asociado finalmente un error menor, de manera general.
En la figura 12 se muestran los epicentros de la familia A (azul) obtenidos en
este trabajo y los calculados por el IGN (verde). Se observa que, tras la relocalizacin,
los epicentros se han concentrado agrupndose en una rea reducida, siguiendo una
lnea de direccin de unos 52 con el norte geogrfico y abarcando una extensin de
unos 2 Km. La localizacin del IGN tambin mostraba una distribucin NE-SW para la
mayora de los sismos pero en una rea mucho ms extensa. Las mayores diferencias
corresponden a los sismos 1 y 2, muy alejados en la distribucin del IGN, mientras
que en la obtenida en este trabajo aparecen muy prximos a los restantes sismos.
La figura 13 muestra la distribucin de hipocentros en profundidad segn un
corte en direccin 52 N. Se observa que los focos se distribuyen entre los 8 y 10 Km
de profundidad, siguiendo una lnea con un buzamiento de 37 hacia el noreste.

30

Los epicentros de los sismos relocalizados de la familia B (azul) y los calculados

por el IGN (verde) se muestran en la figura 14. Se observa una distribucin de
epicentros prcticamente N-S similar a la obtenida por el IGN, si bien los terremotos
abarcan una extensin de unos 2 Km en este trabajo, frente a los cerca de 14 Km en el
caso del IGN.
En la figura 15 se muestra la distribucin en profundidad de los focos de la
familia B, para un corte N-S. Salvo el sismo 5 (con profundidad 5.6 Km), los dems se
distribuyen a lo largo de una lnea con un buzamiento de 56 hacia el sur. No
obstante, este sismo es el que tiene un mayor error asociado, por lo que, si se tiene en
cuenta este error, se puede considerar que el sismo sigue la tendencia general de un
buzamiento de 56 hacia el sur.
-2.60

-2.55

-2.50

-2.45

-2.40

A Gergal
7

37.10

7
5 2
3
4
6
1

A'

37.10

Principal

37.05

37.05

37.00

37.00

36.95
-2.60

-2.55

-2.50

-2.45

36.95
-2.40

Figura 12. Familia A. En verde se representan los epicentros de los sismos a partir
de las coordenadas del IGN. En azul estn representados los epicentros tras la
relocalizacin con correlacin cruzada. La estrella indica la localizacin del sismo
principal, segn la localizacin del IGN.
Familia A

A'

Profundidad (km)

2
4
6
8
1
10
12
12

5
4
6
3

13

2
7
14

15

16

17

Figura 13. Corte en profundidad en direccin (37.07N, 2.06W) (37.12,2.52W)

de los hipocentros de la familia A.

31

-2.60

-2.55

-2.50

-2.45

-2.40

Gergal
A

37.10

37.10
Principal

2
6
5
4

6
3

37.05

37.05
1
4

37.00

37.00

36.95
-2.60

-2.55

-2.50

36.95
-2.40

-2.45

Figura 14. Familia B. En verde se representan los epicentros de los sismos a partir
de las coordenadas del IGN. En azul estn representados los epicentros tras la
relocalizacin con correlacin cruzada. La estrella indica la localizacin del sismo
principal, segn la localizacin del IGN.
Familia B

A'

Profundidad (km)

2
4
5

1
2

6
3

10
12

12

13

14

15

16

17

Figura 15. Corte en profundidad en direccin norte-sur de los hipocentros de la

familia B.
Cesca (2005) obtuvo para el sismo principal de Grgal un mecanismo de falla
normal con un plano de falla orientado en direccin NW-SE, con buzamiento hacia el
NE y el foco a 9 Km de profundidad. El otro plano, orientado en direccin NE-SW, es
prcticamente vertical. Por tanto no es posible relacionar los alineamientos de las
familias A y B con el mecanismo del sismo principal. No obstante, hay que recordar
que el sismo principal y sus rplicas no tienen por qu estar asociados a la misma
falla, pudiendo ocurrir que el sismo principal desencadene el movimiento de fallas
secundarias en las que se producen las rplicas.

32

Esta metodologa ha sido aplicada tambin a la serie de rplicas del sismo de

Alhucemas (Marruecos) del 24 de febrero de 2004 (mb=6.2). Sin embargo, para estos
sismos, pese a que las rplicas son de mayor magnitud (4.0 5.0), no se han obtenido
resultados satisfactorios ya que la estacin ms cercana se encontraba a 40 Km de
distancia, estando las restantes a distancias entre 100 y 500 Km. Ya se ha visto que el
mtodo es muy sensible a poseer registros a distancias epicentrales pequeas, lo que
permite disponer de buenas formas de onda y en consecuencia poder realizar una
buena correlacin.

33

5. CONCLUSIONES
La relocalizacin de terremotos por medio de formas de onda presenta claras
ventajas frente a los mtodos tradicionales, ya que
a) permite disminuir los errores tanto de la profundidad como en la hora origen,
al no utilizar los tiempos de llegada absolutos de la onda P,
b) la menor influencia del modelo de Tierra permite una mejor determinacin de
la profundidad. Debido a la cercana entre s de los hipocentros, los rayos
recorrern trayectos similares hasta la estacin, por lo que el modelo de Tierra
no tiene una influencia directa sobre la relocalizacin.
Los resultados de la aplicacin de esta metodologa a la serie de Grgal
(Almera) de 2002, formada por 39 sismos, ha permitido concluir que
1. es de gran importancia que los registros de cada sismo se hayan producido en
las mismas estaciones para poder llevar a cabo el mayor nmero de
correlaciones

posible,

con

el

fin

de

lograr

una

mejor

relocalizacin

especialmente en profundidad,
2. se han obtenido dos familias de terremotos formadas por 7 sismos (familia A) y
6 sismos (familia B) de entre los 39 sismos de la serie ssmica,
3. la familia A presenta una distribucin de epicentros alineados segn una
direccin N52E, con un extensin de unos 2 Km. Los focos se distribuyen en
profundidad entre 8 10 Km, siguiendo una lnea de buzamiento de 37 hacia
el noreste,
4. por su parte, los epicentros de la familia B se distribuyen segn una lnea con
orientacin N-S con una extensin de aproximadamente 1.5 Km. En
profundidad, los hipocentros se sitan entre los 5.6 Km y los 8 Km segn una
distribucin con un buzamiento de 56 hacia al sur.
5. las soluciones de ambas familias presentan menores errores en profundidad y
RMS frente a los valores calculados por el IGN. En particular, se ha logrado
determinar la profundidad de varios sismos que el IGN no haba determinado,

34

6. ninguna de las familias puede ser asociada al mecanismo focal propuesto por
Cesca (2005) para el sismo principal de la serie de Grgal,
7. esta metodologa permite obtener detalles de la ocurrencia de series de
rplicas, posibilitando el estudio ms detallado del proceso de generacin y
ocurrencia de sismos, y
8. se ha empleado la transformada wavelet, frente a la transformada de Fourier,
como viene siendo habitual, lo que permite un anlisis simultneo en el
dominio del tiempo y de las frecuencias.
Como futura lnea de trabajo se propone:
1. la utilizacin de la onda S en la correlacin cruzada, para obtener una
relocalizacin ms precisa al aumentar el nmero de correlaciones disponibles
para cada par de sismos, y
2. una nueva lnea de investigacin, abierta por este trabajo, que va a permitir
estudiar con mayor profundidad la relocalizacin tanto de series ms
numerosas de sismos, como enjambres de terremotos, con el fin de conocer las
estructuras sismognicas de diversas zonas de estudio, como es el caso
particular del mar de Alborn.

35

6. BIBLIOGRAFA

Cesca, S. (2005). Inversin del tensor momento ssmico de terremotos

superficiales a distancias regionales. Tesis doctoral. Universidad Complutense
de Madrid. Madrid. 289 pp.

Daobeitia, J.J., V. Sallares y J. Gallart. (1998). Local earthquakes seismic

tomography in the Betic Cordillera (southern Spain). Tectonophysics, 160,
225-239.

Deichmann, N. y M. Garca-Fernndez. (1992). Rupture geometry from highprecision relative hypocentre locations of microearthquake clusters. Geophys.
J. Int., 110 (3), 501-517.

Frechet, J. (1985). Sismogense et doublets sismiques. Tesis doctoral.

Universit Scientifique et Mdicale de Grenoble. 206 pp.

Gaucher, E. (1998). Comportement Hydromchanique d'un Massif Fractur:

Apport de la Microsismicit Induite. PhD thesis, Universit de Paris VII.

Geiger, L.: Herdbestimmung bei Erdbeben aus den Ankunftszeiten. Nachrichten

von der Kniglichen Gesellschaft der Wissenschaften zu Gttingen 4 (1910),
331-349.

Geller, R. J. and S. Mueller. (1980). Four similar earthquakes in central

California. Geophys. Res. Lett., 7(10), 821-824.

Goded Milln, T. (2006). Reevaluacin de daos de los terremotos de Mlaga de

1494 y 1680. Trabajo de investigacin para la obtencin del Diploma de
Estudios Avanzados. Universidad Complutense de Madrid. 152 pp.

Goldstein, P. (1987). Seismic Analysis Code (SAC). Universidad de California.

Got, J.-L., J. Frechet, and F. Klein. (1994). Deep fault plane geometry inferred
from multiplet relative relocation beneath the south flank of Kilauea. J.
Geophys. Res., 99,15.375-15.386.

Goupillaud, P., A. Grossman, y J. Morlet. (1984). Cycle-Octave and Related

Transforms in Seismic Signal Analysis. Geoexploration, 23:85-102.

36

Haukkson, E. y P. Shearer (2005). Southern California hypocenter relocation

with waveform cross-correlation, part 1: results using the double-difference
method. Bull. Seism. Soc. Am. 95, 3, 896903.

Ito, A. (1985). High resolution relative hypocenters of similar earthquakes by

cross-spectral analysis method. J. Phys. Earth, 33, 279294.

Ito, A. (1990). Earthquake swarm activity revealed from high-resolution

relative hypocenters clustering of microearthquakes. Tectonophysics, 175,
47-66.

Klein,

F.W.

(1978).

Hypocenter

location

program

HYPOINVERSE,

U.S.

Geological Survey. Open-File Report 78-694, 113 pp.

Lee, W. H. K. y J. C. Lahr (1972). HYP071: A computer program for

determining hypocenter, magnitude, and first motion pattern of local
earthquakes, Open File Report, U. S. Geological Survey, 100 pp.

Levenstein, V. I. (1965) Binary codes with use of deletions, insertions and

substitutions of symbols. Dokl. A. N. SSSR, 163,845-848.

Lyon-Caen, H., P. Papadimitriou, A. Deschamps, P. Bernard, K. Makropoulos,

F. Pacchiani y G. Patau (2004). First results of the CRLN seismic array in the
western Corinth Gulf: Evidence for old fault reactivation. C. R. Geoscience,
336:343-351.

Menke, W. (1999). Using waveform similarity to constrain earthquake

locations. Bull. Seism. Soc. Am. 89, 1143-1146.

Pacchiani, Francesco. (2006). tude sismologique des failles normales actives

du rift de Corinthe. Thse de doctorat. Universit Paris XI. 383 pp.

Paige, C.C. y M.A. Saunders (1982). LSQR: Sparse linear equations and least
squares problems. ACM Transactions on mathematical software 8/2, 195-209.

Pavlis, G.L. (1986) Appraising earthquake hypocenter location errors: a

complete, practical approach for single-event locations. Bull. Seism. Soc. Am.,
76, 1699-1717.

37

Polikar, Robi. (1994). The wavelet tutorial.

http://www-ccrma.stanford.edu/~unjung/mylec/WTpart1.html
http://www-ccrma.stanford.edu/~unjung/mylec/WTpart2.html
http://www-ccrma.stanford.edu/~unjung/mylec/WTpart3.html

Poupinet, G., W. L. Ellsworth, y J. Frechet (1984). Monitoring velocity

variations in the crust using earthquake doublets: An application to the
Calaveras Fault, California. J. Geophys. Res., 89(B7), 5719-5731.

Rodrguez Abad, I. (2004) Estudio de la determinacin hipocentral e influencia

de los errores. Trabajo de investigacin para la obtencin del Diploma de
Estudios Avanzados. Universidad Complutense de Madrid. 133 pp.

Shearer, P., E. Hauksson y G. Lin (2005). Southern California Hypocenter

Relocation with Waveform Cross-Correlation, Part 2: Results Using SourceSpecific Station Terms and Cluster Analysis. Bull. Seism. Soc. Am., 95, 904915.

Waldhauser, F. y W. L. Ellsworth. (2000). A double-difference earthquake

location algorithm: method and application to the Northern Hayward Fault.
Bull. Seismol. Soc. Am., 90(6): 1353-1368.

Waldhauser, F. y W. L. Ellsworth. (2002). Fault structure and mechanics of

the Hayward Fault, California, from double-difference earthquake locations. J.
Geophys. Res., 107(B3), doi: 10.1029/2000JB000084.

Zhizhin, M., A. Gvishiani, S. Bottard, B. Mohammadioun y J. Bonnin (1992).

Classification of strong motion waveforms from different geological regions
using syntactic pattern recognition scheme. Cah. C. Europ. Geodyn. Sismol., 6,
33-42.

38

7. ANEXO 1

Cuentas digitales

Ampliacin de la figura 4.

t (s)

39