Journal of Volcanology and Geothermal Research最新文献

{"title":"Comprehensive orientation method using cross-line laser: A case study of paleomagnetic analysis on pyroclastic deposits from Asama Volcano (Japan)","authors":"Tatsuo Kanamaru ,&nbsp;Kuniyuki Furukawa ,&nbsp;Akari Maki","doi":"10.1016/j.jvolgeores.2025.108293","DOIUrl":"10.1016/j.jvolgeores.2025.108293","url":null,"abstract":"<div><div>Accurate orientation of terrestrial materials is essential in earth sciences, particularly in structural geology and paleomagnetic studies. In this study, we present a novel and versatile sample orientation method utilizing commercially available cross-line lasers. This method address challenges in collecting oriented geological samples caused by factors such as sample size, fragility, rough surfaces, and local magnetic anomalies that complicate the use of magnetic compasses. To evaluate the effectiveness of this method, paleomagnetic measurements were performed on juvenile clasts from pyroclastic density current deposits and the pyroclastic fall deposit of the 1783 eruption of Asama Volcano. The result coincides within error limits with the geomagnetic field predicted by a paleosecular variation curve, demonstrating the potential of our orienting method. However, further experiments are required to quantify the precision of our orienting method and to assess its dependence on factors such as clast size and shape. Preliminary paleomagnetic analyses of the pyroclastic fall deposit suggest a complex emplacement history and indicate a surprisingly high-temperature emplacement. These findings offer new insights into the emplacement processes and thermal history of pyroclastic deposits. This innovative method enhances the capability to collect and analyze oriented geological samples, providing significant contributions to volcanological and structural geological research.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"461 ","pages":"Article 108293"},"PeriodicalIF":2.4,"publicationDate":"2025-02-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143421907","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}

{"title":"The spatiotemporal development of downsag and trapdoor structures during caldera subsidence with 1–10 km in diameter in analogue sandbox experiments","authors":"Ryuhei Sanjo ,&nbsp;Toshihiko Sugai","doi":"10.1016/j.jvolgeores.2025.108294","DOIUrl":"10.1016/j.jvolgeores.2025.108294","url":null,"abstract":"<div><div>Calderas are volcanic depressions formed when the roof of a magma chamber collapses due to the depletion of magma within the chamber. Investigating the development of caldera boundary faults that accommodate chamber-roof subsidence is crucial for understanding caldera collapse events and predicting contemporaneous and subsequent volcanic eruptions. Because subsurface structures in natural calderas are difficult to observe, analogue sandbox and numerical experiments, and inversions of geodetic data are often used to reconstruct caldera structures. Recent experimental and analytical studies have revealed localized downsag during the upward propagation of caldera boundary faults and that the fault structures of trapdoor calderas are highly variable. However, how downsag becomes localized during upward boundary fault propagation, and when and why trapdoor structures originate during caldera collapse remain open questions. Here, we performed analogue sandbox experiments in a transparent sandbox and used image analysis techniques to provide insights into these questions. Calderas with 1–10 km in diameter are more favorable than larger ones (tens of kilometers in diameter) for assessing downsagging and trapdoor subsidence, because these volcanic deformations may be masked by crustal deformation due to regional tectonics controlling large caldera subsidence. We therefore focused on calderas with 1–10 km in diameter, and excluded some factors such as pre-existing stress regimes, faults and magma flow, in order to evaluate the development of downsag and trapdoor subsidence themselves. Despite these experimental limitations, our experiments are in agreement with natural calderas and show that as caldera subsidence increases, downsag and trapdoor faulting evolve in four stages. (1) At the beginning of downsag, the horizontal displacement velocity is greatest along the edge of the downsagging region. (2) As the outward-dipping reverse faults nucleate at depth and propagate upward, the peak of horizontal displacement velocity localizes along the surface projection of the concealed faults. (3) Caldera collapse blocks then undergo trapdoor subsidence when the outward-dipping reverse faults reach the surface. (4) Finally, a second episode of trapdoor subsidence of the collapse block occurs when inward-dipping normal faults nucleate beyond the outward-dipping reverse faults. The first and second stages suggest that the existence and location of concealed caldera boundary faults can be predicted from the downsag deformation pattern. The third and fourth stages explain temporal and structural variations observed in trapdoor calderas in the context of caldera boundary fault development. Our model will improve the resolution of caldera structural reconstructions and associated inferences regarding magma chamber dynamics.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"461 ","pages":"Article 108294"},"PeriodicalIF":2.4,"publicationDate":"2025-02-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143421908","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Microearthquakes identification based on convolutional neural networks and clustering techniques","authors":"Fernando Lara ,&nbsp;Román Lara-Cueva ,&nbsp;Felipe Grijalva ,&nbsp;Ana Zambrano","doi":"10.1016/j.jvolgeores.2025.108282","DOIUrl":"10.1016/j.jvolgeores.2025.108282","url":null,"abstract":"<div><div>Microearthquakes are critical for understanding volcanic activity, leading to monitoring many volcanoes worldwide with seismic sensor networks. These networks generate a substantial amount of data, making visual analysis challenging. Consequently, researchers have focused on developing automatic microearthquake recognition systems over the past decades. A primary challenge with these systems is their reliance on labeled databases for training supervised learning models, where the output labels depend on the database labels. We propose using clustering algorithms in conjunction with a Fine-tuned Convolutional Neural Network (CNN) as a feature extractor to identify overlapping microearthquakes, and other types of microearthquakes withoutneeding labeled datasets. This methodology has two stages: The First stage relies on Transfer Learning, to specialize the CNN in microearthquake recognition. The Second stage uses the Fine-tuned CNN as a feature extractor. This methodology is applied to the Cotopaxi Volcano and validated in the Llaima Volcano. It uses unsupervised databases to find clusters of isolated events with similar characteristics to Long Period (LP), Volcano Tectonic (VT), Tremor (TRE), among others. Additionally, it identifies a cluster with overlapping microearthquakes. In the validation stage, 79 % of the VT events are associated to the same cluster without the need to adjust the Fine-tuned CNN. This test is performed on a dataset of a volcano never seen by CNN or Clustering algorithms. Normalized Entropy is used as a metric to verify the generalization of knowledge, the proposed work is compared with Uniform Manifold Approximation and Projection for Dimension Reduction (UMAP). The proposed work obtains 0.04 lower uncertainty with respect to UMAP.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"460 ","pages":"Article 108282"},"PeriodicalIF":2.4,"publicationDate":"2025-02-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143378816","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}

{"title":"Development of permeable networks by viscous-brittle deformation in a shallow rhyolite intrusion. Part 2: Microstructural analysis","authors":"Taylor Witcher ,&nbsp;Steffi Burchardt ,&nbsp;Tobias Mattsson ,&nbsp;Michael J. Heap ,&nbsp;Anne Pluymakers ,&nbsp;Kai Li ,&nbsp;Peter Lazor","doi":"10.1016/j.jvolgeores.2025.108278","DOIUrl":"10.1016/j.jvolgeores.2025.108278","url":null,"abstract":"<div><div>Volcanic and magmatic outgassing mechanisms can determine eruptive behavior of shallow silicic magma bodies. Most outgassing mechanisms proposed take place along conduit margins, where the highest strain rates drive ascending magma to brittle failure. However, these mechanisms do not account for outgassing large volumes of magma away from the conduit walls. Here, we present a continuum of porosity preserved in the microcrystalline rhyolitic Sandfell laccolith, Eastern Iceland. Three stages in the continuum are described: porous flow bands, pore channels, and fracture bands. These deformation features are present throughout the entire exposed volume of the Sandfell laccolith in meter-long band geometries, ranging from mm- to dm-scale thickness, and interlayered with coherent, undeformed rhyolite. Using microstructural analytical methods and drawing on the result of previous experimental studies, we show that emplacement-related deformation induced strain partitioning around a crystal content of 45 % that resulted in the segregation of melt-rich and melt-poorer flow bands. Subsequent deformation induced by continued magma emplacement caused strain partitioning in the melt-rich flow bands. Depending on strain rate, different types of deformation features developed, through dilation or porosity redistribution (porous flow bands), cavitation (pore channels), or tensile fracture (fracture bands). Porous flow bands have permeability values ∼4 orders of magnitude higher than undeformed rhyolite. Pore channels and fracture bands have much larger length scales, and so permeability increases dramatically in those systems. Hence, the abundance and interconnectivity of deformation features preserved in the Sandfell laccolith provided an efficient outgassing mechanism for the bulk of the intrusion. Outgassing due to viscous-brittle magma deformation during magma emplacement should therefore be considered for crystal-rich magmas, e.g., during effusive lava dome extrusion.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"461 ","pages":"Article 108278"},"PeriodicalIF":2.4,"publicationDate":"2025-02-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143421909","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Multiple dikes make eruptions easy","authors":"Agust Gudmundsson","doi":"10.1016/j.jvolgeores.2025.108284","DOIUrl":"10.1016/j.jvolgeores.2025.108284","url":null,"abstract":"<div><div>Dikes supply magma to most volcanic eruptions. Understanding how propagating dikes may, or may not, reach the surface is thus one of the fundamental tasks for volcanology. Many, perhaps most, dike segments injected from magma sources do not reach the surface to feed volcanic eruptions. Instead, the dike segments become arrested (stop their propagation), commonly at or close to contacts between mechanically dissimilar layers/units, at various crustal depths. This means that many and perhaps most volcanic unrest periods with dike injections do not result in eruptions. There are several conditions that make dike arrest likely, but the main one is layering where the layers have contrasting mechanical properties. Such layering means that local stresses are heterogeneous and anisotropic and, therefore, in some layers unfavourable for dike propagation – hence the dike arrest. Here I show that once a dike has formed, however, its very existence tends to make the local stress field along the dike homogeneous (with invariable orientation of principal stresses) and favourable (with dike-parallel orientation of the maximum compressive principal stress) for later dike injections. This means that subsequence dikes may use an earlier dike as a path, either along the margin or the centre of the earlier dike, thereby generating a multiple dike. Because earlier feeder-dikes form potential paths for later-injected dikes to the surface, many volcanic eruptions are fed by multiple dikes. Examples include recent eruptions in the volcanoes Etna (Italy) and Kilauea (Hawaii), and the Icelandic volcanoes Krafla, Hekla, Fagradalsfjall, and the Sundhnukur crater row. Thus, multiple dikes favour dike propagation to the surface; thereby making dike-fed eruptions easier.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"460 ","pages":"Article 108284"},"PeriodicalIF":2.4,"publicationDate":"2025-02-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143388290","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"The influence of graben geometry on dike propagation","authors":"Khaled Rahimi ,&nbsp;Marcus Bursik ,&nbsp;Janine L. Kavanagh","doi":"10.1016/j.jvolgeores.2024.108254","DOIUrl":"10.1016/j.jvolgeores.2024.108254","url":null,"abstract":"<div><div>Understanding dike behavior under varying topographical stresses is essential for predicting volcanic vent locations. This study addresses how graben geometry influences dike propagation by employing an experimental approach to simulate dike propagation, using fluids as analogs for magmas. We inject these fluids into the center of the base of an elastic gelatin medium, molded with an elongated topographic low representing a graben structure in the Earth's crust. Experiments involving different graben geometries were conducted to assess how local stress fields influence dike propagation. The results demonstrate that shallow grabens facilitate vertical to sub-vertical dike propagation within the graben, while deeper grabens cause dike deflection, often leading to off-graben venting or the formation of canoe-like sheet intrusions. These findings introduce new variables into the understanding of volcanic vent locations and eruption forecasts, particularly in rift zones with complex topography. The results offer critical implications for improving volcanic hazard assessments by incorporating the influence of topographical stress on dike propagation. Overall, the study advances our understanding of the mechanisms controlling dike propagation and provides a basis for further research into mitigating risks associated with volcanic eruptions in a graben setting.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"458 ","pages":"Article 108254"},"PeriodicalIF":2.4,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143180896","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}

{"title":"Magmatic differentiation and plumbing system beneath Nyamulagira volcano (Virunga Volcanic Province, East African Rift)","authors":"Ephrem Kamate Kaleghetso ,&nbsp;Olivier Namur ,&nbsp;Benoît Smets ,&nbsp;Jacqueline Vander Auwera ,&nbsp;François Lubala ,&nbsp;Thomas Van Gerve ,&nbsp;Sander M. Molendijk","doi":"10.1016/j.jvolgeores.2024.108264","DOIUrl":"10.1016/j.jvolgeores.2024.108264","url":null,"abstract":"<div><div>Nyamulagira volcano (eastern Democratic Republic of the Congo) is one of Africa's most active volcanoes and poses a significant threat to the nearby rapidly expanding population centers. Situated in the Virunga Volcanic Province in the western branch of the East African Rift System, Nyamulagira's frequent eruptions offer a valuable opportunity to study effusive rift volcanism. However, understanding its subsurface processes remains challenging, impeding monitoring efforts.</div><div>Here, we report a comprehensive mineralogical, petrological, and geochemical study of rocks from historical eruptions at Nyamulagira. Mineral textures and compositions show strong spatial and temporal variations. We identify multiple active magma reservoirs in which magma compositions and textures evolve via crystal accumulation, fractional crystallization, magma recharge, mixing, and convection.</div><div>We distinguish three distinct reservoirs: a deep reservoir at ∼ 22–30 km depth, the main storage region at ∼ 13–18 km depth, and a shallow reservoir at ∼ 2–9 km depth. We propose that differentiation at Nyamulagira integrates olivine and clinopyroxene crystallization and accumulation in the deep reservoir, and fractional crystallization and magma mixing/homogenization in the interconnected intermediate and shallow reservoirs. Primitive magmas from the deep reservoir are predominantly emitted via distal eruptions, whereas more evolved magmas from the shallower reservoirs have been erupted on the flanks or at the summit caldera in recent decades.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"458 ","pages":"Article 108264"},"PeriodicalIF":2.4,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143180871","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}

{"title":"Corrigendum to “Evolution of a large Quaternary monogenetic field; the multifaceted volcanism of the Serdán-Oriental basin, México” [Journal Of Volcanology and Geothermal Research, 446 (2024) 107999].","authors":"Jaime Alberto Cavazos-Alvarez ,&nbsp;Gerardo Carrasco-Núñez ,&nbsp;Giovanni Sosa-Ceballos ,&nbsp;Federico Lucci","doi":"10.1016/j.jvolgeores.2024.108239","DOIUrl":"10.1016/j.jvolgeores.2024.108239","url":null,"abstract":"","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"458 ","pages":"Article 108239"},"PeriodicalIF":2.4,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143180350","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}

{"title":"Modeling magma recharge dynamics during the 2016 Nevados de Chillan eruption: An interacting two-chamber system evidenced by petrology and geodesy","authors":"Camila Novoa Lizama ,&nbsp;D. Remy ,&nbsp;J.C. Baez ,&nbsp;A. Oyarzun ,&nbsp;S. Bonvalot ,&nbsp;A. Hooper","doi":"10.1016/j.jvolgeores.2024.108253","DOIUrl":"10.1016/j.jvolgeores.2024.108253","url":null,"abstract":"<div><div>The Nevados de Chillan is a large composite stratovolcanic complex located in the Southern Andean Volcanic Zone, at the south of Chile. Its last eruption spanned six years from January 8, 2016 until January 2023. After three years without any deformation, exhibiting phreatic and phreatomagmatic activity, an uplift episode was reported in June 2019, marking the start of its magmatic phase. How this geodetic activity evolved in the following three years is still not clear, neither how it correlates with its superficial activity. Here we analyze InSAR time series, as well as daily GNSS time series from 2015 to early 2022, highlighting the wide range of ground surface displacements observed during the eruption. In-depth analysis of these displacements together with recent petrological and geochemical results leads us to consider a recharge mechanism involving a double-reservoir model to explain the observed geodetic activity. We develop an analytical model of dynamic magma flow, which we coupled with a boundary element method to account for any geometry of the reservoirs and the topography of the volcano. The model consists of a shallow elongated source located at 5.8 km depth below the volcanic edifice connected by an incompressible magma-filled hydraulic pipe to a deeper sill like source centered at 15 km depth. We propose that the activation of the system started with a small magma intrusion of one month of duration, which was sufficient to overheat the hydrothermal system and re-mobilize magma in the shallow chamber, thus explaining the non-deformation during its phreatic phase and the slight subsidence observed during its phreato-magmatic stage. Then, we propose that a new and larger magma intrusion occurred in June 2019, explaining the uplift episode observed, which continued for the following three years of the eruption, decaying exponentially. Our model indicates that this intrusion was triggered by magma coming from the crust-mantle boundary to the deep reservoir at constant rate of 0.016 km³y⁻¹ from June 2019 to Jan. 2022, with small changes to this rate that would explain the small fluctuations observed during this uplift episode. We show that the conduit-dominated magma transfer between both reservoirs controls the dynamics of the system. A deep mafic reservoir recharging an evolved shallow reservoir would explain the mafic enclaves found in the dacites in the latter eruption and thus offers a physical model to jointly explain the observations we get from petrology, geochemistry and geophysics, bridging the disciplines. The model presented here can be used to study potential recharge mechanism occurring in the different stages of an eruptive cycle including pre-eruptive, co-eruptive and post-eruptive stages.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"458 ","pages":"Article 108253"},"PeriodicalIF":2.4,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143180868","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Development of the Late Oligocene Zákupy Diatreme, Czech Republic: New insights into protracted diatreme emplacement","authors":"Vladislav Rapprich ,&nbsp;Michael S. Petronis ,&nbsp;Sarah F. Trevino ,&nbsp;Jennifer Lindline ,&nbsp;Benjamin van Wyk de Vries ,&nbsp;Ondřej Pour ,&nbsp;Zuzana Tasáryová ,&nbsp;Matt Heizler","doi":"10.1016/j.jvolgeores.2024.108257","DOIUrl":"10.1016/j.jvolgeores.2024.108257","url":null,"abstract":"<div><div>The Zákupy diatreme is well exposed by an abandoned quarry on Kamenický Vrch Hill (Czech Republic) located in the Lusatian (Lausitz) Volcanic Field in the eastern part of the NE-SW trending Ohře Rift. The abandoned quarry's walls extend upward over 50 m exposing the lower and upper diatreme facies. Petrographic analyses revealed two main compositionally distinct suites of intrusive and extrusive bodies: dikes, sills, and lava-lakes with relics of fine-grained basanite (to limburgite) and a set of rhönite-bearing picrobasalt dikes and sills. The parallel, almost linear, and relatively fractionated (La_N/Yb_N 21.98–26.21) chondrite-normalized REE patterns suggest the erupted magmas originated through low-degree partial melting and did not undergo any significant fractionation. Close relationships between basanites and picrobasalts is confirmed also with mineral chemistry with almost identical compositional trends of clinopyroxene and spinel-group phases. The negligible difference between bulk-rock Mg# (0.71–0.72 for basanites and 0.74 for picrobasalts) and melt Mg# calculated from olivine composition (0.70 for basanites and 0.68–0.70 for picrobasalts) reflects rather a load of mantle-derived (high Mg#) xenoliths and xenocrysts rather than the presence of cummulitic olivine. Anisotropy of magnetic susceptibility data reveals magnetic fabrics of variable orientations that reflect upward, downward, and subhorizontal magma flow within the diatreme complex. Paleomagnetic data fall into three general groupings that we interpret to reflect magma emplacement during at least three phases. These data also support tilting, rotation, and slumping of diatreme host rocks causing physical displacement of the intrusions following remanence acquisition. Four new ⁴⁰Ar/³⁹Ar age determinations from the picrite (<em>n</em> = 1) and basinite (<em>n</em> = 3) yield relatively undisturbed age spectra that constrain the age of emplacement between 31.24 and 32.40 Ma, indicating that magmatic evolution of the Zákupy diatreme spanned 1.2 million years. This multidisciplinary study reveals that the Zákupy diatreme evolved from several similarly sourced magmatic pulses. Repeated molten fuel coolant interactions produced a complex diatreme setting consisting of multiple eruptive events leading to coherent intrusions, various pyroclastic and hyaloclastite materials, local contact metamorphism, and other evidence of explosive hot rock-fluid interactions. The complex evolution and eruptive behavior, together with a protracted time of the diatreme activity, illuminate that small, seemingly simple volcanoes experience a very dynamic evolution involving multiple intrusive events, changing eruptive styles, and surface morphologies.</div></div>","PeriodicalId":54753,"journal":{"name":"Journal of Volcanology and Geothermal Research","volume":"458 ","pages":"Article 108257"},"PeriodicalIF":2.4,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143180869","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}