Mantle Temperature Research Articles

Magnetotelluric image of the Patagonian slab window: Constraints on upper mantle physical properties and sources of intraplate magmatism

The Patagonian slab window (PSW) is a region of the Southamerican subduction zone where the absence of subducted slabs is interpreted, due to the subduction of the Chile Mid-Ocean Ridge at the Chile Triple Junction. Here we report the results of a long-period magnetotelluric (MT) study conducted in two 300 km-long trench-parallel transects crossing the northern boundary of the PSW in the proximal backarc. We modeled the MT data using 3-D inversion, obtaining an electrical resistivity model of the continental crust and upper mantle up to a depth of ∼150 km. Our model shows a heterogeneous resistivity structure in the uppermost mantle, dominated by resistivities >300 Ωm below the array of sites even within the PSW, and some low-resistivity zones (LRZs, <10 Ωm) mainly at the edge of the array. Using petrophysical models, we estimated the mantle temperature, water content, melt fraction, and viscosity based on obtained resistivity values and a preexistent model of P-wave velocity (Vp) at 50 km and 100 km depth. These estimates suggest that the uppermost mantle within the PSW region is heterogeneous and dominated by high-viscosity blocks, compatible with the continental mantle lithosphere or even subducted slabs. Based on relatively hot and low-viscosity zones estimated in the periphery of LRZs, we interpret the presence of asthenospheric mantle in areas where LRZs coincide with relatively low Vp. According to this interpretation, asthenospheric upwelling in the study area at depths ≤150 km would be localized rather than ubiquitous over the interpreted extent for the PSW. Such localized asthenosphere upwelling processes in the past could explain the scattered distribution of Neogene basaltic lavas in the southern Patagonia backarc. The continental crust exhibits LRZs in the upper and lower crust. Remarkably, ensembles of LRZs at different crustal depths within the presumable area of the PSW were found below the General Carrera Lake, and towards the North Patagonian Icefield, likely indicating the presence of hidden intraplate magmatic and/or hydrothermal systems.

Estimating Primary Magmas From Mars With PRIMARSMELT: Implications for the Petrogenesis of Some Martian Rocks and the Thermal Evolution of Mars

AbstractPrimary magmas form by partial melting in the mantle of a terrestrial planet and represent the starting material for building its crust. The compositions of primary magmas are critical for understanding the thermal history of planetary interiors, as they can be used to estimate mantle potential temperatures (TP) and track changes in the conditions of mantle partial melting over time. Here, we introduce PRIMARSMELT, a new member of the PRIMELT software family, calibrated to estimate the composition of Martian primary magmas and their formation conditions. We applied PRIMARSMELT to a comprehensive database of basaltic compositions from Mars. Our results are consistent with their petrology, requiring olivine addition to restore fractionated compositions to their primary parents and olivine subtraction from cumulate rocks. Individual primary magma solutions provide insights into the petrogenesis of specific Martian meteorites, with implications for the near‐primary nature of some primitive meteorites and the relationship between lithologies A and B in meteorite EETA 79001. Taken together, our results suggest nearly constant or potentially increasing mantle potential temperatures throughout the geological history of Mars. The average TP for young shergottite meteorites is ∼1,442 ± 40°C, similar to ambient mantle temperatures inferred from geophysical models. In contrast, older basaltic rocks record potential temperatures as low as ∼1,320 ± 48°C for igneous clasts in meteorites NWA 7034/7533. We suggest that, rather than plume‐related magmatism, shergottite meteorites record ambient mantle temperatures, with the thermal evolution trend possibly resulting from inefficient heat loss, as expected for a planet in stagnant‐lid mode.

New constraints to subduction zone arc and backarc mantle temperatures: a test of the corner flow model

The heat source for the hot volcanic arc of subduction zones and the uniformly warm backarc is not well understood, particularly where backarc extension is absent. In the widely studied corner flow model, the traction of the underthrusting plate drags down the overlying viscous mantle, and the induced shallower seaward flow brings in heat from the landward and deeper asthenosphere. However, earlier comparisons with observations, especially backarc heat flow, gave poor agreement. There are several new constraints to the temperatures and structures above the underthrusting plate and in the backarc, especially those in western North America, for comparison with model predictions. Seismic receiver functions map horizontal boundaries, including the lithosphere-asthenosphere boundary (LAB). Seismic shear wave tomography gives mantle velocity-depth that provides an approximate proxy for temperature. Attenuation of seismic wave speed indicates high temperature or partial melting. The presence of recent backarc volcanics and their geochemistry give the depth and temperature of the LAB. From these constraints in western North America, the LAB reflects ponding of partial melt, is consistently at a depth of 65±3 km and temperature of 1350±25°C, and overlies nearly constant temperature (adiabatic) asthenosphere across the backarc from Mexico to Alaska. These depth and temperature estimates are greatly different from corner flow model predictions. An alternative hypothesis more consistent with these constraints is vigorous small-scale convection in the backarc asthenosphere. Whether or how this convection coexists with the large-scale flow patterns inferred from the interpretation of seismic anisotropy analyses, and to what degree it prevails in the arc-forearc region, deserve further research.

A shallow mantle seismic discontinuity beneath northeast India: Evidence from receiver function analyses

The crust and shallow upper mantle structure beneath the Upper Brahmaputra Valley, Indo-Burma Ranges, and Bengal Basin of Northeast India have been investigated based on receiver function (RF) analysis of teleseismic earthquakes recorded by 11 seismological stations. The study reveals a thin crust (∼35 km) beneath the Brahmaputra Valley (at JORH station) with a surface sedimentary layer of ∼4 km thick. The crustal thickness is observed to increase towards the north in the Himalaya (∼40 km at ZIRO and ITAN) and to the south (up to ∼46 km at KOHI). The crustal thickness near the Tripura fold-belt and Bengal Basin varies within ∼36–40 km. The study reveals the existence of a shallow mantle discontinuity (Hales discontinuity) at a variable depth range of ∼54–78 km characterized by a step increase (∼7.5–11 %) in shear wave velocity observed in the inverted models. The mineralogical phase transformation from spinel to garnet is considered as the origin of this discontinuity. The shallow depth of the discontinuity indicates an increase in upper mantle temperature which conforms to the high geothermal gradient reported in the region. The variation of depth of the discontinuity can be interpreted in terms of the addition of Cr+3 that shifts the spinel-garnet stability field to higher depths whereas Fe+2 shifts it to lower depths. Despite the high temperature in the upper mantle, the observed low Vp/Vs ratio (1.65–1.75) below the Hales discontinuity can be explained by the presence of a high fraction of orthopyroxene.

High pressure and high temperature Brillouin scattering measurements of pyrope single crystals using flexible CO2 laser heating systems

Single-crystal Brillouin scattering measurements are important for interpreting seismic velocities within the Earth and other planetary interiors. These measurements are rare, however, at temperatures above 1000 K, due to the fact that the transparent samples cannot be heated by common laser heating systems operating at a wavelength on the order of 1 μm. Here we present Brillouin scattering data on pyrope collected at pressures up to 23.8 GPa and temperatures between 850 and 1900 K using a novel CO2 laser heating system confined in either a flexible hollow silica waveguide or an articulated arm with mirrors mounted in each junction to direct the laser to the exit point. Pyrope has been chosen because it has been extensively studied at high pressures and moderate temperatures and therefore it is an excellent sample for bench-marking the CO2 laser heating system. The new high-temperature velocity data collected in this study allow the room pressure thermal parameters of pyrope to be constrained more tightly, resulting in values that reproduce the temperature dependence of the unit-cell volume of pyrope measured in recent studies at ambient pressure. Aggregate wave velocities of pyrope calculated along an adiabat using the thermoelastic parameters determined in this study are larger than those obtained using published values, implying that velocities for many mantle components may be underestimated at mantle temperatures because high temperature experimental data are lacking.

Seismic evidence of an extremely thin crust along a 70 Ma slow-spreading crustal segment in the equatorial Atlantic Ocean

The oceanic crust at slow-spreading ridges is formed by a combination of magmatic and tectonic processes. Seismic tomographic studies, commonly used to determine crustal structures, indicate that mature slow-spreading crust has an average thickness of 6.1 ± 0.9 km, with a lower to upper crustal thickness ratio of 2.3. However, tomographic methods are based on high-frequency approximations and can only provide information about large-scale crustal structures. Using seismic full-waveform inversion applied to ocean bottom seismometer data, we report the presence of nearly uniform and extremely thin crust (4 ± 0.2 km) with a lower to upper crustal thickness ratio of 0.74 along an ∼100 km, 70 Ma crustal segment formed at the slow-spreading Mid-Atlantic Ridge in the equatorial Atlantic Ocean. This thin crust is underlain by a Moho transition zone that is 1 ± 0.5 km thick, resulting in a combined crustal and Moho transition zone thickness of 5 ± 0.6 km. The crustal thinning is primarily due to the thinning of the lower crust, suggesting that a significant part of the crust is formed by lava flows and dike intrusions. The presence of an extremely thin crust indicates a decrease in mantle temperature by at least 50 °C, suggesting a vast extent of the upper mantle thermal minimum in the equatorial Atlantic Ocean, and the absence of influence of the mantle from the Siera Leone mantle plume 70 m.y. ago.

Tracking subduction-related metasomatism of the subcontinental lithospheric mantle using Ca-, O-, and H-isotopes

Mantle xenoliths provide effective records of the metasomatic processes that affect continental lithosphere evolution, such as interaction with subducted components or modification via small-degree melts. Correlations between major/trace element geochemistry with stable and radiogenic isotope compositions can help constrain the source and timing of this metasomatism. We report new δ18O, δ44/40Ca, and δD values for twelve kimberlite-hosted mantle xenoliths from the Slave Craton (NWT, Canada), which show varying degrees of metasomatism. The δ18O values of olivine (δ18Ool = +5.33 ± 0.13‰; 1σ; n = 12) overlap average mantle values. Clinopyroxene and garnet δ18O values (δ18Ocpx = +5.31 ± 0.10‰; δ18Ogrt = +5.37 ± 0.23‰; 1σ) extend below those reported in most mantle peridotites and are strongly correlated with clinopyroxene δ44/40Ca (avg. = +1.00 ± 0.10‰; 1σ) and garnet δ44/40Ca (avg. = +1.18 ± 0.19‰; 1σ) respectively, extending from typical mantle values to low δ18O and high δ44/40Ca values. In general, Δ18Ocpx-ol and Δ18Ogrt-ol (ranging from −0.19‰ to +0.19‰ and from −0.56‰ to +0.35‰, respectively) are lower than expected equilibrium values at mantle temperatures. Strong negative correlations are found between δ18Ogrt and Δ18Ogrt-ol and garnet major and trace element composition (Na2O, H2O, La/YbN). Furthermore, phlogopite-bearing kelyphitic rims have δD values (avg. = −126 ± 13‰; 1σ) lower than typical mantle values. Whole rock Sm-Nd model ages and oxygen isotope diffusion modeling suggest that metasomatism occurred during the Mesozoic, shortly before kimberlite entrainment, consistent with indications from diamond-forming fluids from the Slave craton. The combined low δ18O, δD, and high δ44/40Ca signature of the mantle peridotite xenoliths, along with the age constraints, suggest the metasomatic fluid/melt is sourced from a recycled oceanic crust component related to Mesozoic subduction in western North America.

Plate tectonics through Earth’s history: constraints from the thermal evolution of Earth’s upper mantle

ABSTRACT Temporal changes in Earth’s tectonic style play a crucial role in understanding the planet’s evolution. Modern-style tectonics is characterized by the formation of basaltic crust at divergent plate boundaries and its subsequent recycling at subduction zones, accompanied by wedge mantle formation and arc magmatism. It is commonly believed that the secular cooling of the mantle modified the tectonic style from stagnant lid or heat pipe on an early hotter Earth to horizontal tectonics during Meso-Neoarchean. However, various field, petrographic, and geochemical studies suggest that the onset of plate tectonics ranges from Hadean to Neoproterozoic. In this study, we re-evaluate the primary magma temperature (mantle potential temperature, Tp) of the upper ambient mantle, spanning from Eoarchean to Neoproterozoic. We used basalts from several Archean and Proterozoic greenstone belts worldwide, along with Proterozoic ophiolites, to (re)calculate the Tp using the FRACTIONATE-PT method. We observed a Tp range of 1444–1639°C during the Archean and 1414–1611°C during the Proterozoic. These findings indicate a strong correlation with previously estimated Tp values obtained from PRIMELT3 method. This further highlights strong internal consistency among different methods and supports models of a hot ambient mantle during the Archean and Proterozoic. We further reviewed numerical models regarding the effect of mantle temperature on the viability of early Earth plate tectonics. Such model results, composition of Archean continental crust, recent crustal growth models, and field evidence are consistent with the operation of plate tectonics on an early hotter Earth. The transition from a hotter mantle to a colder one from Eo-Neoarchean resulted in thinner oceanic lithosphere and a less depleted lithospheric mantle. Subduction of this thinner oceanic lithosphere led to modern-style tectonics. The intensity and style of plate tectonics on the early Earth differed from modern tectonics, which emerged during the Neoarchean.

Highly variable magmatic accretion at the ultraslow-spreading Gakkel Ridge

Crustal accretion at mid-ocean ridges governs the creation and evolution of the oceanic lithosphere. Generally accepted models1–4 of passive mantle upwelling and melting predict notably decreased crustal thickness at a spreading rate of less than 20 mm year−1. We conducted the first, to our knowledge, high-resolution ocean-bottom seismometer (OBS) experiment at the Gakkel Ridge in the Arctic Ocean and imaged the crustal structure of the slowest-spreading ridge on the Earth. Unexpectedly, we find that crustal thickness ranges between 3.3 km and 8.9 km along the ridge axis and it increased from about 4.5 km to about 7.5 km over the past 5 Myr in an across-axis profile. The highly variable crustal thickness and relatively large average value does not align with the prediction of passive mantle upwelling models. Instead, it can be explained by a model of buoyant active mantle flow driven by thermal and compositional density changes owing to melt extraction. The influence of active versus passive upwelling is predicted to increase with decreasing spreading rate. The process of active mantle upwelling is anticipated to be primarily influenced by mantle temperature and composition. This implies that the observed variability in crustal accretion, which includes notably varied crustal thickness, is probably an inherent characteristic of ultraslow-spreading ridges.

Double subduction controls on long-lived continental tectonics and subcontinental mantle temperatures

The complexities of convergent margins commonly include the interactions of subduction zones, with many geological records of “double” subduction. Here, we build two-dimensional numerical models to explore the evolution of complex subduction systems by systematically testing single and inward-dipping double subduction beneath a continental upper plate and the impact of continental collision on these systems. When compared to single subduction models, the inward-dipping double subduction shows hindered trench migrations and larger volumes of upwelling mantle enhanced by excess sinking slab mass. Double subduction draws larger volumes of hotter mantle beneath the continent in an area much broader than the marginal basins of single subductions, contributing to subcontinental heating by ∼200 °C. As collision jams one margin of a double subduction system, the other margin follows the evolution of migrating single subduction zones, although characterized by persisting higher mantle temperatures and strong upwellings, inherited from the double subduction stage, and large-scale upper plate extension. The modeling outcomes are compared to scaling arguments to test the viability of the mechanism proposed for tectonics of the Cenozoic South China Sea and Neoproterozoic Yangtze Block (southeastern China), where the inward-dipping double subduction provides a context for protracted large-scale continental extension, hotter subcontinental temperatures, and channeled mantle flow not easily reconciled with the dynamics of single subduction zones.

Coupled fates of Earth's mantle and core: Early sluggish-lid tectonics and a long-lived geodynamo.

Conventional Earth evolution models are unable to simultaneously reproduce two fundamental observations: the mantle's secular temperature record and a long-lived geodynamo before inner core nucleation. Today, plate tectonics efficiently cools the mantle, but if assumed to operate throughout Earth's history, past mantle temperature and plate motion become unrealistically high. Through coupled core-mantle modeling that self-consistently predicts multiple mantle convection regimes, we show that over most of the Precambrian, Earth likely operated in a distinct "sluggish-lid" tectonic mode, characterized by partial decoupling between the lithosphere and mantle. This dominant early regime is due to a hotter Earth and the presence of the asthenosphere. This mode regulates the core-mantle boundary heat flow, which powers the geodynamo before inner core nucleation. Both sluggish-lid tectonics and a long-lived dynamo demonstrate the inextricably connected paths of the core-mantle system. Moreover, our simulations simultaneously satisfy diverse geological observations and are consistent with emerging interpretations of such records.

Stabilizing inverse ringwoodite with defects, and a possible origin for the 560-km seismic discontinuity

Ringwoodite is an important mineral in the mantle transition zone, and its cationic disorder can profoundly affect its physicochemical properties, but there is currently much controversy about this disorder. In this study, we investigate the cation disorder states of pure Mg2SiO4-ringwoodite and defective ringwoodite under mantle transition zone conditions through DFT calculations and thermodynamic models. Two stable endmembers are seen, one with normal ringwoodite structure and the other with inverted structure (its Si atoms and half of its Mg atoms have swapped sites). Our results indicate that pure ringwoodite does not invert (swap Mg and Si cations) under normal mantle temperatures but the introduction of a Si-excess, Mg-deficient defect induces a swap at normal mantle temperatures and this swap is likely induced by a wide range of defects including water. Thus, in the presence of such a defect or similar defects the olivine phase transition sequence may then go from olivine to wadsleyite to inverse ringwoodite, and then normal ringwoodite. We calculate the seismic properties of normal and inverse ringwoodite and find significantly slower wave speeds in inverted ringwoodite. Due to this difference the presence of inverse ringwoodite may provide a potential explanation for the discontinuous interface of seismic waves at the depth of ∼560 km.

Factors Contributing to Slab Locations and Geometries in Reconstructions of Past Mantle Flow

AbstractIndividual sinking slabs present markedly different geometries between 410 and 660 km depths, from vertical slabs penetrating the lower mantle to slabs stagnating above the lower mantle. The proposed factors determining these contrasted geometries include mantle viscosity and the magnitude and evolution of trench retreat. Here, we assess the success of paleo‐geographically driven global mantle flow models in matching slabs in tomographic models between 400 km and 1,000 km depth. We quantify the spatial match between predicted present‐day mantle temperature anomalies and vote maps of tomographic models. We investigate the sensitivity of the spatial match to input parameters of the mantle flow model: imposed tectonic reconstruction, model start age, and viscosity contrast between the upper and lower mantle. We evaluate the visual match between model slabs and tomographic vote maps for three circum‐Pacific regions with contrasted slab dip angles between 400 km and 1,000 km depth. Predicted model slabs better match slabs inferred from tomography when there is an increase in viscosity at 660 km depth. The temporal evolution of the models and the global match at present day suggest that the subduction history could be refined in the global tectonic reconstructions that we considered. For example, we suggest that the subduction to the east of Japan should be offset by approximately 100 km to the west at ∼80 Ma to match the anchoring of a continuous slab into the lower mantle suggested by tomography.

The distribution and generation of carbonatites

The physio-chemical framework that generates carbonatites and, ultimately, the associated rare earth element deposits remains contentious. This primarily reflects the diverse tectonic settings in which carbonatites occur: large igneous provinces, continental rifts and major extensional terranes, syn- to post-collisional settings, or ocean islands. There is, however, a broad consensus that carbonatites (or their parental melts) originate in the mantle. These exotic melts have small volumes that make them ideal probes of conditions in their underlying source regions. We combine the carbonatite locations with global maps of lithospheric thickness, derived from seismic tomography, and show that post-Neoproterozoic carbonatites occur preferentially above the margins of thick cratonic lithosphere (e.g., adjacent to the South Atlantic and Indian Oceans or in North America, Greenland, and Asia) and where once thick lithosphere has undergone stretching (e.g., eastern Asia). Our thermal modeling reveals that lateral and vertical heat conduction on rifted craton margins, or rapid stretching of cratonic lithosphere, can mobilize carbonated peridotite at the temperatures (950−1250 °C) and pressures (2−3 GPa) required to form primary carbonatites or their parental alkali silicate melts. Importantly, our models show that heat conduction from upwelling mantle plumes or ambient mantle on rifted cratonic margins may sufficiently modify the temperature of the lithospheric mantle to cause melting of carbonated peridotite, settling the long-standing debate on the role of rifting and heating in the generation of carbonatites.

Effect of Al‐Incorporation on the Sound Velocities of Superhydrous Phase B at High Pressure and High Temperature

AbstractSound velocities and densities of Al‐bearing superhydrous phase B (Al‐bearing SuB) were investigated up to 24 GPa and 1300 K by synchrotron X‐rays combined with ultrasonic interferometry techniques in a multi‐anvil apparatus. We found that Al + H incorporation decreases the adiabatic bulk modulus and shear modulus of SuB. Our results however suggest that this effect is less important than that of temperature. The presence of hydrous pyrolite with ∼10 wt.% Al‐bearing SuB in cold subducting slabs could explain up to ∼2% high velocity anomalies at the bottom of the mantle transition zone (500–660 km) while it turns into ∼7.7% low velocity anomalies below 660 km. Al‐bearing SuB is reportedly stable at mantle temperatures, where it could account for ∼12.2% low velocity anomalies beneath subduction zones in the uppermost lower mantle.

Determining mid-ocean ridge geography from upper mantle temperature

In this study, we examine the influence of the mantle and large-scale tectonics on the global mid-ocean ridge (MOR) system. Using solely seismically-inferred upper mantle temperatures below the melting zone (260-600 km) and an interpretable machine learning model (Random Forest and Principal Component Analysis), we can predict a priori the location (ocean basin and ridge system) of global MOR with up to 90% accuracy. Two features provide > 50% of the discriminative power: the temperature difference between the mid-layer (340-500 km) and other depths, and the depth-averaged temperature of the upper mantle. We suggest long-term (100s Myr) tectonic and large-scale convective processes primarily driven by deep subduction left ample, distinct but hidden fingerprints in the mantle that allow us to separate regions at ∼1000 km scale. Our result implies that the large-scale geophysical and geochemical differences observed along the MOR system are reflective, not primarily of shallow processes associated with melting, but of the integrated long-term tectonic, subduction, and convective flow history that determines the present-day upper mantle temperature structure.

Generation of Archean TTGs via sluggish subduction

Abstract The trondhjemite-tonalite-granodiorite (TTG) suite of rocks prominent in Earth’s Archean continents is thought to form by melting of hydrated basalt, but the specific tectonic settings of formation are unclear. Models for TTG genesis range from melting of downgoing mafic crust during subduction into a hotter mantle to melting at the base of a thick crustal plateau; while neither uniquely defines a global tectonic regime, the former is consistent with mobile lid tectonics and the latter a stagnant lid. One major problem for a subduction model is slabs sinking too quickly and steeply in a hotter mantle to melt downgoing crust. I show, however, that grain size reduction in the lithosphere leads to relatively strong plate boundaries on the early Earth, which slow slab sinking. During this “sluggish subduction,” sinking plates can heat up enough to melt when the mantle temperature is ≳1600 °C. Crustal melting via sluggish subduction can thus explain TTG formation during the Archean due to elevated mantle temperatures and the paucity of TTG production since due to mantle cooling.

Sound velocities in lunar mantle aggregates at simultaneous high pressures and temperatures: Implications for the presence of garnet in the deep lunar interior

Recent experimental and theoretical studies on lunar magma ocean crystallisation have suggested the presence of significant proportions of garnet in the deep lunar interior. While phase relation studies indicate a deep lunar mantle consisting of olivine, pyroxene, and garnet, the compatibility of such an assemblage with seismic models of the lunar interior is yet untested. In this study we report compressional and shear wave velocities in an iron-rich assemblage consisting of olivine, orthopyroxene, clinopyroxene, and garnet up to ∼8 GPa and 1300 K, by means of ultrasonic interferometry measurements combined with synchrotron techniques using the multi-anvil press apparatus. Sound velocity and density models of lunar mantle rocks along a selenotherm based on our experimental results find good agreement with the seismic and density profiles at lunar interior depths of 740–1260 km. Further models are constructed, allowing for the variation of chemical composition, phase proportion, and temperature; these suggest that a garnet-rich deep lunar mantle is compatible with present-day lower lunar mantle temperatures of between 1400–1800 K. Our results show that lunar mantle rocks with up to 33 wt.% garnet may provide an explanation for the observed high velocities of the lower lunar mantle. The presence of garnet in the lowermost part of the Moon's mantle has significant implications for the depth and temperature of the Moon's magma ocean as well as the composition, structure and internal dynamics of the solid Moon.

Alumina solubility in periclase determined to lower mantle conditions and implications for ferropericlase inclusions in diamonds

A series of high-pressure multi-anvil experiments were conducted in the MgO-Al2O3 system between 15 and 50 GPa and at temperatures up to 2623 K, to determine the solubility of Al2O3 in periclase in the stability fields of coexisting spinel, corundum, the Mg2Al2O5 modified ludwigite phase, the MgAl2O4 calcium ferrite phase, and the MgAl2O4 calcium titanate phase. The Al2O3 solubility in periclase is strongly temperature-dependent over the conditions investigated. Conversely, periclase Al2O3 solubility exhibits a negative relationship with pressure. Experiments with up to 40 mol.% FeO in ferropericlase show Al2O3 solubilities that are near identical to those of periclase, within experimental uncertainties. A simple thermodynamic model incorporating significant literature on ambient pressure measurements is able to describe periclase Al2O3 solubility up to 40 GPa. This model is used to investigate the Al2O3 contents of ferropericlase inclusions observed in natural diamonds, which vary up to 0.35 mol. %. Pressure–temperature curves along which particular inclusions formed can be produced if equilibrium can be assumed with Al-bearing minerals found in the same diamond. Alternatively, the maximum possible Al2O3 content in ferropericlase can be determined for a certain pressure–temperature profile and inclusions with Al2O3 contents that exceed this curve can be excluded from those conditions. To obtain such solubility curves, calculations are performed for periclase coexisting with the Al-rich phases spinel, garnet, bridgmanite and the MgAl2O4 calcium ferrite phase. The calculations indicate that periclase, and by inference ferropericlase, Al2O3-contents cannot be greater than 0.5 mol.% under present day adiabatic mantle temperatures and go through a minimum at mantle transition zone conditions. This excludes a number of Al-rich ferropericlase inclusions found in natural diamonds from being formed in the transition zone, unless temperatures were super-adiabatic. This subset of inclusions likely formed either at the base of the upper mantle or the top of the lower mantle, but must have formed at near adiabatic temperatures. The majority of ferropericlase inclusions have Al2O3 contents that would be consistent with formation in the transition zone at near slab temperatures, but could still have been formed at higher temperatures if Al2O3 activities were low.

Onset of Plate Motion in the Presence of Chemical Heterogeneities in the Mantle and the Effect of Mantle Temperature

AbstractChemical heterogeneities of various origins have been observed in the Earth's mantle. Their assumed higher density acts on the style of mantle convection and therefore, as tectonic plates form the highly viscous upper boundary layer of mantle convection, the onset of surface motion is also affected by chemical heterogeneities. We perform 2D basally heated thermochemical mantle convection models which initially include layers of dense material at different mantle depths. In comparison to purely thermal convection models, the onset of first plate motion is delayed. This delay is even more pronounced, the deeper the location of dense material, the larger its volume and the higher its density. Additionally, we varied the starting temperature of our models. In an initially hot mantle, a strong temperature‐dependent viscosity contrast between the mantle and cold surface leads to a cold, highly viscous plate (stagnant lid). For an initially cold mantle, rising plumes first need to transport basal heat upwards to create a sufficient viscosity contrast. Consequently, the formation and subsequent breaking of the stagnant lid is delayed. Considering a hotter mantle after Earth's magma ocean phase, we find that plate motion can occur within approximately the first 0.5 Gyr of solid‐state convection if no chemical structures are present or for dense material situated at the surface. Deep chemical heterogeneities delay the onset by at least one order of magnitude. Furthermore, the early surface motion will have been more erratic than todays stable plate motion, probably with a few single subduction events.