Estudio de la fuente sísmica con estaciones en tierra y mar Study of the seismic source using onshore and ocean stations

Resumen   Abstract   Índice   Conclusiones


Cabieces Díaz, Roberto

2023-A
Descargar PDF  



Resumen

Los parámetros focales de los terremotos con foco en mar están generalmente determinados con incertidumbres muy grandes, debido a la ausencia de estaciones sísmicas en tierra cercanas al epicentro. El uso de redes con estaciones en tierra en combinación con sismómetros de fondo marino, ayuda a reducir la distancia epicentral y la brecha acimutal, mejorando de este modo la estimación de los parámetros focales de los terremotos. Un segundo factor muy importante que mejora la precisión de las localizaciones hipocentrales, es el uso de un modelo de Tierra específico para la zona de estudio y en especial modelos 3D. Aunque, actualmente existen numerosos modelos de Tierra 1D globales, es difícil encontrar modelos 3D para determinadas regiones y que además tengan una resolución valida, para poder ser utilizados como modelos en las localizaciones hipocentrales de terremotos. Además de estos factores, es importante utilizar un método de búsqueda hipocentral que sea robusto, como es el caso de los métodos actuales de localización no lineal, y que se pueda aplicar con modelos de Tierra 3D.

El objetivo de esta tesis es, estudiar como la utilización de redes sísmicas con estaciones en tierra y una pequeña red de sismómetros de fondo marino junto con modelos de Tierra 3D y métodos de búsqueda hipocentral no lineales, influyen en la estimación de los parámetros focales de terremotos ocurridos en la región del Cabo San Vicente. La región de estudio de esta tesis es el Sur-Oeste de Iberia, cerca del Cabo San Vicente. En esta región han sucedido algunos de los terremotos y tsunamis más grandes en Europa, como el terremoto de Lisboa de 1755 (Imax = X). La red submarina utilizada en esta tesis está formada por cinco sismómetros de fondo marino, de banda ancha, situados a 200 km de la costa del Cabo San Vicente y que servirán para estudiar la sismicidad de la región durante el período que duro el despliegue (7 meses), junto con las estaciones de las redes sísmicas permanentes del Instituto Geográfico Nacional, del Instituto Português do Mar e Atmosfera y la red Western Mediterranean.

Se han estudiado 38 terremotos, con magnitudes 2 < mbLg < 5, que ocurrieron durante el periodo de despliegue de la red submarina (septiembre 2015 – abril 2016). Estas localizaciones se han llevado a cabo usando tres modelos de tierra diferentes. Un modelo de tierra 1D y dos modelos 3D.

En la tesis se distinguen tres partes bien diferenciadas: La primera parte, está basada en la localización de terremotos con modelos de tierra 3D y usando un método de búsqueda y localización no lineal, en la segunda parte, se realiza la determinación y el análisis del mecanismo focal de siete terremotos, también con foco en mar y en la tercera parte, se expone y analiza la nueva metodología diseñada en esta tesis que, sin la necesidad de utilizar estaciones en tierra, utiliza la red submarina como una antena sísmica para determinar los epicentros de los terremotos.

Los resultados de esta tesis demuestran que el uso de la red sísmica submarina, aunque sea de pequeño tamaño, y de modelos de Tierra 3D de alta resolución, mejoran la determinación de los parámetros focales de todos los terremotos analizados en comparación con las localizaciones del catálogo del Instituto Geográfico Nacional. Concretamente, la reducción de incertidumbre asociada a la profundidad y al epicentro mejora en hasta un 50%.

Respecto a los mecanismos focales, se muestra cómo los resultados obtenidos utilizando la red combinada de estaciones en tierra y mar, concuerdan con el patrón de estrés en la región de esfuerzos regionales en la zona objeto de estudio de esta tesis. Además, los resultados de los mecanismos focales de los terremotos más grandes, con localizaciones próximas a los terremotos analizados, también son coherentes.

En lo referente a la nueva metodología desarrollada en esta tesis, se mostrara como esta es un robusto algoritmo capaz de facilitar la detección de señales sísmicas, que están altamente contaminadas en ruido oceánico. Además, la metodología es eficaz incluso siendo utilizada en redes submarinas que cuenten con pocos sensores, como es el caso de esta tesis, y además es flexible en cuanto a las dimensiones geométricas de la red utilizada. Los resultados de las localizaciones epicentrales obtenidas utilizando esta metodología y con la red submarina están en acuerdo con las determinadas mediante métodos de localización convencional, mostrando que no existen diferencias significativas (distancias medias entre epicentros < 12 km), entre ambas determinaciones.

Finalmente, se obtiene que la utilización de una red sísmica submarina junto con modelos de Tierra 3D y métodos de búsqueda hipocentral no lineales, mejora significativamente la localización de los terremotos en la región de estudio, y además que la red submarina se puede utilizar como una antena sísmica capaz de localizar de forma precisa sismos regionales.



Abstract

The focal parameters of earthquakes with an offshore focus are generally determined with very large uncertainties, due to the absence of seismic stations on land, close to the epicenter. The use of networks with ground stations in combination with ocean bottom seismometers helps to reduce the epicentral distance and the azimuthal gap, thus improving the estimation of the focal parameters of earthquakes. A second important factor that improves the precision of hypocentral locations is the use of a specific Earth model for the study area and especially 3D models. Although there are currently numerous global 1D Earth models, it is difficult to find 3D models for certain regions and that also have a valid resolution, to be able to be used as models in hypocentral earthquake locations. In addition to these factors, it is important to use a hypocentral search method that is robust, as is the case with current non-linear location methods, and can be applied with 3D models.

The objective of this thesis is to study how the use of seismic networks with stations on land and a small network of ocean bottom seismometers together with 3D Earth models, influence the estimation of the focal parameters of earthquakes that occurred in the study region (the South- West of Iberia, near Cape Saint Vincent). Some of the largest earthquakes and tsunamis in Europe have occurred in this region, such as the Lisbon earthquake of 1755 (Imax = X). The submarine network used in this thesis is made up of five broadband ocean bottom seismometers located 200 km from the coast of Cape Saint Vincent and which will serve to study the seismicity of the region during the period that the deployment lasted (8 months), together with the stations of the permanent seismic networks of the Instituto Geográfico Nacional, the Instituto Português do Mar e Atmosfera and the Western Mediterranean network. It has been studied 38 earthquakes, with magnitudes 2 <mbLg <5, which occurred during the submarine network deployment period (September 2015 – April 2016). These locations have been carried out using three different land models. One 1D earth model and two 3D models.

In the thesis three well differentiated parts are distinguished: The first part, is based on the location of earthquakes with 3D earth models and using a non-linear search and location method, in the second part the determination and the analysis of the focal mechanism of seven earthquakes, also with a focus on the sea and in the third part, the new methodology designed in this thesis is exposed and analyzed which, without the need to use land stations, uses the submarine network as a seismic array to determine the epicenters of earthquakes.

The results of this thesis show that the use of the submarine seismic network, even if it is small, and of high-resolution 3D Earth models, improve the determination of the focal parameters of all the earthquakes analyzed compared to the locations in the catalog of the Instituto Geográfico Nacional. Specifically, the reduction of uncertainty associated with depth and epicenter improves by up to 50%. Regarding the focal mechanisms, it will be shown how the results obtained using the combined network of land and sea stations, agree with the stress pattern in the region of regional efforts in the area under study of this thesis. Furthermore, the results of the focal mechanisms of the largest earthquakes, with locations close to the earthquakes analyzed, are also consistent.

It is also shown how the methodology is a robust algorithm capable of detecting seismic signals, which are highly polluted in oceanic noise. Furthermore, the methodology is effective even when used in underwater networks that have few sensors, as is the case in this thesis, and it is flexible in terms of the geometric dimensions of the network used. The results of the epicentral locations obtained using this methodology and with the submarine network are in agreement with those determined by conventional location methods, showing non-significant differences (mean distance between epicenters <12 km), between both determinations.

Finally, the use of a submarine seismic network together with 3D Earth models and non- linear hypocentral location methods, provides a significant improvement in the location of earthquakes, and also that the submarine network can be used as a seismic array of precise epicentral location.

 



Índice

Resumen / Summary, 1
Capítulo 1. Introducción / Introduction, 5
Capítulo 2. Metodología, 11
2.1 Localización hipocentral, 11
2.1.1 Localización hipocentral por métodos lineales, 11
2.2.2 Localización hipocentral por métodos no lineales, 18
2.2 Mecanismo focal, 27
2.2.1 Determinación de la orientación del plano de fractura, 27
2.3. Técnica de Antena Sísmica para tratamiento de señales incoherentes basado en la Transformada Wavelet Continua, 29
2.3.1 Técnicas de array, 29
2.3.2 Respuesta de la Antena Sísmica, 31
2.3.3 Desarrollo de nueva metodología para la estimación del vector lentitud por medio de la Continuos Wavelet Transform (CWT), 34
2.3.4 Método de Amplitudes (AM), 38
2.3.5 Broad-Band Frequency-Wavenumber Beamforming (BB-FK), 40
2.3.6. Estimación de las coordenadas epicentrales por medio de las CFs, 42
Capítulo 3. Procesado y análisis de las observaciones, 43
3.1. Observaciones, 43
3.2 Ocean Bottom Seismometer (OBS), 45
3.3 Determinación de la respuesta instrumental, 47
3.4 Detección y selección de eventos sísmicos, 52
3.5 Corrección de la señal de tiempo de los OBSs, 53
3.6 Corrección de la orientación de las componentes horizontales de los OBSs, 59
3.7 Estimación del ruido sísmico en OBSs, 63
3.8 Modelos de Tierra, 68
3.9 Software, 69
3.9.1 CWT-Seismology, 69
3.9.2 Software para la localización de terremotos con método no lineal: NonLinLoc, 70
3.9.3 Software para la localización de terremotos con método lineal: SeisComP3, 70
3.9.4 Librerías de código adicionales, 70
Capítulo 4. Resultados, 71
4.1 Localizaciones hipocentrales, 72
4.1.1 Selección de terremotos, 72
4.1.2 Resultados de la localización hipocentral NLL, 75
4.1.3 Comparación con localizaciones IGN, 85
4.1.4 Comparación de resultados en función del modelo de Tierra, 89
4.1.5 Casos especiales: terremotos localizados fuera de la red de OBSs, 91
4.1.6 Comparación de resultados con trabajos previos, 93
4.2 Mecanismo focal: orientación del plano de falla, 95
4.3 Resultados vector lentitud, 100
4.3.1 Estimación del vector lentitud mediante la técnica CWT, 100
4.3.2 Eficacia y robustez del método CWT, 106
5. Conclusiones / Conclusions, 117
6. Referencias, 121
Anexo I. Coordenadas de los OBSs del array, 133
Anexo II. Curvas de sincronización de los relojes de los OBSs, 135
Anexo III. Terremotos utilizados para la orientación de las componentes horizontales de los OBSs, 137
Anexo IV. Probability Density Functions (PDFs) del array de OBSs, 139
Anexo V. Modelos de Tierra 3D, 141
Anexo VI. Parametrización del método CWT, 143
A6.1 Algoritmo y pseudocódigo, 143
A6.2 Parametrización del algoritmo CWT, 144
Anexo VII. Resultados de las localizaciones hipocentrales, 147
A7.1 Parametrización, 147
A7.2 Tablas de relocalizaciones hipocentrales, 148
Anexo VIII. Mecanismos Focales, 155


Conclusiones

Los resultados de esta tesis permiten obtener las siguientes conclusiones:

• Es uso de un array de OBSs para la localización de terremotos con foco en mar (aunque
sea de solo unos pocos OBSs,) junto con estaciones en tierra, mejora la calidad de los
hipocentros, disminuyendo la incertidumbre de la solución. La mejora más notable
corresponde a la profundidad.
• Los modelos de Tierra 3D específicos para la región, son altamente recomendables frente
a los modelos 1D, ya que reducen las incertidumbres de las localizaciones de terremotos
hasta en un 50%.
• Los métodos de localización no lineal son preferibles frente a los métodos lineales, y en
especial en terremotos con una mala cobertura azimutal de estaciones y ausencia de
estaciones a distancias pequeñas. El parámetro que mejor refleja esto es la profundidad.
• El uso conjunto de estaciones en tierra y OBSs, junto con un método de localización no
lineal y un modelo de Tierra 3D específico para la zona (Cabo S.Vicente – Estrecho de
Gibraltar) ha facilitado la obtención de buenas localizaciones hipocentrales con los datos
de la campaña 2015-2016 del Cabo S. Vicente.
• Los resultados de los datos de esta campaña han permitido confirmar la ocurrencia de focos
activos entre 30-50 km de profundidad en la región B, zona entre en Cabo S. Vicente y el
array de OBSs.
• Se ha propuesto y desarrollado una metodología, basada en la obtención de funciones
características por medio de la Transformada Continua Wavelet (CWT), para poder estimar
el vector lentitud y el epicentro de terremotos regionales con arrays de gran apertura y
pocas estaciones.
• La aplicación del método CWT a los datos de la campaña marina 2015-2016 ha permitido
comprobar la eficacia del método CWT para localizar epicentros de terremotos regionales.
Como trabajo futuro derivado de los resultados obtenidos en esta tesis, se proponen las
siguientes líneas de investigación:
• Mejorar el algoritmo de búsqueda y localización de terremotos basados en la inversión de
los tiempos de recorrido de las ondas sísmicas explorando nuevas posibilidades en las
nuevas y emergentes técnicas de inversión de la forma de onda y stacking que han
demostrado resultados muy prometedores en los últimos años (Li et al., 2020).
• Estudiar la posible incorporación al Sistema de Alerta Temprana del sur de la Península
Ibérica (Carranza et al., 2017) una red de OBSs usando las nuevas tecnologías de telemetría
en tiempo cuasi real de redes de fibra óptica y/o redes OBSs (Berger et al., 2016; Williams
et al., 2019; Kohler et al., 2020).
• Mejorar la capacidad de detección del algoritmo CWT de señales emergentes y con baja
SNR a través de la optimización del compromiso entre la resolución en tiempo y frecuencia
de la wavelet seleccionada, ya que actualmente únicamente se utiliza la Morlet Wavelet
compleja.
• Explorar técnicas que permitan la reducción del ruido de las componentes horizontales de
los OBSs (Negi et al., 2020).
• Procesar de forma conjunta los datos procedentes de la campaña marina del proyecto
NEAREST (Integrated observations from NEAR shore sourcES of Tsunamis: towards an
early warning system ,2007-2008 (Geissler et al., 2010)), junto con los de la campaña
2015-2016 de esta tesis para determinar el espesor de la corteza bajo la zona de despliegue
de OBSs, con inversión de funciones receptoras y/o autocorrelaciones (Langston, 1979;
Phạm y Tkalčić, 2018).
• Realización de nuevas campañas marinas con OBSs en otras zonas sísmicamente activas
del sur de Iberia (v.g., Mar de Alborán).