Magmatic Conditions Research Articles

Mineralogy and petrology of the mafic Panjal Traps, Kashmir, India: constraints on magmatic conditions and post-emplacement metamorphism

Mineralogy and petrology of the mafic Panjal Traps, Kashmir, India: constraints on magmatic conditions and post-emplacement metamorphism

Direct observation of degassing during decompression of basaltic magma.

Transitions in eruptive style during volcanic eruptions strongly depend on how easily gas and magma decouple during ascent. Stronger gas-melt coupling favors highly explosive eruptions, whereas weaker coupling promotes lava fountaining and lava flows. The mechanisms producing these transitions are still poorly understood because of a lack of direct observations of bubble dynamics under natural magmatic conditions. Here, we combine x-ray radiography with a novel high-pressure/high-temperature apparatus to observe and quantify in real-time bubble growth and coalescence in basaltic magmas from 100 megapascals to surface. For low-viscosity magmas, bubbles coalesce and recover a spherical shape within 3 seconds, implying that, for lava fountaining activity, gas and melt remain coupled during the ascent up to the last hundred meters of the conduit. For higher-viscosity magmas, recovery times become longer, promoting connected bubble pathways. This apparatus opens frontiers in unraveling magmatic/volcanic processes, leading to improved hazard assessment and risk mitigation.

Extreme geochemical fractionation during mantle melting: Insights from Hf-Nd isotopically ultra-depleted eclogite

Extreme fractionation of elements and isotopes in mafic igneous rocks is abnormal in deciphering the source nature and melting conditions of mafic magmatism. Especially, identification of geochemically ultra-depleted mafic melts and their mantle sources has great bearing on the property of crust-mantle differentiation at plate margins. This is illustrated by extreme Hf-Nd isotope fractionation in ultrahigh-pressure eclogites from the Sulu orogen in east-central China. In addition to the previous finding of ultrahigh εNd(0) values, we report here new data of whole-rock trace elements and Lu-Hf isotopes in eclogites and related rocks from the same region. The present results show extremely high Lu/Hf ratios and abnormally high εHf(0) values of up to 576 for the eclogites, significantly different from the garnet amphibolites and other eclogite-facies metamorphic rocks in the same orogen. This feature is coupled with the ultrahigh εNd(0) values as well as the severe depletion of light rare earth elements (LREE) and high field strength elements (HFSE). Because HFSE and LREE are immobile in aqueous solutions and the effect of melt extraction is insignificant during the continental deep subduction, the extreme fractionation of Lu/Hf and Sm/Nd indicate their origination from a geochemically ultra-depleted mantle source. These eclogites have the depleted mantle Hf model ages of 1.27 Ga to 1.61 Ga, similar to the depleted mantle Nd model ages of 1.39 Ga to 1.67 Ga as previously reported. This suggests that the protolith of the extremely high εHf-εNd eclogites was a kind of mafic igneous rocks derived from fractional crystallization of geochemically ultra-depleted mafic melts, which were produced by partial melting of the highly refractory lithospheric mantle during a series of seafloor spreading initiation-failure cycles at a divergent plate margin after the breakup of supercontinent Columbia in the Early Mesoproterozoic. The mafic igneous rocks were located in a passive continental margin in the Late Paleozoic and experienced deep subduction and exhumation in the Triassic, giving rise to the presently studied eclogites. The ancient geochemical signatures were retained without considerable influence by mantle convection, providing insights into the nature of crust-mantle differentiation during the tectonic transition from supercontinental breakup to seafloor spreading beneath the sub-ridge lithospheric mantle.

On the origin of alkali feldspar megacrysts in granitoids. Part 2: evidence for nucleation and growth under magmatic conditions from crystal size distributions of the Cathedral Peak Granodiorite, California, USA

The mechanisms whereby alkali feldspar megacrysts form have been debated for several decades; yet, we do not understand well the processes that lead to their formation. We take advantage of glacially polished outcrop surfaces from the Cathedral Peak Granodiorite in the Tuolumne Intrusive Complex, CA to quantitatively characterize alkali feldspar textures, to provide better insight into their origin. On the glacially polished surfaces, we traced alkali feldspar crystals > 10 mm in the field. From the same localities, we also collected large slabs and stained them to reveal feldspar textures for crystals < 20 mm in size. We scaned the resulting field tracings and rock slabs to quantify CSDs using image processing techniques with the software ImageJ. The CSDs from glacially polished outcrop surfaces and complementary polished and stained rock slabs reveal two stages of crystallization. Crystals > 20 mm show log-linear CSDs with shallow slopes, suggesting magmatic nucleation and growth on timescales of thousands of years. Crystals < 20 mm define a second stage of crystallization, with much steeper slopes, suggesting a period of enhanced nucleation leading to formation of a groundmass during the final stages of solidification on timescales of decades to centuries. We do not find any evidence for CSDs affected by textural coarsening, or any effects of subsolidus processes. Our data suggest that these megacrysts form in large, slowly cooling magma, where low nucleation rates dominate. These crystals are not special in their magmatic formation—only in their size. A change in solidification conditions led to the formation of a groundmass, which warrants further study to better understand this crystallization stage in a plutonic environment.

Constraints on the source of Siberian Trap magmas from Mo isotope evidence

The relative contributions of asthenospheric mantle, lithospheric mantle, and continental crust in the genesis of the Siberian Traps Large Igneous Province (ST-LIP) remain poorly constrained. Most models invoke partial melting of asthenospheric mantle within a mantle plume with an inventory of recycled crustal material, with or without melting of subcontinental lithospheric mantle, and crustal contamination during ascend through the continental lithosphere. A greater understanding of this topic is of fundamental importance because the ST-LIP basalts are associated with the large Permian-Triassic extinction and the world’s largest magmatic Ni-Cu-platinum group element sulfide resources (Ni-Cu-PGE), the Norilsk-Talnakh mining camp. The ∼250 ± 2 Ma Siberian Traps at Norilsk contain a classic sequence of basaltic rocks that provide a spatial and stratigraphic context for changing chemistry of the eruptive products. We present a detailed geochemical and Sr-Nd-Hf-Mo isotopic investigation of ST-LIP volcanic rocks and associated sedimentary rocks, including coal and anhydrite, from the Norilsk area. The 98Mo/95Mo isotope ratios (reported as δ98Mo ratio relative to NIST SRM 3134) vary significantly from −0.62 to 0.07 ‰ for the older basalt formations and −0.41 to 0.03 ‰ for the younger basalt formations. The range of δ98Mo and its correlation with the other geochemical tracers cannot be explained by post-magmatic alteration, magmatic differentiation, or sulfide fractionation as Mo behaves as a lithophile element under these magmatic conditions. We suggest that the δ98Mo range can be explained by interaction of the plume with subcontinental lithospheric mantle modified by subduction processes. This is particularly prominent in the earlier Ivakinsky to Gudchikhinsky formations where light δ98Mo values coupled with low Mo/Ce can be explained by contributions from a dehydrated eclogitic component, whereas more rare heavy δ98Mo values and high Mo/Ce likely require contributions from a fluid metasomatized mantle source. Later magmas of the Nadezhdinsky formation show clear evidence of crustal contamination in their combined Mo-Sr-Nd-Hf isotope and trace element systematics, while the later more voluminous Morongovsky type magmas are shallower, large degree melts with limited crustal interaction. Our data shows the usefulness of Mo isotopes in deciphering magma sources of large igneous provinces eruptions.

Multi-stage magmatism and Sn-polymetallic mineralization in the Shuicheng Region, SW China

The economically important tin‑tungsten mineralization in southwestern China has attracted much interest in recent years, yet its age, tectonic setting, and genetic relationship to the host granites remain uncertain in most regions. Our focus in this study is the Shuicheng tin-polymetallic deposit, situated in the western sector of the Diantan Batholith, within the northern Tengchong Block. The granites in this area consist of two intrusive units, i.e., 1) biotite monzogranite (T1), with LA-ICP-MS zircon UPb ages ranging from 76.88 ± 0.43 Ma to 74.69 ± 0.57 Ma, whole-rock εNd(t) values ranging from −9.75 to −8.11, and zircon εHf(t) values from −12.19 to −8.85; and 2) alkali-feldspar granite (T2), with an LA-ICP-MS zircon UPb age of 52.16 ± 0.23 Ma, whole-rock εNd(t) values of −10.15 to −8.25, and zircon εHf(t) values of −9.64 to −8.06. These two intrusive events have been linked to prolonged low-angle subduction and retreat of the subducting Neo-Tethys Ocean Plate, respectively. The Shuicheng granites exhibit high concentrations of SiO2 and Al2O3, contrasted with low levels of FeOT, CaO, MgO, TiO2, and P2O5, suggesting that they are highly fractionated S-type granites. High abundances of Rb, K, U, and Hf, along with low contents of Ba, Sr, and Ti, further indicate significant fractionation and crystallization processes during magma evolution. Importantly, a transition from purely magmatic conditions to a magmatic-hydrothermal state occurred during the latter stages of magmatic evolution. This transition, associated with an amplified potential for tin mineralization, provides fresh insight into the role of magmatic processes in ore genesis, enriching our understanding of the tin‑tungsten mineralization history of the Tengchong region.

Decoding degassing modes of magma chamber of arc volcanoes: Insights from CO2/Cl and Cl/H2O ratios of magmatic fluids in groundwater

The H2O-CO2-Cl composition of the fluids released from the magma chambers for arc volcanoes was calculated based on the solubilities of H2O and CO2 in the silicate melt and the partition coefficient of Cl between the aqueous fluid and melt for the stages defined by a simple evolution model of the magma chamber. The degassing modes consisted of mafic bubbling (MB: first boiling) and felsic bubbling (FB: second boiling) releasing fluids as separate bubbles due to the mafic magma supply and felsic magma formation by differentiation, respectively, and felsic solidification (FS) releasing fluid due to felsic pluton formation. The variations in the CO2/Cl and Cl/H2O ratios of the fluids were examined for each degassing mode at pressures ranging from 150 to 400 MPa and were found to have very large variations at 0.00024–845 for CO2/Cl and 0.00043–0.041 for Cl/H2O. The CO2/Cl is concluded to be a good indicator of magma chamber conditions because each mode has a specific ratio distinguished from that of the other modes; i.e., CO2/Cl > 43 for MB, 0.027–23 for FB, and < 0.62 for FS.The alteration of magmatic composition of fluid during ascent was discussed considering the phase separation of fluid, CO2 bubbling in groundwater, dissolution, and deposition of minerals; furthermore, the limitation of the method when applied to natural groundwater systems were outlined. This method was tested for a natural volcanic system, the Kuju volcano in central Kyusyu, Japan, which is an active volcano complex composed of various magma types from rhyolitic to basaltic, replacing the eruption centers from west to east. Groundwater samples from springs and boreholes were collected and analyzed for chemical and isotopic compositions to determine the magmatic H2O, dissolved inorganic carbon (DIC), and Cl concentrations. The new approach to determine the magmatic CO2 concentration in groundwater using 3He concentration proposed herein cancelled the effect of carbonate deposition/dissolution, addition of DIC from organic matter, and bubbling of CO2. The estimated magmatic CO2/Cl ratios were 0.0009–14, which overlapped consistently with the modeled variation in the CO2/Cl ratio. The spatial variation of magmatic CO2/Cl in groundwater showed that higher-CO2/Cl are found near the central part of the volcanic complex which composed of felsic to mafic volcanoes, and low-CO2/Cl are only found around old felsic volcanoes. The CO2 bubbling springs were found to be associated with younger mafic volcanoes. The condition of the magma chambers, as shown by the spatial variation in CO2/Cl, were consistent with the record of the eruptive activities of the Kuju volcano. The CO2/Cl method proposed in this study using groundwater is useful for estimating magmatic conditions in chambers.

The influence of olivine settling on the formation of basaltic cumulates revealed by micro-tomography and numerical simulations

The speed at which crystals settle in magmatic reservoirs affects the solidification rate of magmas and their differentiation. Despite extensive prior work on the subject, most of our quantitative understanding of the process is still restricted to treating crystals as spherical particles and does not address the geometric complexities of natural crystals. Here, we use three-dimensional (3D) X-ray microcomputed tomography (X-μCT) observations on olivine crystals from Kīlauea Volcano (Hawai‘i) to document their highly intricate and variable geometries, and the textural growth relationships between the olivine crystals and their inclusions (melt, spinel, and fluid/vapor bubbles). Olivine crystals generally have clustered polyhedral or skeletal shapes, which reflect variable magmatic conditions (or growth rates) during their formation. The cumulative presence of spinel, melt, fluid, and vapor inclusions affects the density of the host crystals by up to 6% relative, and thus plays a limited role on modifying crystal settling rate. In contrast, the overall crystal shape plays a major role. We performed numerical simulations employing a finite element method to investigate the effect of crystal morphology on settling rate and the evolution of the particle volume fraction in a magmatic convective layer. We show that for all olivine geometries investigated, the settling velocity is highest when the long axis of the crystal is aligned with the flow direction of the melt. Increasing the aspect ratio of the olivine tends to decrease its settling velocity and results in an increase in the influence of its orientation on the terminal velocity (UT). Extrapolation of the simulation results to variable particle volume fractions (Φ = 0–0.5) indicates high crystal settling rates (UT = ∼9.7 × 10−6–1.4 × 10−5 m/s) that are used to estimate the timescales for the formation of olivine cumulates in natural melt-dominated basaltic systems. The formation of olivine cumulates is therefore rapid, potentially leading to the accumulation of a crystal layer at the bottom, where the frictional contacts between the crystals exert a rheological lock-up acting against further convection. Crystal accumulation in the locked layer (parametrized with the solid/liquid volume ratio in the reservoir) is a function of the reservoir size and crystal fraction, and takes a few years in small reservoirs (<1 km thick) and a few decades in larger reservoirs (several kms thick). We propose a calibration of olivine suspension timescales for mafic magma reservoirs (based on the knowledge of the particle volume fraction, reservoir height, and olivine morphology). This calibration is used to estimate the rate of cumulate build-up, and can help interpret crystal size distributions in the framework of crystal suspension times.

Late Ediacaran post–collisional granite in Babouri–Figuil (Northwest Cameroon): Implication for crustal–mantle contributions in an extensional setting

The Babouri–Figuil Magmatic Complex (BFMC, adjacent to the central-south Chad domain) is located on the southern flank of the Saharan Metacraton. It is amongst the least geologically studied regions in Cameroon; only a few key research papers on geochronology and isotopes have been completed. The emplacement of granitic bodies is poorly constrained between the early and late–Neoproterozoic with inherited zircons of older Paleoproterozoic–Neoproterozoic ages. In this study, undeformed biotite granite massifs (Lombel and Badesi) were studied, and their geochronologic, geochemical and mineral chemical data, including Nd–Sr–Hf isotopes, are presented. The massifs are enriched in pink K-feldspar (orthoclase and microcline; 60–65 wt%) and plagioclase (8–10 wt%) anastomosing to intergranular quartz (18–20 wt%) and biotite (3–5 wt%) crystals. The Lombel massif yielded a Concordia age of 566 ± 2.5 Ma (MSWD = 1.6, zircon U–Pb), slightly older than the Concordia age of 556 ± 2.7 Ma (MSWD = 2.1, zircon U–Pb) obtained for the Badesi massif. These are interpreted as emplacement ages, which post–date arc magmatism by ∼ 45 Ma. The geochemistry of the massifs showed sub- to mid–alkaline, high–K alkali–calcic, and weakly peraluminous compositional characteristics. The calculated biotite temperatures ranged from 613 to 647 °C, pressures from 100 to 300 MPa, and H2O from 3.0 to 3.5 wt%. Whole–rock composition geothermobarometry revealed low pressure of 100–500 MPa, depths of 3.5–18 km and low oxygen fugacity (ΔQFM = −0.6), indicating reducing crystallization conditions. The estimated radiogenic isotopic values are similar for both massifs, falling in the range from initial εHf(t) = +2.6 to +8.4 and initial εNd(t) = +1.3 to +2.1, with a corresponding average initial Sr of 0.70212–0.70458. These radiogenic isotope values suggest a mantle–derived origin for this granite or the melting of a juvenile mafic crustal protolith, as confirmed by the two–stage model age obtained (TDM2 = ∼ 1.2 Ga). Given that biotite is interstitial in alkali feldspar and plagioclase in both massifs, this suggests that biotite chemistry records conditions of magma evolution. Overall, these results confirm that the studied massifs are post–collisional products emplaced at shallow crustal levels under evolved magmatic conditions. They are coeval with a few post–magmatic intrusions in the central–south Chad and emplaced after the Sao–Francisco Congo Craton collided with the Saharan Metacraton. The ages obtained enabled a possible correlation with several Neoproterozoic post–collisional granites in Cameroon, south-central Chad, Central African Republic and Borborema Province of Brazil. The ages also revealed that continental assembly coupled with regional late–magmatism occurred during the Proterozoic.

Chemistry, Growth, and Fate of the Unique, Short‐Lived (2019–2020) Water Lake at the Summit of Kīlauea Volcano, Hawaii

AbstractLess than a year after the 2018 Kīlauea caldera collapse and eruption, water appeared in newly deepened Halemaʻumaʻu crater. The lake—unprecedented in the written record—grew to a depth of ∼50 m before lava from the December 2020 eruption boiled it away. Surface water heightened concerns of potential phreatic or phreatomagmatic explosions but also offered a new means of possibly identifying eruption precursors. The U.S. Geological Survey Hawaiian Volcano Observatory (HVO) monitored the lake via direct visual observation, webcams, thermal imaging, colorimetry, and laser rangefinders. HVO also employed uncrewed aircraft systems to sample the water and measure near‐lake gas composition. The lake's δD and δ18O indicate a groundwater source with substantial evaporation. The initial sample had a salinity (total dissolved solids concentration) of 71,000 mg/L and was rich in sulfate (∼53,000 mg/L), iron (∼500 mg/L), and magnesium (∼10,000 mg/L). Subsequent samples were slightly more dilute. The water's pH (∼4), δ34S (+4.3‰), and surface temperatures (up to 85°C) suggest, rather than significant scrubbing of magmatic volatiles, leaching of basalt and reactions with sulfate minerals resulted in high concentrations of sulfate and other solutes. Thermodynamic modeling and precipitate mineralogy indicate that water composition was controlled by iron oxidation and sulfate dissolution. Although the lake exhibited no detectable precursors before the next eruption, and phreatic or phreatomagmatic explosions did not materialize, our multi‐parameter approach to monitoring yielded an enhanced understanding of the hydrologic, geologic, and magmatic conditions that led to the formation of the unique and short‐lived lake.

Volatile-rich and oxidized ore-forming magmas of the Zijinshan porphyry Cu deposit, China: A mineral composition perspective

The Zijinshan orefield (Fujian Province) is a representative Cu-Au-Mo-Ag porphyry-epithermal mineralization system in eastern China. It contains both Cu ore-bearing and ore-barren granitoids, which provide an opportunity to study the unique properties of fertile magmas. This study focuses on chemistry of minerals such as zircon, apatite, amphibole, and plagioclase in ore-related (zircon U–Pb ages of 103–101 Ma) and ore-barren (zircon U–Pb ages of 157–156 Ma and ∼105 Ma) granitoids to constrain the characteristics of mineralized magmas. The results of whole-rock geochemistry and zircon Hf isotopes suggest that both ore-barren and ore-related granitoids were derived from partial melting of the lower crust, with the ore-related magma being influenced by mantle materials. The analysis of mineral chemistry reveals that the Zijinshan Cu ore-bearing porphyry magmas have high oxidation states (ΔFMQ +0.23–+2.25), high water contents (4.4–4.9 wt%), and shallow emplacement depths (1.9–4.0 km). These magmatic conditions are favorable for fluid exsolution and mineralization, similar to the global Cu-fertile porphyry magmas, but distinct from regional pre-mineralization ore-barren igneous rocks. Apatites from the ore-related granitoids have significantly lower Cl contents (0.03–0.83 wt%) and higher F/Cl ratios (1.52–377) compared to coeval ore-barren monzogranites (Cl = 0.14–2.42 wt%; F/Cl = 0.07–2.82). Additionally, the apatites’ volatile budgets, including XOH, XCl and XF values, and their covariation relationship, in the ore-related intrusions are consistent with apatite crystallization in H2O-saturated conditions. This suggests that extensive fluid exsolution from the parent magma in the chamber plays a crucial role in the formation of porphyry Cu deposits from a low-Cl magma. These low-Cl-apatite intrusions have the potential to form porphyry deposits. Together, we propose that the high oxygen fugacity, high water contents, and fluid exsolution of the parent magma may have triggered the Zijinshan porphyry mineralization.

GEODYNAMICS OF JOINT ZONES OF REYKJANES AND KOLBEINSEY SPREADING RIDGES WITH ICELAND RIFT ZONES

Joint zones of Reykjanes and Kolbeinsey spreading ridges with Iceland large igneous province considerably differ in structure from adjacent spreading segments despite of similar kinematics. Tjörnes transform zone is a complicated system. It comprises several volcanic and amagmatic structures. In contrast, Reykjanes rift zone has relatively simple structure and homogeneous type of tectonic and magmatic activity. The causes of those differences and their modern dynamics are not fully explained hitherto. Basing on morphometric analysis of normal fault scarps parameters it was concluded that such significant differences of transform zones are result of spatial and temporal stability of adjacent structures. In turn, the latter is controlled by periodic increase of Iceland plume magmatic activity. Rift structures development within transform zones directly correlates with their position over adjacent spreading segments and their magmatic conditions. Modern development of both transform zones is caused by Iceland rift zones instability and migration impacted by Iceland plume thermal pulses. Consequently, transform zones undergo kinematic changes. For Tjörnes transform zone it is expressed in its structure gradual simplification: the western branch and block structures cease their activity. In Reykjanes rift zone rift axis gradually migrates southwards that probably results in its intensive volcanism.

Reduced post-collisional magmas for Tengchong volcanoes, southeastern Tibetan Plateau

Rapidly cooled and well-preserved magnetite and ilmenite in young volcanic rocks provide good samples for constraining the magmatic redox conditions. We report chemical compositions of magnetite and ilmenite from the youngest Tengchong volcanoes at the southeastern margin of the Tibetan Plateau, to constrain the oxygen fugacity and temperatures for the Tengchong magmas. Oxygen fugacity values (△FMQ) relative to the fayalite-magnetite-quartz (FMQ) buffer are −1.8 ∼ −0.7 (averaging −1.1 ± 0.4) for Dayingshan; −1.3 ∼ −0.7（averaging −1.1 ± 0.3） for Maanshan; and −1.3 for Heikongshan. The temperatures for magnetite-ilmenite pairs of Tengchong lavas range from 789 to 1090 °C. The relatively reducing oxygen fugacity values (△FMQ = −1.8 ∼ −0.7) of the Tengchong magmas are significantly lower than those oxidizing arc magmas (△FMQ = +0.5 ∼ +3.5). The reduced magmas at Tengchong may exemplify other Quaternary post-collisional lavas of the Tibetan Plateau. The crystallization temperatures from mineral-liquid equilibrium are 1160–1192 °C for olivine, 1049–1110 °C for plagioclase, 978–1122 °C for clinopyroxene, 789–1090 °C for Fe-Ti oxides, 743–792 °C for zircon. The high Fe-Ti oxide temperature and Ti-in-zircon temperature indicate warm magmatic conditions (>750 °C) prior to eruptions of Tengchong magmas.

Zircon and the role of magmatic petrogenesis in the formation of felsic-hosted volcanogenic massive sulfide (VMS) deposits: a case study from the mid-Paleozoic Yukon-Tanana terrane, northern Canadian Cordillera

Magmatism is a critical component in sustaining hydrothermal convection and metal transport during the formation of volcanogenic massive sulfide (VMS) deposits. Previous studies of magmatic petrogenesis in VMS systems have demonstrated that ore-related volcanic rocks have distinct whole-rock geochemical and isotopic signatures (i.e., high HFSE, REE, Th, εHf-Nd, zircon saturation T) relative to barren volcanic rocks, which supports models of elevated crustal heat flow during periods of ore deposition; however, the petrologic characteristics and intrinsic parameters (e.g., T, fO2) related to these magmatic events in VMS districts remain poorly understood. Arc–back-arc assemblages from the mid-Paleozoic Yukon-Tanana terrane are well-characterized and include the Finlayson Lake VMS district, which is host to several felsic-hosted deposits (e.g., Kudz Ze Kayah, GP4F, Wolverine) that were generated in a peri-Laurentian continental back-arc tectonic setting. In this study, zircon from back-arc and coeval arc rocks in the Yukon-Tanana terrane was used as a proxy for primary magma formation conditions that generated VMS-proximal and VMS-distal stratigraphy. Our results indicate that zircon grains in VMS-proximal environments have unique textural, geochemical, and isotopic characteristics (e.g., low-aspect ratios, greater abundance of zircon-phosphate intergrowths, Th/U > 1, Zr/Hf > 80, Tzrc > 780 °C, εHfi > –7) that are clearly distinguished from zircon in VMS-distal rocks in both the back-arc and arc settings (Th/U < 1, Zr/Hf < 80, Tzrc < 780 °C, εHfi < –7). These signatures correlate to VMS-proximal magmas that were hotter, less fractionated, and contained greater juvenile melt contributions compared to VMS-distal magmas and reflect a series of high-flux magmatic events that directly correspond to the early tectonic development of Yukon-Tanana terrane. Moreover, this study underscores the importance of mineral-scale petrology, geochemistry, and geochronology in defining the primary magmatic conditions that generated VMS-related felsic rocks and highlights the utility of zircon as a prospectivity tool in both grassroots and brownfields VMS exploration.

Hot subduction in the southern Paleo-Asian Ocean: Insights from clinopyroxene chemistry and Sr–Nd–Hf–Pb isotopes of Carboniferous volcanics in West Junggar

The chemical evolution and pressure–temperature conditions of subduction zone magmatism along ancient suture zones in orogenic belts can provide important information regarding plate convergence processes in paleo-oceans. Carboniferous magmatism in West Junggar is key to understanding the tectonothermal and subduction history of the Junggar Ocean, which was a branch of the Paleo-Asian Ocean, as well as the accretionary processes in the southwestern Central Asian Orogenic Belt (CAOB). We undertook a geochronological, mineralogical, geochemical, and Sr–Nd–Hf–Pb isotopic study of volcanic rocks from the Baikouquan area of West Junggar. We used these data to determine the petrogenesis, mantle source, and pressure–temperature conditions of these magmas, and further constrain the subduction and tectonic history of the Junggar Ocean. The studied volcanic rocks yielded zircon U–Pb ages of 342–337 Ma and are characterized by enrichments of large-ion lithophile elements (LILEs), and depletions in high-field-strength elements (HFSEs), indicative of an island arc affinity. The volcanic rocks have positive ƐNd(t) (5.83–7.04) and ƐHf(t) (13.47–15.74) values, 87Sr/86Sr(t) ratios of 0.704023–0.705658, and radiogenic 207Pb/204Pb(t) and 208Pb/204Pb(t) ratios at a given 206Pb/204Pb(t) ratio, indicative of a depleted mantle source contaminated by subduction-related materials. Geochemical modeling calculations indicate that ≤1% of a subduction component comprising fluid and sediment melt could have generated the source of the parental melts of the Baikouquan volcanic rocks. Clinopyroxene phenocrysts in the volcanic rocks are classified as high- and low-Ti clinopyroxene, and pressure–temperature calculations suggest the host rocks formed at high temperatures (∼1300 °C) and shallow to moderate depths (<2 GPa). The magma was probably generated by hot and hydrous melting in a mantle wedge in response to subduction of young, hot oceanic lithosphere. The present results, combined with published data, suggest that the Baikouquan volcanic rocks record a transition in tectonic setting from normal cold to anomalous hot subduction of young oceanic lithosphere close to a mid-ocean ridge. This indicates ridge subduction began shortly after 337 Ma. Our results provide new insights into the tectonomagmatic evolution during intra-oceanic subduction prior to ridge subduction.

Natural Ti‐Rich Mineral (Ilmenite, Titanite and Perovskite) Reference Materials for In Situ Ti Isotopic Measurement by LA‐MC‐ICP‐MS

Iron‐titanium oxides such as ilmenite (FeTiO3), titanite (CaTiSiO5) and perovskite (CaTiO3) are the common Ti‐rich mineral phases crystallised during magmatic and metamorphic processes on Earth. Depending on magma types or conditions of phase equilibria, formation of these Ti‐rich minerals can result in Ti isotopic fractionation in the range of ‐1.52‰ to +2.90‰ on δ49TiOL‐Ti (i.e., the per mil difference of 49Ti/47Ti ratio relative to the OL‐Ti reference material), making the Ti isotope ratios of these minerals potential tracers for conditions of magmatism and metamorphism. Due to their resistance to aqueous alteration, these Ti‐rich accessory minerals are also commonly present as pristine, detrital phases in sedimentary rocks, which offer an opportunity to study the evolution of magmatism and metamorphism throughout the Earth's geological history. Here we have developed a novel technique for in situ Ti isotopic measurement in ilmenite, titanite and perovskite using femtosecond laser ablation multi‐collector inductively coupled plasma‐mass spectrometry (fs‐LA‐MC‐ICP‐MS) under wet plasma conditions. Samples were ablated with different laser spot sizes (15–50 μm) and different laser energy densities (0.6–4.2 J cm‐2) to obtain adequate Ti signal intensity at a fixed laser repetition rate of 2 Hz. When 49Ti signal intensity of samples ranged from 0.3 to 3.7 V, no significant signal‐dependent Ti isotopic fractionation was observed under wet plasma conditions. Repeated measurements on twelve Ti‐rich minerals using different analytical protocols provided comparable δ49TiOL‐Ti values within uncertainties, confirming the accuracy of the proposed fs‐LA‐MC‐ICP‐MS method. With one exception (i.e., RUS1), all analysed minerals are homogeneous in their Ti isotopic compositions between individual chips with an intermediate precision of ±0.13‰ to ±0.17‰ (2s) on δ49TiOL‐Ti. Collectively, these minerals record a significant δ49TiOL‐Ti variation ranging from ‐0.46‰ to +2.12‰. These warrant the suitability of these materials as Ti isotopic reference materials for in situ Ti isotopic measurement.

Priming and eruption of andesite magmas at Taranaki volcano recorded in plagioclase phenocrysts

Plagioclase phenocrysts in pre-historic andesites provide insight to the dynamics of magma formation and eruption at Taranaki volcano, New Zealand. The phenocryst population has a diversity of relic cores and a total in situ 87Sr/86Sr range of 0.70440–0.70486. Within-sample 87Sr/86Sr variations of 0.00018 to 0.00043 indicate that many phenocrysts are antecrysts and/or xenocrysts, derived from multiple crystal mush bodies. The Sr-isotopic differences in the phenocrysts of consecutive eruptions indicate that different magmas were tapped or formed on a centennial timescale. Most phenocrysts have multiple resorption/calcic regrowth zone(s) with elevated FeO* but invariant MgO zonation profiles. They likely record mafic melt inputs, and subsequent storage at elevated temperature caused re-equilibration of the Mg gradient. However, distinct rim types record different final pre/syn-eruptive magmatic conditions. Those in magmas erupted at 1030–1157 CE, 1290–1399 CE and 1780–1800 CE are characterised by resorption and calcic regrowth with sharp MgO and FeO* gradients. They record the entry of mafic melt into the system a few days or less before eruption based on Mg diffusion chronometry. In contrast, most phenocrysts erupted at 1755 CE, 1655 CE, and a few pre-1 ka events, have texturally uniform rims, compositionally consistent with closed-system crystallisation. This suggests alternating external and internal eruption triggers. Alternatively, the rate of magma reactivation via intrusion may dictate whether there was sufficient time for a mineralogical response to be recorded in part or all of the system. With respect to anticipating future eruptions, the plagioclase phenocrysts suggest multi-stage magma priming but rapid onset of the final trigger.

Aso volcano, Japan: assessing the 100-year probability of a new caldera-forming eruption based on expert judgements with Bayes Net and Importance Sampling uncertainty analysis

The Aso-4 explosive eruption on Kyushu, Japan, 89,500 years ago was one of the biggest eruptions in the last one hundred millennia, with a magnitude of approximately M8. Modern society requires the likelihood of natural events with potentially disastrous consequences to be evaluated, even if probabilities of occurrence are diminishingly small. For some situations, it is not satisfactory to assert an event scenario probability is “negligible” or can be “ignored”. Judicial hearings or litigation may require risk levels to be quantified, in which case, statements of scientific confidence could be decisive. Internationally, e.g., for nuclear site safety evaluations, event likelihoods on order of 10–7/year are often considered for quantitative assessment. At such hazard levels, this might include evaluating the proposition that a particular volcano can deliver a future super-eruption, a supposition that could be attached to Aso volcano. But, simplistically taking the average recurrence interval between past caldera-forming eruptions at a given volcano is an unreliable guide to the likelihood of a future repeat: each past event represented a unique set of tectonic and magmatic conditions within a continually evolving volcanic system. Such processes are not temporally stationary nor statistically uniform. To evaluate the probability of a new M8 event at Aso, within the next 100 years, we performed a comprehensive stochastic probability uncertainty analysis using a model implemented with advanced computational Bayes Net (BN) software. Our eruption process model is informed by multiple strands of evidence from volcanology, petrology, geochemistry and geophysics, together with estimates of epistemic (knowledge) uncertainty, adduced from reviews of published data, modelling and from expert judgement elicitation. Several lines of evidence characterise the likely structure, magmatic composition and eruptive state of the present-day Aso volcano, which has had numerous smaller eruptions since Aso-4. To calculate the probability of another M8 eruption of Aso, we implemented probabilistic ‘Importance Sampling’ in our model. With this approach, we find the chance of an Aso-4 scale eruption (characterised by mean volume 500 km3 DRE and approximate 90% credible interval [210 ‥ 1200] km3 DRE) is less than 1–in–1 billion in the next 100 years (i.e., < 10–9 probability). Based on current volcanological understanding and evidence, we believe this probability estimate is robust to within an order of magnitude.

Iron Isotope Systematics of the Skaergaard Intrusion and Implications for its Liquid Line of Descent

Abstract The Skaergaard intrusion is one of the most thoroughly studied layered mafic intrusions on Earth and an exceptional example of (near) closed-system magmatic differentiation. We report new Fe isotope data for whole rocks, and magnetite and ilmenite separates through the layered series (LS) and upper border series (UBS) of the intrusion. δ56Fe values for gabbroic rocks range from 0.033‰ to 0.151‰ with an abrupt step increase at the base of Lower Zone c (LZc) within LS with the appearance of cumulus magnetite and subsequent decline accompanying FeTi oxide fractionation. The lowest δ56Fe values are found near the Upper Zone b (UZb) to c (UZc) boundary followed by a sharp rise across UZc approaching the Sandwich Horizon. Magnetite–ilmenite separates straddle bulk rock compositions with fractionation factors (Δ56Femt-ilm) of 0.081‰ to 0.239‰, consistent with subsolidus equilibration. Granophyric rocks occurring as pods, sheets and wispy layers from the upper zone and UBS equivalents and having unradiogenic Sr similar to gabbroic rocks of Skaergaard, are isotopically heavier than their host ferrodiorites (Δ56Fegranophyre-ferrodiorite ≥ 0.1‰) reaching a maximum δ56Fe of 0.217‰ in UBS. A fused xenolith from UBS has δ56Fe = 0.372‰. This range in δ56Fe spans much of that reported for terrestrial igneous rocks, and like the global dataset, shows a pronounced increase in δ56Fe with inferred silica content of modeled Skaergaard liquids. Forward modeling of closed system fractional solidification was undertaken to account for Fe isotope systematics, first by testing published liquid lines of descent (LLD), and then by exploring improvements and considering the impacts of liquid immiscibility, crustal contamination, fluid exsolution and diffusional processes. Our modeling relies on published Fe+2 and Fe+3 force constants for magmatic minerals and silicate glasses, and the most reliable estimates of the average bulk composition and mass proportions of the well-defined subzones of the intrusion. We show that the increase in δ56Fe across the LZb–LZc boundary is readily explained by the increased incorporation of Fe+3 into the crystallizing solid including magnetite. We further demonstrate that the classic Fenner LLD, involving strong Fe enrichment at nearly constant silica, does not lead to a rise in δ56Fe toward the end stages of evolution, while a Bowen-like LLD, with little Fe enrichment and strong Si enrichment, also underestimates enrichment in heavy Fe isotopes in the ferrodiorites of UZc. A LLD following an intermediate path involving modest Fe and Si enrichment, followed by Fe depletion best explains the observations. We predict ~3.5% (by mass) residual liquid after crystallization of UZc having a composition similar to felsic segregations in pegmatitic bodies found in the intrusion. While liquid immiscibility may have been encountered within fractionating mush at the margins of the intrusion, the Fe isotope systematics do not support liquid phase separation of the bulk magma. Crustal contamination, fluid exsolution, hydrothermal alteration and thermal diffusion are also shown to have no resolvable effect on the Fe isotope composition of the gabbroic and granophyric rocks. We conclude that the Fe isotope systematics documented in the Skaergaard intrusion reflect the dominant role of fractionating Fe-rich minerals from gabbroic through ferrodioritic to rhyolitic liquids. The success of our model to account for the observed Fe isotope systematics for Skaergaard demonstrates the utility of Fe+2 and Fe+3 force constants determined at ambient conditions to model magmatic conditions and gives critical insights into plutonic processes fractionating Fe isotopes complementary to the volcanic record.

Diffusion of halogens (F, Cl, Br, I) in silicic melt

Chemical diffusion of the halogens F, Cl, Br, and I in silica-rich natural melts was experimentally investigated by the diffusion couple technique. Experiments were conducted under anhydrous conditions at atmospheric pressure and hydrous conditions (∼1.5 wt% H2O) at 160 MPa, over a temperature range of 750–1000 °C and 1000–1200 °C, respectively. Quenched trachytic melt samples were analyzed using an electron microprobe (EPMA) and secondary ion mass spectrometry (SIMS).All halogens exhibit Arrhenian behaviour during diffusion in the investigated melt compositions with F always diffusing fastest. The other halogens show progressively slower diffusion (F > Cl > Br > I) correlated to their ionic radii. In anhydrous melt a diffusivity range of 3–4 orders of magnitude is covered among the halogens with DF(1000 °C) ∼5 × 10−13 m2s−1 and DI(1000 °C) ∼1 × 10−16 m2s−1. The diffusivities of all halogens increase in hydrous melt with the largest increase being observed for the slowest-diffusing halogens, e.g., resulting in an increase of up to 2 orders of magnitude for iodine at 1000 °C compared to the anhydrous case. This behavior yields a narrower overall diffusive range of only 1–2 orders of magnitude among all halogens. Activation energies (EA) of diffusion consistently range from ∼200 to 390 kJ mol−1 in anhydrous melts. In hydrous melt EA generally decreases, with the highest decrease determined for F (∼131 kJ mol−1) and only slight changes for the other halogens (∼201–222 kJ mol−1).Our diffusivity data of the anhydrous series exhibit a pronounced correlation of diffusivity with the ionic radii, however, this correlation is attenuated in hydrous melt. While in anhydrous melt, halogen diffusion is closely related to ionic porosity, in hydrous experiments, the process of ionic detachment becomes more important as a rate-limiting diffusion mechanism, e.g., comparable to the case of diffusion of divalent/trivalent cations.The results of this study provide the first consistent diffusion dataset including all halogens under naturally relevant magmatic conditions and highlight the pronounced compositional effect of both, major element and dissolved H2O on halogen diffusion. The derived diffusion parameters may be readily used for modelling of diffusive fractionation in silicic melts or determining the timescales of natural silicic volcanic processes based on halogen concentration measurements. Furthermore, these data emphasize the potential of diffusive fractionation among the halogens, which may be applied as a monitoring tool for volcanic unrest on actively degassing volcanoes.