Magma Batches Research Articles

Mixing of magmas from multiple sources in the petrogenesis of multi‐stage Early Cretaceous syenites in the Wulingshan alkaline complex, northern North China Craton: Evidence from enclaves

There are several different models for the origin of syenites, but the role of magma mixing in the formation of syenites remains unclear. The Wulingshan alkaline complex in the northern North China Craton consists mainly of porphyritic syenite and syenite with abundant enclaves. These enclaves may provide new insights into the petrogenesis of syenites. We obtained zircon U–Pb age, mineral chemistry, whole‐rock major and trace elements, and Sr–Nd isotopic data for the enclaves and their host rocks to constrain the petrogenesis and identify the role of magma mixing during the formation of the different syenites in the Wulingshan alkaline complex. The results of zircon U–Pb dating indicate that the enclaves and host rocks crystallized contemporaneously at ca. 133 Ma. The enclaves contain abundant clinopyroxene, amphibole, and biotite, and their average Nb/Ta (18.46) and Th/Ce (0.04) ratios are similar to those of the mantle. The enclaves have relatively high Fe2O3T, MgO, and CaO contents, and their SiO2contents are equivalent to those of intermediate rocks, indicating that they formed by mixing the mafic and felsic magmas to some extent. The different mineral assemblages, major and trace elements and isotopic compositions of the enclaves in the porphyritic syenite and syenite indicate that they are two batches of parental magma with different properties derived from the mantle. The porphyritic syenite and syenite have high SiO2and low Fe2O3T, MgO, and CaO contents, as well as low V, Cr, Co, and Ni contents. These major and trace element characteristics of the host rocks require the involvement of crustal components. The Sr–Nd isotopic compositions of the enclaves and the host syenitoids plot on a mixing curve between enriched lithospheric mantle and lower crust, indicating that these rocks are the products of magma mixing and crust–mantle interaction. Considering the contact relationships, the geochemistry data, and the mixing model presented in the paper, we propose that the upwelling of the hot asthenosphere heated the overlying enriched lithospheric mantle and triggered low‐degree partial melting. The alkalic mafic magmas derived from enriched mantle sources were mixed with felsic magmas generated by partial melting of the lower (or upper) crust to form the different syenites.

Consecutive underplating of cognate magmas: Contributions to the enclaves and large-volume plagiogranites in the North Qaidam orogen (NW China)

Magma mingling between mafic and felsic magmas is a plausible mechanism for the formation of mafic microgranular enclaves (MMEs) which are common worldwide. Through dissecting the MMEs with parentage of magma mingling, it can be ascertained how two types of magmas with contrasting physicochemical properties interplay and what implicit information magma mingling records in crust-mantle interaction. Here we present petrography, geochronology, major and trace elements, Sr-Nd-Hf isotopes, and mineral geochemistry for rock associations from the North Qaidam orogen, which provides new insights into the processes of magma mingling, petrogenesis of large-volume host plagiogranites and crust-mantle magmatic interaction. MORB (mid-ocean ridge basalt)-like and arc-like magmas can be distinguished from mantle-derived mafic samples with zircon UPb ages ranging from 480 to 465 Ma. MORB-like magmas originated from a MORB-type asthenospheric mantle, whereas arc-like magmas were derived from an ultra-depleted asthenospheric mantle as a result of consecutive extraction of MORB-like magmas and late-stage metasomatism by subduction-related aqueous fluids. Given the content discrepancies for major oxides, gabbro-dioritic MMEs are formed by incipient mechanical mingling between gabbro-diorite magmas and felsic plagiogranite magmas. MMEs and host plagiogranites yield coincident crystallization ages of 478 Ma and overlapped whole-rock SrNd ((87Sr/86Sr)i = 0.7044–0.7052; εNd(t) = +4.6 ∼ +5.6) and zircon Hf isotopes (εHf(t) = +13 ∼ +16). The host plagiogranite magmas are most likely derived by partial melting of pre-existing analogues of gabbro-diorite intrusions under the condition of consecutive underplating of different batches of cognate magmas over short time scales, which in turn determine the consistent isotope compositions between MMEs and host plagiogranites. As for mineral geochemistry, the plagioclases with the highest An% value and the hornblendes implying the highest melt water content, appear within the MME-bearing plagiogranites next to the chilled margin of MMEs, illustrating that aqueous fluid and heat transfer play a pioneering but finite role in incipient magma mingling. The consecutive underplating of cognate mantle-derived magmas as an ultimate engine leads to the partial melting of newly emplaced basites to form massive plagiogranite magmas and geochemically stratified crust. With magma ascent into a shallow level of crust, magma mingling between underplated mafic magmas and resultant plagiogranite magmas is finally realized.

Compositional Evolution of Polygenetic Fissure Volcanic Systems: Insights From the Latest Eruptions at Craters of the Moon Volcanic Field

AbstractThe Craters of the Moon (COM) Volcanic Field in Idaho is the largest mostly Holocene lava field in the conterminous U.S. and the site of the most recent volcanism in the Yellowstone‐Snake River Plain province. This prominent example of a polygenetic fissure field has produced over 60 eruptions over the past 15,000 years from the Great Rift fissure system. The eight most recent lava flows, known as Period A, were erupted between about 2,500 and 2,000 years ago. In this study, major and trace element and Pb‐Sr‐Nd isotope data were obtained, with several samples collected from each flow from widely spaced locations to evaluate heterogeneity in their source magmas. The results reveal three previously unrecognized compositional groups that may represent separate magma batches, and that assimilation and fractional crystallization (AFC) controlled the evolution of the entire Period A suite. The three oldest flows were periodic withdrawals from a body of highly evolved and poorly homogenized magma, and four later flows represent successive samplings from a stratified body of more primitive magma. These two flow groups are related via variable AFC modeled with Magma Chamber Simulator software, with Neogene rhyolite as a plausible crustal assimilant. A third compositional type that requires a different parent magma and AFC history is represented by a single flow erupted from vents a few km away from the others. The compositional evolution and variability observed in this sequence of flows provide new insights into how magmas are generated, modified, and delivered in such polygenetic fissure systems.

Variations in Volcanism and Tectonics Along the Hotspot‐Influenced Reykjanes Ridge

AbstractMapping and sampling four sections of the slow‐spreading Reykjanes Ridge provide insight into how tectonic and volcanic activity varies with distance from the Iceland plume. The studied areas are characterized by significant variations in water depth, lava chemistry, crustal thickness, thermal structure, and ridge morphology. For each study area, fault pattern and dimension, tectonic strain, seamount morphology, and density are inferred from 15 m‐resolution bathymetry. These observations are combined with geochemical analysis from glass samples and sediment thickness estimations along Remotely Operated Vehicle‐dive videos. They reveal that (a) tectonic and volcanic activity along the Reykjanes Ridge, do not systematically vary with distance from the plume center. (b) The tectonic geometry appears directly related to the deepening of the brittle/ductile transition and the maximum change in tectonic strain related to the rapid change in crustal thickness and the transition between axial‐high and axial valley (∼59.5°N). (c) Across‐axis variations in the fault density and sediment thickness provide similar widths for the neo‐volcanic zone except in regions of increased seamount emplacement. (d) The variations in seamount density (especially strong for flat‐topped seamounts) are not related to the distance from the plume but appear to be correlated with the interaction between the V‐shape ridges (VSR) flanking the ridge and the ridge axis. These observations are more compatible with the buoyant upwelling melting instability hypothesis for VSR formation and suggest that buoyant melting instabilities create many small magma batches which by‐pass the normal subaxial magmatic plumbing system, erupting over a wider‐than‐normal area.

Genesis of the Questa Mo Porphyry Deposit and Nearby Polymetallic Mineralization, New Mexico, USA

Abstract The Oligocene Latir magmatic center in northern New Mexico is an exceptionally well exposed volcano-plutonic complex that hosts a variety of magmatic-hydrothermal deposits, ranging from relatively deep, F-rich porphyry Mo mineralization to shallower epithermal deposits. We present new whole-rock chemical and isotopic data for plutonic rocks from the Latir magmatic center, including extensive sampling of drill core samples of intrusive rocks from the Questa porphyry Mo deposit. These data document temporal chemical trends of porphyry-related mineralization that occurred after caldera-forming magmatism and during postcaldera batholith assembly. Silicic magmas were generated multiple times throughout the history of the Latir magmatic center, but few are associated with the formation of a mineral deposit. Whole-rock trace element ratios and Sr, Nd, and Pb isotope compositions vary throughout the protracted history of silicic magmatism. The caldera-forming ignimbrite and early phase of postcaldera intrusions are unmineralized, more enriched in high field strength elements, and generally contain less radiogenic Sr and Pb and more radiogenic Nd than later intrusions. The Questa porphyry Mo deposit formed immediately after the most isotopically primitive phase of the batholith was assembled, ruling out simple reworking of juvenile mantle-derived crust as the source for mineralizing magmas. Rhyolite dikes associated with polymetallic sulfide deposits intruded ~800 k.y. after Mo mineralization, and Nd isotope data indicate that these dikes are associated with different batches of magma and are unrelated to the Mo-mineralizing intrusions at the Questa mine. Together, these data indicate that the source of magmas changed significantly throughout the 10-m.y. history of the magmatic center. We assess multiple genetic models for porphyry-related magmatism against this data set, favoring models with discrete periods of magma genesis from a deep hybridized zone in the lower crust giving rise to the punctuated periods of mineralization. These observations suggest that the formation of mineral deposits within a central magmatic locus is likely the result of the piecemeal assembly of individual hydrothermal-magmatic systems, and that distal and younger polymetallic mineralization commonly observed near known porphyry deposits represents decoupled processes.

The 2013–2020 seismic activity at Sabancaya Volcano (Peru): Long lasting unrest and eruption

Sabancaya volcano is the youngest and second most active volcano in Peru. It is part of the Ampato-Sabancaya volcanic complex which sits to the south of the ancient Hualca Hualca volcano and several frequently active faults, thus resulting in complex volcano-tectonic interactions. After 15 years of repose, in 2013, a series of 4 earthquakes with magnitude >4.5 occurred within 24 h, marking the beginning of a new episode of unrest. Several additional swarms of earthquakes occurred in the following years until magmatic eruptive activity started on 6 November 2016. This activity is ongoing as of this writing, with an average of 50 explosions per day. In this study, we present results of multiparametric monitoring of Sabancaya's activity observed during 2013–2020. Seismic data are used to create a one-dimensional seismic velocity model, to catalog, locate, and characterize earthquakes, to detect repeating earthquake families, and to monitor seismic velocity variations by ambient noise cross-correlation. These analyses are complemented by visual and remote sensing observations and ground deformation measurements. All monitored parameters showed significant changes on 6 November 2016, the day of eruption onset, thus dividing the eruptive activity into pre-eruptive and eruptive stages.The unrest is characterized by high levels of seismic activity with hundreds of events detected per day. Volcano-tectonic (VT) earthquakes were dominant during the pre-eruptive period while long-period (LP) events and explosions have been most numerous since the eruption onset. Earthquake locations highlight long-lasting seismogenic zones along multiple previously active regional faults, as well as along newly identified faults. This VT seismicity is mainly distributed in a sector from the northwest to the east of the volcanic complex at distances of up to 30 km from the crater. We focus our analysis on two eruptive episodes: the eruption onset and subsequent crater migration from south to north, and the increase of lava dome extrusion rate in 2019. Both episodes are accompanied by seismic velocity decreases of up to 0.2% and are preceded by a few weeks by bursts of distal VT activity, including numerous repeating earthquakes. These repeated events were located on several remote tectonic faults (5–25 km from the vent). We suggest that these phenomena could be due to the injection of a batch of magma in the deep reservoir and/or conduit, which would generate 1) a pressure wave propagating in the hydrothermal system, triggering the bursts of seismic activity and 2) slow rising of magma by melting old material filling the conduit that eventually produced the eruptive and dome growth acceleration events.

Geochronology and Mapping Constraints on the Time-Space Evolution of the Igneous and Hydrothermal Systems in the Taurus Cu-Mo District, Eastern Alaska

Abstract The Taurus porphyry Cu-Mo district contains four mineralized porphyry centers in the eastern interior of Alaska. All four centers were emplaced during a magmatic episode that spanned from ca. 72 to 67 Ma, with seven distinct igneous suites. Each igneous suite resulted in hydrothermal alteration and mineralization, with younger pulses overprinting older pulses. Each magmatic-hydrothermal system is not present at all four mineralized centers. Apart from the Dennison occurrence, each mineralized center records pulses of repeated intermediate-silicic magmatism and associated alteration and mineralization. Laser ablation-inductively coupled plasma-mass spectrometry U-Pb zircon crystallization ages indicate that an early quartz porphyry dike swarm ranges in age from ca. 71 to 70 Ma and is associated with potassic, sericitic, and propylitic alteration. Quartz latite intrusions were emplaced at ca. 69 Ma and exhibit early sodiccalcic alteration overprinted by potassic, sericitic, and propylitic alteration. The Taurus monzonite suite is cut by quartz latite but yielded an ca. 70 Ma emplacement age and exhibits the largest footprint of potassic and sericitic alteration. Feldspar porphyry dikes were emplaced ca. 69 Ma and have significant tourmaline-bearing potassic and sericitic alteration. This suite was followed by development of an igneous breccia with a monzonitic igneous matrix. Sodic-calcic alteration was associated with the igneous brecciation. A small stock of monzonite was emplaced at ca. 68 Ma causing locally pervasive sericite-tourmaline-pyrite alteration. The youngest suite of magmatism dated in the district is a series of granodiorite porphyry dikes with weak sodic-calcic and propylitic alteration that truncates earlier alteration assemblages. Mineralization in the district consists of chalcopyrite and molybdenite associated with sugary quartz veins with potassium feldspar and biotite alteration envelopes (A veins). Less common banded quartz-molybdenite veins (B veins) occur with potassium feldspar envelopes. Gold occurs throughout the district and is strongly correlated with copper grade. Sericitic alteration contains lower copper contents and is predominantly associated with quartz-pyrite veins with sericite envelopes (D veins). Pyrrhotite and local arsenopyrite are present in sericitic assemblages. Pyrrhotite also occurs as inclusions in pyrite within D veins. Magmas across the district exhibit oxidized characteristics, evidenced by the presence of abundant magnetite, rare titanite, and elevated Eu/Eu* and Ce/Ce* in zircon. Zircon Th/U and Yb/Gd compositions suggest a fractionation path controlled by apatite, titanite, and hornblende. Zircon rare earth element ratios and trace element data indicate two distinct batches of magma evolved from mafic parental compositions to monzonite and granodioritic compositions via fractional crystallization. In the early pulse of magma (ca. 72–69 Ma), fractional crystallization was key to ore formation. Earlier, better mineralized suites evolve to less negative Eu anomalies (Eu/Eu* > 0.7), indicating more oxidized and higher-water-pressure conditions evidenced by the suppression of plagioclase crystallization, compared to later, more poorly mineralized suites. The temporal and spatial evolution of the district was determined from mapping and U-Pb and Re-Os geochronology. Mapping of igneous and hydrothermal assemblages indicates that the locus of the intrusive suites and hydrothermal systems shifted spatially over time, based on the presence of high-temperature (K-silicate–dominant) alteration, which is coincident with the highest Cu and Au grades. The earliest hydrothermal system was centered at Bluff and East Taurus and transitioned to West Taurus during emplacement of the second magmatic suite. Emplacement of the third magmatic suite was centered back at East Taurus, and the fourth and fifth suites were centered at West Taurus. The latest suites were widespread without a core of high-temperature alteration marking a central locus. East Taurus contains the overlap of six of the seven magmatic and hydrothermal suites and has the highest intersected grades and tonnages in the district. The Bluff and Dennison occurrences exhibit fewer igneous suites and hydrothermal assemblages with weak mineralization. Sodic-calcic alteration, common on the deep and distal flanks of porphyry systems, is only present at West Taurus and is indicative of a localized source of high-salinity nonmagmatic fluids.

The volcanic activity changes occurred in the 2021–2022 at Vulcano island (Italy), inferred by the abrupt variations of soil CO2 output

The active cone of La Fossa is a close conduit volcano characterized by solphataric activity, manifested by discharging fluids through fumaroles and soil degassing. Since 1978 several degassing crises have been observed and interpreted as early signals of volcanic unrests. Recently, from June 2021 to May 2022, we measured the changes in soils CO2 release to evaluating the level and duration of the actual exhaling crises. The CO2 output has been evaluated by surveys carried out in anomalously degassing areas, located both in the La Fossa cone summit area and in other peripheral zones, coupled to near-real time monitoring data acquired by three automated stations. The strong and deep input of volatiles released from an underlying magma batch modified the chemical composition of the shallow plumbing system, bringing the system to a higher level of CO2 total pressure. This work highlights that a geochemical networks of stations, located at some distance from the fumaroles release and/or from eruptive conduits, is useful and can be applied to characterizing and monitoring any other active volcanic system. This type of studies can be useful to contribute to forecast the next evolution of the studied systems.

Zircon Hf-O isotopic constraints on the origin and temporal evolution of the Toba volcanic system, Indonesia

Combined U-Pb-Hf-O isotopic analyses of zircons are useful probes to investigate the origin and temporal evolution of long-lived magma reservoirs feeding large eruptions. This study presents a complete dataset of zircon Hf-O isotopes of the latest five Toba eruptions in Quaternary. The new isotopic data differ slightly within and between eruptive units, with the majority of zircons having ε Hf values between −8 to −5 and δ 18 O values between 6 and 6.8. These Hf-O isotopic compositions of zircons were considered to inherit from an enriched and slightly heterogenous magma source. Specifically, zircons of the Oldest Toba Tuff (OTT) have rather restricted and constant Hf-O isotopes, which may crystallize from a reservoir protractedly supplied from the same magma source. Zircons of the Middle Toba Tuff (MTT) have a large range of δ 18 O values varying from ~6 to 8, with low ε Hf values extending to −14, indicating shallow-level crustal contamination. The largely negative correlations between ε Hf and δ 18 O for zircons from the MTT is indicative of contamination of high-δ 18 O crustal materials, whereas those of the positive correlations for zircons from the Youngest Toba Tuff (YTT) implies contamination of hydrothermally altered low-δ 18 O crustal materials. It is important that the Hf-O isotopes of all the zircons do not record any single batch of magma that was derived from partial melting of depleted mantle peridotites, but are consistent with an enriched and heterogenous mantle source hybridized by the subducted sediments, followed by contamination of low-δ 18 O or high-δ 18 O materials during their ascent to the surface or storage in the shallow magma reservoirs. The increase of zircon δ 18 O values towards eruptions in both the MTT and YTT suggests that the two eruptions may be accompanied by roof collapse that added high-δ 18 O materials to the magma. Thus, assimilation of high-δ 18 O crustal materials could be a possible trigger for the MTT and YTT eruptions. • Toba volcanic rocks were derived from an isotopically enriched and heterogenous mantle source with significant crustal contamination. • Hf-O isotopes of MTT zircon crystals exhibit characteristics of shallow-level crustal contamination. • The YTT magma reservoir was contaminated by both high- and low-δ 18 O crustal materials. • Assimilation of high-δ 18 O materials by roof collapse could be a likely eruption trigger for the MTT and YTT eruptions.

UNDERSTANDING Cu DEFICIENCY AND Mo ENRICHMENT IN THE JURASSIC ZHANGGUANGCAI-LESSER XING’AN CONTINENTAL ARC (NE CHINA): INSIGHTS FROM THE LUMING PORPHYRY MO DEPOSIT

Abstract Continental arcs are favorable sites for porphyry Cu ± Mo deposits. In contrast, the Jurassic Zhangguangcai-Lesser Xing’an continental arc, northeastern China, hosts numerous porphyry Mo deposits with only minor total Cu endowment. The reasons for this remain elusive. Here, we address the issue by tracking the origin and evolution of a long-lived (ca. 12 m.y.) magmatic suite associated with one of the largest Mo deposits in the Zhangguangcai-Lesser Xing’an range continental arc—the Luming porphyry Mo deposit. In combination with previous geochronological data, our results indicate that the pre- (187.5–186.5 Ma) and syn-ore (178.6–175.6 Ma) intrusions at Luming represent separate magma batches from isotopically similar parental magmas with the signature of subduction-metasomatized mantle-derived magmas mixed with crustal melts. The Sr/Y, La/YbN, Dy/YbN, and Eu/Eu* values, together with fO2 and apatite S and Cl concentrations of the syn-ore intrusions, however, are systematically higher than pre-ore intrusions. These data indicate that there is a distinct change in magma chemistry during the period of magmatism quiescence between the pre- and syn-ore magmatism. The marked change in magma chemistry is coincident with an episode of major compression and crustal thickening in the Zhangguangcai-Lesser Xing’an range. We, therefore, interpret that the chemical changes most likely reflect a deepening of the locus of lower crustal magma evolution, linked to increased crustal thickness and/or melt H2O contents accompanying an increase in orogenic stress. Considering the tectonic setting, magmatic evolution in the deep crust, emplacement depth, fractionation degree, and oxygen fugacity of the porphyry Mo-related intrusions at Luming, they are very similar to porphyry Cu deposits worldwide. Accordingly, we propose that the deficiency in Cu and enrichment in Mo of the porphyry deposits in the Zhangguangcai-Lesser Xing’an range most likely formed due to mantle-derived magmas interacting with reduced crustal materials or melts in the deep crust. This resulted in substantial sequestration of Cu and S, decreasing porphyry Cu potential. Such an inference is supported by the distinctly lower S and Cl concentrations of the Mo-related intrusions in the Zhangguangcai-Lesser Xing’an range compared to those associated with porphyry Cu deposits. Our results confirm the fundamental role of increasing compressive stress, crustal thickening, and transcrustal petrologic processes in the formation of porphyry deposits and highlight the role of crustal components in modulating the Mo/Cu ratio of porphyry deposits in a continental arc setting. In addition, we tentatively propose that magmatic apatite SO3 and/or Cl contents can be used to discriminate porphyry Mo-related intrusions from those associated with porphyry Cu deposits in continental arcs.

Timing and Petrogenesis of the Permo-Carboniferous Larvik Plutonic Complex, Oslo Rift, Norway: New Insights from U–Pb, Lu-Hf, and O Isotopes in Zircon

Abstract The Permo-Carboniferous Oslo Rift is a narrow, aborted continental rift with comparatively little extension but voluminous magmatism, developed at a Proterozoic lithospheric discontinuity. On the southern flank of the onshore part of the rift, the multiple intrusive phase Larvik plutonic complex (LPC) encompasses almost the width of the rift with a conspicuous over- to undersaturated assemblage of monzonite (larvikite, tønsbergite) and nepheline syenite (lardalite, foyaite). We present new single-crystal U–Pb–, Lu–Hf–, and O–in–zircon isotope data for the intrusive centers of the complex. U–Pb ages of slightly over- to undersaturated monzonites show a shift in igneous activity from 299 Ma on the eastern edge of the complex to 289 Ma in the west. The western part of the complex is built of two consanguineous magmatic systems with a northward shift in the locus of magmatism, from 296 to 289 Ma (undersaturated monzonite) and 296 to ~288 Ma (nepheline syenite). Moreover, an early oversaturated monzonite was emplaced in the eastern part of the complex at ~302 Ma. The Hf and O isotope composition of the monzonites and nepheline syenites shows little variation with zircon εHf (295 Ma) values of +5.5 to +8 and δ18Ο values of +4.79‰ to +5.49‰. No change in isotope values is observed with varying alkalinity and the Hf and O isotope compositions are interpreted as mantle source values. The source of the precursors of the monzonitic and nepheline syenitic magmas was probably a mildly depleted, sublithospheric peridotite in the lithosphere-asthenosphere boundary zone (at 3–4 GPa). Rhyolite–MELTS modeling implies polybaric fractionation of weakly alkaline basaltic melts from the source having led to an oversaturated/saturated liquid line of descent at ~0.5 GPa and an undersaturated one at ~1 GPa. Magmatism had an imperative role in the evolution of the rift by localizing strain and weakening the lithosphere along the discontinuity. A fractionating crustal melt column was tapped at varying depths in the course of advancing rupturing with no significant mixing of the over- and undersaturated melts, and magma batches were emplaced as sequential nested plutons in response to differential opening of the rift. Among continental saturated-undersaturated alkaline complexes, the LPC may be unique in the emplacement of successive magma batches across a continental rift in the direction of opening, tapping of increasingly deeper parts of the melt column with advancing rifting, and intrusion of two contemporaneous, contrasting magmatic lineages at the end of igneous activity. The marked lithospheric step at the rift locus in the Precambrian basement of southeastern Norway was probably the driving force for the inception and evolution of the Larvik magmatic system.

Revisiting the origin of the Carboniferous Oytag pluton in West Kunlun orogenic belt, northwest China

The Oytag pluton in the West Kunlun orogenic belt of China represents the rare exposed Carboniferous granitoids associated with the early-stage evolution of the Paleo-Tethys Ocean, and its origin is controversial. New zircon U-Pb data suggest an emplacement age of 327–324 Ma for trondhjemite, which slightly predated tonalite (∼314 Ma). These results, along with available published ages, indicate several episodes of magmatism to build the Oytag pluton. Both trondhjemite and tonalite have positive whole-rock Nd isotopic compositions (ε Nd (t) = +4.56 to +8.01 for trondhjemite; ε Nd (t) = +4.68 to +6.15 for tonalite), which positively correlate with whole-rock MgO values. Given the presence of Ordovician inherited zircon grains, variations of whole-rock Nd isotopic compositions are attributed to addition of <10% crustal contaminants during assimilation and fractional crystallization (AFC) processes. Interestingly, trondhjemite and tonalite follow distinct AFC trends, consistent with their different ages. Zircons from two trondhjemite samples have overlapping U-Pb ages but show distinctly different trace elements and internal structures (sector-zoned vs. oscillatory-zoned). The sector-zoned zircons have a less evolved signature of lower Hf contents, higher Eu/Eu* and Th/U ratios than the oscillatory-zoned ones, indicating significant crystal fractionation processes for their parental melts. The more variable Hf contents of the oscillatory-zoned zircons than those of the sector-zoned ones reflect a more dynamic environment where they crystallized. We propose that the sector-zoned zircons were crystallized from locked, interstitial melt in a mush and the oscillatory-zoned zircons from extracted melt. In situ crystal-liquid segregation in a mush chamber thereby could be an efficient approach to the differentiation of Oytag trondhjemite. On the other hand, zircons from tonalite have similar Hf contents to the sector-zoned zircons from trondhjemite, but display slightly higher Eu/Eu* ratios and much higher Th/U ratios. This indicates that tonalite should have experienced much different AFC processes from trondhjemite. Both the episodic magmatism feature and the distinct differentiation processes argue for rises of different batches of magmatism to generate trondhjemite and tonalite in Oytag. Along with the large variations of Hf isotopes (5–6 epsilon units) in zircons, we propose a partial melting of juvenile mafic crust model for the Oytag pluton. Heat from upwelling of asthenosphere may have triggered the partial melting of wet, mafic underplates during the opening of a Paleo-Tethyan back-arc basin. • The Oytag pluton comprises predominant trondhjemite with minor tonalite. • Oytag trondhjemites (327–324 Ma) formed earlier than tonalites (∼314 Ma). • The Oytag pluton was built by different batches of magmas via partial melting.

Generation of magmatism under active continental margins: A thermodynamic study of subduction and translithospheric diapirs

A recent model of continental arc magmatism states that subducted sediments mix physically with the oceanic slab basalts in the subduction channel and produce diapirs that ascend through the mantle wedge, undergo melting and relaminate to the base of the continental crust. From here, melt or magma batches ascend through the crust, producing the Cordilleran batholiths and associated volcanism ranging in composition from 0 to 25 wt% maficity (FeO t + MgO) and 44–79 wt% SiO 2 , with a gap around 11.5 wt% maficity and 57–60 wt% SiO 2 . This model has been predicted by thermomechanical numerical studies and later supported by phase equilibrium experiments; however, thermodynamic modelling to verify whether the composition of the Cordilleran batholiths is reproduced by this model has not yet been carried out to complement the experimental approach. In this article, the evolution of the mélange along the subduction path and during the ascent of the diapirs is investigated, with focus on the conditions of generation of the diapirs, the composition and proportion of melt at the P-T conditions of relamination of the diapirs, the relationship between the basalt:sediment proportion (the composition) of the diapir to this melt composition and proportion and whether the composition of the melts match the Cordilleran batholiths. The key point to test is the process or processes that produce the compositional variability observed in these granitoids. Our findings indicate that (1) parental melts and/or magmas range in maficity from ∼1 to ∼11.5 wt%, with the compositional gap of the Cordilleran trend representing the most mafic composition possible of the parent, (2) the diapirs are restricted to having 50–60% of sediment component: diapirs with lower sediment contents may be too dense to detach from the slab and mélanges with higher sediment contents may be too less dense to be subducted to the depth of formation of diapirs, (3) melts at 1100 °C and 1.5 GPa from these diapirs have the bulk granodioritic composition of the batholiths, (4) restite unmixing occurs when magma batches segregate from the relaminated diapirs, (5) the lower the temperature of the region of the diapir from which the batch segregates, the higher the restite content, (6) the temperature, by controlling the restite content, is the ultimate factor controlling the compositional trends of Cordilleran batholiths, and (7) fractional crystallization or cotectic evolution is still needed to account for the full compositional range of these batholiths, from the parental 1–11.5 wt% to the observed 0–25 wt% in maficity. • Cordilleran granitoids can be formed by melting of subducted basalt and sediment. • Parental magmas of Cordilleran granitoids derive from diapirs relaminated to the crust. • The main process controlling the composition of granitoids is restite entrainment. • The higher the temperature of the diapir, the lower the amount of restites entrained. • Composition of diapirs is restricted to having 50–60% sediment component.

The petrogenesis of chemically zoned, phonolitic, Plinian and sub-Plinian eruptions of Somma-Vesuvius, Italy: Role of accessory phase removal, independently filled magma reservoirs with time, and transition from slightly to highly silica undersaturated magmatic series in an ultrapotassic stratovolcano

A mineralogical, major, LA-ICP-MS trace element mineral chemistry and bulk-rock geochemical study of juvenile samples of the Mercato, Avellino, Pompeii and Pollena eruptions, collected in stratigraphically and volcanologically well-characterized sections of the Somma-Vesuvius stratovolcano (Roman Magmatic Province), along with reference data on the 1944 CE and the Pomici di Base eruptions, highlights the compositional variability of bulk-rock and glass from leucite phonotephrite to garnet-bearing phonolite . The latter products have extreme fractionation of trace elements (e.g., La/Yb n = 126, Zr/Y = 89, Zr/Hf =78, Nb/Ta = 40; Th/U = 2.3), very low Sc, V, Y, HREE and very high As, Tl, Cs, Pb, Th and U. The Pomici di Base products, older than the eruptions described above, range from leucite-bearing shoshonites to trachytes , are devoid of garnet and belong to an independent liquid-line-of-descent, having also different fractionation between trace elements (e.g., La/Yb n = 15; Zr/Y = 12.4, Zr/Hf = 50, Nb/Ta = 15.6; Th/U = 3 in the Pomici di Base trachytes). A marked chemical variability of the observed phases is found. The geochemistry of garnet, amphibole , clinopyroxene and other phases shows wide variations of concentrations and elemental ratios (e.g., La/Yb n up to 520 in the sadanagaite coexisting with garnet). The magmatic evolution is dominated by low-pressure, oxidized, nearly closed-system fractional crystallization of clinopyroxene, plagioclase , leucite, ±magnetite, ±biotite, ±olivine and apatite in the transition tephrite-phonotephrite, and of potassic sanidine (±hyalophane), Fe-clinopyroxene, melanite garnet±Fe-amphibole in the transition tephriphonolite-phonolite. Mineralogic and geochemical evidence and modelling points out the existence of independent, zoned magma batches throughout the activity of the stratovolcano, which possibly started to crystallize at similar depths. The Somma-Vesuvius magmas thus evolved in shallow independent reservoirs with respect to those of the neighbouring volcanic complexes (Campi Flegrei, Ischia) as shown, for instance, by the contrasting compositional trends of clinopyroxene and amphibole, and have very limited evidence of crustal contamination (and/or carbonate assimilation). • Pollena, Pompeii, Avellino and Mercato are zoned eruptions with phonolites • The magmas which formed the various eruptions evolved independently • Different garnets dominate HREE, Y and HFSE; amphibole is strongly HREE-poor • The phonolites and their mineral compositions have no tendency to become peralkaline • F, Cl and H 2 O dominate the volatile budget

Assembly and development of large active calderas hosting geothermal systems: Insights from Los Humeros volcanic complex (Mexico)

Caldera volcanoes are complex geological systems that show, during and after their formation, a wide variability in terms of eruptive styles, magmatic and geochemical evolution, and volcanic structures. Los Humeros Volcanic Complex (LHVC) is a key area where to study these factors. It is located in the easternmost sector of the Trans-Mexican Volcanic Belt and hosts an active geothermal system. LHVC exemplifies the complex nature of calderas, recording periods of alternated effusive and explosive eruptive phases, a heterogenous magmatic source, diverse basaltic to rhyolitic products, and the overlapping of caldera collapse events. This study provides an up-to-date interpretation of the caldera framework. Our work is based on a detailed analysis of literature and novel data aimed at the morpho-structural caldera configuration, and supported by geophysical modeling, and petrology of the volcanic products. The new results confirm a more complex evolution of LHVC involving geometric configurations related to multiple caldera collapse events producing both asymmetric trapdoor, and piecemeal styles. The present configuration of the caldera framework involved the initial formation (at 164 ka) of a large caldera (Los Humeros) of ca. 15–17 km diam (shorter than previously reported), which was overlapped asymmetrically by the younger (69 ka) and smaller (ca. 8.5–10 km diam) Los Potreros caldera over the west side, as a result of multiple sequential collapsing events. The final configuration of the caldera rims (scarps) was promoted by incremental growth. A post-caldera resurgence phase was induced by the injection of multiple shallow intrusions of silicic bodies, causing a more localized deformation of the caldera floor. The Holocene volcanism records the recurrent injection of compositionally distinct magma batches uprising from different depths of the LHVC transcrustal magmatic plumbing system. This latter indicates the transition from the caldera stage magmatic system dominated by a single, large, and shallow magmatic body to a more complex, polybaric post-caldera stage plumbing system, made up of a lower crust mafic reservoir feeding smaller magma batches vertically distributed in the whole crust. The integration of all the existing data constitutes the archive to better understand how large active calderas hosting exploitable geothermal systems assemble through time.

Inferences on the magmatic plumbing system at Stromboli volcano (Italy) from trace element geochemistry of matrix glasses and minerals in different types of explosive eruptions

Stromboli (Italy) is a basaltic volcano characterized by persistent, mild strombolian activity, occasionally interrupted by lava effusion and more violent explosive events, named major explosions and paroxysms depending on their intensity and magnitude. The normal activity is fed by a shallow and degassed highly porphyritic (HP) shoshonitic basalt carrying about 50 vol.% crystals settled in a shoshonitic glassy matrix (K2O > 3.8 wt.%). The more energetic explosions erupt a deep, volatile-rich, low-porphyritic (LP) magma with < 10 vol.% crystals in a shoshonitic basaltic glassy matrix (K2O < 2.4 wt.%). Products with intermediate glass composition are also found in the more violent explosive events. In this study, we present a new data set of major and trace element contents in matrix glasses and minerals performed in products from different types of explosive activity that occurred at Stromboli between 1998 and 2020. This large data set is used to put constraints on the evolution and architecture of the intermediate plumbing system, where the transformation from LP to HP occurs. Results indicate that, compared to paroxysms, the glassy matrices of the LP pumices from major explosions are richer in incompatible trace elements (and K2O wt.%) due to < 15 wt.% fractionation of clinopyroxene and olivine. This points to a chemical zoning of the deep reservoir and suggests that major explosions are fed by magmas residing in its upper part. Among the major explosions, the homogeneous intermediate glasses in the products from the 19 July 2020 event originate from the interplay of mixing and crystal fractionation processes. The crystallization of euhedral microphenocrysts of An-rich plagioclase suggests that batches of magma can pond and crystallize for few days (< 11) at the base of the intermediate zone of the plumbing system, at pressure coinciding with the entering of plagioclase into the system (< 100 MPa). As a relevant point for understanding the pre- and syn-eruptive magma dynamics, data indicate a positive correlation between the magnitude of the explosions and the depth of the supply magma.

Rare-Metal Enrichment and Nb–Ta Fractionation during Magmatic–Hydrothermal Processes in Rare-Metal Granites: Evidence from Zoned Micas from the Yashan Pluton, South China

Abstract Rare-metal mineralization in aluminous granites (i.e. rare-metal granites, RMGs) is usually associated with high degree of magma fractionation, strong enrichment in fluxing elements (Li, B, P, and F), and the operation of both magmatic and hydrothermal processes. Experimental data have established the theoretical basis for melt- or fluid-driven rare-metal mineralization. To investigate these mechanisms in natural samples, this paper presents detailed textural and compositional variations for zoned micas and rare-metal minerals from rocks of the Late Jurassic Yashan rare-metal granitic pluton, South China. This pluton preserves an evolutionary sequence from ore-barren rocks to those containing ores formed by Ta–Nb–Li mineralization during magmatic–hydrothermal processes. Three main units of the Yashan pluton are exposed: from bottom to top, these comprise protolithionite–muscovite granite (Unit I), Li-mica granite (Unit II), and topaz–lepidolite granite (Unit III), representing crystallization from successive magma batches from a deep-level magma chamber. The gradual decrease in Nb/Ta and K/Rb in both whole rocks and micas from units I to II to III illustrates fractionation within the magma chamber, accompanied by successive enrichment of incompatible elements in the residual melt. Furthermore, the strong enrichment of fluxing elements such as Li, P, and F in Unit III likely led to lower magma viscosity, which may have helped to expulse interstitial residual melt from the magma chamber. Zoned micas from the Yashan pluton are composed of Li-phengite or lepidolite inner cores and muscovite rims. Cores show increasing Li, F, and rare-metal elements from units I and II to Unit III, consistent with the trend of magmatic fractionation. In contrast, rims that formed in the residual hydrosilicate melt show marked decrease in F and rare metals (i.e. Li, Cs, Rb, Nb, Ta, Sn, and W). Chemical variations of the zoned micas and rare-metal minerals were used together with Rayleigh fractionation modeling calculations to reconstruct the contribution of melt–fluid immiscibility to the crystallization of rare metals, which took place via the combination of melts enriched in high-field-strength elements (e.g. Nb and Ta) and fluids enriched in mobile elements (e.g. Mn, Fe, W, and Cs). Strong fractionation of Nb and Ta, along with extreme enrichment of Ta, in RMGs is caused mainly by the crystallization of micas and columbite-group minerals, and also depends on the degree of fractional crystallization. The subsolidus alteration of micas by acidic fluids may have leached some Nb rather than Ta during chloritization, which could have partly contributed to the enrichment of rare metals and fractionation of Nb–Ta. It is concluded that magmatic–hydrothermal processes, including magmatic evolution and fluid exsolution, are critical for rare-metal enrichment and Nb–Ta fractionation in RMGs.

Genesis of the Wilsons Promontory batholith of southeast Victoria, Australia: heterogeneities in the Proterozoic basement

The monzogranites and syenogranites of the Early Devonian (Emsian), post-tectonic, high-level, S-type Wilsons Promontory batholith contain magmatic garnet and cordierite, in addition to biotite, as well as biotite–quartz pseudomorphs after early orthopyroxene. The batholith forms the southern tip of the Australian mainland and consists of four penecontemporaneous plutons, emplaced at shallow crustal depth, intruding Ordovician to Lower Devonian lower greenschist-facies, quartz-rich metaflysch rocks in the Bassian tectonic zone. The magmas were initially H2O-undersaturated and at T > 800 °C, having been generated by partial melting of greywacke-dominated, Paleoproterozoic to Mesoproterozoic (ca 1.6 Ga) basement rocks in the Selwyn Block (the northern tip of the VanDieland microcontinent). However, the specific source rocks differ from those of the S-type granitic plutons in the neighbouring Melbourne tectonic zone. This probably reflects compositional variations within the Selwyn Block, which has been proposed to form the basement terrane beneath both the Melbourne and Bassian zones. The compositional zonation in the batholith is due to sequential intrusion of separate magma batches rather than in situ differentiation, and the main mechanism responsible for primary chemical variations in the rocks appears to have been peritectic assemblage entrainment. Some local differentiation also occurred, probably by filter pressing, to produce aplites and tourmaline-bearing pegmatites. Normal granitic rocks (i.e. not aplites or pegmatites) with SiO2 >76 wt% were not produced through differentiation of other granitic magmas in the batholith. At the global scale, this suggests that fractionation models for the origins of high-silica granitic and rhyolitic magmas may need revision. KEY POINTS S-type magmas of the Wilsons Promontory batholith were H2O-undersaturated and at T > 800 °C. They were generated by partial melting of greywacke-dominated, Paleoproterozoic to Mesoproterozoic (ca 1.6 Ga) basement. The sources differ from those of the plutons in the Melbourne Zone, suggesting heterogeneity in the Selwyn Block. Most compositional variation is due to intrusion of separate magma batches and peritectic assemblage entrainment, not differentiation.

Petrology, geochemistry and new U–Pb ages of the Epembe syenites and carbonatite, northwest Namibia

The Epembe Alkaline Carbonatite Complex (EACC) in northwestern Namibia, was emplaced along a fault zone into medium-to high-grade Palaeoproterozoic basement rocks of the Epupa Metamorphic Complex (EMC) and extends over a distance of 9 km in a south-easterly direction with a width of 1 km. Alkaline rocks constitute the main lithologies and are cross-cut by a calcite-carbonatite dyke. The alkaline rocks can be classified as syenite and nepheline-bearing syenites, with alkalic, metaluminous and ferroan affinities. An emplacement age of 1220 ± 3 Ma (2 SE, MSWD = 1.3) is determined for the syenite, whereas two nepheline syenite samples give similar magmatic ages of 1209 ± 3 Ma (2 SE, MSWD = 1.1) and 1205 ± 13 Ma (2 SE, MSWD = 2). The nepheline syenite ages correspond with the 1198 ± 5 Ma (2 SE, MSWD = 1.1) age of the carbonatite interpreted as magmatic. The nepheline-normative syenites define broadly linear trends in Harker plots consistent with evolution by fractional crystallization involving pyroxene and apatite. The rare-earth element (REE) pattern of the syenite shows a negative Eu anomaly, consistent with prior plagioclase formation and fractionation. The syenite is characterized by the absence of a negative Nb anomaly on a Primitive Mantle (PRIMA)-normalized diagram with Ce/Pb and Nb/U ratios of 12 and 19, respectively, suggesting that it was not affected by crustal contamination during ascent and emplacement. The absence of a crustal signature in the evolution of the syenite indicates that the syenite is not related to the nepheline syenite by combined assimilation and fractional crystallization and is likely to have been emplaced as a distinct magma batch derived from the mantle. REE patterns of the nepheline syenite and carbonatite are coherent showing LREE enrichment relative to HREE. On a PRIMA-normalized multi-element diagram both the carbonatite and syenite display depletions in Zr, Hf and Pb with positive Sr. Low Mg number and Cr and Ni contents of the carbonatite and nepheline syenite suggest that these rocks crystallized from an evolved parental magma rather than a primary mantle derived-melt. Geochemical similarities between the carbonatite and nepheline syenite also suggest that the origin of the Epembe carbonatite is not a product of immiscible separation between a carbonate and silicate melt. Instead, the carbonate may represent a melt fraction that formed after partial crystallization of the nepheline syenite. Magma assimilation of the country rocks by the carbonatite magma is evident from the occurrence of biotite, feldspars and syenitic xenoliths in the carbonatite.

Homogenization of zircon Hf isotopes during late-stage granite crystallization

Zircon Hf isotope compositions provide valuable constraints on granite petrogenesis and the evolution of continental crust. However, the causes of variations in zircon Hf isotope ratios during magma evolution are poorly understood. Here we compiled zircon UPb ages and LuHf isotope compositions of the Permian-Triassic granites in South China to demonstrate that Hf isotope variability is common in the early stages of granitic magma evolution, but Hf isotope ratios are homogenized during the final stages of granite crystallization. Late-grown zircons have relatively high concentrations of U and have εHf(t) values approximating the median values of multiple analyses for the low-U zircon in the same samples. Hafnium isotope variations measured in low-U zircon could be related to incremental assembly of multiple isotopically distinct magma batches derived from partial melting of heterogeneous crustal rocks; these discrete magma batches may have remained isolated from each other during segregation and extraction in the source and during ascent through the crust. After emplacement, these magma batches may have experienced variable degrees of mixing resulting in homogenization of Hf isotopes in late crystallized zircon. Coupled with whole-rock geochemical data of the Permian-Triassic granites, we argue that the Hf isotope homogenization is related to crystallization of zircon from late-stage magmas or precipitation of hydrothermal zircon from exsolved hydrothermal fluids during the final stages of magma solidification. Hafnium isotope homogenization in late-stage melts/fluids may have been driven by an increased concentration of fluorine and the efficient diffusion of hafnium-fluoride complexes through silicic melt. These results provide a deeper understanding of the mechanisms that contribute to Hf isotope heterogeneity in zircon during granite crystallization.