Evolution of lithospheric mantle in the north of Nain‐Baft oceanic crust (Neo‐Tethyan ophiolite of Ashin,Central Iran)

This study is focused on a plagioclase‐bearing spinel lherzolite from Chah Loqeh area in the Neo‐Tethyan Ashin ophiolite. It is exposed along the west of left‐lateral strike‐slip Dorouneh Fault in the northwest of Central‐East Iranian Microcontinent. Mineral chemistry (Mg#^olivine < ~ 90, Cr#^{clinopyroxene} < ~ 0.2, Cr#^spinel < ~ 0.5, Al₂O₃^{orthopyroxene} > ~ 2.5 wt%, Al₂O₃^{clinopyroxene} > ~ 4.5 wt%, Al₂O₃^spinel > ~ 41.5 wt%, Na₂O^{clinopyroxene} > ~ 0.11 wt%, and TiO₂^{clinopyroxene} > ~ 0.04 wt%) confirms Ashin lherzolite was originally a mid‐oceanic ridge peridotite with low degrees of partial melting at spinel‐peridotite facies in a lithospheric mantle level. However, some Ashin lherzolites record mantle upwelling and tectonic exhumation at plagioclase‐peridotite facies during oceanic extension and diapiric motion of mantle along Nain‐Baft suture zone. This mantle upwelling is evidenced by some modifications in the modal composition (i.e. subsolidus recrystallization of plagioclase and olivine between pyroxene and spinel) and mineral chemistry (e.g. increase in TiO₂ and Na₂O of clinopyroxene, and TiO₂ and Cr# of spinel and decrease in Mg# of olivine), as a consequence of decompression during a progressive upwelling of mantle. Previous geochronological and geochemical data and increasing the depth of subsolidus plagioclase formation at plagioclase‐peridotite facies from Nain ophiolite (~ 16 km) to Ashin ophiolite (~ 35 km) suggest a south to north closure for the Nain‐Baft oceanic crust in the northwest of Central‐East Iranian Microcontinent.     

Formation time of the big mantle wedge beneath eastern China and a new lithospheric thinning mechanism of the North China craton—Geodynamic effects of deep recycled carbon

High-resolution P wave tomography shows that the subducting Pacific slab is stagnant in the mantle transition zone and forms a big mantle wedge beneath eastern China. The Mg isotopic investigation of large numbers of mantle-derived volcanic rocks from eastern China has revealed that carbonates carried by the subducted slab have been recycled into the upper mantle and formed carbonated peridotite overlying the mantle transition zone, which becomes the sources of various basalts. These basalts display light Mg isotopic compositions(δ26 Mg = –0.60‰ to –0.30‰) and relatively low87 Sr/86 Sr ratios(0.70314–0.70564) with ages ranging from 106 Ma to Quaternary, suggesting that their mantle source had been hybridized by recycled magnesite with minor dolomite and their initial melting occurred at 300-360 km in depth. Therefore, the carbonate metasomatism of their mantle source should have occurred at the depth larger than 360 km, which means that the subducted slab should be stagnant in the mantle transition zone forming the big mantle wedge before 106 Ma. This timing supports the rollback model of subducting slab to form the big mantle wedge. Based on high P-T experiment results, when carbonated silicate melts produced by partial melting of carbonated peridotite was raising and reached the bottom(180–120 km in depth) of cratonic lithosphere in North China, the carbonated silicate melts should have 25–18 wt% CO2 contents, with lower Si O2 and Al2 O3 contents, and higher Ca O/Al2 O3 values, similar to those of nephelinites and basanites, and have higher εNdvalues(2 to 6). The carbonatited silicate melts migrated upward and metasomatized the overlying lithospheric mantle, resulting in carbonated peridotite in the bottom of continental lithosphere beneath eastern China. As the craton lithospheric geotherm intersects the solidus of carbonated peridotite at 130 km in depth, the carbonated peridotite in the bottom of cratonic lithosphere should be partially melted, thus its physical characters are similar to the asthenosphere and it could be easily replaced by convective mantle. The newly formed carbonated silicate melts will migrate upward and metasomatize the overlying lithospheric mantle. Similarly, such metasomatism and partial melting processes repeat, and as a result the cratonic lithosphere in North China would be thinning and the carbonated silicate partial melts will be transformed to high-Si O2 alkali basalts with lower εNdvalues(to-2). As the lithospheric thinning goes on,initial melting depth of carbonated peridotite must decrease from 130 km to close 70 km, because the craton geotherm changed to approach oceanic lithosphere geotherm along with lithospheric thinning of the North China craton. Consequently, the interaction between carbonated silicate melt and cratonic lithosphere is a possible mechanism for lithosphere thinning of the North China craton during the late Cretaceous and Cenozoic. Based on the age statistics of low δ26 Mg basalts in eastern China, the lithospheric thinning processes caused by carbonated metasomatism and partial melting in eastern China are limited in a timespan from 106 to25 Ma, but increased quickly after 25 Ma. Therefore, there are two peak times for the lithospheric thinning of the North China craton: the first peak in 135-115 Ma simultaneously with the cratonic destruction, and the second peak caused by interaction between carbonated silicate melt and lithosphere mainly after 25 Ma. The later decreased the lithospheric thickness to about70 km in the eastern part of North China craton.     

Diapirism of depleted peridotite — a model for the origin of hot spots

A model is proposed for the origin of hot spots that depends on the existence of major-element heterogeneities in the mantle. Generation of basaltic crust at spreading centers produces a layer of residual peridotite ～20–25 km thick directly beneath the crust which is depleted in Fe/Mg, TiO₂, CaO, Al₂O₃, Na₂O and K₂O, and which has a slightly lower density than undepleted peridotite beneath it. Upon recycling of this depleted peridotite back into the deep mantle at subduction zones, it becomes gravitationally unstable, and tends to rise as diapirs through undepleted peridotite. For a density contrast of 0.05 g cm^?3, a diapir 60 km in diameter would rise at roughly 8 cm y^?1, and could transport enough heat to the base of the lithosphere to cause melting and volcanism at the surface. Hot spots are thus viewed as a passive consequence of mantle convection and fractionation at spreading centers rather than a plate-driving force.It is suggested that depleted diapirs exist with varying amounts of depletion, diameters, upward velocities and source volumes. Such variations could explain the occurrence of hot spots with widely varying lifetimes and rates of lava production. For highly depleted diapirs with very low Fe/Mg, the diapir would act as a heat source and the asthenosphere and lower lithosphere drifting across the diapir would serve as the source region of magmas erupted at the surface. For mildly depleted diapirs with Fe/Mg only slightly less than in normal undepleted mantle, the diapir could provide not only the source of heat but also most or all of the source material for the erupted magmas. The model is consistent with isotopic data that require two separate and ancient source regions for mid-ocean ridge and oceanic island basalts. The source for mid-ocean ridge basalts is considered to be material upwelling at spreading centers from the deep mantle. This material forms the oceanic lithosphere. Oceanic island basalts are considered to be derived from varying mixtures of sublithospheric and lower lithospheric material and the rising diapir itself.     

The role of lithospheric phase transitions on seafloor flattening at old ages

Uplift of old ocean seafloor relative to the familiar t^1/2 curve is generally presumed to be a consequence of an additional heat source at the base of the oceanic lithosphere. Several geodynamical mechanisms have been invoked to explain this behavior. We have taken a petrological rather than geodynamical approach and considered whether the causes of the flattening can be found within the conductive portion of the lithosphere, where viscous creep is insignificant. Accordingly, self-consistent calculations of the phase assemblages of several candidate mantle compositions have been performed, subject to the constraint of experimentally determined calorimetric measurements. We have found that, in general, for temperatures below around 800°C and pressures greater than 10 kbar, the phase change from spinel peridotite to garnet peridotite has a significantly negative Clapeyron slope, which could cause some amount of uplift of old seafloor. The transition, for representative thermal profiles, should be depressed from about 35 km (at 40 Ma) to 45 km (at 160 Ma). Since spinel peridotite is of lower density the net effect is to raise the seafloor topography. The extent of the uplift, which depends primarily on the Al₂O₃ content of the peridotite and on the effective thermal diffusivity of the lithosphere, should be on the order of 140–220 m.     

Geodynamic significance of ultramafic xenoliths from Eastern Serbia: Relics of sub-arc oceanic mantle?

A suite of highly depleted peridotite xenoliths in East Serbian Palaeogene basanites represents the lithospheric mantle beneath the Balkan Peninsula. The xenoliths are harzburgites, clinopyroxene-poor lherzolites and rare dunites. They contain mostly <5 vol.% of modal clinopyroxene and are characterized by high Mg# in silicates (>91), high Cr# in spinel (mostly 0.5–0.7), and by distinctively low Al₂O₃ contents in orthopyroxene (mostly 1–2 wt.%). They have experienced some mantle metasomatism which has slightly obscured their original composition. Nevertheless, the general characteristics of the xenoliths imply a composition which is significantly more depleted than most non-cratonic sub-continental mantle xenolith suites, as well as orogenic peridotites and abyssal peridotites. Geological and compositional evidence suggests that the xenoliths do not represent Archean mantle. The existence of Proterozoic mantle cannot be entirely excluded, although it is in disagreement with geological evidence. On the other hand, the studied xenoliths are compositionally very similar to peridotites of modern oceanic sub-arc settings. The existence of such a depleted lithospheric mantle segment is also inferred from the presence of rare orthopyroxene-rich xenoliths in the same suite. These are interpreted to have originated as lithospheric precipitates of high-Mg, SiO₂-saturated magmas that require a highly depleted mantle source. Such source is typically required by boninitic-like magmas of intraoceanic suprasubduction settings. A proposed geodynamic model to explain these observations involves accretion or underplating of the lower parts of the Tethyan oceanic lithosphere during the Upper Jurassic closure of the eastern branch of the Vardar ocean.     

Composition,temperature, and thickness of the lithosphere of the Archean Kaapvaal craton

A new model is proposed for the structure of the Kaapvaal craton lithosphere. Based on chemical thermodynamics methods, profiles of the chemical composition, temperature, density, and S wave velocities are constructed for depths of 100–300 km. A solid-state zone of lower velocities is discovered on the S velocity profile in the depth interval 150–260 km. The temperature profiles are obtained from absolute values of P and S velocities, taking into account phase transformations, anharmonicity, and anelastic effects. The examination of the sensitivity of seismic models to the chemical composition showed that relatively small variations in the composition of South African xenoliths result in lateral temperature variations of ～200°C. Inversion of some seismic profiles (including IASP91) with a fixed bulk composition of garnet peridotites (the primitive mantle material) leads to a temperature inversion at depths of 200–250 km, which is physically meaningless. It is supposed that the temperature inversion can be removed by gradual fertilization of the mantle with depth. In this case, the craton lithosphere should be stratified in chemical composition. The depleted lithosphere composed by garnet peridotites exists to depths of 175–200 km. The lithospheric material at depths of 200–250 km is enriched in basaltoid components (FeO, Al₂O₃, and CaO) as compared with the material of garnet peridotites but is depleted in the same components as compared with the fertile substance of the underlying primitive mantle. The material composing the craton root at a depth of ～275 km does not differ in its physical and chemical characteristics from the composition of the normal mantle, and this allows one to estimate the thickness of the lithosphere at 275 km. The results of this work are compared with data of seismology, thermal investigations, and thermobarometry.     

The role of magma mixing in the petrogenesis of mafic alkaline lavas,Rockeskyllerkopf Volcanic Complex,West Eifel,Germany

The quaternary Rockeskyllerkopf Volcanic Complex (RVC) comprises three spatially and temporally distinct volcanic centers that can also be distinguished on the basis of their geochemical signatures. All the volcanic products in the complex are olivine basanites whose major and trace element compositions span almost the entire range defined for the West Eifel field as a whole. The RVC lavas have lower Al₂O₃, Na₂O and Y contents and higher TiO₂, CaO, K₂O, Sc, V, Co, Rb, and Ba than the Tertiary lavas in nearby Hocheifel volcanic field. Within the complex, the oldest South East Lammersdorf Center (SEL) comprises primitive lavas with an average MgO content of ～11 wt.% and La_N/Yb_N of 29?±?2. The second center, Mäuseberg, has similar MgO to SEL but is distinct in its much higher La_N/Yb_N of 42?±?2. The Rockeskyllerkopf Center, which was erupted after a break in activity, comprises lavas similar in composition to the SEL Center but with distinctly higher Al₂O₃ and lower MgO contents. Given the lack of evidence for significant fractionation or assimilation in the RVC lavas, we attribute the compositional variations within and between the centers of the RVC to be due to variations in the composition of the source region in combination with magma mixing. Our preferred model involves 1–5% partial melting of LREE-enriched mantle in the garnet stability field, likely within the thermal boundary layer at the base of the lithospheric mantle. These melts mixed to variable degrees with 2–4% partial melts of phlogopite-spinel peridotite formed at higher levels in the modally metasomatised lithospheric mantle.     

The mantle origins of Karoo picrites

Karoo picrites are divided into a volumetrically dominant high Ti-Zr (HTZ) group and a less abundant low Ti-Zr (LTZ) group. Distinguished in this way, it is then possible to use the comparison of major element compositions with those of experimentally determined partial melts of mantle peridotites to constrain the depth of segregation for each magma type. This approach also identifies the nature of source and residual mantle materials from which such magma segregation may have occurred. The LTZ group show uniformly shallow pressures (ca. 13–15 kb) of magma segregation and the HTZ group show a range of pressures (10 kb to > 30 kb) of segregation. In the HTZ group a range of K₂O, TiO₂ and Zr contents is observed, with greater concentrations in picrites which have segregated at higher pressures (high-NaK# HTZ picrites) relative to a low-NaK# HTZ picrite group (NaK#=[Na/₂O+K₂O][Na₂O+K₂O+CaO]). If we measure the refractory character of a mantle source peridotite by its CaO and Al₂O₃ content then it is shown that the LTZ picrites have segregated from sources which are more refractory than MORB or OIB sources and which lack the strong K, Ti and Zr enrichments of HTZ sources. Mixing of at least two mantle components is required to fully explain the major element and isotopic data in the HTZ group and is combined with variations in composition produced by different depths of magma segregation. Consideration of the major element phase relations coupled with trace element abundances and isotope data leads us to suggest that one component is the Kaapvaal lithospheric mantle (incompatible element enriched but relatively refractory with respect to basaltic components) and the other is asthenospheric mantle. This lithospheric mantle component is most evident in the high-NaK# HTZ picrite end-member which has segregated at a greater depth (ca. 18–32 kb). The asthenospheric component is most closely represented by the low-NaK# HTZ picrites which have segregated at pressures of 10–22 kb. The LTZ source, more refractory mantle than MORB, is inferred to be the shallow sub-continental lithospheric mantle which may also be the source of the low-MgO tholeiites related to these picrites, and therefore the source for the main phase of Karoo basaltic eruptions.     

Major element chemistry of ocean island basalts — Conditions of mantle melting and heterogeneity of mantle source

We estimate average compositions of near-primary, ‘reference’ ocean island basalts (OIBs) for 120 volcanic centers from 31 major island groups and constrain the depth of lithosphere–asthenosphere boundary (LAB) at the time of volcanism and the possible depth of melt–mantle equilibration based on recently calibrated melt silica activity barometer. The LAB depth versus fractionation corrected OIB compositions (lava compositions, X, corrected to Mg# 73, X_OIB^#73, i.e., magmas in equilibrium with Fo₉₀, if olivine is present in the mantle source) show an increased major element compositional variability with increasing LAB depths. OIBs erupted on lithospheres < 40 km thick approach the compositions (e.g. SiO₂^#73, TiO₂^#73, [CaO/Al₂O₃]^#73) of primitive ridge basalts and are influenced strongly by depth and extent of shallow melting. However, X_OIB^#73 on thicker lithospheres cannot be explained by melt–mantle equilibration as shallow as LAB. Melt generation from a somewhat deeper (up to 50 km deeper than the LAB) peridotite source can explain the OIB major element chemistry on lithospheres ≤ 70 km. However, deeper melting of volatile-free, fertile peridotite is not sufficient to explain the end member primary OIBs on ≥ 70 km thick lithospheres. Comparison between X_OIB^#73 and experimental partial melts of fertile peridotite indicates that at least two additional melt components need to be derived from OIB source regions. The first component, similar to that identified in HIMU lavas, is characterized by low SiO₂^#73, Al₂O₃^#73, [Na₂O/TiO₂]^#73, and high FeO?^#73, CaO^#73, [CaO/Al₂O₃]^#73. The second component, similar to that found in Hawaiian Koolau lavas, is characterized by high SiO₂^#73, moderately high FeO?^#73, and low CaO^#73 and Al₂O₃^#73. These two components are not evenly sampled by all the islands, suggesting a heterogeneous distribution of mantle components that generate them. We suggest that carbonated eclogite and volatile-free, silica-excess eclogite are the two most likely candidates, which in conjunction with fertile mantle peridotite, give rise to the two primitive OIB end members.     

Petrology of the Green Knobs diatreme and implications for the upper mantle below the Colorado Plateau

Peridotite inclusions, crystal fragments, and kimberlite breccia at Green Knobs, New Mexico, have been studied to evaluate compositions and processes in the upper mantle below the Colorado Plateau. Most peridotite inclusions are spinel lherzolites and harzburgites, or their partly hydrated equivalents, in the Cr-diopside group. Orthopyroxene-rich websterites and olivine websterites comprise 3% of the peridotites and formed as cumulates. Typical anhydrous or slightly hydrated peridotites contain aluminous, calcic diopside (5–7% Al₂O₃), aluminous orthopyroxene (3–6% Al₂O₃), spinel, and olivine (near Fa₉). Geothermometers based on different mineral pairs yield temperatures from above 1100°C to below 700°C in single rocks. High values, derived from pyroxenes with included exsolution lamellae, may approximate temperatures of primary crystallization. Low values, based on olivine-spinel and olivine-clinopyroxene pairs, approach upper mantle temperatures before eruption. In rare samples, some spinel grains are rimmed by garnet while others are not rimmed; garnet formation was controlled by nucleation kinetics. About one-third of the peridotites were deformed shortly before eruption, with effects ranging from mild cataclasis to the production of ultramylonites.Discrete crystals of garnet, olivine (near Fa₈), and Cr-diopside represent garnet peridotite. Eclogites were not found. The garnet peridotite is more depleted than overlying spinel peridotite, and it is not a likely source for the minettes associated with the kimberlites.The mantle below Green Knobs consists of spinel peridotite from 45 to perhaps 60 km depth immediately underlain by more-depleted garnet peridotite. The position of the spinel-garnet transition may be fixed by kinetics. The kimberlite may have been produced when heat from ascending minette magma released volatiles from otherwise depleted garnet peridotite. Resulting gas-solid mixtures erupted along zones of deformation associated with Colorado Plateau monoclines. Sheared lherzolites formed during renewed movement along these zones.     

Thermal structure and thickness of the lithospheric mantle underlying the Siberian Craton from the kraton and kimberlit superlong seismic profiles

A self-consistent approach is proposed for the investigation of the thermal conditions, chemical composition, and internal structure of the upper mantle of the Earth. Using this approach, the thermal state of the lithospheric mantle beneath the Siberian Craton (SC) is reconstructed from P velocities, taking into account the phase transitions, anharmonicity, and the effects of anelasticity. The velocities of seismic waves are more sensitive to temperature than to the composition of the mantle rocks, which allows the velocity models to be effectively used for reconstruction of the thermal regime of the mantle. The temperature at depths 100–300 km is reconstructed by inversion of the Kraton and Kimberlit superlong seismic profiles for compositions of the garnet harzburgite, lherzolite, and intermediate composition of garnet peridotite. The averaged temperature in the normal continental mantle is reconstructed by inversion of the IASP91 reference model for depleted and fertile substance. One-dimensional models and two-dimensional thermal fields undergo a substantial fall in temperature (～300–600°C) beneath the Siberian Craton as compared to the temperatures of the continental mantle and paleotemperatures inferred from the thermobarometry of xenoliths. Temperature profiles of the Siberian Craton deduced from seismic data lie between the conductive geotherms of 32.5–40.0 mW/m² and below the P(H)-T values obtained for low- and high-temperature xenoliths from the Mir, Udachnaya, and Obnazhennaya kimberlite pipes. The thickness of the thermal lithosphere estimated from the intersection with the potential adiabat is 300–320 km, which is consistent with the data on heat flows and seismotomographic observations. This provides grounds for the assumption that the low-temperature anomalies (thermal roots of continents) penetrate down to a depth of 300 km. The analysis of the sensitivity of seismic velocity and density to the variations in temperature, pressure, and chemical and phase composition of petrological models shows that recognition of fine differences in chemical composition of the lithospheric rocks by seismic methods is impossible.     

Effective elastic thicknesses of the lithosphere in the Central Iberian Peninsula from heat flow: Implications for the rheology of the continental lithospheric mantle

The traditional view of the rheology of the continental lithosphere, sometimes known as the “jelly sandwich model”, consists of a strong upper crust, a weak lower crust, and a strong upper lithospheric mantle. Some authors argue, however, that the lithospheric mantle is weak and contributes little to the total strength and the effective elastic thickness of the lithosphere; this weakness is claimed to be due to the mantle being wet or subjected to temperatures higher than usually believed. This paper uses the relationship between rheology of the lithosphere and heat flow to calculate theoretical effective elastic thicknesses for three regions of the central Iberian Peninsula (the Duero Basin, the Spanish Central System and the Tajo Basin), taking into account the contribution of the crust and the lithospheric mantle, for dry and wet rheologies. We found that a wet peridotite rheology for the lithospheric mantle is generally consistent with independent (based on Bouguer coherence or flexural modeling) estimates of the effective elastic thickness for the study area, whereas a dry peridotite rheology cannot be reconciled with them. Moreover, the contribution of the mantle to the bending moment of the lithosphere, and therefore to both the effective elastic thickness and the total strength of the lithosphere, is important, and it may even be the dominant contribution. Therefore, the jelly sandwich model may be considered valid for the central Iberian Peninsula.     

A Nd and Sr isotopic study of the Trinity peridotite; implications for mantle evolution

Field evidence indicates that the Trinity peridotite was partially melted during its rise as a part of the upwelling convecting mantle at a spreading center. A SmNd mineral isochron for a plagioclase lherzolite yields an age,T = 427 ± 32 Ma and initialε_Nd = + 10.4 ? 0.4 which is distinctly higher than that expected for typical depleted mantle at this time. This age is interpreted as the time of crystallization of trapped melt in the plagioclase lherzoliteP-T field. This time of crystallization probably represents the time when the massif was incorporated as a part of the oceanic lithosphere. The SmNd model age of the plagioclase lherzolite totalrock isT_CHUR^Nd = 3.4 AE. This suggests that the Trinity peridotite was derived from a mantle that was depleted rather early in earth history. The peridotite contains many generations of pyroxenite dikes and some microgabbro dikes. We report data for two dikes that clearly crosscut the main metamorphic fabric of the peridotite. A microgabbro dike yields a SmNd mineral isochron age ofT = 435 ± 21 Ma andε_Nd = + 6.7 ? 0.3. A pyroxenite dike yields an initialε_Nd = + 7.3 ± 0.4. The initialε_Nd values for the pyroxenite and gabbro dikes are fairly similar to those for the depleted mantle at this time and are distinct from the lherzolite—demonstrating that they are not genetically related. RbSr data do not give any coherent pattern. However, some bounds can be put on initial Sr values ofε_Sr ? ?21 for the plagioclase lherzolite andε_Sr ? ?8.7 for the microgabbro dike. It is plausible that the dikes represent cumulates left behind from island arc magmas that rose through the the oceanic lithosphere within the vicinity of a subduction zone. Major and trace elements and SmNd isotopic data indicate a multiple stage history for the Trinity peridotite; a small melt fraction was extracted from an undepleted source ～ 3.4 AE or more ago to produce the proto-lherzolite; a large fraction of melt (～ 12 to 23%) was extracted from the proto-lherzolite to produce the present rock; the lherzolite was then crosscut by dikes from average depleted mantle ～ 0.44 AE ago. The data are compatible with the depleted mantle source being formed very early in earth history. Although most available data indicate that the depleted upper mantle has been relatively well stirred through time, the Trinity data suggest that very ancient Nd isotopic values are preserved and thus chemical and physical heteorgeneities are sometimes preserved in the depleted source of mid-ocean ridge basalts as well as the oceanic lithosphere which they intrude.     

Is upper-mantle phosphorus contained in sodic garnet?

Garnets crystallized experimentally from within the anhydrous melting ranges of an olivine tholeiite, a tholeiitic andesite and an augite leucitite at pressures between 18 and 45 kbars contain up to 0.4% Na₂O and 0.6% P₂O₅. The Na and P are thought to form a substitution couple, replacing Ca and Si in the garnet structure; representing limited solid solution between grossular (Ca₃Al₂Si₃O₁₂) and the phosphate Na₃Al₂P₃O₁₂. This substitution is enhanced by increasing pressure and by falling temperature (increasing degree of crystallization) at constant pressure.Current knowledge of the crystalline site of P in the upper mantle is hampered by lack of data on the stability of apatite and other phosphates at appropriate pressures and temperatures. If all samples of garnetiferous upper mantle brought to the surface by magmatic processes have been depleted to some extent by previous escape of a partial-melt fraction, P₂O₅ concentrations below 0.1% in their garnets could nevertheless signify that this phase was the sole predepletion host for P in the upper mantle, at the depths from which such inclusions are derived. If garnet and apatite are the principal minerals containing P in the upper mantle, it may be possible to use covariances between P and rare-earth elements in mafic liquids to detect which of these phases was the dominant host for P at the site of magma genesis. This approach confirms the widely-held opinion that strongly alkalic mafic magmas are products of upper-mantle partial fusion in the presence of residual garnet. It also leads to a contrasting proposal that mid-ocean ridge basalts may be generated by upper-mantle partial fusion at comparatively small depths, in the presence of residual apatite.     

Si- and alkali-rich melt inclusions in minerals of mantle peridotites from eastern China: Implication for lithospheric evolution

CO2 fluid inclusions in mantle minerals are an im-portant source for us to get the information of mantle fluids. Fluid inclusions are mainly composed of CO2, with minor CO, H2O, CH4, N2, H2S, SO2, F, etc., which were demonstrated by lots of Raman spec-trometer analyses in recent years. In contrast, there are very few researches on CO2-bearing melt inclusions since it is more difficult to do so. The available studies have found that the primary CO2-bearing melt inclu-sions are basaltic …     

Garnet-ilmenite-perovskite transitions in the system Mg4Si4O12-Mg3Al2Si3O12 at high pressures and high temperatures: phase equilibria, calorimetry and implications for mantle structure

Phase relations in the system Mg₄Si₄O₁₂-Mg₃Al₂Si₃O₁₂ were examined at pressures of 19-27 GPa and relatively low temperatures of 800-1000 °C using a multianvil apparatus to clarify phase transitions of pyroxene-garnet assemblages in the mantle. Both of glass and crystalline starting materials were used for the experiments. At 1000 °C, garnet solid solution (s.s.) transforms to aluminous ilmenite s.s. at 20-26 GPa which is stable in the whole compositional range in the system. In Mg₄Si₄O₁₂-rich composition, ilmenite s.s. transforms to a single-phase aluminous perovskite s.s., while Mg₃Al₂Si₃O₁₂-rich ilmenite s.s. dissociates into perovskite s.s. and corundum s.s. These newly determined phase relations at 1000 °C supersede preliminary phase relations determined at about 900 °C in the previous study. The phase relations at 1000 °C are quite different from those reported previously at 1600 °C where garnet s.s. transforms directly to perovskite s.s. and ilmenite is stable only very close to Mg₄Si₄O₁₂. The stability field of Mg₃Al₂Si₃O₁₂ ilmenite was determined at 800-1000 °C and 25-27 GPa by reversed phase boundaries. In ilmenite s.s., the a-axis slightly increases but the c-axis and molar volume decrease substantially with increasing Al₂O₃ content. Enthalpies of ilmenite s.s. were measured by differential drop-solution calorimetry method using a high-temperature calorimeter. The excess enthalpy of mixing of ilmenite s.s. was almost zero within the errors. The measured enthalpies of garnet-ilmenite and ilmenite-perovskite transitions at 298 K were 105.2±10.4 and 168.6±8.2 kJ/mol, respectively, for Mg₄Si₄O₁₂, and 150.2±15.9 and 98.7±27.3 kJ/mol, respectively, for Mg₃Al₂Si₃O₁₂. Thermodynamic calculations using these data give rise to phase relations in the system Mg₄Si₄O₁₂-Mg₃Al₂Si₃O₁₂ at 1000 and 1600 °C that are generally consistent with those determined experimentally, and confirm that the single-phase field of ilmenite expands from Mg₄Si₄O₁₂ to Mg₃Al₂Si₃O₁₂ with decreasing temperature. The earlier mentioned phase relations in the simplified system as well as those in the Mg₂SiO₄-Fe₂SiO₄ system are applied to estimate mineral proportions in pyrolite as a function of depth along two different geotherms: one is a horizontally-averaged temperature distribution in a normal mantle, and the other being 600 °C lower than the former as a possible representative geotherm in subducting slabs. Based on the previously described estimated mineral proportions versus depth along the two geotherms, density and compressional and shear wave velocities are calculated as functions of depth, using available mineral physics data. Along a normal mantle geotherm, jumps of density and velocities at about 660 km corresponding to the post-spinel transition are followed by steep gradients due to the garnet-perovskite transition between 660 and 710 km. In contrast, along a low-temperature geotherm, the first steep gradients of density and velocities are due to the garnet-ilmenite transition between 610 and 690 km. This is followed by abrupt jumps at about 690 km for the post-spinel transition, and steep gradients between 700 and 740 km that correspond to the ilmenite-perovskite transition. In the latter profile along the low-temperature geotherm, density and velocity increases for garnet-ilmenite and ilmenite-perovskite transitions are similar in magnitude to those for the post-spinel transition. The likely presence of ilmenite in cooler regions of subducting slabs is suggested by the fact that the calculated velocity profiles along the low-temperature geotherm are compatible with recent seismic observations indicating three discontinuities or steep velocity gradients at around 600-750 km depth in the regions of subducting slabs.     

Approximate solutions for the formation of the lithosphere

The lithosphere is interpreted as a thermal boundary layer. Approximate solutions of the boundary layer cooling problem are developed which include mantle radioactivity, partial melt in the asthenosphere, a temperature gradient in the asthenosphere, and a non-zero lithospheric thickness at the ridge crests. The cooling history of oceanic lithosphere is found to be remarkably insensitive to assumptions about the amount of radioactivity in the upper mantle and the extent of melting in the asthenosphere. Determinations of the thickness of oceanic lithosphere and the depths of oceans as a function of age are in excellent agreement with boundary layer predictions which include a heat flux from the asthenosphere. However, the determinations do not resolve how much of the total asthenospheric heat flux might be caused by a temperature gradient in the asthenosphere. Simple thermal arguments indicate that the initial lithospheric thickness, L₀, at ridge crests should depend on the local half-spreading rate, V, as L₀ = 3 km/V(cm/year).     

Softening of the subcontinental lithospheric mantle by asthenosphere melts and the continental extension/oceanic spreading transition

The majority of ophiolitic peridotites in the Alpine–Apennine system show evidence of extensive interaction between subcontinental lithospheric mantle and fractional melts of asthenospheric origin. This interaction led to petrological, structural, and geochemical changes in the lithospheric mantle, and was accompanied by a temperature increase to near-asthenospheric values, resulting in the thermomechanical erosion of the lithosphere. We term the parts of mantle lithosphere thus affected the asthenospherized lithospheric mantle or ALM.The thermal and rheological consequences of thermomechanical erosion are explored by modelling the temperature and rheological properties of the thinned lithosphere as a function of thickness of ALM and time since asthenospherization (i.e., since the beginning of thermal relaxation). Results are given both in terms of rheological profiles (strength envelopes) and total lithospheric strength (TLS) for different lower crustal rheologies. The TLS decreases as a consequence of thermomechanical erosion. This decrease is a non-linear function of the thickness of the ALM. While practically negligible if less than 50% of lithospheric mantle is affected, it becomes significant (up to almost one order of magnitude) if thermomechanical erosion approaches the Moho. The maximum decrease in TLS is achieved within a short time span (1–2 Ma) after the end of the heating episode.As a working hypothesis, we propose that thermomechanical erosion of the lithospheric mantle, related to lithosphere/asthenospheric melts interaction, can be an important factor in a geologically rapid decrease in TLS. This softening could lead to whole lithospheric failure and consequently to a transition from continental extension to oceanic spreading.     

地幔热导率的选取对动力学数值模拟的影响——以岩石圈张裂过程为例

目前存在有多种地幔热导率模型,不同模型在数值和随温压变化的特征上有明显的差异.为探究不同热导率模型对动力学数值模拟结果的影响,本文对不同模型下的岩石圈张裂过程进行模拟研究,探讨地幔热导率对岩石圈热传输、变形和熔融过程的影响及其作用机理.结果显示,不同热导率模型下,岩石圈的变形和熔融特征表现出明显差异.高热导率模型下,岩石圈破裂较晚,形成陆缘较为宽阔,地壳熔融强烈而地幔熔融较弱;低热导率模型下,岩石圈破裂较早,形成陆缘较为狭窄,地幔熔融强烈而地壳熔融较弱.这种差异源于不同地幔热导率下岩石圈和地幔热状态的变化及相应力学性质的改变.高热导率下,热传导的增温效应显著,岩石圈呈现较热的状态,其强度整体较低,壳幔耦合减弱;而低热导率下,热对流的增温效应显著,岩石圈呈较冷的状态,其强度整体较高,壳幔耦合增强.基于模拟结果,本文认为地幔热导率的选取对动力学模拟的结果有着较为显著的影响,相对于随温压的变化,热导率数值的差异对动力学数值模拟的结果影响更大,尤其是对于地幔熔融过程的影响.     

The system enstatite-pyrope at high pressures and temperatures and the mineralogy of the earth's mantle

In a diamond-anvil pressure cell coupled with laser heating, the system enstatite (MgSiO₃)-pyrope (3 MgSiO₃ · Al₂O₃) has been studied in the pressure region between about 100 and 300 kbar at about 1000°C using glass starting materials. The high-pressure phase behavior of the intermediate compositions of the system contrasts greatly with that of the two end-members. Differences between MgSiO₃ and 95% MgSiO₃ · 5% Al₂O₃ are especially remarkable. The phase assemblages β-Mg₂SiO₄ + stishovite and γ-Mg₂SiO₄ (spinel) + stishovite displayed by MgSiO₃ were not observed in 95% MgSiO₃ · 5% Al₂O₃, and the garnet phase, which was observed in 95% MgSiO₃ · 5% Al₂O₃ at high pressure, was not detected in MgSiO₃. These results suggest that the high-pressure phase transformations found in pure MgSiO₃ would be inhibited under mantle conditions by the presence even of small amounts of Al₂O₃ (?4% by weight). On the other hand, pyrope displays a wide stability field, finally transforming at 240–250 kbar directly to an ilmenite-type modification of the same stoichiometry. The two-phase region, within which orthopyroxene and garnet solid solutions coexist, is very broad. The structure of the earth's mantle is discussed in terms of the phase transformations to be expected in a simple mixture of 90% MgSiO₃ · 10% Al₂O₃ and Mg₂SiO₄. The seismic discontinuity at a depth of 400 km in the earth's mantle is probably due entirely to the olivine → β-phase transition in Mg₂SiO₄, with the progressive solution of pyroxene in garnet (displayed in 90% MgSiO₃ · 10% Al₂O₃) occurring at shallower depths. The inferred discontinuity at 650 km is due to the combination of the phase changes spinel → perovskite + rocksalt in Mg₂SiO₄ and garnet → ilmenite in 90% MgSiO₃ · 10% Al₂O₃. The 650-km discontinuity is thus characterized by an increase in the primary coordination of silicon from 4 to 6. A further discontinuity in the density and seismic wave velocities at greater depth associated with the ilmenite-perovskite phase transformation in 90% MgSiO₃ · 10% Al₂O₃ is expected.