Magma Eruption Research Articles

Diffusion of Sr and Ba in plagioclase: Composition and silica activity dependencies, and application to volcanic rocks

Strontium and barium diffusion chronometry in plagioclase is routinely applied to mafic and felsic magmatic systems. This technique can be used to determine the timescales of magma reservoir assembly and the cooling rates of plutons and volcanic rocks, which has emerged as a useful method to assess volcanic hazards. Here we report diffusion experiments that aim to constrain the diffusivities of Sr and Ba in oligoclase and labradorite at 1 atm pressure, between 900 and 1,200°C, and assessing diffusion in different crystallographic orientations. In all of the reported experiments, silica activity (aSiO2) is buffered by varying stable phase assemblages in the diffusant source powder. The experimental products were analysed by secondary ion mass spectrometry (SIMS) depth profiling and laser ablation inductively coupled mass spectrometry (LA-ICP-MS) line scanning. There is no resolvable dependence of Sr and Ba diffusion in plagioclase upon aSiO2 or crystal orientation. However, Sr and Ba diffusivities are found to vary as functions of the plagioclase anorthite content. As such, we parameterise the diffusivity of Sr and Ba in plagioclase as a function of temperature and anorthite content as follows:log10DSr(m2s−1)=−1.65(±0.24)XAn−3.03(±1.16)−[368,142(±27,141)2.303RT],log10DBa(m2s−1)=−1.43(±0.20)XAn−4.65(±0.96)−[337,037(±22,969)2.303RT],where XAn is the plagioclase anorthite content (mole fraction), T is the temperature (K), and R is the gas constant (J mol-1 K-1). The diffusion rate of Sr in plagioclase determined in this study is ∼1.5–2 orders of magnitude slower than previously determined, whereas Ba diffusion is similar to previous studies. This is most likely due to the Ba-feldspar stability at the experimental conditions employed by previous studies, whereas Sr-feldspar was absent from the source powder assemblage. By applying the diffusivities determined in this study to plagioclase crystals from the Cerro Galán ignimbrite (Argentina) and Santorini caldera (Greece), we find timescales of ∼105 years, with a good agreement between results from Sr and Ba diffusion modelling. Therefore, our data reconcile experimental diffusion data with measured Sr and Ba profiles in plagioclase and suggest that, at least regarding the Cerro Galán ignimbrite and Santorini caldera, plagioclase records the time needed to differentiate magma reservoirs and assemble large volumes of eruptible magma.

Volcanic-intrusive connections and crystal-melt segregation in the Dulan tilted crustal section: Insights from accessory mineral evolution

The volcanic-intrusive connection is critical for understanding the formation of transcrustal magma systems and the mechanisms driving magma eruption. Yet, the genetic links between abyssal intrusive rocks and acidic volcanic rocks remain poorly understood. Here, we present new data from the Dulan volcanic-intrusive complex, a tilted section of the East Kunlun orogenic belt, integrating zircon UPb dating, trace elements, apatite and zircon Hf isotope geochemistry, and whole-rock geochemistry to investigate volcanic-intrusive connections and track crystal-melt segregation processes. The Dulan complex comprises rhyolite, granite porphyry, and granodiorite, all showing consistent zircon UPb ages, rare earth element (REE) patterns, and εHf(t) values, suggesting a shared magma source. Glomerocrysts of plagioclase, alkali feldspar, and hornblende, along with antecrystic zircon and apatite, imply a cumulate affinity in granodiorite and granite porphyry. Compositional variations in whole-rock and accessory minerals indicate that plagioclase, apatite, and zircon fractionation significantly influenced magma evolution. We propose a two-stage crystal segregation process in a ∼ 2 Myr-long magma reservoir system. In the early stage, at ∼33 km depth, aqueous fluid influx initiated crystal-melt segregation, forming the Dulan granodiorite from residual mush. In the later stage, mafic melt and fluid injections caused crystal-melt segregation in a shallower reservoir at ∼12 km depth, with high-silica magma producing rhyolite, while residual mush consolidated into granite porphyry. Our findings suggest that volcanic-intrusive connections are more complex than mere compositional equivalence or complementarity, requiring multidimensional analysis through accessory minerals and whole-rock geochemistry.

Unique seismic and eruption precursors to the 1996 and ongoing magmatic eruptions of Popocatépetl: Coupled and fluidized bed events

We describe three unique types of seismicity at Popocatépetl volcano that accompanied the initial vent-clearing eruptive activity in December 1994 through the eruption of the first two domes in 1996. We identify and describe two types of coupled events, 1) spasmodic burst coupled events, a burst of volcano tectonic (VT) events coupled with a large eruptive explosion (21 December 1994 and 5 March 1996), and 2) explosion-couplet coupled events, a pair of events with a deeper first event and whose second event correlates with a small gas emission or explosion. Explosion-couplets occurred with the onset of magmatic ash dominated eruptions and their properties are very useful for forecasting magmatic eruptions at Popocatépetl volcano. Measurable quantities including the time between the first and second phase, daily numbers of events, and the amplitude ratio between the second and first phase systematically changed as the first two domes approached the surface between March and June 1996. Interevent times decreased from many tens of seconds to a few seconds, while amplitude ratios increased from about 2 to 10 or more. Event numbers increased prior to and during initial dome extrusion. These changes were used to forecast the eruption of the first two domes. We also report on another unusual type of low frequency seismicity that accompanied emissions called fluidized bed events. Fluidized bed events occurred following the initial eruption on 21 December 1994, beginning mid-January 1995, when SO2 emissions were initially elevated, and the conduit was open. Fluidized bed events began to wane in late March through the end of May 1995 as did SO2 emissions. Consistent gas, visual and seismic observations by Centro Nacional de Prevención de Desastres (CENAPRED) enabled the direct correlation of the gas, magma, and seismic phenomena.

Hydraulically linked reservoirs simultaneously fed the 1975–1984 Krafla Fires eruptions: Insights from petrochemistry

The 1975–1984 Krafla Fires in northeast Iceland was the first plate-boundary rifting episode to be tracked using seismic and geodetic monitoring. Geophysical observations from this episode have inspired conceptual models of magma transport during plate spreading, but a lack of complementary petrologic insights has hindered a holistic understanding of the events. To address this knowledge gap, we studied the petrochemistry of all nine Krafla Fires basaltic eruptions. Our large dataset of new whole-rock, matrix glass and mineral analyses from samples collected during or shortly after each eruption reveal a clear compositional bimodality in the erupted magmas that persisted across the episode, with evolved quartz tholeiite (MgO = 5.7–6.4 wt.%) erupted inside Krafla caldera, and more primitive (usually olivine-normative) tholeiite (MgO = 6.4–8.7 wt%) erupted north of the caldera margin. Barometric calculations indicate tapping of these magmas from distinct reservoirs: a primitive lower-crustal reservoir at a most probable depth of ∼14–19 km, and a more evolved, shallower reservoir at a most probable depth of ∼7–9 km beneath the caldera. These reservoirs were tapped simultaneously in several of the nine eruptions, and in three events the two magma types mixed near the northern caldera margin. Varying levels of trace element depletion in the deep-sourced primitive melts reflect incomplete mixing of diverse mantle-derived melts at depth; the most enriched of these melts could be parental to evolved inside-caldera magma via fractional crystallization. Clinopyroxene rims on gabbroic nodules from primitive September 1984 lavas record lower crustal pressures, while diffusion models suggest that these rims grew up to within a few months before eruption. Ascent of the primitive magma from the lower crust thus occurred over timescales much shorter than eruptive repose periods, without prolonged stalling at shallow depths. These observations are inconsistent with the view that the eruptions were entirely fed by lateral magma outflow from the shallow reservoir. They instead require some decoupling of the flow paths of the two magma types: the primitive magma either bypassed the sub-caldera reservoir laterally or ascended vertically beneath the northern vents. The two reservoirs nonetheless shared a hydraulic connection and jointly responded to rifting. Comparison of the Krafla Fires with other rifting events and eruptions highlights the complexity and diversity of magma transport during plate boundary rifting events, which is not yet captured by a generalizable model. Integration of petrologic, geochemical and geophysical data is essential to provide a holistic view of future rifting events in Iceland and at other spreading centres.

The 2019 pumice raft forming eruption of Volcano-F (Volcano 0403–091) and implications for hazards posed by submerged calderas

Low volcanic explosivity index (VEI) eruptions are common occurrences in the Southwest Pacific but, as demonstrated by the 2021/2022 eruption of Hunga Volcano, submerged calderas in the region are also capable of producing much larger and more hazardous eruptions. As such, characterising smaller events from potentially hazardous systems is essential. The 2019 eruption of Volcano-F, a submerged caldera, would likely have gone totally undetected had it not produced a pumice raft that inundated beaches in Fiji and eventually washed up in Australia. New data, acquired 5 months after the eruption, reveal the development of a new vent and the accumulation of at least 3.1*107 m3 bulk volume (dense rock equivalent of 5.6*106 m3) of material on the seafloor. Between 30 and 70% of erupted material entered the raft, while the rest remained near to or was dispersed down-current of the vent. This previously unaccounted for material increases the volume estimate for the eruption, confirming it as a VEI 3 event and highlights the importance of considering not just the floating component of a pumice raft forming eruption for VEI estimation. Geochemical analysis reveals the eruption comprised a homogenous batch of dacitic magma, with compositional characteristics similar to that erupted from the same volcano in 2001, and an until-now-unidentified pumice raft in the Coral Sea in 1964. Volcano-F therefore appears to have had at least three explosive eruptions in the last 60 years, indicating it is significantly at unrest. Repeated eruptions of similar composition and low crystal content magma over decadal to centennial scales indicate the existence of a melt-dominant magma body beneath the volcano. Submerged calderas, like Volcano-F, are common in the wider Southwest Pacific region, with many such calderas producing regular eruptions, implicating active magmatic recharge. Our findings motivate a need to carefully monitor and characterise even apparently small eruptions at this volcano, and others along the Tonga-Kermadec Arc. This is because such eruptions have the potential to subsequently prime or trigger more explosive eruptions and provide critical geochemical evidence about the plumbing system and evolution of the volcano, essential for understanding the diverse hazards they pose.

The color systematics of volcanic ashfall samples in estimating eruption sequences: a case study of the 2017–2018 eruption at Shinmoe-dake, Kirishima volcano, Southwest Japan

The color of pyroclasts is fundamental, because it reflects various magma properties and eruption processes, including particle morphology, chemistry, and petrological characteristics. However, deriving the componentry ratio (CR) of pyroclasts for ongoing eruption monitoring remains challenging due to the lack of a robust quantitative standard. The derivation of the CR, as well as other petrological analyses, is too laborious and time-consuming to introduce as a sustainable monitoring method. To address this, we employed spectroscopic colorimetry to rapidly and quantitatively describe eruptive product colors, enabling CR derivation based on clear, objective standards for ash particle classification. Through color spectroscopy of bulk and sieved ash samples, we analyzed the major size fraction for time-series samples during the waxing stage of the 2017–2018 Shinmoe-dake eruption in Kirishima volcano, Southwest Japan. Our findings reveal that the color changes in bulk ash systematically changed with the evolution of componentry. This temporal color change was due to an increase in the amount of vesicular particles with clear glass against dark angular lava particles, as well as a grain size change, which we interpret as an indication of a transition from phreatic/phreatomagmatic to magmatic eruption. Subsequently, the color of the ash changed when the amount of different lava particles increased gradually, coinciding with a shift toward a more dominant effusion of lava. As the lava effusion continued, a slight reddening of the ash, indiscernible to the naked eye, was clearly detected by the spectrometer before the onset of intermittent Vulcanian eruptions. We interpreted this reddening as oxidation resulting from decreased ascent speed and caprock formation, which accumulates overpressure for Vulcanian explosions. These results highlight the potential of rapid, objective color value and componentry derivation for sustainable quasi-real-time monitoring of ongoing eruption materials.Graphical

Volcanic Eruption in the Nanoworld: Efficient Oxygen Exchange at the Si/SnO2 Interface.

Here, a phenomenon of efficient oxygen exchange between a silicon surface and a thin layer of tin dioxide during chemical vapor deposition is presented, which leads to a unique Sn:SiO2 layer. Under thermodynamic conditions in the temperature range of 725-735°C, the formation of nanostructures with volcano-like shapes in "active" and "dormant" states are observed. Extensive characterization techniques, such as electron microscopy, X-ray diffraction, synchrotron radiation-based X-ray photoelectron, and X-ray absorption near-edge structure spectroscopy, are applied to study the formation. The mechanism is related to the oxygen retraction between tin(IV) oxide and silicon surface, leading to the thermodynamically unstable tin(II)oxide, which is immediately disproportionate to metallic Sn and SnO2 localized in the SiO2 matrix. The diffusion of metallic tin in the amorphous silicon oxide matrix leads to larger agglomerates of nanoparticles, which is similar to the formation of a magma chamber during the natural volcanic processes followed by magma eruption, which here is associated with the formation of depressions on the surface filled with metallic tin particles. This new effect contributes a new approach to the formation of functional composites but also inspires the development of unique Sn:SiO2 nanostructures for diverse application scenarios, such as thermal energy storage.

Petrogenesis of flood basalts and shield basanite from Mertolemariam- Abamineos area, northwestern Ethiopian plateau: Assessments for mantle source variations

Petrographic and geochemical data (major and trace elements) are presented for Oligocene flood basalts and Miocene shield lavas from the Mertolemariam-Abamineos area, northwestern Ethiopian Plateau to examine the petrogenesis of the erupted magmas and the nature of mantle source compositions. The Mertolemariam flood basalts have mainly aphyric and plagioclase phyric textures. The Abamineos shield lavas have sparsely clinopyroxene phyric to highly plagioclase-clinopyroxene phyric textures. Geochemical classification shows that the Mertolemariam Oligocene flood basalts are sub-alkaline in composition, whereas the Abamineos Miocene shield lavas are highly alkaline in composition. The Abamineos shield has a unique composition (i.e., basanites) compared with all other shield basalts (transitional to alkaline) in the northwestern Ethiopian Plateau. Major and trace element compositions display two distinct trends between Mertolemariam flood basalts and Abamineos shield basanites. Fractional crystallization, partial melting, and crustal contamination of a homogeneous mantle source cannot explain the compositional variations between the Mertolemariam flood basalts and the Abamineos shield basanites. The trace element composition suggests that the Mertolemariam Oligocene flood basalts were more likely generated from the mixing of OIB (mantle plume) and E-MORB (enriched asthenosphere) mantle components. The Abamineos Miocene shield basanites were derived from the mantle plume (OIB) component. LREE/MREE and LREE/HREE indicate that both groups possibly originated from a mantle source within the stability field of spinel and garnet. In comparison, the Mertolemariam flood basalts were formed by a higher degree of partial melting from relatively shallow depths than the Abamineos shield basanites. We propose a scenario that explains the magmatic genesis in the northwestern Ethiopian Plateau: volcanism initiated by the initial arrival of the mantle plume (OIB-like) beneath the lithosphere, which comes across a geochemically fertile and enriched MORB (E-MORB) mantle component in the upper asthenosphere. The hot mantle plume triggered melting of the fertile and enriched MORB, and then melting occurred in the plume (OIB-like) at depth in the stability field of spinel and garnet. Melts from the mantle plume (OIB-like) and E-MORB components mixed to produce sub-alkaline flood basalts during the Oligocene in the northwestern Ethiopian Plateau Subsequently, melts from the advanced upwelling mantle plume (OIB-like) produced Miocene shield lavas in the northwestern Ethiopian Plateau.

Resolving the “megacryst paradox”: Feldspar orientation relationships record crystal mobility in granites

K-feldspar megacrysts are common in silicic plutons, but there is a long-running debate around how they form and what their presence tells us about magmatic systems. Field, textural, and geochemical evidence supports growth in a melt-rich environment, but experimental evidence and phase-equilibria modeling indicate that K-feldspar grows late in the crystallization sequence, when the magma is highly crystalline. We provide a new perspective on this problem by examining the arrangement of plagioclase inclusions within megacrysts to test whether they exhibit the systematic low-energy crystallographic relationships expected from attachment by synneusis in melt-rich environments where crystals have space to rotate. We use electron backscatter diffraction to quantify the crystal orientations and find that the megacrysts’ plagioclase inclusions do occupy these preferred orientations and therefore were incorporated in a melt-rich environment. K-feldspar is also present as an interstitial network, but plagioclase crystals hosted within this network have non-systematic orientations. This transition from systematic to non-systematic plagioclase orientations marks the point at which the crystals formed a rigid, interconnected framework that impeded rotation into low-energy orientations. Phase-equilibria modeling indicates that this transition occurred when the magma was ∼55% crystalline. The remaining ∼45% melt crystallized at the eutectic, forming the interstitial phases. Thus, we resolve the “megacryst paradox”; the megacrysts grew freely in melt, and the groundmass K-feldspar formed after crystal lock-up. Megacrysts therefore provide a detailed textural and chemical record of a critical period in the system’s evolution: the transition from a mobile and potentially eruptible magma to an immobile mush.

Newly identified small vulcanian eruptions during the caldera-forming stage of Towada Volcano, Northeast Japan

Understanding the detailed eruptive history and conditions during the buildup to catastrophic caldera-forming eruptions is essential for understanding the evolution of caldera volcanoes and predicting eruptive hazards. Towada Volcano is an active caldera volcano in the northern part of the Northeast Japan Arc. At least two catastrophic caldera-forming eruptions (eruptive episodes N and L) occurred during its caldera-forming stage (61–15.7 ka). A detailed geological survey identified three small vulcanian tephra layers that were erupted during the caldera-forming stage. The tephra layers are blue–gray ash fall deposits that consist mainly of fresh blocky dacite–rhyolite fragments. In ascending order, these deposits are assigned to eruptive episodes O′, N′, and M′. Each eruption had a volume of <0.11 km3, which is much smaller than that of other eruptions during the caldera-forming stage. The eruptive ages of episodes O′, N′, and M′ are estimated to be ca. 40, 23.1, and 17.2 ka, respectively, based on new 14C ages and paleosol thicknesses data. At least three cycles of a medium- to large-scale (3–20 km3) explosive eruption, preceded by a small vulcanian eruption (<0.11 km3), occurred during the latter half of the caldera-forming stage. The small vulcanian eruptions (episodes O′, N′, and M′) occurred after relatively long periods of quiescence (4000–14,000 years) and were followed by medium- to large-scale explosive eruptions (episodes N, M, and L) 1500–4000 years later. This cyclic activity probably reflects changes in overpressure in the magma reservoir. As the overpressure increased, episodes O′, N′, and M′ probably occurred when the overpressure had not increased sufficiently to trigger a medium- to large-scale explosive eruption, and episodes N, M, and L occurred at higher overpressures. These cycles characterize the fermentation phase of Towada Volcano, and terminated with the formation of the Towada Caldera during episode L at 15.7 ka. The magma composition and eruption frequency changed abruptly after episode L, and a small stratovolcano was formed through intermittent eruptions of mafic magma. This change suggests that the caldera collapse during episode L changed the entire shallow magmatic system that existed during the fermentation phase, and the system shifted to a recovery (post-caldera) phase. Based on eruptive volumes and frequencies, the present Towada Volcano has not yet reached the conditions that existed during the late caldera-forming stage and is therefore unlikely to produce a catastrophic caldera-forming eruption in the near future.

Submarine volcanism in the Sicilian Channel revisited

The origin and role of volcanism in continental rifts remains poorly understood in comparison to other volcano-tectonic settings. The Sicilian Channel (central Mediterranean Sea) is largely floored by continental crust and represents an area affected by pronounced crustal extension and strike-slip tectonism. It hosts a variety of volcanic landforms closely associated with faults, which can be used to better understand the nature and distribution of rift-related volcanism. A paucity of appropriate seafloor data in the Sicilian Channel has led to uncertainties regarding the location, volume, sources and timing of submarine volcanism. To improve on this situation, we use newly acquired geophysical data (multibeam echosounder and magnetic data, sub-bottom profiles) and dredged seafloor samples to: (i) re-assess the evidence for submarine volcanism in the Sicilian Channel and define its spatial pattern, (ii) infer the relative age and style of magmatism, and (iii) relate this to the dominant tectonic structures in the region. Quaternary rift-related volcanism has been focused at Pantelleria and Linosa, at the northwest boundaries of their respective NW-SE trending grabens. Subsidiary and older volcanic sites potentially occur at the Linosa III and Pantelleria SE seamounts, collectively representing the only sites of recent volcanism that can be directly related to the main rift process. These long-lived polygenetic volcanic landforms have been shaped by magmatism that is directly correlated with extensional faulting and buried igneous bodies. Older volcanic landforms, sharing a similar scale and alignment, occur to the north at Nameless Bank and Adventure Bank. These deeply eroded volcanoes have likely been inactive since the Pliocene and are probably related to earlier stages of crustal thinning and underlying feeder structures in the northern region of the Sicilian Channel. Along a similar alignment, Pinne Bank, SE Pinne Bank and Cimotoe in the northern Sicilian Channel lack a surface volcanic signature but are associated with intrusive bodies or deeply buried volcanic rock masses. Terrible Bank, in the same region, also shows evidence of ancient, polygenetic magmatism, but was subject to significant erosion and lacks a prominent alignment. The much younger volcanism at Graham Volcanic Field and along the northern Capo-Granitola-Sciacca Fault Zone differs markedly from that observed in the other study areas. Here, the low-volume and scattered volcanic activity is driven by shallow-water mafic magma eruptions, which gave rise to small individual cones. These sites are associated with large fault structures away from the main rift axis and may have a distinct magmatic origin. Dispersed active fluid venting occurs across both ancient and young volcanic sites in the region and is directly associated with shallow magmatic bodies within tectonically-controlled basins. Our study provides the foundation for an updated tectonic and magmatic framework for the Sicilian Channel, and for future detailed chronological and geochemical assessment of the sources and evolution of magmatic processes in the region.

Improved artificial bee colony algorithm for pressure source parameter inversion of Sakurajima volcano from InSAR data

A novel artificial bee colony algorithm was introduced for the eruption event of the Sakurajima volcano on August 9, 2020, to invert the magma source characteristics below the volcano based on the point source Mogi model. Considering that the Sakurajima volcano is surrounded by sea, all the deformation data are used to obtain the location and magma eruption volume of the volcano. In response to the weak local search capability of the artificial swarm algorithm, the difference between the global optimal individual and the un-roulette screened individual is introduced as the variance component in the onlooker stage. Detailed simulation experiments verify the improvement of the algorithm in terms of convergence speed. In real experiments, the Sakurajima volcano inversion shows closer fitting results and smaller residuals compared to the existing literature. Meanwhile, the convergence speed of the algorithm echoes with the simulation experiments.

Changes in the magmatic plumbing system associated with the Gotemba sector collapse at Mount Fuji, Japan

It is well known that magmatic plumbing systems change over time, but there is much debate as to why and how. We studied volcanic ejecta continuously deposited in an outcrop at Kagosaka Pass at the eastern base of Mount Fuji to investigate the factors responsible for changes in the magmatic plumbing system. The sample consisted of pyroclastic sediments from explosive eruptions for approximately 3000 y, which sandwiched the time of the Gotemba sector collapse at approximately 2500 BP. Chemical analyses of whole rocks, minerals, and matrix glasses, as well as mode measurements of glass and bubbles, were performed on samples collected from approximately 30 layers; significant changes were observed before and after the collapse. For example, before and after the collapse, matrix glass area increased around 60% to over 80% and anorthite content (Ca / (Ca + Na) * 100) of phenocryst plagioclase decreased from over 80 to below 65. For a period after the collapse, possibly hundreds of years, the plagioclase and olivine phenocrysts exhibited characteristics indicative of crystallization at low temperatures and pressures, and the pyroclast matrix became highly vitreous. Eruptions with ejecta of these characteristics continued more than a dozen times, lasting about 500 years. In addition, the trend in the distribution of the bulk rock chemical composition changed significantly, showing a differentiation trend with only plagioclase and clinopyroxene crystal separation. An investigation using the MELTS software revealed that the phenomenon of direct eruptions from deep magma chambers to the surface, bypassing shallow magma reservoirs, continued for several hundred years after the collapse. This can be interpreted as a decrease in confining pressure associated with the collapse, facilitating the eruption of magma from the depths. Furthermore, based on an examination of the water content in the magma during this period, we posit that the trigger for the rise of magma from the deep magma chamber of Mount Fuji is the acquisition of excess pressure by the injection of magma from a deeper level.

Petrological and geochemical insights into the magma plumbing system of the Daliuchong dacite eruption, Tengchong Volcanic Field, SW China

The formation of highly evolved, dacitic magmas has been attributed to various processes, including crystal fractionation, partial melting of overlying crust, and/or assimilation of crustal material into an evolving magma chamber. These processes are undoubtedly primary processes involved in the formation of dacites, but they may not be the only mechanism involved in the formation of high-silica dacites. For instance, mafic magma replenishment has been proposed as an additional mechanism but has not been assessed, and thus, its role has not been well-constrained. The Daliuchong volcano is the result of one of the largest eruptive events within the Tengchong Volcanic Field (TVF) in southwest China during the Early-Middle Pleistocene. Here, we conducted detailed mineral textures, mineral chemistry, and geochemical studies on Daliuchong pyroclastic rocks to explore the pre-eruptive storage conditions and evolution processes of the magma. The Daliuchong pyroclastic rocks are dacitic in composition. The samples show porphyritic textures characterized by phenocrysts of plagioclase, amphibole, clinopyroxene, orthopyroxene, and Fe-Ti oxides. Additionally, two distinct types of glomerocryst are identified: a gabbroic glomerocryst containing plagioclase + clinopyroxene + orthopyroxene ± Fe-Ti oxides assemblage and a dioritic glomerocryst containing plagioclase + amphibole ± pyroxene ± Fe-Ti oxides assemblage. Both phenocryst and glomerocryst show rich micro-textures. The Daliuchong dacite exhibits bulk compositional heterogeneity. Analysis of bulk-rock data suggests that this heterogeneity may arise from both the differentiation of the dacite itself and the injection of mafic magma. The compositional similarity between the Daliuchong dacite and experimentally produced partial melts of metamorphic basalt supports that the Daliuchong dacite was predominantly formed through the partial melting of the mafic lower crust. Thermobarometry estimation indicates that clinopyroxenes with high Mg# crystallized at 560–870 MPa, whereas amphibole and clinopyroxenes with low Mg# crystallized at 185–300 MPa. Based on the observed phase relations and the calculated crystallization conditions, we propose that during the differentiation of the Daliuchong dacite, heterogeneous dacitic magma formed by partial melting accumulated in a deep magma reservoir (21–32 km) before subsequently ascending toward shallower depths. Crystallization of plagioclase, amphibole, Fe-Ti oxides, and small amounts of pyroxene and apatite occurred at a shallower depth (7–10 km). The presence of coarse-sieve texture, fine-sieve texture, and oscillatory zoning with high amplitude in plagioclase suggests intermittent injection of mafic magma into the shallow magma reservoir, with the eruption of dacitic magma occurring after the final mafic magma replenishment. The petrological evidence above advocates that primitive magma replenishment could have been involved in the formation and triggered the eruption of dacite in the Daliuchong volcano.

Effect of crustal stress state on magmatic stalling and ascent: case study from Puyehue-Cordón Caulle, Chile

The Southern Volcanic Zone (SVZ) in Chile is an active continental arc with a complex history of volcanism, where a range of magmatic compositions have been erupted in a variety of styles. In the Central SVZ, both monogenetic and polygenetic volcanoes exist, in close proximity to the Liquiñe-Ofqui Fault System (LOFS), but with variable local stress states. Previous studies have inferred varying crustal storage timescales, controlled by the orientation of volcanic centres relative to the N-S striking LOFS and σHMax in this region. To assess the relationship between volcanism and crustal stress states affected by large-scale tectonic structures and edifice controls, we present whole rock geochemical data, to ensure consistency in source dynamics and crustal processing, mineral-specific compositional data, thermobarometry, and Fe–Mg diffusion modelling in olivine crystals from mafic lavas, to assess ascent timescales, from the stratovolcanic edifice of Puyehue-Cordón Caulle and proximal small eruptive centres. Textural observations highlight differences in crystal maturation timescales between centres in inferred compression, transpression, and extension, yet source melting dynamics remain constant. Only samples from the stratovolcanic edifice (in regional compression) preserve extensive zonation in olivine macrocrysts; these textures are generally absent from proximal small eruptive centres in transtension or extension. The zonation in olivines from stratovolcanic lavas yields timescales on the order of a few days to a few weeks, suggesting that even in environments which inhibit ascent, timescales between unrest and eruption of mafic magmas may be short. Significantly, high-resolution compositional profiles from olivine grains in the studied samples record evidence for post-eruptive growth and diffusion, highlighting the importance of careful interpretation of diffusion timescales from zoned minerals in more slowly cooled lavas when compared with tephra samples.

The Nieve volcanic cluster: A Pliocene - Pleistocene lava dome cluster in the Michoacán-Guanajuato volcanic field (México)

The Nieve monogenetic volcanic cluster is located in the central–eastern region of the Michoacán–Guanajuato volcanic field, along the Huiramba fault zone, a relay ramp in the Morelia–Acambay fault system produced by oblique north-northwest transtension. This volcanic cluster includes at least 17 middle Pliocene to late Pleistocene lava domes, two small shield volcanoes, and two scoria cones. Between 4 and 3.8 Ma, two effusive eruptions built two small shield volcanoes overlying one another, with a magma volume of 3.93 km3. Between 2.9 Ma and 21.4 ka, 17 lava domes and two scoria cones were emplaced on the flanks of these volcanoes. The entire cluster resulted in a total erupted volume of 17 km3, covering an area of 326 km2 and reaching a thickness of emplaced volcanic material of 1200 m, resulting in a magma eruption rate equivalent to 0.004 km3/ka. All the rocks associated with this cluster are within a relatively restricted range in composition, between 53.9 and 64.2 wt% SiO₂, from andesite enriched in silica to basaltic andesite. The presence of intrusive-rock xenoliths and xenocrysts with dissolution textures reveals that assimilation processes modified the magmas. Based on the regional geological record, we suggest that the establishment of the Nieve volcanic cluster has been controlled by tectonic structures and the basement of the region, which has allowed the chemical evolution of these magma batches that probably had sources in at least two deep reservoirs as reflected by the Nb/Th versus Ta/U ratio.

Early arc crust formation preserved in the Grenadines archipelago, southern Lesser Antilles arc.

Intra-arc diversity in volcanic activity and composition is ubiquitous, but its underlying causes remain largely unresolved in many settings. In this work, we examine such variability in the Grenadines archipelago, southern Lesser Antilles arc. Here, juxtaposed volcanic centres exhibit eruptive longevities and chemistries distinct from northern counterparts in the same arc. Our goal is to explain this deviation by investigating variations in magmatic processes using petrological data from erupted crustal xenoliths and lavas, and interpreting these findings within the context of the archipelago's tectonic history and geophysical structure. Textural analyses of xenoliths reveal crystallization over a wide range of pressure-temperature-melt composition conditions in the crust. Mineral phases display discrete compositional trends pointing towards significant inter-island variability in underlying plumbing systems. The geochemical variety of erupted magmas is reminiscent of the entire arc. We speculate that the Grenadines represents the early onset of subduction forming the modern-day Lesser Antilles arc. Extrusive volcanism initiated as submarine activity. Subsequent uplift eroded the original topography of these volcanic centres following the eventual cessation of volcanism in the Neogene. The positioning of the Grenadines on an elevated platform provides rare modern insight into early arc crust formation not commonly preserved in established active arcs.

A multi-methodological approach to record dynamics and timescales of the plumbing system of Zaro (Ischia Island, Italy)

Determining the time spans of processes related to the assembly of eruptible magma at active volcanoes is fundamental to understand magma chamber dynamics and assess volcanic hazard. This information can be recorded in the chemical zoning of crystals. Nevertheless, this kind of study is still poorly employed for the active volcanoes of the Neapolitan area (Southern Italy), in particular, for Ischia island where the risk is extremely high and this information can provide the basis for probabilistic volcanic hazard assessment. For these reasons, we acquired chemical composition on clinopyroxene crystals erupted at Ischia during the Zaro eruption (6.6 ± 2.2 ka) and performed numerical simulations of the input of mafic magma into a trachytic reservoir, in order to investigate various aspects of pre-eruptive dynamics occurring at different timescales. This event emplaced a ~ 0.1 km3 lava complex, in which the main trachytic lava flows host abundant mafic to felsic enclaves. Previous petrological investigation suggested that mafic magma(s) mixed/mingled with a trachytic one, before the eruption. In this work, the clinopyroxene zoning patterns depict the growth of crystals in different magmatic environments, recording sequential changes occurred in the plumbing system before the eruption. The evolution of the plumbing system involved a hierarchy of timescales: a few hours for magma mingling caused by mafic recharge(s) and likely occurred multiple times over a decade during which a dominant magmatic environment was sustained before the eruption. Such timescales must be considered in volcanic hazard assessment at Ischia and similar active volcanoes in densely populated areas.

Mesozoic–Cenozoic tectonic–palaeogeographical evolution of Bayingobi Basin: response to subduction and collision of the Mongol–Okhotsk Ocean plate

The structures, sedimentary fill and magmatic rocks of sedimentary basins are the products of regional tectonic evolution and are crucial to studying the basin's tectonic–palaeogeographical environment and regional tectonic evolution. The present study discusses the zircon U–Pb dating of volcanic rock samples, in situ Hf isotope and whole-rock geochemistry analysis, and the palaeogeographical, sedimentary filling and tectonic background of the Bayingobi basin through outcrop and core descriptions of Meso-Cenozoic strata, stratigraphic correlation and lithofacies analysis. The U–Pb age of the volcanic rocks is 132–102 Ma, and the in situ Hf isotopic values range from −20.99 to +29.48. The 87 Sr/ 86 Sr and 143 Nd/ 144 Nd isotopic ratios of the volcanic rocks are 0.707049–0.879761 and 0.511846–0.512540, respectively, and the Nd isotopic values range from −0.62 to −6.83. The volcanic rocks in the basin have the characteristics of continental margin island arc volcanic rocks (CAA) with loss of Nb and Ta and enrichment of Pb, originating mainly from dehydration melting of the subducting plate. The lower Cretaceous Bayingobi, Suhongtu and Yingen Formations developed specifically delta–lacustrine deposits. The lower and upper member of the Bayingobi Formation was deposited in the Berriasian to Valanginian and Valanginian to late Aptian, respectively. The Suhongtu Formation was deposited during the late Aptian to early Albian and was controlled by the strike of the Engeer Us fault. The Yingen Formation was deposited in the late faulted depression stage in the late Albian. Because of the subduction of the Mongol–Okhotsk Ocean plate in the Early Cretaceous, several depressions (sags) were formed in the Bayingobi basin, accompanied by the eruption of continental plate margin island arc magmas. With the closure of the Mongol–Okhotsk Ocean in the late Early Cretaceous, the basin was uplifted as a whole and the Upper Cretaceous Ulansuhai Formation was deposited. Supplementary material : Supplementary figures and tables are available at https://doi.org/10.6084/m9.figshare.c.7105807

Magma intrusion process during pre-magmatic period (2010−2013) of Sinabung volcano as revealed by seismicity of volcano-tectonic and hybrid earthquakes

Sinabung volcano in Indonesia reactivated, with a series of phreatic eruptions, in August and September 2010. These eruptions were followed by various types of magmatic eruption at the summit, including lava dome growth, pyroclastic density currents, lava flow, and vulcanian eruptions, starting in December 2013. Prior to the magmatic eruptions, Sinabung exhibited two sequences of phreatic eruption and an increase in the seismicity of volcanic earthquakes. This study clarifies the progress of magma intrusion during the pre-magmatic period based on the hypocenter distribution, waveform similarity, seismic moment, rupture length, and stress drop for volcano-tectonic (VT) and hybrid earthquakes. It was found that the hypocenters of VT earthquakes are distributed north and northwest of the summit at a depth range of 1–10 km below sea level. Deep seismicity (4–10 km) alternated with shallow seismicity (2–4 km). The last shallow seismicity, from July to mid-December 2013, was different from previous seismicity in terms of higher intensity, as determined from the seismic moment, a temporary increase in the rupture length, and the migration of hypocenters towards the summit. The seismicity of VT earthquakes was altered by a swarm of hybrid earthquakes in mid-December. The hybrid earthquakes were smaller, in terms of seismic moment (mostly <3 × 1010 Nm), than the VT earthquakes. Their hypocenters were concentrated in the shallowest depth range (−0.5 to 1.5 km below sea level) directly below the summit, their source process was repeatable (they could be grouped into six earthquake families), their dominant peak was in a low-frequency range (2.5–4.5 Hz), and they had a relatively low stress drop (<0.15 MPa). This suggests that the swarm of hybrid earthquakes was induced by a frequent repeated fracture of fluid-filled cracks due to the intrusion of magma up to directly below the summit. The transition of the earthquake family and the change in its source parameters, namely an increase in the stress drop prior to the appearance of a lava dome and then a slight decrease, may reflect a gradual change in internal pressure in the hypocentral zone through the magma intrusion process.