Shallow marine volcanism can produce violent phreatomagmatic eruptions that pose a significant hazard to surrounding coastal communities. Kolumbo volcano, located 7 km northeast of Santorini, is one of the most hazardous volcanic centers in the Aegean Sea. Kolumbo last erupted explosively in 1650, causing over 70 casualties and forming the 2500 m wide and 500 m deep present-day crater. While 2D seismic reflection data have provided a general overview of its the temporal evolution, the full complexity of Kolumbo’s internal structure can only be reconstructed using three-dimensional (3D) data. Here, we use high-resolution 3D seismic reflection data covering most of Kolumbo’s edifice, and present a detailed reconstruction of its internal architecture and the processes that have controlled its evolution. We show that Kolumbo’s edifice was formed by more than the previously assumed five major explosive eruption. Furthermore, we present evidence for at least one additional explosive eruption that formed a crater larger than the present-day crater from the 1650 eruption. Together, the increased explosive potential of past eruptions and the higher eruption frequency indicate a greater than was previously appreciated. Considering its proximity to Santorini, one of the most visited islands in the eastern Mediterranean, our results emphasise the need for more comprehensive natural hazard monitoring strategies at the central Aegean Sea.
How to cite: Blanch Jover, M., Karstens, J., Kopp, H., Berndt, C., Crutchley, G. J., Preine, J., and Nomikou, P.: Reconstructing Kolumbo’s explosive history: insights from 3D seismic data (Santorini), EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-10907, https://doi.org/10.5194/egusphere-egu25-10907, 2025.
The Bathymetrists Seamounts (BSM) are an elongated intra-plate volcanic province in the equatorial Eastern Atlantic Ocean. They are situated north of the Sierra Leone Rise, a smooth, aseismic seafloor elevation, believed to have formed above a hot spot at the Mid Atlantic Ridge. The arrangement of the 40 elongated seamounts, some of which are over 100 km long, is enigmatic. The strike direction of the seamounts is either parallel to the intersecting fracture and transcurrent zones or tilted by about 30° to 60°, which has been interpreted as evidence of Riedel shears. Previous age models are based on the geochemistry of a few dredged (surface) rock samples.
With the first set of ca. 4000 km of high-resolution multichannel seismic reflection data, we provide new insights into the structural and temporal evolution of the BSMs. Subsurface images reveal stacked seamount structures, offering direct evidence for their relative order of formation. Based on the DSDP site 366 findings, we establish a chrono-stratigraphic model for the sediment basins between the seamounts. This allows us to approximate the age of the underlying volcanic flanks and to refine the spatial and temporal evolution of BSM edifices through the Paleocene to Eocene. Furthermore, our data reveal signs of recent hydrothermal and magmatic activity, including intrusions, mud diapirism, and fluid chimneys penetrating the volcanic flanks. These findings highlight the BSMs as a dynamic system with both ancient origins and ongoing activity, offering new perspectives on understudied intra-plate volcanism and its associated processes.
How to cite: Hartge, M., Hübscher, C., Seidel, E., and Preine, J.: Volcanic evolution of the Bathymetrists Seamounts (equatorial Atlantic Ocean) since the Paleocene, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-2894, https://doi.org/10.5194/egusphere-egu25-2894, 2025.
Campi Flegrei is an active caldera in a densely populated area, currently experiencing significant ground uplift and seismicity. Leveraging precise relocations of extensive seismicity since 2014, we determined high resolution (250 m), 3D P- and S-wave seismic images of the inner caldera which we combine with a novel rock-physics experiment to characterize the primary features of the caldera’s 3D structure: a gas-rich reservoir below 2 km depth, a deformed caprock at 1 to 2 km depth, and a funnel-shaped, (thermo-metamorphic) basement below 3.5 km depth. Seismicity migrates downwards from the caprock to the reservoir, and, following reservoir depletion, stress loading triggers deeper, larger magnitude events along the inner-caldera boundary faults. The reservoir extent and the seismicity distribution closely correlate with the area of maximum uplift, where accelerating deformation due to pore-fluid pressure is corroborated by laboratory experiments using site-relevantin-situ rock samples. These findings suggest that coupling between gas-reservoir pressure and the fibrous microstructure of the confining caprock drives the ground uplift. This structural-dynamic reconstruction of interconnected seismic and ground deformation processes provides a framework for forecasting the evolution of unrest, which is crucial for enhancing medium- and short-term multi-hazard assessment and mitigation strategies. Our results indicate that seismic activity and the potential for a phreatic explosion should be considered as plausible scenarios, prompting a reevaluation of the hazard assessment for the area.
How to cite: De Landro, G., Vanorio, T., Titouan, M., Russo, G., Lomax, A., Virieux, J., and Zollo, A.: 3D Structure and Dynamics of Campi Flegrei Enhance Multi-Hazard Assessment, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-11536, https://doi.org/10.5194/egusphere-egu25-11536, 2025.
Integrating seismic and petrological analyses is crucial for comprehensively understanding volcanic systems. In this study, we used the Receiver Function (RF) technique to investigate the crust and upper mantle beneath the Campi Flegrei Caldera (CFc), one of the world’s most active and complex volcanic systems. Over the past decades, this caldera has exhibited signs of unrest, including significant ground deformation, degassing, and seismicity.
RF analysis is highly susceptible to seismic discontinuities, enabling the detection of significant structural features and providing critical insights into P and S-wave velocity distributions. This study extends previous research by incorporating petrological constraints to understand better the relationship between seismic velocity anomalies and the magmatic system beneath CFc.
Using data from thirteen seismic stations located within the caldera, we applied a multi-taper deconvolution method to derive RFs. We then employed Bodin et al.'s (2012) transdimensional inversion approach to retrieve 1D velocity profiles and determine the probability of seismic discontinuities. To connect geophysics and petrology, we used petrological modelling tools to estimate the liquid fraction and the rock temperature in the identified velocity anomalies, offering a robust interpretation of the magmatic system's physical state.
Our analysis revealed two significant Low-Velocity Zones beneath CFc. The first is a crustal reservoir located beneath Nisida island, extending from 8 to 16 km depth with dimensions of approximately 4 x 5 x 12 km. The Vs values are between 3.5 – 4 Km/s, and a melt fraction ranges between 0 to 5%. Below this anomaly, a deeper magma source zone was identified at depths ranging from 16 to 33 km. This layer is characterised by Vs values between 2.3 – 3.3 km/s and melt fraction ranging from 15% to 30%. These seismic models and petrological data suggest that the shallow reservoir contains only partially molten material. At the same time, the deeper zone represents a more significant magma source, potentially feeding the volcanic system. Our results also indicate the presence of an older, dense intrusive complex within the crust, which may influence the migration and storage of magma.
This study demonstrates the value of integrating seismic and petrological analyses to enhance our understanding of volcanic plumbing systems. The identified velocity anomalies provide critical evidence of the magmatic system's geometry and physical properties, highlighting the interplay between crustal and upper mantle structure and magmatic processes. These findings contribute to the ongoing monitoring efforts and hazard assessment of Campi Flegrei.
How to cite: Ortega-Ramos, V., D'Auria, L., Granja-Bruña, J. L., Cabrera-Pérez, I., Pappalardo, L., Buono, G., and Pérez, N.: Integrating seismic and petrological data unveils the magmatic system beneath Campi Flegrei Caldera., EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-14066, https://doi.org/10.5194/egusphere-egu25-14066, 2025.
Terceira is a volcanic Island of the Azores Archipelago, located in a tectonic triple junction expression of the intersection of the Mid-Atlantic Ridge with the complex western sector of the Eurasian-Africa Boundary. The island sits on the so-called Terceira rift, the Azorean segment of the triple junction, oblique to the main spreading direction with ultra-slow spreading rate.
The main active volcanic structure is the Santa Barbara volcano to the west of the island, geologically the younger part. The island is divided by a roughly WNW-ESE fissural volcanic system, that intersects Santa Bárbara volcano to the north, crosses the active Guilherme Moniz volcano, in the mid of the island, and the extinct Serra do Cume volcano to the south. Pico Alto, another active volcano, is located north of the fissural system, whereas to the NE lies the significant tectonic structure of the Lajes Graben.
We will present a 3D tomographic model of Terceira Island, obtained from coupled inversion of local earthquake and ambient noise tomography, using a dense temporary seismic deployment (2019-2020) coupled with data from the permanent IPMA Seismic Network (2000-2019).
The Santa Bárbara volcano presents a distinct tomographic signature, with a high Vp and Vp/Vs anomalies located beneath the caldera and extending ~4-6 km, and significant seismicity up to ~6 km depths. On the other hand, although with the same signal the Vp and Vp/Vs anomalies beneath the Pico Alto volcano are much weaker, with few seismicity associated. The shallow layers beneath Guilherme Moniz present weak perturbations on both Vp and Vp/Vs, but a stronger deep high Vp anomaly is clear around 7 km depths with no significant Vp/Vs variation associated. The fissural system signature is usually associated with shallow anomalies, with high Vp north of Santa Bárbara volcano, but low Vp and very low Vp/Vs ~1.6 beneath the area between Pico Alto and Guilherme Moniz volcanoes; the low Vp/Vs anomaly in the center of the island is very shallow <3 km and is probably due to the presence and migration of geofluids in the vicinity of the geothermal powerplant.
The fissural system in this area presents a strong seismic activity, with the deepest events recorded by the networks ~9 km; a relatively high number of very shallow and low magnitude events were recorded around the area of the geothermal powerplant.
The SW area of the Serra do Cume volcano presents high Vp values and usually a normal Vp/Vs ~1.77, suggesting an already frozen volcanic system. It has some significant seismicity associated probably due to the intersection of the fissural system that splits the island. To the NE, a bit surprisingly, no significant seismicity was recorded along the Lajes Graben; the lack of seismicity also gave a low tomographic coverage.
This is a contribution to RESTLESS (DOI:10.54499/PTDC/CTA-GEF/6674/2020) and GEMMA (DOI: 10.54499/PTDC/CTA-GEO/2083/2021) and funded by the Portuguese Fundação para a Ciência e a Tecnologia (FCT) I.P./MCTES through national funds (PIDDAC): UID/50019/2025 and LA/P/0068/2020 https://doi.org/10.54499/LA/P/0068/2020), and UIDB/04683 and UIDP/04683 – Instituto de Ciências da Terra.
How to cite: Dias, N. A., Fontiela, J., Silveira, G., Moreira, M., and Matias, L.: Terceira island, Azores: crustal imaging and correlation of volcano-tectonic structures with seismicity, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-12383, https://doi.org/10.5194/egusphere-egu25-12383, 2025.
The Soufrière Hills Volcano on the island of Montserrat in the Lesser Antilles was activated in 1995 after centuries of dormancy. Since then and until 2010, five major phases of activity led to explosive eruptions, dome building and collapse events, pyroclastic flows and ash clouds, which buried the capital, Plymouth, caused loss of life and extensive destruction of infrastructure, rendering the southern part of the island uninhabitable. Montserrat is currently rebuilding, with efforts focusing among others in the importance of better understanding the volcanic system and its associated geothermal field for both risk management and potential resource development.
In this work we jointly invert 7,112 P-wave and 1,376 S-wave arrival times recorded on 18 seismic stations from 1039 local events, classified by the Montserrat Volcano Observatory as tectonic or volcano-tectonic earthquakes, to determine the three-dimensional compressional and shear wave velocities and to simultaneously improve the event locations. Finite frequency sensitivity kernels are used instead of geometrical rays-paths to account for finite frequency effects, providing a more realistic representation of data’s sensitivity. To evaluate the quality and the uncertainty of the tomographic images, we utilize high performance direct sparse algorithms and graph-theory based techniques that allow the efficient calculation of the model’s resolution and posterior covariance matrices, allowing to assess the robustness of the features revealed in the tomographic models.
Our results show significant velocity anomalies associated with the volcanic centres and the associated geothermal field. Fast seismic P- and S- wave velocities are imaged beneath the active volcano of Soufrière Hills and the older dormant Central Hills volcanic centre. The fast velocity anomaly becomes more prominent and wider beneath Soufrière Hills at depths approximately between 1 and 3 km below sea level. These regions correspond to present and older volcanic cores and are possibly comprised from andesitic crystalized rocks of dome cores and intrusive magmatic bodies such as dikes and sills. The surrounding regions appear slower, possibly associated with deposits of volcanic ash, lava fragments, pyroclastic flows and lahars. A significant P- and S-wave low velocity anomaly is observed to the western side of the Island to the area of a high-temperature geothermal field that was inferred from previous geophysical studies and confirmed by three geothermal wells. In our models this feature is dipping towards E/SE, and extends down to ~3.5 km. These results are in good agreement with previous studies in the region. The Vp/Vs ratio map derived from the two velocity models reveals low Vp/Vs ratio at the volcanic centre, and positive anomalies at the flanks of Soufrière Hills. The low-velocity region that is associated with the geothermal field is characterized by low Vp/Vs ratio values too.
How to cite: Bogiatzis, P., Kendall, J.-M., Baird, A. F., Blundy, J. D., Stewart, R. C., and Ryan, G. A.: Insights on the structure of the Soufrière Hills volcano and the associated geothermal field using joint inversion of passive P- and S-wave travel time data, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-19072, https://doi.org/10.5194/egusphere-egu25-19072, 2025.
In this study, we develop a multiparameter ambient noise adjoint tomography method, and apply it for the first time to image the crustal magmatic system and regional deformation of the Toba region. Using Rayleigh and Love waves at periods of 5-20 s extracted from ambient noise, we construct a new multiparameter 3D crustal model that includes shear-wave velocity, radial and azimuthal anisotropy. The isotropic component of our model reveals 1) over 30% Vs reductions beneath the Toba caldera with a melt fraction ranging from 14.5% to 18.5%, and 2) two low Vs bodies located in the middle crust (10-20 km) beneath the Helatoba volcano and the upper crust beneath the Lubukraya volcano, suggesting a large transcrustal magmatic mush model beneath this volcanic arc region. Our anisotropic model shows > 10% positive radial anisotropy (Vsh > Vsv) in the middle crust of the volcanic regions, indicating the presence of horizontally layered melt sills. In the upper crust, we find predominantly weak negative radial anisotropy and significant azimuthal anisotropy, suggesting subvertical rock fabrics dominate upper crustal anisotropy. The orientation of fast velocity directions (FVDs) mostly aligns with the Sumatran Fault due to fault fabrics resulting from shearing deformation along the plate boundary. In the Tarutung region with rich geothermal resources, FVDs shift to being fault-perpendicular probably due to the alignment of stress-induced, fluid-rich microcracks. Our study provides new insights into crustal magmatic architecture and deformation regimes of the Toba region shaped by regional tectonics and magmatic processes.
How to cite: Wang, K., Tong, P., Zhang, Z., and Su, L.: Mapping the crustal magmatic system and regional deformation of the Toba region by multiparameter ambient noise adjoint tomography, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-5325, https://doi.org/10.5194/egusphere-egu25-5325, 2025.
Changbaishan volcano (CBV), located on the border between China and North Korea, has experienced multiple large-scale eruptions in the past. The well-known millennial eruption occurred in 946 and is one of the largest eruptions in the world. In 2002~2005, the CBV experienced an episode of unrest with intensive seismicity and some ground deformations, which has raised great concerns for the public. Many seismic tomography studies have been conducted in the CBV region, but due to the sparse distribution of seismic stations or smaller spatial coverage of some dense temporary stations, the crustal magma plumbing system for the CBV is still not well resolved.
To characterize the fine crustal structure of the CBV area and characterize its magmatic system, we deployed a dense seismic array over a one-month period, which consists of 277 short-period stations. The array covers an area of 200 km in the east-west direction, and 190 km in the north-south direction, with an average station spacing of 10 km. We also assembled continuous data from 14 permanent volcanic monitoring stations operated by the Jilin Earthquake Agency.
To extract high-quality empirical Green’s functions (EGFs) from one month of continuous ambient noise data, we first followed the conventional ambient noise data processing flow to compute hourly cross-correlation functions (CCFs) for each station pair. Subsequently, the template-matching-based selection method was applied to select CCFs with distinct surface wave signals for some time segments. These selected CCFs were then denoised using the SVD-based wiener filter (SVDWF) to further enhance the signal-to-noise ratios (SNRs) of CCFs. Finally, phase-weighted stacking (PWS) was employed to obtain the final CCFs for each station pair. This processing workflow significantly improved the SNRs of EGFs, enhanced the quality of the dispersion spectra and extended the surface wave dispersion frequency band from ~10 s to ~20 s.
After extracting high-quality dispersion curves, we employed the direct surface wave tomography method to invert the crustal velocity structure beneath the CBV region. The inverted Vs model reveals a prominent low-velocity anomaly in the mid-to-upper crust beneath the Tianchi crater, and the presence of a wide-spread low-velocity layer approximately 10 km thick in the middle crust. In addition, our Vs model also indicates low-velocity anomalies in the upper crust beneath the Wangtian’e crater. A unified low-velocity body connects the Tianchi crater with the Namphothe crater in the mid-to-lower crust, which suggests that these volcanic systems may originate from a common deeper source.
How to cite: Pan, D., Gao, J., and Zhang, H.: Ambient Noise Tomography Reveals Magmatic Plumbing System Beneath the Changbaishan Volcano Field, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-16194, https://doi.org/10.5194/egusphere-egu25-16194, 2025.
Posters on site: Wed, 30 Apr, 16:15–18:00 | Hall X1
In 2023 a comprehensive seismic monitoring experiment across the Aeolian volcanic archipelago was deployed as part of a multinational collaborative effort to investigate the relationship between seismicity and volcanic activity. Covering the three southernmost islands in the southern Tyrrhenian Sea, this study included 120 nodal seismic sensors operating in Fall and Winter 2023. Later, 16 broadband stations were deployed over the southern Aeolian Islands of Alicudi, Filicudi, Salina, Panarea, Lipari and Vulcano, and in the vicinity of Milazzo town, in the northern Sicilian coast yield high-resolution seismic data over a wide frequency range. Preliminary results from spectral analysis reveal distinct seismic signatures associated with volcanic processes, as supported by geochemical observations of subsurface activity (Federico et al., 2023; Inguaggiato et al., 2023). These findings emphasize the importance of dense, multiscale seismic networks and interdisciplinary approaches in advancing volcanic hazard assessment and early warning capabilities.
This study is funded by the INGV Pianeta Dinamico project 2023-2025 CAVEAT (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: Di Luccio, F. and the The CAVEAT Team: The Tyrrhenian Sea under watch: multinational, multiscale seismic monitoring across the Aeolian Islands, southern Italy, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-18022, https://doi.org/10.5194/egusphere-egu25-18022, 2025.
Since magma and gas uprise and accumulation into the crust differently affect seismic velocities, seismic tomography is revealing a potential tool to detect how melt versus exsolved magmatic volatile phase reservoirs are distributed at depth into the volcanic plumbing systems. 4D tomography (in space and time) can detect temporal changes in seismic velocities, enabling to follow gas and magma accumulation, a crucial aspect in monitoring open vent volcanoes and unresting calderas. To achieve this goal we combined fully non-linear tomography with improved and fast seismic phases detection allowed by machine-learning. The use of a fully non-linear approach allows overcoming some limits of standard linearized methods, which can obscure details because of damping, smearing and blurring of seismic anomalies, due to the need of global regularization of the inverse problem. We are developing a Bayesian approach to local earthquake tomography that erases the dependence on arbitrary starting velocity models, providing more reliable absolute velocity values and model uncertainties. This is an important aspect when observed changes in seismic velocity have to be compared with theoretical predictions from petrophysical modeling. Furthermore, a data-driven self-adapting parameterizations of the earth structure strongly enhances the resolution capability in regions where high gradients in velocity are expected and in regions poorly illuminated by seismic rays, permitting to reveal seismic anomalies not detectable by standard approaches because of pre-determined model parametrization. Here, we present our novel approach and the application to the Campi Flegrei caldera, where the discrimination between gas and magma injection at shallow depth is crucial for unraveling the causal process of the unrest and in multirisk assessment and forecasting.
How to cite: Giacomuzzi, G., Chiarabba, C., De Gori, P., Fonzetti, R., and Piana Agostinetti, N.: 4D non-linear tomography of the Campi Flegrei caldera, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-3825, https://doi.org/10.5194/egusphere-egu25-3825, 2025.
Two months of continuous records from 40 three-component nodal stations deployed from late October to late December 2023 are used to reveal the subsurface structure down to 2 km depth beneath the Salina Island, in the Tyrrhenian Sea (southern Italy). We present the first 3D shear wave velocity model of a relatively small island of ~26 km2 areal extension derived by ambient noise tomography. We calculate Green’s functions for the vertical components and the anti-causal and causal components were manually and separately inspected to obtain the group velocity dispersion curves via frequency-time analysis (FTAN). Using the Noisy Dispersion Curve Picking program (Granados-Chavarría et al., 2019) we obtain a total of 616 dispersion curves, from 0.25 to 4 Hz. Between each station-pair, we considered the distance along the relief instead of the great-circle distance, to take into account the abrupt topographic gradients crossing the two cone-shaped volcanoes of Salina, Monte dei Porri to the west and Monte Fossa delle Felci to the east, the highest peak of the entire archipelago.
Salina Island hosts a widely variability of volcanic activities (volcanic fissure, collapses, stratovolcanoes, dyke intrusions and diatremes) in a small area. Our model allowed us to image the most representative volcanic episodes that shaped the island: 1) the oldest and highly eroded Pizzo di Corvo in the westernmost region, 2) the volcanic fissure-type volcanic episodes of Pizzo Capo (northeastern region) and Monte Rivi in central Salina (which transits to a central-type activity), 3) the two main stratovolcanoes, Monte Fossa delle Felci and Monte dei Porri, and 4) the youngest activity at the northwestern region, the collapse of the Pollara depression.
This study is funded by the INGV Pianeta Dinamico project 2023-2025 CAVEAT (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., and Ventura, G.: Ambient Noise Tomography of the Salina volcanic island, southern Italy, from nodal seismic data, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-11567, https://doi.org/10.5194/egusphere-egu25-11567, 2025.
We present a new 3D overall model of Vp, Vs and Vp/Vs for the Mount Etna (southern Italy), the largest and most active volcano in Europe. We applied the LOTOS code (Koulakov, BSSA 2009) to a dataset of ~4600 crustal earthquakes that occurred in the study area during the last 26 years (Totaro et al., SciRep. 2024). The selected dataset, representing the longest time-interval ever analyzed for Mt. Etna, allowed us to characterize the volcano velocity structure getting over possible singularities due to specific eruptive phases. We estimated and jointly interpreted P- and S-wave velocity patterns together with the Vp/Vs ratio, particularly effective to discriminate the presence of groundwater, gas, and melts and thus very precious for volcano investigations (Kuznetsov et al., Geosciences 2017; Vargas et al., SciRep. 2017; Totaro et al., SciRep. 2022). The obtained 3D seismic velocity patterns allowed us to add further details on already known anomalies and to identify new previously undetected ones. Focusing on the latter, at the shallowest layer we highlight the presence of two high Vp/Vs volumes, located in close correspondence with low resistivity areas (Siniscalchi et al., JVGR 2010, JGR-SE 2012), that can be associated to underground aquifers generated by meteoric water penetrating the volcano edifice. Moreover, a high Vp/Vs anomaly characterized by intense seismic activity has been clearly detected along the eastern flank of Mt. Etna representing a volume of strongly fractured sedimentary rocks through which a large amount of fluids may rise. Finally, on the western side, a high Vp/Vs area with very low seismicity is detectable. The achieved velocity patterns may suggest fluid accumulation, probably not associated to the volcanic activity, even if further investigations are necessary to better solve and understand this previously unknown anomalous region. In conclusion, our study furnished a comprehensive velocity model that, encompassing specific volcanic phases and allowing a joint interpretation of Vp, Vs and Vp/Vs patterns, provides a more complete modelling of the main features of Mt. Etna.
How to cite: Totaro, C., Aloisi, M., Ferlito, C., Orecchio, B., Presti, D., and Scolaro, S.: 3D Local Earthquake Tomography provides new insights on Mount Etna velocity structure, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-16041, https://doi.org/10.5194/egusphere-egu25-16041, 2025.
The creation of seismic event catalogs has been revolutionized by the availability of continuous waveform data and the integration of deep learning algorithms for tasks such as event detection, phase picking, association, event localization, and classification. In this study, we showcase the application of these advanced methodologies to generate comprehensive "deep" seismic catalogs for volcanic regions, focusing on the Canary Islands. We also demonstrate how these enhanced catalogs contribute to seismic tomography studies.
Our first analysis evaluates the performance of deep learning-based phase pickers when applied to volcano-tectonic events. These pickers, characterized by minimal parameter tuning requirements (typically only a probability threshold for valid picks) offer a significant advantage. However, as they are primarily trained on datasets lacking volcanic events, their sensitivity to such earthquakes may be reduced, and false positives could be more frequent. To address these challenges, we propose a robust workflow combining deep learning-based phase picking, event association, and relocation. This approach yields seismic catalogs that are more complete and accurate compared to those generated using conventional methods.
Finally, we utilize these improved seismic catalogs to construct 3D P- and S-wave velocity models for regions within the Canary Islands, including La Palma, as well as the central archipelago's regional structure. These models provide new insights into the subsurface dynamics of this volcanic system.
How to cite: Villaseñor, A., Díaz Suárez, E., Domínguez-Cerdeña, I., del Fresno, C., and Bartolomé, R.: Seismic Catalogs and Tomographic Velocity Models for the Canary Islands: A Deep Learning Approach, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-10257, https://doi.org/10.5194/egusphere-egu25-10257, 2025.
Popocatépetl, one of Mexico's most active and hazardous volcanoes, threatens millions living in its shadow. Since its reactivation in 1994, several seismic tomography studies have been conducted. However, their results show significant inconsistencies, fueling debate about the volcano's internal structure. This study presents new findings from an enhanced local earthquake tomography with data from 2019 to 2024.
We used a seismic network comprising up to 18 broadband stations. Ten of these stations belong to the National Center for Disaster Prevention (CENAPRED) and form the volcano’s permanent monitoring network. The remaining eight are managed by us at the Institute of Geophysics at UNAM. The spatial distribution of these stations ensures adequate azimuthal coverage of the volcanic edifice and its seismicity, enabling the highest depth resolution achieved to date.
The local earthquake database was created using a machine learning-based workflow. This involved automatic phase picking and association, followed by post-processing to refine the VT catalog. This approach identified additional events absent from CENAPRED’s official catalog, improving ray coverage across the region. We then applied Enhanced Seismic Tomography, combining Double-Difference Tomography with the Weighted Average Model post-processing to minimize bias from initial parameters.
Over 1,000 events were jointly inverted to produce three-dimensional P- and S-wave velocity models and relocations. Resolution tests indicate the ability to resolve structures of approximately 1×1×0.5 km, offering unprecedented detail of the volcano's features. Finally, we compare our results with previous models to refine our understanding of Popocatépetl’s interior.
This research was supported by the UNAM-PAPIIT Program: IN103823.
How to cite: Bernal-Manzanilla, K. and Calò, M.: Enhanced Local Earthquake Tomography at the Popocatépetl Volcano (Mexico), EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-734, https://doi.org/10.5194/egusphere-egu25-734, 2025.
Understanding the complex subsurface structures in volcanic regions is crucial for effective hazard assessment and resource management. In such areas, traditional active seismic studies often require substantial human and material resources, while operating in densely populated areas. In the contrary, passive seismic data offers a valuable alternative by making use of the extensive seismic records accumulated over time, but introduces complexities such as unknown source location, magnitude and fault mechanism. This study aims to develop a workflow for processing passive-source seismic data with irregular source-receiver distributions, identifying reflection/converted phases to produce a seismic profile suitable for seismic reflection imaging. To achieve this, we simulated the Campi Flegrei caldera case-study with 400 micro-earthquakes, randomly distributed between 1 and 3 km in depth, and recorded by 13 stations. We explore different gather configurations and evaluate their effectiveness for seismic imaging. We apply several preprocessing steps to both synthetic and real seismic data to ensure high-quality imaging results. Reflection phases are carefully identified by utilizing the predicted theoretical arrival times of both reflected and converted phases. Finally, through this processing procedure, we ultimately generate a passive-source seismic profile with irregular source-receiver distributions, that is exploitable for Kirchoff migration imaging.
How to cite: Zhou, X., Nazeri, S., Zand, T., Virieux, J., and Zollo, A.: Optimizing Passive Seismic Data Processing in Volcanic Areas with Irregular Source-Receiver Distributions , EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-21232, https://doi.org/10.5194/egusphere-egu25-21232, 2025.
The Christiana-Santorini-Kolumbo volcanic field (CSKVF) in the Aegean Sea represents one of Europe's most active volcanic centers. Over the past several hundred thousand years, the volcanic evolution of Santorini has encompassed a wide range of eruptive styles, ranging from catastrophic caldera-forming events (exemplified by the Minoan and Cape Riva eruptions) to explosive submarine episodes (such as the 726 CE Kameni and 1650 CE Kolumbo eruptions) and predominantly effusive activity that has characterized the Kameni islands in recent centuries. Due to this remarkable variability, Santorini represents one of the most prominent natural laboratories for volcanological research and education worldwide. However, the majority of deposits from these eruptions are preserved offshore around Santorini, making them inaccessible for direct examination and sampling. In 2019 and 2024, we acquired over 1800 km of high-resolution 2D seismic reflection profiles within and around Santorini to study the volcano-tectonic evolution of the CSKVF, complemented by drilling results from IODP Expedition 398. This unique dataset allows us to directly correlate seismic units with specific eruptions, enabling us to identify characteristic seismic signatures for various volcanic deposits and to constrain their volumes and emplacement dynamics. This integrated approach allows us to fill critical gaps in the eruptive record and develop a comprehensive catalogue linking seismic facies to volcanoclastic deposit types, which can serve as an analogue for interpreting marine seismic data from less well-documented volcanic regions.
How to cite: Karstens, J., Preine, J., Blanch Jover, M., Berndt, C., Crutchley, G. J., Nomikou, P., and Kutterolf, S.: From reflections to eruptions – reconstructing volcanic eruptions using marine seismic data from offshore Santorini, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-9391, https://doi.org/10.5194/egusphere-egu25-9391, 2025.
Submarine explosive volcanism poses significant hazards to coastal communities, infrastructures, and marine and air traffic. However, our understanding of the mechanisms, frequencies, and distributions of submarine explosive eruptions is limited due to their inaccessibility. This observational gap is particularly critical for shallow submarine environments, where the interplay between magma and seawater can lead to violent explosive eruptions that may trigger destructive tsunamis, pyroclastic surges, extensive pumice rafts, and large airborne ash plumes. In this study, we use high-resolution seismic, bathymetric, and seafloor imagery to investigate the formation mechanism of submarine volcanoes at the Northern Reykjanes Ridge. Our seismic images reveal distinctive volcanic edifices characterized by low width-height ratios, stratified outward-dipping reflections, and extensive volcanoclastic aprons. These overly a glacial erosion unconformity, indicating Holocene formation. We show that the volcanoes of the northern Reykjanes Ridge formed predominantly by shallow submarine eruptions, some of which breached the sea-surface forming short-lived islands that were historically observed.
By comparing the seismic reflection patterns from these volcanoes with submarine volcanoes from the Azores and Aegean, we establish key characteristic to distinguish three distinct submarine volcanic formation modes: (1) deep-water explosive eruptions, (2) shallow-water explosive eruptions, and (3) Surtseyan eruptions. Our study highlights the potential hazards of future eruptions along the Reykjanes Ridge, which may include tsunamis, ash plumes, or extensive pumice rafts. By establishing seismic signatures of submarine explosive eruptions, our study highlights the potential of seismic imaging as a powerful tool for understanding submarine volcanic processes, providing key insights into the formation and evolution of submarine volcanoes in marine environments.
How to cite: Preine, J., Hübscher, C., Pałgan, D., van der Zwan, F. M., Dittmers, C., Friedrich, A., Beethe, S., Ehlies, V., Ford, J., Haimerl, B., Ischebeck, L., Lackner, M., Schmidt, M., Eisermann, J. O., Budke, L., Óðinsson, D. Þ., and Augustin, N.: Seismic Imaging of Submarine Volcanoes at the Northern Reykjanes Ridge, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-1433, https://doi.org/10.5194/egusphere-egu25-1433, 2025.
