Hello Guest User, You are visiting this website from a computer with an IP address of 172.70.131.62 with the name of '?' since Tue Apr 23, 2024 at 9:08:53 AM PT for approx. 0 minutes now.
The Sheahan Diamond Literature Reference Compilation - Scientific and Media Articles based on Major Keyword - Subduction
The Sheahan Diamond Literature Reference Compilation is compiled by Patricia Sheahan who publishes on a monthly basis a list of new scientific articles related to diamonds as well as media coverage and corporate announcements called the Sheahan Diamond Literature Service that is distributed as a free pdf to a list of followers. Pat has kindly agreed to allow her work to be made available as an online digital resource at Kaiser Research Online so that a broader community interested in diamonds and related geology can benefit. The references are for personal use information purposes only; when available a link is provided to an online location where the full article can be accessed or purchased directly. Reproduction of this compilation in part or in whole without permission from the Sheahan Diamond Literature Service is strictly prohibited. Return to Diamond Keyword Index
Sheahan Diamond Literature Reference Compilation - Scientific Articles by Author for all years
Each article reference in the SDLRC is tagged with one or more key words assigned by Pat Sheahan to highlight the main topics of the article. In an effort to make it easier for users to track down articles related to a specific topic, KRO has extracted these key words and developed a list of major key words presented in this Key Word Index to which individual key words used in the article reference have been assigned. In most of the individual Key Word Reports the references are in crhonological order, though in some such as Deposits the order is first by key word and then chronological. Only articles classified as "technical" (mainly scientific journal articles) and "media" (independent media articles) are included in the Key Word Index. References that were added in the most recent monthly update are highlighted in yellow.
Subduction is a subset of "plate tectonics" which deals with what happens when oceanic and continental plates converge, namely the denser oceanic plate descends beneath the lighter continental crust.
Early to mid-Tertiary inverted metamorphic gradient under the ColoradoPlateau: evidence from eclogite xenoliths in ultramafic microbreccias, Navajo volcanic field
Journal of Geophysical Research, Vol. Paper # 91JB00284
Nappes, tectonics of oblique plate convergence, and metamorphic evolution related to 140 million years of continuous subduction, Franciscan Complex, California
Journal of Geology, Vol. 100, No. 1, January pp. 19-40
Nd isotope systematics of 2.7 Ga adakites, magnesian andesites and arc basalts, Superior Province: evidence for shallow crustal recycling at Archean subduction zones
Earth and Planetary Science Letters, Vol. 202, 2, pp. 345-60.
Long tern memory of subduction processes in the lithospheric mantle: evidence from the geochemistry of basic dykes in the Gardar Province of South Greenland.
Journal of the Geological Society, Vol. 159, 6, pp. 705-714.
Computation of phase equilibration temperatures by linear programming: a tool for geodynamic modeling and its application to subduction zone decarbonation.
Earth and Planetary Science Letters, Advanced in press,
Computation of phase equilibration temperatures by linear programming: a tool for geodynamic modeling and its application to subduction zone decarbonation.
Earth and Planetary Science Letters, Advanced in press,
In situ x-ray observations of phase assemblages in peridotite and basalt compositions at lower mantle conditions: implications for density of subducted...
Recycling of continental crust at modern subduction zones, implications for Precambrian crustal growth, supercontinent constructions, and littering the mantle
Geological Society of America Annual Meeting ABSTRACTS, Nov. 7-10, Paper 82-3, Vol. 36, 5, p. 205.
Evidence from ophiolites, blueschists, and ultrahigh pressure metamorphic terranes that the modern episode of subduction tectonics began in Neoproterozoic time.
Wang, Q., Ji, S., Salisbury, M.H., Xia, B., Pan, M., Xu, Z.
Pressure dependence and anisotropy of P wave velocities in ultrahigh pressure metamorphic rocks from the Dabie Sulu orogenic belt: implications for seismic properties
High pressure intermediate temperature metamorphism in the southern Barbarton granitoid greenstone terrain, South Africa: a consequence of subduction driven ...
Benn, K., Mareschal, J-C., Condie, K.C. Archean Geodynamics and Environments, AGU Geophysical Monograph, No. 164, pp. 239-254.
Melting of subducted continent: element and isotopic evidence for a genetic relationship between Neoproterozoic and Mesozoic granitoids in the Sulu orogen.
Chemical Geology, Vol. 229, 4, May 30, pp. 227-256.
Evolution from oceanic subduction to continental collision: a case study from the northern Tibetan Plateau based on geochemical and geochronological data.
A subduction wedge origin for Paleoarchean peridotitic diamonds and harzburgites from the PAnd a kimberlite, Slave Craton: evidence from Re Os isotope systematics.
Contributions to Mineralogy and Petrology, Vol. 152, 3, pp. 275-294.
Mainprice, D., Le Page, Y., Rodgers, J., Jouanna, P.
Predicted elastic properties of hydrous D phase at mantle pressures: implications for the anisotropy of subducted slabs near 670 km discontinuity and in the lower mantle.
Earth and Planetary Science Letters, Vol. 259, 3-4, pp. 283-296.
Mainprice, D., Le Page, Y., Rodgers, J., Jouanna, P.
Predicted elastic properties of hydrous D phase at mantle pressures: implications for the anisotropy of subducted slabs near 670 km discontinuity and in the lower mantle.
Earth and Planetary Science Letters, Vol. 259, 3-4, pp. 283-296.
Mo, X., Zhao, Z., Deng, J., Flower, M., Yu, X., Luo, Z., Li, Y., Zhou, S., Deng, G., Zhu, D.
Petrology and geochemistry of post collisional volcanic rocks from the Tibetan plateau: implications for lithosphere heterogeneity and collision induced mantle
Geological Society of America, Special Paper, No. 409, pp. 507-530.
Crustal recycling at modern subduction zones applied to the past-issues of growth and preservation of continental basement crust, mantle geochemistry reconstruction.
Geological Society of America, Memoir, No. 200, pp. 9-32.
Conditions for a crustal block to be sheared off from subducted continental lithosphere: what is an essential factor to cause features associated with collision?
Journal of Geophysical Research, Vol. 113, B4, B04414
Conditions for a crustal block to be sheared off from the subducted continental lithosphere: what is an essential factor to cause features associated with collision?
Journal of Geophysical Research, Vol. 113, B004414.
Polyphase formation and exhumation of high to ultrahigh pressure rocks in continental subduction zones: numerical modeling and application to the Sulu ultrahigh pressure terrane in eastern China.
Journal of Geophysical Research, Vol. 114. B9, B09406
Evolution of the slab bending radius and the bending dissipation in three dimensional subduction models with a variable slab to upper mantle viscosity ratio.
Earth and Planetary Science Letters, Vol. 288, 1-2, pp. 309-319.
Metamorphic evolution of the Naga Hills eclogite and blueschist northeast India: implications for early subduction of the Indian Plate under Burma microplate
Journal of Metamorphic Geology, Vol. 28, 2, pp. 209-225.
Geochemistry of mafic dikes in the Singhbhum Orissa craton: implications for subuction related metasomatism of the mantle beneath the eastern Indian craton.
International Geology Review, Vol. 52, 1, pp. 79-94.
Antigorite peridotite, metaserpentinite, and other inclusions within diatremes on the Colorado Plateau, SW USA: implications for the mantle wedge during low angle
Variation of mineral composition, fabric and oxygen fugacity from massive to foliated eclogites during exhumation of subducted ocean crust in North Qiilian suture
Journal of Metamorphic Geology, Vol. 29, 7, pp. 699-720.
Role of iron and reducing conditions on the stability of dolomite + coesite between 4.25 and 6 GPa - a potential mechanism for diamond formation during subduction
European Journal of Mineralogy, Vol. 23, 1, pp. 5-16.
The alkaline ultrabasic magmatism of the Onega peninsula Nenoksa fields - reflection (display) of the plume and subduction processes in Belomorsky region.
Vladykin, N.V., Deep Seated Magmatism: its sources and plumes, pp. 145-163.
Delamination and ultradeep subduction of continental crust: constraints from elastic wave velocity and density measurement in ultrahigh pressure metamorphic rocks
Journal of Metamorphic Geology, Vol. 29, 7, pp. 781-801.
High pressure phase transitions and subduction behaviour of continent crust at pressure temperature conditions up to the upper part of the lower mantle.
Earth and Planetary Science Letters, Vol. 357-358, pp. 31-41.
Martin, A.M., Laporte, D., Koga, K.T., Kawamoto, T., Hammouda, T.
Experimental stidy of the stability of a dolomite + coesite assembalge in contact with peridotite: implications for sediment-mantle interaction and diamond formation during subduction.
Carbonatites out of a subducted altered oceanic crust? New experimental evidences for "low temperature" carbonatitic melts in COH bearing gabbros at 3.8-4.2 Gpa.
The importance of crystal chemistry on REE partitioning between mantle minerals ( garnet, clinopyroxene, orthopyroxene, and olivine) and basaltic melts.
Armytage, R.M.G., Brandon, A.D., Peslier, A.H., Lapen, T.J.
Osmium isotope evidence for Early to Middle Proterozoic mantle lithosphere stabilization and concommitant production of juvenile crust in Dish Hill, CA peridotite xenoliths.
Geochimica et Cosmochimica Acta, Vol. 137, pp. 113-133.
Generation of hydrous plumes in the mantle transition zone linked to the tectonic erosion of continental crust.
V.S. Sobolev Institute of Geology and Mineralogy Siberian Branch Russian Academy of Sciences International Symposium Advances in high pressure research: breaking scales and horizons ( Courtesy of N. Poikilenko), Held Sept. 22-26, 3p. Abstract
Russian Geology and Geophysics, Vol. 56, pp. 978-995.
Mantle
Subduction
Abstract: The Hadean and Archean geologic history of the Earth is discussed in the context of available knowledge from different sources: space physics and comparative planetology; isotope geochronology; geology and petrology of Archean greenstone belts (GB) and tonalite-trondhjemite-granodiorite (TTG) complexes; and geodynamic modeling review to analyse plate-tectonic, plume activity, and impact processes. Correlation between the age peaks of terrestrial Hadean-Early Archean zircons and late heavy bombardment events on the Moon, as well as the Hf isotope composition of zircons indicating their mostly mafic sources, hint to an important role of impact processes in the Earth’s history between 4.4 and 3.8 Ga. The earliest continental crust (TTG complexes) formed at 4.2 Ga (Acasta gneisses), while its large-scale recycling left imprint in Hf isotope signatures after 3.75 Ga. The associations and geochemistry of rocks suggest that Archean greenstone belts formed in settings of rifting, ocean floor spreading, subduction, and plume magmatism generally similar to the present respective processes. The Archean history differed in the greater extent of rocks derived from mantle plumes (komatiites and basalts), boninites, and adakites as well as in shorter subduction cycles recorded in alternation of typical calc-alkaline andesite-dacite-rhyolite and adakite series that were generated in a hotter mantle with more turbulent convection and unsteady subduction. The Archean is interpreted as a transient period of small plate tectonics.
Abstract: The infiltration of fluids into continental lithospheric mantle is a key mechanism for controlling abrupt changes in the chemical and physical properties of the lithospheric root1, 2, as well as diamond formation3, yet the origin and composition of the fluids involved are still poorly constrained. Such fluids are trapped within diamonds when they form4, 5, 6, 7 and so diamonds provide a unique means of directly characterizing the fluids that percolate through the deep continental lithospheric mantle. Here we show a clear chemical evolutionary trend, identifying saline fluids as parental to silicic and carbonatitic deep mantle melts, in diamonds from the Northwest Territories, Canada. Fluid–rock interaction along with in situ melting cause compositional transitions, as the saline fluids traverse mixed peridotite–eclogite lithosphere. Moreover, the chemistry of the parental saline fluids—especially their strontium isotopic compositions—and the timing of host diamond formation suggest that a subducting Mesozoic plate under western North America is the source of the fluids. Our results imply a strong association between subduction, mantle metasomatism and fluid-rich diamond formation, emphasizing the importance of subduction-derived fluids in affecting the composition of the deep lithospheric mantle.
Abstract: The lithospheric mantle beneath West Antarctica has been characterized using petrology, whole-rock and mineral major element geochemistry, whole-rock trace element chemistry and Mössbauer spectroscopy data obtained on a suite of peridotite (lherzolite and harzburgite) and pyroxenite xenoliths from the Mount Morning eruptive centre, Southern Victoria Land. The timing of pyroxenite formation in Victoria Land overlaps with subduction of the Palaeo-Pacific plate beneath the Gondwana margin and pyroxenite is likely to have formed when fluids derived from, or modified by, melting of the subducting, eclogitic, oceanic crustal plate percolated through peridotite of the lithospheric mantle. Subsequent melting of lithospheric pyroxenite veins similar to those represented in the Mount Morning xenolith suite has contributed to the enriched trace element (and isotope) signatures seen in Cenozoic volcanic rocks from Mount Morning, elsewhere in Victoria Land and Zealandia. In general, the harzburgite xenoliths reflect between 20 and 30% melt depletion. Their depleted element budgets are consistent with Archaean cratonization ages and they have mantle-normalized trace element patterns comparable with typical subcontinental lithospheric mantle. The spinel lherzolite mineral data suggest a similar amount of depletion to that recorded in the harzburgites (20-30%), whereas plagioclase lherzolite mineral data suggest <15% melt depletion. The lherzolite (spinel and plagioclase) xenolith whole-rocks have compositions indicating <20% melt depletion, consistent with Proterozoic to Phanerozoic cratonization ages, and have mantle-normalized trace element patterns comparable with typical depleted mid-ocean ridge mantle. All peridotite xenoliths have undergone a number of melt-rock reaction events. Melting took place mainly in the spinel peridotite stability field, but one plagioclase peridotite group containing high-sodium clinopyroxenes is best modelled by melting in the garnet field. Median oxygen fugacity estimates based on Mössbauer spectroscopy measurements of spinel and pyroxene for spinel-facies conditions in the rifted Antarctic lithosphere are -0·6 ?log fO2 at Mount Morning and –1·0 ± 0·1 (1?) ?log fO2 for all of Victoria Land, relative to the fayalite-magnetite-quartz buffer. These values are in good agreement with a calculated global median value of -0·9 ± 0·1 (1?) ?log fO2 for mantle spinel-facies rocks from continental rift systems.
Abstract: Volcanic arcs above subduction zones are enriched in volatiles and fluid-mobile elements with respect to mid-oceanic ridge basalts. There is general consensus that this particular subduction zone signature is generated by fluid-induced extraction of these elements from subducted oceanic crust and its sedimentary cover. However, how these fluids are transferred through the mantle wedge to the locus of partial melting and what modification the fluids will experience is unresolved. Here we investigate the interaction of slab fluids with the mantle wedge through a series of high-pressure experiments. We explore two end-member processes of focused and porous reactive flow of hydrous slab melts through the mantle. Transfer by porous flow leads to the formation of hydrous minerals that sequester fluid-mobile elements and residual fluids characterized by trace element patterns inconsistent with typical arc lavas. In contrast, no hydrous minerals are formed in the reaction zone of experiments mimicking focused flow, and the typical trace element signature acquired during fluid extraction from the slab is preserved, indicating that this is an efficient process for element transfer through the mantle wedge.
Abstract: The formation mechanisms for early Archaean continental crust are controversial. Continental crust may have accumulated via horizontal accretion in modern-style subduction zones or via vertical accretion above upper mantle upwelling zones. However, the characteristics of the continental crust changes at the transition between the Archaean and Proterozoic eons, suggesting that continental crust did not form in subduction zones until at least the late Archaean. Here I use seismic receiver function data to analyse the bulk properties of continental crust in Western Australia, which formed and stabilized over a billion years in the Archaean. I find that the bulk seismic properties of the crust cluster spatially, with similar clusters confined within the boundaries of tectonic terranes. I use local Archaean crustal growth models to show that both plume and subduction processes may have had a role in creating crust throughout the Archaean. A correlation between crustal age and the bulk seismic properties of the crust reveals a trend: from about 3.5?Gyr ago (Ga) to the end of the Archaean, the crust gradually thickened and simultaneously became more evolved in composition. I propose that this trend reflects the transition between crust dominantly formed above mantle plumes, to crust formed in subduction zones-a transition that may reflect secular cooling of Earth’s mantl
Journal of Metamorphic Geology, Vol. 33, 6, pp. 649-670.
Mantle
Subduction
Abstract: Subduction of hydrated lithospheric mantle introduces Hinline imageO, ferric iron, oxidized carbon and sulphur to the subduction zone system. The fate of these components is poorly known, but is intimately linked to the global geochemical cycles of iron, carbon and sulphur, the genesis of arc-related ore deposits, the temporal evolution of mantle redox state and subduction-related earthquakes and magmatism. thermocalc is used to provide first-order constraints on the effect of subduction zone metamorphism on metamorphic redistribution of iron, carbon, sulphur and water in ultramafic rocks via construction of P?T and T-X(O) pseudosections with open system calculation of the effect of fluid loss. The calculations replicate observed mineral assemblages in high-P to low-T ultramafic rocks at P?T conditions consistent with those suggested by other workers. The results are consistent with open system fluid loss without significant fluid infiltration. Water loss is complete by 850 inline imageC, the corresponding depth of fluid loss being consistent with that inferred for earthquakes in subducting slabs. Losses of carbon and sulphur are relatively minor, at around <5% and <1%, respectively, so it is envisaged that most carbon and sulphur subducted in ultramafic lithologies is transported to >5 GPa, below the depths of the source zone for arc volcanoes. Oxygen activity for rocks in closed systems that evolve with a fixed redox budget is calculated to change from ?FMQ ?1 at 350 inline imageC to over ?FMQ +3 at 850 inline imageC. This result emphasizes the need to consider redox budget as well as oxygen activity when the results of experiments performed at fixed oxygen activity relative to some buffer are interpreted in the context of natural systems. In open systems, devolatilization is calculated to increase the redox budget and oxygen activity of the residue via loss of methane and Hinline imageS at the brucite-out and serpentine-out reactions respectively. No fluid-induced mechanism for oxidation of sub-arc mantle by transfer of redox budget from hydrated ultramafic lithologies to the overlying sub-arc mantle was identified, although further thermodynamic data on fluid species such as Sinline image are required to confirm this.
Eakin, C.M., Long, M.D., Scire, A., Beck, S.L., Wagner, L.S., Zandt, G., Tavera, H.
Internal deformation of the subducted Nazca slab inferred from seismic anisotropy. ..new study suggests that the Earth's rigid tectonic plates stay strong when they slide under another plate, known as subduction, may not be universal.
Abstract: Interactions between crustal and mantle reservoirs dominate the surface inventory of volatile elements over geological time, moderating atmospheric composition and maintaining a life-supporting planet. While volcanoes expel volatile components into surface reservoirs, subduction of oceanic crust is responsible for replenishment of mantle reservoirs. Many natural, 'superdeep' diamonds originating in the deep upper mantle and transition zone host mineral inclusions, indicating an affinity to subducted oceanic crust. Here we show that the majority of slab geotherms will intersect a deep depression along the melting curve of carbonated oceanic crust at depths of approximately 300 to 700 kilometres, creating a barrier to direct carbonate recycling into the deep mantle. Low-degree partial melts are alkaline carbonatites that are highly reactive with reduced ambient mantle, producing diamond. Many inclusions in superdeep diamonds are best explained by carbonate melt-peridotite reaction. A deep carbon barrier may dominate the recycling of carbon in the mantle and contribute to chemical and isotopic heterogeneity of the mantle reservoir.
Abstract: Forty-one diamonds sourced from the Juina-5 kimberlite pipe in Southern Brazil, which contain optically identifiable inclusions, have been studied using an integrated approach. The diamonds contain <20 ppm nitrogen (N) that is fully aggregated as B centres. Internal structures in several diamonds revealed using cathodoluminescence (CL) are unlike those normally observed in lithospheric samples. The majority of the diamonds are composed of isotopically light carbon, and the collection has a unimodal distribution heavily skewed towards ?13C ~ ?25 ‰. Individual diamonds can display large carbon isotope heterogeneity of up to ~15 ‰ and predominantly have isotopically lighter cores displaying blue CL, and heavier rims with green CL. The light carbon isotopic compositions are interpreted as evidence of diamond growth from abiotic organic carbon added to the oceanic crust during hydrothermal alteration. The bulk isotopic composition of the oceanic crust, carbonates plus organics, is equal to the composition of mantle carbon (?5 ‰), and we suggest that recycling/mixing of subducted material will replenish this reservoir over geological time. Several exposed, syngenetic inclusions have bulk compositions consistent with former eclogitic magnesium silicate perovskite, calcium silicate perovskite and NAL or CF phases that have re-equilibrated during their exhumation to the surface. There are multiple occurrences of majoritic garnet with pyroxene exsolution, coesite with and without kyanite exsolution, clinopyroxene, Fe or Fe-carbide and sulphide minerals alongside single occurrences of olivine and ferropericlase. As a group, the inclusions have eclogitic affinity and provide evidence for diamond formation at pressures extending to Earth’s deep transition zone and possibly the lower mantle. It is observed that the major element composition of inclusions and isotopic compositions of host Juina-5 diamonds are not correlated. The diamond and inclusion compositions are intimately related to subducted material and record a polybaric growth history across a depth interval stretching from the lower mantle to the base of the lithosphere. It is suggested that the interaction of slab-derived melts and mantle material combined with subsequent upward transport in channelised networks or a buoyant diapir explains the formation of Juina-5 diamonds. We conclude that these samples, despite originating at great mantle depths, do not provide direct information about the ambient mantle, instead, providing a snapshot of the Earth’s deep carbon cycle.
Abstract: The formation of the Earth's continents is enigmatic. Volcanic arc magmas generated above subduction zones have geochemical compositions that are similar to continental crust, implying that arc magmatic processes played a central role in generating continental crust. Yet the deep crust within volcanic arcs has a very different composition from crust at similar depths beneath the continents. It is therefore unclear how arc crust is transformed into continental crust. The densest parts of arc lower crust may delaminate and become recycled into the underlying mantle. Here we show, however, that even after delamination, arc lower crust still has significantly different trace element contents from continental lower crust. We suggest that it is not delamination that determines the composition of continental crust, but relamination. In our conceptual model, buoyant magmatic rocks generated at arcs are subducted. Then, upon heating at depth, they ascend and are relaminated at the base of the overlying crust. A review of the average compositions of buoyant magmatic rocks — lavas and plutons — sampled from the Aleutians, Izu-Bonin-Marianas, Kohistan and Talkeetna arcs reveals that they fall within the range of estimated major and trace elements in lower continental crust. Relamination may thus provide an efficient process for generating lower continental crust.
Abstract: Geological-geochemical evidence point towards higher mantle potential temperature and a different type of tectonics (global plume-lid tectonics) in the early Earth (>3.2 Ga) compared to the present day (global plate tectonics). In order to investigate tectono-magmatic processes associated with plume-lid tectonics and crustal growth under hotter mantle temperature conditions, we conduct a series of 3D high-resolution magmatic-thermomechanical models with the finite-difference code I3ELVIS. No external plate tectonic forces are applied to isolate 3D effects of various plume-lithosphere and crust-mantle interactions. Results of the numerical experiments show two distinct phases in coupled crust-mantle evolution: (1) a longer (80-100 Myr) and relatively quiet ‘growth phase’ which is marked by growth of crust and lithosphere, followed by (2) a short (?20 Myr) and catastrophic ‘removal phase’, where unstable parts of the crust and mantle lithosphere are removed by eclogitic dripping and later delamination. This modelling suggests that the early Earth plume-lid tectonic regime followed a pattern of episodic growth and removal also called episodic overturn with a periodicity of ?100 Myr.
Abstract: Whether subduction initiation is a forced or spontaneous process is a fundamental geological problem. Using drill core data from the Amami Sankaku basin, Arculus et al.1 suggest that subduction in the Izu–Bonin–Mariana arc began spontaneously. Here we argue that the evidence presented favours neither spontaneous nor forced
Abstract: We reconstructed the spatial distribution of eclogites in the cratonic mantle based on thermobarometry for ~ 240 xenoliths in 4 kimberlite pipes from different parts of the Slave craton (Canada). The accuracy of depth estimates is ensured by the use of a recently calibrated thermometer, projection of temperatures onto well-constrained local peridotitic geotherms, petrological screening for unrealistic temperature estimates, and internal consistency of all data. The depth estimates are based on new data on mineral chemistry and petrography of 148 eclogite xenoliths from the Jericho and Muskox kimberlites of the northern Slave craton and previously reported analyses of 95 eclogites from Diavik and Ekati kimberlites (Central Slave). The majority of Northern Slave eclogites of the crustal, subduction origin occurs at 110-170 km, shallower than in the majority of the Central Slave crustal eclogites (120-210 km). The identical geochronological history of these eclogite populations and the absence of steep suture boundaries between the central and northern Slave craton suggest the lateral continuity of the mantle layer relatively rich in eclogites. We explain the distribution of eclogites by partial preservation of an imbricated and plastically dispersed oceanic slab formed by easterly dipping Proterozoic subduction. The depths of eclogite localization do not correlate with geophysically mapped discontinuities. The base of the depleted lithosphere of the Slave craton constrained by thermobarometry of peridotite xenoliths coincides with the base of the thickened lithospheric slab, which supports contribution of the recycled oceanic lithosphere to formation of the cratonic root. Its architecture may have been protected by circum-cratonic subduction and shielding of the shallow Archean lithosphere from the destructive asthenospheric metasomatism.
Progress in Earth and Planetary Science, Vol. 3, Apr. 28, P. 9-
Mantle
Subduction
Abstract: In the forearc region, aqueous fluids are released from the subducting slab at a rate depending on its thermal state. Escaping fluids tend to rise vertically unless they meet permeability barriers such as the deformed plate interface or the Moho of the overriding plate. Channeling of fluids along the plate interface and Moho may result in fluid overpressure in the oceanic crust, precipitation of quartz from fluids, and low Poisson ratio areas associated with tremors. Above the subducting plate, the forearc mantle wedge is the place of intense reactions between dehydration fluids from the subducting slab and ultramafic rocks leading to extensive serpentinization. The plate interface is mechanically decoupled, most likely in relation to serpentinization, thereby isolating the forearc mantle wedge from convection as a cold, potentially serpentinized and buoyant, body. Geophysical studies are unique probes to the interactions between fluids and rocks in the forearc mantle, and experimental constrains on rock properties allow inferring fluid migration and fluid-rock reactions from geophysical data. Seismic velocities reveal a high degree of serpentinization of the forearc mantle in hot subduction zones, and little serpentinization in the coldest subduction zones because the warmer the subduction zone, the higher the amount of water released by dehydration of hydrothermally altered oceanic lithosphere. Interpretation of seismic data from petrophysical constrain is limited by complex effects due to anisotropy that needs to be assessed both in the analysis and interpretation of seismic data. Electrical conductivity increases with increasing fluid content and temperature of the subduction. However, the forearc mantle of Northern Cascadia, the hottest subduction zone where extensive serpentinization was first demonstrated, shows only modest electrical conductivity. Electrical conductivity may vary not only with the thermal state of the subduction zone, but also with time for a given thermal state through variations of fluid salinity. High-Cl fluids produced by serpentinization can mix with the source rocks of the volcanic arc and explain geochemical signatures of primitive magma inclusions. Signature of deep high-Cl fluids was also identified in forearc hot springs. These observations suggest the existence of fluid circulations between the forearc mantle and the hot spring hydrothermal system or the volcanic arc. Such circulations are also evidenced by recent magnetotelluric profiles.
Earth and Planetary Science Letters, Vol. 449, pp. 118-126.
Gondwana, Rodinia
Subduction
Abstract: Periodic assembly and dispersal of continental fragments has been a characteristic of the solid Earth for much of its history. Geodynamic drivers of this cyclic activity are inferred to be either top-down processes related to near surface lithospheric stresses at plate boundaries or bottom-up processes related to mantle convection and, in particular, mantle plumes, or some combination of the two. Analysis of the geological history of Rodinian crustal blocks suggests that internal rifting and breakup of the supercontinent were linked to the initiation of subduction and development of accretionary orogens around its periphery. Thus, breakup was a top-down instigated process. The locus of convergence was initially around north-eastern and northern Laurentia in the early Neoproterozoic before extending to outboard of Amazonia and Africa, including Avalonia-Cadomia, and arcs outboard of Siberia and eastern to northern Baltica in the mid-Neoproterozoic (?760 Ma). The duration of subduction around the periphery of Rodinia coincides with the interval of lithospheric extension within the supercontinent, including the opening of the proto-Pacific at ca. 760 Ma and the commencement of rifting in east Laurentia. Final development of passive margin successions around Laurentia, Baltica and Siberia was not completed until the late Neoproterozoic to early Paleozoic (ca. 570-530 Ma), which corresponds with the termination of convergent plate interactions that gave rise to Gondwana and the consequent relocation of subduction zones to the periphery of this supercontinent. The temporal link between external subduction and internal extension suggests that breakup was initiated by a top-down process driven by accretionary tectonics along the periphery of the supercontinent. Plume-related magmatism may be present at specific times and in specific places during breakup but is not the prime driving force. Comparison of the Rodinia record of continental assembly and dispersal with that for Nuna, Gondwana and Pangea suggests grouping into two supercycles in which Nuna and Gondwana underwent only partial or no break-up phase prior to their incorporation into Rodinia and Pangea respectively. It was only after this final phase of assembly that the supercontinents then underwent full dispersal.
Geophysical Research Letters, Vol. 43, 10, pp. 4945-4953.
Mantle
Subduction
Abstract: Advances in global seismic tomography have increasingly motivated identification of subducted lithosphere in Earth’s deep mantle, creating novel opportunities to link plate tectonics and mantle evolution. Chief among those is the quest for a robust subduction reference frame, wherein the mantle assemblage of subducted lithosphere is used to reconstruct past surface tectonics in an absolute framework anchored in the deep Earth. However, the associations heretofore drawn between lower mantle structure and past subduction have been qualitative and conflicting, so the very assumption of a correlation has yet to be quantitatively corroborated. Here we show that a significant, time-depth progressive correlation can be drawn between reconstructed subduction zones of the last 130 Myr and positive S wave velocity anomalies at 600 -2300 km depth, but that further correlation between greater times and depths is not presently demonstrable. This correlation suggests that lower mantle slab sinking rates average between 1.1 and 1.9 cmyr 1.
Earth and Planetary Science Letters, Vol. 446, pp. 137-148.
Mantle
Coesite, subduction
Abstract: The ability to accurately constrain the secular record of high- and ultra-high pressure metamorphism on Earth is potentially hampered as these rocks are metastable and prone to retrogression, particularly during exhumation. Rutile is among the most widespread and best preserved minerals in high- and ultra-high pressure rocks and a hitherto untested approach is to use mineral inclusions within rutile to record such conditions. In this study, rutiles from three different high- and ultrahigh-pressure massifs have been investigated for inclusions. Rutile is shown to contain inclusions of high-pressure minerals such as omphacite, garnet and high silica phengite, as well as diagnostic ultrahigh-pressure minerals, including the first reported occurrence of exceptionally preserved monomineralic coesite in rutile from the Dora -Maira massif. Chemical comparison of inclusion and matrix phases show that inclusions generally represent peak metamorphic assemblages; although rare prograde phases such as titanite, omphacite and corundum have also been identified implying that rutile grows continuously during prograde burial and traps mineralogic evidence of this evolution. Pressure estimates obtained from mineral inclusions, when used in conjunction with Zr-in-rutile thermometry, can provide additional constraints on the metamorphic conditions of the host rock. This study demonstrates that rutile is an excellent repository for high- and ultra-high pressure minerals and that the study of mineral inclusions in rutile may profoundly change the way we investigate and recover evidence of such events in both detrital populations and partially retrogressed samples.
Abstract: The theory of plate tectonics describes how the surface of Earth is split into an organized jigsaw of seven large plates1 of similar sizes and a population of smaller plates whose areas follow a fractal distribution2, 3. The reconstruction of global tectonics during the past 200 million years4 suggests that this layout is probably a long-term feature of Earth, but the forces governing it are unknown. Previous studies3, 5, 6, primarily based on the statistical properties of plate distributions, were unable to resolve how the size of the plates is determined by the properties of the lithosphere and the underlying mantle convection. Here we demonstrate that the plate layout of Earth is produced by a dynamic feedback between mantle convection and the strength of the lithosphere. Using three-dimensional spherical models of mantle convection that self-consistently produce the plate size -frequency distribution observed for Earth, we show that subduction geometry drives the tectonic fragmentation that generates plates. The spacing between the slabs controls the layout of large plates, and the stresses caused by the bending of trenches break plates into smaller fragments. Our results explain why the fast evolution in small back-arc plates7, 8 reflects the marked changes in plate motions during times of major reorganizations. Our study opens the way to using convection simulations with plate-like behaviour to unravel how global tectonics and mantle convection are dynamically connected.
Geochimica et Cosmochimica Acta, Vol. 186, pp. 207-225.
Mantle
Subduction, melting
Abstract: The subducted continental crust material will be gravitationally trapped in the deep mantle after having been transported to depths of greater than ?250 -300 km (the “depth of no return”). However, little is known about the status of this trapped continental material as well as its contribution to the mantle heterogeneity after achieving thermal equilibrium with the surrounding mantle. Here, we conduct an experimental study over pressure and temperature ranges of 9 -16 GPa and 1300 -1800 °C to constrain the fate of these trapped upper continental crust (UCC). The experimental results show that partial melting will occur in the subducted UCC along normal mantle geotherm to produce K-rich melt. The residual phases composed of coesite/stishovite + clinopyroxene + kyanite in the upper mantle, and stishovite + clinopyroxene + K-hollandite + garnet + CAS-phase in the mantle transition zone (MTZ), respectively. The residual phases achieve densities greater than the surrounding mantle, which provides a driving force for descent across the 410-km seismic discontinuity into the MTZ. However, this density relationship is reversed at the base of the MTZ, leaving the descended residues to be accumulated above the 660-km seismic discontinuity and may contribute to the “second continent”. The melt is ?0.6 -0.7 g/cm3 less dense than the surrounding mantle, which provides a buoyancy force for ascent of melt to shallow depths. The ascending melt, which preserves a significant portion of the bulk-rock rare earth elements (REEs), large ion lithophile elements (LILEs), and high-filed strength elements (HFSEs), may react with the surrounding mantle. Re-melting of the metasomatized mantle may contribute to the origin of the “enriched mantle sources” (EM-sources). Therefore, the deep subducted continental crust may create geochemical/geophysical heterogeneity in Earth’s interior through subduction, stagnation, partial melting and melt segregation.
Abstract: Two areas with different types of hydration (serpentinization), which occurred in two settings distinct in temperatures, pressures, and stresses, are spatially individualized in the ophiolitic ultramafic massifs of the Polar Urals. The high-temperature hydration of ultramafic rocks occurred in the lithosphere of the mantle wedge directly above the subducted slab. The initial conditions of hydration are limited to 1.2-2 GPa and 650-700°C; a stable assemblage of olivine + antigorite + magnetite ? amphibole ? talc ? chlorite was formed at 0.9-1.2 GPa and 550-600°C. The low-temperature mesh lizardite-chrysotile serpentinization occurred in the crustal, near-surface conditions. Both types of hydration were accompanied by release of hydrogen, which participates in abiogenic CH4 synthesis in the presence of CO2 dissolved in water.
Abstract: Convergent margins, being the boundaries between colliding lithospheric plates, form the most disastrous areas in the world due to intensive, strong seismicity and volcanism. We review global geophysical data in order to illustrate the effects of the plate tectonic processes at convergent margins on the crustal and upper mantle structure, seismicity, and geometry of subducting slab. We present global maps of free-air and Bouguer gravity anomalies, heat flow, seismicity, seismic Vs anomalies in the upper mantle, and plate convergence rate, as well as 20 profiles across different convergent margins. A global analysis of these data for three types of convergent margins, formed by ocean-ocean, ocean-continent, and continent-continent collisions, allows us to recognize the following patterns. (1) Plate convergence rate depends on the type of convergent margins and it is significantly larger when, at least, one of the plates is oceanic. However, the oldest oceanic plate in the Pacific ocean has the smallest convergence rate. (2) The presence of an oceanic plate is, in general, required for generation of high-magnitude (M > 8.0) earthquakes and for generating intermediate and deep seismicity along the convergent margins. When oceanic slabs subduct beneath a continent, a gap in the seismogenic zone exists at depths between ca. 250 km and 500 km. Given that the seismogenic zone terminates at ca. 200 km depth in case of continent-continent collision, we propose oceanic origin of subducting slabs beneath the Zagros, the Pamir, and the Vrancea zone. (3) Dip angle of the subducting slab in continent-ocean collision does not correlate neither with the age of subducting oceanic slab, nor with the convergence rate. For ocean-ocean subduction, clear trends are recognized: steeply dipping slabs are characteristic of young subducting plates and of oceanic plates with high convergence rate, with slab rotation towards a near-vertical dip angle at depths below ca. 500 km at very high convergence rate. (4) Local isostasy is not satisfied at the convergent margins as evidenced by strong free air gravity anomalies of positive and negative signs. However, near-isostatic equilibrium may exist in broad zones of distributed deformation such as Tibet. (5) No systematic patterns are recognized in heat flow data due to strong heterogeneity of measured values which are strongly affected by hydrothermal circulation, magmatic activity, crustal faulting, horizontal heat transfer, and also due to low number of heat flow measurements across many margins. (6) Low upper mantle Vs seismic velocities beneath the convergent margins are restricted to the upper 150 km and may be related to mantle wedge melting which is confined to shallow mantle levels.
Journal of Geophysical Research,, Vol. 121, 7, pp. 5013-5030.
Canada
Subduction
Abstract: The cratonic cores of the continents are remarkably stable and long-lived features. Their ability to resist destructive tectonic processes is associated with their thick (?250 km), cold, chemically depleted, buoyant lithospheric keels that isolate the cratons from the convecting mantle. The formation mechanism and tectonic stability of cratonic keels remains under debate. To address this issue, we use P wave and S wave relative arrival-time tomography to constrain upper mantle structure beneath southeast Canada and the northeast USA, a region spanning three quarters of Earth's geological history. Our models show three distinct, broad zones: Seismic wave speeds increase systematically from the Phanerozoic coastal domains, through the Proterozoic Grenville Province, and to the Archean Superior craton in central Québec. We also recover the NW-SE trending track of the Great Meteor hot spot that crosscuts the major tectonic domains. The decrease in seismic wave speed from Archean to Proterozoic domains across the Grenville Front is consistent with predictions from models of two-stage keel formation, supporting the idea that keel growth may not have been restricted to Archean times. However, while crustal structure studies suggest that Archean Superior material underlies Grenvillian age rocks up to ?300 km SE of the Grenville Front, our tomographic models show a near-vertical boundary in mantle wave speed directly beneath the Grenville Front. We interpret this as evidence for subduction-driven metasomatic enrichment of the Laurentian cratonic margin, prior to keel stabilization. Variable chemical depletion levels across Archean-Proterozoic boundaries worldwide may thus be better explained by metasomatic enrichment than inherently less depleted Proterozoic composition at formation.
Abstract: Subduction of heterogeneous lithologies (sediments and altered basalts) carries a mixture of volatile components (H2O ± CO2) into the mantle, which are later mobilized during episodes of devolatilization and flux melting. Several petrologic and thermodynamic studies investigated CO2 decarbonation to better understand carbon cycling at convergent margins. A paradox arose when investigations showed little to no decarbonation along present day subduction geotherms at subarc depths despite field based observations. Sediment diapirism is invoked as one of several methods for carbon transfer from the subducting slab. We employ high-resolution 2D petrological-thermomechanical modeling to elucidate the role subduction dynamics has with respect to slab decarbonation and the sediment diapirism hypothesis. Our thermodynamic database is modified to account for H2O-CO2 binary fluids via the following lithologies: GLOSS average sediments (H2O: 7.29 wt.% & CO2: 3.01 wt.%), carbonated altered basalts (H2O: 2.63 wt.% & CO2: 2.90 wt.%), and carbonated peridotites (H2O: 1.98 wt.% & CO2: 1.50 wt.%). We include a CO2 solubility P-x[H2O wt.%] parameterization for sediment melts. We parameterize our model by varying two components: slab age (20, 40, 60, 80 Ma) and convergence velocity (1, 2, 3, 4, 5, 6 cm year? 1). 59 numerical models were run and show excellent agreement with the original code base. Three geodynamic regimes showed significant decarbonation. 1) Sedimentary diapirism acts as an efficient physical mechanism for CO2 removal from the slab as it advects into the hotter mantle wedge. 2) If subduction rates are slow, frictional coupling between the subducting and overriding plate occurs. Mafic crust is mechanically incorporated into a section of the lower crust and undergoes decarbonation. 3) During extension and slab rollback, interaction between hot asthenosphere and sediments at shallow depths result in a small window (~ 12.5 Ma) of high integrated CO2 fluxes (205 kg m? 3 Ma? 1).
Abstract: The growth and recycling of continental crust has resulted in the chemical and thermal modification of Earth's mantle, hydrosphere, atmosphere, and biosphere for ?4.0 b.y. However, knowledge of the protolith that gave rise to the first continents and whether the environment of formation was a subduction zone still remains unknown. Here, tonalite melts are formed in high P-T experiments in which primitive oceanic plateau starting material is used as an analogue for Eoarchean (3.6-4.0 Ga) oceanic crust generated at early spreading centers. The tonalites are produced at 1.6-2.2 GPa and 900-950 °C and are mixed with slab-derived aqueous fluids to generate melts that have compositions identical to that of Eoarchean continental crust. Our data support the idea that the first continents formed at ca. 4 Ga and subsequently, through the subduction and partial melting of ?30-45-km-thick Eoarchean oceanic crust, modified Earth's mantle and Eoarchean environments and ecosystems.
Earth and Planetary Science Letters, Vol. 451, pp. 114-124.
Mantle
Subduction
Abstract: The elastic properties of two single crystals of majoritic garnet (Mg3.24Al1.53Si3.23O12 and Mg3.01Fe0.17Al1.68Si3.15O12), have been measured using simultaneously single-crystal X-ray diffraction and Brillouin spectroscopy in an externally heated diamond anvil cell with Ne as pressure transmitting medium at conditions up to ?30 GPa and ?600 K. This combination of techniques makes it possible to use the bulk modulus and unit-cell volume at each condition to calculate the absolute pressure, independently of secondary pressure calibrants. Substitution of the majorite component into pyrope garnet lowers both the bulk (KsKs) and shear modulus (G ). The substitution of Fe was found to cause a small but resolvable increase in KsKs that was accompanied by a decrease in ?Ks/?P?Ks/?P, the first pressure derivative of the bulk modulus. Fe substitution had no influence on either the shear modulus or its pressure derivative. The obtained elasticity data were used to derive a thermo-elastic model to describe VsVs and VpVp of complex garnet solid solutions. Using further elasticity data from the literature and thermodynamic models for mantle phase relations, velocities for mafic, harzburgitic and lherzolitic bulk compositions at the base of Earth's transition zone were calculated. The results show that VsVs predicted by seismic reference models are faster than those calculated for all three types of lithologies along a typical mantle adiabat within the bottom 150 km of the transition zone. The anomalously fast seismic shear velocities might be explained if laterally extensive sections of subducted harzburgite-rich slabs pile up at the base of the transition zone and lower average mantle temperatures within this depth range.
Earth and Planetary Science Letters, Vol. 441, pp. 200-210.
Mantle
Subduction
Abstract: Subduction zones are one of the most striking feature on Earth. They represent one of the two types of convergent plate boundaries, in which one tectonic plates sinks underneath another into the Earth’s mantle. Soon after the advent of the theory of plate tectonics scientists recognized that subduction zones are one of the main drivers of plate motion and mantle convection [Elsasser, 1971]. With trench motion during progressive subduction, overriding plates incorporated in subduction zones may follow the trench and/or deform internally. Such deformation is often characterized by backarc extension, which leads to opening of backarc basins, such as the Tyrrhenian Sea, the Scotia Sea, the Aegean Sea, the North Fiji Basin, and the Lau Basin.
Abstract: The wealth of seismic observations collected over the past 20 years has raised intriguing questions about the three-dimensional (3D) nature of the mantle flow field close to subduction zones and provided a valuable constraint for how the plate geometry may influence mantle flow proximal to the slab. In geodynamics, there has been a new direction of subduction zone modelling that has explored the 3D nature of slab-driven mantle flow, motivated in part by the observations from shear wave splitting, but also by the observed variations in slab geometries worldwide. Advances in high-performance computing are now allowing for an unprecedented level of detail to be incorporated into numerical models of subduction. This paper summarizes recent advances from 3D geodynamic models that reveal the complex nature of slab-driven mantle flow, including trench parallel flow, toroidal flow around slab edges, mantle upwelling at lateral slab edges, and small scale convection within the mantle wedge. This implies slab-driven mantle deformation zones occur in the asthenosphere proximal to the slab, wherein the mantle may commonly flow in a different direction and rate than the surface plates, implying laterally variable plate-mantle coupling. The 3D slab-driven mantle flow can explain, in part, the lateral transport of geochemical signatures in subduction zones. In addition, high-resolution geographically referenced models can inform the interpretation of slab structure, where seismic data are lacking. The incorporation of complex plate boundaries into high-resolution, 3D numerical models opens the door to a new avenue of research in model construction, data assimilation, and modelling workflows, and gives 3D immersive visualization a new role in scientific discovery.
Abstract: Geological observations indicate that there are only a few rocks of Archean Earth and no Hadean rocks on the surface of the present-day Earth. From these facts, many scientists believe that the primordial continents never existed during Hadean Earth, and the continental volume has kept increasing. On the other hand, recent studies reported the importance of the primordial continents on the origin of life, implying their existence. In this paper, we discussed the possible process that could explain the loss of the primordial continents with the assumption that they existed in the Hadean. Although depending on the timing of the initiation of plate tectonics and its convection style, subduction erosion, which is observed on the present-day Earth, might have carried the primordial continents into the deep mantle.
Abstract: Numerical experiments are used to investigate the thermo-mechanical controls for inducing flat subduction and why flat subduction is rare relative to normal/steep subduction. Our modeling results demonstrate that flat subduction is an end-member of a steady state subduction geometry and is characterized by a curved slab with a nearly-horizontal slab section. Intermediate cases between normal/steep and flat subduction appear to be transient in origin and evolve toward one of the stable end-members. Physical parameters inducing flat subduction can be classified into four categories: buoyancy of the subducting oceanic lithosphere (e.g., slab age, oceanic crustal thickness), viscous coupling between the overriding and downgoing plates (e.g., initial subduction angle), external kinematic conditions, and rheological properties of the subduction zone. On the basis of parameter sensitivity tests and the main characteristics of present-day flat subduction zones, positive buoyancy from either the young slab or the thickened oceanic crust is considered as the primary controlling parameter. Our results show that the possibility of flat subduction is directly proportional to oceanic crustal thickness and inversely proportional to the slab age. Furthermore, oceanic crust must be thicker than 8 km to induce flat subduction, when the slab is older than 30 Ma with an initial subduction angle of ? 20° and without absolute trenchward motion of the overriding plate. The lower the initial subduction angle or the thicker the overriding continental lithosphere, the more likelihood for flat subduction. The initial subduction angle is more influential for the development of flat subduction than the overriding lithospheric thickness, and a thick overriding lithosphere induces flat subduction only under the condition of an initial subduction angle of ? 25°, with a slab age of ? 30 Ma and without absolute trenchward motion of the overriding plate. However, when the initial subduction angle is increased to > 25°, no flat subduction is predicted. All the parameters are evaluated within the constraints of a mechanical framework in which the slab geometry is regarded as a result of a balance between the gravitational and hydrodynamic torques. Any factor that can sufficiently reduce gravitational torque or increase hydrodynamic torque will exert a strong effect on flat subduction development. Our results are consistent with the observations of modern flat subduction zones on Earth.
Abstract: The primordial crust on the Earth formed from the crystallization of the surface magma ocean during the Hadean. However, geological surveys have found no evidence of rocks dating back to more than 4 Ga on the Earth's surface, suggesting the Hadean crust was lost due to some processes. We investigated the subduction of one of the possible candidates for the primordial crust, anorthosite and KREEP crust similar to the Moon, which is also considered to have formed from the crystallization of the magma ocean. Similar to the present Earth, the subduction of primordial crust by subduction erosion is expected to be an effective way of eliminating primordial crust from the surface. In this study, the subduction rate of the primordial crust via subduction channels is evaluated by numerical simulations. The subduction channels are located between the subducting slab and the mantle wedge and are comprised of primordial crust materials supplied mainly by subduction erosion. We have found that primordial anorthosite and KREEP crust of up to ?50 km thick at the Earth's surface was able to be conveyed to the deep mantle within 0.1-2 Gy by that mechanism.
Abstract: Continental crust is buoyant compared with its oceanic counterpart and resists subduction into the mantle. When two continents collide, the mass balance for the continental crust is therefore assumed to be maintained. Here we use estimates of pre-collisional crustal thickness and convergence history derived from plate kinematic models to calculate the crustal mass balance in the India-Asia collisional system. Using the current best estimates for the timing of the diachronous onset of collision between India and Eurasia, we find that about 50% of the pre-collisional continental crustal mass cannot be accounted for in the crustal reservoir preserved at Earth’s surface today—represented by the mass preserved in the thickened crust that makes up the Himalaya, Tibet and much of adjacent Asia, as well as southeast Asian tectonic escape and exported eroded sediments. This implies large-scale subduction of continental crust during the collision, with a mass equivalent to about 15% of the total oceanic crustal subduction flux since 56 million years ago. We suggest that similar contamination of the mantle by direct input of radiogenic continental crustal materials during past continent-continent collisions is reflected in some ocean crust and ocean island basalt geochemistry. The subduction of continental crust may therefore contribute significantly to the evolution of mantle geochemistry.
Contributions to Mineralogy and Petrology, Vol. 171, pp. 107-
Mantle
Subduction
Abstract: This experimental study simulates the interaction of hotter, deeper hydrous mantle melts with shallower, cooler depleted mantle, a process that is expected to occur in the upper part of the mantle wedge. Hydrous reaction experiments (~6 wt% H2O in the melt) were conducted on three different ratios of a 1.6 GPa mantle melt and an overlying 1.2 GPa harzburgite from 1060 to 1260 °C. Reaction coefficients were calculated for each experiment to determine the effect of temperature and starting bulk composition on final melt compositions and crystallizing assemblages. The experiments used to construct the melt-wall rock model closely approached equilibrium and experienced <5% Fe loss or gain. Experiments that experienced higher extents of Fe loss were used to critically evaluate the practice of “correcting” for Fe loss by adding iron. At low ratios of melt/mantle (20:80 and 5:95), the crystallizing assemblages are dunites, harzburgites, and lherzolites (as a function of temperature). When the ratio of deeper melt to overlying mantle is 70:30, the crystallizing assemblage is a wehrlite. This shows that wehrlites, which are observed in ophiolites and mantle xenoliths, can be formed by large amounts of deeper melt fluxing though the mantle wedge during ascent. In all cases, orthopyroxene dissolves in the melt, and olivine crystallizes along with pyroxenes and spinel. The amount of reaction between deeper melts and overlying mantle, simulated here by the three starting compositions, imposes a strong influence on final melt compositions, particularly in terms of depletion. At the lowest melt/mantle ratios, the resulting melt is an extremely depleted Al-poor, high-Si andesite. As the fraction of melt to mantle increases, final melts resemble primitive basaltic andesites found in arcs globally. An important element ratio in mantle lherzolite composition, the Ca/Al ratio, can be significantly elevated through shallow mantle melt-wall rock reaction. Wall rock temperature is a key variable; over a span of <80 °C, reaction with deeper melt creates the entire range of mantle lithologies from a depleted dunite to a harzburgite to a refertilized lherzolite. Together, the experimental phase equilibria, melt compositions, and reaction coefficients provide a framework for understanding how melt-wall rock reaction occurs in the natural system during melt ascent in the mantle wedge.
Abstract: Oceanic lithospheric mantle is generally interpreted as depleted mantle residue after mid-ocean ridge basalt extraction. Several models have suggested that metasomatic processes can refertilize portions of the lithospheric mantle before subduction. Here, we report mantle xenocrysts and xenoliths in petit-spot lavas that provide direct evidence that the lower oceanic lithosphere is affected by metasomatic processes. We find a chemical similarity between clinopyroxene observed in petit-spot mantle xenoliths and clinopyroxene from melt-metasomatized garnet or spinel peridotites, which are sampled by kimberlites and intracontinental basalts respectively. We suggest that extensional stresses in oceanic lithosphere, such as plate bending in front of subduction zones, allow low-degree melts from the seismic low-velocity zone to percolate, interact and weaken the oceanic lithospheric mantle. Thus, metasomatism is not limited to mantle upwelling zones such as mid-ocean ridges or mantle plumes, but could be initiated by tectonic processes. Since plate flexure is a global mechanism in subduction zones, a significant portion of oceanic lithospheric mantle is likely to be metasomatized. Recycling of metasomatic domains into the convecting mantle is fundamental to understanding the generation of small-scale mantle isotopic and volatile heterogeneities sampled by oceanic island and mid-ocean ridge basalts.
Geochimica et Cosmochimica Acta, Vol. 198, pp. 70-91.
Mantle
Subduction
Abstract: The equilibrium between aqueous fluids and allanite-bearing eclogite has been investigated to constrain the effect of temperature (T) and fluid composition on the stability of allanite and on the mobility of major and trace elements during the dehydration of eclogites. The experiments were performed at 590-800 °C and 2.4-2.6 GPa, and fluids were sampled as synthetic fluid inclusions in quartz using an improved entrapment technique. The concentrations and bulk partition coefficients were determined for a range of major (Mg, Ca, Na, Fe, Al, Ti) and 16 trace elements as a function of T and fluid composition. The results reveal a significant effect of T on element partitioning between the fluids and the solid mineral assemblage. The partition coefficients increase by more than an order of magnitude for most of the major and trace elements, and several orders of magnitude for light rare-earth elements (LREE) from 590 to 800 °C. The addition of various ligand species into the fluid at 700 °C results in distinctive trends on element partitioning. The concentrations and corresponding partition coefficients of most of the elements are enhanced upon addition of NaF to the fluid. In contrast, NaCl displays a nearly opposite effect by suppressing the solubilities of major elements and consequently affecting the mobility of trace elements that form stable complexes with alkali-(alumino)-silicate clusters in the fluid, e.g. high field strength elements (HFSE). The results further suggest that fluids in equilibrium with orthopyroxene and/or diopsidic clinopyroxene are peralkaline (ASI ?0.1-0.7), whereas fluids in equilibrium with omphacitic pyroxene are more peraluminous (ASI ?1.15). Therefore, natural aqueous fluids in equilibrium with eclogite at about 90 km depth will be slightly peraluminous in composition. Another important finding of this study is the relatively high capacity of aqueous fluids to mobilize LREE, which may be even higher than that of hydrous melts.
Abstract: When deeply buried in subduction zones, rocks undergo mineral transformations that record the increase of pressure and temperature. The fact that high-pressure metamorphic parageneses are found at the Earth’s surface proves that rock burial is followed by exhumation. Here we use analysis of available data sets from high-pressure metamorphic rocks worldwide to show that the peak pressure is proportional to the subsequent decompression occurring during the initial stage of retrogression. We propose, using a simple mechanical analysis, that this linear relationship can be explained by the transition from burial-related compression to extension at the onset of exhumation. This major switch in orientation and magnitude of principal tectonic stresses leads to a catastrophic pressure drop prior to actual rock ascent. Therefore, peak pressures are not necessarily, as commonly believed, directly dependent on the maximum burial depth, but can also reflect a change of tectonic regime. Our results, which are in agreement with natural data, have significant implications for rock rheology, subduction zone seismicity, and the magnitudes of tectonic pressures sustained by rocks. Current views of subduction dynamics could be reconsidered in that perspective.
Science China Earth Sciences, Vol. 60, 2, pp. 207-217.
Technology
Subduction
Abstract: Diamonds are renowned as the record of Earth’s evolution history. Natural diamonds on the Earth can be distinguished in light of genetic types as kimberlitic diamonds (including peridotitic diamonds and eclogitic diamonds), ultrahigh-pressure metamorphic diamonds and ophiolitic diamonds. According to the inclusion mineralogy, most diamonds originated from continental lithospheric mantle at depths of 140-250 km. Several localities, however, yield ultradeep diamonds with inclusion compositions that require a sublithospheric origin (>~250 km). Ultradeep diamonds exhibit distinctions in terms of carbon isotope composition, N-concentration, mineral inclusions and so on. The present study provides a systematic compilation concerning the features of ultradeep diamonds, based on which to expound their genesis affinity with mantle-carbonate melts. The diamond-parental carbonate melts are proposed to be stemmed from the Earth’s crust through subduction of oceanic lithosphere. Ultradeep diamonds are classified into a subgroup attaching to kimberlitic diamonds grounded by formation mechanism, and present connections in respect of carbon origin to eclogitic diamonds, ultrahigh-pressure metamorphic diamonds and ophiolitic diamonds.
Abstract: There exist three mainstream opinions regarding the timing of the initial collision between the Indian and Eurasian continents, namely, 65±5, 45±5, and 30±5 Ma. Five criteria are proposed for determining which tectonic event was related to the initial collision between India and Asia: the rapid decrease in the rate of plate motion, the cessation of magmatic activity originating from the subduction of oceanic crust, the end of sedimentation of oceanic facies, the occurrence of intracontinental deformation, and the exchange of sediments sourced from two continents. These criteria are used to constrain the nature of these tectonic events. It is proposed that the 65±5 Ma tectonic event is consistent with some of the criteria, but the upshot of this model is that the magmatic activity originating from the Tethyan subduction since the Mesozoic restarted along the southern margin of the Asian continent in this time after a brief calm, implying that the subduction of the Neotethys slab was still taking place. The magmatic activity that occurred along the southern margin of the Asian continent had a 7-Myr break during 72-65 Ma, which in this study is interpreted as having resulted from tectonic transformation from subduction to transform faulting, indicating that the convergence between the Indian and Asian continents was once dominated by strike-slip motion. The 30±5 Ma tectonic event resulted in the uplift of the Tibetan Plateau, which was related to the late stage of the convergence between these two continents, namely, a hard collision. The 45±5 Ma tectonic event is in accordance with most of the criteria, corresponding to the initial collision between these two continents.
Earth and Planetary Science Letters, Vol. 464, pp. 10-23.
Mantle
Subduction
Abstract: Seismic tomography shows that subducting slabs can either sink straight into the lower mantle, or lie down in the mantle transition zone. Moreover, some slabs seem to have changed mode from stagnation to penetration or vice versa. We investigate the dynamic controls on these modes and particularly the transition between them using 2D self-consistent thermo-mechanical subduction models. Our models confirm that the ability of the trench to move is key for slab flattening in the transition zone. Over a wide range of plausible Clapeyron slopes and viscosity jumps at the base of the transition zone, hot young slabs (25 Myr in our models) are most likely to penetrate, while cold old slabs (150 Myr) drive more trench motion and tend to stagnate. Several mechanisms are able to induce penetrating slabs to stagnate: ageing of the subducting plate, decreasing upper plate forcing, and increasing Clapeyron slope (e.g. due to the arrival of a more hydrated slab). Getting stagnating slabs to penetrate is more difficult. It can be accomplished by an instantaneous change in the forcing of the upper plate from free to motionless, or a sudden decrease in the Clapeyron slope. A rapid change in plate age at the trench from old to young cannot easily induce penetration. On Earth, ageing of the subducting plate (with accompanying upper plate rifting) may be the most common mechanism for causing slab stagnation, while strong changes in upper plate forcing appear required for triggering slab penetration.
Earth and Planetary Science Letters, Vol. 465, pp. 155-167.
Mantle
Subduction
Abstract: Subduction is integral to mantle convection and plate tectonics, yet the role of the subslab mantle in this process is poorly understood. Some propose that decoupling from the slab permits widespread trench parallel flow in the subslab mantle, although the geodynamical feasibility of this has been questioned. Here, we use the source-side shear wave splitting technique to probe anisotropy beneath subducting slabs, enabling us to test petrofabric models and constrain the geometry of mantle fow. Our global dataset contains 6369 high quality measurements - spanning ?40,000?km?40,000?km of subduction zone trenches - over the complete range of available source depths (4 to 687?km) - and a large range of angles in the slab reference frame. We find that anisotropy in the subslab mantle is well characterised by tilted transverse isotropy with a slow-symmetry-axis pointing normal to the plane of the slab. This appears incompatible with purely trench-parallel flow models. On the other hand it is compatible with the idea that the asthenosphere is tilted and entrained during subduction. Trench parallel measurements are most commonly associated with shallow events (source depth <50?km<50?km) - suggesting a separate region of anisotropy in the lithospheric slab. This may correspond to the shape preferred orientation of cracks, fractures, and faults opened by slab bending. Meanwhile the deepest events probe the upper lower mantle where splitting is found to be consistent with deformed bridgmanite.
Abstract: Many carbonatites are associated both spatially and temporally with large igneous provinces (LIPs), and considered to originate from a mantle plume source lacking any contribution from recycled crustal materials. Here, we report an occurrence of carbonatite enriched in rare-earth elements (REE) and associated with the Tarim LIP in northwestern China. The Tarim LIP comprises intrusive and volcanic products of mantle plume activity spanning from ~ 300 to 280 Ma. The carbonatites at Wajilitage in the northwestern part of Tarim are dominated by calcite and dolomite varieties, and contain abundant REE minerals (principally, monazite and REE-fluorcarbonates). Th-Pb age determination of monazite yielded an emplacement age of 266 ± 5.3 Ma, i.e. appreciably younger than the eruption age of flood basalts at ~ 290 Ma. The carbonatites show low initial 87Sr/86Sr (0.7037-0.7041) and high ?Nd(t) (1.2-4) values, which depart from the isotopic characteristics of plume-derived basalts and high-Mg picrites from the same area. This indicates that the Wajilitage carbonatites derived from a mantle source isotopically distinct from the one responsible for the voluminous (ultra)mafic volcanism at Tarim. The carbonatites show ?26MgDSM3 values (? 0.99 to ? 0.65‰) that are significantly lower than those in typical mantle-derived rocks and rift carbonatites, but close to marine sediments and orogenic carbonatites. We propose that the carbonatites in the Tarim LIP formed by decompressional melting of recycled sediments mixed with the ambient mantle peridotite. The enriched components in the Tarim plume could be accounted for by the presence of recycled sedimentary components in the subcontinental mantle.
Abstract: Some of Earth's coldest mantle is found in subduction zones at the tip of the mantle wedge that lies between the subducting and overriding plates. This forearc mantle is isolated from the flow of hot material beneath the volcanic arc, and so is inferred to reach temperatures no more than 600 to 800 °C — conditions at which hydrous mantle minerals should be stable. The forearc mantle could therefore constitute a significant reservoir for water if sufficient water is released from the subducting slab into the mantle wedge. Such a reservoir could hydrate the plate interface and has been invoked to aid the genesis of megathrust earthquakes and slow slip events. Our synthesis of results from thermal models that simulate the conditions for subduction zones globally, however, indicates that dehydration of subducting plates is too slow over the life span of a typical subduction zone to hydrate the forearc mantle. Hot subduction zones, where slabs dehydrate rapidly, are an exception. The hottest, most buoyant forearcs are most likely to survive plate collisions and be exhumed to the surface, so probably dominate the metamorphic rock record. Analysis of global seismic data confirms the generally dry nature of mantle forearcs. We conclude that many subduction zones probably liberate insufficient water to hydrate the shallower plate boundary where great earthquakes and slow slip events nucleate. Thus, we suggest that it is solid-state processes and not hydration that leads to weakening of the plate interface in cold subduction zones.
Abstract: Diamonds are renowned as the record of Earth’s evolution history. Natural diamonds on the Earth can be distinguished in light of genetic types as kimberlitic diamonds (including peridotitic diamonds and eclogitic diamonds), ultrahigh-pressure metamorphic diamonds and ophiolitic diamonds. According to the inclusion mineralogy, most diamonds originated from continental lithospheric mantle at depths of 140-250 km. Several localities, however, yield ultradeep diamonds with inclusion compositions that require a sublithospheric origin (>~250 km). Ultradeep diamonds exhibit distinctions in terms of carbon isotope composition, N-concentration, mineral inclusions and so on. The present study provides a systematic compilation concerning the features of ultradeep diamonds, based on which to expound their genesis affinity with mantle-carbonate melts. The diamond-parental carbonate melts are proposed to be stemmed from the Earth’s crust through subduction of oceanic lithosphere. Ultradeep diamonds are classified into a subgroup attaching to kimberlitic diamonds grounded by formation mechanism, and present connections in respect of carbon origin to eclogitic diamonds, ultrahigh-pressure metamorphic diamonds and ophiolitic diamonds.
Abstract: Some of Earth's coldest mantle is found in subduction zones at the tip of the mantle wedge that lies between the subducting and overriding plates. This forearc mantle is isolated from the flow of hot material beneath the volcanic arc, and so is inferred to reach temperatures no more than 600 to 800 °C - conditions at which hydrous mantle minerals should be stable. The forearc mantle could therefore constitute a significant reservoir for water if sufficient water is released from the subducting slab into the mantle wedge. Such a reservoir could hydrate the plate interface and has been invoked to aid the genesis of megathrust earthquakes and slow slip events. Our synthesis of results from thermal models that simulate the conditions for subduction zones globally, however, indicates that dehydration of subducting plates is too slow over the life span of a typical subduction zone to hydrate the forearc mantle. Hot subduction zones, where slabs dehydrate rapidly, are an exception. The hottest, most buoyant forearcs are most likely to survive plate collisions and be exhumed to the surface, so probably dominate the metamorphic rock record. Analysis of global seismic data confirms the generally dry nature of mantle forearcs. We conclude that many subduction zones probably liberate insufficient water to hydrate the shallower plate boundary where great earthquakes and slow slip events nucleate. Thus, we suggest that it is solid-state processes and not hydration that leads to weakening of the plate interface in cold subduction zones.
Abstract: The lower plate is the dominant agent in modern convergent margins characterized by active subduction, as negatively buoyant oceanic lithosphere sinks into the asthenosphere under its own weight. This is a strong plate-driving force because the slab-pull force is transmitted through the stiff sub-oceanic lithospheric mantle. As geological and geochemical data seem inconsistent with the existence of modern-style ridges and arcs in the Archaean, a periodically-destabilized stagnant-lid crust system is proposed instead. Stagnant-lid intervals may correspond to periods of layered mantle convection where efficient cooling was restricted to the upper mantle, perturbing Earth's heat generation/loss balance, eventually triggering mantle overturns. Archaean basalts were derived from fertile mantle in overturn upwelling zones (OUZOs), which were larger and longer-lived than post-Archaean plumes. Early cratons/continents probably formed above OUZOs as large volumes of basalt and komatiite were delivered for protracted periods, allowing basal crustal cannibalism, garnetiferous crustal restite delamination, and coupled development of continental crust and sub-continental lithospheric mantle. Periodic mixing and rehomogenization during overturns retarded development of isotopically depleted MORB (mid-ocean ridge basalt) mantle. Only after the start of true subduction did sequestration of subducted slabs at the core-mantle boundary lead to the development of the depleted MORB mantle source. During Archaean mantle overturns, pre-existing continents located above OUZOs would be strongly reworked; whereas OUZO-distal continents would drift in response to mantle currents. The leading edge of drifting Archaean continents would be convergent margins characterized by terrane accretion, imbrication, subcretion and anatexis of unsubductable oceanic lithosphere. As Earth cooled and the background oceanic lithosphere became denser and stiffer, there would be an increasing probability that oceanic crustal segments could founder in an organized way, producing a gradual evolution of pre-subduction convergent margins into modern-style active subduction systems around 2.5 Ga. Plate tectonics today is constituted of: (1) a continental drift system that started in the Early Archaean, driven by deep mantle currents pressing against the Archaean-age sub-continental lithospheric mantle keels that underlie Archaean cratons; (2) a subduction-driven system that started near the end of the Archaean.
Abstract: The breakup of continents and their subsequent drifting plays a crucial role in the Earth's periodic plate aggregation and dispersal cycles. While continental aggregation is considered the result of oceanic closure during subduction, what drives sustained divergence in the following stages remains poorly understood. In this study, thermo-mechanical numerical experiments illustrate the single contribution of subduction and coupled mantle flow to the rifting and drifting of continents. We quantify the drag exerted by subduction-induced mantle flow along the basal surface of continental plates, comparing models of lithospheric slab stagnation above the upper-lower mantle boundary with those where slabs penetrate into the lower mantle. When subduction is upper-mantle confined, divergent basal tractions localise at distances comparable to the effective upper mantle thickness (~ 500 km), causing the opening of a marginal basin. Instead, subduction of lithosphere in the lower mantle reorganises the flow into a much wider cell localising extensional stresses at greater distances from the trench (~ 3000 km). Sub-continental tractions are higher and more sustained over longer time periods in this case, and progressively increase as the slab sinks deeper. Although relatively low, basal-shear stresses when integrated over large plates, generate tension forces that may exceed the strength of the continental lithosphere, eventually leading to breakup and opening of a distal basin. The models illustrate the emergence of a similar mechanism, which results in the formation of back-arc basins above upper-mantle confined subduction, and scales to much larger distances for deeper subduction. Examples include the Atlantic Ocean formation and drifting of the South and North American plates during the Mesozoic-Cenozoic Farallon plate subduction.
Abstract: The net flux of carbon between the Earth’s interior and exterior, which is critical for redox evolution and planetary habitability, relies heavily on the extent of carbon subduction. While the fate of carbonates during subduction has been studied, little is known about how organic carbon is transferred from the Earth’s surface to the interior, although organic carbon sequestration is related to sources of oxygen in the surface environment. Here we use high pressure–temperature experiments to determine the capacity of rhyolitic melts to carry carbon under graphite-saturated conditions in a subducting slab, and thus to constrain the subduction efficiency of organic carbon, the remnants of life, through time. We use our experimental data and a thermodynamic model of CO2 dissolution in slab melts to quantify organic carbon mobility as a function of slab parameters. We show that the subduction of graphitized organic carbon, and the graphite and diamond formed by reduction of carbonates with depth, remained efficient even in ancient, hotter subduction zones where oxidized carbon subduction probably remained limited. We suggest that immobilization of organic carbon in subduction zones and deep sequestration in the mantle facilitated the rise (~103–5 fold) and maintenance of atmospheric oxygen since the Palaeoproterozoic and is causally linked to the Great Oxidation Event. Our modelling shows that episodic recycling of organic carbon before the Great Oxidation Event may also explain occasional whiffs of atmospheric oxygen observed in the Archaean.
Abstract: As subducting plates reach the base of the upper mantle, some appear to flatten and stagnate, while others seemingly go through unimpeded. This variable resistance to slab sinking has been proposed to affect long-term thermal and chemical mantle circulation. A review of observational constraints and dynamic models highlights that neither the increase in viscosity between upper and lower mantle (likely by a factor 20–50) nor the coincident endothermic phase transition in the main mantle silicates (with a likely Clapeyron slope of –1 to –2 MPa/K) suffice to stagnate slabs. However, together the two provide enough resistance to temporarily stagnate subducting plates, if they subduct accompanied by significant trench retreat. Older, stronger plates are more capable of inducing trench retreat, explaining why backarc spreading and flat slabs tend to be associated with old-plate subduction. Slab viscosities that are ?2 orders of magnitude higher than background mantle (effective yield stresses of 100–300 MPa) lead to similar styles of deformation as those revealed by seismic tomography and slab earthquakes. None of the current transition-zone slabs seem to have stagnated there more than 60 m.y. Since modeled slab destabilization takes more than 100 m.y., lower-mantle entry is apparently usually triggered (e.g., by changes in plate buoyancy). Many of the complex morphologies of lower-mantle slabs can be the result of sinking and subsequent deformation of originally stagnated slabs, which can retain flat morphologies in the top of the lower mantle, fold as they sink deeper, and eventually form bulky shapes in the deep mantle.
Abstract: The ophiolites that crop out discontinuously along the ?2000 km Yarlung Zangbo Suture zone (YZSZ) between the Nanga Parbat and Namche Barwa syntaxes in southern Tibet represent the remnants of Neotethyan oceanic lithosphere (Fig. 1a). We have investigated the internal structure and the geochemical makeup of mafic-ultramafic rock assemblages that are exposed in the westernmost segment of the YZSZ where the suture zone architecture displays two distinct sub-belts of ophiolitic and mélange units separated by a continental Zhongba terrane (Fig. 1b). These two sub-belts include the Daba – Xiugugabu in the south (Southern sub-belt, SSB) and the Dajiweng – Saga in the north (Northern sub-belt, NSB). We present new structural, geochemical, geochronological data from upper mantle peridotites and mafic dike intrusions occurring in these two sub-belts and discuss their tectonomagmatic origin. In-situ analysis of zircon grains obtained from mafic dikes within the Baer, Cuobuzha and Jianabeng massifs in the NSB, and within the Dongbo, Purang, Xiugugabu, Zhaga and Zhongba in the SSB have yielded crystallization ages ranging between130 and 122 Ma. Dike rocks in both sub-belts show N-MORB REE patterns and negative Nb, Ta and Ti anomalies, reminiscent of those documented from SSZ ophiolites. Harzburgitic host rocks of the mafic dike intrusions mainly display geochemical compositions of abyssal peridotites (Fig. 2), with the exception of the Dajiweng harzburgites, which show the geochemical signatures of forearc peridotites (Lian et al., 2016). Extrusive rocks that are spatially associated with these peridotite massifs in both sub-belts also have varying compositional and geochemical features. Tithonian to Valanginian (150 – 135 Ma) basaltic rocks in the Dongbo massif have OIB-like geochemistry and 138 Ma basaltic lavas in the Purang massif have EMORB-like geochemistry (Liu et al., 2015). Tuffaceous rocks in the Dajiweng massif are 140 Ma in age and show OIB-like geochemistry. We interpret these age and geochemical data to reflect a rifted continental margin origin of the extrusive rock units in both sub-belts. These data and structural observations show that the western Yarluang Zangbo ophiolites represent fragments of an Ocean-Continent Transition (OCT) peridotites altered by fluids in an initial supersubduction setting. We infer that mafic-ultramafic rock assemblages exposed in the SSB and NSB initially formed in an ocean – continent transition zone (OCTZ) during the late Jurassic, and that they were subsequently emplaced in the forearc setting of an intraoceanic subduction zone within a Neotethyan seaway during 130 to 122 Ma. The NSB and SSB are hence part of a single, S-directed nappe sheet derived from a Neotethyan seaway located north of the Zhongba terrane.
Abstract: Fluids are important for mass transfer at the slab–mantle interface in subduction zones. However, it is usually difficult to trace fluids from specific sources in a subducting slab, especially those derived from dehydration of serpentinite. Coesite-bearing whiteschist at Dora-Maira in the Western Alps is characterized by strong Mg enrichment relative to the country rocks, which requires infiltration of Mg-rich fluids into the supracrustal rock. In order to constrain the origin of such Mg-rich fluids, we have performed an integrated study of whole-rock Mg and O isotopes, zircon U–Pb ages and O isotopes for the whiteschist and related rocks. Zircons in the whiteschist show two groups of U–Pb ages at ?262 Ma and ?34 Ma, respectively, for relict and newly grown domains. The Permian U–Pb ages of relict magmatic domains are consistent with the protolith age of host metagranite, suggesting that their common protolith is the Permian granite. The Tertiary U–Pb ages occur in coesite-bearing metamorphic domains, consistent with the known age for ultrahigh-pressure metamorphism. The metamorphic domains have ?18O values of
‰5.8–6.8‰, whereas the relict magmatic domains have high ?18O values of ‰?10‰. Such high ?18O values are also characteristic of the metagranite, indicating that the whiteschist protolith underwent metasomatism by metamorphic fluids with low ?18O value of f ‰?2–4‰. The whiteschist mostly has whole-rock ?26Mg values of ?0.07 to 0.72‰, considerably higher than country-rock ?26Mg values of ?0.54 to ‰?0.11‰. Thus, the metamorphic fluids are not only rich in Mg but also heavy in Mg isotopes. They were probably derived from the breakdown of Mg-rich hydrous minerals such as talc and antigorite in serpentinite at the slab–mantle interface in the subduction channel. Therefore, the dehydration of mantle wedge serpentinite during the subduction and exhumation of continental crust can provide the Mg-rich fluids responsible for the metasomatism of crustal rocks at subarc depths.
Physics of the Earth and Planetary Interiors, Vol. 269, 1, pp. 112-132.
Mantle
subduction
Abstract: The present study aims at better deciphering the different mechanisms involved in the functioning of the subduction interplate. A 2D thermo-mechanical model is used to simulate a subduction channel, made of oceanic crust, free to evolve. Convergence at constant rate is imposed under a 100 km thick upper plate. Pseudo-brittle and non-Newtonian behaviours are modelled. The influence of the subduction channel strength, parameterized by the difference in activation energy between crust and mantle (?Ea) is investigated to examine in detail the variations in depth of the subduction plane down-dip extent, zcoup . First, simulations show that numerical resolution may be responsible for an artificial and significant shallowing of zcoup if the weak crustal layer is not correctly resolved. Second, if the age of the subducting plate is 100 Myr, subduction occurs for any ?Ea . The stiffer the crust is, that is, the lower ?Ea is, the shallower zcoup is (60 km depth if ?Ea = 20 kJ/mol) and the hotter the fore-arc base is. Conversely, imposing a very weak subduction channel (?Ea > 135 J/mol) leads there to an extreme mantle wedge cooling and inhibits mantle melting in wet conditions. Partial kinematic coupling at the fore-arc base occurs if ?Ea = 145 kJ/mol. If the incoming plate is 20 Myr old, subduction can occur under the conditions that the crust is either stiff and denser than the mantle, or weak and buoyant. In the latter condition, cold crust plumes rise from the subduction channel and ascend through the upper lithosphere, triggering (1) partial kinematic coupling under the fore-arc, (2) fore-arc lithosphere cooling, and (3) partial or complete hindrance of wet mantle melting. zcoup then ranges from 50 to more than 250 km depth and is time-dependent if crust plumes form. Finally, subduction plane dynamics is intimately linked to the regime of subduction-induced corner flow. Two different intervals of ?Ea are underlined: 80-120 kJ/mol to reproduce the range of slab surface temperature inferred from geothermometry, and 10-40 kJ/mol to reproduce the shallow hot mantle wedge core inferred from conditions of last equilibration of near-primary arc magmas and seismic tomographies. Therefore, an extra process controlling mantle wedge dynamics is needed to satisfy simultaneously the aforementioned observations. A mantle viscosity reduction, by a factor 4-20, caused by metasomatism in the mantle wedge is proposed. From these results, I conclude that the subduction channel down-dip extent, zcoup , should depend on the subduction setting, to be consistent with the observed variability of sub-arc depths of the subducting plate surface.
Abstract: There is much discussion of the timing of the onset of plate tectonics, yet there is increasing evidence that magma types similar to those from recent within plate and subduction related settings were generated in different areas at broadly similar times in the early Archaean. It may therefore be helpful to consider when plate tectonics became the dominant mechanism associated with the generation of continental crust, rather than just when it started. To do this we discuss the geochemical and mechanical characteristics of the lithosphere in the search to discern differences between Early Archaean and younger tectonic environments. Seismic tomography allows increasingly detailed mapping of the lithosphere, and it provides some evidence that the degree of anisotropy is different in different Archaean terrains. Structural styles also appear to vary from basin and swell, or vertical tectonics, as in the Australia Pilbara and southern Africa, to those with more strongly developed regional fabrics and greater seismic anisotropy, as in North America. These terrains tend to be characterized by inferred within-plate and subduction-related magmatism respectively, and we consider possible links between the degree of crustal and mantle anisotropy and the nature of the magmatic record. At least in some areas, terrains with stronger regional fabrics may be younger than those in which such fabrics are less well developed. A model is developed for the generation and stabilization of continental lithosphere in the Archaean. It seeks to reconcile evidence for hot shallow melting with melt fractions up to 40% to generate residual peridotites now preserved as mantle xenoliths, and the lower degrees of melting required to generate the mafic sources of TTGs (fractionated Lu/Hf and Sm/Nd, and perhaps not Rb/Sr).
Earth and Planetary Science Letters, Vol. 474, pp. 516-526.
Mantle
subduction
Abstract: Higher mantle potential temperatures characterized the early Earth, resulting in thicker, more mafic oceanic crust entering subduction systems. This change in the composition of subducted slabs, combined with the enhanced temperature contrast between the slab and ambient mantle, would have altered the buoyancy forces driving subduction in the early Earth. Here we investigate this “compositional effect” through a combination of petrologic and thermal modeling. Specifically, we construct density profiles for sinking slabs under modern and early Earth conditions based on a range of mafic crust and mantle compositions. Slab and mantle densities are then determined from mineral assemblages calculated using the thermodynamic modeling program Perple_X along slab geotherms estimated from an analytic thermal model. Consistent with previous studies, we find that modern MORB compositions are typically less dense than the ambient mantle in the basalt barrier zone, located immediately beneath the mantle transition zone. By contrast, possible early Earth oceanic crust compositions are denser than ambient mantle at all depths down to 1000 km. This compositional effect results in slabs that would have more readily penetrated the transition zone, promoting single-layered convection and effective mantle mixing in the early Earth.
Earth and Planetary Science Letters, Vol. 474, pp. 59-67.
Mantle, India-Eurasia
subduction
Abstract: We present three-dimensional numerical models to investigate the dynamics of continental collision, and in particular what happens to the subducted continental lithosphere after oceanic slab break-off. We find that in some scenarios the subducting continental lithosphere underthrusts the overriding plate not immediately after it enters the trench, but after oceanic slab break-off. In this case, the continental plate first subducts with a steep angle and then, after the slab breaks off at depth, it rises back towards the surface and flattens below the overriding plate, forming a thick horizontal layer of continental crust that extends for about 200 km beyond the suture. This type of behaviour depends on the width of the oceanic plate marginal to the collision zone: wide oceanic margins promote continental underplating and marginal back-arc basins; narrow margins do not show such underplating unless a far field force is applied. Our models show that, as the subducted continental lithosphere rises, the mantle wedge progressively migrates away from the suture and the continental crust heats up, reaching temperatures >900?°C. This heating might lead to crustal melting, and resultant magmatism. We observe a sharp peak in the overriding plate rock uplift right after the occurrence of slab break-off. Afterwards, during underplating, the maximum rock uplift is smaller, but the affected area is much wider (up to 350 km). These results can be used to explain the dynamics that led to the present-day crustal configuration of the India–Eurasia collision zone and its consequences for the regional tectonic and magmatic evolution.
Journal of Geophysical Research: Solid Earth, Vol. 122, 8, pp. 6599-6624.
Mantle
subduction
Abstract: Due to the multidisciplinary nature of combined geodynamics and shear wave splitting studies, there is still much to be understood in terms of isolating the contributions from mantle dynamics to the shear wave splitting signal, even in a two-dimensional (2-D) mantle flow framework. This paper investigates the viscous flow, lattice preferred orientation (LPO) development, and predicted shear wave splitting for a suite of buoyancy-driven subduction models using a non-linear rheology to shed light on the nature of the slab-driven asthenospheric flow and plate-mantle coupling. The slab-driven zone of influence in the mantle, LPO fabric, and resulting synthetic splitting are sensitive to slab strength and slab initial slab dip. The non-linear viscosity formulations leads to dynamic reductions in asthenospheric viscosity extending over 600 km into the mantle wedge and over 300 km behind the trench, with peak flow velocities occurring in models with a weaker slab and moderate slab dip. The olivine LPO fabric in the asthenosphere generally increases in alignment strength with increased proximity to the slab but can be transient and spatially variable on small length scales. The results suggest that LPO formed during initial subduction may persist into the steady state subduction regime. Vertical flow fields in the asthenosphere can produce shear wave splitting variations with back azimuth that deviate from the predictions of uniform trench-normal anisotropy, a result that bears on the interpretation of complexity in shear wave splitting observed in real subduction zones. Furthermore, the models demonstrate the corner flow paradigm should not be equated with a 2-D subduction framework.
Abstract: Impact cratering was a dominant geologic process in the early Solar System that probably played an active role in the crustal evolution of the young terrestrial planets. The Earth’s interior during the Hadean, 4.56 to 4 billion years ago, may have been too hot to sustain plate tectonics. However, whether large impacts could have triggered tectonism on the early Earth remains unclear. Here we conduct global-scale tectonic simulations of the evolution of the Earth through the Hadean eon under variable impact fluxes. Our simulations show that the thermal anomalies produced by large impacts induce mantle upwellings that are capable of driving transient subduction events. Furthermore, we find that moderate-sized impacts can act as subduction triggers by causing localized lithospheric thinning and mantle upwelling, and modulate tectonic activity. In contrast to contemporary subduction, the simulated localized subduction events are relatively short-lived (less than 10?Myr) with relatively thin, weak plates. We suggest that resurgence in subduction activity induced by an increased impact flux between 4.1 and 4.0 billion years ago may explain the coincident increase in palaeointensity of the magnetic field. We further suggest that transient impact-driven subduction reconciles evidence from Hadean zircons for tectonic activity with other lines of evidence consistent with an Earth that was largely tectonically stagnant from the Hadean into the Archaean.
Abstract: Boninites are widely distributed along the western margin of the Pacific Plate extruded during the incipient stage of the subduction zone development in the early Paleogene period. This paper discusses the genetic relationships of boninite and antecedent protoarc basalt magmas and demonstrates their recycled ancient slab origin based on the T-P conditions and Pb-Hf-Nd-Os isotopic modeling. Primitive melt inclusions in chrome spinel from Ogasawara and Guam islands show severely depleted high-SiO2, MgO (high-silica) and less depleted low-SiO2, MgO (low-silica and ultralow-silica) boninitic compositions. The genetic conditions of 1?346?°C at 0.58?GPa and 1?292?°C at 0.69?GPa for the low- and ultralow-silica boninite magmas lie on adiabatic melting paths of depleted mid-ocean ridge basalt mantle with a potential temperature of 1?430?°C in Ogasawara and of 1?370?°C in Guam, respectively. This is consistent with the model that the low- and ultralow-silica boninites were produced by remelting of the residue of the protoarc basalt during the forearc spreading immediately following the subduction initiation. In contrast, the genetic conditions of 1?428?°C and 0.96?GPa for the high-silica boninite magma is reconciled with the ascent of more depleted harzburgitic source which pre-existed below the Izu-Ogasawara-Mariana forearc region before the subduction started. Mixing calculations based on the Pb-Nd-Hf isotopic data for the Mariana protoarc basalt and boninites support the above remelting model for the (ultra)low-silica boninite and the discrete harzburgite source for the high-silica boninite. Yb-Os isotopic modeling of the high-Si boninite source indicates 18-30?wt% melting of the primitive upper mantle at 1.5-1.7?Ga, whereas the source mantle of the protoarc basalt, the residue of which became the source of the (ultra)low-Si boninite, experienced only 3.5-4.0?wt% melt depletion at 3.6-3.1?Ga, much earlier than the average depleted mid-ocean ridge basalt mantle with similar degrees of melt depletion at 2.6-2.2?Ga.
Earth and Planetary Science Letters, Vol. 478, pp. 203-213.
Mantle
subduction
Abstract: After the onset of continental collision, magmatism often persists for tens of millions of years, albeit with a different composition, in reduced volumes, and with a more episodic nature and more widespread spatial distribution, compared to normal arc magmatism. Kinematic modelling studies have suggested that slab breakoff can account for this post-collisional magmatism through the formation of a slab window and subsequent heating of the overriding plate and decompression melting of upwelling asthenosphere, particularly if breakoff occurs at depths shallower than the overriding plate. To constrain the nature of any melting and the geodynamic conditions required, we numerically model the collision of two continental plates following a period of oceanic subduction. A thermodynamic database is used to determine the (de)hydration reactions and occurrence of melt throughout this process. We investigate melting conditions within a parameter space designed to generate a wide range of breakoff depths, timings and collisional styles. Under most circumstances, slab breakoff occurs deeper than the depth extent of the overriding plate; too deep to generate any decompressional melting of dry upwelling asthenosphere or thermal perturbation within the overriding plate. Even if slab breakoff is very shallow, the hot mantle inflow into the slab window is not sustained long enough to sufficiently heat the hydrated overriding plate to cause significant magmatism. Instead, for relatively fast, shallow breakoff we observe melting of asthenosphere above the detached slab through the release of water from the tip of the heating detached slab. Melting of the subducted continental crust during necking and breakoff is a more common feature and may be a more reliable indicator of the occurrence of breakoff. We suggest that magmatism from slab breakoff alone is unable to explain several of the characteristics of post-collisional magmatism, and that additional geodynamical processes need to be considered when interpreting magmatic observations.
Abstract: During continental collision, a part of the lower-plate material can be subducted, emplaced at the base of the upper plate, and eventually incorporated into its crust. This mechanism of continental-crust transformation is called relamination and it has been invoked to explain occurrences of high-pressure felsic rocks in different structural positions of several orogenic systems. In the present study we reproduced relamination during continental collision in a thermo-mechanical numerical model. We performed a parametric study and distinguished three main types of evolution regarding the fate of the subducted continental crust: (i) return along the plate interface in a subduction channel or wedge, (ii) flow at the bottom of the upper-plate lithosphere and subsequent trans-lithospheric exhumation near the arc or in the back-arc region ("sub-lithospheric relamination"), and (iii) nearly horizontal flow directly into the upper-plate crust ("intra-crustal relamination"). Sub-lithospheric relamination is preferred for relatively quick convergence of thin continental plates. An important factor for the development of sub-lithospheric relamination is melting of the subducted material, which weakens the lithosphere and opens a path for the exhumation of the relaminant. In contrast, a thick and strong overriding plate typically leads to exhumation near the plate interface. If the overriding plate is too thin or weak, intra-crustal relamination occurs. We show that each of these evolution types has its counterpart in nature: (i) the Alps and the Caledonides, (ii) the Himalayan-Tibetan system and the European Variscides, and (iii) pre-Cambrian ultra-hot orogens.
Abstract: How new subduction zones form is an emerging field of scientific research with important implications for our understanding of lithospheric strength, the driving force of plate tectonics, and Earth's tectonic history. We are making good progress towards understanding how new subduction zones form by combining field studies to identify candidates and reconstruct their timing and magmatic evolution and undertaking numerical modeling (informed by rheological constraints) to test hypotheses. Here, we review the state of the art by combining and comparing results coming from natural observations and numerical models of SI. Two modes of subduction initiation (SI) can be identified in both nature and models, spontaneous and induced. Induced SI occurs when pre-existing plate convergence causes a new subduction zone to form whereas spontaneous SI occurs without pre-existing plate motion when large lateral density contrasts occur across profound lithospheric weaknesses of various origin. We have good natural examples of 3 modes of subduction initiation, one type by induced nucleation of a subduction zone (polarity reversal) and two types of spontaneous nucleation of a subduction zone (transform collapse and plumehead margin collapse). In contrast, two proposed types of subduction initiation are not well supported by natural observations: (induced) transference and (spontaneous) passive margin collapse. Further work is therefore needed to expand on and understand the implications of these observations. Our future advancing understanding of SI will come from better geologic insights, laboratory experiments, and numerical modeling, and with improving communications between these communities.
Abstract: Water transported by deep subduction to the mantle transition zone (MTZ) that is eventually released and migrates upwards is invoked as a likely cause for hydroweakening and cratonic lithosphere destruction. The destruction of the North China Craton (NCC) during the Mesozoic has been proposed to be related to hydroweakening. However, the source of water related to large-scale craton destruction in the NCC is poorly constrained. Some suggest that the water was mainly released from a flat-lying (or stagnating) slab in the MTZ, whereas others posit that most water was released from a previously existing strongly hydrous MTZ then perturbed by the stagnating subduction in the MTZ layer. In this study, we use numerical modeling to evaluate the water carrying ability of flat-lying slabs in the MTZ with different slab ages and water contents to simulate its maximum value and discuss its potential role on large-scale hydroweakening and craton destruction. Our results reveal that a single flat-lying slab in the MTZ cannot provide enough water for large-scale cratonic lithosphere hydroweakening and thinning. Water estimates invoked for craton destruction as experienced by the NCC can only be the result of long-term piling of multiple slabs in the MTZ or penetrating deeper into the lower mantle.
Geochemistry, Geophysics, Geosystems, Vol. 18, 12, pp. 4469-4486.
Mantle
subduction
Abstract: The classical Wilson Cycle concept, describing repeated opening and closing of ocean basins, hypothesizes spontaneous conversion of passive continental margins into subduction zones. This process, however, is impeded by the high strength of passive margins, and it has never occurred in Cenozoic times. Here using thermomechanical models, we show that additional forcing, provided by mantle flow, which is induced by neighboring subduction zones and midmantle slab remnants, can convert a passive margin into a subduction zone. Models suggest that this is a long-term process, thus explaining the lack of Cenozoic examples. We speculate that new subduction zones may form in the next few tens of millions of years along the Argentine passive margin and the U.S. East Coast. Mantle suction force can similarly trigger subduction initiation along large oceanic fracture zones. We propose that new subduction zones will preferentially originate where subduction zones were active in the past, thus explaining the remarkable colocation of subduction zones during at least the last 400 Myr.
Abstract: Mineral water contents, together with the major and trace element compositions of minerals and whole-rock, were determined for garnet pyroxenites enclosed by ultrahigh-pressure (UHP) metamorphic gneiss at Hujialin in the Sulu orogen. The garnet pyroxenites have low SiO2 contents of 40.25 to 46.68 wt% and MgO contents of 10.99 to 14.79 wt%. They are characterized by enrichment in LREE and LILE (Ba, Sr, Pb) but depletion in HFSE (Nb, Zr) and HREE. They were generated in the Triassic by metasomatic reaction of the mantle wedge peridotite with hydrous felsic melts derived from partial melting of the deeply subducted continental crust. Measured water contents vary from 523 to 1213 ppm for clinopyroxene, and 55 to 1476 ppm for garnet. These mineral water contents are not only correlated with mineral major and trace element abundances but also relatively homogenous within single mineral grains. Such features preclude significant disturbance of the mineral water contents during pyroxenite exhumation from the mantle depth to the surface and thus indicate preservation of the primary water contents for the UHP metasomatites. The garnet pyroxenites are estimated to have bulk water contents of 424-660 ppm, which are higher than those for the MORB source, similar to or higher than those for the OIB sources and close to the lower limit for the arc magma source. The relationships between contents of mineral water and some elements suggest that the high water contents of garnet pyroxenites are primarily determined by the abundance of water-rich clinopyroxene. Garnet also has the high water contents, suggesting its importance in hosting water at mantle depths. Calculated whole-rock H2O/Ce ratios are 63-145, higher than those for Hawaiian garnet pyroxenites and SWIR abyssal pyroxenites. These observations suggest that metasomatic pyroxene-rich lithologies have the capacity to contribute high H2O concentrations and variable H2O/Ce ratios to the mantle. This lends support to the interpretation that the source of some intraplate basalts may be a heterogeneous mixture of peridotite and pyroxenite. On the other hand, the high water contents of garnet pyroxenites suggest that the presence of ultramafic metasomatites in the mantle wedge would enhance its water storage and thus reduce the water transport into deeper mantle by subduction.
Abstract: How new subduction zones form is an emerging field of scientific research with important implications for our understanding of lithospheric strength, the driving force of plate tectonics, and Earth's tectonic history. We are making good progress towards understanding how new subduction zones form by combining field studies to identify candidates and reconstruct their timing and magmatic evolution and undertaking numerical modeling (informed by rheological constraints) to test hypotheses. Here, we review the state of the art by combining and comparing results coming from natural observations and numerical models of SI. Two modes of subduction initiation (SI) can be identified in both nature and models, spontaneous and induced. Induced SI occurs when pre-existing plate convergence causes a new subduction zone to form whereas spontaneous SI occurs without pre-existing plate motion when large lateral density contrasts occur across profound lithospheric weaknesses of various origin. We have good natural examples of 3 modes of subduction initiation, one type by induced nucleation of a subduction zone (polarity reversal) and two types of spontaneous nucleation of a subduction zone (transform collapse and plumehead margin collapse). In contrast, two proposed types of subduction initiation are not well supported by natural observations: (induced) transference and (spontaneous) passive margin collapse. Further work is therefore needed to expand on and understand the implications of these observations. Our future advancing understanding of SI will come from better geologic insights, laboratory experiments, and numerical modeling, and with improving communications between these communities.
Abstract: Water transported by deep subduction to the mantle transition zone (MTZ) that is eventually released and migrates upwards is invoked as a likely cause for hydroweakening and cratonic lithosphere destruction. The destruction of the North China Craton (NCC) during the Mesozoic has been proposed to be related to hydroweakening. However, the source of water related to large-scale craton destruction in the NCC is poorly constrained. Some suggest that the water was mainly released from a flat-lying (or stagnating) slab in the MTZ, whereas others posit that most water was released from a previously existing strongly hydrous MTZ then perturbed by the stagnating subduction in the MTZ layer. In this study, we use numerical modeling to evaluate the water carrying ability of flat-lying slabs in the MTZ with different slab ages and water contents to simulate its maximum value and discuss its potential role on large-scale hydroweakening and craton destruction. Our results reveal that a single flat-lying slab in the MTZ cannot provide enough water for large-scale cratonic lithosphere hydroweakening and thinning. Water estimates invoked for craton destruction as experienced by the NCC can only be the result of long-term piling of multiple slabs in the MTZ or penetrating deeper into the lower mantle.
Abstract: Carbonatitic magmatism plays a significant role in Earth's carbon cycle, which is also a lithoprobe of crust-mantle interaction, mantle metasomatism and partial melting. Due to different mineral assemblages and geochemical compositions, and diverse tectonic settings, the origin of carbonatite has long been debated. At subduction zones, sediments (including carbonates) are subducted into the mantle with the downgoing oceanic slab. However, the detailed mechanism of how subducted carbonates contribute to carbonatitic magmatism remains unclear. Here we present geochronological, geochemical and isotopic study on the Taohuala Mountain carbonatite at the southern margin of the Alxa Block, North China Craton. The classification of carbonatite from the Taohuala Mountain relies strongly on the observations of obvious intrusion contact relationships and flow structures in field outcrop. The Taohuala Mountain carbonatite has SiO2 ranging from 2.37 wt.% to 11.45 wt%, high CaO (45.93-53.86 wt%) and low MgO (0.51-4.39 wt%), and is characterized by enrichment of LILE (Ba, Sr), depletion of HFSE (Nb, Ta, Zr, Hf), and slightly negative Ce and Eu anomalies. Carbonates in the samples have high 87Sr/86Sr (0.70686-0.70694) and low 143Nd/144Nd (0.511635-0.511924). Remarkably, the highly fractionated ?18OVSMOW (11.83-25.92‰) indicates components of both sedimentary and mantle origin. Detailed zircon in situ U-Pb dating and oxygen isotope analysis exhibit contrast ages and ?18OVSMOW from core to rim, i.e., old ages (mainly > 800 Ma), high Th/U (mainly > 0.5) and low ?18OVSMOW (6.37-11.44‰) in cores (inherited), whereas young ages (~ 400 Ma), low Th/U (mainly < 0.01) and high ?18OVSMOW (20.04-24.54‰) in rims, suggesting that the Taohuala Mountain carbonatite may have been generated from melting of subducted sedimentary carbonates. Considering all these evidences, and that the collision along Qilian Mountains was older than the carbonatite, we propose that a large volume of sedimentary carbonates subducted and remained in the lithospheric mantle under the Alxa block during the closure of the Paleo-Qilian Ocean. Subsequently, the carbonatite was formed by melting of carbonates with minor contributions from the mantle during the breakoff or rollback of the Paleo-Asian oceanic slab.
Abstract: Subducting serpentinized lithosphere has distinct ?D and ?18O values compared to normal mantle. Slab-derived fluids that infiltrate the mantle wedge can alter its oxygen and hydrogen isotope composition, raising or lowering the ?18O and ?D values depending on the nature of the subducted components. Hydrous minerals in peridotite xenoliths from the Colorado Plateau (southwestern USA) have ?D values (up to ?33‰) much higher than average mantle (?80‰), but similar to ?D values of olivine-hosted melt inclusions within arc basalts, suggesting a slab-derived fluid source. Oxygen isotope ratios of olivine from these xenoliths are similar to average mantle, yet display a strong negative correlation with clinopyroxene Ce/Sm, a proxy of metasomatism. This correlation is most simply explained by metasomatism from fluids derived from the serpentinized portion of the Farallon slab. Although ?18O values of mantle minerals span a narrow range, integration of stable isotope data with other geochemical tracers can provide new constraints on modern and ancient subduction-related processes, potentially providing a method for probing Archean lithospheric mantle for evidence of early subduction.
Earth and Planteray Science Letters, Vol. 489, pp. 258-266.
Mantle
subduction
Abstract: The impact of remotely forced mantle flow on regional subduction evolution is largely unexplored. Here we investigate this by means of 3D thermo-mechanical numerical modeling using a regional modeling domain. We start with simplified models consisting of a 600 km (or 1400 km) wide subducting plate surrounded by other plates. Mantle inflow of ?3 cm/yr is prescribed during 25 Myr of slab evolution on a subset of the domain boundaries while the other side boundaries are open. Our experiments show that the influence of imposed mantle flow on subduction evolution is the least for trench-perpendicular mantle inflow from either the back or front of the slab leading to 10-50 km changes in slab morphology and trench position while no strong slab dip changes were observed, as compared to a reference model with no imposed mantle inflow. In experiments with trench-oblique mantle inflow we notice larger effects of slab bending and slab translation of the order of 100-200 km. Lastly, we investigate how subduction in the western Mediterranean region is influenced by remotely excited mantle flow that is computed by back-advection of a temperature and density model scaled from a global seismic tomography model. After 35 Myr of subduction evolution we find 10-50 km changes in slab position and slab morphology and a slight change in overall slab tilt. Our study shows that remotely forced mantle flow leads to secondary effects on slab evolution as compared to slab buoyancy and plate motion. Still these secondary effects occur on scales, 10-50 km, typical for the large-scale deformation of the overlying crust and thus may still be of large importance for understanding geological evolution.
Abstract: Alkaline lavas occur globally in subduction-related volcanic arcs. Conventional models for the origin of these lavas typically invoke a multi-stage process in which mantle wedge peridotite, enriched in phlogopite and/or amphibole due to prior metasomatism, partially melts during infiltration by fluids and melts derived from subducted oceanic lithosphere. However, geochemical systematics in the majority of subduction-related alkaline lavas require physical mixing of subducted components and peridotite prior to partial melting. This can be explained by the mélange diapir model, which predicts the generation of arc magmas during advection of buoyant material from the slab-wedge interface into the mantle wedge below arcs. Here we report results from experiments in which natural mélange materials were partially melted at upper mantle conditions to produce alkaline magmas. Partial melts produced in our experiments have trace-element abundance patterns that are typical of alkaline arc lavas, such as enrichment in large ion lithophile elements (LILEs) and depletion in Nb and Ta. These results favor generation of alkaline magmas in the arc and backarc regions of subduction zones by partial melting of mélange materials rather than previously metasomatized peridotite.
Geochemistry, Geophysics, Geosystems, Vol. 19, 3, pp. 865-881.
Mantle
subduction
Abstract: The buoyancy of continental crust opposes its subduction to mantle depths, except where mineral reactions substantially increase rock density. Sluggish kinetics limit such densification, especially in dry rocks, unless deformation and hydrous fluids intervene. Here we document how hydrous fluids in the subduction channel invaded lower crustal granulites at 50-60 km depth through a dense network of probably seismically induced fractures. We combine analyses of textures and mineral composition with thermodynamic modeling to reconstruct repeated stages of interaction, with pulses of high-pressure (HP) fluid at 650-6708C, rehydrating the initially dry rocks to micaschists. SIMS oxygen isotopic data of quartz indicate fluids of crustal composition. HP growth rims in allanite and zircon show uniform U-Th-Pb ages of 65 Ma and indicate that hydration occurred during subduction, at eclogite facies conditions. Based on this case study in the Sesia Zone (Western Italian Alps), we conclude that continental crust, and in particular deep basement fragments, during subduction can behave as substantial fluid sinks, not sources. Density modeling indicates a bifurcation in continental recycling: Chiefly mafic crust, once it is eclogitized to >60%, are prone to end up in a subduction graveyard, such as is tomographically evident beneath the Alps at 550 km depth. By contrast, dominantly felsic HP fragments and mafic granulites remain positively buoyant and tend be incorporated into an orogen and be exhumed with it. Felsic and intermediate lithotypes remain positively buoyant even where deformation and fluid percolation allowed them to equilibrate at HP.
Mineralogy and Petrology, 10.1007/ s00710-018- 0577-8, 19p.
Africa
metasomatism, subduction, geobarometry
Abstract: The Kaapvaal craton and its underlying mantle is probably one of the best studied Archean entity in the world. Despite that, discussion is still vivid on important aspects. A major debate over the last few decades is the depth of melting that generated the mantle nuclei of cratons. Our new evaluation of melting parameters in peridotite residues shows that the Cr2O3/Al2O3 ratio is the most useful pressure sensitive melting barometer. It irrevocably constrains the pressure of melting (melt separation) to less than 2 GPa with olivine (ol), orthopyroxene (opx) and spinel (sp) as residual phases. Garnet (grt) grows at increasing pressure during lithosphere thickening and subduction via the reaction opx?+?sp ? grt?+?ol. The time of partial melting is constrained by Re-depletion model ages (TRD) mainly to the Archean (Pearson and Wittig 2008). However, only 3% of the ages are older than 3.1 Ga while crustal ages lie mainly between 3.1 to 2.8 Ga for the W- and 3.7 to 2.8 Ga for the E-block. Many TRD-ages are probably falsified by metasomatism and the main partial melting period was older than 3.1 Ga. Also, Nd- and Hf- model ages of peridotitic lithologies from the W-block are 3.2 to 3.6 Ga old. The corresponding very negative ?Nd (?40) and ?Hf values (?65) signal the presence of subducted crustal components in these old mantle portions. Subducted components diversify the mantle in its chemistry and thermal structure. Adjustment towards a stable configuration occurs by fluid transfer, metasomatism, partial melting and heat transfer. Ages of metasomatism from the Lu-Hf isotope system are 3.2 Ga (Lace), 2.9 Ga (Roberts Victor) and 2.62 Ga (Finsch) coinciding with the collision of cratonic blocks, the growth of diamonds, metamorphism of eclogites and of Ventersdoorp magmatism. The cratonic lithosphere was stabilized thermally by the end of the Archean and cooled since then with a rate of 0.07 °C/Ma.
Geophysical Research Letters, Vol. 45, 11, pp. 5336-5343.
Mantle
subduction
Abstract: Rocks on the Earth's surface are cooled, hardened, eventually forming rigid plates that move around relative to one another. When two plates converge, one plate overrides the other, which sinks into the Earth's deep mantle. The sinking plate carries water, which softens rocks and also affects the behavior of the sinking/overriding plates and surrounding mantle flows (“subduction dynamics”). To investigate the role of water in subduction dynamics, 2?D fluid dynamical simulations were performed. The simulations suggest that subduction dynamics change significantly with the level of hydration of the sinking plate, which is represented by the thickness of a hydrous layer. When the hydrous layer is thin, the plate sinks rapidly with a shifting boundary and stagnates above the lower mantle. In contrast, when the hydrous layer is thick, plate convergence is sluggish, the plate boundary remains stationary, and the sinking plate penetrates into the lower mantle. These results indicate that a small amount of water is expected for the northwest part of the Pacific Plate, characterized by the rapid convergence, plate boundary shifting, and stagnation of the sinking plate.
Abstract: Water is carried down into Earth at subduction zones by the process of plate tectonics. Much of the water escapes close behind the subduction zone, promoting melting of the mantle and giving rise to the volcanic chains in the Pacific Ocean basin known as the Ring of Fire, and many other volcanoes elsewhere. But can water be carried even further into the mantle? How would we even know? Why is it important, and what are the effects of such deep water storage? Diamonds can give us the answers to these questions. Recent discoveries of water-containing mineral inclusions and even free water held at high pressures in diamonds tell us that water is carried into Earth’s deep interior—perhaps as deep as 700 km.
Abstract: Type IIb diamonds, such as the Hope diamond, contain trace amounts of boron and are prized for their blue colors. Since boron is a quintessential crustal element, it is completely unexpected in diamond-forming fluids at mantle depths. Despite the mineralogical/geochemical interest in type IIb diamonds, almost nothing is known about how they form chiefly because of their rarity (?0.02% of all diamonds) and high gem value. To investigate the type of mantle host rock, the depth of origin (lithospheric vs. convecting mantle), and the source of boron, the high-volume diamond grading stream of the Gemological Institute of America was systematically screened to find type IIb diamonds with inclusions. Over a period of about two years, 46 prospective samples were identified and examined optical microscopy, X-ray diffraction, and infrared/Raman spectroscopy; a few diamonds were also analyzed for carbon isotopic composition and polished for electron probe microanalysis of inclusions. The examined inclusions represent retrogressed highpressure minerals, from metabasic to metaperidotitic hosts in the lowermost mantle transition zone (MTZ) to lower mantle (LM). These include former CaSiO3-perovskite, majorite, bridgmanite, stishovite, calcium-ferrite-type phase, and ferropericlase. The variably light carbon isotope compositions and inclusion mineralogy indicate diamond growth in deeply subducted oceanic lithosphere (crust and mantle). Some inclusions are found to have coexisting fluid (CH4 ± H2) that suggests the original high-pressure minerals interacted with hydrous media. We propose that the boron resided in serpentinized oceanic lithosphere. During subduction, the serpentine was metamorphosed to dense hydrous magnesium silicates (DHMS) that retained some boron. Upon breakdown in the MTZ/LM, these DHMS yielded boron-bearing hydrous fluids conducive to diamond growth.
Abstract: The H/C ratio in earth’s exosphere is higher than it is in the source region of primitive basalts, suggesting an enriched carbon reservoir in the mantle[1]. A plausible explanation is that subduction of carbon may have enriched the mantle in recycled carbon over time. Average basaltic crust contains ~ 2 wt.% CO2 [2], and modeling of slab devolatilisation suggests that subducted carbonate may survive to be transported deeper into the mantle [3]. Carbonated oceanic crust should melt in the transition zone along most subduction geotherms due to a deep trough in the carbonated basalt solidus, and mineral inclusions in superdeep diamonds testify to carbonate melt in their formation [4]. Along cool subduction geotherms carbonate may subduct into the lower mantle, potentially enriching the deep mantle in carbon. Here we report on laser-heated diamond anvil cell experiments in the CaO-MgO-SiO2-CO2 and FeO-MgO-SiO2-CO2 systems at lower mantle pressures where we investigate the stability of carbonate in oceanic crust, and test for decarbonation and diamond forming reactions involving carbonate and coexisiting free silica. We find that carbonate reacts with silica to form bridgmanite ± Ca-perovskite + CO2 at pressures in the range of ~50 to 70 GPa. These decarbonation reactions form an impenetrable barrier to subduction of carbonate into the deeper lower mantle, however, slabs may carry solid CO2 (Phase V) into the deeper lower mantle. We also identify reactions where carbonate or CO2 dissociate to form diamond plus oxygen. We suggest that the deep lower mantle may become enriched in carbon in the form of diamond over time due to subduction of carbonate and solid CO2 and its eventual dissociation to form diamond plus oxygen. Release of oxygen during diamond formation may also provide a mechanism for locally oxidizing the deep mantle.
Abstract: The continental basement in the Eastern Mediterranean represents the northern edge of Gondwana, which has been the site of repeated crustal accretion and has subsequently been modified by consecutive rifting events. We investigated the geologic and thermal history of the North Gondwana lower crust by examining the U-Pb-Hf isotope systematics in zircons within 6 mafic granulite xenoliths from Pliocene lava cone in North Israel. The lava cone protrudes through the platform cover that seals the late Neoproterozoic junction between the Arabian-Nubian basement to the South and the Cadomian basement exposed in the Taurides to the North. The mafic granulite xenoliths are composed of plagioclase + orthopyroxene + clinopyroxene ± garnet ± spinel ± secondary amphibole. U-Pb zircon ages from the granulites vary among the different samples with distinct zircon age populations at 400-1200 Ma, 170-350 Ma, and 3.6-4.2 Ma, attesting the lower crust preserves a prolonged thermal and igneous history. While 400-550 Ma U-Pb ages are interpreted to be the result of Pb loss, the wide scatter of zircon grains aged between 550 and 1200 Ma, alongside their diverse ?Hf(t) values (?25-+10), is an extraordinary evidence for the accretion of Neoproterozoic sediments into the North Gondwana arc root lower crust. The U-Pb-Hf signature of these zircons resembles Cadomian sediments of the Tauride block to the north, indicating southward (present coordinates) subduction under North Gondwana and possible accretion of fore-arc sediments to the lower crust through relamination in the latest Neoproterozoic. One xenolith contained metamorphic-shaped zircons aged 170-350 Ma with positive ?Hf values and Hf-TDM of 0.85 Ga interpreted to reflect Paleozoic recycling of the Neoproterozoic juvenile Arabian basement, which we consider to form a major component of the lower crust in the region. An overwhelming cluster of Carboniferous zircons concentrating at 305 Ma with exclusively negative ?Hf values around ?6, was retrieved from three xenoliths. Some of these zircons portrayed igneous textures and shape. While Carboniferous igneous activity is the hallmark of Western Europe's Variscan orogeny, the latter did not affect the southern rifted edge of Neo-Tethys where our xenoliths were retrieved. The Paleozoic age-Hf composition in our xenoliths is therefore interpreted to result from syn-Variscan recycling of Neoproterozoic sedimentary remains in the lower crust, and some degree of melting in a non-orogenic environment. Rather than with horizontal plate motions and orogeny, the Carboniferous zircon ages in the xenoliths appear to coalesce with significant vertical movements that created continental scale unconformities and a broad basin and swell architecture known to develop over the entire North Gondwana margin at that time. The Carboniferous aged zircons in northern Israel lower crustal xenoliths are therefore a unique gauge of the thermal perturbation that accompanied the large-scale mantle dynamics below the then passive North African margin of Gondwana, while Variscan orogenic accretion occurred on the Eurasian margin. These lower crustal granulites xenoliths therefore contain important information with respect to the nature of the lower crust under Israel, with implications on the geodynamic setting during the Cadomian and Variscan cycles.
Physics of the Earth and Planetary Interiors, Vol. 269, 1, pp. 112-132.
Mantle
subduction
Abstract: The present study aims at better deciphering the different mechanisms involved in the functioning of the subduction interplate. A 2D thermo-mechanical model is used to simulate a subduction channel, made of oceanic crust, free to evolve. Convergence at constant rate is imposed under a 100 km thick upper plate. Pseudo-brittle and non-Newtonian behaviours are modelled. The influence of the subduction channel strength, parameterized by the difference in activation energy between crust and mantle () is investigated to examine in detail the variations in depth of the subduction plane down-dip extent, . First, simulations show that numerical resolution may be responsible for an artificial and significant shallowing of if the weak crustal layer is not correctly resolved. Second, if the age of the subducting plate is 100 Myr, subduction occurs for any . The stiffer the crust is, that is, the lower is, the shallower is (60 km depth if kJ/mol) and the hotter the fore-arc base is. Conversely, imposing a very weak subduction channel ( J/mol) leads there to an extreme mantle wedge cooling and inhibits mantle melting in wet conditions. Partial kinematic coupling at the fore-arc base occurs if kJ/mol. If the incoming plate is 20 Myr old, subduction can occur under the conditions that the crust is either stiff and denser than the mantle, or weak and buoyant. In the latter condition, cold crust plumes rise from the subduction channel and ascend through the upper lithosphere, triggering (1) partial kinematic coupling under the fore-arc, (2) fore-arc lithosphere cooling, and (3) partial or complete hindrance of wet mantle melting. then ranges from 50 to more than 250 km depth and is time-dependent if crust plumes form. Finally, subduction plane dynamics is intimately linked to the regime of subduction-induced corner flow. Two different intervals of are underlined: 80-120 kJ/mol to reproduce the range of slab surface temperature inferred from geothermometry, and 10-40 kJ/mol to reproduce the shallow hot mantle wedge core inferred from conditions of last equilibration of near-primary arc magmas and seismic tomographies. Therefore, an extra process controlling mantle wedge dynamics is needed to satisfy simultaneously the aforementioned observations. A mantle viscosity reduction, by a factor 4-20, caused by metasomatism in the mantle wedge is proposed. From these results, I conclude that the subduction channel down-dip extent, , should depend on the subduction setting, to be consistent with the observed variability of sub-arc depths of the subducting plate surface.
Earth and Planetary Science Letters, Vol. 502, pp. 166-173.
India
subduction
Abstract: Tectonic plate velocities predominantly result from a balance between the potential energy change of the subducting slab and viscous dissipation in the mantle, bending lithosphere, and slab-upper plate interface. A range of observations suggest that slabs may be weak, implying a more prominent role for plate interface dissipation than previously thought. The shallow thrust interface is commonly assumed to be weak due to an abundance of fluids and near-lithostatic pore fluid pressures, but little attention has been paid to the influence of the deeper, viscous interface. Here we show that the deep interface viscosity in subduction zones is strongly affected by the relative proportions of sedimentary to mafic rocks that are subducted to depth. Where sediments on the down-going plate are sparse, the deep interface is dominated by mafic lithologies that metamorphose to eclogites, which exhibit viscosities 1-2 orders of magnitude higher than the asthenospheric mantle, and reduce subduction plate speeds. In contrast, where sediments are abundant and subducted to depth, the deep interface viscosity is 1-2 orders of magnitude lower than the asthenospheric mantle, thus allowing significantly faster plate velocities. This correlation between subduction plate speed and deep sediment subduction may help explain dramatic accelerations (or decelerations) in convergence rates, such as the acceleration documented for India-Asia convergence during the mid-Cenozoic.
Nature Geoscience, 10.1038/s41561-018-0215-4 pp. 682-687.
Russia, Siberia
subduction
Abstract: Magmatic volatile release to the atmosphere can lead to climatic changes and substantial environmental degradation including the production of acid rain, ocean acidification and ozone depletion, potentially resulting in the collapse of the biosphere. The largest recorded mass extinction in Earth’s history occurred at the end of the Permian, coinciding with the emplacement of the Siberian large igneous province, suggesting that large-scale magmatism is a key driver of global environmental change. However, the source and nature of volatiles in the Siberian large igneous province remain contentious. Here we present halogen compositions of sub-continental lithospheric mantle xenoliths emplaced before and after the eruption of the Siberian flood basalts. We show that the Siberian lithosphere is massively enriched in halogens from the infiltration of subducted seawater-derived volatiles and that a considerable amount (up to 70%) of lithospheric halogens are assimilated into the plume and released to the atmosphere during emplacement. Plume-lithosphere interaction is therefore a key process controlling the volatile content of large igneous provinces and thus the extent of environmental crises, leading to mass extinctions during their emplacement.
Mineralogy and Petrology, Vol. 112, 5, pp. 603-616.
Europe, Austria
subduction
Abstract: Mafic layers displaying transition between clinopyroxenite and eclogite within peridotite from felsic granulite in the Bohemian Massif (Lower Austria) have been investigated. The mafic-ultramafic bodies shared a common granulite facies metamorphism with its hosting felsic rocks, but they still preserve evidence of eclogite facies metamorphism. The selected mafic layer for this study is represented by garnet with omphacite in the core of coarse-grained clinopyroxene, while fine-grained clinopyroxene in the matrix is diopside. In addition, garnet contains inclusions of omphacite, alkali feldspars, hydrous and other phases with halogens and/or CO2. Textural relations along with compositional zoning in garnet from the clinopyroxenite-eclogite layers favour solid-state recrystallization of the precursor minerals in the inclusions and formation of garnet and omphacite during subduction. Textures and major and trace element distribution in garnet indicate two stages of garnet growth that record eclogite facies and subsequent granulite facies overprint. The possible model explaining the textural and compositional changes of minerals is that the granulite facies overprint occurred after formation and exhumation of the eclogite facies rocks.
Journal of Structual Geology, Vol. 115, pp. 284-296.
Mantle
subduction
Abstract: During the last decades, processes taking place in subduction zones have been a topic of extensive debate. There are many models on a tectonic scale, but knowledge of deformation at microscopic levels is essential for an understanding of the processes involved. Crystallographic preferred orientation (CPO) investigations of minerals in eclogites yield the possibility to study deformation in the crustal part of subducted oceanic lithosphere, as well as subducted continental units containing basic intrusions, which are frequently exhumed during continental collision. From CPO data of omphacite and garnet, as well as other constituent minerals in eclogites, conclusions on stress and strain at depth in currently active subduction zones can be drawn, as the pressure-temperature and strain path of the now exhumed rocks can be unraveled. This contribution provides an overview of CPO studies of eclogites concerning slip systems, deformation, strain and other possible CPO forming mechanisms in all major constituent mineral phases in eclogites and the implications of these microstructural data on subduction zone processes are discussed. Nevertheless, there are still many open questions and future research is essential.
Earth and Planetary Science Letters, Vol. 501, pp. 67-77.
Mantle
subduction
Abstract: The mechanisms underpinning the formation of a focused volcanic arc above subduction zones are debated. Suggestions include controls by: (i) where the subducting plate releases water, lowering the solidus in the overlying mantle wedge; (ii) the location where the mantle wedge melts to the highest degree; and (iii) a limit on melt formation and migration imposed by the cool shallow corner of the wedge. Here, we evaluate these three proposed mechanisms using a set of kinematically-driven 2D thermo-mechanical mantle-wedge models in which subduction velocity, slab dip and age, overriding-plate thickness and the depth of decoupling between the two plates are systematically varied. All mechanisms predict, on the basis of model geometry, that the arc-trench distance, D, decreases strongly with increasing dip, consistent with the negative D-dip correlations found in global subduction data. Model trends of sub-arc slab depth, H, with dip are positive if H is wedge-temperature controlled and overriding-plate thickness does not exceed the decoupling depth by more than 50 km, and negative if H is slab-temperature controlled. Observed global H-dip trends are overall positive. With increasing overriding plate thickness, the position of maximum melting shifts to smaller H and D, while the position of the trenchward limit of the melt zone, controlled by the wedge's cold corner, shifts to larger H and D, similar to the trend in the data for oceanic subduction zones. Thus, the limit imposed by the wedge corner on melting and melt migration seems to exert the first-order control on arc position.
Earth and Planetary Science Letters, Vol. 502, pp. 96-103.
Mantle
subduction
Abstract: Since several decades, the processes allowing for the subduction of the continental lithosphere less dense than the mantle in a collision context have been widely explored, but models that are based upon the premise that slab pull is the prominent driver of plate tectonics fail. The India-Asia collision, where several episodes of continental subduction have been documented, constitute a case study for alternative views. One of these episodes occurred in the early collision time within the Asian plate where continental lithosphere not attached to any oceanic lithosphere subducted southward in front of the Indian lithosphere during its northward subduction that followed the oceanic subduction of the Tethys ocean. This process, known as collisional subduction, has a counter-intuitive behavior since the subduction is not driven by slab pull. It has been speculated that the mantle circulation can play an important role in triggering collisional subduction but a detailed, qualitative analysis of it is not available, yet. In this work we explore the southward subduction dynamics of the Asian lithosphere below Tibet by means of analogue experiments with the aim to highlight how the mantle circulation induces or responds to collisional subduction. We found that during the northward oceanic subduction (analogue of Tethys subduction) attached to the indenter (Indian analogue), the main component of slab motion is driven vertically by its negative buoyancy, while the trench rolls back. In the mantle the convective pattern consists in a pair of wide convective cells on both sides of the slab. But when the indenter starts to bend and plunge in the mantle, trench motion reverses. Its advance transmits the far field forces to two upper plates (Asian analogues). The more viscous frontal plate thickens, and the less viscous hinterland plate, which is attached to the back wall of the box, subducts. During this transition, a pair of sub-lithospheric convective cells is observed on both sides of the Asian analogue slab, driven by the shortening of the frontal plate. It favors the initiation of the backwall plate subduction. Such subduction is maintained during the entire collision by a wide cell with a mostly horizontal mantle flow below Tibet, passively advecting the Asian analogue slab. Experimental results suggest that once the tectonic far-field force related to the forward horizontal motion becomes dominant upon the buoyancy forces, trench advancing and the transmission of the tectonic force to the upper and backwall plates are promoted. This peculiar condition triggers the subduction of the backwall plate, despite it is light and buoyant.
Journal of African Earth Sciences, Vol. 149, pp. 215-234.
Africa, Cameroon
subduction
Abstract: High- to ultrahigh-pressure metamorphic assemblages consisting of garnet-omphacitic clinopyroxene bearing mafic rocks have been identified within the Paleoproterozoic Nyong Group in SW Cameroon, at the northwestern margin of the Archean Congo craton. These rocks were investigated in detail and for the first time evidence for eclogite facies metamorphism at ca 25?kbar and 850?°C is provided. A clockwise P-T path with nearly isothermal decompression (ITD) is deduced from mineral zoning and textural relationships characterized by mineral recrystallization and multi-layered coronitic overgrowths of plagioclase and clinopyroxene surrounding garnet porphyroblasts. These P-T conditions imply a burial depth greater than 90?km, at lower geothermal gradient of ca 10?°C/km. The geochemical signature of ten representative rock samples show that two groups of eclogite facies rocks genetically originate from mostly basaltic and basaltic andesite compositions, with a characteristic upper mantle-derived tholeiitic trend. Moreover, their chondrite and MORB normalized REE and trace element concentrations are characterized by nearly flat REE patterns with very little to no Eu anomaly, (La/Sm)N???1 and Zr/Nb???10, as well as a gradual depletion from LREE to HREE with also very little to no Eu anomaly, but (La/Sm)N < 1, Zr/Nb > 10 and negative anomalies in Th, K, Nb, Ta, Sr, Zr and Ti consistent with mid-ocean ridge basalt (MORB) contaminated by a subduction component or by a crustal component. Previous available geochronological data coupled with our new petrological, mineralogical and geochemical findings clearly indicate that the eclogite facies metabasites from the Eburnean Nyong Group between 2100 and 2000 Ma represent one of the oldest subducted oceanic slab or trace of a suture zone so far recorded within the West Central African Fold Belt (WCAFB). The geodynamic implications of these eclogites suggest a subduction-related process followed by a rapid exhumation of their protoliths, therefore, providing critical information corroborating that plate tectonic processes operated during the Paleoproterozoic.
Contributions to Mineralogy and Petrology, Vol. 173, 21p. Doi.org/10.1007/s00410-018-1515-9
Mantle
subduction
Abstract: The role of aqueous fluid in fracturing in subducting slabs was investigated through a series of deformation experiments on dunite that was undersaturated (i.e., fluid-free) or saturated with water (i.e., aqueous-fluid bearing) at pressures of 1.0-1.8 GPa and temperatures of 670-1250 K, corresponding to the conditions of the shallower regions of the double seismic zone in slabs. In situ X-ray diffraction, radiography, and acoustic emissions (AEs) monitoring demonstrated that semi-brittle flow associated with AEs was dominant and the creep/failure strength of dunite was insensitive to the dissolved water content in olivine. In contrast, aqueous fluid drastically decreased the creep/failure strength of dunite (up to ~ 1 GPa of weakening) over a wide range of temperatures in the semi-brittle regime. Weakening of the dunite by the aqueous fluid resulted in the reduction of the number of AE events (i.e., suppression of microcracking) and shortening of time to failure. The AE hypocenters were located at the margin of the deforming sample while the interior of the faulted sample was aseismic (i.e., aseismic semi-brittle flow) under water-saturated conditions. A faulting (slip rate of ~ 10?³ to 10?? s?¹) associated with a large drop of stress (?? ~ 0.5 to 1 GPa) and/or pressure (?P ~ 0.5 GPa) was dominant in fluid-free dunite, while a slow faulting (slip rate < 8 × 10?? s?¹) without any stress/pressure drop was common in water-saturated dunite. Aseismic semi-brittle flow may mimic silent ductile flow under water-saturated conditions in subducting slabs.
Abstract: Short- and long-term processes at or close to the subduction plate interface (e.g.,mineral transformations, fluid release, seismicity and more generally deformation) might be more closely related than previously thought. Increasing evidence from the fossil rock record suggests that some episodes of their long geological evolution match or are close to timescales of the seismic cycle. This contribution uses rocks recovered (episodically) from subduction zones, together with insights from thermomechanical modelling, to provide a new dynamic vision of the nature, structure and properties of the plate interface and to bridge the gap between the mechanical behavior of active subduction zones (e.g.,coupling inferred from geophysical monitoring) and fossil ones (e.g.,coupling required to detach and recover subducted slab fragments). Based on critical observations and an exhaustive compilation of worldwide subducted oceanic units (for which the presence near the plate interface, rock types, pressure, temperature, T/P gradients, thickness and timing of detachment can be assessed), the present study demonstrates how long-term mechanical coupling exerts a key control on detachment from the slab and potential rock recovery. Critical assessment of rock T/P characteristics indicates that these fragments can indeed be used as natural probes and provide reliable information on subduction interface dynamics down to ~2.8?GPa. Rock clusters are identified at depths of 30, 5560 and 80?km, with some differences between rock types. Data also reveal a first-order evolution with subduction cooling (in the first ~5?Myr), which is interpreted as reflecting a systematic trend from strong to weak mechanical coupling, after which subduction is lubricated and mostly inhibits rock recovery. This contribution places bounds on the plate interface constitution, regular thickness (<300?m; i.e. where/when there is no detachment), changing geometry and effective viscosity. The concept of ‘coupled thickness' is used here to capture subduction interface dynamics, notably during episodes of strong mechanical coupling, and to link long- and short-term deformation. Mechanical coupling depends on mantle wedge rheology, viscosity contrasts and initial structures (e.g.,heterogeneous lithosphere, existence of décollement horizons, extent of hydration, asperities) but also on boundary conditions (convergence rates, kinematics), and therefore differs for warm and cold subduction settings. Although most present-day subduction zone segments (both along strike and downdip) are likely below the detachment threshold, we propose that the most favorable location for detachment corresponds to the spatial transition between coupled and decoupled areas. Effective strain localization involves dissolution-precipitation and dislocation creep but also possibly brittle fractures and earthquakes, even at intermediate depths.
Geochemistry, Geophysics, Geosystems, Vol. 19, 9, pp. 3140-3163.
United States, California
subduction
Abstract: Mantle convection shapes Earth's surface by generating dynamic topography. Observational constraints and regional convection models suggest that surface topography could be sensitive to mantle flow for wavelengths as short as 1,000 and 250 km, respectively. At these spatial scales, surface processes including sedimentation and relative sea?level change occur on million?year timescales. However, time?dependent global mantle flow models do not predict small?scale dynamic topography yet. Here we present 2?D spherical annulus numerical models of mantle convection with large radial and lateral viscosity contrasts. We first identify the range of Rayleigh number, internal heat production rate and yield stress for which models generate plate?like behavior, surface heat flow, surface velocities, and topography distribution comparable to Earth's. These models produce both whole?mantle convection and small?scale convection in the upper mantle, which results in small?scale (<500 km) to large?scale (>104 km) dynamic topography, with a spectral power for intermediate scales (500 to 104 km) comparable to estimates of present?day residual topography. Timescales of convection and the associated dynamic topography vary from five to several hundreds of millions of years. For a Rayleigh number of 107, we investigate how lithosphere yield stress variations (1050 MPa) and the presence of deep thermochemical heterogeneities favor small?scale (200500 km) and intermediate?scale (500104 km) dynamic topography by controlling the formation of small?scale convection and the number and distribution of subduction zones, respectively. The interplay between mantle convection and lithosphere dynamics generates a complex spatial and temporal pattern of dynamic topography consistent with constraints for Earth.
Abstract: The sinking remnant of a surface plate crosses and interacts with multiple boundaries in Earth's interior. Here, we specifically investigate the prominent dynamic interaction of the sinking plate portion with the upper-mantle transition zone and its corresponding surface elevation signal. We unravel, for the first time, that the collision of the sinking slab with the transition zone induces a sudden, dramatic downward tilt of the upper plate towards the subduction trench. Unraveling this crucial interaction was only possible thanks to state-of-the-art numerical modelling and post-processing. The new model that is introduced here to study the dynamically self-consistent temporal evolution of subduction features accurate subduction-zone topography, robust single-sided plate sinking, stronger plates close to laboratory values, an upper-mantle phase transition, and simple continents at a free surface. To distinguish the impact of the new physical model features, three different setups are used: the simplest model setup includes a basic high-viscosity lower mantle, the second adds a 660-km phase transition, and the third includes, additionally, a continental upper plate. Common to all models is the clear topographic signal upon slab-transition-zone interaction: the upper plate tilts abruptly towards the subduction trench by about 0.05° and over around 10 Ma. This dramatic increase in upper-plate tilt can be related to the slab-induced excitation of the high-viscosity lower mantle, which introduces a wider flow pattern. A large change in horizontal extent of inundation of up to 900 km is observed as a direct consequence of the upper-plate tilting. Such an abrupt variation in surface topography and inundation extent should be clearly visible in temporal records of large-scale surface elevation and might explain continental tilting as observed in Australia since the Eocene and North America during the Phanerozoic.
Abstract: As subducting plates reach the base of the upper mantle, some appear to flatten and stagnate, while others seemingly go through unimpeded. This variable resistance to slab sinking has been proposed to affect long-term thermal and chemical mantle circulation. A review of observational constraints and dynamic models highlights that neither the increase in viscosity between upper and lower mantle (likely by a factor 20-50) nor the coincident endothermic phase transition in the main mantle silicates (with a likely Clapeyron slope of -1 to -2 MPa/K) suffice to stagnate slabs. However, together the two provide enough resistance to temporarily stagnate subducting plates, if they subduct accompanied by significant trench retreat. Older, stronger plates are more capable of inducing trench retreat, explaining why backarc spreading and flat slabs tend to be associated with old-plate subduction. Slab viscosities that are ?2 orders of magnitude higher than background mantle (effective yield stresses of 100-300 MPa) lead to similar styles of deformation as those revealed by seismic tomography and slab earthquakes. None of the current transition-zone slabs seem to have stagnated there more than 60 m.y. Since modeled slab destabilization takes more than 100 m.y., lower-mantle entry is apparently usually triggered (e.g., by changes in plate buoyancy). Many of the complex morphologies of lower-mantle slabs can be the result of sinking and subsequent deformation of originally stagnated slabs, which can retain flat morphologies in the top of the lower mantle, fold as they sink deeper, and eventually form bulky shapes in the deep mantle.
Abstract: During flat subduction, material is scraped off the base of the continental mantle lithosphere, building a migrating keel. This testable mechanism for flat subduction recreates features of the Laramide orogeny.
Abstract: The fraction of radioactive heat production in Earth’s mantle to convective heat loss has decreased with the aging of Earth, as more of its nuclear fuel became spent and more of its heat was lost to space. Earth was therefore hotter in its past, but there is no consensus as to how much higher the mantle’s temperature was in early Earth compared to the present. This is an important problem to understand because it is expected to have imposed secular changes in the formation of oceanic lithosphere at ridges and its cycling at subduction zones (Herzberg and Rudnick, 2012; Foley, 2018). In a hotter early Earth, the ambient mantle melted more extensively, to make thicker basaltic oceanic crust and residual mantle peridotite, the latter of which was depleted in chemical elements that entered the magmas. Sometime later, the basaltic oceanic crust became hydrated by seawater, and it in turn melted to make silicic continental crust. As discussed in more detail below, this transformation led to the juxtaposition of continental crust on top of oceanic lithospheric mantle (Herzberg and Rudnick, 2012). The original “oceanic mantle lithosphere” is now called "continental mantle lithosphere" because it is located below continental crust in cratons. This hypothesis is explored by Servali and Korenaga (2018, p. 1047 in this issue of Geology), and is the reason why they entitle their paper an "oceanic origin of continental mantle lithosphere".
Abstract: Mechanisms of Precambrian orogeny and their contribution to the origin of ultrahigh temperature granulites, granite-greenstone terranes and net crustal growth remain debatable. Here, we use 2D numerical models with 150 °C higher mantle temperatures compared to present day conditions to investigate physical and petrological controls of Precambrian orogeny during forced continental plates convergence. Numerical experiments show that convergence between two relatively thin blocks of continental lithosphere with fertile mantle creates a short-lived cold collisional belt that later becomes absorbed by a long-lived thick and flat ultra-hot accretionary orogen with Moho temperatures of 700-1100 °C. The orogen underlain by hot partially molten depleted asthenospheric mantle spreads with plate tectonic rates towards the incoming lithospheric block. The accretionary orogeny is driven by delamination of incoming lithospheric mantle with attached mafic lower crust and invasion of the hot partially molten asthenospheric wedge under the accreted crust. A very fast convective cell forms atop the subducting slab, in which hot asthenospheric mantle rises against the motion of the slab and transports heat towards the moving orogenic front. Juvenile crustal growth during the orogeny is accompanied by net crustal loss due to the lower crust subduction. Stability of an ultra-hot orogeny is critically dependent on the presence of relatively thin and warm continental lithosphere with thin crust and dense fertile mantle roots subjected to plate convergence. Increased thickness of the continental crust and subcontinental lithospheric mantle, pronounced buoyancy of the lithospheric roots, and decreased mantle and continental Moho temperature favor colder and more collision-like orogenic styles with thick crust, reduced magmatic activity, lowered metamorphic temperatures, and decreased degree of crustal modification. Our numerical modeling results thus indicate that different types of orogens (cold, mixed-hot and ultra-hot) could be created at the same time in the Early Earth, depending on compositional and thermal structures of interacting continental blocks.
Geophysical Research Letters, Vol. 45, 24, pp. 13,298-13,305.
Mantle
subduction
Abstract: A subduction zone's mantle wedge can have a complex pattern of seismic anisotropy where the fast direction often rotates from trench?parallel close to the trench to trench?normal in the backarc. This pattern can be interpreted as induced by either 3?D trench?parallel flow or by the presence of water close to the trench. Almost all models so far favored the trench?parallel flow hypothesis, usually based on indirect or complementary indicators such as the evolution of geochemical signatures of volcanoes along the arc. Here we examine a seismic anisotropy observational signature that can be used to discriminate between the two explanations. The concept is defined using an interdisciplinary approach linking a direct modeling of the flow in the subduction wedge and a computation of seismic wave propagation in anisotropic media. We define a unique water?induced signature that is the presence of a “morph zone” characterized by a weak anisotropy and a decrease of seismic velocities. We apply the model to the Lau Basin where we find this predicted signature, demonstrating for the first time that water rather than trench?parallel flow is responsible for the observed anisotropy pattern there.
Abstract: The supercontinent cycle of episodic assembly and breakup of almost all continents on Earth is commonly considered the longest period variation to affect mantle convection. However, global zircon Hf isotopic signatures and seawater Sr isotope ratios suggest the existence of a longer-term variation trend that is twice the duration of the supercontinent cycle. Here we propose that since ?2 billion years ago the superocean surrounding a supercontinent, as well as the circum-supercontinent subduction girdle, survive every second supercontinent cycle. This interpretation is in agreement with global palaeogeography and is supported by variations in passive margin, orogen, and mineral deposit records that each exhibits both ?500-700 million years periodic signal and a 1000-1500 million years variation trend. We suggest that the supercontinent cycle is modulated by an assembly that alternates between dominantly extroversion after a more complete breakup, and dominantly introversion after an incomplete breakup of the previous supercontinent.
Abstract: Serpentinized oceanic mantle lithosphere is considered an important carrier of water and fluid-mobile elements, including halogens, into subduction zones. Seafloor serpentinite compositions indicate Cl, Br and I are sourced from seawater and sedimentary pore fluids, while F may be derived from hydrothermal fluids. Overall, the heavy halogens are expelled from serpentinites during the lizardite-antigorite transition. Fluorine, on the other hand, appears to be retained or may be introduced from dehydrating sediments and/or igneous rocks during early subduction. Mass balance calculations indicate nearly all subducted F is kept in the subducting slab to ultrahigh-pressure conditions. Despite a loss of Cl, Br and I from serpentinites (and other lithologies) during early subduction, up to 15% of these elements are also retained in the deep slab. Based on a conservative estimate for serpentinite thickness of the metamorphosed slab (500 m), antigorite serpentinites comprise 37% of this residual Cl, 56% of Br and 50% of I, therefore making an important contribution to the transport of these elements to the deep mantle.
Abstract: Across the entire mantle we interpret 94 positive seismic wave-speed anomalies as subducted lithosphere and associate these slabs with their geological record. We document this as the Atlas of the Underworld, also accessible online at www.atlas-of-the-underworld.org, a compilation comprising subduction systems active in the past ~ 300 Myr. Deeper slabs are correlated to older geological records, assuming no relative horizontal motions between adjacent slabs following break-off, using knowledge of global plate circuits, but without assuming a mantle reference frame. The longest actively subducting slabs identified reach the depth of ~ 2500 km and some slabs have impinged on Large Low Shear Velocity Provinces in the deepest mantle. Anomously fast sinking of some slabs occurs in regions affected by long-term plume rising. We conclude that slab remnants eventually sink from the upper mantle to the core-mantle boundary. The range in subduction-age versus - depth in the lower mantle is largely inherited from the upper mantle history of subduction. We find a significant depth variation in average sinking speed of slabs. At the top of the lower mantle average slab sinking speeds are between 10 and 40 mm/yr, followed by a deceleration to 10-15 mm/yr down to depths around 1600-1700 km. In this interval, in situ time-stationary sinking rates suggest deceleration from 20 to 30 mm/yr to 4-8 mm/yr, increasing to 12-15 mm/yr below 2000 km. This corroborates the existence of a slab deceleration zone but we do not observe long-term (> 60 My) slab stagnation, excluding long-term stagnation due to compositional effects. Conversion of slab sinking profiles to viscosity profiles shows the general trend that mantle viscosity increases in the slab deceleration zone below which viscosity slowly decreases in the deep mantle. This is at variance with most published viscosity profiles that are derived from different observations, but agrees qualitatively with recent viscosity profiles suggested from material experiments.
Geophysical Research Letters, Vol. 46, 1, pp. 179-188.
Mantle
subduction
Abstract: Eclogite mainly consists of pyrope?almandine?grossular garnet and sodium?rich pyroxene (omphacite) and is a key component of the Earth's upper mantle and oceanic crust. It plays an important role in the mantle convection. The lack of thermoelastic parameters of eclogitic garnets and omphacites hampers accurate modeling of eclogite density at deep?Earth pressure?temperature conditions. In this study, we obtained the thermoelastic parameters of natural eclogitic garnets and omphacites and then modeled the densities of high?Fe and low?Fe eclogites in the subducted oceanic crust and the normal upper mantle. In the upper mantle, eclogite enhances the slab subduction into the transition zone; however, the presence of the metastable low?Fe eclogite would promote the slab stagnation within the upper range of the transition zone. Additionally, eclogite can explain positive density anomalies at depths of 100-200 km of the upper mantle of Asia identified by seismic observations.
Geochemistry, Geophysics, Geosystems, Vol. 20, 2, pp. 1189-1201.
Mantle
subduction
Abstract: Tectonic plates plunge into the mantle at trenches, carrying water from the oceans. Some of this water may go down to the mantle transition zone between 410? and 660?km depth, where minerals have a large water storage capacity. In this study, we use observations of seismic anisotropy, the directional dependency of seismic wave speed, which is sensitive to the water content in the mantle transition zone. We find that the mantle transition zone beneath some subduction zones is drier than previously thought.
Abstract: Seismic anisotropy provides key information to map the trajectories of mantle flow and understand the evolution of our planet. While the presence of anisotropy in the uppermost mantle is well established, the existence and nature of anisotropy in the transition zone and uppermost lower mantle are still debated. Here we use three-dimensional global seismic tomography images based on a large dataset that is sensitive to this region to show the ubiquitous presence of anisotropy in the lower mantle beneath subduction zones. Whereas above the 660?km seismic discontinuity slabs are associated with fast SV anomalies up to about 3%, in the lower mantle fast SH anomalies of about 2% persist near slabs down to about 1,000-1,200?km. These observations are consistent with 3D numerical models of deformation from subducting slabs and the associated lattice-preferred orientation of bridgmanite produced in the dislocation creep regime in areas subjected to high stresses. This study provides evidence that dislocation creep may be active in the Earth’s lower mantle, providing new constraints on the debated nature of deformation in this key, but inaccessible, component of the deep Earth.
Abstract: Subduction-Transform Edge Propagator (STEP) faults are the locus of continual lithospheric tearing at slab edges, resulting in sharp changes in the lithospheric and crustal thickness and triggering lateral and/or near-vertical mantle flow. However, the mechanisms at the lithospheric mantle scale are still poorly understood. Here, we present the microstructural study of olivine-rich lherzolite, harzburgite and wehrlite mantle xenoliths from the Oran volcanic field (Tell Atlas, NW Algeria). This alkali volcanic field occurs along a major STEP fault responsible for the Miocene westward slab retreat in the westernmost Mediterranean. Mantle xenoliths provide a unique opportunity to investigate the microstructures in the mantle section of a STEP fault system. The microstructures of mantle xenoliths show a variable grain size ranging from coarse granular to fine-grained equigranular textures uncorrelated with modal variations. The major element composition of the mantle peridotites provides temperature estimates in a wide range (790-1165?°C) but in general, the coarse-grained and fine-grained peridotites suggest deeper and shallower provenance depth, respectively. Olivine grain size in the fine-grained peridotites depends on the size and volume fraction of the pyroxene grains, which is consistent with pinning of olivine grain growth by pyroxenes as second phase particles. In the coarse-grained peridotites, well-developed olivine crystal preferred orientation (CPO) is characterized by orthorhombic and [100]-fiber symmetries, and orthopyroxene has a coherent CPO with that of olivine, suggesting their coeval deformation by dislocation creep at high-temperature. In the fine-grained microstructures, along with the weakening of the fabric strength, olivine CPO symmetry exhibits a shift towards [010]-fiber and the [010]- and [001]-axes of orthopyroxene are generally distributed subparallel to those of olivine. These data are consistent with deformation of olivine in the presence of low amounts of melts and the precipitation of orthopyroxenes from a melt phase. The bulk CPO of clinopyroxene mimics that of orthopyroxene via a topotaxial relationship of the two pyroxenes. This observation points to a melt-related origin of most clinopyroxenes in the Oran mantle xenoliths. The textural and geochemical record of the peridotites are consistent with interaction of a refractory harzburgite protolith with a high-Mg# melt at depth (resulting in the formation of coarse-grained clinopyroxene-rich lherzolite and wehrlite), and with a low-Mg# evolved melt in the shallow subcontinental lithospheric mantle (forming fine-grained harzburgite). We propose that pervasive melt-peridotite reaction - promoted by lateral and/or near-vertical mantle flow associated with lithospheric tearing - resulted in the synkinematic crystallization of secondary lherzolite and wehrlite and played a key effect on grain size reduction during the operation of the Rif-Tell STEP fault. Melt-rock reaction and secondary formation of lherzolite and wehrlite may be widespread in other STEP fault systems worldwide.
Abstract: The burial and exhumation of continental crust during collisional orogeny exert a strong control on the dynamics of mountain belts and plateaus. Constraining the rates and style of exhumation of deeply buried crust has proven difficult due to complexities in the local geology and thermochronometric methods typically used. To advance this field, we applied trace?element and U?Pb laser ablation inductively coupled plasma mass spectrometry analyses to rutile from eclogite and amphibolite samples from the Western Gneiss Complex of Norway—an archetypal continental (ultra)high?pressure (UHP) terrane. Peak temperature and timing of midcrustal cooling were constrained for samples collected along a subduction? and exhumation?parallel transect, using Zr?in?rutile thermometry and U?Pb rutile geochronology, respectively. Peak temperatures decrease from 830 °C in the UHP domain to 730 °C at the UHP?HP transition, remain constant at 730 °C across most of the terrane, and decrease to 620 °C at the eclogite?out boundary. U?Pb results show that most of the terrane cooled through 500 °C at 380-375 Ma except for the lowest grade region, where cooling occurred approximately 20 million years earlier. The results indicate that exhumation was a two stage process, involving (1) flexural rebound and slab flattening at depth combined with foreland?directed extrusion, followed by (2) synchronous cooling below 500 °C across the, by then, largely flat?lying Western Gneiss Complex. The latter implies and requires relatively homogeneous mass removal across a large area, consistent with erosion of an overlying orogenic plateau. The Caledonides were at near?equatorial latitudes at the time. A Caledonian paleo?plateau thus may represent a so far unrecognized factor in Devonian and Carboniferous atmospheric circulation and climate forcing.
Geophysical Research Letters, Vol. 46, 7, pp. 3663-3675.
Mantle
subduction
Abstract: Separation of microcontinental blocks from their parent continent is usually attributed to abrupt relocation of concentrated extension from the mid?oceanic ridge to the adjacent continental margin. In the context of extensional passive margin evolution, previous extensive numerical and analog studies have revealed that hot upwelling mantle flow plays a key role in the mechanical weakening of the passive margin lithosphere needed to initiate a ridge jump. This, in turn, results in continental breakup and subsequent microcontinent isolation. However, the consequences of mantle plume impingement on the base of a moving lithospheric plate that is already involved into subduction are still unexplored quantitatively. Here we present the results of 3?D thermo?mechanical models showing that even in the context of induced plate motion (contractional boundary conditions), which are necessary to sustain continuous convergence, thermal and buoyancy effects of the mantle plume emplaced at the bottom of the continental part of the subducting plate are sufficient to initiate continental breakup and the subsequent opening of a new oceanic basin that separates the microcontinental block from the main body of the continent. With these models, we show that it is physically possible to form microcontinents in a convergent setting without the cessation of subduction.
Journal of the Geological Society, Vol. 176, pp. 348-374.
Mantle
carbon, subduction
Abstract: Carbon is subducted to depths where metamorphism liberates water-bearing fluids. The C-bearing fluids facilitate partial melting of the upper mantle, generating magmas that may erupt as arc volcanics. Degassing of the magmas releases CO2 and other volatile species to the atmosphere. Over geological time, this process contributes to the composition of the atmosphere and planetary habitability. Here I summarize the background needed to carry out theoretical geochemical modelling of fluids and fluid-rock interactions from surficial conditions into the upper mantle. A description of the general criteria for predicting equilibrium and non-equilibrium chemical reactions is followed by a summary of how the thermodynamic activities of species are related to measurable concentrations through standard states and activity coefficients. Specific examples at ambient conditions involving dilute water are detailed. The concept of aqueous speciation and how it can be calculated arises from this discussion. Next, I discuss how to calculate standard Gibbs free energies and aqueous activity coefficients at elevated temperatures and pressures. The revised Helgeson-Kirkham-Flowers equations of state are summarized and the revised predictive correlations for the estimation of equation of state coefficients in the Deep Earth Water (DEW) model are presented. Finally, the DEW model is applied to the solubility and speciation of aqueous aluminium.
Abstract: Inherited structural, compositional, thermal, and mechanical properties from previous tectonic phases can affect the deformation style of lithosphere entering a new stage of the Wilson cycle. When continental crust jams a subduction zone, the transition from subduction to extension can occur rapidly, as is the case following slab breakoff of the leading subducted oceanic slab. This study explores the extent to which geometric and physical properties of the subduction phase affect the subsequent deformation style and surface morphology of post subduction extensional systems. We focus on regions that transition rapidly from subduction to extension, retaining lithospheric heterogeneities and cold thermal structure inherited from subduction. We present numerical models suggesting that following failed subduction of continental crust (with or without slab breakoff), the extensional deformation style depends on the strength and dip of the preexisting subduction thrust. Our models predict three distinct extensional modes based on these inherited properties: (1) reactivation of the subduction thrust and development of a rolling?hinge detachment that exhumes deep crustal material in a domal structure prior to onset of an asymmetric rift; (2) partial reactivation of a low?angle subduction thrust, which is eventually abandoned as high?angle, “domino”?style normal faults cut and extend the crust above the inherited thrust; and (3) no reactivation of the subduction fault but instead localized rifting above the previous subduction margin as new rift?bounding, high?angle normal faults form. We propose that the first mode is well exemplified by the young, rapidly exhumed Dayman?Suckling metamorphic core complex that is exhuming today in Papua New Guinea.
Abstract: On contemporary Earth, subduction recycles mafic oceanic crust and associated volatile elements, creating new silicic continental crust in volcanic arcs. However, if the mantle was hotter in the Precambrian, the style of subduction, the depth of devolatilization and the formation of silicic continental crust may have been different. Consequently, the generation of the tonalite-trondhjemite-granodiorite (TTG) suite, which is characteristic of Archean crust, may not have been related to subduction. Here, we use a two-dimensional numerical magmatic-thermomechanical model to investigate intraoceanic subduction for contemporary mantle conditions and at higher mantle temperatures, as appropriate to the Precambrian. In each case, we characterize the thermal structure of the subducting plate and investigate magma compositions and production rates. We use these results to assess the potential growth of silicic continental crust associated with intraoceanic subduction at different mantle temperatures. For the Precambrian, in a set of experiments with ?T?=?150?K and decreasing subducting plate velocity, we find that the contemporary style of subduction was preceded by an arc-free regime dominated by rapid trench rollback and vigorous upwelling of asthenospheric mantle into the space created above the retreating slab. In this regime, formation of magmas by fluid-fluxed melting of the mantle is suppressed. Instead, decompression melting of upwelling asthenospheric mantle results in the widespread development of voluminous plateau-like basaltic magmas. In addition, retreating subduction at higher mantle temperature causes faster descent of the downgoing slab, leads to colder thermal gradients, similar to those associated with active subduction in the western Pacific today, and suppresses melting of the basaltic crust, limiting production of silicic (adakite-like) magmas. With increasing maturity of the subduction system, retreat of the subducting plate ceases, the role of decompression melting strongly decreases and fluid-fluxed melting of the mantle coupled with melting of the hydrated slab begins to produce basaltic and felsic arc volcanic rocks similar to those formed during contemporary subduction. In an additional series of individual experiments at various ?T, an increase of the mantle temperature above ?T?=?150?K leads to episodic and short-lived subduction accompanied by limited production of silicic continental crust. The results of our experimental study demonstrate that a hotter mantle in the Precambrian changes dramatically both the slab dynamics and the processes of magma generation and crustal growth associated with intraoceanic subduction zones. These changes may preclude growth of the early Precambrian silicic continental crust by processes that were dominantly similar to those associated with contemporary subduction.
Abstract: Kimberlite magmas are volatile-rich, potassic, and ultramafic, and they are host to most of the world’s diamond deposits. A continental-scale kimberlite magmatic belt (the central Cretaceous kimberlite corridor [CCKC]) is found in the interior of the North American continent. Parallel to and coeval with the CCKC, the Cretaceous Omineca magmatic belt (OMB) is located in the Cordilleran orogen. Cordilleran magmatism, including the OMB, is commonly explained through long-lived east-dipping subduction beneath the western margin of the continent. However, this does not explain the temporal and spatial relationships between the OMB and CCKC. We suggest that west-dipping subduction of North American lithosphere beneath the eastern side of the Cordillera explains both. In this model, subduction resulted in arc magmatism of the OMB. The contemporaneous CCKC was formed by extensional stress acting on the continent as it flexed upon entry into the trench. Using a semi-infinite elastic beam model, we show that flexure of a subducting continental plate (elastic thickness = 120 km) produces tensile stresses in the deep continental lithosphere, coincident with the location of the CCKC.
Abstract: The Earth’s crust-mantle boundary, the Mohorovi?i? discontinuity, has been traditionally considered to be the interface between the magnetic crust and the non-magnetic mantle1. However, this assumption has been questioned by geophysical observations2,3 and by the identification of magnetic remanence in mantle xenoliths4, which suggest mantle magnetic sources. Owing to their high critical temperatures, iron oxides are the only potential sources of magnetic anomalies at mantle depths5. Haematite (?-Fe2O3) is the dominant iron oxide in subducted lithologies at depths of 300 to 600 kilometres, delineated by the thermal decomposition of magnetite and the crystallization of a high-pressure magnetite phase deeper than about 600 kilometres6. The lack of data on the magnetic properties of haematite at relevant pressure-temperature conditions, however, hinders the identification of magnetic boundaries within the mantle and their contribution to observed magnetic anomalies. Here we apply synchrotron Mössbauer source spectroscopy in laser-heated diamond anvil cells to investigate the magnetic transitions and critical temperatures in Fe2O3 polymorphs7 at pressures and temperatures of up to 90 gigapascals and 1,300 kelvin, respectively. Our results show that haematite remains magnetic at the depth of the transition zone in the Earth’s mantle in cold or very cold subduction geotherms, forming a frame of deep magnetized rocks in the West Pacific region. The deep magnetic sources spatially correlate with preferred paths of the Earth’s virtual geomagnetic poles during reversals8 that might not reflect the geometry of the transitional field. Rather, the paths might be an artefact caused by magnetized haematite-bearing rocks in cold subducting slabs at mid-transition zone depths. Such deep sources should be taken into account when carrying out inversions of the Earth’s geomagnetic data9, and especially in studies of planetary bodies that no longer have a dynamo10, such as Mars.
Annual Reviews of Earth and Planetary Sciences, Vol. 47, pp. 41-66.
Mantle
subduction
Abstract: In this review, we address the current status of numerical modeling of the mantle transition zone and uppermost lower mantle, focusing on the hydration mechanism in these areas. The main points are as follows: (a) Slab stagnation and penetration may play significant roles in transporting the water in the whole mantle, and (b) a huge amount of water could be absorbed into the deep mantle to preserve the surface seawater over the geologic timescale. However, for further understanding of water circulation in the deep planetary interior, more mineral physics investigations are required to reveal the mechanism of water absorption in the lower mantle and thermochemical interaction across the core-mantle boundary region, which can provide information on material properties to the geodynamics community. Moreover, future investigations should focus on determining the amount of water in the early planetary interior, as suggested by the planetary formation theory of rocky planets. Moreover, the supplying mechanism of water during planetary formation and its evolution caused by plate tectonics are still essential issues because, in geodynamics modeling, a huge amount of water seems to be required to preserve the surface seawater in the present day and to not be dependent on an initial amount of water in Earth's system.
Abstract: Tectonic plate velocities predominantly result from a balance between the potential energy change of the subducting slab and viscous dissipation in the mantle, bending lithosphere, and slab-upper plate interface. A range of observations suggest that slabs may be weak, implying a more prominent role for plate interface dissipation than previously thought. The shallow thrust interface is commonly assumed to be weak due to an abundance of fluids and near-lithostatic pore fluid pressures, but little attention has been paid to the influence of the deeper, viscous interface. Here we show that the deep interface viscosity in subduction zones is strongly affected by the relative proportions of sedimentary to mafic rocks that are subducted to depth. Where sediments on the down-going plate are sparse, the deep interface is dominated by mafic lithologies that metamorphose to eclogites, which exhibit viscosities 1-2 orders of magnitude higher than the asthenospheric mantle, and reduce subduction plate speeds. In contrast, where sediments are abundant and subducted to depth, the deep interface viscosity is 1-2 orders of magnitude lower than the asthenospheric mantle, thus allowing significantly faster plate velocities. This correlation between subduction plate speed and deep sediment subduction may help explain dramatic accelerations (or decelerations) in convergence rates, such as the acceleration documented for India-Asia convergence during the mid-Cenozoic.
Abstract: Archean cratons consist of crustal granite-greenstone terrains (GGTs) coupled to roots of strong, buoyant cratonic lithospheric mantle (CLM). Although this association is unique to the Archean and formed from ca. 4.0 to 2.5 Ga, the origins of terrestrial cratons are debated. I propose that crustal plateaus, quasi-circular craton-like features (?1400-2400 km diameter, 0.5-4 km high), on Earth’s sister planet Venus might serve as analogs for Archean cratons. Crustal plateaus, which are isostatically supported by a compositionally controlled low-density root, host a distinctive surface called ribbon-tessera terrain. Ribbon-tessera also occurs as arcuate-shaped inliers in the Venus lowlands, widely interpreted as remnants of rootless crustal plateaus. Within each crustal plateau, surface ribbon-tessera terrain comprises a vast igneous province analogous to terrestrial GGTs, and the plateau root is analogous to CLM. Crustal plateaus and ribbon-tessera terrain collectively represent Venus’ oldest preserved features and surfaces, and they formed during an ancient period of globally thin lithosphere. To explain the linked features of crustal plateaus, a bolide impact hypothesis has been proposed in which a large bolide pierces ancient thin lithosphere, leading to massive partial melting in the sublithospheric mantle. In this model, melt escapes to the surface, forming an enormous lava pond, which evolves to form ribbon-tessera terrain; mantle melt residue forms a strong, resilient buoyant root, leading to plateau support and long-term stability of an individual crustal plateau. Building on the similarity of GGT-CLM and Venus crustal plateaus, I propose an exogenic hypothesis for Archean craton formation in which a large bolide pierces thin Archean lithosphere, causing localized high-temperature, high-fraction partial melting in the sublithospheric mantle; melt rises, forming an igneous province that evolves to form a GGT, and melt residue develops a complementary CLM. By this mechanism, Archean cratons may have formed in a spatially and temporally punctuated fashion at a time when large bolides showered Archean Earth.
Nature Communications, doi.org/10.1038/ s41467-019-11072-5 10p. Pdf
Europe, Italy
subduction
Abstract: Recycling of Earth’s crust through subduction and delamination contributes to mantle heterogeneity. Melt inclusions in early crystallised magmatic minerals record greater geochemical variability than host lavas and more fully reflect the heterogeneity of magma sources. To date, use of multiple isotope systems on small (300 ?m) melt inclusions was hampered by analytical limitations. Here we report the first coupled Sr-Nd-Pb isotope data on individual melt inclusions from potassium-rich lavas from neighbouring Quaternary volcanoes in central Italy and infer the presence of a previously unidentified ancient lower crustal component in the mantle. We suggest derivation from Variscan or older basement included in the upper mantle by either delamination, sediment recycling, subduction erosion and/or slab detachment processes during Cenozoic subduction and collision of the western Mediterranean. The capability to determine isotope ratios in individual melt inclusions permits the detection of distinctive mantle contaminants and can provide insights into how geodynamic processes affect subduction recycling.
Abstract: Many of the world's largest known REE deposits are associated with post-collisional alkaline-carbonatite magmatic complexes (e.g., the Minanning-Dechang belt, China). These systems are potassic to ultrapotassic in composition and contain LREE-dominated mineralisation associated with F and Ba-rich carbonatite breccias, carbonatite dykes and carbo-hydrothermal veins. They are typically emplaced through major shear zones during a period of 'relaxation' that postdates continental collision by up to 75 Ma. The subduction of sediment during continental collision is potentially a key control on the 'fertility' of the mantle source, and understanding the role of sediment is a crucial step towards better exploration models. However, the identification of sediment source components to alkaline systems has not been straightforward because their petrological complexity precludes traditional methods such as trace-element ratios and major-element modelling of crystal fractionation. We use a global database of Sr, Nd and Hf isotope compositions for alkaline and carbonatite systems, alongside geodynamic reconstructions to identify favourable source components for mineralisation and to provide direct information about the origin of the metals of interest. Subduction of shale and carbonate sequences is likely to introduce REE + HFSE and potentially mineralising ligands (F-, CO3 2-) into the mantle source for post-collisional alkaline systems; clastic sediments are poorer in these vital components. This research provides a framework through which the mineral exploration industry can identify tectonic environments that are predisposed to form REE mineralisation, providing regional-scale (100-1000 km) guidance especially for systems hidden beneath sedimentary cover.
Abstract: Subduction is a component of plate tectonics, which is widely accepted as having operated in a manner similar to the present-day back through the Phanerozoic Eon. However, whether Earth always had plate tectonics or, if not, when and how a globally linked network of narrow plate boundaries emerged are matters of ongoing debate. Earth's mantle may have been as much as 200-300 °C warmer in the Mesoarchean compared to the present day, which potentially required an alternative tectonic regime during part or all of the Archean Eon. Here we use a data set of the pressure (P), temperature (T), and age of metamorphic rocks from 564 localities that vary in age from the Paleoarchean to the Cenozoic to evaluate the petrogenesis and secular change of metamorphic rocks associated with subduction and collisional orogenesis at convergent plate boundaries. Based on the thermobaric ratio (T/P), metamorphic rocks are classified into three natural groups: high T/P type (T/P > 775 °C/GPa, mean T/P ~1105 °C/GPa), intermediate T/P type (T/P between 775 and 375 °C/GPa, mean T/P ~575 °C/GPa), and low T/P type (T/P < 375 °C/GPa, mean T/P ~255 °C/GPa). With reference to published thermal models of active subduction, we show that low T/P oceanic metamorphic rocks preserving peak pressures >2.5 GPa equilibrated at P-T conditions similar to those modeled for the uppermost oceanic crust in a wide range of active subduction environments. By contrast, those that have peak pressures <2.2 GPa may require exhumation under relatively warm conditions, which may indicate subduction of young oceanic lithosphere or exhumation during the initial stages of subduction. However, low T/P oceanic metamorphic rocks with peak pressures of 2.5-2.2 GPa were exhumed from depths where, in models of active subduction, the slab and overriding plate change from being decoupled (at lower P) to coupled (at higher P), possibly suggesting a causal relationship. In relation to secular change, the widespread appearance of low T/P metamorphism in the Neoproterozoic represents a “modern” style of cold collision and deep slab breakoff, whereas rare occurrences of low T/P metamorphism in the Paleoproterozoic may reveal atypical localized regions of cold collision. Low T/P metamorphism is not known from the Archean geological record, but the absence of blueschists in particular is unlikely to reflect secular change in the composition of the oceanic crust. In addition, the premise that the formation of lawsonite requires abnormally low thermal gradients and the postulate that oceanic subduction-related rocks register significantly lower maximum pressures than do continental subduction-related rocks, and imply different mechanisms of exhumation, are not supported. The widespread appearance of intermediate T/P and high T/P metamorphism at the beginning of the Neoarchean, and the subsequent development of a clear bimodality in tectono-thermal environments are interpreted to be evidence of the stabilization of subduction during a transition to a globally linked network of narrow plate boundaries and the emergence of plate tectonics.
Abstract: The subduction of sediment connects the surface nitrogen cycle to that of the deep Earth. To understand the evolution of nitrogen in the atmosphere, the behavior of nitrogen during the subduction and melting of subducted sediments has to be estimated. This study presents high-pressure experimental measurements of the partitioning of nitrogen during the melting of sediments at sub-arc depths. For quantitative analysis of nitrogen in minerals and glasses, we calibrated the electron probe micro-analyzer on synthetic ammonium feldspar to measure nitrogen concentrations as low as 500??g?g?1. Nitrogen abundances in melt and mica are used together with mass balance calculations to determine DN(Mica/Melt), DN(Fluid/Mica), and DN(Fluid/Melt). Calculated partition coefficients correspond to expected values for NH4+, which behaves similarly to Rb+ due to its nearly identical size. Nitrogen partitioning between fluid and melt (DN(Fluid/Melt)) and fluid and bulk residue (melt+mica) (DN(Fluid/Bulk)) increase linearly with temperature normalized to pressure. This linear relationship can be used to calculate DN(Fluid/Melt) and DN(Fluid/Bulk) for slab melts from 800 to 1200?°C following: and [nasty equation that did not copy]. We used these partition coefficients to quantify the amount of N recycled into the mantle as 50?±?6% of today's atmospheric N. Depending on the rate of mantle N degassing we calculated 4 different scenarios for atmospheric pN2 evolution. All 4 scenarios estimate pN2 to be 8-12% higher at the beginning of the Phanerozoic. These estimates diverge towards the past due to uncertainties in the mechanism and magnitude of N degassing from the mantle. Assuming degassing of N in the past was close to modern degassing rates from MORB, pN2 was up to 40% higher at the onset of plate tectonics at 3-4?Ga. However, degassing rates were probably higher than this: assuming 10× and 20× times higher rates at the onset of plate tectonics leads to pN2 within 20% of modern values. If N degassing from the mantle is increased to 40× the modern MORB rate, pN2 in the Archean would have been 50% lower than today's, which is in accordance with observations from paleoatmospheric studies.
American Mineralogist, Vol. 104, 8, pp. 1065-1082.
Mantle
subduction
Abstract: Subduction is a component of plate tectonics, which is widely accepted as having operated in a manner similar to the present-day back through the Phanerozoic Eon. However, whether Earth always had plate tectonics or, if not, when and how a globally linked network of narrow plate boundaries emerged are matters of ongoing debate. Earth's mantle may have been as much as 200-300 °C warmer in the Mesoarchean compared to the present day, which potentially required an alternative tectonic regime during part or all of the Archean Eon. Here we use a data set of the pressure (P), temperature (T), and age of metamorphic rocks from 564 localities that vary in age from the Paleoarchean to the Cenozoic to evaluate the petrogenesis and secular change of metamorphic rocks associated with subduction and collisional orogenesis at convergent plate boundaries. Based on the thermobaric ratio (T/P), metamorphic rocks are classified into three natural groups: high T/P type (T/P > 775 °C/GPa, mean T/P ~1105 °C/GPa), intermediate T/P type (T/P between 775 and 375 °C/GPa, mean T/P ~575 °C/GPa), and low T/P type (T/P < 375 °C/GPa, mean T/P ~255 °C/GPa). With reference to published thermal models of active subduction, we show that low T/P oceanic metamorphic rocks preserving peak pressures >2.5 GPa equilibrated at P-T conditions similar to those modeled for the uppermost oceanic crust in a wide range of active subduction environments. By contrast, those that have peak pressures <2.2 GPa may require exhumation under relatively warm conditions, which may indicate subduction of young oceanic lithosphere or exhumation during the initial stages of subduction. However, low T/P oceanic metamorphic rocks with peak pressures of 2.5-2.2 GPa were exhumed from depths where, in models of active subduction, the slab and overriding plate change from being decoupled (at lower P) to coupled (at higher P), possibly suggesting a causal relationship. In relation to secular change, the widespread appearance of low T/P metamorphism in the Neoproterozoic represents a “modern” style of cold collision and deep slab breakoff, whereas rare occurrences of low T/P metamorphism in the Paleoproterozoic may reveal atypical localized regions of cold collision. Low T/P metamorphism is not known from the Archean geological record, but the absence of blueschists in particular is unlikely to reflect secular change in the composition of the oceanic crust. In addition, the premise that the formation of lawsonite requires abnormally low thermal gradients and the postulate that oceanic subduction-related rocks register significantly lower maximum pressures than do continental subduction-related rocks, and imply different mechanisms of exhumation, are not supported. The widespread appearance of intermediate T/P and high T/P metamorphism at the beginning of the Neoarchean, and the subsequent development of a clear bimodality in tectono-thermal environments are interpreted to be evidence of the stabilization of subduction during a transition to a globally linked network of narrow plate boundaries and the emergence of plate tectonics.
Abstract: We show that the peripheral Pangea subduction zone closely followed a polar great circle. We relate it to the band of faster?than?average velocities in lowermost mantle. Both structures have an axis of symmetry in the equatorial plane. Assuming geologically long term stationarity of the deep mantle structure, we propose to use the axis of symmetry of Pangea to define an absolute reference frame. This reference frame is close to the slab remnants and NNR frames of reference but disagrees with hot spots based frames. We apply this model to the last 400 Myr. We show that a hemispheric supercontinent appeared as early as 400 Ma. However, at 400 Ma, the axis of symmetry was situated quite far south and progressively migrated within the equatorial plane that it reached at 300 Ma. From 300 to 110?100 Ma, it maintained its position within the equatorial plane. We propose that the stationarity of Pangea within a single hemisphere surrounded by subduction zones led to thermal isolation of the underlying asthenosphere and consequent heating as well as a large accumulation of hot plume material. We discuss some important implications of our analysis concerning the proposition that the succession of supercontinents and dispersed continents is controlled by an alternation from a degree one to a degree two planform.
Journal of Geophysical Research: Solid Earth, Vol. 124, 7, pp. 6504-6522.
Mantle
subduction
Abstract: Mantle wedge hybridization by crust?derived melt is a crucial mechanism responsible for arc lavas. However, how the melt?rock reactions proceed in the mantle wedge and affect melt compositions is poorly understood. Garnet peridotites from Jiangzhuang in the Sulu orogen (eastern China) host garnetite and pyroxenite veins formed by slab?mantle interactions at different melt/rock ratios. The Jiangzhuang peridotites consist mainly of garnet lherzolites and minor harzburgites and represent a fragment of the mantle wedge influenced by ultrahigh?pressure metamorphism (5.2-6.1 GPa) in the subduction channel. Petrography, major and trace element geochemistry, and in situ clinopyroxene Sr isotope values of the garnetite and pyroxenite veins reveal their derivation from interactions between mantle wedge peridotites and deeply subducted crust?derived melts. The two veins share a common metamorphic and metasomatic history and have similar mineral assemblages and compositions, enriched isotope signatures, and formation P?T conditions, indicating the same source for their reacting melts. The different mineral proportions and microtextures between the garnetite and pyroxenite veins are ascribed to different melt/rock ratios. The garnetite vein formed at relatively high melt/rock ratios (>1:1), which would likely produce hybrid slab melts with Mg?rich, high?silica adakitic signatures. In contrast, the pyroxenite vein formed at low melt/rock ratios (<1:1), and the expected hybrid slab melts would evolve into high?Mg andesites. Moreover, recycled heterogeneous garnetite and pyroxenite could contribute to the mantle sources of intraplate magmas. Therefore, slab?mantle interactions at different melt/rock ratios could be an important crustal input to lithological and geochemical heterogeneities in the mantle.
Contributions to Mineralogy and Petrology, Vol. 174, 22p. Pdf
Mantle
subduction, melting
Abstract: To investigate the first melts of the mantle wedge in subduction zones and their relationship to primitive magmas erupted at arcs, the compositions of low degree melts of hydrous garnet lherzolite have been experimentally determined at 3.2 GPa over the temperature range of 925-1150 °C. Two starting compositions with variable H2O contents were studied; a subduction-enriched peridotite containing 0.61% Na2O, 0.16 K2O% (wt%) with 4.2 wt% H2O added (Mitchell and Grove in Contrib Mineral Petrol 170:13, 2015) and an undepleted mantle peridotite (Hart and Zindler in Chem Geol 57:247-267, 1986) with 14.5% H2O added (Till et al. in Contrib Mineral Petrol 163:669-688, 2012). Saturating phases include olivine, orthopyroxene, clinopyroxene, garnet and rutile. Melting extent is tracked from near solidus (~?5 wt%) to 25 wt%, which is close to or beyond the point where clinopyroxene and garnet are exhausted. The beginning of melting is a peritectic reaction where 0.54 orthopyroxene?+?0.17 clinopyroxene?+?0.13 garnet react to produce 1.0 liquid?+?0.88 olivine. The melt production rate near the solidus is 0.1 wt% °C?1 and increases to 0.3 wt% °C?1 over the experimentally studied interval. These values are significantly lower than that observed for anhydrous lherzolite (~?1 wt% °C?1). When melting through this reaction is calculated for a metasomatized lherzolite source, the rare earth element characteristics of the melt are similar to melts of an eclogite, as well as those observed in many subduction zone magmas. Moreover, since rutile is stable up to?~?8 wt% melting, the first melts of a hydrous lherzolite source could also show strong high field strength element depletions as is observed in many subduction zone lavas. The silicate melts measured at the lowest temperatures and melting extents (10 wt%) are high silica andesites (56-60 wt% SiO2) and contain very low Ca/Al and high alkalis. These deep low degree andesitic melts, if added to experimentally produced hydrous liquids from melting (20-25 wt%) of harzburgite residues at shallow pressures (1.0-1.2 GPa, Mitchell and Grove in Contrib Mineral Petrol 170:13, 2015), can match the compositional characteristics of primitive natural basaltic andesite and magnesian andesite lavas found globally in arcs. In addition to a silicate melt phase, there is a small amount of silicate dissolved in the H2O supercritical fluid that coexists with the silicate liquid and solids in our experiments. The composition of this fluid is in equilibrium with the Mg-rich minerals and it is granitic. The results presented here are used to develop a model for producing hydrous arc magmas. We hypothesize that mantle wedge melting produced by the flux of hydrous fluid from the slab occurs over a range of depths that begins at the base of the mantle wedge and ends at shallow mantle depths. These melts ascend and remain isolated until they mix in the shallow, hottest part of the mantle wedge. In this melting scenario, the metasomatic “slab melt” contributions to arc magmas is small (~?5 wt%), but its effect on the alkali, REE and incompatible trace element budget of the derivative magmas is large and able to reproduce the trace elemental characteristics of the primitive andesites. Higher proportions of slab or sediment melt do not resemble primitive high magnesian arc andesites and basaltic andesites.
Geochimica et Cosmochimica Acta, in press available. 45p.
Mantle
subduction
Abstract: Archean peridotite xenoliths in the ?2.52 Ga Zhulagou diorite (Yinshan Block, North China Craton) show chemical and Li isotopic evidence for metasomatism above an ancient subduction zone. The peridotite xenoliths are composed of olivine + orthopyroxene + amphibole + phlogopite + serpentine. The peridotite xenoliths have low whole-rock Mg# (80-81) and low Mg# (81-84) in olivine, indicating that they are cumulates that formed near the crust-mantle boundary. Petrological observations, mineral trace element data and isotopic ages show that the sequence of hydrous minerals is amphibole-serpentine-phlogopite. SIMS U-Pb dating of zircon from peridotites yielded an upper intercept age at ?2.53 Ga, and a U-Pb lower intercept age at ?1.8 Ga. The age of ?2.53 Ga is interpreted to date the crystallization of zircon from the metasomatized mantle melt that formed the Zhulagou cumulate peridotite. Rb-Sr mineral isochrons date phlogopite formation at ?1760 Ma, consistent with the lower intercept age of zircon. Pargasitic amphibole from the Zhulagou peridotites has fractionated REE, pronounced depletions of Nb, Ta, Zr and Ti, and heavy ?7Li (?+14‰) and light ?11B (?-11‰). Combined with slightly depleted mantle whole rock ?Nd (?+1.3) and high zircon ?18O (+5.6 to +7.0‰), the amphibole composition reflects that the peridotite xenoliths formed from melts that carried the geochemical and isotopic fingerprint typical for a metasomatized mantle wedge above a subduction zone. The Zhulagou peridotite xenoliths have the highest ?7Li values (?+12‰) recorded in Archean peridotites. Isotopically heavy Li and light B in olivine, orthopyroxene, and amphibole from the peridotite xenoliths show that Li and B may decouple during partial melting or fluid release from the subducted slab. The decoupling of Li and B may have a variety of reasons, including different host minerals for Li and B in the source and different protoliths in the subducted slab. The Li and B isotopic record on the recycling of ancient material demonstrates that modern-style subduction operated already in the late Archean.
Abstract: At subduction zones, most diamonds form by carbon saturation in hydrous fluids released from lithospheric plates on equilibration with mantle rocks. Although organic molecules are predicted among dissolved species which are the source for carbon in diamonds, their occurrence is not demonstrated in nature, and the physical model for crustal diamond formation is debated. Here, using Raman microspectroscopy, I determine the structure of carbon-based phases inside fluid inclusions in diamond-bearing rocks from the Alps. The results provide direct evidence that diamond surfaces are coated by sp2-, and sp3-bonded amorphous carbon and functional groups of carboxylic acids (e.g., carboxyl, carboxylate, methyl, and methylene), indicating the geosynthesis of organic compounds in deep hydrous fluids. Moreover, this study suggests diamond nucleation via metastable molecular precursors. As a possible scenario, with carbon saturation by reduction of carboxylate groups, I consider tetrahedral H-terminated C groups as templates for the growth of sp3-structured carbon.
Geochemistry, Geophysics, Geosystems, Vol. 20, pp. 3268-3288.
Mantle
subduction
Abstract: Molnar and England (1990, https://doi.org/10.1029/JB095iB04p04833) introduced equations using a semianalytical approach that approximate the thermal structure of the forearc regions in subduction zones. A detailed new comparison with high?resolution finite element models shows that the original equations provide robust predictions and can be improved by a few modifications that follow from the theoretical derivation. The updated approximate equations are shown to be quite accurate for a straight?dipping slab that is warmed by heat flowing from its base and by shear heating at its top. The approximation of radiogenic heating in the crust of the overriding plate is less accurate but the overall effect of this heating mode is small. It is shown that the previous and updated approximate equations become increasingly inaccurate with decreasing thermal parameter and increasing variability of slab dip. It is also shown that the approximate equations cannot be extrapolated accurately past the brittle?ductile transition. Conclusions in a recent paper (Kohn et al., 2018, https://doi.org/10.1073/pnas.1809962115) that modest amount of shear heating can explain the thermal conditions of past subduction from the exhumed metamorphic rock record are invalid due to a number of compounding errors in the application of the Molnar and England (1990, https://doi.org/10.1029/JB095iB04p04833) equations past the brittle?ductile transition. The use of the improved approximate equations is highly recommended provided their limitations are taken into account. For subduction zones with variable dip and/or low thermal parameter finite element modeling is recommended.
Nature Communications, doi.org/10.1038/ s41467-019-12984-y 8p. Pdf
Mantle
subduction
Abstract: At subduction zones, most diamonds form by carbon saturation in hydrous fluids released from lithospheric plates on equilibration with mantle rocks. Although organic molecules are predicted among dissolved species which are the source for carbon in diamonds, their occurrence is not demonstrated in nature, and the physical model for crustal diamond formation is debated. Here, using Raman microspectroscopy, I determine the structure of carbon-based phases inside fluid inclusions in diamond-bearing rocks from the Alps. The results provide direct evidence that diamond surfaces are coated by sp2-, and sp3-bonded amorphous carbon and functional groups of carboxylic acids (e.g., carboxyl, carboxylate, methyl, and methylene), indicating the geosynthesis of organic compounds in deep hydrous fluids. Moreover, this study suggests diamond nucleation via metastable molecular precursors. As a possible scenario, with carbon saturation by reduction of carboxylate groups, I consider tetrahedral H-terminated C groups as templates for the growth of sp3-structured carbon.
Nature Communications, doi.org/10.1038/s41467-019-11065-4 pdf 10p.
Mantle
subduction
Abstract: Ore deposits are loci on Earth where energy and mass flux are greatly enhanced and focussed, acting as magnifying lenses into metal transport, fractionation and concentration mechanisms through the lithosphere. Here we show that the metallogenic architecture of the lithosphere is illuminated by the geochemical signatures of metasomatised mantle rocks and post-subduction magmatic-hydrothermal mineral systems. Our data reveal that anomalously gold and tellurium rich magmatic sulfides in mantle-derived magmas emplaced in the lower crust share a common metallogenic signature with upper crustal porphyry-epithermal ore systems. We propose that a trans-lithospheric continuum exists whereby post-subduction magmas transporting metal-rich sulfide cargoes play a fundamental role in fluxing metals into the crust from metasomatised lithospheric mantle. Therefore, ore deposits are not merely associated with isolated zones where serendipitous happenstance has produced mineralisation. Rather, they are depositional points along the mantle-to-upper crust pathway of magmas and hydrothermal fluids, synthesising the concentrated metallogenic budget available.
Abstract: Understanding the dominant force responsible for supercontinent breakup is crucial for establishing Earth's geodynamic evolution that includes supercontinent cycles and plate tectonics. Conventionally, two forces have been considered: the push by mantle plumes from the sub-continental mantle which is called the active force for breakup, and the dragging force from oceanic subduction retreat which is called the passive force for breakup. However, the relative importance of these two forces is unclear. Here we model the supercontinent breakup coupled with global mantle convection in order to address this question. Our global model features a spherical harmonic degree-2 structure, which includes a major subduction girdle and two large upwelling (superplume) systems. Based on this global mantle structure, we examine the distribution of extensional stress applied to the supercontinent by both sub-supercontinent mantle upwellings and subduction retreat at the supercontinent peripheral. Our results show that: (1) at the center half of the supercontinent, plume push stress is ?3 times larger than the stress induced by subduction retreat; (2) an average hot anomaly of no higher than 50 K beneath the supercontinent can produce a push force strong enough to cause the initialization of supercontinent breakup; (3) the extensional stress induced by subduction retreat concentrates on a ?600 km wide zone on the boundary of the supercontinent, but has far less impact to the interior of the supercontinent. We therefore conclude that although circum-supercontinent subduction retreat assists supercontinent breakup, sub-supercontinent mantle upwelling is the essential force.
Abstract: Reactions involving carbon in the deep Earth have limited manifestations on Earth's surface, yet they have played a critical role in the evolution of our planet. The metal-silicate partitioning reaction promoted carbon capture during Earth's accretion and may have sequestered substantial carbon in Earth's core. The freezing reaction involving iron-carbon liquid could have contributed to the growth of Earth's inner core and the geodynamo. The redox melting/freezing reaction largely controls the movement of carbon in the modern mantle, and reactions between carbonates and silicates in the deep mantle also promote carbon mobility. The 10-year activity of the Deep Carbon Observatory has made important contributions to our knowledge of how these reactions are involved in the cycling of carbon throughout our planet, both past and present, and has helped to identify gaps in our understanding that motivate and give direction to future studies.
Geochemistry, Geophysics, Geosystems, Vol. 20, 1, pp. 5667-5690.
Mantle
subduction
Abstract: The amount of H2O and CO2 that is carried into deep mantle by subduction beyond subarc depths is of fundamental importance to the deep volatile cycle but remains debated. Given the large uncertainties surrounding the spatio?temporal pattern of fluid flow and the equilibrium state within subducting slabs, a model of H2O and CO2 transport in slabs should be balanced between model simplicity and capability. We construct such a model in a two?part contribution. In this Part I of our contribution, thermodynamic parameterization is performed for the devolatilization of representative subducting materials—sediments, basalts, gabbros, peridotites. The parameterization avoids reproducing the details of specific devolatilization reactions, but instead captures the overall behaviors of coupled (de)hydration and (de)carbonation. Two general, leading?order features of devolatilization are captured: (1) the released volatiles are H2O?rich near the onset of devolatilization; (2) increase of the ratio of bulk CO2 over H2O inhibits overall devolatilization and thus lessens decarbonation. These two features play an important role in simulation of volatile fractionation and infiltration in thermodynamically open systems. When constructing the reactive fluid flow model of slab H2O and CO2 transport in the companion paper Part II, this parameterization can be incorporated to efficiently account for the open?system effects of H2O and CO2 transport.
Abstract: The strength of subducted slabs in the mantle transition zone influences the style of mantle convection. Intense deformation is observed particularly in relatively old subducted slabs in the deep mantle transition zone. Understanding the cause of this regional and depth variation in slab deformation requires constraint of the rheological properties of deep mantle materials. Here, we report results of in situ deformation experiments during the olivine to ringwoodite phase transformation, from which we infer the deformation process under the conditions of cold slabs deep in the mantle transition zone. We find that newly transformed fine-grained ringwoodite deforms by diffusion creep and that its strength is substantially smaller than that of coarser-grained minerals but increases with time. Scaling analysis, based on a model of transformation kinetics and grain-size evolution during a phase transformation, suggests that a cold slab will be made of a mixture of weak, fine-grained and strong, coarse-grained materials in the deep transition zone, whereas a warm slab remains strong because of its large grain size. We propose that this temperature dependence of grain size may explain extensive deformation of cold slabs in the deep transition zone but limited deformation of relatively warm slabs.
Annual Review of Earth and Planetary Sciences, Vol. 47, pp. 41-66.
Mantle
subduction, water
Abstract: In this review, we address the current status of numerical modeling of the mantle transition zone and uppermost lower mantle, focusing on the hydration mechanism in these areas. The main points are as follows: (a) Slab stagnation and penetration may play significant roles in transporting the water in the whole mantle, and (b) a huge amount of water could be absorbed into the deep mantle to preserve the surface seawater over the geologic timescale. However, for further understanding of water circulation in the deep planetary interior, more mineral physics investigations are required to reveal the mechanism of water absorption in the lower mantle and thermochemical interaction across the core-mantle boundary region, which can provide information on material properties to the geodynamics community. Moreover, future investigations should focus on determining the amount of water in the early planetary interior, as suggested by the planetary formation theory of rocky planets. Moreover, the supplying mechanism of water during planetary formation and its evolution caused by plate tectonics are still essential issues because, in geodynamics modeling, a huge amount of water seems to be required to preserve the surface seawater in the present day and to not be dependent on an initial amount of water in Earth's system.
Abstract: This paper presents and discusses petrological observations from high- to ultrahigh-pressure (U)HP metamorphic terrains in relation to existing geophysical and numerical models for subduction and exhumation processes in orogenic belts. The interpretations are mostly based on observations from gneiss terrains bearing abundant bodies mafic (meta-)eclogite and ultramafic garnet peridotite and pyroxenite, exposed in collisional orogens. The inclusions and compositional zoning of minerals are considered to be first order information that is needed to constrain PT paths of HP-UHP rocks and reconstruct the related geodynamic models for subduction and exhumation of crustal and mantle rocks. The Bohemian Massif of the European Variscides is used as the basis for a model example to explain these processes, but (U)HP rocks from various other terrains are taken into consideration to discuss available PT paths in relation to proposed subduction and exhumation rates of (U)HP rocks based on geophysical and geochronological data. Primarily information used in this respect include textural relations and preserved prograde zoning in minerals from many (U)HP rocks, which reveal that a relatively cool geothermal gradient typical of subduction zones tended to prevail during the prograde and peak pressure segments of PT paths prior to initiation of exhumation and may have continued, even with cooling, if exhumation rates were rapid. The commonly applied interpretation of isothermal decompression during exhumation is critically appraised, considering whether a simple thermal relaxation (and radiogenic heating) during exhumation is responsible for formation of post-peak pressure, retrograde mineral assemblages and textures observed in (U)HP rocks. We go on to consider whether this can satisfactorily explain the often pervasive medium-pressure, high-temperature metamorphic re-equilibration of (U)HP rocks or whether an additional, external source of heat is a better explanation. We conclude that the commonly observed high-temperature metamorphic overprint exhibited by (U)HP rocks occurs mostly after rocks have been exhumed from the subduction channel and have reached normal crustal positions, when mantle upwelling resulting from slab breakoff (delamination) or slab rollback takes place at the onset of continent-continent collision. We also explore contrasting PT trajectories for mantle rocks that have been entrained into crustal material during their subduction or exhumation; PT paths of mantle and subducted crustal rocks tend to converge as mantle rocks impinge upon the cooler subduction zone and, once entrained, share a common evolution that depends on the exhumation mechanism and rate. Considering all of the data presented in this work we conclude that the diverse, polyphase metamorphic evolution exhibited by (U)HP terrains, embodied in the PT paths of HP and UHP rocks, has important consequences for reconstructing their changing thermal regimes and provides important constraints for geodynamic models involving subduction and the transition to collision.
Annual Review of Earth and Planetary Sciences, Vol. 48, 30p. pdf
Mantle
subduction, metamorphism
Abstract: If we accept that a critical condition for plate tectonics is the creation and maintenance of a global network of narrow boundaries separating multiple plates, then to argue for plate tectonics during the Archean requires more than a local record of subduction. A case is made for plate tectonics back to the early Paleoproterozoic, when a cycle of breakup and collision led to formation of the supercontinent Columbia, and bimodal metamorphism is registered globally. Before this, less preserved crust and survivorship bias become greater concerns, and the geological record may yield only a lower limit on the emergence of plate tectonics. Higher mantle temperature in the Archean precluded or limited stable subduction, requiring a transition to plate tectonics from another tectonic mode. This transition is recorded by changes in geochemical proxies and interpreted based on numerical modeling. Improved understanding of the secular evolution of temperature and water in the mantle are key targets for future research. 1) Higher mantle temperature in the Archean precluded or limited stable subduction, requiring a transition to plate tectonics from another tectonic mode. 2) Plate tectonics can be demonstrated on Earth since the early Paleoproterozoic (since c. 2.2 Ga), but before the Proterozoic Earth's tectonic mode remains ambiguous. 3) The Mesoarchean to early Paleoproterozoic (3.2-2.3 Ga) represents a period of transition from an early tectonic mode (stagnant or sluggish lid) to plate tectonics. 4) The development of a global network of narrow boundaries separating multiple plates could have been kick-started by plume-induced subduction.
Annual Review of Earth and Planetary Sciences, Vol. 47, pp. 41-66. pdf
Mantle
water, subduction
Abstract: In this review, we address the current status of numerical modeling of the mantle transition zone and uppermost lower mantle, focusing on the hydration mechanism in these areas. The main points are as follows: (a) Slab stagnation and penetration may play significant roles in transporting the water in the whole mantle, and (b) a huge amount of water could be absorbed into the deep mantle to preserve the surface seawater over the geologic timescale. However, for further understanding of water circulation in the deep planetary interior, more mineral physics investigations are required to reveal the mechanism of water absorption in the lower mantle and thermochemical interaction across the core-mantle boundary region, which can provide information on material properties to the geodynamics community. Moreover, future investigations should focus on determining the amount of water in the early planetary interior, as suggested by the planetary formation theory of rocky planets. Moreover, the supplying mechanism of water during planetary formation and its evolution caused by plate tectonics are still essential issues because, in geodynamics modeling, a huge amount of water seems to be required to preserve the surface seawater in the present day and to not be dependent on an initial amount of water in Earth's system. 1) Slab stagnation and penetration of the hydrous lithosphere are essential for understanding the global-scale material circulation. 2) Thermal feedback caused by water-dependent viscosity is a main driving mechanism of water absorption in the mantle transition zone and uppermost lower mantle. 3) The hydrous state in the early rocky planets remains to be determined from cosmo- and geochemistry and planetary formation theory. 4) Volatile cycles in the deep planetary interior may affect the evolution of the surface environment.
Lithos, doi.org/10.101016/j. lithos.2020.105460, 55p. Pdf
United States, Colorado
subduction
Abstract: The zonation of trace elements in Cr-pyrope has been analyzed to investigate histories recorded by garnet in cool mantle and implications for processes in the mantle wedge during low-angle subduction. The garnets are from the Navajo Volcanic Field (NVF) of the Colorado Plateau in the southwestern United States. Their host rocks were emplaced between 30 and 24 Ma, not long after low-angle subduction of the Farallon plate and the Laramide orogeny. Twenty-seven of the 31 garnets are discrete fragments from diatremes of serpentinized ultramafic microbreccia (SUM). Many Ni and Mn values in those grains yield temperatures in and below the range 800 °C to 600 °C, consistent with inclusions of chlorite and other hydrous minerals. Gradients of Na, Ti, Mn, Ni, Y, REE, and other elements are common in the discrete grains. Diffusion calculations constrain possible histories of a garnet with growth zonation of Y and REE and retrograde zonation of Mn and Ni. The garnet must have formed and resided at temperatures below about 900 °C, and that zonation is unlikely to have been preserved from the Proterozoic. The Mn gradient is consistent with cooling during the time period of low-angle Farallon subduction. Many of the discrete garnets grew from sources metasomatically enriched in Li, Na, Ti, and Zr, and they have distinctive high Li/Na. The high Li/Na is a likely consequence of Li metasomatism in the mantle wedge. Nb values in the Cr-pyrope fragments from northern SUM diatremes are lower than those of most garnets from kimberlites, but the low values are similar to and lower than those of pyrope in peridotite emplaced in orogenic belts; those northern diatremes also host lawsonite eclogite and garnetite. Most but not all relatively Cr-rich discrete garnets (4 to 8 wt% Cr2O3) record enrichment in LREE but little enrichment in Ti and Zr; LREE enrichment by an aqueous fluid is consistent with chlorite inclusions in two of these garnets. Oriented lamellae of rutile and other phases are common, and one texture may record fluid interactions associated with interface-coupled dissolution-reprecipitation. Garnets from the four minette-hosted peridotites record markedly hotter temperatures than do the SUM-hosted grains, and they lack the high Li/Na and low Nb. Gradients of trace elements in minette-hosted garnets record metasomatism shortly before crustal emplacement, but that metasomatism resulted in a decrease of Nb and had little effect on Li/Na. The minette-hosted garnet peridotites may represent mantle residual from Proterozoic accretion of the Colorado Plateau. In contrast, at least some of the SUM-hosted garnets probably grew in rock eroded from the forearc wedge by low-angle subduction and emplaced below the NVF.
Journal of Geophysical Research: Solid Earth, Vol. 125, e2019JB017884.
Mantle, South America
subduction
Abstract: Nazca subduction beneath South America is one of our best modern examples of long?lived ocean?continent subduction on the planet, serving as a foundation for our understanding of subduction processes. Within that framework, persistent heterogeneities at a range of scales in both the South America and Nazca plates is difficult to reconcile without detailed knowledge of the subducted Nazca slab structure. Here we use teleseismic travel time residuals from >1,000 broadband and short?period seismic stations across South America in a single tomographic inversion to produce the highest?resolution contiguous P wave tomography model of the subducting slab and surrounding mantle beneath South America to date. Our model reveals a continuous trench?parallel fast seismic velocity anomaly across the majority of South America that is consistent with the subducting Nazca slab. The imaged anomaly indicates a number of robust features of the subducted slab, including variable slab dip, extensive lower mantle penetration, slab stagnation in the lower mantle, and variable slab amplitude, that are incorporated into a new, comprehensive model of the geometry of the Nazca slab surface to ~1,100 km depth. Lower mantle slab penetration along the entire margin suggests that lower mantle slab anchoring is insufficient to explain along strike upper plate variability while slab stagnation in the lower mantle indicates that the 1,000 km discontinuity is dominant beneath South America.
Abstract: A novel structure of calcium orthocarbonate Ca2CO4-Pnma has been discovered using ab-initio crystal structure prediction methods (AIRSS and USPEX) based on the density functional theory. This phase appears above 13 GPa and remains stable to at least 50 GPa and 2000 K according to the calculations within quasi-harmonic approximation. Thus, the discovered phase can be stable at Earth's transition zone and lower mantle P-T conditions. The carbon atoms in the discovered phase are presented in 4-fold coordination, and its structure is similar to the high-pressure and high-temperature ?'H-Ca2SiO4 phase.
Abstract: Subduction diamonds represent the sequestration of carbon from fluids released from lithospheric plates at mantle depths. In deep fluids, besides reactive molecular species (e.g., CO2), inorganic, and organic aqueous ionic species have been proposed as a source of carbon in diamonds (Sverjensky et al., 2014). Unequivocal signatures of organic species, however, have not been found, neither a unified model for diamond nucleation and growth has been proposed. Here, I use Raman microspectroscopy to determine the structure and composition of carbon-based phases precipitated inside diamond-bearing fluid inclusions in metamorphic rocks from the Alps to reveal the spontaneous products on carbon saturation in deep Earth’s aqueous fluids. I find that microand nano-sized diamonds are coated by sp2-, and sp3-bonded amorphous carbon that shows Raman modes of attached organic functional group structures (e.g., carboxyl, carboxylate, methyl, and methylene). Present data suggest that decomposition of complex carboxylic acids can induce diamond nucleation on the reduction of the carboxyl groups, whereas sp3-bonded radicals can create structural intermediates allowing diamond growth (Frezzotti, 2019). This formation mechanism is consistent with nucleation models via metastable molecular precursors (Gebbie et al., 2018). The present study provides direct evidence that, deep within the Earth, dissolved inorganic carbon can spontaneously evolve to organic species in the absence of biologically catalyzed processes. Results suggest that the Earth’s interior should be considered as a favorable environment for the origin of prebiotic organic compounds.
Physics of the Earth and Planetary Interiors, Vol. 305, 13p. Pdf
Mantle
subduction
Abstract: The composition of Archean volcanic crust can be characterized by a higher Mg/Si ratio than modern mid-ocean ridge basalt (MORB), because of the higher degree melting from the warmer mantle in the Archean. Although modern MORB may become less dense than the surrounding mantle beneath the mantle transition zone (MTZ), the Mg-rich composition of Archean volcanic crust may result in the different density, and therefore different sinking behavior near the MTZ. In order to understand the compositional effect of Archean volcanic crust on the sinking behaviors and the scale of mantle mixing in the Archean, we investigated the mineralogy and density of Archean volcanic crust near the MTZ (470-910 km-depth). We conducted experiments at 19-34 GPa and 1400-2400 K using the laser-heated diamond anvil cell (LHDAC) combined with in-situ X-ray diffraction (XRD). The in-situ XRD and the chemical analysis revealed that Archean volcanic crust forms garnet and ringwoodite (84 and 16 vol%, respectively), which gradually transforms to Brg and CaPv (82 and 18 vol%, respectively) at 23-25 GPa and 1800 K. Our in-situ XRD experiments allowed us to measure the volumes of stable phases and to estimate their densities at high pressure and temperature. The results suggest that Archean volcanic crust maintains greater density than the pyrolitic mantle in the Archean regardless of temperature at 20-34 GPa (570-850 km-depth), promoting further sinking into the deeper mantle in the Archean. We also considered the density of the subducting slab in the Archean. The density model showed that the subducting slab is still denser or at least equally dense as the surrounding pyrolitic mantle in the Archean.
Abstract: A novel structure of calcium orthocarbonate Ca2CO4-Pnma has been discovered using ab-initio crystal structure prediction methods (AIRSS and USPEX) based on the density functional theory. This phase appears above 13 GPa and remains stable to at least 50 GPa and 2000 K according to the calculations within quasi-harmonic approximation. Thus, the discovered phase can be stable at Earth's transition zone and lower mantle P-T conditions. The carbon atoms in the discovered phase are presented in 4-fold coordination, and its structure is similar to the high-pressure and high-temperature ?'H-Ca2SiO4 phase.
Tectonics, in press available, e2020TC006063 39p. Pdf
Global
subduction
Abstract: Microcontinents and continental fragments are small pieces of continental crust that are surrounded by oceanic lithosphere. Although classically associated with passive margin formation, here we present several preserved microcontinents and continental fragments associated with subduction systems. They are located in the Coral Sea, South China Sea, central Mediterranean and Scotia Sea regions and a ‘proto?microcontinent’, in the Gulf of California. Reviewing the tectonic history of each region and interpreting a variety of geophysical data allows us to identify parameters controlling the formation of microcontinents and continental fragments in subduction settings. All these tectonic blocks experienced long, complex tectonic histories with an important role for developing inherited structures. They tend to form in back?arc locations and separate from their parent continent by oblique or rotational kinematics. The separated continental pieces and associated marginal basins are generally small and formation is quick (<50 Myr). Microcontinents and continental fragments formed close to large continental masses tend to form faster than those created in systems bordered by large oceanic plates. A common triggering mechanism for formation is difficult to identify, but seems to be linked with rapid changes of complex subduction dynamics. The young ages of all contemporary pieces found in situ suggest that microcontinents and continental fragments in these settings are short lived. Although presently the amount of in?situ subduction?related microcontinents is meagre (an area of 0.56% and 0.28% of global, non?cratonic, continental crustal area and crustal volume respectively), through time microcontinents contributed to terrane amalgamation and larger continent formation.
Abstract: The chemical and physical processes operating during subduction-zone metamorphism can profoundly influence the cycling of elements on Earth. Deep-Earth carbon (C) cycling and mobility in subduction zones has been of particular recent interest to the scientific community. Here, we present textural and geochemical data (CO, Sr isotopes and bulk and in-situ trace element concentrations) for a suite of ophicarbonate rocks (carbonate-bearing serpentinites) metamorphosed over a range of peak pressure-temperature (P-T) conditions together representing a prograde subduction zone P-T path. These rocks, in order of increasing peak P-T conditions, are the Internal Liguride ophicarbonates (from the Bracco unit, N. Apennines), pumpellyite- and blueschist-facies ophicarbonates from the Sestri-Voltaggio zone (W. Ligurian Alps) and the Queyras (W. Alps), respectively, and eclogite-facies ophicarbonates from the Voltri Massif. The Bracco oceanic ophicarbonates retain breccia-like textures associated with their seafloor hydrothermal and sedimentary origins. Their trace element concentrations and ?18OVSMOW (+15.6 to +18.2‰), ?13CVPDB (+1.1 to +2.5‰) and their 87Sr/86Sr (0.7058 to 0.7068), appear to reflect equilibration during Jurassic seawater-rock interactions. Intense shear deformation characterizes the more deeply subducted ophicarbonates, in which prominent calcite recrystallization and carbonation of serpentinite clasts occurred. The isotopic compositions of the pumpellyite-facies ophicarbonates overlap those of their oceanic equivalents whereas the most deformed blueschist-facies sample shows enrichments in radiogenic Sr (87Sr/86Sr?=?0.7075) and depletion in 13C (with ?13C as low as ?2.0‰). These differing textural and geochemical features for the two suites reflect interaction with fluids in closed and open systems, respectively. The higher-P-metamorphosed ophicarbonates show strong shear textures, with coexisting antigorite and dolomite, carbonate veins crosscutting prograde antigorite foliation and, in some cases, relics of magnesite-nodules enclosed in the foliation. These rocks are characterized by lower ?18O (+10.3 to 13.0‰), enrichment in radiogenic Sr (87Sr/86Sr up to 0.7096) and enrichment in incompatible and fluid-mobile element (FME; e.g., As, Sb, Pb). These data seemingly reflect interaction with externally-derived metamorphic fluids and the infiltrating fluids likely were derived from dehydrating serpentinites with hybrid serpentinite-sediment compositions. The interaction between these two lithologies could have occurred prior to or after dehydration of the serpentinites elsewhere. We suggest that decarbonation and dissolution/precipitation processes operating in ancient subduction zones, and resulting in the mobilization of C, are best traced by a combination of detailed field and petrographic observations, C, O and Sr isotope systematics (i.e., 3D isotopes), and FME inventories. Demonstration of such processes is key to advancing our understanding of the influence of subduction zone metamorphism on the mobilization of C in subducting reservoirs and the efficiency of delivery of this C to depths beneath volcanic arcs and into the deeper mantle.
Abstract: Microcontinents and continental fragments are small pieces of continental crust that are surrounded by oceanic lithosphere. Although classically associated with passive margin formation, here we present several preserved microcontinents and continental fragments associated with subduction systems. They are located in the Coral Sea, South China Sea, central Mediterranean and Scotia Sea regions, and a “proto?microcontinent,” in the Gulf of California. Reviewing the tectonic history of each region and interpreting a variety of geophysical data allows us to identify parameters controlling the formation of microcontinents and continental fragments in subduction settings. All these tectonic blocks experienced long, complex tectonic histories with an important role for developing inherited structures. They tend to form in back?arc locations and separate from their parent continent by oblique or rotational kinematics. The separated continental pieces and associated marginal basins are generally small and their formation is quick (<50 Myr). Microcontinents and continental fragments formed close to large continental masses tend to form faster than those created in systems bordered by large oceanic plates. A common triggering mechanism for their formation is difficult to identify, but seems to be linked with rapid changes of complex subduction dynamics. The young ages of all contemporary pieces found in situ suggest that microcontinents and continental fragments in these settings are short lived. Although presently the amount of in?situ subduction?related microcontinents is meager (an area of 0.56% and 0.28% of global, non?cratonic, continental crustal area and crustal volume, respectively), through time microcontinents contributed to terrane amalgamation and larger continent formation.
Doklady Earth Sciences, Vol. 493, 1, pp. 513-516. pdf
Russia, Yakutia
subduction
Abstract: The data available indicate the complex evolution of deformed peridotites of mantle xenoliths, the P-T parameters of which indicate that they are fragments of the metasomatized lower part of the cratonic lithosphere. The zoning established in garnets from xenoliths in kimberlite pipes is interpreted as a result of metasomatism that occurred shortly before xenoliths reached the surface. Metasomatic alterations in xenoliths of deformed harzburgites were manifested not only in the development of zoning of minerals. The study results show that there is a discrepancy between the data calculated based on the contents of incompatible elements in minerals of xenoliths and those obtained due to direct measurements of the bulk composition of xenoliths. To determine the balance of incompatible elements, we have carried out experiments on leaching xenoliths of deformed lherzolites from the Udachnaya kimberlite pipe. It was established that a significant part of LREEs in the studied xenoliths occurs in the intergranular space. The distribution pattern of incompatible elements and, in particular, the presence of a positive Eu anomaly indicate that the appearance of the intergranular component is not associated with contamination of xenoliths by the kimberlite melt. Quite a few xenoliths demonstrate a positive Eu anomaly, which indicates the influence of the subducted crustal component at one of the modification stages of xenoliths.
Abstract: Understanding water transport by the subducting slab and the corner flow of the mantle wedge is a crucial topic because it is a prime control on seismic tremors, arc-to-intraplate volcanoes as well as on global water distribution in the mantle. However, most of previous studies focused on water transport by the subducting slab and did not quantitatively evaluated the amount of water carried by the corner flow into the deep mantle. Using two-dimensional numerical experiments, we model both the dehydration of the subducting slab and (de)hydration of the mantle wedge and quantify the amount of water transported by both of them. We use the water solubilities of basalt and peridotite derived from laboratory measurements and from thermodynamic calculations, and compare the implications of their differences. Our calculations show that the two models for the water solubilities of basalt result in either abundant or scarce free water through extensive or negligible dehydration of the sub-forearc oceanic crust, leading to a hydrated or a dry cold nose of the mantle wedge, respectively. Further, the oceanic crust of the subducting slab is almost dehydrated prior to reaching a depth of 250 km, regardless of subduction parameters and the models for the water solubilities of basalt. The dehydration depth of the lithospheric mantle of the subducting slab deepens with decreasing slab temperature. The lithospheric mantle of cold subducting slab (e.g., Northeast Japan) experiences partial dehydration at sub-backarc depths and transports the remaining bound water beyond a depth of 250 km, regardless of the models for the water solubilities of peridotite. Deep water transport by the corner flow of the mantle wedge is negligible regardless of the models for the water solubilities of peridotite. The water carried by the lithospheric mantle may be the cause of backarc and intraplate volcanoes in Northeast Asia.
Abstract: Serpentinization greatly affects the physical and chemical properties of lithospheric mantle. Here we address the fate of serpentinized peridotites and their influence over an entire Wilson cycle. We document the near-surface journey of serpentinized subcontinental peridotites exhumed during rifting and continental breakup, reactivated as buoyant material during subduction, and ultimately emplaced as "ophiolite-like" fragments within orogenic belts. This life cycle is particularly well documented in former Tethys margins, where recent studies describe the ongoing incorporation of Mesozoic serpentinized subcontinental peridotites that diapirically rise from a subducting lower plate’s mantle to be emplaced into the accretionary prism in front of a continental arc. This newly recognized mode of subduction-linked serpentine diapirism from the downgoing lithospheric slab is consistent with the origin of some exhumed serpentinized subcontinental peridotites in the Apennines (Italy), these assemblages reaching their present locations during Alpine orogenesis. Transfer of serpentinized subcontinental peridotites from the downgoing to the overriding plate motivates the concept of a potentially "leaky" subduction channel. Weak serpentine bodies may in fact rise into, preferentially migrate within, and eventually leave the intraplate shear zone, leading to strong lateral heterogeneities in its composition and mechanical strength.
Lithos, doi. 10.1016/j. lithos.2020.105949 1p. Pdf
Mantle
subduction
Abstract: Crystal morphologies are essential for deciphering the reaction history of igneous and metamorphic rocks because they often record the interplay between nucleation and growth rates controlled by the departure from equilibrium. Here, we report an exceptional record of the morphological transition of olivine formed during subduction metamorphism and high-pressure dehydration of antigorite-serpentinite to prograde chlorite-harzburgite in the Almirez ultramafic massif (Nevado-Filábride Complex, Betic Cordillera, SE Spain). In this massif, rare varied-textured chlorite-harzburgite (olivine+enstantite+chlorite+oxides) —formed after high-P dehydration of antigorite-serpentinite— exhibits large olivine porphyroblasts made up of rounded cores mantled by coronas of tabular olivine grains, similar to single tabular olivines occurring in the matrix. The correlative X-ray ?-CT and EBSD study of two varied-textured chlorite-harzburgite samples show that tabular olivine in coronas is tabular on (100)Ol with c > b >> a, and grew in nearly the same crystallographic orientation as the rounded olivine cores of the porphyroblast. Quantitative textural analysis and mass balance indicate that varied-textured chlorite-harzburgite is the result of a two-stage nucleation and growth of olivine during the progress of the high-P dehydration of antigorite-serpentinite to chlorite-harzburgite reaction. The first stage occurred under a low affinity (?Gr) and affinity rate (?Gr/dt) of the antigorite dehydration reaction that resulted in a low time-integrated nucleation rate and isotropic growth of olivine, forming rounded olivine porphyroblasts. With further progress of the dehydration reaction, a second stage of relatively higher affinity and affinity rate resulted in a higher time-integrated nucleation rate of olivine coeval with a shift from isotropic to anisotropic olivine growth, leading to tabular olivines. The two-stage evolution resulted in olivine porphyroblasts made up of rounded cores mantled by coronas of tabular olivine grains characteristic of varied-texture chlorite-harzburgite. Although a switch to anisotropic tabular olivine in the second stage is consistent with the relative increase in the affinity and affinity rate, these changes cannot solely account for the growth of Almirez olivine tabular on (100). Tabular olivines in komatiites and other igneous rocks are tabular on (010)Ol with either a > c >> b, or a ? c > > b, in agreement with experimentally determined growth rates of olivine phenocrysts under moderate to high undercooling and cooling rates. On the other hand, olivine tabular on (100) is expected in the presence of highly polymerized fluids where inhibited growth of the olivine (100) and (010) interfaces occurs, respectively, due to dissociative and molecular adsorption of water monolayers. Rounded and tabular olivines in Almirez varied-textured chlorite-harzburgite show differing trace element compositions that we interpret as due to the infiltration of external fluids during antigorite dehydration. Isothermal infiltration of highly polymerized fluids would explain the shift in the affinity and affinity rate of the antigorite dehydration reaction, as well as the olivine morphology tabular on (100) due to the inhibited growth on the (100) and, to a lesser extent, (010). Our study shows that surface-active molecules may play an essential role in shaping the morphology of growing crystals during fluid-present metamorphic crystallization.
Abstract: We propose that subducting slabs may cause lithospheric removal by directing mantle flow along the craton margin. This process could carve and shape the cratons, leading to conditions that impact the overall (in)stability of the lithosphere. We use three-dimensional geodynamic models to investigate how subduction-driven directed flow interacts with cratonic lithosphere of differing shape, concluding that the margin shape controls both channelization of flow around the craton as well as the potential for destruction. While the simulations show that all craton shapes aid in channelization, the cratons with straight vertical margins are the most resistant to deformation, and the cratons with gradually thickening margins are less resistant to deformation. The dependence on shape could contribute to the progressive removal of cratonic lithosphere along its margin in a runaway process until a more stable vertical margin shape evolves.
Abstract: The Trans-Hudson Orogen (THO), formed from the convergence between the Superior craton and the composite Churchill Upper Plate (CUP), is one of the best-preserved examples of a collisional orogen in the Paleoproterozoic. Similar to modern collision systems such as the Himalayan orogen, it is characterized by a composite upper plate in which terrane accretion established a continental plateau that was tectonically and magmatically active for >100 myr. Our study presents new petrological and geochronological data for four samples collected in three lithotectonic domains of the south Rae craton (one of the CUP terranes). The results presented here allow us to re-define the previously proposed extent of THO reworking in the CUP and afford the opportunity to study and compare the evolution of various fragments that illustrate differing levels of a collapsed plateau in the CUP hinterland. The new data indicate that the south Rae craton locally preserves evidence for burial at 1.855-1.84 Ga with peak metamorphic conditions at approximately 790 °C and 9.5-12.5 kbar followed by rapid cooling and decompression melting (P < 6 kbar) at ca. 1.835-1.826 Ga. These results, which provide important and so far missing Pressure-Temperature-time (P-T-t) constraints on the evolution of the south Rae craton in the Northwest Territories at Trans-Hudson time, coupled with existing regional geochronological and geochemical data, are used to propose an updated model for the post-1.9 Ga THO collision and extensional collapse. Our results reveal that: i) initial thickening in the upper plate started at Snowbird time (ca. 1.94 Ga), then continued via Sask collision (with high-grade metamorphism recorded in the south Rae craton, ca. 1.85 Ga), and ended with Superior collision (ca. 1.83 Ga); ii) the extent of the THO structural and metamorphic overprint in the SW CUP is much broader across strike than previously recognized, and iii) T-t data in the south Rae are indicative of relatively fast cooling rates (8-25 °C/Ma) compared to other known Precambrian orogens. We suggest that the Paleoproterozoic THO represents the first record of a major ‘modern-style’ orogenic plateau collapse in Earth’s history.
Abstract: When one tectonic plate dives beneath another at a subduction zone, it recycles huge amounts of water and other chemicals into Earth’s mantle. The sinking plate carries seawater trapped in sediments and crust or chemically bound in minerals like serpentine. Later release of this water in the mantle contributes to key geological processes, such as earthquakes and the formation of volcano-feeding magma.
By volume, the largest portion of a subducting plate is its bottom layer, which comprises upper mantle material. Estimates of the amount of water in downgoing slabs of upper mantle vary widely: Some suggest that worldwide, subduction zones have swallowed more than two oceans’ worth of water in the past 540 million years. However, new research by Miller et al. suggests that water transport at the Middle America Trench subduction zone is an order of magnitude less than previously estimated. As a plate approaches a subduction zone, it bends downward, causing faults to form. Models and earlier observations have suggested that this bending and faulting allow seawater to infiltrate into the upper mantle, where it fills cracks in fault zones, reacts with olivine to produce serpentine, and is later carried deeper into the subduction zone. Previous estimates of how much water reaches the upper mantle along bending faults have relied on measurements of the speed of seismic waves as they pass through a subducting plate. However, those measurements and estimates could not discern whether the upper mantle layer is uniformly hydrated or whether water is confined to bending fault zones. To address that limitation, the new study accounted for seismic anisotropy characterizing how the speed of seismic waves depends on the direction they travel through a material. The researchers used data collected by seafloor seismometers to measure seismic anisotropy along the Middle America Trench near Nicaragua, which enabled a much more detailed picture of upper mantle hydration. The data revealed that in the region studied, water storage in the upper mantle is limited to serpentinized fault zones that thin rapidly with depth, suggesting that fault dynamics and serpentinization reaction kinetics prevent seawater from hydrating the mantle between bending faults. New estimates of water transport that incorporate this finding are an order of magnitude lower than previous estimates for the Middle America Trench. Because the same processes occur at other subduction zones, the researchers report that far less water may be transported worldwide than previously estimated. (Journal of Geophysical Research: Solid Earth, https://doi.org/10.1029/2020JB020982, 2021)
Abstract: Subducting tectonic plates carry water and other surficial components into Earth’s interior. Previous studies suggest that serpentinized peridotite is a key part of deep recycling, but this geochemical pathway has not been directly traced. Here, we report Fe-Ni-rich metallic inclusions in sublithospheric diamonds from a depth of 360 to 750 km with isotopically heavy iron (?56Fe = 0.79 to 0.90‰) and unradiogenic osmium (187Os/188Os = 0.111). These iron values lie outside the range of known mantle compositions or expected reaction products at depth. This signature represents subducted iron from magnetite and/or Fe-Ni alloys precipitated during serpentinization of oceanic peridotite, a lithology known to carry unradiogenic osmium inherited from prior convection and melt depletion. These diamond-hosted inclusions trace serpentinite subduction into the mantle transition zone. We propose that iron-rich phases from serpentinite contribute a labile heavy iron component to the heterogeneous convecting mantle eventually sampled by oceanic basalts.
Annual Review of Earth and Planetary Sciences, Vol. 49, pp. 37-70.
Mantle
subduction
Abstract: Volatile elements (water, carbon, nitrogen, sulfur, halogens, and noble gases) played an essential role in the secular evolution of the solid Earth and emergence of life. Here we provide an overview of Earth's volatile inventories and describe the mechanisms by which volatiles are conveyed between Earth's surface and mantle reservoirs, via subduction and volcanism. Using literature data, we compute volatile concentration and flux estimates for Earth's major volatile reservoirs and provide an internally balanced assessment of modern global volatile recycling. Using a nitrogen isotope box model, we show that recycling of N (and possibly C and S) likely began before 2 Ga and that ingassing fluxes have remained roughly constant since this time. In contrast, our model indicates recycling of H2O(and most likely noble gases) was less efficient in the past. This suggests a decoupling of major volatile species during subduction through time, which we attribute to the evolving thermal regime of subduction zones and the different stabilities of the carrier phases hosting each volatile. This review provides an overview of Earth's volatile inventory and the mechanisms by which volatiles are transferred between Earth reservoirs via subduction. The review frames the current thinking regarding how Earth acquired its original volatile inventory and subsequently evolved through subduction processes and volcanism.
Geochimica et Cosmochimica Acta, Vol. 306, pp. 143-170. pdf
Europe, Italy
subduction
Abstract: The volatile transfer in subduction zones and the role of sulfate as a vector for the mobilization of oxidized components from down-going slabs remain hotly debated issues. Orogenic spinel and garnet peridotite lenses from the Ulten Zone (Eastern Alps, Italy), exhumed as part of felsic metamorphic terranes in continental collision zones, bear witness to mass transfer processes in these pivotal environments. In this study, we carried out a multi-method investigation of mantle sulfides coexisting with four generations of carbonates, indicating coupled sulfur and carbon mobility throughout the peridotites’ metamorphic evolution as part of the Variscan subduction architecture. Detailed petrography, bulk rock measurements, in situ chemical and geochemical analyses of sulfides as well as Sr isotope analyses of associated clinopyroxene and amphibole are combined with the aim to constrain the origin, nature and effect of multiple C-O-H-S-bearing fluids and melts the peridotites interacted with. The first, pre-peak, metasomatic pulse (Stage 1) is represented by an H2S-CO2-bearing melt from the subduction-modified hot mantle wedge, which formed a pyroxenite layer hosting matrix pentlandite with ?34S of +2.77‰. Matrix carbonates occasionally occur in the coarse-grained peridotite under eclogite-facies conditions (Stage 2), with heavier ?34S (up to +3.43‰), radiogenic Sr (87Sr/86Srclinopyroxene > 0.7052) and elevated Pb abundances. These are ascribed to interaction with isotopically heavy melts carrying recycled crustal component, permissive of, but not requiring, involvement of oxidized S species. Conversely, isotopically lighter matrix pentlandite (?34S = ?1.62 to +0.67‰), and radiogenic Sr in amphibole (87Sr/86Sr = 0.7056) and associated dolomite (published data) from fine-grained garnet-amphibole peridotites may point to involvement of H2S-CO2-bearing crustal fluids, which variably equilibrated with the mantle before interacting with the peridotites. The post-peak Stage 3 marks the entrapment of peridotites into a tectonic mélange. Here, kelyphitization of garnet is catalyzed by further ingress of a S-bearing fluid (?34S = ?0.38‰), while carbonate veining with occasional sulfides bear witness to channelized fluid flow. Sulfide and amphibole grains in retrogressed spinel peridotites reveal the highest contents of fluid-mobile elements (As, Sb) and 87Sr/86Sramphibole up to 0.7074, suggesting late interactions with isotopically heavy crustal fluids at high fluid-rock ratios. Textural observations indicate that, during Stage 4, serpentinization of peridotites at low ƒS2 played an active role not only in CO2 release by conversion of dolomite to calcite + brucite intergrowths, but also in local removal of 32S during the final exhumation stage. Late channelized sulfur remobilization is evidenced by the serpentine + magnetite (±millerite ± calcite) vein carrying > 300 ppm S. Overall, the relatively narrow range of sulfur isotope composition (?34S = ?1.62 to +3.76‰) is indicative of limited interaction with isotopically heavy crustal liquids, and points to a subordinate role of subduction-derived sulfate throughout the extended fluid(melt)/rock evolution of the Ulten Zone peridotites, first in the mantle wedge and then as part of a tectonic mélange.
Abstract: The thermal evolution of subducting slabs controls a range of subduction processes, yet we lack a robust understanding of how thermal structure develops over a subduction zone's lifetime. We investigate the time-dependence of slab thermal structure using dynamically consistent, time evolving models. Pressure-temperature (P-T) conditions along the slab Moho and slab top exhibit substantial variability throughout the various phases of subduction: initiation, free sinking, and mature subduction. This variability occurs in response to time-dependent subduction properties (e.g., fast vs. slow convergence) and thermal structure inherited from previous phases (e.g., due to upper plate aging). At a given depth, the slab cools rapidly during initiation, after which slower cooling occurs. In the case of the Moho, additional cooling occurs during the free sinking phase. We explore the implications of time-dependent thermal structure on exhumed rocks and slab dehydration. Modeled slab top P-T paths span much of the P-T space associated with exhumed rocks, suggesting a significant component of recorded variability may have dynamic origins. Coupling our P-T profiles with thermodynamic models of oceanic lithosphere, we show that dehydrating ultramafic rocks at the slab Moho provide the bulk of hydrous fluid at subarc depths during the earliest phases. Over subsequent phases, these rocks carry fluids into the deeper mantle, and it is mafic crust along the slab top that releases water at subarc depths. We conclude that varying subduction conditions, and non-steady-state thermal structure, challenge the utility of kinematically driven models with constant subduction parameters, particularly for investigating thermal structure in the geological past.
Mantle Convection and Surface Expressions, Geophysical Monograph, No. 263, Chapter 7, pp. 179- 14p. Pdf
Mantle
inclusions, subduction, CLIPPER
Abstract: Some rare diamonds originate below the lithosphere, from depths of 300-800 km and perhaps deeper. Ongoing sublithospheric or super-deep diamond research is providing new insight into the mantle and the hidden consequences of plate tectonics. Here we highlight several advances in the past decade, stemming from the discovery of inclusions from oceanic crust at lower mantle depths; inclusions having geochemical imprints of low-degree car-bonatitic melt, possibly from subducted slabs; hydrous ringwoodite and other signs of deep water; major mantle minerals preserved in their original crystal structure, including ringwoodite and CaSiO 3-perovskite; additional diamond varieties with a super-deep origin (CLIPPIR and type IIb diamonds), greatly increasing the known prevalence and diversity of super-deep diamonds; and consistent, recurring Fe-Ni-C-S metallic melt inclusions from depths of 360-750 km. Redox freezing of oxidized, slab-derived fluid/melt upon interaction with ambient metal-saturated mantle appears to be a phenomenon broadly recorded by many super-deep diamonds. Melting of car-bonate, as well as dehydration reactions, from subducted slabs are relevant mechanisms that may generate fluid/ melt contributing to diamond growth. Fe-Ni metal, with dissolved carbon, sulfur, and other elements is also indicated as a possible diamond-forming melt. These mobile and dynamic phases are agents of chemical mass-transfer in the deep mantle.
Abstract: Impingement of a hot buoyant mantle plume head on the lithosphere is one of the few scenarios that can initiate a new subduction zone without requiring any pre-existing weak zones. This mechanism can start subduction and plate tectonics on a stagnant lid and can also operate during active plate tectonics where plume-lithosphere interactions is likely to be affected by plate motion. In this study, we explore the influence of plate motion on lithospheric response to plume head-lithosphere interaction including the effect of magmatic weakening of lithosphere. Using 3d thermo-mechanical models we show that the arrival of a new plume beneath the lithosphere can either (1) break the lithosphere and initiate subduction, (2) penetrate the lithosphere without subduction initiation, or (3) spread asymmetrically below the lithosphere. Outcomes indicate that lithospheric strength and plume buoyancy control plume penetration through the lithosphere whereas the plate speed has a subordinate influence on this process. However, plate motion may affect the geometry and dynamics of plume-lithosphere interaction by promoting asymmetry in the subduction zone shape. When a sufficiently buoyant plume hits a young but subductable moving lithosphere, a single-slab modern-style subduction zone can form instead of multiple subduction zones predicted by stagnant lid models. In the case of subduction initiation of older moving oceanic lithosphere, asymmetrical cylindrical subduction is promoted instead of more symmetrical stagnant lid subduction. We propose that the eastward motion of the Farallon plate in Late Cretaceous time could have played a key role in forming one-sided subduction along the southern and western margin of the Caribbean and NW South America.
Abstract: Subduction of oceanic crust buries an average thickness of 300-500 m of sediment that eventually dehydrates or partially melts. Progressive release of fluid/melt metasomatizes the fore-arc mantle, forming serpentinite at low temperatures and phlogopite-bearing pyroxenite where slab surface reaches 700-900 °C. This is sufficiently high to partially melt subducted sediments before they approach the depths where arc magmas are formed. Here, we present experiments on reactions between melts of subducted sediments and peridotite at 2-6 GPa/750-1100 °C, which correspond to the surface of a subducting slab. The reaction of volatile-bearing partial melts derived from sediments with depleted peridotite leads to separation of elements and a layered arrangement of metasomatic phases, with layers consisting of orthopyroxene, mica-pyroxenite, and clinopyroxenite. The selective incorporation of elements in these metasomatic layers closely resembles chemical patterns found in K-rich magmas. Trace elements were imaged using LA-ICP-TOFMS, which is applied here to investigate the distribution of trace elements within the metasomatic layers. Experiments of different duration enabled estimates of the growth of the metasomatic front, which ranges from 1-5 m/ky. These experiments explain the low contents of high-field strength elements in arc magmas as being due to their loss during melting of sedimentary materials in the fore-arc.
Africa, South Africa, Europe, Greenland, China, Russia, Siberia, Canada, South America, Brazil
subduction, metasomatism
Abstract: To investigate halogen heterogeneity in the subcontinental lithospheric mantle (SCLM), we measured the concentrations of Cl, Br, and I in kimberlites and their mantle xenoliths from South Africa, Greenland, China, Siberia, Canada, and Brazil. The samples can be classified into two groups based on halogen ratios: a high-I/Br group (South Africa, Greenland, Brazil, and Canada) and a low-I/Br group (China and Siberia). The halogen compositions were examined with the indices of crustal contamination using Sr and Nd isotopes and incompatible trace elements. The results indicate that the difference between the two groups was not due to different degrees of crustal contamination but from the contributions of different mantle sources. The low-I/Br group has a similar halogen composition to seawater-influenced materials such as fluids in altered oceanic basalts and eclogites and fluids associated with halite precipitation from seawater. We conclude that the halogens of the high-I/Br group are most likely derived from a SCLM source metasomatized by a fluid derived from subducted serpentinite, whereas those of the low-I/Br group are derived from a SCLM source metasomatized by a fluid derived from seawater-altered oceanic crust. The SCLM beneath Siberia and China could be an important reservoir of subducted, seawater-derived halogens, while such role of SCLM beneath South Africa, Greenland, Canada, and Brazil seems limited.
Earth and Planetary Science Letters, in press available 8p. Pdf
Mantle
subduction
Abstract: Long-term patterns of mantle convection are illustrated by the locations of large low-shear-velocity provinces (LLSVPs) in the lowermost mantle, as well as an enclosing girdle of subduction zones. These structures, stable since Pangea's breakup, have been proposed to provide an absolute reference frame to anchor plate motions in deep time. Simple conceptual models of the supercontinent cycle (introversion, extroversion, orthoversion) predict differing loci of global subduction zones and inferred LLSVPs. We develop a tool to statistically compare idealized supercontinent cycle models with paleolatitude distributions of global subduction zones from paleogeographic reconstructions. We find that subduction zone locations younger than 250 Ma are moderately well described by an idealized girdle around the LLSVPs, but more robust conceptual models must take into account a Tethyan locus or “arm” within the girdle. Between 540 and 250 Ma, such an orthogonal arm is not needed to generate robust correlations with paleogeography; but the global subduction girdle is found to rotate progressively by approximately 90°. Our results suggest that planetary degree-two mantle structures are long-lived but not eternal: they reorganize their absolute locations from one supercontinent cycle to the next.
Earth and planetary Science Letters, Vol. 578, 117293, 8p.
Pangea
subduction
Abstract: Long-term patterns of mantle convection are illustrated by the locations of large low-shear-velocity provinces (LLSVPs) in the lowermost mantle, as well as an enclosing girdle of subduction zones. These structures, stable since Pangea's breakup, have been proposed to provide an absolute reference frame to anchor plate motions in deep time. Simple conceptual models of the supercontinent cycle (introversion, extroversion, orthoversion) predict differing loci of global subduction zones and inferred LLSVPs. We develop a tool to statistically compare idealized supercontinent cycle models with paleolatitude distributions of global subduction zones from paleogeographic reconstructions. We find that subduction zone locations younger than 250 Ma are moderately well described by an idealized girdle around the LLSVPs, but more robust conceptual models must take into account a Tethyan locus or “arm” within the girdle. Between 540 and 250 Ma, such an orthogonal arm is not needed to generate robust correlations with paleogeography; but the global subduction girdle is found to rotate progressively by approximately 90°. Our results suggest that planetary degree-two mantle structures are long-lived but not eternal: they reorganize their absolute locations from one supercontinent cycle to the next.
Science Advances, doi.10.1126/sciadv.abj1325 1p. Pdf
Mantle
subduction
Abstract: Earth’s carbon cycle is strongly influenced by subduction of sedimentary material into the mantle. The composition of the sedimentary subduction flux has changed considerably over Earth’s history, but the impact of these changes on the mantle carbon cycle is unclear. Here, we show that the carbon isotopes of kimberlite magmas record a fundamental change in their deep-mantle source compositions during the Phanerozoic Eon. The 13C/12C of kimberlites before ~250 Ma preserves typical mantle values, whereas younger kimberlites exhibit lower and more variable ratios-a switch coincident with a recognized surge in kimberlite magmatism. We attribute these changes to increased deep subduction of organic carbon with low 13C/12C following the Cambrian Explosion when organic carbon deposition in marine sediments increased significantly. These observations demonstrate that biogeochemical processes at Earth’s surface have a profound influence on the deep mantle, revealing an integral link between the deep and shallow carbon cycles.
Earth and Planetary Science Letters, Vol. 584, 14p. 117475
Mantle
subduction
Abstract: Melting experiments undertaken with finely ground powder of phengite-bearing eclogite yielded solidus temperatures of about 970 °C at 4 GPa and 1250 °C at 9 GPa. Additional experiments with a rock powder of psammopelitic composition established a solidus at 9 GPa at a temperature of 1350 °C. Initial melts produced from both rocks are rich in potassium. The melts generated from eclogite tend to become richer in Na and Ca with rising temperature due to increasing decomposition of clinopyroxene. At the maximum temperatures of the experiments with eclogite, up to 450 °C above the solidus at 4 GPa, this phase is still present in the restite together with abundant garnet. In the temperature interval of 1100-1300 °C, when 22-30% of the studied eclogite was melted, the melts are quartz monzonitic in composition. According to the reported experimental results, we suggest that partial melting of oceanic crust is unlikely in a subduction zone. However, ascending melange diapirs, composed of material from the upper portion of a deep-seated subducted oceanic slab, can partially melt in the hot mantle wedge. The thus generated melts further ascend to contribute to lavas of magmatic arc systems.
Earth and Planetary Science Letters, Vol. 584, 13p. 117424
Mantle
subduction
Abstract: Hydrous arc magmas are produced when water-bearing materials from subducted oceanic plates are transported to the mantle beneath volcanic arcs, though the mechanism of mass transport remains debated. The geochemical characteristics of the slab component have important implications for the thermal structures of down-going plates and the fluxes of elements into the deep mantle. If slab temperatures are low, then elemental fluxes from the slab will be carried in a dilute fluid. If temperatures are high, the slab may melt instead. While a long-standing paradigm for arc volcanism has been that sediments melt and ocean crust dehydrates, a growing body of evidence from arc geochemistry and experimental petrology suggests both sediment and ocean crust melt. The low solubility of many elements in aqueous fluids prevents them from making a substantial contribution to arc mass-balance. Constraints from Sr concentrations and 87Sr/86Sr ratios require a large flux of Sr from the ocean crust, which is only possible if the crust melts. H2O/Sr ratios of arc volcanics are also inconsistent with slab fluids. These conclusions are supported by thermo-mechanical models indicating that slab temperatures exceed the hydrous solidus for both ocean crust and sediments. Examination of experimental data shows a likely strong effect of oxygen fugacity on residual phases during slab melting. Arc data are best explained if the ocean crust melts beneath all arcs under oxidizing conditions somewhere between FMQ and NNO+2. Experimental constraints on sediments also require melting and that sediment melt compositions depend on bulk composition as well as temperature. If these experiments serve as analogs to sediment melting beneath arcs, then sediment bulk compositions are a necessary input for any rare earth element-based slab thermometer. We present compositions for ocean crust partial melts and partition coefficients for sediment melting based on existing experiments, physical models, and arc data, that can be used in geochemical models of arc volcanism.
Abstract: Cubic CaSiO3 perovskite is a major phase in subducted oceanic crust, where it forms at a depth of about 550?kilometres from majoritic garnet1,2,28. However, its rheological properties at temperatures and pressures typical of the lower mantle are poorly known. Here we measured the plastic strength of cubic CaSiO3 perovskite at pressure and temperature conditions typical for a subducting slab up to a depth of about 1,200?kilometres. In contrast to tetragonal CaSiO3, previously investigated at room temperature3,4, we find that cubic CaSiO3 perovskite is a comparably weak phase at the temperatures of the lower mantle. We find that its strength and viscosity are substantially lower than that of bridgmanite and ferropericlase, possibly making cubic CaSiO3 perovskite the weakest lower-mantle phase. Our findings suggest that cubic CaSiO3 perovskite governs the dynamics of subducting slabs. Weak CaSiO3 perovskite further provides a mechanism to separate subducted oceanic crust from the underlying mantle. Depending on the depth of the separation, basaltic crust could accumulate at the boundary between the upper and lower mantle, where cubic CaSiO3 perovskite may contribute to the seismically observed regions of low shear-wave velocities in the uppermost lower mantle5,6, or sink to the core-mantle boundary and explain the seismic anomalies associated with large low-shear-velocity provinces beneath Africa and the Pacific.