Asthenospheric Temperature Research Articles

The Effect of Intracrystalline Water on the Mechanical Properties of Olivine at Room Temperature

AbstractThe effect of small concentrations of intracrystalline water on the strength of olivine is significant at asthenospheric temperatures but is poorly constrained at lower temperatures applicable to the shallow lithosphere. We examined the effect of water on the yield stress of olivine during low‐temperature plasticity using room‐temperature Berkovich nanoindentation. The presence of water in olivine (1,600 ppm H/Si) does not affect hardness or yield stress relative to dry olivine (≤40 ppm H/Si) outside of uncertainty but may slightly reduce Young’s modulus. Differences between water‐bearing and dry crystals in similar orientations were minor compared to differences between dry crystals in different orientations. These observations suggest water content does not affect the strength of olivine at low homologous temperatures. Thus, intracrystalline water does not play a role in olivine deformation at these temperatures, implying that water does not lead to weakening in the coldest portions of the mantle.

Carboniferous back-arc bimodal rocks in West Kunlun during northward subduction of the Paleo-Tethys

ABSTRACT The widespread occurrence of Carboniferous Wuluate Formation volcanic rocks in the West Kunlun is closely related to the evolution of the Paleo-Tethys but has received little attention. In this study, we conducted detailed fieldwork on the Wuluate Formation in West Kunlun and performed zircon U‒Pb dating, whole-rock major and trace element analysis, whole-rock Sr-Nd isotope analysis, and in situ zircon Hf isotope analyses. The Wuluate Formation is characterized by a bimodal volcanic suite. Zircon U‒Pb dating indicates that the rhyolites erupted at ca. 348–353 Ma, which is coeval with the zircon U‒Pb dating result of 354 Ma for the basalt. The basalts show a tholeiitic magmatic evolution trend with flattened REE patterns and depleted whole-rock Sr-Nd and zircon Hf isotopes, indicating geochemical affinity to back-arc basin basalt (BABB). Trace element modelling results indicate that the primitive basaltic compositions fit well with the products of the DMM + sediments + H2O source under melting conditions of asthenospheric temperature and thin lithosphere thickness, and the REE modelling results show that the basalts underwent 40% to 95% crystallographic differentiation of olivine, monazite, hornblende, plagioclase, and magnetite. The silicic rocks are calcic to calcic-alkaline and meta/peraluminous. They exhibit flattened REE patterns and slightly enriched zircon Hf isotopes, which could originate from Precambrian metaigneous sources of West Kunlun. The Kungai bimodal suite is likely to have formed in a back-arc setting corresponding to the northward subduction of Paleo-Tethys. The subduction of Paleo-Tethys had already occurred, and a back-arc basin had developed before the Early Carboniferous.

Lateral variation in slab window viscosity inferred from global navigation satellite system (GNSS)–observed uplift due to recent mass loss at Patagonia ice fields

AbstractThe geographic coincidence of the Chile Ridge slab window and the Patagonia ice fields offers a unique opportunity for assessing the effects of slab window rheology on glacial isostatic adjustment (GIA). Mass loss of these ice fields since the Little Ice Age causes rapid but variable crustal uplift, 12–24 mm/yr around the North Patagonia ice field, increasing to a maximum of 41 mm/yr around the South Patagonia ice field, as determined from newly collected or processed geodetic data. We used these observational constraints in a three-dimensional Maxwell viscoelastic finite element model of GIA response above both the subducting slab and slab window in which the upper-mantle viscosity was parameterized to be uniform with depth. We found that the viscosity of the northern part of the slab window, ~2 × 1018 Pa·s, is lower than that of the southern part by approximately an order of magnitude. We propose that this along-strike viscosity contrast is due to late Cenozoic ridge subduction beneath the northern part of the slab window, which increases asthenospheric temperature and reduces viscosity.

Global variation of seismic energy release with oceanic lithosphere age

Variations in Mid Ocean Ridge seismicity with age provide a new tool to understand the thermal evolution of the oceanic lithosphere. The sum of seismic energy released by earthquakes during a time, and for an area, is proportional to its lithospheric age. Asthenospheric temperatures emerge on ridge centers with new crust resulting in high seismic activity; thus, the energy released sum is highest on the young lithosphere and decreases with age. We propose a general model that relates the systematic variation of seismic energy released with the lithospheric age. Our analysis evaluates the main physical factors involved in the changes of energy released sum with the oceanic lithosphere age in MOR systems of different spreading rates. These observations are substantiated based on three cross-sections of the East Pacific Rise, six sections in the Mid Atlantic Ridge, and three profiles in the Central Indian Ridge. Our global model provides an additional tool for understanding tectonic processes, including the effects of seismicity and mid-plate volcanism, and a better understanding of the thermal evolution for the young oceanic lithosphere.

Quantifying the Relationship Between Short‐Wavelength Dynamic Topography and Thermomechanical Structure of the Upper Mantle Using Calibrated Parameterization of Anelasticity

AbstractOceanic residual depth varies on 5,000 km wavelengths with amplitudes of ±1 km. A component of this short‐wavelength signal is dynamic topography caused by convective flow in the upper ∼300 km of the mantle. It exerts a significant influence on landscape evolution and sea level change, but its contribution is often excluded in geodynamic models of whole‐mantle flow. Using seismic tomography to resolve buoyancy anomalies in the oceanic upper mantle is complicated by the dominant influence of lithospheric cooling on velocity structure. Here, we remove this cooling signal from global surface wave tomographic models, revealing a correlation between positive residual depth and slow residual velocity anomalies at depths <300 km. To investigate whether these anomalies are of sufficient amplitude to account for short‐wavelength residual depth variations, we calibrate an experimentally derived parameterization of anelastic deformation at seismic frequencies to convert shear wave velocity into temperature, density, and diffusion creep viscosity. Asthenospheric temperature anomalies reach +150°C in the vicinity of major magmatic hot spots and correlate with geochemical and geophysical proxies for potential temperature along mid‐ocean ridges. Locally, we find evidence for a ∼150 km‐thick, low‐viscosity asthenospheric channel. Incorporating our revised density structure into models of whole‐mantle flow yields reasonable agreement with residual depth observations and suggests that ±30 km deviations in local lithospheric thickness account for a quarter of total amplitudes. These predictions remain compatible with geoid constraints and substantially improve the fit between power spectra of observed and predicted dynamic topography. This improvement should enable more accurate reconstruction of the spatiotemporal evolution of Cenozoic dynamic topography.

Quantifying Asthenospheric and Lithospheric Controls on Mafic Magmatism Across North Africa

AbstractAfrican basin‐and‐swell morphology is often attributed to the planform of subplate mantle convection. Across North Africa, the coincidence of Neogene and Quaternary (i.e., <23 Ma) magmatism, topographic swells, long wavelength gravity anomalies, and slow shear wave velocity anomalies within the asthenosphere provides observational constraints for this hypothesis. Admittance analysis of topographic and gravity fields corroborates the existence of subplate support. To investigate quantitative relationships between intraplate magmatism, shear wave velocity, and asthenospheric temperature, we collected and analyzed a suite of 224 lava samples from Tibesti, Jabal Eghei, Haruj, Sawda/Hasawinah, and Gharyan volcanic centers of Libya and Chad. Forward and inverse modeling of major, trace, and rare Earth elements was used for thermobarometric studies and to determine melt fraction as a function of depth. At each center, mafic magmatism is modeled by assuming adiabatic decompression of dry peridotite with asthenospheric potential temperatures of 1300‐1360 °C. Surprisingly, the highest temperatures are associated with the low‐lying Haruj volcanic center rather than with the more prominent Tibesti swell. Our results are consistent with earthquake tomographic models which show that the slowest shear wave anomalies within the upper mantle occur directly beneath the Haruj center. This inference is corroborated by converting observed velocities into potential temperatures, which are in good agreement with those determined by geochemical inverse modeling. Our results suggest that North African volcanic swells are primarily generated by thermal anomalies located beneath thinned lithosphere.

Quantitative Relationships Between Basalt Geochemistry, Shear Wave Velocity, and Asthenospheric Temperature Beneath Western North America

AbstractWestern North America has an average elevation that is ∼2 km higher than cratonic North America. This difference coincides with a westward decrease in average lithospheric thickness from ∼240 to <100 km. Tomographic models show that slow shear wave velocity anomalies lie beneath this region, coinciding with the pattern of basaltic magmatism. To investigate relationships between magmatism, shear wave velocity, and temperature, we analyzed a suite of >260 basaltic samples. Forward and inverse modeling of carefully selected major, trace, and rare earth elements were used to determine melt fraction as a function of depth. Basaltic melt appears to have been generated by adiabatic decompression of dry peridotite with asthenospheric potential temperatures of 1340 ± 20 °C. Potential temperatures as high as 1365 °C were obtained for the Snake River Plain. For the youngest (i.e., <5 Ma) basalts with a subplate geochemical signature, there is a positive correlation between shear wave velocities and trace element ratios such as La/Yb. The significance of this correlation is explored by converting shear wave velocity into temperature using a global empirical parameterization. Calculated temperatures agree with those determined by inverse modeling of rare earth elements. We propose that regional epeirogenic uplift of western North America is principally maintained by widespread asthenospheric temperature anomalies lying beneath a lithospheric plate, which is considerably thinner than it was in Late Cretaceous times. Our proposal accounts for the distribution and composition of basaltic magmatism and is consistent with regional heat flow anomalies.

Neogene Uplift and Magmatism of Anatolia: Insights From Drainage Analysis and Basaltic Geochemistry

AbstractIt is agreed that mantle dynamics have played a role in generating and maintaining the elevated topography of Anatolia during Neogene times. However, there is debate about the relative importance of subduction zone and asthenospheric processes. Key issues concern onset and cause of regional uplift, thickness of the lithospheric plate, and the presence/absence of temperature and/or compositional anomalies within the convecting mantle. Here, we tackle these interlinked issues by analyzing and modeling two disparate suites of observations. First, a drainage inventory of 1,844 longitudinal river profiles is assembled. This database is inverted to calculate the variation of Neogene regional uplift through time and space by minimizing the misfit between observed and calculated river profiles subject to independent calibration. Our results suggest that regional uplift commenced at 20 Ma in the east and propagated westward. Second, we have assembled a database of geochemical analyses of basaltic rocks. Two different approaches have been used to quantitatively model this database with a view to determining the depth and degree of asthenospheric melting across Anatolia. Our results suggest that melting occurs at depths as shallow as 60 km in the presence of mantle potential temperatures as high as 1400°C. There is evidence that temperatures are higher in the east, consistent with the pattern of subplate shear wave velocity anomalies. Our combined results are consistent with isostatic and admittance analyses and suggest that elevated asthenospheric temperatures beneath thinned Anatolian lithosphere have played a first‐order role in generating and maintaining regional dynamic topography and basaltic magmatism.

Constraints on dynamic topography from asymmetric subsidence of the mid-ocean ridges

Stresses from mantle convection deflect Earth's surface vertically, producing dynamic topography that is important for continental dynamics and sea-level change but difficult to observe due to overprinting by isostatic topography. For long wavelengths (∼104 km), the amplitude of dynamic topography is particularly uncertain, with mantle flow models typically suggesting larger amplitudes (>1000 m) than direct observations. Here we develop a new constraint on the amplitude of long-wavelength dynamic topography by examining asymmetries in seafloor bathymetry across mid-ocean ridges. We compare bathymetric profiles across the Mid-Atlantic Ridge (MAR) and the East Pacific Rise (EPR) and we find that the South American flank of both ridges subsides faster than its opposing flank. This pattern is consistent with dynamic subsidence across South America, supported by downwelling in the lower mantle. To constrain the amplitude of dynamic topography, we compare bathymetric profiles across both ridges after correcting bathymetry for several different models of dynamic topography with varying amplitudes and spatial patterns. We find that long-wavelength dynamic topography with an amplitude of only ∼500 m explains the observed asymmetry of the MAR. A similar model can explain EPR asymmetry but is complicated by additional asymmetrical topography associated with tectonic, crustal thickness, and/or asthenospheric temperature asymmetries across the EPR. After removing 500 m of dynamic topography, both the MAR and EPR exhibit a slower seafloor subsidence rate (∼280–290 m/Myr1/2) than previously reported. Our finding of only ∼500 m of long-wavelength dynamic topography may indicate the importance of thermochemical convection and/or large viscosity variations for lower mantle dynamics.

Continental margin subsidence from shallow mantle convection: Example from West Africa

Spatial and temporal evolution of the uppermost convecting mantle plays an important role in determining histories of magmatism, uplift, subsidence, erosion and deposition of sedimentary rock. Tomographic studies and mantle flow models suggest that changes in lithospheric thickness can focus convection and destabilize plates. Geologic observations that constrain the processes responsible for onset and evolution of shallow mantle convection are sparse. We integrate seismic, well, gravity, magmatic and tomographic information to determine the history of Neogene-Recent (<23 Ma) upper mantle convection from the Cape Verde swell to West Africa. Residual ocean-age depths of +2 km and oceanic heat flow anomalies of +16 ± 4 mW m−2 are centered on Cape Verde. Residual depths decrease eastward to zero at the fringe of the Mauritania basin. Backstripped wells and mapped seismic data show that 0.4–0.8 km of water-loaded subsidence occurred in a ∼500 × 500 km region centered on the Mauritania basin during the last 23 Ma. Conversion of shear wave velocities into temperature and simple isostatic calculations indicate that asthenospheric temperatures determine bathymetry from Cape Verde to West Africa. Calculated average excess temperatures beneath Cape Verde are >+100 °C providing ∼103 m of support. Beneath the Mauritania basin average excess temperatures are <−100 °C drawing down the lithosphere by ∼102 to 103 m. Up- and downwelling mantle has generated a bathymetric gradient of ∼1/300 at a wavelength of ∼103 km during the last ∼23 Ma. Our results suggest that asthenospheric flow away from upwelling mantle can generate downwelling beneath continental margins.

Seismic High Attenuation Region Observed Beneath Southern New England From Teleseismic Body Wave Spectra: Evidence for High Asthenospheric Temperature Without Melt

AbstractSeismic attenuation exhibits strong geographic variability in northeastern North America, with the highest values associated with the previously recognized Northern Appalachian Anomaly (NAA) in southern New England. The shear wave quality factor at 100 km depth is 14 &lt; QS &lt; 25, the ratio of P wave and S wave quality factors is QP/QS = 1.2 ± 0.03 (95%), and the frequency dependence parameter is α = 0.39 ± 0.025 (95%). The high values of QP/QS and α are compatible with laboratory measurements of unmelted rock and, in the case of α, incompatible with widespread melting. The low QS implies high mantle temperatures (~1,550–1,650°C) at 100 km depth (assuming no melt). Small‐scale variations in attenuation suggest structural heterogeneity within the NAA, possibly due to lithospheric delamination caused by asthenospheric flow.

Oriented grain growth and modification of ‘frozen anisotropy’ in the lithospheric mantle

Seismic anisotropy throughout the oceanic lithosphere is often assumed to be generated by fossilized texture formed during deformation at asthenospheric temperatures close to the ridge. Here we investigate the effect of high-temperature and high-pressure static annealing on the texture of previously deformed olivine aggregates to simulate residence of deformed peridotite in the lithosphere. Our experiments indicate that the orientation and magnitude of crystallographic preferred orientation (CPO) will evolve due to the preferential growth of grains with low dislocation densities. These observations suggest that texture and stored elastic strain energy promote a style of grain growth that modifies the CPO of a deformed aggregate. We demonstrate that these microstructural changes alter the orientation distributions and magnitudes of seismic wave velocities and anisotropy. Therefore, static annealing may complicate the inference of past deformation kinematics from seismic anisotropy in the lithosphere.

Asthenospheric kimberlites: Volatile contents and bulk compositions at 7 GPa

During ascent, kimberlites react with the lithospheric mantle, entrain and assimilate xenolithic material, loose volatiles and suffer from syn- and post-magmatic alteration. Consequently, kimberlite rocks deviate heavily from their primary melt. Experiments at 7 GPa, 1300–1480 °C, 10–30 wt% CO2 and 0.46 wt% H2O on a proposed primitive composition from the Jericho kimberlite show that saturation with a lherzolitic mineral assemblage occurs only at 1300–1350 °C for a carbonatitic melt with <8 wt% SiO2 and >35 wt% CO2. At asthenospheric temperatures of >1400 °C, where the Jericho melt stays kimberlitic, this composition saturates only in low-Ca pyroxene, garnet and partly olivine.We hence forced the primitive Jericho kimberlite into multiple saturation with a lherzolitic assemblage by adding a compound peridotite. Saturation in olivine, low- and high-Ca pyroxene and garnet was obtained at 1400–1650 °C (7 GPa), melts are kimberlitic with 18–29 wt% SiO2 + Al2O3, 22.1–24.6 wt% MgO, 15–27 wt% CO2 and 0.4–7.1 wt% H2O; with a trade-off of H2O vs. CO2 and temperature. Melts in equilibrium with high-Ca pyroxene with typical mantle compositions have ≥2.5 wt% Na2O, much higher than the commonly proposed 0.1–0.2 wt%.The experiments allow for a model of kimberlite origin in the convective upper mantle, which only requires mantle upwelling that causes melting at the depth where elemental carbon (in metal, diamond or carbide) converts to CO2 (at ∼250 km). If primary melts leading to kimberlites contain a few wt% H2O, then adiabatic temperatures of 1400–1500 °C would yield asthenospheric mantle melts that are kimberlitic (>18 wt% SiO2 + Al2O3) but not carbonatitic (<10 wt% SiO2 + Al2O3) in composition, carbonatites only forming 100–200 °C below the adiabat. These kimberlites represent small melt fractions concentrating CO2 and H2O and then acquire part of their chemical signature by assimilation/fractionation during ascent in the subcratonic lithosphere.

Constraints on early Franciscan subduction rates from 2-D thermal modeling

Mesozoic subduction along the western margin of North America in California jumped west in mid-Jurassic time to form the Franciscan accretionary wedge outboard of the Coast Range Ophiolite and the Great Valley forearc basin. Based on Mesozoic plate motion analysis for the Pacific region, overall rates of subduction are likely to have been on the order of 50 km/m.y. New and previously published Lu–Hf garnet growth ages on high-grade blocks in the Franciscan, however, suggest that early Franciscan metamorphism took place in a slowly cooling environment, characteristic of slow subduction rates. In order to assess rates of early Franciscan subduction, 2-D thermal modeling was carried out based on the following assumptions: (1) Subduction was initiated at 170 Ma, based on the age of the oldest tectonic block, and on the age of the Coast Range Ophiolite, believed to have formed in a forearc setting behind the newly initiated subduction zone. (2) At the start of subduction the hanging wall mantle wedge was at asthenospheric temperatures, consistent with initiation of subduction at the ophiolite spreading center. (3) The earliest subducted lithosphere was 10 m.y. old at the start of subduction. (4) Tectonic blocks were metamorphosed at an average pressure of 1.4 GPa, corresponding to a depth of 43 km beneath mantle rocks (assuming ρ = 3.3 g cm − 3 ). Subduction rates were investigated between 1 and 300 km/m.y., and the best fit between the ages and estimated metamorphic temperatures of the blocks was found to be 10 ± 5 km/m.y. This is five times slower than the ~ 50 km/m.y. determined from plate motion analysis for that time, and suggests that some reappraisal of the configuration and velocities of the oceanic plates may be needed.

Small-scale convection beneath oceans and continents

Small-scale convection supplies heat flow of ∼17 mW m−2 to the base of stable continents where xenolith studies resolve the geotherm. However, effects of small-scale convection are difficult to resolve in ocean basins. On first pass, most seafloor appears to subside to an asymptote compatible with ∼40 mW m−2 convective heat flow. These common regions are tracked by hotspots so uplift associated with ponded mantle material is an attractive alternative. Unaffected seafloor in the North and South Atlantic continues to subside with the square root of age as expected from pure conduction. The theory of stagnant-lid convection provides good scaling relationships for heat flow. For linear viscosity, heat flow is proportional to the underlying “half-space” viscosity to the −1/3 power and the temperature to change viscosity by a factor of e to the 4/3 power. The formalism is easily modified to represent convection beneath a lid of highly viscous and buoyant cratonal lithosphere and to represent transient convection beneath thickening oceanic lithosphere. Asthenospheric mantle with linear, strongly temperature-dependent, and weakly depth-dependent viscosity is compatible with both oceanic and continental data. More complicated rheology may allow vigorous small-scale convection under most but not all old ocean basins. Still viable hypotheses require poorly understood global features, including lateral variations of asthenospheric temperature. Seismological studies have the potential to resolve the lithosphere-asthenosphere boundary, including local variations of its depth associated with small-scale convection.

A Neogene chronology of Iceland plume activity from V-shaped ridges

The V-shaped ridges, which straddle the mid-oceanic ridges on either side of Iceland, are recognized as a window into transient convective circulation of the mantle. The geometry of these ridges places important constraints on the nature of asthenospheric flow within the head of the Icelandic plume. Associated changes in the bathymetry of the Greenland–Scotland Ridge may have modulated the overflow of North Atlantic Deep Water. In order to improve the Neogene chronology of these V-shaped ridges, we have re-analyzed a set of legacy seismic reflection profiles which cross the Reykjanes Ridge south of Iceland. After converting these profiles to depth, the sediment–basement interface was carefully picked and corrected for sediment loading. Ages were assigned using magnetic anomalies which were correlated to an astronomically tuned timescale. We then calculated the residual height of oceanic crust along each profile by removing components of long-wavelength subsidence which are caused by lithospheric plate cooling and by dynamic support of the plume. By correlating the pattern of residual height with free-air gravity anomalies, we can discriminate between the effects of normal faulting and crustal thickness variations. The pattern of residual height primarily reflects asthenospheric temperature fluctuations which flow away from the center of the Icelandic plume. If we assume that asthenospheric flow is radial, our results can be used to estimate the mass flux of the Icelandic plume which varies from 1–3 × 10 14 kg/yr over the last 17 million years. Earlier, during Paleogene times, there is indirect evidence for higher mass flux. If the plume conduit beneath Iceland has a radius of 150 km, our inferred mass flux corresponds to an upwelling rate of 270 mm/yr. Finally, we have compared the Neogene record of residual height predicted for the Denmark Straits with estimates of the overflow of Northern Component Water, the Neogene precursor of North Atlantic Deep Water. There is a striking correlation between the two records over the last 5 million years which suggests that regional, plume-driven, vertical motion controlled deep-water overflow at the Denmark Straits.

Evaluating hot spot–ridge interaction in the Atlantic from regional‐scale seismic observations

We probe variations in mantle temperature, composition, and fabric along hot spot–influenced sections of the Mid‐Atlantic Ridge (MAR), using surface waves from nearby ridge earthquakes recorded on broadband island‐based seismic stations. We invert frequency‐dependent phase delays from these events to estimate one‐dimensional mean shear velocity and radial shear anisotropy profiles in the upper 200 km of the mantle within two seafloor age intervals: 5–10 Ma and 15–20 Ma. Mean shear velocity profiles correlate with apparent hot spot flux: lithosphere formed near the low‐flux Ascension hot spot is characterized by high mantle velocities, while the MAR near the higher‐flux Azores hot spot has lower velocities. The impact of the high‐flux Iceland hot spot on mantle velocities along the nearby MAR is strongly asymmetric: the lithospheric velocities near the Kolbeinsey ridge are moderately slow, while velocities near the Reykjanes ridge estimated in previous studies are much slower. Within each region the increase in shear velocity with age is consistent with a half‐space cooling model, and the velocity variations observed between Ascension, the Azores, and Kolbeinsey are consistent with approximately ±75° potential‐temperature variation among these sites. In comparison, the Reykjanes lithosphere is too slow to result purely from half‐space cooling of a high‐temperature mantle source. We speculate that the anomalously low shear velocities within the lithosphere produced at the Reykjanes ridge result from high asthenospheric temperatures of +50–75 K combined with ∼12% (by volume) gabbro retained in the mantle due to the imbalance between high hot spot–influenced melt production and relatively inefficient melt extraction along the slow spreading Reykjanes. Radial shear anisotropy in the upper 150 km also indicates an apparent hot spot influence: mantle fabric near Ascension is quite weak, consistent with previous models of anisotropy produced by corner flow during slow seafloor spreading. The fabric near the Azores and the Kolbeinsey ridge is stronger, suggesting that the hot spot increases mantle deformation beyond that produced by slow seafloor spreading in these regions.

Evidence for thin oceanic crust on the extinct Aegir Ridge, Norwegian Basin, NE Atlantic derived from satellite gravity inversion

Satellite gravity inversion incorporating a lithosphere thermal gravity correction has been used to map crustal thickness and lithosphere thinning factor for the NE Atlantic. Predicted oceanic crustal thicknesses in the Norwegian Basin are between 4 and 7 km on the extinct Aegir Ridge, increasing to 9 – 14 km at the margins, consistent with volcanic margin continental breakup at the end of the Paleocene. The observation (from gravity inversion and seismic refraction studies) of thin oceanic crust produced by the Aegir Ridge in the Oligocene may have implications for the temporal evolution of asthenosphere temperature under the NE Atlantic during the Tertiary. Thin Oligocene oceanic crust may imply cool (normal) asthenosphere temperatures during the Oligocene in contrast to elevated asthenosphere temperatures in the Paleocene and Miocene‐Recent as indicated by the formation of volcanic margins and Iceland respectively.

Petrochemistry of peridotites from North China: Significance for lithospheric mantle evolution

A petrochemical analysis was undertaken of peridotitic xenoliths in volcanic rocks that erupted from the Paleozoic to the Cenozoic within the eastern part of the North China craton, and the peridotites as tectonic intrusion in the Early Mesozoic from the Sulu orogen. The results show that the cratonic mantle, which was refractory and existed when the kimberlites intruded in the Paleozoic, had almost been replaced by the newly accreted fertile lithospheric mantle during the Mesozoic-Cenozoic. The erosion, metasomatism, and intermingling caused by the accreted asthenospheric material acting on the craton mantle along the weak zone and deep fault (such as the Tanlu fault) in the existing lithosphere resulted in the lithospheric thinning at a larger scale 100 Ma ago (but later than 178 Ma). The largest thinning would be in the Eogene. The upwelling asthenospheric material transformed into accreted lithospheric mantle due to the asthenospheric temperature falling in the Neogene (leading to relatively slight lithospheric incrassation), and finally accomplished mantle replacement. The peridotitic body in the Sulu orogen represents the products spreading from the modified cratonic lithospheric mantle.

Fractionation of the noble metals by physical processes

During partial melting in the earth’s mantle, the noble metals become fractionated. Os, Ir, Ru, and Rh tend to remain in the mantle residue whereas Pt, Pd, and Re behave mildly incompatible and are sequestered to the silicate melt. There is consensus that sulfide plays a role in the fractionation process; the major noble metal repository in the mantle is sulfide, and most primitive mantle melts are sulfide-saturated when they leave their mantle sources. However, with sulfide–silicate partitioning, the fractionation cannot be modeled properly. All sulfide–silicate partition coefficients are so extremely high that a silicate melt segregating from a mantle source with residual sulfide should be largely platinum-group elements free. We offer a physical alternative to sulfide–silicate chemical partitioning and provide a mechanism of generating a noble metal-rich melt from a sulfide-saturated source: Because sulfide is at least partially molten at asthenospheric temperature, it will behave physically incompatible during melt segregation, and a silicate melt segregating from a mantle residue will entrain molten residual sulfide in suspension and incorporate it in the basaltic pool melt. The noble metal abundances of a basalt then become independent of sulfide–silicate chemical partitioning. They reflect the noble metal abundances in the drained sulfide fraction as well as the total amount of sulfide entrained. Contrary to convention, we suggest that a fertile, sulfide-rich mantle source has more potential to generate a noble metal-enriched basaltic melt than a refractory mantle source depleted by previous partial melting events.