SM7.2 | Volcano seismology: observations and modelling
EDI
Volcano seismology: observations and modelling
Convener: Ivan Lokmer | Co-conveners: Chris Bean, Kristín Jónsdóttir, Arthur Jolly
Orals
| Thu, 18 Apr, 16:15–18:00 (CEST)
 
Room -2.47/48
Posters on site
| Attendance Fri, 19 Apr, 10:45–12:30 (CEST) | Display Fri, 19 Apr, 08:30–12:30
 
Hall X1
Orals |
Thu, 16:15
Fri, 10:45
Volcanic seismicity is fundamental for monitoring and investigating volcanic systems, their structure and their underlying processes. Volcanoes are very complex objects, where both the pronounced heterogeneity and topography can strongly modify the recorded signals for a wide variety of source types. In source inversion work, one of the challenges is to capture the effect of small-scale heterogeneities in order to remove complex path effects from seismic data. This requires high-resolution imagery, which is a significant challenge in heterogeneous volcanoes. In addition, the link between the variety of physical processes beneath volcanoes and their seismic response (or lack of) is often not well known, leading to large uncertainties in the interpretation of volcano dynamics based on the seismic observations. Taking into account all of these complexities, many standard techniques for seismic analysis may fail to produce breakthrough results.

In order to address the outlined challenges, this session aims to bring together seismologists, volcano and geothermal seismologists, wave propagation and source modellers, working on different aspects of volcano seismology including: (i) seismicity catalogues, statistics and spatio-temporal evolution of seismicity, (ii) seismic wave propagation and scattering, (iii) new developments in volcano imagery, (iii) seismic source inversions, and (iv) seismic time-lapse monitoring. Contribution on controlled geothermal systems in volcanic environments are also welcome. Contributions on developments in instrumentation and new methodologies (e.g. Machine Learning) are particularly welcome.
By considering interrelationships in these complementary seismological areas, we aim to build up a coherent picture of the latest advances and outstanding challenges in volcano seismology.

Orals: Thu, 18 Apr | Room -2.47/48

Chairpersons: Ivan Lokmer, Kristín Jónsdóttir, Arthur Jolly
16:15–16:20
16:20–16:30
|
EGU24-4796
|
ECS
|
Virtual presentation
Felix Rodrigo Rodriguez Cardozo, Jochen Braunmiller, Taha Sadeghi Chorsi, Gerardo Mendo Pérez, Vala Hjörleifsdóttir, Kristín Jónsdóttir, Yesim Cubuk Sabuncu, Glenn Thompson, Stephen McNutt, Jacqueline Dixon, Timothy Dixon, and Rocco Malservisi

In March 2021, a fissure eruption in Fagradalsfjall marked the onset of the first volcanic activity in over 800 years on the Reykjanes Peninsula, Iceland. Since then, three more fissure eruptions occurred in August 2022, July 2023, and December 2023. Intense seismic swarms that included Mw≥ 4.0  earthquakes preceded all eruptions. Concurrently with swarm activity, large surface deformations related to magma intrusions were observed by Interferometric Synthetic Aperture Radar (InSAR), and Global Navigation Satellite System (GNSS) data. Ground displacements exceeded the deformation expected from earthquakes by far, suggesting the intrusion process was primarily aseismic. However, the intrusions may have influenced the prevalent earthquake faulting style in the pre-eruptive swarms. For instance, seismic moment tensors before and during the early 2021 swarm indicate that right-lateral bookshelf faulting along roughly north-south trending strike-slip faults dominated seismic deformation. This changed to more northeast-southwest oriented oblique-normal to normal faulting mechanisms for the later swarms consistent with graben formation after a dyke intrusion. Seismic moment release during the 2021 seismic swarm, which included ten relatively large Mw 5+ earthquakes, was larger than for subsequent swarms, where only a few events reached Mw 5. This may indicate that the 2021 intrusion, which marked a reawakening of volcanic activity, may have triggered bookshelf faults close to failure that might have otherwise ruptured in the near term due to the underlying oblique rifting.  The transition of source mechanisms toward oblique and normal faulting during the later swarms, though, may reflect an active role of the intrusion processes on fault orientation of triggered seismicity rather than simply inducing seismicity on pre-existing faults. 

How to cite: Rodriguez Cardozo, F. R., Braunmiller, J., Sadeghi Chorsi, T., Mendo Pérez, G., Hjörleifsdóttir, V., Jónsdóttir, K., Cubuk Sabuncu, Y., Thompson, G., McNutt, S., Dixon, J., Dixon, T., and Malservisi, R.: The role of magma intrusions on faulting mechanisms during pre-eruptive seismic swarms in the Reykjanes Peninsula, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-4796, https://doi.org/10.5194/egusphere-egu24-4796, 2024.

16:30–16:40
|
EGU24-10597
|
ECS
|
On-site presentation
Alea Joachim, Eva P. S. Eibl, Daniel Müller, and Thomas R. Walter

On March 19, 2021, an effusive eruption lasting for six months began in the Geldingadalir valley on the Reykjanes peninsula, in the southwest of Iceland. This eruption was characterised by episodic lava effusion from 2 May to 18 September and changing eruptive behaviour. Here, we analyse five of such effusion episodes of 8 June 2021 by using drone video data acquired over the active crater lake together with volcanic tremors that were recorded by a seismometer located at 5.5 km distance from the active vent. We are thus able to study each of the five episodes in terms of tremor amplitude evolution and its frequency spectrum, and compare it with the timing, height and dynamics of the lava lake, its bubbling and crater overflow. We observe a slow rise of the lava lake by 25.1 to 26.2 metres within 6.15 ± 2.35 minutes, followed by a rapid fall of the lava lake surface to its previous level within 1.6 ± 0.12 minutes. In contrast, the quiescence period in the tremor lasts ~ 10.12 ± 0.68 minutes followed by another 2.55 ± 0.2 minutes of tremor. Thus, the duration of tremor generation is shorter than the time required for the lava lake to reach its maximum height. Furthermore, the tremor amplitude reaches its maximum after the lava lake has started to sink. We discuss the volcanic tremor generation in relation to lava lake elevation and related processes.

How to cite: Joachim, A., Eibl, E. P. S., Müller, D., and Walter, T. R.: Linking volcanic tremor amplitude to drone records of lava lake level changes during the 2021 Geldingadalir eruption, Iceland, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-10597, https://doi.org/10.5194/egusphere-egu24-10597, 2024.

16:40–17:00
|
EGU24-21112
|
solicited
|
On-site presentation
Anne Obermann, Bettina Goertz-Allmann, Pilar Sanchez-Pastor, Peidong shi, Sin-Mei Wu, and Vala Hjórleifsdóttir

Hengill volcano and its associated geothermal fields represent Iceland's most productive harnessed high-temperature geothermal fields, where energy is provided by cooling magmatic intrusions connected to three volcanic systems. The crustal structure in this area is highly heterogeneous and shaped by the intricate interplay between tectonic forces and magmatic/hydrothermal activities, making detailed subsurface characterization challenging. In the Northern part of the Hengill geothermal field, super-hot geothermal resources have been spotted that are currently considered for geothermal exploration.

Over the past years, we have studied the site in great detail, and acquired high-quality datasets from a 40+broadband seismic station array, a dense 500+ station nodal array and distributed acoustic sensing data from a fibre line crossing the area. We compare the results that we obtained from various seismic imaging methods e.g., earthquake tomography, ambient noise methods and discuss their potential and limitation to enable high-resolution seismic methods for exploration and monitoring of geothermal plays in such complex volcanic environments.

How to cite: Obermann, A., Goertz-Allmann, B., Sanchez-Pastor, P., shi, P., Wu, S.-M., and Hjórleifsdóttir, V.: High-Resolution Seismic Methods for Geothermal Exploration and Monitoring across the Hengill volcanic area, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-21112, https://doi.org/10.5194/egusphere-egu24-21112, 2024.

17:00–17:10
|
EGU24-7488
|
ECS
|
On-site presentation
Sophie Butcher, Andrew Bell, Stephen Hernandez, Peter La Femina, James Grannell, and Mario Ruiz

The 2018 eruption at Sierra Negra volcano, Galapagos Islands, was accompanied by 8.5 metres of caldera subsidence and intense seismicity, resulting from deflation of a shallow sill-like magma reservoir at ~2 km depth. High-precision hypocentre locations from manually picked phase arrivals show that earthquake sources are tightly constrained within a complex, multi-stranded trapdoor fault system (TDF) above the sill. However, the incompleteness of this high-precision catalogue leaves outstanding questions about the spatio-temporal evolution of seismicity through the eruption, and how it relates to deformation and magma efflux.

Here we present the results of an automated workflow to streamline the production of a ‘temporally-enriched’ seismic catalogue for 2018 eruption at Sierra Negra. We initially utilise PhaseNet, a deep-neural-network-based automatic phase picker, to identify events missing from the initial manually picked catalogue, expanding the detections from 1,618 to 9,871 events. Our catalogue identifies more events in the immediate aftermath of the Mw5.4 earthquake that initiated the eruption, and new small magnitude events (< ML2.0) in the period more than 72 hours after the eruption onset. We then use a template matching approach to further supplement these detections. Specifically, these events fill gaps in the catalogue where tremor amplitudes make manual event detection more difficult. Hypocentre locations for newly detected events are also constrained to the TDF zone, however there is more variety in depth estimates. This has implications for how we consider the TDF with depth, and allows us to consider other potential sources of seismicity in the system.

Our workflow offers an efficient method of producing ‘temporally-enriched’ catalogues at Sierra Negra, and can be readily adapted for the sparse seismic network that remains during the current inter-eruptive phase. However, our experience at Sierra Negra suggests that applying automated earthquake detection and location methods can be challenging in volcanic settings, and requires careful parameterization and quality control.

How to cite: Butcher, S., Bell, A., Hernandez, S., La Femina, P., Grannell, J., and Ruiz, M.: Temporal enrichment of the seismic record of the 2018 eruption at Sierra Negra using deep neural networks, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-7488, https://doi.org/10.5194/egusphere-egu24-7488, 2024.

17:10–17:20
|
EGU24-13020
|
On-site presentation
Karina Bernal-Manzanilla, Marco Calò, Karina Eloisa Rodríguez García, Daniel Martínez Jaramillo, and Sébastien Valade

Popocatépetl, one of Mexico's most active volcanoes, poses significant risks to the dense populations in its vicinity. Effective monitoring of its seismic activity is crucial for understanding and mitigating these hazards. This study employs data collected with a network of 19 seismic stations surrounding the volcano, combined with machine learning techniques and spatial coherence methods, to generate comprehensive seismic catalogs spanning from 2019 to the present. Our automated workflow includes the identification and localization of long period (LP) events, tremors, and volcano-tectonic (VT) earthquakes.

For this purpose, an improved classification model based on Support Vector Machines was developed to distinguish LP events and tremors within continuous recordings. Their locations were determined using a cross-correlation-based method. Additionally, the VT earthquake catalog was compiled using deep learning-based models for phase picking, followed by standard location methods. Our findings not only corroborate trends observed in manual analyses at the volcano's observatory but also uncover additional events, highlighting trends in the volcano’s dynamics not observed before.

To showcase the use of these catalogs, we will present a multiparametric analysis integrating seismic data with thermal anomalies, SO2 emissions, and GPS measurements. This research not only deepens our comprehension of volcanic processes but also underscores the transformative role of technology in geophysical research.

 

Research supported by the program UNAM-DGAPA-PAPIIT: IN103823.

How to cite: Bernal-Manzanilla, K., Calò, M., Rodríguez García, K. E., Martínez Jaramillo, D., and Valade, S.: Dynamics of the Popocatépetl Volcano, Mexico, revealed by Machine Learning-Based Seismic Catalogs, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-13020, https://doi.org/10.5194/egusphere-egu24-13020, 2024.

17:20–17:30
|
EGU24-15278
|
Virtual presentation
Clément Hibert, Joachim Rimpot, Lise Retailleau, Jean-Marie Saurel, Jean-Philippe Malet, Germain Forestier, Jonathan Weber, Tord S. Stangeland, Antoine Turquet, and Pascal Pelleau

Continuous seismological observations provide valuable information to deepen our understanding of processes occurring in both aerial and submarine volcanoes. However, the wealth of the seismicity recorded near volcanoes makes exhaustive exploration of these seismological chronicles very complex and time-consuming. In this study, we present a systematic analysis of two months of seismological records using a self-supervised learning (SSL) approach for the unsupervised clustering of continuous seismic data acquired by ocean bottom seismometers deployed in the vicinity of the Fani Maoré volcano (Mayotte). The proposed clustering process allows the identification of individual seismic events, seismic crisis and tremors that would be challenging to observe using conventional approaches. We show that our approach detects and classifies both known and new events, including two eruptive sequences previously unknown. We also demonstrate the potential of self-supervised methods for the analysis of seismological records, providing a synoptic view and facilitating the discovery of insightful yet rare events. This approach has numerous applications in exploring various seismological datasets, simplifying analysis while making it more comprehensive.

How to cite: Hibert, C., Rimpot, J., Retailleau, L., Saurel, J.-M., Malet, J.-P., Forestier, G., Weber, J., Stangeland, T. S., Turquet, A., and Pelleau, P.: Self-supervised learning for the exploration of continuous seismic records at the Fani Maoré submarine volcano (Mayotte), EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-15278, https://doi.org/10.5194/egusphere-egu24-15278, 2024.

17:30–17:40
|
EGU24-4438
|
On-site presentation
Zacharie Duputel, Lucile Costes, Valérie Ferrazzini, and Olivier Lengliné

Volcanoes often exhibit very long period (VLP) signals, with periods ranging from 2 to 100 seconds. Due to their very long wavelengths, these waveforms are not strongly affected by volcano structural heterogeneities and provide invaluable insights into dike and magma properties that are not easily accessible though other observations. Historically, only a few VLP events have been recorded at the Piton de la Fournaise, primarily associated with collapses in the summit caldera in 1986, 2002 and 2007. However, since 2010, the monitoring network has evolved significantly and it is now equipped with several broadband stations, allowing a wide range of signals to be recorded. In this study, we show that VLP events are quite common at Piton de la Fournaise. Specifically, we identify swarms of VLP events during eruptions and during magma injections preceding eruptions. Source analysis of VLP events during eruptions indicates the resonance of the dike during sudden decreases in magma flow. Pre-eruptive VLP waveforms exhibit significant differences from those observed during eruptions, pointing to a distinct source process as the dike propagates laterally toward the volcano flank.

How to cite: Duputel, Z., Costes, L., Ferrazzini, V., and Lengliné, O.: Very long-period observations at Piton de la Fournaise volcano, La Réunion, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-4438, https://doi.org/10.5194/egusphere-egu24-4438, 2024.

17:40–17:50
|
EGU24-16604
|
On-site presentation
Andres Barajas, Nikolai Shapiro, and German Prieto

We present observations showing that during episodes of volcanic tremors, the phase of inter-station cross-correlations becomes stable. We propose a new quantity, the phase coherence, to identify the differential phase stability in recordings obtained from a single pair of stations, which is extrapolated to the seismic network. Then, we present a new approach based on the estimation of differential travel times through the differential phase measurements, to locate the sources of tremors occurring at the end of 2015 at the Klyuchevskoy Volcanic Group in Kamchatka, Russia. We present evidence supporting the existence of two types of activity happening simultaneously during the tremor episode: the main tremor source, originating from a region located between 7 and 9 km depth under the main volcanoes, and the widespread occurrence of weak low-frequency earthquakes occurring at random locations. We show how the phase coherence and the differential phases can be used to provide information on the stability of the tremor source position and to estimate its location.

How to cite: Barajas, A., Shapiro, N., and Prieto, G.: Differential phase analysis for volcanic tremor detection and source location., EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-16604, https://doi.org/10.5194/egusphere-egu24-16604, 2024.

17:50–18:00
|
EGU24-20804
|
On-site presentation
Georg Rümpker and Fabian Limberger

In this study we utilize 3D numerical models to simulate seismic resonances in a volcanic edifice, arising from the interaction between an externally excited wavefield and a magma chamber-conduit system. The resultant wavefield holds the potential to provide significant insights into the properties of the magmatic system. Contrary to previous assumptions that required an internal source, our findings show that the magma chamber and conduit efficiently capture the incident wavefield of both P- and S-waves, excited by a high-frequency (~10 Hz) earthquake located within the edifice. Due to multiple internal reflections off the boundaries of the chamber and the conduit, prolonged reverberations occur, which are guided along the conduit. Temporal and spectral analyses of synthetic seismograms illustrate that the size of the magma chamber and the width of the conduit are critical in determining the magnitude and dominant frequencies of the seismic resonances. Specifically, models with larger magma chambers and wider conduits consistently yield larger resonance amplitudes at distinct frequencies. At greater distances from the conduit, an intensified scattered wavefield with a broad frequency range indicates the presence of a substantial magma chamber within the volcanic edifice. Resonance frequencies reach up to 23 Hz, underscoring significant frequency shifts. In general, these externally initiated resonances may appear as tremor-like signals at seismic stations on the edifice, accompanying more conventional seismic events in its proximity.

How to cite: Rümpker, G. and Limberger, F.: Numerical modeling predicts seismic resonances in the magma chamber-conduit system due to wavefield capturing, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-20804, https://doi.org/10.5194/egusphere-egu24-20804, 2024.

Posters on site: Fri, 19 Apr, 10:45–12:30 | Hall X1

Display time: Fri, 19 Apr 08:30–Fri, 19 Apr 12:30
X1.47
|
EGU24-6379
|
ECS
Emmanuel Caballero-Leyva, Nikolai Shapiro, Cyril Journeau, Léonard Seydoux, Jean Soubestre, and Andrés Barajas

Continuous seismic signals recorded in the vicinity of active volcanoes are composed of seismic waves generated by a variety of internal and environmental sources and propagating through different parts of the plumbing system. This implies that these signals are very sensitive to the state of the plumbing system. A change in the volcanic activity affects the properties of the seismo-volcanic sources while a change in the plumbing structure affects the media through which the seismic waves propagate. Network-based analysis of continuous seismic records has been developed to incorporate information from multiple stations simultaneously. Here we use an approach based on the network covariance matrix that combines an ensemble of inter-station cross-correlations. We compute the width of the eigenvalue distribution of this matrix at a given frequency in a moving time window, resulting in a compact time-frequency representation of continuously recorded seismic wavefield.

 We apply this analysis to ten years (2013-2023) of continuous seismic data from the Piton de la Fournaise volcano located in la Réunion, France. The resulting spectral width distributions indicate that continuous signals are characterized by multiple narrow spectral peaks, which are observed during co-eruptive tremors as well as during periods without visible volcanic activity. We propose a normalization process to enhance these peaks in both the frequency and time domains. We observe numerous spectral peaks in the 1-3 Hz frequency band that remain nearly constant for extended periods (weeks to months). We observe a distinct difference in the spectral peak distribution between co-eruptive and quiet periods, as well as significant variations during long-standing eruptions. Locations of sources of the co-eruptive signals correlate well with the eruption sites. The inter-eruptive signals seem to originate from a combination of environmental and weak internal sources, and changes in their spectral properties might reflect the medium changes after major eruptions.

How to cite: Caballero-Leyva, E., Shapiro, N., Journeau, C., Seydoux, L., Soubestre, J., and Barajas, A.: Long-term evolution of continuous seismic signals at Piton de la Fournaise volcano inferred from the network covariance matrix., EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-6379, https://doi.org/10.5194/egusphere-egu24-6379, 2024.

X1.48
|
EGU24-13932
|
ECS
Yuly Paola Rave-Bonilla, Mel Rodgers, Félix Rodríguez-Cardozo, Jochen Bruanmiller, and John Makario Londoño

Volcano seismology is a key contributor for assessing the activity status of volcanoes and forecasting future eruptive behaviour. One measurement that may indicate a change in the volcanic behaviour is a temporal change in the spectral content, or frequency distribution, of seismic waveforms. Significant changes in frequency distribution have been observed at Soufrière Hills, Montserrat; Telica Volcano, Nicaragua; Redoubt Volcano, Alaska USA, and at other volcanoes prior to volcanic eruptions. Identification of such changes in spectral energy could indicate changes in volcanic activity, and hence be used in forecasting, as well as to investigate the physical processes behind eruptive processes. With the objective of identifying possible changes in the spectral energy distribution of the seismic sources of the Nevado del Ruiz Volcano (NRV) prior to eruptions, we analysed the spectral energy of volcano tectonic (VT), low frequency (LF) and hybrid seismic events during 2012, when the NRV had two Volcanic Explosivity Index (VEI) 2 eruptions on May 29th and June 30th. We analysed nearly 27,000 events and implemented a semi-automatic pre-processing pipeline in ObsPy for selecting the optimal stations and seismograms based on the signal-to-noise ratio and proximity to seismic clusters. Then, we cut the seismograms to isolate the seismic signals and bandpass-filtered the data before calculating metric such as dominant frequency, ratio of high to low spectral energy. In this work, we present preliminary results on the temporal changes in spectral energy of the seismic events at NRV and whether this could be linked, along with other geophysical measurements, with changes to eruptive behaviour.

How to cite: Rave-Bonilla, Y. P., Rodgers, M., Rodríguez-Cardozo, F., Bruanmiller, J., and Londoño, J. M.: Spectral energy distribution of Nevado del Ruiz Volcano seismicity near eruptive events, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-13932, https://doi.org/10.5194/egusphere-egu24-13932, 2024.

X1.49
|
EGU24-11355
|
ECS
Maurice Weber, Christopher Bean, Patrick Smith, Ivan Lokmer, Luciano Zuccarello, Silvio De Angelis, Jean Soubestre, and Vittorio Minio

When it comes to volcanic tremor, low frequency signals (below 5 Hz) are well investigated. Such tremor signals can usually be linked to magma movement or gas fluctuations. However, little is known about seismic tremor signals on Mount Etna above 10 Hz. Hence, a large field campaign targeting high frequencies was undertaken in the summer of 2022. It consisted of the deployment of six dense circular arrays ranging from 30 to 200 m apertures of seismic nodes installed around the summit craters. It led to the detection of tremor bands between 10 and around 20 Hz as well as the typical tremor signals below 5 Hz.

The tremor is detected with good coherency at stations within one array (despite an extreme level of scattering) in good agreement with the energy distribution in the average amplitude spectra of the array. The high frequency tremor varies strongly in intensity over time periods of one hour and re-occurs several times throughout the deployment period of almost a week. In contrast the tremor below 5 Hz is relatively constant. This suggests that the high frequency tremor could be a separate signal due to a process that may not yet be fully understood.

Localisations of these tremor episodes point to or near the Bocca Nuova Summit Crater which was actively degassing at the time. Interestingly, high frequency seismic tremor is matched in time very well by a narrow 3.5-5 Hz acoustic band. While the match in time clearly suggests a connection between the two signals, the different frequencies indicate two different but linked processes happening simultaneously. The acoustic signal implies degassing processes. Later during the deployment tremor episodes are found which are accompanied by much weaker acoustic signals (if at all present) suggesting gases might not necessarily be involved in generating the detected seismic tremor at all.

In summer 2023 we undertook a complementary second deployment of seismic, acoustic and optical camera data in the Bocca Nuova summit area. Once again, we find tremor below 5 Hz, however high frequency characteristics are different to the previous year with tremor bands less dominant than before and much more constant over time. More than one acoustic band is found as well, also constant over time. In this second data set we use camera recordings of the crater activity as a proxy for degassing activity to try and understand the precise origin of these seismic and acoustic volcanic signals.

How to cite: Weber, M., Bean, C., Smith, P., Lokmer, I., Zuccarello, L., De Angelis, S., Soubestre, J., and Minio, V.: Characterisation and locations of volcanic high frequency tremor above 10 Hz on Mount Etna, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-11355, https://doi.org/10.5194/egusphere-egu24-11355, 2024.

X1.50
|
EGU24-11934
|
ECS
|
Elisabeth Glück, Titouan Muzellec, Stephane Garambois, Jean Vandemeulebrouck, Þorbjörg Ágústsdóttir, Egill Árni Guðnason, and Anette Mortensen

Krafla, one of the five central volcanoes of the Northern Volcanic Zone in NE-Iceland, last erupted during the Krafla Fires in the 70s and 80s. During the same period, a geothermal power plant was built within Krafla caldera, first operated in 1978. Both scientific and industrial interest led to an increase of knowledge of the complex system through systematic exploration with a wide variety of geophysical methods including seismic and electromagnetics coupled with borehole information.
Among them, a local seismic network operated by Landsvirkjun and Iceland GeoSurvey, comprising 12 permanent broadband stations, has been continuously recording seismic data since 2013. We supplemented this network in June 2022 with a dense network of 98 nodes, resulting in two arrays, one large-scale, the other small-scale, operating in parallel.
Here we present multi-scale 3-D velocity models for P-, S- and surface waves, independently derived for both networks through local earthquake and ambient noise tomographies. These models offer a glimpse into the subsurface structures of the volcanic system by utilizing various types of waves that are responsive to distinct rock/fluid properties and depths. The relocated and clustered seismic activity, documented by both permanent and temporary networks, underscores active structures pinpointed through tomography. With this we hope to strengthen the understanding of the connected volcanic and geothermal systems. Indeed, both the seismicity and strong velocity anomalies are located at similar depths as the magma batch that was drilled into with the IDDP1.

How to cite: Glück, E., Muzellec, T., Garambois, S., Vandemeulebrouck, J., Ágústsdóttir, Þ., Guðnason, E. Á., and Mortensen, A.: Multi-scale seismic imaging and related seismicity patterns of Krafla volcano and its geothermal system, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-11934, https://doi.org/10.5194/egusphere-egu24-11934, 2024.

X1.51
|
EGU24-14227
|
ECS
Ivan Granados-Chavarria, Francesca Di Luccio, Marco Calò, Matteo Lupi, Mimmo Palano, Daniela Famiani, Federica Magnoni, Antonio Scaltrito, Laura Scognamiglio, Anna Tramelli, Andrea Ursino, Tullio Ricci, and Alessandro Marchetti

One of the tasks of the multidisciplinary project CAVEAT “Central-southern Aeolian islands: Volcanism and tEAring in the Tyrrhenian subduction system” was to deploy a dense seismic network composed by 120 wirelesss nodes, which were deployed for approximately 2 months on three islands of the Aeolian volcanic archipielago.

The main goals of this data acquisition are to:

  • detect and locate the shallow weak seismicity related to both volcanic and tectonic activity,
  • perform regional and local tomographic studies based on both passive methods and ambient noise cross-correlations (previously done for Lipari, Calò et al., 2023) to constrain the crust structure
  • characterize the source mechanisms.

In this work we show the details of the data acquisition strategy and the preliminary analyses of the continuous ambient noise records, to reconstruct the shallow structure of the Aeolian Arc. This study is part of the INGV Pianeta Dinamico project 2023-2025 CAVEAT “Central-southern Aeolian islands: Volcanism and tEAring in the Tyrrhenian subduction system” (grant no. CUP D53J19000170001) supported by the Italian Ministry of University and Research “Fondo finalizzato al rilancio degli investimenti delle amministrazioni centrali dello Stato e allo sviluppo del Paese, legge 145/2018.

How to cite: Granados-Chavarria, I., Di Luccio, F., Calò, M., Lupi, M., Palano, M., Famiani, D., Magnoni, F., Scaltrito, A., Scognamiglio, L., Tramelli, A., Ursino, A., Ricci, T., and Marchetti, A.: The ongoing seismological research for the Central-Southern Aeolian Islands (CAVEAT project), Italy, from dense terrestrial seismic arrays, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-14227, https://doi.org/10.5194/egusphere-egu24-14227, 2024.

X1.52
|
EGU24-22030
|
ECS
|
Andrea Di Benedetto, Anna Figlioli, Antonino D’Alessandro, and Giosue’ Lo Bosco

The collection of a significant catalog of seismo-volcanic data involves the selection of relevant parts of raw signals, that can be automatized by using the Short-term over Long-term Average (STA/LTA) method. Since it is parametric, the common approach to the choice is the adoption of literature-suggested parameters. To overcome these limitations, we propose a methodology for the automatic selection of STA/LTA parameters able to optimize the extraction of local events from a seismo-volcanic raw signal. The parameters are found by a grid search over an index named Quality-Numerosity Index (QNI) that measures the accordance in the automatic cuts and the consequent quantity of triggered seismo-volcanic events with the ones suggested by a human expert. The method was applied in the volcano domain, for the specific application of Explosion Quake signals extraction in Stromboli Volcano. Experiments have been conducted selecting a subset of the dataset as training where to search for the best parameters, which were subsequently adopted in a test set. The results demonstrate that the selected parameters significantly improve the quality of the extraction when compared to those extracted by adopting the parameters indicated in the literature.

How to cite: Di Benedetto, A., Figlioli, A., D’Alessandro, A., and Lo Bosco, G.: Grid-search method for STA/LTA parameters tuning: an application to Stromboli Explosion Quakes , EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-22030, https://doi.org/10.5194/egusphere-egu24-22030, 2024.

X1.53
|
EGU24-4959
Marine Menager, Zacharie Duputel, Lise Retailleau, Valérie Ferrazzini, and Ian McBrearty

Eruptions of Piton de la Fournaise volcano (Reunion Island, France) are preceded by intense pre-eruptive seismicity swarms characterized by hundreds, or even thousands of micro earthquakes (magnitude < 2). These volcano-tectonic events are triggered by the upward migration of magma toward the surface and their location provides important information regarding the future eruption location. Yet, regarding the large number of earthquakes, it is difficult to locate them all during seismicity swarms. Hence, we have implemented an approach at the Piton de la Fournaise Volcano Observatory (OVPF-IPGP) based on machine learning to automatically detect and locate these events.. First, we use PhaseNet to pick P and S waves from 17 seismic stations installed on and around the volcano. Then, phase association and source location are done using a Graph Neural Network (GNN) approach called GENIE (Graph Earthquake Neural Interpretation Engine). To implement GENIE specifically at Piton de la Fournaise, we trained the code with seismic stations and velocity models used by OVPF-IPGP to monitor the volcano.. After phase association, we perform a final hypocenter localization using the probabilistic approach of NonLinLoc. To study the results quality, we compare origin time and source location to the OVPF manual catalog as well as a catalog resulting from template matching and double-difference relocation. In particular, we focus on pre- and syn-eruptive time-periods for multiple eruptions since 2014 in order to investigate the effect of elevated seismicity rate and eruptive tremor on the performance of the workflow. We also assess to what extent the quality of the resulting automatic locations are sufficient to provide an indication of the future eruption site without expert manual input.

Applying this approach allows us to improve the monitoring of the seismicity at Piton de la Fournaise volcano. A work in progress is to implement the same approach over the entire island of La Réunion, which will enable the monitoring of other active areas in the region.

How to cite: Menager, M., Duputel, Z., Retailleau, L., Ferrazzini, V., and McBrearty, I.: Implementation of machine learning approaches to monitor pre-eruptive swarms at Piton de la Fournaise volcano, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-4959, https://doi.org/10.5194/egusphere-egu24-4959, 2024.

X1.54
|
EGU24-12228
|
ECS
Eoghan Totten, Chris Bean, and Gareth O'Brien

Recent advances in machine learning offer a new way for Earth Scientists to make predictions about geological subsurface properties. In particular, methods based on Fourier Neural Operators (FNOs) are increasingly being used as a substitute for conventional approaches based on numerical forward modelling and inversion, at a fraction of the computational cost. Most importantly, FNOs have been shown to predict accurate 2D and 3D forward modelling simulations of seismic waves up to several hundred times faster than physics-based solvers after training.

In synthetic volcanic settings to date, FNOs have been applied successfully to both the forward and inverse problem, capturing the fine-scale velocity structure of heterogeneous models. However, transferring the successful performance of simulation-trained FNOs to make accurate predictions from field-gathered seismic data is yet to be achieved. In order to accomplish this for volcanological data, training models would need to contain representative small-scale velocity heterogeneities and topography in order to produce highly scattered codas in the synthetic seismograms.

This research presents work in progress on simulation-to-real applications of FNOs using field-gathered seismic data from offshore sedimentary basin settings as a testbed environment. Historical seismic survey datasets from Atlantic sedimentary basins are often supplemented with alternative geophysical surveys and site-specific geological constraints. Combining seismic borehole and stratigraphic logs with regional seismic datasets provides a link between field-gathered seismic waveforms, stratigraphy and depth-dependent, small-scale fluctuations in seismic velocity. This in turn enables the creation of synthetic velocity models and seismograms with field-derived properties, centring the collation of data for real-world machine learning applications in the numerical domain. We aim to bring insights gained from training FNOs on a better understood seismic environment to volcanic contexts in future work.

How to cite: Totten, E., Bean, C., and O'Brien, G.: Predicting near-surface seismic data and velocity models using synthetically-trained deep learning methods: applications in data-rich environments, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-12228, https://doi.org/10.5194/egusphere-egu24-12228, 2024.

X1.55
|
EGU24-17156
|
ECS
|
Louisa Murray-Bergquist, Ayon Garcia Pina, Martin Thorwart, Christopher Ulloa, Janneke van Ginkel, Lisanne van Huisstede, Richard Wessels, and Anouk Beniest

The Ojos del Salado Volcano in Chile is the highest altitude volcano in the world, but has not been visibly active since 1993, when steam was observed rising from near the summit. Since then the volcano has been considered dormant, however, it is unclear if it could become active again, and if so, what the ramifications would be. Little is known about the size and structure of the magma chamber, or its potential interaction with hydrothermal fluids and whether this is the heat source of the nearby hydrothermal lake, Laguna Verde.

The region surrounding the volcano is very arid, but due to the cold climate at such high-altitude there are some remnant glaciers and permafrost from the last glacial maximum. Summer meltwater from these cryosphere features contributes to the water budget of the valley below, which adds to the importance of understanding the extent of remaining permafrost and the interactions between local volcanism and the cryosphere that could increase the rate of melting.

In this study we combine InSAR data with local seismicity to investigate local crustal deformation and seismic activity at the Ojos del Salado Volcano. A small seismic array, deployed in 2022, measured seismicity in the vicinity of the volcano and acted as a pilot for the 2024 campaign and network design. A dense network of geophones was deployed on the volcano’s flanks in early 2024. We present the first results of the analysis of the data from this denser network, and the pilot study, which already showed seismicity in the vicinity of the volcano. From this analysis we gain insight into the level of activity of the Ojos del Salado Volcano, local deformation patterns and the style of faulting which is also an indication of potential fluid pathways that could link the volcano to the hydrothermal lake. With this information we can better understand the interaction between this unique volcano and the local cryosphere.

How to cite: Murray-Bergquist, L., Garcia Pina, A., Thorwart, M., Ulloa, C., van Ginkel, J., van Huisstede, L., Wessels, R., and Beniest, A.: Seismic Investigation of the Ojos del Salado Volcano, Chile: The Highest Altitude Volcano in the World, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-17156, https://doi.org/10.5194/egusphere-egu24-17156, 2024.

X1.56
|
EGU24-11254
|
ECS
|
Eleanor Dunn, Chris Bean, Ivan Lokmer, and Andrew Bell

Dynamic earthquake triggering refers to the phenomenon where local seismic activity is induced by dynamic stress disturbances, originating from teleseismic earthquakes. An understanding of dynamic triggering on volcanoes offers a window into volcano stress states and seismicity initiation. Sierra Negra, a basaltic shield volcano situated on Isabela Island, Galápagos, has been the site of recurring episodes of dynamic triggering. Sierra Negra features a large elliptical summit caldera with a trap-door fault system and a magma reservoir extending 2km below the surface. Sierra Negra experienced an eruption in June 2018, characterized by a sequence of pre-eruption inflation, co-eruption deflation, and post-eruption inflation. The occurrence of dynamic earthquake triggering at Sierra Negra was observed in response to high magnitude teleseismic events from 2010 to 2018. The frequency of dynamically triggered earthquakes correlates with the increasing inflation of the magma reservoir. In this study, we aim to answer two questions: 1) How confident are we that the seismicity on Sierra Negra is dynamically triggered? And, 2) What is the location of these dynamically triggered events? Random simulations are used to calculate the likelihood that triggered events are related to teleseismic arrivals rather than being representative of local seismic activity. Results show that for the pre-2018 eruption, the likelihood that events are dynamic triggering is very high, compared to post-2018 eruption where events are more likely to be representative of local seismic activity. We only have access to a single station (VCH1) on Sierra Negra meaning the single-station location method must be used to locate all dynamically triggered events. To test and refine this method, 79 known seismic events are located using a full network from April 2018 – December 2018. Rotation of the 3-component VCH1 into the RTZ (radial-transverse-vertical) coordinate system is used to calculate the back-azimuth and the P-wave to S-wave delay is used to calculate the distance between event and station. 21 unknown dynamically triggered events are located in and around the caldera using this method. Looking forward we hope to understand the relationship between the location and timing of dynamic triggering, and its potential use in understanding volcano unrest state.

How to cite: Dunn, E., Bean, C., Lokmer, I., and Bell, A.: Understanding unrest and dynamic triggering processes on Sierra Negra, Galápagos Islands., EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-11254, https://doi.org/10.5194/egusphere-egu24-11254, 2024.