Pub Date : 2026-01-01Epub Date: 2025-10-25DOI: 10.1016/j.jvolgeores.2025.108480
K.M. Williams , A. Geyer , C. Annen , J. Kavanagh
Sill geometry and magma flow indicators preserved within fossil sills are used to determine magma source locations, understand economic potential of magmatic ore deposits, and forecast potential volcanic eruption sites. However, existing models struggle to incorporate complex flow dynamics and quantify flow variability, thus inhibiting their potential to explain spatially variable magma flow within sills spanning up to hundreds of kilometres. We present results of new 2D finite element numerical simulations coupling fluid and thermal dynamics within a dyke-fed sill using a multiphysics approach. As magma enters the model sill from below via one (or several) feeding dyke(s), magma jets of variable height develop within the sill depending on dyke thickness and inlet velocity. Low-velocity zones occur near the feeding dyke(s), with recirculation present between multiple feeding dyke(s) when they are present. These findings demonstrate the significant impact that intrusion geometry has on the magma flow dynamics within the sill, and we postulate that the presence of magma jets at the dyke-to-sill transition may be one source of sill lobes. They suggest random crystal orientations could be expected close to feeder dykes (low-velocity recirculating flow and rapid solidification), but where strain rates are high crystal alignment may still occur. These results potentially explain complex magma flow as interpreted from field observations and petrographic analysis of sill. Our results show that incorporating intrusion geometries, flow dynamics and thermal processes into models is crucial for bridging the gap between field observations and the underlying processes that govern natural systems.
{"title":"Modelling magma flow within dyke-fed sill geometries: A coupled thermal and fluid dynamics approach","authors":"K.M. Williams , A. Geyer , C. Annen , J. Kavanagh","doi":"10.1016/j.jvolgeores.2025.108480","DOIUrl":"10.1016/j.jvolgeores.2025.108480","url":null,"abstract":"<div><div>Sill geometry and magma flow indicators preserved within fossil sills are used to determine magma source locations, understand economic potential of magmatic ore deposits, and forecast potential volcanic eruption sites. However, existing models struggle to incorporate complex flow dynamics and quantify flow variability, thus inhibiting their potential to explain spatially variable magma flow within sills spanning up to hundreds of kilometres. We present results of new 2D finite element numerical simulations coupling fluid and thermal dynamics within a dyke-fed sill using a multiphysics approach. As magma enters the model sill from below via one (or several) feeding dyke(s), magma jets of variable height develop within the sill depending on dyke thickness and inlet velocity. Low-velocity zones occur near the feeding dyke(s), with recirculation present between multiple feeding dyke(s) when they are present. These findings demonstrate the significant impact that intrusion geometry has on the magma flow dynamics within the sill, and we postulate that the presence of magma jets at the dyke-to-sill transition may be one source of sill lobes. They suggest random crystal orientations could be expected close to feeder dykes (low-velocity recirculating flow and rapid solidification), but where strain rates are high crystal alignment may still occur. These results potentially explain complex magma flow as interpreted from field observations and petrographic analysis of sill. Our results show that incorporating intrusion geometries, flow dynamics and thermal processes into models is crucial for bridging the gap between field observations and the underlying processes that govern natural systems.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"469 ","pages":"Article 108480"},"PeriodicalIF":2.3,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145520729","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}
This study investigated the products of successive Vulcanian eruptions in 1976 at the Minamidake summit crater of the Sakurajima Volcano, namely pumice clasts from an explosion on May 13 and dense juvenile fragments from May 17. Effective undercooling (ΔTeff) on May 11–13 and 13–17 were estimated at 53–86 and 33–45 °C, respectively. Plagioclase phenocrysts exhibited broken surfaces with regrowth rims, which were attributed to decompression by the immediately preceding explosions. Based on the regrowth rim widths, the average growth rates were estimated to be 0.4–1.7 × 10−11 and 1.7–4.7 × 10−12 m s−1 on May 11–13 and 13–17, respectively, suggesting silicon diffusion-controlled growth. These growth rates were consistent with those from previous decompression experiments at similar ΔTeff under the same conditions (i.e., temperature and melt SiO2 and H2O content) as Sakurajima. Assuming that plagioclase nanolites and fine-grained microlites with lengths shorter than the regrowth rims widths nucleated in the corresponding periods, the nucleation rate measured on May 11–13 agreed well with the recalculated experimental rates. In contrast, the nucleation rate measured on May 13–17 was higher by at least 1.9 orders of magnitude, suggesting nucleation acceleration at a pressure of <40 MPa. The high nucleation rates increased the groundmass total crystallinity by a factor of 2–5 compared with the experiments at similar ΔTeff, enabling the groundmass viscosity to reach the stress fragmentation threshold for a Vulcanian explosion. A precise understanding of these shallow magma crystallization kinetics thus paves the way for prediction of eruptive styles.
本研究调查了1976年樱岛火山Minamidake山顶火山口连续火山喷发的产物,即5月13日爆炸的浮石碎屑和5月17日密集的幼体碎片。5月11-13日和13-17日的有效过冷(ΔTeff)估计分别为53-86°C和33-45°C。斜长石斑晶表现出破碎的表面和再生的边缘,这是由于之前的爆炸造成的减压。根据再生环宽度估算,5月11 - 13日和13-17日的平均生长率分别为0.4-1.7 × 10−11和1.7-4.7 × 10−12 m s−1,表明硅的扩散控制生长。这些生长速率与先前在类似ΔTeff条件下(即温度和熔体SiO2和H2O含量)与Sakurajima相同的减压实验一致。假设长度小于再生棱宽的斜长石纳米岩和细粒微岩在相应时期成核,则5月11-13日测量的成核速率与重新计算的实验速率吻合较好。而5月13日至17日的成核速率至少提高了1.9个数量级,表明在40 MPa的压力下成核加速。与类似ΔTeff条件下的实验相比,高成核率使地面物质的总结晶度提高了2-5倍,使地面物质粘度达到了火山爆发的应力破碎阈值。因此,对这些浅层岩浆结晶动力学的精确理解为预测喷发风格铺平了道路。
{"title":"Rates of plagioclase growth, nanolite nucleation, and viscosity increase during Vulcanian activity of Sakurajima Volcano, Japan","authors":"Aulia Syafitri , Michihiko Nakamura , Naoki Araya , Mayumi Mujin , Daisuke Miki , Masato Iguchi","doi":"10.1016/j.jvolgeores.2025.108483","DOIUrl":"10.1016/j.jvolgeores.2025.108483","url":null,"abstract":"<div><div>This study investigated the products of successive Vulcanian eruptions in 1976 at the Minamidake summit crater of the Sakurajima Volcano, namely pumice clasts from an explosion on May 13 and dense juvenile fragments from May 17. Effective undercooling (ΔT<sub>eff</sub>) on May 11–13 and 13–17 were estimated at 53–86 and 33–45 °C, respectively. Plagioclase phenocrysts exhibited broken surfaces with regrowth rims, which were attributed to decompression by the immediately preceding explosions. Based on the regrowth rim widths, the average growth rates were estimated to be 0.4–1.7 × 10<sup>−11</sup> and 1.7–4.7 × 10<sup>−12</sup> m s<sup>−1</sup> on May 11–13 and 13–17, respectively, suggesting silicon diffusion-controlled growth. These growth rates were consistent with those from previous decompression experiments at similar ΔT<sub>eff</sub> under the same conditions (i.e., temperature and melt SiO<sub>2</sub> and H<sub>2</sub>O content) as Sakurajima. Assuming that plagioclase nanolites and fine-grained microlites with lengths shorter than the regrowth rims widths nucleated in the corresponding periods, the nucleation rate measured on May 11–13 agreed well with the recalculated experimental rates. In contrast, the nucleation rate measured on May 13–17 was higher by at least 1.9 orders of magnitude, suggesting nucleation acceleration at a pressure of <40 MPa. The high nucleation rates increased the groundmass total crystallinity by a factor of 2–5 compared with the experiments at similar ΔT<sub>eff</sub>, enabling the groundmass viscosity to reach the stress fragmentation threshold for a Vulcanian explosion. A precise understanding of these shallow magma crystallization kinetics thus paves the way for prediction of eruptive styles.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"469 ","pages":"Article 108483"},"PeriodicalIF":2.3,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145520681","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 : 2026-01-01Epub Date: 2025-11-21DOI: 10.1016/j.jvolgeores.2025.108497
J.J. Rawlings , T.M. Gernon , M.J. Stock , M.R. Palmer , C.M. Petrone , R.J. Brown , E. Humphreys-Williams
Kimberlites are a class of alkaline igneous rock which sometimes contain diamonds and derive from mantle depths typically exceeding 150 km. We investigate the Upper Pleistocene to Holocene (12 ka) Igwisi Hills kimberlites to characterise the olivine grains and better understand their origins and evolution. We observe four distinct olivine populations. Macrocrysts are >1500 μm in diameter, rounded, monocrystalline grains with a Mg-rich xenocrystic core (Fo 90.5–92.8) and elevated Ni but depleted Ca and Mn contents, typical of granular mantle peridotite. Nodules are >1500 μm in diameter, rounded, polycrystalline grains with an Mg-rich core (Fo 91–92.4), also derived from granular peridotite. Microcryst grains are <1500 μm, subhedral-euhedral and can be subdivided into two types. Microcryst-a grains have a Mg-rich (Fo 90–92.5) xenocrystic core indicative of a mantle peridotite origin and microcryst-b grains have a Fe-rich (Fo 89–91) core that likely originated from disaggregated neoblasts. Olivine in all four populations displays four magmatic zones: internal zone, rim, rind and outermost rind, typical of kimberlites worldwide. Fe-rich internal zones are likely derived from a primitive kimberlite melt which entrained the cores. Mg-rich rims formed as the kimberlite melt evolved and crystallised around the internal zones. Fe-rich rinds formed as the kimberlite melt continued to ascend and further crystallised around the rims. Mg-rich outermost rinds crystallised from the kimberlite melt during the final stages of ascent. Based on these observations we present a new ascent model for the Igwisi Hills magmas.
{"title":"A classification scheme for olivine populations in the Upper Pleistocene to Holocene Igwisi Hills (Tanzania) kimberlite lavas and investigation of grain origins and evolution","authors":"J.J. Rawlings , T.M. Gernon , M.J. Stock , M.R. Palmer , C.M. Petrone , R.J. Brown , E. Humphreys-Williams","doi":"10.1016/j.jvolgeores.2025.108497","DOIUrl":"10.1016/j.jvolgeores.2025.108497","url":null,"abstract":"<div><div>Kimberlites are a class of alkaline igneous rock which sometimes contain diamonds and derive from mantle depths typically exceeding 150 km. We investigate the Upper Pleistocene to Holocene (12 ka) Igwisi Hills kimberlites to characterise the olivine grains and better understand their origins and evolution. We observe four distinct olivine populations. Macrocrysts are >1500 μm in diameter, rounded, monocrystalline grains with a Mg-rich xenocrystic core (Fo 90.5–92.8) and elevated Ni but depleted Ca and Mn contents, typical of granular mantle peridotite. Nodules are >1500 μm in diameter, rounded, polycrystalline grains with an Mg-rich core (Fo 91–92.4), also derived from granular peridotite. Microcryst grains are <1500 μm, subhedral-euhedral and can be subdivided into two types. Microcryst-a grains have a Mg-rich (Fo 90–92.5) xenocrystic core indicative of a mantle peridotite origin and microcryst-b grains have a Fe-rich (Fo 89–91) core that likely originated from disaggregated neoblasts. Olivine in all four populations displays four magmatic zones: internal zone, rim, rind and outermost rind, typical of kimberlites worldwide. Fe-rich internal zones are likely derived from a primitive kimberlite melt which entrained the cores. Mg-rich rims formed as the kimberlite melt evolved and crystallised around the internal zones. Fe-rich rinds formed as the kimberlite melt continued to ascend and further crystallised around the rims. Mg-rich outermost rinds crystallised from the kimberlite melt during the final stages of ascent. Based on these observations we present a new ascent model for the Igwisi Hills magmas.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"469 ","pages":"Article 108497"},"PeriodicalIF":2.3,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145618215","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}
Hydrothermal fluids are well known to circulate through fracture zones (e.g., faults) in geothermal areas. Impermeable clay layers can play a role of caprock, retaining geothermal reservoirs underneath. Hydrothermal alteration layers are known to include clay minerals and to have lower hydraulic permeability and higher electrical conductivity than the surrounding rocks. The magnetotelluric (MT) method of geophysical exploration can estimate the resistivity structure below the surface. Therefore, MT method can reveal relations among hydrothermal fluid circulation, impermeable clay layers, and faults. Nevertheless, because of sparseness of observation sites and the resulting low spatial resolution, earlier studies using the MT method have rarely revealed such detailed relations. For this study, high-density audio-frequency MT (AMT) soundings were conducted at 83 sites in the Unzen hot spring area, Japan, with site spacing of approximately 50–150 m to estimate details of the three-dimensional (3-D) resistivity structure near the surface (shallower than 1 km). Consequently, the distribution of conductive bodies was clearly identified in the shallow subsurface. Some conductive bodies are exposed on the surface, where hot springs and fumaroles are distributed. The conductive bodies are regarded as representing hydrothermal alteration layers. However, no shallow conductive body was found on the north side of an active fault (Oshidori-no-ike fault; OF). This feature suggests sharp variation of subsurface temperatures at OF, possibly because of recharging of meteoric water along OF. We conclude that high-density 3-D AMT surveys can reveal shallow subsurface geothermal systems.
众所周知,热液流体在地热区的断裂带(如断层)中循环。不透水的粘土层可以起到盖层的作用,保留地下的地热储层。已知热液蚀变层含有粘土矿物,并且比围岩具有较低的水力渗透率和较高的导电性。地球物理勘探中的大地电磁法可以估计地表以下的电阻率结构。因此,MT方法可以揭示热液流体循环、不渗透粘土层和断层之间的关系。然而,由于观测点的稀疏性和由此产生的低空间分辨率,使用MT方法的早期研究很少揭示这种详细的关系。本研究在日本云禅温泉地区的83个测点进行高密度音频MT (AMT)测深,测点间距约为50-150 m,以估计地表附近(浅于1 km)的三维电阻率结构细节。从而明确了浅层地下导电体的分布。一些导电体暴露在地表,地表分布着温泉和喷气孔。导电体被认为是热液蚀变层的代表。而在活动断层(oshidori -no- like fault; of)的北侧未发现浅层导电体。这一特征表明of的地下温度变化剧烈,可能是由于沿of的大气水补给所致。我们认为高密度三维AMT测量可以揭示浅层地下地热系统。
{"title":"Detailed resistivity structure by high-density AMT surveys in the geothermal area: A case study near Unzen volcanoes, Japan","authors":"Nagi Yamashita , Tada-nori Goto , Keiichi Ishizu , Kodo Umakoshi , Hiroshi Sasaki","doi":"10.1016/j.jvolgeores.2025.108466","DOIUrl":"10.1016/j.jvolgeores.2025.108466","url":null,"abstract":"<div><div>Hydrothermal fluids are well known to circulate through fracture zones (e.g., faults) in geothermal areas. Impermeable clay layers can play a role of caprock, retaining geothermal reservoirs underneath. Hydrothermal alteration layers are known to include clay minerals and to have lower hydraulic permeability and higher electrical conductivity than the surrounding rocks. The magnetotelluric (MT) method of geophysical exploration can estimate the resistivity structure below the surface. Therefore, MT method can reveal relations among hydrothermal fluid circulation, impermeable clay layers, and faults. Nevertheless, because of sparseness of observation sites and the resulting low spatial resolution, earlier studies using the MT method have rarely revealed such detailed relations. For this study, high-density audio-frequency MT (AMT) soundings were conducted at 83 sites in the Unzen hot spring area, Japan, with site spacing of approximately 50–150 m to estimate details of the three-dimensional (3-D) resistivity structure near the surface (shallower than 1 km). Consequently, the distribution of conductive bodies was clearly identified in the shallow subsurface. Some conductive bodies are exposed on the surface, where hot springs and fumaroles are distributed. The conductive bodies are regarded as representing hydrothermal alteration layers. However, no shallow conductive body was found on the north side of an active fault (Oshidori-no-ike fault; OF). This feature suggests sharp variation of subsurface temperatures at OF, possibly because of recharging of meteoric water along OF. We conclude that high-density 3-D AMT surveys can reveal shallow subsurface geothermal systems.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"469 ","pages":"Article 108466"},"PeriodicalIF":2.3,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145520731","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 : 2026-01-01Epub Date: 2025-11-25DOI: 10.1016/j.jvolgeores.2025.108502
Phil Shane , Shanice Mascarenhas , Teresa Ubide , Shane J. Cronin
Understanding the explosivity of future volcanic eruptions is a cornerstone of hazard evaluation, and magma storage days to years leading up to eruption could provide timely information for mitigation strategies. To explore the influence of pre-ascent conditions on eruption style we compared the growth histories of plagioclase and amphibole phenocrysts from four prominent, Holocene-aged, explosive erupted deposits with six effusive eruptions at Taranaki volcano, New Zealand. The explosively erupted andesite-dacite magmas display a range of bulk compositions and crystallinities that overlap with those erupted effusively. Regardless of their eruption style, each phenocryst population displays a high degree of inter-crystal textural, geochemical and Sr-isotope discordance, including in their final rim growth. This can be explained by the periodic intrusive disruption of a redeveloping crystal mush that lacked a unidirectional assimilation-fractional-crystallisation path in the centuries to millennia between eruptions. However, amphiboles from most of the effusive eruptions are distinguished from those erupted explosively by recording late-stage halogen (and likely other volatiles) accumulation in the magma and subsequent decomposition during slow ascent. Also more prevalent in the effusive eruptions are plagioclase textures that record late-stage mingling with a mafic melt. These features support conceptual models that suggest magma compressibility due to volatile accumulation allows more efficient magma mingling during recharge and thus heating. The resulting melt viscosity reduction and enhanced outgassing of the magma then promotes effusive eruption. Testing this hypothesis via mineralogical determination of intensive parameters would need to overcome the high degree of disequilibrium that occurred at all stages of the magma evolution.
{"title":"Does pre-ascent magma storage influence eruption style? An evaluation from textures and geochemistry of phenocrysts at Taranaki volcano","authors":"Phil Shane , Shanice Mascarenhas , Teresa Ubide , Shane J. Cronin","doi":"10.1016/j.jvolgeores.2025.108502","DOIUrl":"10.1016/j.jvolgeores.2025.108502","url":null,"abstract":"<div><div>Understanding the explosivity of future volcanic eruptions is a cornerstone of hazard evaluation, and magma storage days to years leading up to eruption could provide timely information for mitigation strategies. To explore the influence of pre-ascent conditions on eruption style we compared the growth histories of plagioclase and amphibole phenocrysts from four prominent, Holocene-aged, explosive erupted deposits with six effusive eruptions at Taranaki volcano, New Zealand. The explosively erupted andesite-dacite magmas display a range of bulk compositions and crystallinities that overlap with those erupted effusively. Regardless of their eruption style, each phenocryst population displays a high degree of inter-crystal textural, geochemical and Sr-isotope discordance, including in their final rim growth. This can be explained by the periodic intrusive disruption of a redeveloping crystal mush that lacked a unidirectional assimilation-fractional-crystallisation path in the centuries to millennia between eruptions. However, amphiboles from most of the effusive eruptions are distinguished from those erupted explosively by recording late-stage halogen (and likely other volatiles) accumulation in the magma and subsequent decomposition during slow ascent. Also more prevalent in the effusive eruptions are plagioclase textures that record late-stage mingling with a mafic melt. These features support conceptual models that suggest magma compressibility due to volatile accumulation allows more efficient magma mingling during recharge and thus heating. The resulting melt viscosity reduction and enhanced outgassing of the magma then promotes effusive eruption. Testing this hypothesis via mineralogical determination of intensive parameters would need to overcome the high degree of disequilibrium that occurred at all stages of the magma evolution.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"469 ","pages":"Article 108502"},"PeriodicalIF":2.3,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145618217","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-12-01Epub Date: 2025-10-21DOI: 10.1016/j.jvolgeores.2025.108478
Prasanta K. Patro , S. Dhamodharan , Venkata Durga , K.K. Abdul Azeez , Narendra Babu , K. Chinna Reddy , Arvind K. Gupta , M. Shiva Krishna
The Panamik–Changlung hot springs, located along the Karakoram Fault (KF) in the Trans-Himalayan region of Ladakh in India, were investigated for the first time using the Magnetotelluric (MT) method to assess their geoelectrical properties. Situated in the Shyok–Nubra Valley between the Ladakh and Karakoram Batholiths, the region shows abundant geothermal manifestations. We acquired MT data along a 30 km profile parallel to the KF and performed 3D joint inversion of the impedance tensor and tipper data using the ModEM program. The resulting resistivity model reveals an upper crustal conductor (<10 Ω.m) at 4–10 km depth acting as geothermal reservoir with strong out-of-quadrant phases (>90° at periods >10 s) indicate a more complex system, potentially involving current channeling and electrical anisotropy. A shallower conductor (surface to ∼1.5 km) corresponds to sediments of Nubra Formation and fractured zones of the Karakoram Metamorphic Complex saturated with mixed thermal fluids and groundwater.
{"title":"Delineation of a geothermal source beneath the Panamik-Changlung Hot Springs along the Karakoram Fault, Ladakh, India, using magnetotelluric studies","authors":"Prasanta K. Patro , S. Dhamodharan , Venkata Durga , K.K. Abdul Azeez , Narendra Babu , K. Chinna Reddy , Arvind K. Gupta , M. Shiva Krishna","doi":"10.1016/j.jvolgeores.2025.108478","DOIUrl":"10.1016/j.jvolgeores.2025.108478","url":null,"abstract":"<div><div>The Panamik–Changlung hot springs, located along the Karakoram Fault (KF) in the Trans-Himalayan region of Ladakh in India, were investigated for the first time using the Magnetotelluric (MT) method to assess their geoelectrical properties. Situated in the Shyok–Nubra Valley between the Ladakh and Karakoram Batholiths, the region shows abundant geothermal manifestations. We acquired MT data along a 30 km profile parallel to the KF and performed 3D joint inversion of the impedance tensor and tipper data using the ModEM program. The resulting resistivity model reveals an upper crustal conductor (<10 Ω.m) at 4–10 km depth acting as geothermal reservoir with strong out-of-quadrant phases (>90° at periods >10 s) indicate a more complex system, potentially involving current channeling and electrical anisotropy. A shallower conductor (surface to ∼1.5 km) corresponds to sediments of Nubra Formation and fractured zones of the Karakoram Metamorphic Complex saturated with mixed thermal fluids and groundwater.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"468 ","pages":"Article 108478"},"PeriodicalIF":2.3,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145424715","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-12-01Epub Date: 2025-09-15DOI: 10.1016/j.jvolgeores.2025.108444
Susumu Umino
The Fukube Falls in the upper Jadani River exposes successive layers of essentially subaerial ignimbrite of the 70 Ma Futomiyama Group, that intervenes a rheoignimbrite layer of giant pillows. The pillowed layer extends 320 m from east to west and varies in thickness from 15 m to 84 m. The bifurcating pillow lobes indicate that the western, central and eastern pillows flowed southward, eastward and southwestward, respectively. The pillow lobes range from 1 to more than 9 m in height and from 1 to 13 m in width, with an average of 3.4 ± 1.6 m and 4.2 ± 2.1 m, respectively. The height/width ratio is 0.78, which is similar to the ranges for tholeiite pillows of the Recent Loihi Seamount, Hawaii, and larger than those of the Cretaceous Oman Ophiolite. The obsidian glass of the pillow margin has a rhyolitic composition with 77 mass% SiO2. The magma temperature is estimated from 782 °C to 810 °C, depending on the water content, based on the hornblende geothermometer. Analogue experiments on pillow lobe height constrain the upper and lower limits of extrusion rate and viscosity of the rheoignimbrite. In addition, the Devorah number, the ratio of rate of strain accumulation to the rate of strain relaxation, also gives the maximum viscosity and extrusion rate that allows for the formation of coherent pillows without fragmentation. Based on these constraints, the conditions of pillow lobe formation are viscosities ranging from 104.5 Pa·s to 106.8 Pa·s, and extrusion rates from 2.4 × 10−3 m3/s to 3.1 × 103 m3/s for point sources and from 2.4 × 10−3 m2/s to 3.1 × 102 m2/s for line sources. Such viscosities require water contents of the magma more than 1.4 mass%, suggesting the hydrous welding of ignimbrite that was initially extruded on land and was emplaced under water-saturated conditions.
{"title":"Pillow flows of rhyolite rheoignimbrite exposed in the Fukube Falls, Hakusan-Shirakawago White Road","authors":"Susumu Umino","doi":"10.1016/j.jvolgeores.2025.108444","DOIUrl":"10.1016/j.jvolgeores.2025.108444","url":null,"abstract":"<div><div>The Fukube Falls in the upper Jadani River exposes successive layers of essentially subaerial ignimbrite of the 70 Ma Futomiyama Group, that intervenes a rheoignimbrite layer of giant pillows. The pillowed layer extends 320 m from east to west and varies in thickness from 15 m to 84 m. The bifurcating pillow lobes indicate that the western, central and eastern pillows flowed southward, eastward and southwestward, respectively. The pillow lobes range from 1 to more than 9 m in height and from 1 to 13 m in width, with an average of 3.4 ± 1.6 m and 4.2 ± 2.1 m, respectively. The height/width ratio is 0.78, which is similar to the ranges for tholeiite pillows of the Recent Loihi Seamount, Hawaii, and larger than those of the Cretaceous Oman Ophiolite. The obsidian glass of the pillow margin has a rhyolitic composition with 77 mass% SiO<sub>2</sub>. The magma temperature is estimated from 782 °C to 810 °C, depending on the water content, based on the hornblende geothermometer. Analogue experiments on pillow lobe height constrain the upper and lower limits of extrusion rate and viscosity of the rheoignimbrite. In addition, the Devorah number, the ratio of rate of strain accumulation to the rate of strain relaxation, also gives the maximum viscosity and extrusion rate that allows for the formation of coherent pillows without fragmentation. Based on these constraints, the conditions of pillow lobe formation are viscosities ranging from 10<sup>4.5</sup> Pa·s to 10<sup>6.8</sup> Pa·s, and extrusion rates from 2.4 × 10<sup>−3</sup> m<sup>3</sup>/s to 3.1 × 10<sup>3</sup> m<sup>3</sup>/s for point sources and from 2.4 × 10<sup>−3</sup> m<sup>2</sup>/s to 3.1 × 10<sup>2</sup> m<sup>2</sup>/s for line sources. Such viscosities require water contents of the magma more than 1.4 mass%, suggesting the hydrous welding of ignimbrite that was initially extruded on land and was emplaced under water-saturated conditions.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"468 ","pages":"Article 108444"},"PeriodicalIF":2.3,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145107581","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-12-01Epub Date: 2025-10-22DOI: 10.1016/j.jvolgeores.2025.108475
Gaetano Ferrante , Helge Gonnermann , Céline Fliedner , Thomas Giachetti , Amy G. Ryan
{"title":"Corrigendum to “Viscosity of bubbly magmas from torsional experiments on pumice” [Journal of Volcanology and Geothermal Research 461 (2025) 108297]","authors":"Gaetano Ferrante , Helge Gonnermann , Céline Fliedner , Thomas Giachetti , Amy G. Ryan","doi":"10.1016/j.jvolgeores.2025.108475","DOIUrl":"10.1016/j.jvolgeores.2025.108475","url":null,"abstract":"","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"468 ","pages":"Article 108475"},"PeriodicalIF":2.3,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145528405","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-12-01Epub 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-12-01","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}
Pub Date : 2025-12-01Epub 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-12-01","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}