Pub Date : 2025-10-17DOI: 10.1016/j.jvolgeores.2025.108472
Wandono Wandono , Mohamad Ramdhan , Atin Nur Aulia , Arya Bani Pangestu , Edi Hidayat , Syuhada Syuhada , Sri Widiyantoro , Nicholas Rawlinson , Titi Anggono , Bayu Pranata , Muhammad Maruf Mukti , Mufti Putra Maulana
The western Java region is situated along the Sunda Arc, which is formed by subduction of the Indo-Australian Plate beneath the Eurasian Plate. The subduction process causes seismic activity in the Wadati-Benioff zone of the descending slab, produces a chain of active volcanoes, and forms numerous faults in the terrestrial and back-arc region of the overriding plate. This study examines the region's tectonic system via seismic tomography, in which the most recent earthquake catalog from the Indonesian Agency for Meteorology, Climatology, and Geophysics (BMKG) is utilized. Body wave travel time from 4561 local-regional earthquakes recorded by 78 seismic stations produces a 3-D Vp and Vp/Vs ratio model. The subducting Indo-Australian Plate slab is characterized mainly by high Vp and low Vp/Vs anomalies and exhibits steeper subduction than the Slab2 model. Partial melting zones, which serve as magma sources for several volcanoes, are observed at depths of about 90–100 km, characterized by low Vp and high Vp/Vs anomalies. The Lembang and Cimandiri fault lineaments are visible at 10 km depth and are characterized by negative Vp and positive Vp/Vs anomalies. This study also successfully demonstrates that the source of the damaging 2022 Cianjur earthquake occurred in a region dominated by low Vp and high Vp/Vs. These anomalies most likely relate to fluids or molten material in the fault zone. This study again highlights the importance of the BMKG seismic network for understanding tectonic systems on a local-regional scale, in this case, through the lens of seismic tomography.
{"title":"Seismic velocity structure beneath the western Java region, Indonesia, from local earthquake tomography","authors":"Wandono Wandono , Mohamad Ramdhan , Atin Nur Aulia , Arya Bani Pangestu , Edi Hidayat , Syuhada Syuhada , Sri Widiyantoro , Nicholas Rawlinson , Titi Anggono , Bayu Pranata , Muhammad Maruf Mukti , Mufti Putra Maulana","doi":"10.1016/j.jvolgeores.2025.108472","DOIUrl":"10.1016/j.jvolgeores.2025.108472","url":null,"abstract":"<div><div>The western Java region is situated along the Sunda Arc, which is formed by subduction of the Indo-Australian Plate beneath the Eurasian Plate. The subduction process causes seismic activity in the Wadati-Benioff zone of the descending slab, produces a chain of active volcanoes, and forms numerous faults in the terrestrial and back-arc region of the overriding plate. This study examines the region's tectonic system via seismic tomography, in which the most recent earthquake catalog from the Indonesian Agency for Meteorology, Climatology, and Geophysics (BMKG) is utilized. Body wave travel time from 4561 local-regional earthquakes recorded by 78 seismic stations produces a 3-D Vp and Vp/Vs ratio model. The subducting Indo-Australian Plate slab is characterized mainly by high Vp and low Vp/Vs anomalies and exhibits steeper subduction than the Slab2 model. Partial melting zones, which serve as magma sources for several volcanoes, are observed at depths of about 90–100 km, characterized by low Vp and high Vp/Vs anomalies. The Lembang and Cimandiri fault lineaments are visible at 10 km depth and are characterized by negative Vp and positive Vp/Vs anomalies. This study also successfully demonstrates that the source of the damaging 2022 Cianjur earthquake occurred in a region dominated by low Vp and high Vp/Vs. These anomalies most likely relate to fluids or molten material in the fault zone. This study again highlights the importance of the BMKG seismic network for understanding tectonic systems on a local-regional scale, in this case, through the lens of seismic tomography.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"468 ","pages":"Article 108472"},"PeriodicalIF":2.3,"publicationDate":"2025-10-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145333337","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-16DOI: 10.1016/j.jvolgeores.2025.108473
Amador Castro-Colín , Antonio Pola , Hugo Sereno , Juan Daniel Pérez-Orozco , Martha Gabriela Gómez-Vasconcelos , Elia Mercedes Alonso-Guzmán
Ignimbrite welding encompasses processes such as glass sintering, deposit compaction, and pumice fragment flattening, influenced by emplacement temperature, cooling rate, and viscosity. These processes significantly impact the physical properties, mechanical behavior, and textural characteristics of ignimbrites. This study examines the Escalera Ignimbrite in Central Mexico to evaluate the welding degree using physical and mechanical properties alongside ultrasonic wave velocities and dynamic elastic constants. The analysis reveals stratigraphic variations in bulk density, porosity, and oblateness (OB) of pumice fragments, with increased compaction and welding intensity observed at greater depths. Lower stratigraphic sections exhibit higher bulk density (2.26 g/cm3) and reduced porosity (16.08 %), correlating with increased uniaxial compressive strength (UCS) values ranging from 60 to 65 MPa. Similarly, splitting tensile strength (σt) increases with depth, reaching values of 13.3 to 16.3 MPa in highly welded sections. Dynamic elastic modulus (e.g., Ed) varies from 3.97 in upper sections to 8.91 GPa in lower (denser layers), reflecting enhanced compaction and reduced porosity. Additionally, P-wave (Vp) and S-wave (Vs) velocities exhibit trends consistent with welding intensity, with values ranging from 1673 to 2497 m/s for Vp and from 887 to 1216 m/s for Vs, suggesting their reliability as welding indicators compared to density and porosity alone. The results highlight a robust correlation between acoustic and mechanical properties with welding degree, emphasizing the role of post-depositional processes in modifying rock properties. This study contributes to the understanding of the mechanical evolution of ignimbrites and the implications of welding on their behavior.
{"title":"Welding degree assessment of the Escalera Ignimbrite in Central Mexico: Insights from physical and mechanical properties","authors":"Amador Castro-Colín , Antonio Pola , Hugo Sereno , Juan Daniel Pérez-Orozco , Martha Gabriela Gómez-Vasconcelos , Elia Mercedes Alonso-Guzmán","doi":"10.1016/j.jvolgeores.2025.108473","DOIUrl":"10.1016/j.jvolgeores.2025.108473","url":null,"abstract":"<div><div>Ignimbrite welding encompasses processes such as glass sintering, deposit compaction, and pumice fragment flattening, influenced by emplacement temperature, cooling rate, and viscosity. These processes significantly impact the physical properties, mechanical behavior, and textural characteristics of ignimbrites. This study examines the Escalera Ignimbrite in Central Mexico to evaluate the welding degree using physical and mechanical properties alongside ultrasonic wave velocities and dynamic elastic constants. The analysis reveals stratigraphic variations in bulk density, porosity, and oblateness (OB) of pumice fragments, with increased compaction and welding intensity observed at greater depths. Lower stratigraphic sections exhibit higher bulk density (2.26 g/cm<sup>3</sup>) and reduced porosity (16.08 %), correlating with increased uniaxial compressive strength (UCS) values ranging from 60 to 65 MPa. Similarly, splitting tensile strength (<em>σ</em><sub><em>t</em></sub>) increases with depth, reaching values of 13.3 to 16.3 MPa in highly welded sections. Dynamic elastic modulus (<em>e.g.</em>, <em>E</em><sub><em>d</em></sub>) varies from 3.97 in upper sections to 8.91 GPa in lower (denser layers), reflecting enhanced compaction and reduced porosity. Additionally, P-wave (<em>Vp</em>) and S-wave (<em>Vs</em>) velocities exhibit trends consistent with welding intensity, with values ranging from 1673 to 2497 m/s for <em>Vp</em> and from 887 to 1216 m/s for <em>Vs</em>, suggesting their reliability as welding indicators compared to density and porosity alone. The results highlight a robust correlation between acoustic and mechanical properties with welding degree, emphasizing the role of post-depositional processes in modifying rock properties. This study contributes to the understanding of the mechanical evolution of ignimbrites and the implications of welding on their behavior.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"468 ","pages":"Article 108473"},"PeriodicalIF":2.3,"publicationDate":"2025-10-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145333242","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-15DOI: 10.1016/j.jvolgeores.2025.108474
Daniel Díaz , Maximiliano Pavez , Álvaro Amigo
Villarrica volcano is located in the Southern Andes, one of the world's most productive volcanic arcs of the last decades and stands as the most active volcano in South America, in terms of the number of registered eruptive episodes. Besides its large eruptive record, Villarrica has other particularities, such as being currently an open conduit volcano providing a pathway through which volatiles and magma rise to the surface forming a persistent lava lake at its crater.
These features make Villarrica one of the most relevant volcanoes to study along the Andes, and efforts to characterize its internal magmatic structure have arisen during recent years. Along with geochemical and geophysical measurements normally used to image the magmatic system of this volcano, the use of monitoring methods to analyze different changes in time are needed to improve the understanding of active volcanic systems such this one.
Through a research project to test the capabilities of magnetotellurics as a volcano monitoring tool, a first permanent broad-band magnetotelluric station was installed at Villarrica volcano during January 2023, and has been measuring since then. Even considering the limitations of a single station experiment, interesting changes in magnetotelluric parameters measured during 2023 and beginning of 2024 have been estimated, and compared to other geophysical and geological parameters measured during this time, aiming to improve our monitoring capabilities and provide new insights into the magmatic processes at one of the most active volcanoes in the Andes.
{"title":"First year of magnetotelluric observations at Villarrica volcano, Southern Andes","authors":"Daniel Díaz , Maximiliano Pavez , Álvaro Amigo","doi":"10.1016/j.jvolgeores.2025.108474","DOIUrl":"10.1016/j.jvolgeores.2025.108474","url":null,"abstract":"<div><div>Villarrica volcano is located in the Southern Andes, one of the world's most productive volcanic arcs of the last decades and stands as the most active volcano in South America, in terms of the number of registered eruptive episodes. Besides its large eruptive record, Villarrica has other particularities, such as being currently an open conduit volcano providing a pathway through which volatiles and magma rise to the surface forming a persistent lava lake at its crater.</div><div>These features make Villarrica one of the most relevant volcanoes to study along the Andes, and efforts to characterize its internal magmatic structure have arisen during recent years. Along with geochemical and geophysical measurements normally used to image the magmatic system of this volcano, the use of monitoring methods to analyze different changes in time are needed to improve the understanding of active volcanic systems such this one.</div><div>Through a research project to test the capabilities of magnetotellurics as a volcano monitoring tool, a first permanent broad-band magnetotelluric station was installed at Villarrica volcano during January 2023, and has been measuring since then. Even considering the limitations of a single station experiment, interesting changes in magnetotelluric parameters measured during 2023 and beginning of 2024 have been estimated, and compared to other geophysical and geological parameters measured during this time, aiming to improve our monitoring capabilities and provide new insights into the magmatic processes at one of the most active volcanoes in the Andes.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"468 ","pages":"Article 108474"},"PeriodicalIF":2.3,"publicationDate":"2025-10-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145333269","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-12DOI: 10.1016/j.jvolgeores.2025.108470
Ahmed Baamar , Ezzoura Errami , Brahim Karaoui , Zakarya Yajioui , John Browning
The study of dyke swarms is crucial for understanding the processes of magma intrusion and the tectonic environments that influence volcanic activity. Dyke swarms are essential records of tectonic and magmatic events, offering insights into magma ascent, chamber dynamics, and stress conditions during emplacement. In this context, the dyke swarm of the central part of the Saghro Massif (Eastern Anti-Atlas) is investigated for the first time, using structural aspects of the dykes and paleo stress inversion to constrain magma overpressure, the depth of magma reservoir, the state of stress during emplacement, and the tectonic regime responsible for emplacement. We selected seventy-eight mafic and intermediate dykes in the swarm of the Central Part of the Saghro Massif (CPSM), characterized by different orientations N-S to NNE-SSW, NE-SW, E-W and NW-SE. In addition, the size distribution of thicknesses and lengths follows a power-law and a log-normal distribution respectively. We calculate magmatic overpressure using selected dyke aspect ratios to estimate the depth of the magma reservoir. This study suggests that the CPSM dykes form mainly due to the injection of magma from a deep magma reservoir at depths >31 km, hence close to the Moho crust-mantle boundary (31–33 km). Paleo stress reconstruction shows that the CPSM dyke swarm was emplaced when the minimum principal compressive stress (σ3) was oriented WNW-ESE, and the maximum principal compressive stress (σ1) was vertical. The reconstruction also shows that the emplacement of the CPSM swarm occurred during a transtensional tectonic regime, associated with the WACadomian orogeny.
{"title":"Structural parameters of the late Ediacaran dyke swarm in the Saghro Massif (Eastern Anti-Atlas, Morocco): Implications for the depth of magma origin, paleo stress reconstruction and emplacement models","authors":"Ahmed Baamar , Ezzoura Errami , Brahim Karaoui , Zakarya Yajioui , John Browning","doi":"10.1016/j.jvolgeores.2025.108470","DOIUrl":"10.1016/j.jvolgeores.2025.108470","url":null,"abstract":"<div><div>The study of dyke swarms is crucial for understanding the processes of magma intrusion and the tectonic environments that influence volcanic activity. Dyke swarms are essential records of tectonic and magmatic events, offering insights into magma ascent, chamber dynamics, and stress conditions during emplacement. In this context, the dyke swarm of the central part of the Saghro Massif (Eastern Anti-Atlas) is investigated for the first time, using structural aspects of the dykes and paleo stress inversion to constrain magma overpressure, the depth of magma reservoir, the state of stress during emplacement, and the tectonic regime responsible for emplacement. We selected seventy-eight mafic and intermediate dykes in the swarm of the Central Part of the Saghro Massif (CPSM), characterized by different orientations N-S to NN<em>E</em>-SSW, NE-SW, E-W and NW-SE. In addition, the size distribution of thicknesses and lengths follows a power-law and a log-normal distribution respectively. We calculate magmatic overpressure using selected dyke aspect ratios to estimate the depth of the magma reservoir. This study suggests that the CPSM dykes form mainly due to the injection of magma from a deep magma reservoir at depths >31 km, hence close to the Moho crust-mantle boundary (31–33 km). Paleo stress reconstruction shows that the CPSM dyke swarm was emplaced when the minimum principal compressive stress (σ<sub>3</sub>) was oriented WNW-ESE, and the maximum principal compressive stress (σ<sub>1</sub>) was vertical. The reconstruction also shows that the emplacement of the CPSM swarm occurred during a transtensional tectonic regime, associated with the WACadomian orogeny.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"468 ","pages":"Article 108470"},"PeriodicalIF":2.3,"publicationDate":"2025-10-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145333268","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-11DOI: 10.1016/j.jvolgeores.2025.108465
Darren Tan , David Fee , Pavel Izbekov , Taryn Lopez , Társilo Girona , Vanesa Burgos , Stephen McNutt , Matthew Haney , Valerie Wasser , Jessica Larsen , Ronni Grapenthin , Mario Angarita , Pablo Saunders-Shultz , Tara Shreve , Jamshid Moshrefzadeh
Pavlof Volcano, a frequently active volcano in the Aleutian arc, has erupted six times between 2007 and 2022. Many Pavlof Volcano eruptions lack detectable precursory seismicity or ground deformation, making them challenging to forecast. Since 2007, these eruptions have ranged from Strombolian to Subplinian (VEI 2–3), with their characteristics seemingly related to their active vent location. Summit vent eruptions tend to be abrupt and explosive, while southeast flank vent eruptions tend to be more effusive and longer-lived. To better understand Pavlof Volcano’s plumbing system and improve eruption forecasts, we integrate geophysical, petrologic, and satellite-based thermal infrared and gas data to conduct a multidisciplinary data synthesis for eruptions between 2007 and 2022. Notably, a recently developed volcano seismicity detection model reveals unique pre- and syn-eruptive seismic tremor regimes associated with each vent system. The summit vent eruptions show greater tremor diversity and resurgent eruptive behavior, whereas the southeast flank vent eruptions show a steady tremor transition over a single eruptive phase. We use local infrasound data and air-to-ground coupled waves to build explosion catalogs, revealing gas-rich explosions during the 2021–2022 eruption from the summit vent, which was previously thought to be inactive during that eruption. Whole-rock composition and ash analyses indicate a near-uniform basaltic andesite composition spanning multiple eruptions, and a consistent anomalous presence of altered olivines in erupted material. We use these findings to propose a new conceptual model for Pavlof Volcano’s plumbing system: a shallow T-junction outlet controlling seismic tremor diversity and the partitioning of gas-charged magma; and an elongated, heated conduit transporting magma from depth and storing partly molten, residual magma between eruptions. We suggest that magma ascent rate and vent sealing modulate eruption style at Pavlof Volcano, where the T-junction directs rapidly ascending, gas-rich magma to the summit vent, and diverts slower rising, gas-poor magma towards the southeast flank vent.
{"title":"Vent-specific unrest at Pavlof Volcano, Alaska: Insights from multidisciplinary data","authors":"Darren Tan , David Fee , Pavel Izbekov , Taryn Lopez , Társilo Girona , Vanesa Burgos , Stephen McNutt , Matthew Haney , Valerie Wasser , Jessica Larsen , Ronni Grapenthin , Mario Angarita , Pablo Saunders-Shultz , Tara Shreve , Jamshid Moshrefzadeh","doi":"10.1016/j.jvolgeores.2025.108465","DOIUrl":"10.1016/j.jvolgeores.2025.108465","url":null,"abstract":"<div><div>Pavlof Volcano, a frequently active volcano in the Aleutian arc, has erupted six times between 2007 and 2022. Many Pavlof Volcano eruptions lack detectable precursory seismicity or ground deformation, making them challenging to forecast. Since 2007, these eruptions have ranged from Strombolian to Subplinian (VEI 2–3), with their characteristics seemingly related to their active vent location. Summit vent eruptions tend to be abrupt and explosive, while southeast flank vent eruptions tend to be more effusive and longer-lived. To better understand Pavlof Volcano’s plumbing system and improve eruption forecasts, we integrate geophysical, petrologic, and satellite-based thermal infrared and gas data to conduct a multidisciplinary data synthesis for eruptions between 2007 and 2022. Notably, a recently developed volcano seismicity detection model reveals unique pre- and syn-eruptive seismic tremor regimes associated with each vent system. The summit vent eruptions show greater tremor diversity and resurgent eruptive behavior, whereas the southeast flank vent eruptions show a steady tremor transition over a single eruptive phase. We use local infrasound data and air-to-ground coupled waves to build explosion catalogs, revealing gas-rich explosions during the 2021–2022 eruption from the summit vent, which was previously thought to be inactive during that eruption. Whole-rock composition and ash analyses indicate a near-uniform basaltic andesite composition spanning multiple eruptions, and a consistent anomalous presence of altered olivines in erupted material. We use these findings to propose a new conceptual model for Pavlof Volcano’s plumbing system: a shallow T-junction outlet controlling seismic tremor diversity and the partitioning of gas-charged magma; and an elongated, heated conduit transporting magma from depth and storing partly molten, residual magma between eruptions. We suggest that magma ascent rate and vent sealing modulate eruption style at Pavlof Volcano, where the T-junction directs rapidly ascending, gas-rich magma to the summit vent, and diverts slower rising, gas-poor magma towards the southeast flank vent.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"468 ","pages":"Article 108465"},"PeriodicalIF":2.3,"publicationDate":"2025-10-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145363370","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Arc volcanoes often show an alternating behaviour between explosive and effusive activity. Although considered an effusive phenomenon, emplacement of highly evolved lava domes is almost always accompanied by explosive activity and processes that initiate the formation of pyroclastic density currents. Lava dome forming eruptions have distinct hazard patterns and whether a specific volcano has produced lava domes during its history is important. While the presence of a lava dome is self-evident for modern (witnessed) eruptions – the question becomes more difficult for deposits in the stratigraphic record. We here review the most recent eruption of Sanbe Volcano, SW Japan, and highlight challenges in the discrimination between explosive and effusive eruption styles. The Taiheizan eruption produced a large number of deposits from pyroclastic density currents (PDC) of the block and ash flow type, leading to the undisputed interpretation that an active lava dome was present during this eruption about 4000 years ago. However, re-evaluation of the deposits suggests they formed by eruption column collapse during a Plinian eruption. The term “block and ash flow deposit” (a descriptive expression for a deposit with a bimodal grainsize distribution), is often tied to an interpretation (formed by lava dome collapse) in the volcanological literature. There are many processes around active volcanoes that can produce deposits with large blocks in a fine-grained matrix and a dominantly bimodal GSD. This makes recognition of dome forming eruptions in the eruptive record challenging but also gives room for misinterpretation. Beside our re-evaluation of the Taiheizan eruption style, we provide a list of features that can help to identify dome forming eruptions in the stratigraphic record.
{"title":"Discrimination between dome-forming and explosive eruptions in the stratigraphic record – field, textural and petrographic evidence from the Taiheizan eruption, Mt. Sanbe, Japan","authors":"Andreas Auer , Keiko Suzuki-Kamata , Tetsuya Kogure , Daisuke Endo , Hiroshi Kitagawa , Shun Orui , Katsura Kobayashi","doi":"10.1016/j.jvolgeores.2025.108462","DOIUrl":"10.1016/j.jvolgeores.2025.108462","url":null,"abstract":"<div><div>Arc volcanoes often show an alternating behaviour between explosive and effusive activity. Although considered an effusive phenomenon, emplacement of highly evolved lava domes is almost always accompanied by explosive activity and processes that initiate the formation of pyroclastic density currents. Lava dome forming eruptions have distinct hazard patterns and whether a specific volcano has produced lava domes during its history is important. While the presence of a lava dome is self-evident for modern (witnessed) eruptions – the question becomes more difficult for deposits in the stratigraphic record. We here review the most recent eruption of Sanbe Volcano, SW Japan, and highlight challenges in the discrimination between explosive and effusive eruption styles. The Taiheizan eruption produced a large number of deposits from pyroclastic density currents (PDC) of the block and ash flow type, leading to the undisputed interpretation that an active lava dome was present during this eruption about 4000 years ago. However, re-evaluation of the deposits suggests they formed by eruption column collapse during a Plinian eruption. The term “block and ash flow deposit” (a descriptive expression for a deposit with a bimodal grainsize distribution), is often tied to an interpretation (formed by lava dome collapse) in the volcanological literature. There are many processes around active volcanoes that can produce deposits with large blocks in a fine-grained matrix and a dominantly bimodal GSD. This makes recognition of dome forming eruptions in the eruptive record challenging but also gives room for misinterpretation. Beside our re-evaluation of the Taiheizan eruption style, we provide a list of features that can help to identify dome forming eruptions in the stratigraphic record.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"468 ","pages":"Article 108462"},"PeriodicalIF":2.3,"publicationDate":"2025-10-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145333271","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-10DOI: 10.1016/j.jvolgeores.2025.108467
Daniele Tardani , Marco Taussi , Franco Tassi , Jorge E. Romero , Philippe Robidoux , Camila Poblete-González , Fernanda Álvarez-Amado , Lorenzo Manosalva , Claudia Berrios , Daniele L. Pinti , Martin Reich , Diego Morata , Valentina Mura , Gloria Arancibia
Gas geochemistry is key to understanding volcanic processes, offering insights into subsurface magma dynamics and aiding eruption forecasting. We present a 10-year monitoring case study from Andean Southern Volcanic Zone (SVZ), one of the world's most active regions, demonstrating its value in assessing volcanic unrest. The geochemical composition of fumarolic gas emissions from the Nevados de Chillán Volcanic Complex (Chile) strongly varied across different phases of volcanic activity, permitting the evaluation of the dynamic interplay between magmatic and hydrothermal processes. During volcanic quiescence periods, (2013 and 2023), fumarolic gases were predominantly controlled by shallow meteoric-hydrothermal circulation, as suggested by Ar (40Ar/36Ar ∼ 270–290) and water isotopic signatures, low 3He/4He (or Rc) ratios (∼3.5 times the atmospheric value Ra), moderate CO2 levels (between ∼2200 and ∼ 9800 μmol/mol), and the absence of magmatic gaseous species (SO2, HCl, and HF). In contrast, the unrest phase (2016 and 2017) was marked by a rapid and significant gas character shift. Rc/Ra values increased up to >6, CO2 concentrations exceeded 12,000 μmol/mol, and magmatic gaseous species became detectable, with SO2, HCl, and HF reaching 8.5, 4.3, and 0.21 μmol/mol, respectively. Additionally, simultaneous Ar (40Ar/36Ar ∼ 370–410) and water isotopic signature shifts from meteoric origin toward a mixing with deep components were also detected. Gas geothermometry, computed through the H2/Ar* – CH4/CO2 and H2/Ar* - CO/CO2 equilibria systems, revealed a rise from 290 ± 10 °C in quiescence to 340 ± 10 °C during the unrest phase. The long-term geochemical surveillance helped detect a substantial signature of volatile-rich magmatic fluids influx into the hydrothermal system before the eruption phase (2018–2021), offering information on the early detection of unrest.
{"title":"Temporal evolution of fumarolic gas geochemistry at the Nevados de Chillán Volcanic complex (2013−2023): Signals of volcanic unrest and insights into the hydrothermal system","authors":"Daniele Tardani , Marco Taussi , Franco Tassi , Jorge E. Romero , Philippe Robidoux , Camila Poblete-González , Fernanda Álvarez-Amado , Lorenzo Manosalva , Claudia Berrios , Daniele L. Pinti , Martin Reich , Diego Morata , Valentina Mura , Gloria Arancibia","doi":"10.1016/j.jvolgeores.2025.108467","DOIUrl":"10.1016/j.jvolgeores.2025.108467","url":null,"abstract":"<div><div>Gas geochemistry is key to understanding volcanic processes, offering insights into subsurface magma dynamics and aiding eruption forecasting. We present a 10-year monitoring case study from Andean Southern Volcanic Zone (SVZ), one of the world's most active regions, demonstrating its value in assessing volcanic unrest. The geochemical composition of fumarolic gas emissions from the Nevados de Chillán Volcanic Complex (Chile) strongly varied across different phases of volcanic activity, permitting the evaluation of the dynamic interplay between magmatic and hydrothermal processes. During volcanic quiescence periods, (2013 and 2023), fumarolic gases were predominantly controlled by shallow meteoric-hydrothermal circulation, as suggested by Ar (<sup>40</sup>Ar/<sup>36</sup>Ar ∼ 270–290) and water isotopic signatures, low <sup>3</sup>He/<sup>4</sup>He (or Rc) ratios (∼3.5 times the atmospheric value Ra), moderate CO<sub>2</sub> levels (between ∼2200 and ∼ 9800 μmol/mol), and the absence of magmatic gaseous species (SO<sub>2</sub>, HCl, and HF). In contrast, the unrest phase (2016 and 2017) was marked by a rapid and significant gas character shift. Rc/Ra values increased up to >6, CO<sub>2</sub> concentrations exceeded 12,000 μmol/mol, and magmatic gaseous species became detectable, with SO<sub>2</sub>, HCl, and HF reaching 8.5, 4.3, and 0.21 μmol/mol, respectively. Additionally, simultaneous Ar (<sup>40</sup>Ar/<sup>36</sup>Ar ∼ 370–410) and water isotopic signature shifts from meteoric origin toward a mixing with deep components were also detected. Gas geothermometry, computed through the H<sub>2</sub>/Ar* – CH<sub>4</sub>/CO<sub>2</sub> and H<sub>2</sub>/Ar* - CO/CO<sub>2</sub> equilibria systems, revealed a rise from 290 ± 10 °C in quiescence to 340 ± 10 °C during the unrest phase. The long-term geochemical surveillance helped detect a substantial signature of volatile-rich magmatic fluids influx into the hydrothermal system before the eruption phase (2018–2021), offering information on the early detection of unrest.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"468 ","pages":"Article 108467"},"PeriodicalIF":2.3,"publicationDate":"2025-10-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145363369","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-10DOI: 10.1016/j.jvolgeores.2025.108468
Rami Alshembari, James Hickey, Lorenzo Mantiloni
Volcano deformation studies provide crucial insights into magmatic processes, offering a key perspective on the dynamics of magma reservoirs and their eruption potential. Here, we benchmark a Finite Element (FE) poroelastic magma mush model with an analytical model of a spherical melt core and mush shell undergoing melt resupply. In our comparisons, we investigate key reservoir parameters, such as melt core pressure, deformation at the core-shell and shell-rock boundaries, and the evolution of tensile stress. The benchmarking provides excellent agreement between the FE results and the analytical solutions. We extend the analysis by incorporating a previously neglected free surface effect in the numerical models, allowing for the examination of surface deformation beyond the original analytical solution. Comparison of the surface deformation resulting from a core-shell magma mush model to that from a uniform poroelastic magma mush model highlights differences in the predicted surface displacement patterns, emphasizing the importance of understanding the ways in which melt can be stored in the crust. Additionally, we show that melt compressibility plays a more complex role in the core-shell model than in the uniform poroelastic model. A compressible melt in the core stores more pressure during melt resupply, which, upon release, drives post-injection inflation. In contrast, an incompressible melt in the core with a compressible interstitial melt in the shell shows post-injection deflation, as the shell absorbs residual pressure when injection ceases. Our results demonstrate how different reservoir configuration and melt distribution could impact the response of the Earth's crust to magma injections.
{"title":"Benchmarking and testing poroelastic dynamic magma mush volcano deformation models","authors":"Rami Alshembari, James Hickey, Lorenzo Mantiloni","doi":"10.1016/j.jvolgeores.2025.108468","DOIUrl":"10.1016/j.jvolgeores.2025.108468","url":null,"abstract":"<div><div>Volcano deformation studies provide crucial insights into magmatic processes, offering a key perspective on the dynamics of magma reservoirs and their eruption potential. Here, we benchmark a Finite Element (FE) poroelastic magma mush model with an analytical model of a spherical melt core and mush shell undergoing melt resupply. In our comparisons, we investigate key reservoir parameters, such as melt core pressure, deformation at the core-shell and shell-rock boundaries, and the evolution of tensile stress. The benchmarking provides excellent agreement between the FE results and the analytical solutions. We extend the analysis by incorporating a previously neglected free surface effect in the numerical models, allowing for the examination of surface deformation beyond the original analytical solution. Comparison of the surface deformation resulting from a core-shell magma mush model to that from a uniform poroelastic magma mush model highlights differences in the predicted surface displacement patterns, emphasizing the importance of understanding the ways in which melt can be stored in the crust. Additionally, we show that melt compressibility plays a more complex role in the core-shell model than in the uniform poroelastic model. A compressible melt in the core stores more pressure during melt resupply, which, upon release, drives post-injection inflation. In contrast, an incompressible melt in the core with a compressible interstitial melt in the shell shows post-injection deflation, as the shell absorbs residual pressure when injection ceases. Our results demonstrate how different reservoir configuration and melt distribution could impact the response of the Earth's crust to magma injections.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"468 ","pages":"Article 108468"},"PeriodicalIF":2.3,"publicationDate":"2025-10-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145333272","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-10DOI: 10.1016/j.jvolgeores.2025.108469
Lucas Corna , Gert Lube , Daniel H. Uhle , Ermanno Brosch , Jim R. Jones , Michael Manga , Benjamin Andrews
Pyroclastic density currents (PDCs) can cross significant topographic obstacles. The processes that govern the interaction of PDCs with obstacles remain poorly understood leaving uncertainty in hazard planning and mitigation. Here, we report the results of large-scale experimental PDCs comprising hot volcanic particles and gas propagating across ridge-shaped obstacles. Observations from high-speed video and measurements of the velocity, density and temperature structure of the flows are used to identify the flow processes that occur when PDCs propagate across and become partially blocked by hill-shaped topographic obstacles; and how these characteristics are recorded in PDC deposits. The experiments show that the interaction of PDCs with ridges generate strong local perturbations to the internal flow velocity, density and temperature structure. These flow changes are linked to three main processes: the blocking of the lower, concentrated flow region in front of the obstacle; the compression and acceleration of the non-blocked flow regions on the stoss side; and the flow detachment behind the crest and formation of a turbulent wake before flow re-attachment downstream. Flow-topography interactions result in deposition and erosion rates that vary by three and two orders of magnitude, respectively, which explain the strong asymmetry of PDC deposits across topographic obstacles. The facies architecture of experimental deposits across ridges resembles those of natural PDC deposits from Te Maari and Taupō volcanoes (New Zealand). The findings of this study can guide the interpretation of PDC deposits or be taken into consideration in numerical models simulating the propagation of PDCs across complex topography for hazard forecast.
{"title":"Transport and sedimentation of Pyroclastic Density Currents across topographic obstacles","authors":"Lucas Corna , Gert Lube , Daniel H. Uhle , Ermanno Brosch , Jim R. Jones , Michael Manga , Benjamin Andrews","doi":"10.1016/j.jvolgeores.2025.108469","DOIUrl":"10.1016/j.jvolgeores.2025.108469","url":null,"abstract":"<div><div>Pyroclastic density currents (PDCs) can cross significant topographic obstacles. The processes that govern the interaction of PDCs with obstacles remain poorly understood leaving uncertainty in hazard planning and mitigation. Here, we report the results of large-scale experimental PDCs comprising hot volcanic particles and gas propagating across ridge-shaped obstacles. Observations from high-speed video and measurements of the velocity, density and temperature structure of the flows are used to identify the flow processes that occur when PDCs propagate across and become partially blocked by hill-shaped topographic obstacles; and how these characteristics are recorded in PDC deposits. The experiments show that the interaction of PDCs with ridges generate strong local perturbations to the internal flow velocity, density and temperature structure. These flow changes are linked to three main processes: the blocking of the lower, concentrated flow region in front of the obstacle; the compression and acceleration of the non-blocked flow regions on the stoss side; and the flow detachment behind the crest and formation of a turbulent wake before flow re-attachment downstream. Flow-topography interactions result in deposition and erosion rates that vary by three and two orders of magnitude, respectively, which explain the strong asymmetry of PDC deposits across topographic obstacles. The facies architecture of experimental deposits across ridges resembles those of natural PDC deposits from Te Maari and Taupō volcanoes (New Zealand). The findings of this study can guide the interpretation of PDC deposits or be taken into consideration in numerical models simulating the propagation of PDCs across complex topography for hazard forecast.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"468 ","pages":"Article 108469"},"PeriodicalIF":2.3,"publicationDate":"2025-10-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145333273","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The area around Lake Kivu, situated in the western branch of the East African Rift System (WBEARS), which spans the Democratic Republic of Congo, Rwanda, and Burundi, is characterized in its southern part by numerous hot springs. Chemistry and isotopic signatures of hydrothermal fluids have been studied in the northern part of Lake Kivu due to the presence of one of the world's active volcanoes, Nyiragongo. However, the areas along the coast of Lake Kivu, extending south to the shores of Lake Tanganyika, where active volcanism is currently absent, have been largely overlooked, despite the presence of several tens of hydrothermal manifestations with temperatures ranging from 25 to 86 °C, located in a highly faulted region characterized by low to moderate seismicity. This study aimed to fill the gap by collecting samples from 15 hot springs, four lake waters, and one cold spring, located in the Democratic Republic of Congo and Burundi. Stable isotopes of water (δ2H and δ 18O), noble gas isotopes of He, Ne, and Ar, Sr isotopic ratio (87Sr/86Sr), and major ions and trace elements were measured. Most water samples are classified as Na-K-HCO3 or Ca-Mg-HCO3 types. Stable isotopes indicate that the waters are of meteoric origin. The helium isotopic ratios (3He/4He or R), normalized to the atmospheric ratio (Ra = 1.384 × 10−6), range from 0.058 to 1.304, indicating a dominant crustal helium source and a maximum of 13 % of mantle helium. The 87Sr/86Sr ratios range from 0.70362 in the Lac Vert sample, typical of mafic volcanic products, to 0.77520 in hot springs of the Ruzizi basin, a very radiogenic value possibly indicating water-rock interaction with the Proterozoic silicate basement of Central Africa. The 87Sr/86Sr shows a rough trend with the 40Ar/36Ar ratios, indicating a slightly detectable terrigenic 40Ar* excess in hot spring water. These results contrast with those on the northern shore of Lake Kivu, which are characterized by magmatic helium, suggesting that moving further south, the heat is controlled by the local geothermal gradient, with fluids likely circulating deeper into the crust and warming up. This circulation is facilitated by the extensional faults of the African rift, where most of the hot springs are located. Calculated low crustal fluxes of helium in the Lake Kivu region are insufficient to create economically valuable helium reserves, as found in the southern termination of the WBEARS, in the Rukwa Rift Basin.
{"title":"Sources of helium and associated heat in hydrothermal fluids from the central western branch of the East African Rift System (Democratic Republic of Congo and Burundi)","authors":"Wisdom Kambale Kavyavu , Daniele L. Pinti , Bienfait Kambale Simisi , Dario Tedesco","doi":"10.1016/j.jvolgeores.2025.108464","DOIUrl":"10.1016/j.jvolgeores.2025.108464","url":null,"abstract":"<div><div>The area around Lake Kivu, situated in the western branch of the East African Rift System (WBEARS), which spans the Democratic Republic of Congo, Rwanda, and Burundi, is characterized in its southern part by numerous hot springs. Chemistry and isotopic signatures of hydrothermal fluids have been studied in the northern part of Lake Kivu due to the presence of one of the world's active volcanoes, Nyiragongo. However, the areas along the coast of Lake Kivu, extending south to the shores of Lake Tanganyika, where active volcanism is currently absent, have been largely overlooked, despite the presence of several tens of hydrothermal manifestations with temperatures ranging from 25 to 86 °C, located in a highly faulted region characterized by low to moderate seismicity. This study aimed to fill the gap by collecting samples from 15 hot springs, four lake waters, and one cold spring, located in the Democratic Republic of Congo and Burundi. Stable isotopes of water (δ<sup>2</sup>H and δ <sup>18</sup>O), noble gas isotopes of He, Ne, and Ar, Sr isotopic ratio (<sup>87</sup>Sr/<sup>86</sup>Sr), and major ions and trace elements were measured. Most water samples are classified as Na-K-HCO<sub>3</sub> or Ca-Mg-HCO<sub>3</sub> types. Stable isotopes indicate that the waters are of meteoric origin. The helium isotopic ratios (<sup>3</sup>He/<sup>4</sup>He or R), normalized to the atmospheric ratio (Ra = 1.384 × 10<sup>−6</sup>), range from 0.058 to 1.304, indicating a dominant crustal helium source and a maximum of 13 % of mantle helium. The <sup>87</sup>Sr/<sup>86</sup>Sr ratios range from 0.70362 in the Lac Vert sample, typical of mafic volcanic products, to 0.77520 in hot springs of the Ruzizi basin, a very radiogenic value possibly indicating water-rock interaction with the Proterozoic silicate basement of Central Africa. The <sup>87</sup>Sr/<sup>86</sup>Sr shows a rough trend with the <sup>40</sup>Ar/<sup>36</sup>Ar ratios, indicating a slightly detectable terrigenic <sup>40</sup>Ar* excess in hot spring water. These results contrast with those on the northern shore of Lake Kivu, which are characterized by magmatic helium, suggesting that moving further south, the heat is controlled by the local geothermal gradient, with fluids likely circulating deeper into the crust and warming up. This circulation is facilitated by the extensional faults of the African rift, where most of the hot springs are located. Calculated low crustal fluxes of helium in the Lake Kivu region are insufficient to create economically valuable helium reserves, as found in the southern termination of the WBEARS, in the Rukwa Rift Basin.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"468 ","pages":"Article 108464"},"PeriodicalIF":2.3,"publicationDate":"2025-10-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145363368","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}