Continental Lithospheric Mantle Research Articles

Mapping the Structure and Metasomatic Enrichment of the Lithospheric Mantle Beneath the Kimberley Craton, Western Australia

AbstractThe lithology, geochemistry, and architecture of the continental lithospheric mantle (CLM) underlying the Kimberley Craton of north‐western Australia has been constrained using pressure‐temperature estimates and mineral compositions for >5,000 newly analyzed and published garnet and chrome (Cr) diopside mantle xenocrysts from 25 kimberlites and lamproites of Mesoproterozoic to Miocene age. Single‐grain Cr diopside paleogeotherms define lithospheric thicknesses of 200–250 km and fall along conductive geotherms corresponding to a surface heat flow of 37–40 mW/m2. Similar geotherms derived from Miocene and Mesoproterozoic intrusions indicate that the lithospheric architecture and thermal state of the CLM has remained stable since at least 1,000 Ma. The chemistry of xenocrysts defines a layered lithosphere with lithological and geochemical domains in the shallow (<100 km) and deep (>150 km) CLM, separated by a diopside‐depleted and seismically slow mid‐lithosphere discontinuity (100–150 km). The shallow CLM is comprised of Cr diopsides derived from depleted garnet‐poor and spinel‐bearing lherzolite that has been weakly metasomatized. This layer may represent an early (Meso to Neoarchean?) nucleus of the craton. The deep CLM is comprised of high Cr2O3 garnet lherzolite with lesser harzburgite, and eclogite. The peridotite components are inferred to have formed as residues of polybaric partial mantle melting in the Archean, whereas eclogite likely represents former oceanic crust accreted during Paleoproterozoic subduction. This deep CLM was metasomatized by H2O‐rich melts derived from subducted sediments and high‐temperature FeO‐TiO2 melts from the asthenosphere.

Steep subduction of the Indian continental mantle lithosphere beneath the eastern Himalaya revealed by gravity anomalies

Steep subduction of the Indian continental mantle lithosphere beneath the eastern Himalaya revealed by gravity anomalies

A Devonian Shoshonitic Appinite–Granite Suite in the North Qinling Orogenic Belt: Implications for Partial Melting of a Water-Fluxed Lithospheric Mantle in an Extensional Setting

AbstractAppinite–granite suites commonly occur in the final stage of collisional orogenic processes, providing a unique opportunity to reveal the properties of continental lithospheric mantle and crust–mantle geodynamics. In this paper, we present a systematic study of the petrology, mineral chemistry, whole-rock geochemistry and geochronology of the Xiong’erling pluton and adjacent appinite dikes in the northern margin of the North Qinling orogenic belt. The pluton is mainly composed of diorites, quartz monzonites and minor granites. The diorites and appinites have LA-ICP-MS zircon U–Pb ages of c. 389 Ma with variable εHf(t) values of −5.58 to +3.36 and TDM1 model ages peaking at c. 1133 Ma. These rocks belong to the shoshonitic series with high Ba–Sr content and were emplaced in an intraplate extensional environment. The quartz monzonites and granites are oxidized A-type granites and have zircon U–Pb ages of c. 396 and c. 378 Ma with negative εHf(t) values of −11.99 to −1.18 and TDM2 model age peaks at 1434, 1720, and 1888 Ma. The in situ87Sr/86Sr ratios of plagioclase from the mafic and felsic rocks are 0.70565-0.70752 and 0.70695-0.71008, respectively. The Xiong’erling pluton and the coeval appinite dikes represent a Devonian appinite–granite suite in the North Qinling orogenic belt. Asthenospheric upwelling during the intraplate extension triggered partial melting of a phlogopite- and amphibole-bearing garnet lherzolite mantle source that had been previously metasomatized by the subducted oceanic slab, leading to the formation of a primitive hydrous mafic magma. The rapid ascent of the water-rich magma along deep-seated active faults with fast crystallization of amphibole resulted in emplacement of the c. 389 Ma appinite dikes. Replacement of anhydrous minerals by amphibole and biotite via hydrous reaction during magma cooling resulted in the formation of the coeval diorites. The mafic magma underplating triggered episodic remelting of a late Paleoproterozoic to early Mesoproterozoic crustal source, leading to the generation of the c. 396–378 Ma quartz monzonites and granites. Combining our results with existing data, we identify a sequence of (1) northward subduction of the Shangdan oceanic crust beneath the Qinling block at c. 524–438 Ma resulting in island-arc calc-alkaline magmatism, (2) closure of the Shangdan Ocean indicated by collision between the Qinling block and the South Qinling terrane and slab failure magmatism at c. 438–410 Ma, and (3) post-collisional to intraplate extension with alkaline magmatism at c. 410–370 Ma. The Devonian extensive intraplate magmatism marks the end of the Paleozoic orogenesis in the North Qinling belt.

Magmatism during the continent – ocean transition

As continents break apart, the dominant mechanism of extension transitions from faulting and lithospheric stretching to magma intrusion and oceanic crust formation in a new ocean basin. A common feature of this evolution preserved at magmatic rifted margins worldwide are voluminous lava flows that erupted close to sea level during the final stages of development of the continent-ocean transition (COT). The mechanisms responsible for the generation of the melts that contribute to these voluminous flows, the so-called seaward dipping reflectors (SDR), and their significance in the context of COT development, are relatively poorly understood; they lie deep below post-rift strata along submarine rifted margins where they cannot be studied directly. Extensive coring of the Afar Stratoid Series - an areally-extensive sequence of Pliocene-aged basalts and intercalated sediments that lie atop the developing COT in the sub-aerial Afar Depression, northern Ethiopia - offers fresh scope to address this issue. We present a numerical model simulating the formation of enriched metasomes within the continental lithospheric mantle by the passage of magmas resembling modern axial basalts. Thermal destabilization of the metasome, caused by plate stretching, initiates melt formation within the metasome. These melts, when mixed with a depleted lithospheric mantle component, closely match the range of compositions of the Afar Stratoid Series lavas in this study. Metasomatic re-enrichment and subsequent melting of the lithospheric mantle during the COT may contribute to further plate thinning. These results demonstrate a novel mechanism by which large-volume flows may be erupted during the COT.

Geochronology, petrogenesis, and magmatic oxidation state of the Mangling intrusive complex, Northern Qinling Belt, Central China: Implications for magma fertility and tectonic setting

The Mangling intrusive complex has different dioritic to granitic phases and is spatially and temporally related to molybdenum deposits in the Qinling Orogen. Zircon U-Pb dating of the Mangling intrusive complex indicates that dioritic rocks (biotite diorite and biotite diorite enclave; ca. 150−147 Ma) formed earlier than granitic rocks (medium- to fine-grained and fine-grained monzogranite and K-feldspar granite; ca. 145−141 Ma). The Mangling dioritic rocks exhibit large ion lithosphere element (e.g., Rb) and light rare earth element enrichment and high field strength element (e.g., Nb, Ta, and Ti) depletion. They have low to moderate SiO2 (51.33−58.16 wt%), high MgO (3.10−4.75 wt%) and Mg# (48−60), and negative zircon εHf(t) values (−11.6 to −6.8), suggesting origination from the continental lithospheric mantle that may have been metasomatized by previous sediment-derived melts and slab-derived fluids constrained by high Nb/Y, Th/Yb, and Rb/Y ratios. The Mangling granitic rocks are I-type granites and have high SiO2 (67.90−81.88 wt%) and low MgO (0.16−0.74 wt%). They have low and negative zircon εHf(t) values (−18.7 to −1.9) and old zircon Hf two-stage model ages (2334−1287 Ma), as well as similar mineral fractionation (e.g., hornblende, biotite, sphene, and apatite) with the Mangling dioritic rocks, indicating that they were derived from the remelting of old crustal rocks (e.g., Xiong’er and Kuanping groups) by the evolved underplated mafic magma. Compared with the Taoguanping mineralized monzogranite in the Northern Qinling Belt, zircon geochemistry (e.g., Ce4+/Ce3+, Eu/Eu*, and ΔFMQ [relative fayalite-magnetite-quartz buffer]) indicates that magma of the Mangling intrusive complex (except the biotite diorite) has high oxygen fugacity and small fractionated granitic intrusions, which are coeval with the biotite diorite enclave or younger than the Mangling granitic rocks, may have potential for generating porphyry molybdenum mineralization. The combination of this study and previous studies corroborates that the Qinling Orogen underwent an intracontinental orogenic evolution in a post-collisional compression to extension transitional setting during the Late Jurassic to Early Cretaceous, affected by far-field Paleo-Pacific slab subduction.

Reaction between volatile-bearing eclogite and harzburgite as a function of degree of interaction: Experimental constraints at 4 GPa

Abstract The mantle is known to be heterogeneous, mainly composed of peridotite and eclogite. Eclogite-derived hydrous melts may interact with harzburgite at the slab-mantle interface in subduction zones or in the sub-continental lithospheric mantle. In this study, such interactions were simulated by performing hybridization experiments in which a layer of eclogite was juxtaposed to a layer of harzburgite in the presence of H2O-CO2 at 4 GPa and 1200 °C, conditions where eclogite is super-solidus while harzburgite is sub-solidus. A diamond trap was placed in between the two layers to trap the fluid or melt phase, allowing direct determination of their composition. The multi-anvil was rotated at different frequencies to examine the effect of increasing degree of interaction on the melt composition as well as the mineral compositions. The interaction of eclogite-derived hydrous melt and harzburgite results in a reaction layer at the interface between the two lithologies, composed of Opx and garnet. The harzburgite above the reaction layer is metasomatized, containing various amounts of olivine, Opx, Cpx, and garnet. The eclogitic melt is modified during this interaction. With increasing interaction, a thicker reaction layer is formed. Both the eclogitic and the peridotitic garnet compositions approach each other and become intermediate between the composition of the garnet in the eclogite+H2O+CO2 system and the garnet in the harzburgite+H2O+CO2 system at these conditions. The Mg# of the peridotitic olivine and Opx decreases with increasing interaction. The initial basaltic melt in equilibrium with eclogite is metaluminous, turning to a peralkaline melt with increasing interaction with the harzburgite. The metasomatizing effect of the eclogite-derived hydrous melt on the harzburgite is observed by increasing the mode of the peridotitic Opx, Cpx, and garnet at the expense of peridotitic olivine and eclogitic garnet. A slight increase in melt fraction occurs as well. This interaction also results in a gradient in the log fO2. Relatively more oxidizing conditions occur near the reaction layer, becoming more reduced into the peridotite, suggesting that the reaction zones act as partial barriers for the melt to travel through the peridotite. Increased interaction leads to higher log fO2 values. These experiments demonstrate the influence of the degree of interaction on the range of melt compositions found in volcanic arcs as well as the degree of metasomatism in the mantle found in the sub continental lithospheric mantle.

Tungsten in the mantle constrained by continental lithospheric peridotites: Less incompatible and more abundant

Isotope dilution W whole-rock data, coupled with data for other trace elements in 42 peridotite xenoliths from Tariat (Mongolia), Yangyuan and Hannuoba (China), and Letlhakane (Botswana) are used to constrain the W abundance and behavior in the continental lithospheric mantle. Tungsten concentrations are >11 ng g−1, even in the most depleted peridotites, and higher than the W abundance estimate of the bulk silicate Earth (BSE) in all but one sample. Combined W-Th-U systematics reveal fractionation of W from Th and U with increasing degree of melt depletion, while incompatible element-enriched samples display no such fractionation. Integrated melt reactions and degrees of melt depletion are estimated based on either residual Th or residual Yb abundances for each sample, relative to a BSE precursor. This allows modelling of different partial melting scenarios that consider maximized trace element retention (Th-based model; where measured Th contents are decoupled from Al2O3 contents) or correction for metasomatic trace element contributions (Yb-based model; where measured Yb contents are well-coupled with Al2O3 contents). Modelled residual W concentrations reveal W excesses of 11–19 ng g−1, even when correcting for metasomatic W contributions. Combined with available in situ trace element mineral data for the rock-forming silicates and spinel, which all display W abundances much higher than expected, these findings suggest that either (i) W is less incompatible than previously thought, or that (ii) the peridotites experienced W metasomatism after the last melt depletion event with subsequent re-equilibration. The latter is unlikely, as similar geochemical patterns are observed for oceanic lithospheric mantle based on published data, requiring such W metasomatism to be a global phenomenon. We thus conclude that W is not as incompatible as Th and U during upper mantle melting and potential W-Th-U fractionation mechanisms are examined to explain near-constant W/Th and W/U in mantle melts. We propose revised W concentrations for the BSE of 27 ± 16 ng g−1 to 31 ± 14 ng g−1 and for the depleted mantle of 22 ± 13 ng g−1.

Evidence and significance of wehrlitization within continental lithospheric mantle

异剥橄榄岩化是地幔橄榄岩与硅不饱和熔体发生强烈相互作用的地幔过程。这一过程通常伴随地幔物理-化学性质的显著改变，与克拉通岩石圈减薄和破坏密切相关。本文梳理了前人对异剥橄榄岩化的驱动机理、结构和成分特征及其对大陆稳定性影响等方面的定性认识，并在此基础上通过热力学模拟，刻画了不同属性地幔与熔体在不同温度压力条件下的相互作用过程，目标是为定量限定异剥橄榄岩化的发生条件、物理-化学效应及地球动力学意义。研究显示，异剥橄榄岩化通常会出现磷灰石、角闪石、金云母和碳酸盐等标志性交代矿物，并以斜方辉石被单斜辉石替换为典型特征，还可发育海绵边、熔体囊（已冷却结晶为单斜辉石、橄榄石等细粒矿物）和矿物环带等非平衡结构。热力学模拟表明，熔岩反应过程受原岩、熔体性质以及熔/岩比例的共同控制，不同属性地幔（难熔方辉橄榄岩、饱满二辉橄榄岩）与贫硅（如霞石岩）熔体反应都能生成富单斜辉石橄榄岩，且在较高压力（如2.5GPa和4.0GPa）条件下更容易形成异剥橄榄岩；相比之下，不同性质地幔与富硅（如MORB）熔体反应在任何压力条件下均不能产生异剥橄榄岩。此外，两类熔体参与的熔岩反应均能引起围岩密度逐渐增加，并在高压条件下能导致饱满岩石圈地幔密度超过地幔参考值、达到岩石圈失稳的条件。因此，本文从热力学的角度定量约束了异剥橄榄岩化的发生条件和对岩石圈地幔的影响，为前人提出的橄榄岩-熔体反应能够造成岩石圈地幔减薄和破坏这一概念模型提供了理论实证，对评估大陆稳定性具有重要意义。

Mantle thermochemical variations beneath the continental United States through petrologic interpretation of seismic tomography

The continental lithospheric mantle plays an essential role in stabilizing continents over long geological time scales. Quantifying spatial variations in thermal and compositional properties of the mantle lithosphere is crucial to understanding its formation and its impact on continental stability; however, our understanding of these variations remains limited. Here we apply the Whole-rock Interpretive Seismic Toolbox For Ultramafic Lithologies (WISTFUL) to estimate thermal, compositional, and density variations in the continental mantle beneath the contiguous United States from MITPS_20, a joint body and surface wave tomographic inversion for Vp and Vs with high resolution in the shallow mantle (60–100 km). Our analysis shows lateral variations in temperature beneath the continental United States of up to 800–900°C at 60, 80, and 100 km depth. East of the Rocky Mountains, the mantle lithosphere is generally cold (350–850°C at 60 km), with higher temperatures (up to 1000°C at 60 km) along the Atlantic coastal margin. By contrast, the mantle lithosphere west of the Rocky Mountains is hot (typically >1000°C at 60 km, >1200°C at 80–100 km), with the highest temperatures beneath Holocene volcanoes. In agreement with previous work, we find that the chemical depletion predicted by WISTFUL does not fully offset the density difference due to temperature. Extending our results using Rayleigh-Taylor instability analysis, implies the lithosphere below the United States could be undergoing oscillatory convection, in which cooling, densification, and sinking of a chemically buoyant layer alternates with reheating and rising of that layer.

Petrogenesis of the Eocene Yulong potassic intrusion in non‐subduction setting in the Sanjiang Tethys

The giant Yulong porphyry Cu (‐Mo‐Au) deposit was formed in a post‐collisional setting in eastern Tibet, which is associated with Eocene potassic to ultrapotassic felsic intrusive rocks. The origin of the Yulong potassic intrusion remains highly disputed. Hence, we present new zircon U–Pb ages, whole‐rock geochemistry and zircon Lu‐Hf isotopes for the Yulong intrusion. Zircon U–Pb dating from the porphyritic monzogranite and porphyritic quartz monzonite yield ages of 41.1 ± 0.2, 40.8 ± 0.1 and 40.5 ± 0.2 Ma, respectively. Samples from the Yulong intrusion exhibit high K2O (4.13–4.85%), and high K2O/Na2O (1.00–1.25), [La/Yb]n (30.37–40.37) and Sr/Y (57.86–87.68) ratios, low Mg#, MgO, Ni, and Cr contents, enrichment in light rare earth elements (LREEs) and depletion of high‐field‐strength elements (HFSEs, e.g., Nb, Ta, and Ti). Zircon grains from this study have mostly positive εHf(t) values (+1.6 to +2.9), and crustal model ages (TDMC) of 0.90–1.19 Ga. These features suggest that the Yulong intrusion is derived from a thickened juvenile lower crust. During the Permian to late Triassic, Jinshajiang‐Ailaoshan Palaeo‐Tethyan oceanic subduction from the eastern margin of the North Qiangtang produced metasomatic domains within the continental lithospheric mantle and lower crust. During the late Eocene, the continuing India‐Asia continental collision may have preferentially overthickened the crust. Convective removals in the lower part of the thickened lithosphere resulted in the hot asthenosphere upwelling along the Jinshajiang‐Ailaoshan Suture, which induced the partial melting of the residual metasomatized lithospheric mantle, as well as the thickened lower crust in the Eocene, which formed the Yulong intrusion.

Subduction-related mantle accretion and makeover revealed by mantle xenoliths at the Pacific margin of NE Eurasia

Subduction is thought to be a major factor in the origin and makeover of continental lithospheric mantle, but essential aspects of these processes continue to be debated. Here, we report textural, chemical and Re Os isotope data on mantle xenoliths from far eastern Russia located at different distances from the Asian Pacific coast to examine how the proximity to subduction zones affects lithospheric mantle compositions. The xenoliths come from the Khanka block close to the border with China and several sites in the Sikhote-Alin orogenic belt (SAOB) near the Japan Sea coast. The inboard Khanka xenoliths are spinel lherzolites with chemical compositions, suggesting that they are residues of low- to moderate-degrees of melt extraction from fertile mantle and essentially unaffected by later events for 1.9 Ga. Such rocks are also common among peridotites in the SAOB, located closer to the subduction margins, but the SAOB suite also includes: a) heterogeneous rocks with veins and patches of pyroxenites or hydrous minerals, and b) rocks enriched in iron (8.8–13.5 wt%), Al 2 O 3 and incompatible minor elements. The positive correlation between elevated 187 Os/ 188 Os ratios (>0.135) and FeO contents in some Fe-rich SAOB xenoliths is interpreted to indicate that the metasomatic melts were derived from the asthenosphere hybridized with recycled materials from the Pacific slab. The lack of correlation between Re -depletion Os model ages (T RD ) and depletion indices, as well as the similarity of Os isotopic compositions between the SAOB peridotites and abyssal peridotites are interpreted to show that the SAOB mantle possibly formed and accreted to Eurasia from the asthenosphere during eastward migration of the Paleo-Pacific subduction trench. Based these results, we infer that the SAOB lithospheric mantle experienced multi-stage metasomatism after its formation, the most significant of which is localized infiltration of subduction-related basaltic melts with recycled characteristics induced by the extension of Japan Sea and release of volatiles from the slab. By contrast, the inboard Khanka lithospheric mantle shows little, if any, evidence for metasomatism. Such contrasting metasomatic modification implies different degrees and styles of makeover processes at different distances relative to the subduction margins. • The Khanka mantle was little affected by metasomatism since it was formed at ~1.9 Ga. • The Sikhote-Alin lithospheric mantle was formed by subduction and accretion. • Mantle metasomatism beneath Sikhote-Alin was caused by melts with recycled signature. • The nature and degrees of metasomatism are related to the distance from subduction front.

Generation of Continental Lithospheric Mantle by Tectonic Isolation of Oceanic Plate

AbstractContinental formation models invoke subduction or plume‐related processes to create the buoyant, refractory character of continental lithospheric mantle (CLM). From similarities in melt depletion, major element composition, modal clinopyroxene, and Os isotope systematics it has been proposed that oceanic mantle lithosphere is the likely protolith to non‐cratonic CLM, however, a direct link between the two has been difficult to ascertain. Using dredged mantle peridotite xenoliths from the Ferrel Seamount, off the west coast of Baja California, Mexico, we show that tectonic isolation of an oceanic plate may lead to formation of non‐cratonic CLM. Ferrel xenoliths are coarse‐grained spinel lherzolite, or rare harzburgite. Bulk‐rock and clinopyroxene trace element compositions reveal two‐stages of melt refertilization following melt depletion, with infiltration by mid‐ocean ridge basalt‐type melts, followed by melt addition from host alkali basalt. Melt depletion correlations with 187Os/188Os and highly siderophile element abundances indicate preserved melt depletion and refertilization processes are ancient. From these observations, the Ferrel xenoliths represent lithosphere from the abandoned Pacific‐Farallon ridge. The history of melt depletion, followed by MORB‐melt refertilization is consistent with the peridotites representing oceanic mantle lithosphere that was subsequently incorporated into the Baja‐Guadalupe microplate during “ridge jump.” These peridotites demonstrate that isolation of oceanic lithosphere that is rafted onto a continental margin provides a viable means for producing non‐cratonic CLM. We suggest that continuation of late‐stage, low degree melt refertilization may provide a link between oceanic lithosphere and non‐cratonic CLM and propose a tectonic model to preserve and facilitate this continued evolution.

High-C content and CO2/Ba ratio of the Earth’s enriched upper mantle

One of the fundamental constraints on the carbon content in the Earth’s mantle comes from depleted mid-ocean ridge basalts (D-MORB) that experienced minimal CO2 degassing. The CO2/Ba ratio of 120 ± 20 in these D-MORB have been used to estimate the CO2 content in more enriched MORB (e.g., N- and E-MORB) that have lost CO2 through degassing. However, it has been unclear as to how constant the CO2/Ba ratio is in D-, N-, and E-MORB, which are known to derive from partial melting of a chemically and isotopically heterogeneous upper mantle. In this study, we constrained the CO2/Ba ratio of the E-MORB mantle source in the Pacific by studying the compositions of MORB glass samples from the southern East Pacific Rise (on– and off-axial SEPR and Garrett transform fault) as well as volcanic glass samples from the Pukapuka ridges and Rano Rahi seamounts that extend to the west of the SEPR. Some of the submarine glasses from the Rano Rahi seamounts are D-MORB with highly depleted CO2 and incompatible trace element contents, but with isotopic composition pointing to a long-term enriched mantle source. These lavas were generated by re-melting of the residual Pacific E-MORB mantle source that had undergone a recent melting episode responsible for the Pukapuka ridge magmatism. Some Rano Rahi seamount lavas experienced limited CO2 degassing during eruption, and we consider that their most reliable pre-degassing CO2/Ba ratio is the maximum value of 241 ± 30, measured in the most trace element depleted but isotopically enriched sample. While we show that this high CO2/Ba ratio is in part due to fractionation during the earlier melting event that generated the Pukapuka ridges, the estimated initial CO2/Ba ratio of the Pacific E-MORB mantle source is 177-52+71, still higher on average than that of the D-MORB mantle source (120 ± 20). We use this CO2/Ba ratio for the Pacific E-MORB mantle source to provide an estimate of its CO2 content of 391-180+455 ppm. We speculate that the high CO2/Ba and C-enriched nature of the Pacific E-MORB mantle source could originate from either (1) recycling of oceanic lithosphere at subduction zones or (2) foundering of metasomatized continental lithospheric mantle into the convecting mantle.

Ultrahigh-pressure peridotites record Neoarchean collisional tectonics

It is debated when plate tectonics first operated on Earth. One of the arguments against the Archean (> 2.5 Ga) operation of plate tectonics is the lack of rock records that can be best explained to be formed at convergent plate boundaries, such as continental lithosphere metamorphosed at ultrahigh-pressures (> 2.7 GPa or 80–100 km). Here we report Archean ultrahigh-pressure peridotites in Eastern Hebei, the North China Craton. Bulk-rock and mineral compositions suggest that these peridotites are likely cumulates or slivers of metasomatized continental lithospheric mantle. Garnet pseudomorphs and pyroxene exsolution textures are preserved in these Archean peridotites, indicating decompression-induced breakdown of the original garnet and pyroxene from high pressures. We reintegrate the original garnet and clinopyroxene based on mass proportions and compositions of decompression-induced breakdown products. The reconstructed garnet and clinopyroxene compositions of these Archean peridotites indicate that they were brought up from mantle depths of 110–130 km. We propose that these ultrahigh-pressure peridotites are tectonic slivers of a collisional complex, possibly subducted to mantle depths and then exhumed to crustal levels during Neoarchean subduction and subsequent arc/continental collision, similar to those from Phanerozoic continental collisional zones. The Archean ultrahigh-pressure peridotites in the North China Craton provide direct evidence for operation of continental collisional plate tectonics since at least 2.5 billion years ago.

Assessing the role of thermal disequilibrium in the evolution of the lithosphere–asthenosphere boundary: an idealized model of heat exchange during channelized melt transport

Abstract. This study explores how the continental lithospheric mantle (CLM) may be heated during channelized melt transport when there is thermal disequilibrium between (melt-rich) channels and surrounding (melt-poor) regions. Specifically, I explore the role of disequilibrium heat exchange in weakening and destabilizing the lithosphere from beneath as melts infiltrate into the lithosphere–asthenosphere boundary (LAB) in intraplate continental settings. During equilibration, hotter-than-ambient melts would be expected to heat the surrounding CLM, but we lack an understanding of the expected spatiotemporal scales and how these depend on channel geometries, infiltration duration, and transport rates. This study assesses the role of heat exchange between migrating material in melt-rich channels and their surroundings in the limit where advective effects are larger than diffusive heat transfer (Péclet numbers > 10). I utilize a 1D advection–diffusion model that includes thermal exchange between melt-rich channels and the surrounding melt-poor region, parameterized by the volume fraction of channels (ϕ), average relative velocity (vchannel) between material inside and outside of channels, channel spacing (d), and timescale of episodic or repeated melt infiltration (τ). The results suggest the following: (1) during episodic infiltration of hotter-than-ambient melt, a steady-state thermal reworking zone (TRZ) associated with spatiotemporally varying disequilibrium heat exchange forms at the LAB. (2) The TRZ grows by the transient migration of a disequilibrium-heating front at a material-dependent velocity, reaching a maximum steady-state width δ proportional to ϕvchannel(τ/d)n, where n≈2 for periodic thermal perturbations and n≈1 for a single finite-duration thermal pulse. For geologically reasonable model parameters, the spatiotemporal scales associated with establishment of the TRZ are comparable with those inferred for the migration of the LAB based on geologic observations within continental intra-plate settings, such as the western US. The results of this study suggest that, for channelized transport speeds of vchannel=1 m yr−1, channel spacings d≈102 m, and timescales of episodic melt infiltration τ≈101 kyr, the steady-state width of the TRZ in the lowermost CLM is ≈10 km. (3) Within the TRZ, disequilibrium heat exchange may contribute ≈10-5 W m−3 to the LAB heat budget.

Sustained indentation in 2-D models of continental collision involving whole mantle subduction

SUMMARY Continental collision zones form at convergent plate boundaries after the negatively buoyant oceanic lithosphere subducts entirely into the Earth’s mantle. Consequently, orogenesis commences, and the colliding continents are sutured together. During the collision, plate convergence and motion of the sutured boundary towards the overriding plate are manifest in its deformation, as is the case for the long-term (∼50 Ma) and nearly constant convergence rate at the India–Eurasia collisional zone that hosts the Himalaya. However, despite the long history of modelling subduction-collision systems, it remains unclear what drives this convergence, especially in models where subduction is driven solely by buoyancy forces. This paper presents dynamic self-consistent buoyancy-driven 2-D whole-mantle scale numerical models of subduction-and-collision processes to explore variations in density and rheological stratification of the colliding continent and overriding plate (OP) viscosity (a proxy for OP strength) that facilitate post-collisional convergence and collisional boundary migration. In models with a moderately buoyant indenting continent, the collisional boundary advance is comparatively low (0.1–0.6 cm yr–1), and convergence is driven by the dense continental lithospheric mantle that continues to subduct as it decouples from its deforming crust. Conversely, models with a highly buoyant indenting continent show sustained indentation at 0.5–1.5 cm yr–1 until the slab detaches. Furthermore, models with a weaker OP and lower backarc viscosity show an enhanced propensity for indentation by a positively buoyant continent. These models additionally highlight the role of whole mantle flow induced by the sinking of the detached slab in the lower mantle as it sustains slow convergence at an average rate of 0.36 cm yr–1 for ∼25 Myr after break-off as well as prevents the residual slab from educting. In previous buoyancy-driven partial mantle depth models such eduction does generally occur, given that free-sinking of the detached slab in the mantle is not modelled. Although these findings widen the understanding of the long-term convergence of indenting continents, the lower post-collisional advance rates (0.3–1.5 cm yr–1) compared to India’s approximate 1000–2000 km of northward indentation during the last 50 Myr attest to the need for 3-D models.

Thermodynamic modeling of melt addition to peridotite: Implications for the refertilization of the non-cratonic continental mantle lithosphere

In a classic model of evolution of the non-cratonic continental mantle lithosphere, harzburgites represent the refractory (<5% clinopyroxene) residues of high degrees of partial melting of fertile mantle, while lherzolites (>5% clinopyroxene) represent residues of lesser degrees of partial melting. However, partial melting is not the only process that could explain the peridotite compositional variability that ranges from fertile (>2 wt% Al2O3, <45 wt% MgO) to refractory (<2 wt% Al2O3, >45 wt% MgO). In the refertilization process, harzburgite is a refractory protolith (potentially previously formed by partial melting of a fertile mantle) that undergoes reactive percolation of silicate melts derived from the underlying asthenosphere, resulting in the crystallization of a new generation of minerals (mostly clinopyroxene). A simple but critical first step towards understanding the refertilization process is to examine how modal and major element compositions evolve as melts are added to peridotites. Here we use a thermodynamically-constrained two-component mixing model to independently evaluate the roles of five different parameters: pressure, temperature, redox conditions, and compositions of the initial peridotite and the added basaltic melt (hereafter referred to P-T-fO2-Xπ-Xmelt), during melt addition. We compare the results with observed suites of peridotites. The main observations are as follows: (1) the produced model is consistent with the global peridotite database, and (2) T, fO2 and small variations of pressure have almost no impact on the evolution of the system. In contrast, the mineralogy of the percolated harzburgite has a substantial effect on the variation of the modal proportions. The parameter with the most significant impact is Xmelt, which is directly linked to the geodynamic context and melting conditions. This parameter directly controls the refertilization reaction and so, the phase proportions and the bulk-rock composition. Elements that partition preferentially in the melt phase (e.g., Na) display depletions in natural assemblages that are stronger than those predicted from the simple mixing model, consistent with the fact that the natural process occurs in an open system, and that reactive percolation likely results in incompatible element enrichment in the associated melt. Our results corroborate the suggestion that most of the spectrum of compositional variability observed in lithospheric mantle peridotites can be explained by the impregnation of primitive silicate melt in refractory harzburgites.

Petrogenesis of the Yeonhwa ultrapotassic intrusions in the Yeongnam Massif—Evidence for enrichment of the Triassic continental lithospheric mantle beneath the Korean peninsula

The Triassic period is when abundant late-orogenic magmatism occurred in the Korean peninsula, but the characteristics of the lithospheric mantle at that time beneath the Yeongnam Massif have not been well-defined because of the scarcity of mafic intrusions. We report zircon U Pb and phlogopite 40 Ar/ 39 Ar ages, petrography, and geochemistry of the Triassic (227.6 Ma) ultrapotassic igneous rocks in the Yeonhwa I Mine (Yeonhwa Ultrapotassic Intrusions, YUI) emplaced into the Taebaeksan Basin in the Yeongnam Massif. The YUI occur as two stocks with different lithologies: phlogopite-rich and amphibole-rich, where the former underwent carbonate-melt interaction along the margin in contact with the carbonate formation. The two lithologies share similar compositional and isotopic characteristics; however, the phlogopite-rich YUI better retains its mantle-derived features, having idiomorphic high-Mg# (up to 91.9) phlogopite and diopside phenocrysts and bulk concentrations of high MgO (Mg# = molar 100 Mg/[Mg + Fe] > 70), Cr (208–1150 ppm), and Ni (48–219 ppm). The enrichment of the YUI in H 2 O (phlogopite-rich), silica contents (50.2–56.9 wt%), and large ion lithophile elements (LILE; K, Rb, Sr, Ba, Pb), together with depletion of high fields strength elements (HFSE; Ti, Zr, Nb, Hf, Ta), are typical of the ultrapotassic rocks in orogenic settings. The elemental and isotopic data imply that the YUI originated from partial melts of the phlogopite-bearing lithospheric mantle near the spinel-garnet transition zone. The extremely enriched isotopic ratios ( 87 Sr/ 86 Sr(t) = 0.70756 to 0.71227, ε Nd (t) = −8.82 to −11.0) suggest a highly enriched mantle source, possibly more enriched than the contemporary pelagic sediments. Considering that the Yeongnam Massif was in a subduction system during Permian-Triassic, not a continental collisional belt, the enrichment was inherited from a combination of both young, subduction-related and ancient, pre-existing metasomes in the continental lithospheric mantle. Our data confirm the presence of the enriched, volatile-rich lithospheric mantle beneath the Triassic Yeongnam Massif. This enriched origin contrasts with earlier ( c . 240–280 Ma) Yeongnam arc plutons that originated from the depleted asthenospheric mantle (ε Nd (t) ≥ +2). Since the enrichment is also characteristic of most of the c . 220–240 Ma igneous rocks around the Korean peninsula, such lithospheric mantle lithology was widespread during the Late Triassic. • The Yeonhwa Ultrapotassic Intrusions (YUI) intruded into the Yeongnam Massif at around 227.6 Ma. • Reaction with country rocks developed metasomatized margins of the intrusion. • The YUI were derived from the phlogopite-bearing sub-continental lithospheric mantle. • The Triassic Yeongnam Massif had the old lithospheric mantle further modified by the Permina-Triassic subduction zone.

Oxygen isotope (δ18O, Δ′17O) insights into continental mantle evolution since the Archean

Oxygen isotopic ratios are largely homogenous in the bulk of Earth’s mantle but are strongly fractionated near the Earth’s surface, thus these are robust indicators of recycling of surface materials to the mantle. Here we document a subtle but significant ~0.2‰ temporal decrease in δ18O in the shallowest continental lithospheric mantle since the Archean, no change in Δ′17O is observed. Younger samples document a decrease and greater heterogeneity of δ18O due to the development and progression of plate tectonics and subduction. We posit that δ18O in the oldest Archean samples provides the best δ18O estimate for the Earth of 5.37‰ for olivine and 5.57‰ for bulk peridotite, values that are comparable to lunar rocks as the moon did not have plate tectonics. Given the large volume of the continental lithospheric mantle, even small decreases in its δ18O may explain the increasing δ18O of the continental crust since oxygen is progressively redistributed by fluids between these reservoirs via high-δ18O sediment accretion and low-δ18O mantle in subduction zones.

Noble gas isotope systematics in the Canary Islands and implications for refractory mantle components

Noble gas isotope systematics of ocean island basalts (OIB) provide evidence for relatively undegassed and primitive mantle sources. These OIB sources partly derive from the deep mantle by virtue of their distinctiveness from mid-ocean ridge basalts (MORB), which dominantly sample upper mantle. New helium, neon and argon isotope data are presented for Canary Islands lavas, carbonatites and cumulate and mantle xenoliths confirming 3He/4He ratios that are the same or lower than MORB, but that are heterogeneous (∼3–9.5RA) within and between islands in the archipelago. Neon and Ar isotope systematics for lavas are mostly within the range of air compositions. Harzburgite xenoliths from Lanzarote, which are interpreted to represent ancient refractory mantle residues, have distinct He-Ne-Ar isotope systematics from lavas or cumulate xenoliths. The harzburgites are characterized by forsteritic olivine (>Fo91) with uniformly low-3He/4He (6.6 ± 0.2RA), high 40Ar/36Ar (630–4900), and have Ne isotope compositions that range between air values or that are similar or more nucleogenic than depleted MORB mantle (DMM). Similar refractory mantle peridotites have been discovered as xenoliths at other OIB localities and are likely to be distinct from continental lithospheric mantle. Refractory mantle (RM), which by virtue of its high melting temperature is difficult to partially melt, may be a significant component in the convecting mantle and has potential to impart a cryptic noble gas signature to partial melts in intraplate, divergent and convergent margin settings. In this sense, RM may represent a ‘sixth mantle component’ after DMM, high-µ (high 238U/204Pb; HIMU), enriched mantle endmembers (EMI, EMII) and the ‘focus zone’ (FOZO). An RM-type component may be presented in the older eastern Canary Islands and possibly in recent rejuvenated volcanism from Teide (Tenerife). Canary Island intraplate volcanism samples multiple mantle components, confirming that the most gas-rich mantle sources involved in magmatism dominate OIB noble gas compositions.