Melt Extraction Research Articles

The Role of Latent Heat Buffering in the Generation of High-Silica Rhyolites

Abstract The physical process of crystal-melt separation is responsible for the accumulation of small to very large volumes (>100 km3) of eruptible rhyolitic melt in the shallow crust. Granitic intrusions, although providing a terminal, time-integrated image of melt segregation processes, host an unmatched record of the physical properties controlling mechanisms and rates of interstitial melt extraction from a crystal-rich source. We applied mass balance calculations and thermodynamic modeling simulations to an extensive bulk rock geochemistry dataset (>150 samples) collected in a Permian upper-crustal granitoid intrusion of the Italian Southern Alps. Textural and geochemical evidence indicate that this intrusion constituted a single, zoned magma body, with a crystal-rich base and a thick (~2 km), high-silica cap (75–77 wt% SiO₂). The large compositional variability of the crystal-rich materials suggests variable degrees of melt extraction efficiency and corresponding terminal porosities. Specifically, the loosely bimodal distribution of porosity values (φ) indicates that at least two distinct melt segregation mechanisms were operating in this system, which produced both high (0.65–0.45) and low terminal porosities (0.45–0.25) in the crystal-rich, cumulate materials. Modeling of latent heat budget shows that coexistence of cumulate products with differing terminal porosity signature can be explained by melt segregation processes taking place at different depths across a thick, interconnected magmatic reservoir with an initial homogenous water content (~4 wt% H2O). Deep in the mush column, low water activities (aH₂O < 0.5) promoted thermal buffering of cooling magma at high crystallinities, enabling residual melt extraction by percolation through a crystalline framework accompanied by compaction. Instead, at shallower depths, high water activities (aH₂O > 0.5) ensured prolonged magma residence at porosities that promoted crystal melt separation via hindered settling. Distinct melt extraction processes, acting synchronously but at different depths in vertically extensive silicic mush columns, can account for the large volumes of residual, haplogranitic melt mobilized during the relatively short lifespan of upper crustal magma reservoirs (~105 years).

Nature and evolution of lithospheric mantle beneath the Northeast China: Constraints from the Houtuancun mantle xenoliths

Mantle xenoliths entrained in a newly discovered Cenozoic basalt at Houtuancun within the Yilan-Yitong fault zone, situated beneath the Xing'an Mongolian Orogenic Belt in Northeast China, provide a unique opportunity to constrain the nature and evolution of the lithospheric mantle beneath this region. The mantle xenoliths are predominantly spinel lherzolite, with rare occurrences of harzburgite, wehrlite, and websterite. The Houtuancun spinel lherzolite and harzburgite are characterized by high forsterite contents in olivines (89.3–91.1), and their clinopyroxene rare earth element patterns range from slightly light rare earth element (LREE)-depleted to LREE-enriched, reflecting variable degrees of partial melt extraction (up to 20%) overprinted by later metasomatism. Despite this, the Fe isotopic composition of the Houtuancun lherzolite and harzburgite exhibits little variation (δ57/54Fe = −0.05–0.09‰), with an average of 0.02 ± 0.09‰ (2SD, n = 9), close to the suggested δ57/54Fe value of 0.04 ± 0.04‰ for the upper mantle, and there is no significant inter-mineral disequilibria. In contrast, the wehrlite has a slightly lower Mg# [=100 ╳ Mg2+/(Mg2+ + Fe2+)] (88.4) and δ57/54Fe (−0.05‰) compared to the spinel lherzolite and harzburgite, reflecting kinetic iron fractionation during melt-rock interaction. The websterite exhibits the lowest MgO (20.6–24.8 wt%), highest Al2O3 content (5.8–7.3 wt%), and the highest δ57/54Fe (0.25–0.28‰) among the studied xenoliths, suggesting an episode of subduction-related silicate melt circulation through the upper mantle. The results demonstrate compositional heterogeneity of the mantle beneath the Xing'an Mongolian Orogenic Belt in Northeast China. The observed heterogeneity may result from a combination of varying degrees of melt extraction of the subcontinental lithospheric mantle and metasomatism by melts related to the subduction of the (Paleo)-Pacific plate.

Magmatic to solid-state evolution of a shallow emplaced agpaitic tinguaite (the Suc de Sara dyke, Velay volcanic province, France): implications for peralkaline melt segregation and extraction in ascending magmas

Abstract. In the last decades the mush model has been generalized to the complete trans-crustal magmatic system in which differentiation would be driven by segregation and extraction of trapped melts from crystal-rich mushes. Melt extraction processes involved are porous flow and strain localization, the latter being regarded as the main process acting during transfer through dykes and necks along which high differential stresses are acting on. We combine structural measurements together with petrological analyses and textural observations to constrain the model of emplacement and finally emphasize how shear deformation and strain localization structures promoted the residual melt segregation that occurred in a shallow silica-undersaturated peralkaline intrusion (Suc de Sara, Velay volcanic province, French Massif Central). In this study, we demonstrate that segregation and subsequent extraction of the CO2-rich residual melt occurred during magma ascent and final emplacement of the Suc de Sara tinguaite. Contrasting features of shear deformation between the margins that exhibited different permeabilities highlight that melt segregation started by compaction as a loose packing of emerging microlites and continued with melt filling of an anastomosed C/C′ band network developing in the crystal-rich mush subjected to high shear strain. Subsequent melt extraction throughout the country rock was controlled by the permeability of the hanging wall. Along the western hanging wall of the intrusion, extraction of the residual melt was prevented by the 15 cm thick chilled margin. In contrast, segregated melt circulated through the highly porous and permeable eastern margin, causing the fenitization of the country rock.

Granitic mush compaction by impingement of a large felsic enclave: Evidence from the Aztec Wash pluton (Nevada, USA)

Abstract Compaction and melt extraction are key processes in the generation of melt-dominated magma bodies. We analyze magmatic fabrics in the Aztec Wash pluton (15.7 Ma; Nevada, USA) to understand magma dynamics during emplacement of a large (2 m diameter) felsic enclave (LFE) in a silicic magma body. We hypothesize that the LFE moved through crystal-poor magma, settling into crystal-rich magma. We analyze four granite samples collected near the enclave: two 0.25 m below the enclave (E-5 and E-2); one 0.5 m below (E-6); and one 2.5 m to the side (E-3, far field). We use the scanning electron microscope (SEM) to gather data using backscattered electron (BSE) microscopy, energy dispersive spectroscopy (EDS), and electron backscatter diffraction (EBSD) imaging. BSE maps show large (up to 10 mm) euhedral alkali feldspar grains in E-5 and E-2, while E-6 and E-3 have large (up to 6 mm) alkali feldspar with distinct irregular and wavy rims. Dark BSE alkali feldspar rims correlate with Ba depletion in EDS maps. Broken feldspars in E-5 and E-2 have mismatched zones in BSE and orientation differences in EBSD. Grain orientation spread (GOS) maps show that alkali feldspar has the highest degree of intragrain deformation, followed by plagioclase, then quartz. Histograms of the angle between the vertical plane and the pole to the (010) plane demonstrate clustering in E-5, suggesting an enhanced fabric. Differences in rim texture suggest melt extraction beneath the enclave and melt retained in the far field. EBSD data reveal an enhanced foliation beneath the LFE. We conclude that mush under the LFE was deformed during emplacement of the LFE, expelling melt and reorienting grains.

The Weixi High‐silica Granitoids in the Central Sanjiang Orogenic Belt, Southwest China: Implications for Growth of the Continental Crust

AbstractHigh‐silica granitoids record the formation and evolution of the continental crust. A new intrusive complex has been recognized among silicic volcanic rocks of the Weixi arc, Southwest China. The intrusions consist of granites, granitic porphyries, and granodiorites. Zircon U‐Pb age data indicate that the Weixi granitoids formed at 248–240 Ma and were coeval with silicic volcanic rocks of the Weixi arc. The Weixi granitoids are enriched in Rb, Th, and U, depleted in Ba, Sr, Nb, Ta, and Ti, and have high light/heavy rare earth element ratios and slightly negative Eu anomalies. The Weixi granitoids have negative εNd(t) values (–9.8 to –7.8) and negative zircon εHf(t) values (–12.02 to –5.11). The geochemical and isotopic features suggest the Weixi granitoids were derived by partial melting of ancient crustal material. The Weixi granitoids and silicic volcanic rocks were derived from the same magma by crystal accumulation and melt extraction, respectively, and they record the formation of a continental arc in the central Sanjiang orogenic belt.

Chemical constraints and tectonic setting of the Ginebra ophiolite complex

The occurrence and types of ophiolites in the Central Cordillera of the Colombian Andes have not been identified in terms of their tectonic settings and geochemical features. Therefore, the study of the tectonic setting and subduction zones in the ophiolites’ magmatic evolution becomes a significant issue. The Ginebra Ophiolitic Complex (GOC) is located on the western flank of the Central Cordillera in the Valle del Cauca Department, southwest of Colombia, assigned to the Lower Cretaceous period, as indicated by the contact of the GOC with the Buga Batholith. The GOC is formed by three main lithologies: amphibolite, gabbro, and ultramafic rocks (pyroxenite and peridotite), according to Ossa (2007). The rocks in the Puente Piedra section, located in the northeast of Ginebra municipality, are part of the gabbroic cumulates within the ophiolitic sequence. This sequence comprises interspersed layers of gabbroic cumulates, gabbrodiorites, and diorites; interpreted as a result of crystal accumulation and fractional crystallization processes. The mafic rocks of the GOC are sub-alkaline and correspond to the low potassium (K2O wt% 0.03-0.06) tholeiitic series. Geochemically, they have SiO2 wt% ranging from 49 to 61 and an aluminum oxide saturation index of approximately ~0.4 and ~0.8, indicating a metaluminous type. The geochemistry of the studied rocks from Puente Piedra section indicates that the GOC formed through fractional crystallization and accumulation processes from a single magma source. Cluster analysis, used to compare the geochemistry of GOC rocks and the Amaime Complex basalts, suggests a similar magma source, possibly linked to multiple recharge events that underwent fractional crystallization, melt extraction, and accumulation processes. The geochemical parameters are indicative of a suprasubduction zone ophiolite, characterized by a low potassium tholeiitic series affinity, TiO2 values typically <1.2 wt%, Th enrichment typical of subduction zones, high Pb content, and low values of Ti, Y, Yb, Ta, Nb, Zr, and Hf.

Identification of variable contributions to arc rocks by Li[sbnd]Mo and Sr–Nd–Hf–O isotope compositions: A study of the Mafan mafic rocks in the Hong'an orogen, Central China

The geochemistry of arc rocks provides a window to probe compositional variability in mantle wedge and slab-derived inputs. However, quantitative estimation of their contribution is a great challenge. Here we report an integrated study for determination of variable contributions to early Paleozoic arc plutons (the Mafan diorites and gabbros), from the northern Hong'an orogen, central China. Zircon grains from the diorites yielded crystallization age of 445 ± 3 Ma, whereas the gabbros have zircon UPb age of 432 ± 3 Ma. Although both the two mafic rocks suites exhibit typical arc-like trace element distribution patterns, they show a series of differences in stable Li–Mo–O and radiogenic Sr–Nd–Hf isotope compositions, as well as major and trace element features. Geochemically, the diorites are medium-K calc-alkaline series with high and positive whole-rock εNd(t) (+3.8 to +4.5) and zircon εHf(t) values (+12.4 to +14.3), and normal mantle-like zircon δ18O values (5.1 ± 0.2 ‰). Contrarily, the gabbros are tholeiitic and have relatively low whole-rock εNd(t) (+ 2.3 to +3.1) and zircon εHf(t) (+5.0 to +8.5) values with normal mantle-like zircon δ18O values (5.4 ± 0.3 ‰). The diorites have δ98Mo values of −0.15 to −0.03 ‰, while the gabbros have δ98Mo values of −0.28 to −0.12 ‰. These values are slightly higher than those of the mid-ocean ridge basalts (MORB) (−0.20 ± 0.01 ‰). Whereas, δ7Li values of the studied diorites (2.15–4.49 ‰) and gabbros (2.89–5.33 ‰) are mostly consistent within the range of MORB. In combination with the above results, the primary magma of diorites was inferred to be derived from 1 to 5% partial melting of depleted lithospheric mantle metasomatized by slab-derived fluid with minor sediment-derived melts. While the gabbros were suggested to be generated by 5–10% partial melting of more depleted lithospheric mantle that had undergone perhaps extra 0.3% melt extraction, which had further been metasomatized by more subducted materials. In summary, both the slab-derived components and nature of the mantle wedge sources are quantitatively constrained for the mafic rocks in the Mafan plutons. Our new data suggest that the Mafan plutons exhibited the Mariana arc-affinity and represented products of an intra-oceanic subduction system in the northern Hong'an orogen. Combined with previously published data for the early Paleozoic mafic magmatic rocks along the Qinling-Tongbai-Hong'an orogenic belt, there exists an oceanic arc in the east (Hong'an orogen) and an Andean-type continental arc in the west (Qinling-Tongbai orogen).

Peritectic orthopyroxene entrainment during partial melting of garnet peridotite produced the Bushveld Complex chromite deposits

AbstractOne of the largest chromium deposits on Earth occurs in the Rustenburg Layered Suite (RLS) of the Bushveld Complex as laterally continuous chromitite layers. None of the hypotheses proposed for the origin of the chromitites can explain both the abundance of Cr in the RLS and the unusual enrichment in Cr and V over Ni, relative to typical depleted mantle values. This study investigates the possibility that the layering and chromitite formation are consequences of the entrainment of source components into the magmas that formed the RLS. Thermodynamic modelling results reveal a wedge-shaped domain in pressure-temperature space in the subcratonic mantle within which Cr-bearing orthopyroxene forms as a peritectic product of incongruent melting. Entrainment of this orthopyroxene produces magmas that crystallise peritectic olivine and chromite on ascent, due to the consumption of orthopyroxene by melt. The chromite- and olivine-bearing magmas intrude as sills and can produce chromite and dunite layers by density separation. This model, which interprets the RLS Sr-isotopic composition to reflect prior mantle metasomatism by crustal fluids (ideally ancient and of low volume), readily explains the formation of chromitite layers from relatively thin sills, as well as the very high ratios of Cr and V to other compatible elements relative to typical mantle compositions. The special circumstances required to produce the RLS chromitites do not relate to some oddity of repetitive crustal assimilation or magma compositions that allow chromite-only saturation. Rather, they relate to speed of melting and magma extraction which enabled peritectic orthopyroxene entrainment to the magmas.

Calcium isotopic composition of the bulk silicate Earth: A komatiite perspective

Calcium (Ca) isotopic composition of the mantle has been previously estimated as δ44/40Ca. However, δ44/40Ca of the mantle records both its stable isotopic composition and excess/deficit in radiogenic 40Ca that can be produced by the decay of 40K. Here, we present the high-precision Ca isotopic compositions of 16 Archean komatiites from five representative localities (Barberton, Munro, Abitibi, Yilgarn, and Taishan). Due to the high degrees of partial melting that led to their formation, these komatiites are used here to simultaneously estimate ε40/44Ca, δ44/40Ca, and δ44/42Ca of the bulk silicate Earth. The komatiites exhibit a resolvable negative ε40/44Ca relative to SRM915a, with an average of −0.6 ± 0.1 (2SE, 95 c.i.% if not specified). Given the lack of significant differences among the samples and the non-zero value of ε40/44Ca, we suggest using δ44/42Ca to report the stable Ca isotopic composition in this study. Excluding one sample (SD6/400), which may have been affected by alteration, the remaining komatiites display relatively homogeneous δ44/42Ca, ranging from 0.39 ± 0.02‰ (2SE, N = 16) to 0.49 ± 0.02‰ (2SE, N = 10), defining a mean of 0.44 ± 0.02‰ (2SE, N = 15). The absence of correlations between δ44/42Ca and LOI, Sr/Nd, (Dy/Yb)N, and (La/Sm)N among the komatiites indicates that post-eruption alteration, the influence of residual garnet, as well as prior melt extraction from the mantle source have had a negligible influence on their stable Ca isotopic compositions. Based on available experimental petrological data and isotope fractionation factors, it is demonstrated that komatiites were insignificantly fractionated from their mantle sources (i.e., <0.02% in δ44/42Ca). Their δ44/42Ca represents the ambient mantle at the time of segregation. Accordingly, we suggest that the average δ44/42Ca and ε40/44Ca of the komatiites measured here are representative of the bulk silicate Earth. Our new estimate can serve as a reference point for Ca isotopes as a two-dimensional tracer for sources and geological processes.

Reconciling Mars InSight Results, Geoid, and Melt Evolution With 3D Spherical Models of Convection

AbstractWe investigate the geodynamic and melting history of Mars using 3D spherical shell models of mantle convection, constrained by the recent InSight mission results. The Martian mantle must have produced sufficient melt to emplace the Tharsis rise by the end of the Noachian–requiring on the order of 1–3 × 109 km3 of melt after accounting for limited (∼10%) melt extraction. Thereafter, melting declined; however, abundant evidence for limited geologically recent volcanism necessitates some present‐day melt even in the cool mantle inferred from InSight data. We test models with two mantle activation energies and a range of crustal Heat Producing Element (HPE) enrichment factors and initial core‐mantle boundary temperatures. We also test the effect of including a hemispheric (spherical harmonic degree‐1) step in lithospheric thickness to model the Martian dichotomy. We find that a higher activation energy (350 kJ mol−1) rheology produces present‐day geotherms consistent with InSight results, and crustal HPE enrichment factors of 5–10‐times produce localized melting near or up to present‐day. The 10‐times crustal HPE enrichment is consistent with both InSight and geochemical results and also produces present‐day geoid power spectra consistent with Mars. However, calculations that match the present‐day geoid power spectra require more than 60% melt extraction to produce the Tharsis swell. The addition of a degree‐1 hemispheric dichotomy, as an equatorial step in lithospheric thickness, does not significantly improve upon melt production or the geoid.

In-situ trace element and Sr isotope in mantle xenoliths from the West Qinling orogen: Evidence for carbonate metasomatism

Subduction process plays an important role in triggering chemical metasomatism of the lithospheric mantle beneath orogenic belts. The lithospheric mantle beneath the West Qinling orogen is extensively heterogeneous, but the lithospheric modification process and related metasomatic agents remain poorly constrained. In this paper, we report in-situ elemental and isotopic compositions of minerals in mantle xenoliths from the West Qinling orogen. The mantle xenoliths are lherzolites and can be categorized into two types. Lherzolite with refractory compositions represents the residue of lithospheric mantle that have undergone high degrees of melt extraction. In contrast, lherzolites with fertile compositions contain clinopyroxenes variably enriched in incompatible elements and radiogenic Sr isotopes. These fertile lherzolites represent lithospheric mantle modified by metasomatic process in different degrees. High (La/Yb)N, Ca/Al, Nb/Ta, Nd/Yb, 87Sr/86Sr ratios and low Ti/Eu ratios of the clinopyroxenes suggest the metasomatic agent was dominated by carbonate-rich melts. This inference is supported by the identification of carbonate inclusion, reaction structures, and clinopyroxene chemical zoning in the lherzolites. Considering the tectonic background of the West Qinling orogen, the orogenic lithospheric mantle would have been metasomatized by carbonate-rich melts derived from partial melting of subducted carbonate-bearing Paleo-Tethyan oceanic slab. In this regard, petrological observations, along with in-situ elemental and Sr isotopic compositions of clinopyroxenes in mantle xenoliths can be powerful means to identify carbonate metasomatism and constrain the evolution of the lithospheric mantle beneath orogenic belts.

Mineral chemistry and melt evolution of the mantle wedge peridotites in the Late Cretaceous Zagros Belt ophiolites (Iran): clues for the subduction initiation-induced forearc magmatism

We investigate in this paper mineral compositions and geochemical evolution of the mantle wedge peridotites preserved in the Late Cretaceous Zagros ophiolites of SW Iran. Mantle peridotites above subduction zones commonly experience distinct melting, depletion and refertilization processes as a result of the circulation of fluids derived from subducting slabs and flux melting. Our results reveal that the mantle wedge peridotites in the Zagros ophiolites are characterized mainly by residual and impregnated types. Residual peridotites resulted from early depletion and later refertilization processes, whereas impregnated peridotites developed due to episodic melt impregnations within and across the mantle. Mg#s and NiO contents, spinel Cr#, Mg#, and TiO 2 in olivines, Mg# and Al 2 O 3 contents of orthopyroxenes, and Mg#, TiO 2 and Al 2 O 3 contents in clinopyroxenes of dunites, harzburgites and lherzolites indicate the significant role of re-equilibration processes among different mineral phases and interactions with basaltic melts percolating within the host peridotites. The observed geochemical variations in the mineral chemistry of the Zagros peridotites reflect changes in magma chemistry and fluctuations in the degree of melt extraction and melt–rock interactions within the mantle peridotites. However, our data suggest that Mg–Fe distribution in the spinels of some dunites and harzburgites might also have resulted from subsolidus redistribution and exchange with surrounding olivine grains. Spinel and clinopyroxene phases in gabbroic rocks and ultramafic cumulates within the Zagros ophiolites also show significant variations in their compositions, suggesting that their magmas evolved from MORB-like to IAT, calc–alkaline and boninite suites, typical of subduction initiation-generated melts. Hence, the Zagros ophiolites present a case study of time-progressive melt evolution of the forearc oceanic lithosphere. Supplementary material: Datasets of the mineral chemistry and melt evolution of the mantle wedge peridotites in the Late Cretaceous Zagros Belt ophiolites (Iran) are available in tabular and figurative form at https://doi.org/10.6084/m9.figshare.c.7093528 . Thematic collection: This article is part of the Ophiolites, melanges and blueschists collection available at: https://www.lyellcollection.org/topic/collections/ophiolites-melanges-and-blueschists

Multi-stage ultrahigh temperature metamorphism in the lower crust of the Kaapvaal craton recorded by U–Pb ages of garnet

U–Pb ages were determined by split-stream LA-SF/MC-ICPMS in garnets from UHT granulite xenoliths (Star mine, South Africa; 124 Ma). They give a considerable age range of 400 million years with well-defined maximas at 3.09, 3.01 and 2.75 Ga. The oldest peak overlaps with the changeover from tonalites to K-granites at 3.14–3.04 Ga and with zircon ages of the mid-crustal granulites of the Vredefort dome (3.1 Ga) in the wake of the 3.2 Ga collision of three terrains that compose the Witwatersrand block. Subduction (or sagduction) of the uppermost crust in an ultrahot orogen setting brought shales and greenstones to the lower crust. Ultrahigh temperature (UHT) conditions are the result of high mantle potential temperatures and self- heating by the radioactive inventory of the subducted lithologies. Metamorphism, anatexis to very high degrees and melt extraction left UHT granulites as residue. Rejuvenation of UHT conditions was brought about by Dominion Group magmatism between 3.0 and 2.95 Ga. Magmatic uprise caused intense shearing in the lower crust followed by recrystallisation of the shear zones to generate the younger garnet age group. Ventersdorp flood basalt volcanism caused similar processes at around 2.72 Ga and generated the third garnet age group. Zircon gives U–Pb ages mainly around 2.72 Ga (both literature and our own data) i.e. zircon adjusted or newly crystallized at the youngest UHT event. Only few zircon grains retained older ages up to 2.94 Ga. Still unconstrained, but very high closure temperatures (≥ 1100 °C) for the U–Pb system in garnet keep the memory of the oldest ages in UHT granulites. Such ages can only be reset by recrystallization. This way, garnet records a prolonged high-temperature history of the lower crust of the Kaapvaal craton.

An enhanced chemical maturation of juvenile arc crust recorded by the Urad Houqi intrusive suite in Langshan arc, Central Asia

Continental arc magmatism is crucial in producing juvenile continental crust of andesitic to dacitic composition, yet its impact on the formation of the modern mature continental crust, which is more enriched in potassic and granitic components, remains poorly understood. We examined a suite comprised of norite gabbro, gabbroic diorite, high-Mg diorite, dioritic enclave, and medium-K to high-K granites (the Urad Houqi suite) from the Langshan arc in the southern Central Asian Orogenic Belt (CAOB). Zircon U-Pb geochronology indicates that the Urad Houqi suite was formed ca. 278−268 Ma by a continuous magmatic event, marking the final episode of the late Paleozoic continental arc magmatism in the southern CAOB induced by the subduction of the Paleo-Asian Ocean. The suite exhibits a more potassic and fertile composition compared to typical continental arc magmas, and shows a continuous increasing K2O/Na2O ratio and incompatible element concentrations, including K, Rb, Th, and U, with increasing SiO2 content. The geochemical evolution of the suite offers valuable insights into the maturation process of the juvenile arc crust. By combining zircon Hf and bulk-rock Sr-Nd-Pb isotope studies, we demonstrated that the various lithologies within the Urad Houqi suite recorded a three-stage enhancement of the geochemical fertility of the juvenile Langshan arc crust: (1) generation of the initially fertile primary melt from the mantle wedge metasomatized by the recycled subducted sediments; (2) trans-crustal open-system hybridization processes, including peritectic reaction and magma mixing, that further fertilized the derived melts in addition to fractional crystallization; and (3) diverse petrogenetic processes that contributed to the fertility of the granitic magmas, including polybaric fractionation of the primary magma, anatexis of ancient crust, and extraction of high-silica melts from the shallow mush system. The three-stage maturation of the juvenile Langshan arc crust was typically coupled with the accelerated subduction of the Paleo-Asian Ocean since the earliest Permian. Our study finds that the Urad Houqi suite can serve as an example of modern continental crust maturation at continental arc settings.

Magma flux variations triggering shallow-level emplacement of the Takidani pluton (Japan): Insights into the volcanic-plutonic connection

Defining the physico-chemical evolution of felsic magmatic reservoirs in the middle-upper continental crust is crucial for understanding the connection between granitic plutons and silicic volcanic rocks and for research on the trigger mechanisms of volcanic eruptions. In this study, we investigate the formation and evolution of the Hotaka-Takidani volcano-plutonic complex (Japan) using a comprehensive approach that integrates zircon petrochronology with phase-equilibria and thermal modelling. The Takidani pluton is a Pleistocene tilted zoned intrusion of granodiorite-granite composition associated with at least two large caldera-forming eruptions. We have investigated the three uppermost units of the pluton and the oldest and more voluminous Nyukawa eruption, which deposited more than 400 km3 (DRE, dense rock equivalent) of crystal-rich dacitic ignimbrite. Amphibole thermobarometry, U-Pb zircon dates and geochemical data indicate that the Takidani pluton and the Nyukawa dacite were sequentially sourced from a 10–12 km-thick long-lived magmatic reservoir in the middle crust. Thermal evolution of the reservoir led to the extraction of the Nyukawa magma, followed by brief storage in the upper crust and eruption at ca. 1.706 ± 0.091 Ma. For two of the units of the pluton, we combine high-spatial resolution zircon U-Pb geochronology with zircon chemical data. Applying the Ti-in-zircon thermometer, we observe that magma within the reservoir was stored at a relatively constant temperature of ca. 700 °C for 300 kyr and subsequently, starting from 1.4 Ma, the temperature increased progressively reaching 820–850 °C at 1.1 Ma. This temperature shift is inconsistent with a constant magma flux and suggests an increase of the magma injection rate into the reservoir. Using the distribution of zircon ages and thermal simulations, we determined that the reservoir was built by a fairly low average magma flux of ca. 6 × 10−6 km3*km−2*yr−1, shifting to a magma supply rate which is three to four times higher during the last 300 ka of zircon crystallisation. This increase in magma flux caused the final emplacement of the Takidani magma in the upper crust, where it rapidly crystallized. The integration of geochronological data with phase-equilibria and thermo-chemical modelling demonstrates that multiple melt extraction events occurred during the evolution of the Hotaka-Takidani magmatic system. Only the first one triggered a major eruption, while subsequent events resulted in the emplacement of the Takidani pluton in the shallow crust.

Sub-arc Mantle Heterogeneity of the Northern Luzon Volcanic Arc: Mineral and Whole Rock Compositional Variability in Mantle Xenoliths from Lutao Island

Abstract Mantle xenoliths hosted in volcanic rocks from the island of Lutao offer a glimpse into the nature of the mantle beneath the northern Luzon volcanic arc. The xenoliths are spinel-bearing and composed mostly of harzburgite with one lherzolite and one olivine orthopyroxenite. The olivine (Fo92.5–88.9), orthopyroxene (Mg# = 94.6–89.2), and clinopyroxene (Wo49.1–38.1En57.0–45.4Fs3.0–11.0) compositions are similar to those of abyssal peridotites. The spinel compositions are variable and can be principally divided into high-Al (Cr# &amp;lt; 45) and low-Al (Cr# &amp;gt; 45) groupings. The whole rock compositions are similar to abyssal peridotite (Al2O3 = 0.95–2.07 wt %; Mg# = 88.5–90.9) and have U-shaped chondrite normalized rare earth element patterns. The Sr-Nd isotopes of the xenoliths are broadly chondritic (87Sr/86Sri = 0.704400–0.707908; εNd(t) = 0.0 − +1.5). The two-pyroxene equilibrium temperatures range from 900 to 1200 °C with the majority of temperature estimates &amp;gt;1000 °C. The olivine-orthopyroxene-spinel oxygen barometry estimates yielded ΔFMQ values from 0 to +2 and correspond to moderately oxidizing to oxidizing conditions. The xenoliths are likely derived from the Philippine Sea Plate lithospheric mantle that was modified by melt extraction and/or fluid enrichment processes. Trace element and isotopic mixing modeling indicate that 1–2% contamination by subducted South China Sea sediment can explain the Sr-Nd isotopic enrichment and Th and U elemental variability within the xenoliths assuming an initial composition similar to enriched depleted mid-ocean ridge mantle (E-DMM). The anomalously high two-pyroxene equilibrium temperatures of the Lutao xenoliths relative to other regions of the northern Luzon volcanic arc (Iraya &amp;lt;1000 °C) indicate that they were affected by a high-temperature event that was likely a consequence of recent intra-arc rifting that occurred after collision (&amp;lt;6 Ma) between the Luzon arc and the Eurasian margin.

Subduction zone decoupling and rapid forearc spreading at an active continental margin: Evidence from an SSZ-type ophiolite in Qilian Orogen, NW China

The forearc supra-subduction zone (SSZ) ophiolites are of great significance for understanding the geological processes of subduction systems. They are generally considered to mark the initiation of oceanic subduction. Here we present a well-preserved forearc sequence in the Ebo area, North Qilian Accretionary Belt, NW China. Three rock types are identified from bottom to top, including the forearc basalts (FABs), boninites, and calc-alkaline arc rocks, which are similar to those from the Izu-Bonin-Mariana (IBM) area. The FABs display similar compositions to N-MORB but with lower Ti/V ratios (17.16–20.92) and they were derived from 5 to 10 % melting of depleted MORB mantle in the spinel-stability field at P–T conditions of 1.5–2.3 GPa/1442 ± 56 °C after < 2 % melt extraction in the garnet-stability field. The boninites are characterized by significant geochemical characteristics of subduction-zone fluid/melt involvement, which were the products of flux melting of refractory harzburgitic mantle that interacted with slab-derived fluids/ carbonate melt under the P–T conditions of 0.2–1.2 GPa /1316 ± 38 °C. The calc-alkaline rocks show typical arc geochemistry, suggesting the partial melting of the well-established enriched subarc mantle. Zircon U-Pb dating of the boninite and the arc gabbro yielded ages of 507 ± 10 Ma and 483 ± 4 Ma, respectively. Combined with regional data, we infer that the Ebo FABs and boninites are not the products of subduction initiation but formed during subduction-zone decoupling-retreat and rapid forearc spreading following earlier subduction at the North Qilian continental margin, which is subsequently followed by typical continental arc development.

From source to surface: clues from garnet-bearing Carboniferous silicic volcanic rocks, Iberian Pyrite Belt, Portugal

This work investigates the relationships between partial melting, melt extraction, pluton growth and silicic volcanism in garnet-bearing felsic volcanic rocks that were extruded in the Iberian Pyrite Belt, at ca. 345 Ma. The garnets are of peritectic origin, displaying textural and chemical features of disequilibrium crystallization during partial melting reactions involving biotite at high temperatures (up to 870 °C) in the middle-lower crust. Major element composition suggests compositional equilibrium with the entrained and pinitized peritectic cordierite, but reveals some subsequent homogenization by diffusion. Trace element maps and spot analyses of garnet show, nonetheless, significant trace element variations, reflecting biotite and Y-REE-P-rich accessory phase breakdown during partial melting reactions. Peritectic garnet and cordierite growth resulted in the preservation of Th- and Y-rich prograde suprasolidus monazite, which constrains the timing of partial melting of the metapelitic protolith at ca. 356.8 ± 2.4 Ma. The zircon cargo further shows that a significant amount of zircon crystals from previously crystallized felsic melts were also remobilized and erupted. These were likely stored in an upper crustal pluton that grew episodically since ca. 390 Ma during voluminous melt generation periods within the middle to lower crust, which also resulted in voluminous volcanism. The geochemical trends of the felsic volcanic rocks reflect the entrainment of xenoliths of peritectic garnet, cordierite and feldspar, and as such, the garnet-bearing felsic volcanic rocks represent an erupted mixture of a lower-temperature (ca. 770 °C) silicic melt and autocrysts, and peritectic phases and zircon crystals from previously crystallized and stored felsic melts.

A window on the lithospheric mantle beneath the eastern Trans-Mexican volcanic belt: Insights from pargasite-bearing mantle xenocrysts from the Holocene Ocotenco maar volcano

In this short communication, we report and discuss for the first time the occurrence of pargasite-bearing pyroxenite mantle macro-xenocrysts in a basaltic-andesite juvenile material erupted from the Holocene Ocotenco maar at the Serdán-Oriental Basin Volcanic field, situated in the eastern sector of the Trans-Mexican Volcanic Belt. Through an integrated approach based on microtextural observations, mineral chemistry, geothermobarometric modelling, and comparisons with the existing literature, we determined the xenocrysts equilibrated at ca. 0.8 GPa and 1100 °C, compatible with a pargasite-bearing, upper lithospheric mantle setting. The different textures (resorbed vs. unzoned) of xenocrysts, not only provide evidence of a mantle melting event, but also enabled the discrimination between disrupted mantle fragments and peritectic solid phase generated during the pargasite breakdown. Moreover, peculiar characteristics of the xenocrysts such as the dark pargasite rims and the melt embayment in clinopyroxene, together with the lack of phenocrysts from magma fractionation, the swallow-tail texture of groundmass plagioclase, and a vesicle-rich textures, cumulatively suggest a short-term (days to months?) eruptive scenario with a fast melt extraction from the source region to a rapid ascent through the crust and a sudden explosive eruption by the intense interaction of the ascending melt with a source of water close to the surface.

A Machine Learning-Based Thermobarometer for Magmatic Liquids

Abstract Experimentally calibrated models to recover pressures and temperatures of magmas are widely used in igneous petrology. However, large errors, especially in barometry, limit the capacity of these models to resolve the architecture of crustal igneous systems. Here, we apply machine learning to a large experimental database to calibrate new regression models that recover P–T of magmas based on melt composition plus associated phase assemblage. The method is applicable to compositions from basalt to rhyolite, pressures from 0.2 to 15 kbar, and temperatures of 675°C to 1400°C. Testing and optimisation of the model with a filter that removes estimates with standard deviation above the 50th percentile show that pressures can be recovered with root-mean-square-error (RMSE) of 1.1 to 1.3 kbar and errors on temperature estimates of 21°C. Our findings demonstrate that, given constraints on the coexisting mineral assemblage, melt chemistry is a reliable recorder of magmatic variables. This is a consequence of the relatively low thermodynamic variance of natural magma compositions despite their relatively large number of constituent oxide components. We apply our model to two contrasting cases with well-constrained geophysical information: Mount St. Helens volcano (USA), and Askja caldera in Iceland. Dacite whole-rocks from Mount St Helens erupted 1980 to 1986, inferred to represent liquids extracted from cpx–hbl–opx–plag–mt–ilm mush, yield melt extraction source pressures of 5.1 to 6.7 kbar in excellent agreement with geophysical constraints. Melt inclusions and matrix glasses record lower pressures (0.7–3.8 kbar), consistent with magma crystallisation within the upper reaches of the imaged geophysical anomaly and during ascent. Magma reservoir depth estimates for historical eruptions from Askja match the location of seismic wave speed anomalies. Vp/Vs anomalies at 5 to 10 km depth correspond to hot (~990°C) rhyolite source regions, while basaltic magmas (~1120°C) were stored at 7 to 17 km depth under the caldera. These examples illustrate how our model can link petrology and geophysics to better constrain the architecture of volcanic feeding systems. Our model (MagMaTaB) is accessible through a user-friendly web application (https://igdrasil.shinyapps.io/MagmaTaBv4/).