On the importance of advection in determining the local isotopic composition of the mantle

The equations for isotopic evolution in a deforming medium are derived and used to show that the local isotopic composition of the mantle depends primarily on the advection of heterogeneities by the flow field. Various examples of the role of advection are given using a two-scale model of mantle convection. The effect of small-scale flows, which have dimensions smaller than the plates themselves, is to convert any initially localized heterogeneity into thin spiral sheets. The isotopic properties of erupted lavas will depend on the relative size of the zone of partial melt to the spacing between such sheets. An idealized model for⁸⁷Sr/⁸⁶Sr variations from Pacific islands and the covariation of¹⁴³Nd/¹⁴⁴Nd-⁸⁷Sr/⁸⁶Sr is based on the combined effect of isotopic fluxes at subduction zones and advection by the return flow under the moving lithospheric plates.     

Tracking the lithium isotopic evolution of the mantle using carbonatites

Carbonatites are mantle-derived, intraplate magmas that provide a means of documenting isotopic variations of the Earth's mantle through time. To investigate the secular Li isotopic evolution of the mantle and to test whether Li isotopes document systematic recycling of material processed at or near the Earth's surface into the mantle, we analyzed the Li isotopic compositions of carbonatites and spatially associated mafic silicate rocks. The Li isotopic compositions of Archean (2.7 Ga) to Recent carbonatites (δ⁷Li = 4.1 ± 1.3 (n = 23, 1σ)) overlap the range typical for modern mantle-derived rocks, and do not change with time, despite ongoing crustal recycling. Thus, the average Li isotopic composition of recycled crustal components has not deviated greatly from the mantle value (~ + 4) and/or Li diffusion is sufficiently fast to attenuate significant heterogeneities over timescales of 10⁸ years. Modeling of Li diffusion at mantle temperatures suggests that limited δ⁷Li variation in the mantle through time reflects the more effective homogenization of Li in the mantle compared to radiogenic isotope systems. The real (but limited) variations in δ⁷Li that exist in modern mantle-derived magmas as well as carbonatites studied here may reflect isotopic fractionation associated with shallow-level processes, such as crustal assimilation and diffusive isotopic fractionation in magmatic systems, with some of the scatter possibly related to low-temperature alteration.     

Statistical analysis of isotopic ratios in MORB: the mantle blob cluster model and the convective regime of the mantle

A statistical examination of isotopic distributions for MORB from various ocean ridges leads to the “blob cluster model”, in which the oceanic crust accreting at ridges results from the mixing of two components within the ascending mantle. These are (1) upper mantle material and (2) discrete rising blobs of more radiogenic material. The blobs are fractionated to a variable degree and are distributed in the upper mantle circulation in a manner that is related to the spreading rate.(1) Themean values of the isotopic distributions allow us to calculate the probabilities of the two types of material within the mantle. The results show that theproportion of asthenospheric material in the mixtureincreases with the spreading rate, in agreement with the hypothesis of blob dilution within the upper mantle convection.Mass fluxes can be estimated for the rising blobs from these probabilities, which depend on the respective concentrations in the sources of the two types of material. If the blobs originate in the lower mantle, this flux estimation would suggest that a significant part of the lower mantle has been injected into the upper mantle during earth history.(2) Thestandard deviations of the distributions depend on the “efficiency” of the mixing process:the more imbricated are the asthenospheric and blob materials in the mixture,the smaller is theisotopic spread. This efficiency parameter is shown to increase with the spreading rate, as already suggested by previous comparisons between the East Pacific Rise and the Mid-Atlantic Ridge. Moreover, this feature may also be correlated with other data such as ridge bathymetric variations.     

Episodic mantle differentiation: Nd and Sr isotopic evidence

Clinopyroxenes separated from two hydrous and four anhydrous ultramafic nodules, selected from a suit of xenoliths from Dreiser Weiher (DW), West Germany, have been studied for Nd and Sr isotopic composition. Nd exhibits a range of ε_JUV(T) from 0 to +12.4 and ⁸⁷Sr/⁸⁶Sr varies between 0.70185 and 0.70400. T_ICE model ages for anhydrous nodules indicate that the mantle underlying DW was originally depleted ?2 AE ago. Correlation of ¹⁴³Nd/¹⁴⁴Nd with Sm/Nd in this group of samples suggests that a second partial melting event occurred about 560 m.y. ago resulting in LREE enrichment of at least part of the anhydrous mantle. During a later episode, probably contemporaneous with the eruption of the host basalt in Quaternary times, most of the spinel peridotitic upper mantle below DW was modified. This metasomatism led to hydration and incompatible element enrichment of originally anhydrous mantle. The isotopic data for the anhydrous nodules again demonstrate that oceanic-type mantle underlies at least some continental areas. It is apparent that the separation of subcontinental mantle regions from an initially chondritic reservoir may occur in several discrete episodes. However, differing histories of depletion and/or enrichment will produce isotopically distinct mantle reservoirs. Therefore, basalts extracted from these mantle reservoirs will scatter about an average Nd-Sr trend line reflecting the nature of the differentiation in their source regions.     

Osmium-strontium-neodymium-lead isotopic covariations in mid-ocean ridge basalt glasses and the heterogeneity of the upper mantle

Osmium, strontium, neodymium, and lead isotopic data have been obtained for 30 hand picked samples of basaltic glass from the Pacific, Atlantic and Indian mid-oceanic ridges. Large variations in Os isotopic ratios exist in the glasses, from abyssal peridotite-like values to radiogenic compositions similar to oceanic island basalts (¹⁸⁷Os/¹⁸⁶Os and ¹⁸⁷Os/¹⁸⁸Os ratios range from 1.06 to 1.36 and from 0.128 to 0.163, respectively). Os isotopic and elemental data suggest the existence of mixing correlations. This relationship might be ascribed to secondary contamination processes; however, such a hypothesis cannot account for the negative correlation observed between Os and Nd isotopes and the existence of complementary covariations between Os and SrPb isotopes. In this case, OsSrNdPb isotopic variations are unrelated to late post-eruption or shallow level contamination. These relationships provide strong evidence that the Os isotopic composition of the samples are derived from the mantle and thus implies a global chemical heterogeneity of the oceanic upper mantle. The results are consistent with the presence of recycled oceanic crust in the mantle sources of mid-ocean ridge basalts, and indicate that the unique composition of the upper mantle below the Indian ocean results from its contamination by a mantle component characterized by radiogenic Os and particularly unradiogenic Nd and Pb isotopic compositions.     

The NdSr isotopic correlation in mantle materials and geodynamic consequences

The neodymium-strontium isotopic correlation observed in most of the Earth mantle materials is evaluated by means of direct modelling. Several geochemical models are quantitatively developed to explain the observations. The main results of this modelling are that such a correlation is not geochemically trivial and that it corresponds to specific conditions in chemical fractionation. These specific conditions seem to be satisfied by solid-liquid partitioning in magmatic conditions. The discussion of the experimental data supports a continuous convecting-magmatic fractionation model for a large proportion of the mantle.     

Lead concentration and isotopic composition in five peridotite inclusions of probable mantle origin

The lead content of five whole-rock peridotite inclusions (four lherzolites and one harzburgite) in alkali basalt ranges from 82 to 570 ppb (parts per billion). Approximately 30–60 ppb of this amount can be accounted for by analyzed major silicate minerals (olivine ≤ 10 ppb; enstatite 5–28 ppb; chrome diopside ～400 ppb). Through a series of acid leaching experiments, the remainder of the lead is shown to be quite labile and to reside in either glassy or microcrystalline veinlets or accessory mineral phases, such as apatite and mica. The lead isotopic composition of the peridotites (²⁰⁶Pb/²⁰⁴Pb= 18.01–18.90;²⁰⁷Pb/²⁰⁴Pb= 15.52–15.61;²⁰⁸Pb/²⁰⁴Pb= 37.80–38.86) lies within the range of values defined by many modern volcanic rocks and, in particular, is essentially coextensive with the abyssal tholeiite field. In all but one instance, isotopic differences were found between the peridotite and its host alkali basalt. Two of the peridotites clearly demonstrated internal isotopic heterogeneity between leachable and residual fractions that could not simply be due to contamination by the host basalt. However, there is no evidence that these ultramafic rocks form some layer in the mantle with isotopic characteristics fundamentally different from those of the magma sources of volcanic rocks.     

Sr isotopic composition of Afar volcanics and its implication for mantle evolution

Investigations of Rb-Sr systematics of basalts from the Afar depression (Ethiopia) indicate the presence of a heterogeneous mantle source region. The Sr isotopic compositions of the basalts from the Afar axial and transverse ranges identify source regions which are enriched in LIL elements and radiogenic Sr (axial ranges) and others which are relatively depleted (transverse ranges). Sr isotopic composition of basalts from the Red Sea, Gulf of Aden and Gulf of Tadjoura, which range from 0.70300 to 0.70340 are also reported and compared with the more radiogenic Afar region, which is characterized by⁸⁷Sr/⁸⁶Sr ranging from 0.70328 to 0.70410.Available geochemical and isotopic data suggest that a relation exists between magma composition and the advancement of the rifting process through progressive lithosphere attenuation leading to continental break-up. However, the petrogenetic process is not simple and probably implies a vertically zoned mantle beneath the Afar region. Sr isotopic evidence suggests that the vertically zoned mantle is more radiogenic and enriched in LIL elements in its upper part.     

Implications of correlated Nd and Sr isotopic variations for the chemical evolution of the crust and mantle

Published data showing a linear correlation between initial Nd and Sr isotope compositions in young basalts indicate the existence of a spectrum of isotopically distinct reservoirs in the mantle which represent either (1) mixtures of two homogeneous endmember reservoirs, one of which may be undifferentiated material or (2) fractionated reservoirs all derived from a homogeneous initial reservoir with the same ratio of enrichment factors for Sm/Nd and Rb/Sr. The slope of the correlation, which can be described approximately by (⁸⁷Sr/⁸⁶Sr) = ?3.74114 (¹⁴³Nd/¹⁴⁴Nd) + 2.61935orε_Nd = ?2.7 ε_Sr, places constraints on the origin of these reservoirs and hence on the chemical evolution of the crust-mantle system. The reservoirs could be residual regions of the mantle left after ancient partial melting events. If so, the requirement of constant relative fractionation of Sm/Nd and Rb/Sr in refractory residues is a strong constraint on partial melting models. Calculations suggest that batch melting models are more compatible with this constraint than are fractional melting models, but models incorporating currently accepted distribution coefficients and residual phase assemblages cannot reproduce the observed isotope effects except under highly specific conditions. The slope of the correlation is not consistent with the hypotheses that chemical structure in the mantle is due to accretional heterogeneity or variable loss of elements to the core. If the mantle reservoirs are complementary in composition to the continental crust, and if the crust + mantle has ε_Nd = 0andε_Sr = 0 and chondritic Sr/Nd, then Rb/Sr in the crust is calculated to be less than 0.10, suggesting that the crust may be more mafic in composition and contain a smaller proportion of the earth's Rb and heat-producing elements than previously estimated.     

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.     

An approach to the assessment of local trace metal pollution in the Mediterranean marine atmosphere

Assessing the extent of local trace metal pollution in marine atmospheric particulates presents a number of difficulties, and one potential approach to the problem has been evaluated by reference to a series of aerosols from the lower troposphere over the Eastern Mediterranean. This region, which is less remote from the continents than open-ocean areas, has two significantly different atmospheric particulate catchments: being bordered in the north by nations having industrialized, semi-industrialized and rural economies, and in the south by the North African desert belt. The elemental chemistry of the particulates is illustrated in terms of the distributions of Fe and Pb, which are presented in the form of enrichment factor diagrams, and is shown to be essentially controlled by the dilution of a ‘European’ background material, common to ‘remote’ regions of the latitudinal belt, with crust-derived, desert components. Major local perturbations in the dilution relationship can occur when either crust-derived or pollutant components strongly dominate the total particulate population. The extent of these local perturbations can be tentatively assessed using appropriate enrichment factor diagrams, which offer a framework within which to interpret the elemental chemistry of marine atmospheric particulates providing sufficient data is available for their construction. By using such an EF diagram it is suggested that the Eastern Mediterranean atmospheric particulates have not suffered local pollution on a gross scale with respect to Pb.     

Atmospheric vapour and precipitation are not in isotopic equilibrium in a continental mountain environment

Isotopic exchange with atmospheric vapour can strongly influence the isotopic values of evaporating surface water bodies (e.g., lakes), influencing our understanding of hydrological processes across aquatic and terrestrial environments. Rather than measure the isotopic values of the atmosphere directly, it is much more common to estimate values by assuming equilibrium with local precipitation. This assumption may introduce large errors, thereby biasing hydrological inferences and understanding. The pattern and magnitude of this error has been quantified only in a few circumstances. We compared observations of vapour and precipitation isotope values over a four-year record collected in a montane environment in the central Rocky Mountains of North America. We further investigated factors and conditions promoting disequilibrium. Scenario comparisons assessed the impact of theoretical and methodological elements on the accuracy of the equilibrium assumption. We found that the equilibrium assumption was not well supported by direct and continuous observations of vapour isotopes using tower-based laser isotope spectroscopy, particularly during the summer months. Across all scenarios, errors associated with the equilibrium assumption were high, credibly ranging from 14 to 154 ‰ for δ²H and 1.5 to 16.3 ‰ for δ¹⁸O. Environmental covariates (e.g., vapour pressure deficit, air pressure) helped explain some of the apparent disequilibrium. Although the equilibrium assumption for estimating atmospheric vapour isotope values may not be applicable in a continental montane environment, our findings highlight opportunities for using direct vapour isotope measurements to better understand vapour sources, air mass mixing, and phase changes over complex mountainous terrain, which in turn may better constrain regional- to global-scale hydrological process estimates, such as evapotranspiration rates and the water budgets of mountain lakes.     

基于优化参数的地震P、S波震相到时自动拾取及质量评估

准确拾取P、S波震相到时是深入开展地震波研究工作的基础,本文改进了自动拾取参数优化函数算法和质量评估方案,引入了拾取到时优化方案,使用基于参数优化的频带-带宽拾取算法、AICD拾取算法和峰度拾取算法对腾冲地区7个宽频带地震台站记录的地震资料开展了地震P、S波到时自动拾取,对拾取结果进行了优化和质量判定.结果表明：经参数优化、拾取优化后,采用3种方法自动拾取的P、S波到时与人工拾取到时的时差在0.1 s内的记录占比分别达到74.66%、70.98%.这些参数值均优于算法改进前的同类参数,证明了优化方法的可靠性.

     

An integrated assessment of global and regional water demands for electricity generation to 2095

Electric power plants account for approximately half the global industrial water withdrawal. Although continued electric-sector expansion is probable, significant variations in water intensity by electricity technology and cooling system type make its effects on water demands uncertain. Using GCAM, an integrated assessment model of energy, agriculture, and climate change, we establish lower-, median-, and upper-bound estimates for current electric-sector water withdrawals and consumption in 14 geopolitical regions, and compare them with available estimates. We then explore water use for electricity to 2095, focusing on uncertainties in water withdrawal and consumption intensities, power plant cooling system changes, and adoption rates of water-saving technologies. Results reveal a probable decrease in the water withdrawal intensity with capital stock turnover, but a corresponding increase in consumptive use, for which technologies under development may compensate. At a regional scale, water use varies significantly based on the existing capital stock and its evolution over the century.     

A Nd isotopic study of the Kerguelen Islands: Inferences on enriched oceanic mantle sources

The origin of the highly differentiated igneous rocks of the Kerguelen Islands and the nature of their source regions have been investigated by a Nd isotopic study. The Nd isotopic compositions of syenites and granites are identical to those of gabbros and basalts and indicate a common source. The isotopic data preclude the involvement ofold continental crustal material in the genesis of these granitic and alkalic rocks. The data from the Kerguelen samples greatly extend the Nd-Sr isotopic correlation observed for uncontaminated basalts from the oceanic mantle. The large Nd isotopic variations in the Kerguelen samples could be explained by mixing of deep mantle material brought up by a plume and the upper oceanic mantle or by heterogeneities in the lower mantle. An important finding of this study is that there are enriched mantle sources under the oceanic regions. These enriched sources may be ancient in age and are compatible with the 2-b.y. age inferred from the Pb isotope data of these samples. Earth models in future must incorporate this feature of the oceanic mantle in a consideration of mantle-crust evolutionary relationships.     

Strontium and neodymium isotopic evidence for the heterogeneous nature and development of the mantle beneath Afar (Ethiopia)

Neodymium isotope and REE analyses of recent volcanic rocks and spinel lherzolite nodules from the Afar area are reported. The¹⁴³Nd/¹⁴⁴Nd ratios of the volcanic rocks range from 0.51286 to 0.51304, similar to the range recorded from Iceland. However, the⁸⁷Sr/⁸⁶Sr ratios display a distinctly greater range (0.70328–0.70410) than those reported from the primitive rocks of Iceland. Whole rock samples and mineral separates from the spinel lherzolite nodules exhibit uniform¹⁴³Nd/¹⁴⁴Nd ratios (ca. 0.5129) but varied⁸⁷Sr/⁸⁶Sr ratios in the range 0.70427–0.70528.The SrNd isotope variations suggest that the volcanic rocks may have been produced by mixing between two reservoirs with distinct isotopic compositions. Two possible magma reservoirs in this area are the source which produced the “MORB-type” volcanics in the Red Sea and Gulf of Aden and the anomalous source represented by the nodule suite. The isotopic composition of the volcanics is compatible with mixing between these two reservoirs.It is shown that the anomalous source with a high⁸⁷Sr/⁸⁶Sr ratio cannot have been produced by simple processes of partial melting and mixing within normal mantle. Instead the high⁸⁷Sr/⁸⁶Sr is equated with a fluid phase. A primitive cognate fluid, subducted seawater or altered oceanic lithosphere may have been responsible for the generation of the source with a high⁸⁷Sr/⁸⁶Sr ratio.     

Core formation and chemical equilibrium in the Earth—II: Chemical consequences for the mantle and core

We present here a new model of core formation which is based on the current understanding of planetary accretion and discuss its implications for the chemistry of the Earth's mantle and core. Formation of the Earth by hierarchical accretion of progressively larger bodies on a time scale much longer than that of solid body differentiation in the nebula indicates that a significant fraction of metal in the core could be inherited from preterrestrially differentiated planetesimals. An analysis of the segregation of this iron to form the core suggests that most of the metal settles to the core without interaction with silicates; only a small fraction of the metal chemically equilibrates at high temperatures and pressures with the silicates. The siderophile element abundances in the mantle are considered to be a consequence of a two-step equilibration with iron, once preterrestrially in the planetesimals at low temperatures and pressures, and later in the Earth at high temperatures and pressures. The highly siderophile elements such as Re, Au and the platinum group elements in the mantle are essentially excluded from silicates from the preterrestrial equilibration. We attribute the abundances of these elements in the mantle to the later equilibration in the Earth at substantially reduced metal-silicate partition coefficients (D^met/sil), for which there is a considerable experimental evidence now. Mass balance considerations constrain the fraction of core metal involved in such an equilibration at approximately 0.3 – 0.5%. The model accounts for the levels and the near-chondritic ratios of the highly siderophile elements in the mantle. The mantle abundances of the less siderophile elements are largely determined by preterrestrial metal-silicate equilibrium and are not significantly affected by the second equilibration. The extreme depletion of sulfur and the lack of silicate melt-sulfide signature in the noble metal abundances in the mantle are natural consequences of this mode of core formation. Sulfur was added to the magma ocean during the high-T, high-P equilibration in the Earth, not extracted from it by sulfide segregation to the core. Except for Ni and Co, the overall siderophile abundances of the mantle can be well matched in this two-step equilibration model.

The mantle characteristics of Ni and Co are unique to the Earth and hence suggest a terrestrial process as the likely cause. One such process is the flotation and addition of olivine to the primitive upper mantle. In our model of core formation, neither the elemental and isotopic data of Re---Os, nor the low sulfur content of the mantle remains as an objection to the existence of a magma ocean and olivine flotation.

The small fraction of core metal that equilibrates with silicates at high T and P suggests that the light elements O, Si or H are unimportant in the core, leaving S (and possibly C) as prime candidates. Sulfur, as FeS associated with incoming iron metal, is directly sequestered to the core along with the bulk of the iron metal. It appears unlikely that other light elements can be added to the core after its formation. U and Th are excluded from the core but the model allows for entry of some K; however, the extent to which K serves as a heat source in the core remains uncertain.

The model is testable in two ways. One is by investigation of the metal-silicate partitioning at high temperatures and pressures under magma ocean conditions to determine if the (D^met/sil) values are lowered to the levels required in the model. The other is by experiments to determine if a solvus closure between metal and silicate liquids occurs at high temperatures relevant to a magma ocean.     

The southeast Australian lithospheric mantle: isotopic and geochemical constraints on its growth and evolution

Trace elements and isotopic compositions of whole rocks and mineral separates are reported for 15 spinel-bearing harzburgite and lherzolite xenoliths from southeastern Australia. These samples have an exceedingly large range in isotopic compositions, with⁸⁷Sr/⁸⁶Sr ranging from 0.70248 to 0.70834 and ε_Nd values ranging from +12.7 to −6.3. This range in isotopic compositions can be found in xenoliths from a single locality. The isotopic compositions of clinopyroxene separates and their whole rocks were found to be different in some xenoliths. Samples containing small glass pockets, which replace pre-existing hydrous minerals, generally show only small differences in isotopic composition between clinopyroxene and whole rock. In a modally metasomatized peridotite, significant differences in the Sr and Nd isotopic compositions of a coexisting phlogopite-clinopyroxene pair are present. Coexisting clinopyroxenes and orthopyroxenes from an anhydrous lherzolite have Sr isotopic compositions that are significantly different (0.70248 versus 0.70314), and yield an apparent age of 625 Ma, similar to that found previously by Dasch and Green [1]. However, the Nd isotopic compositions of the clinopyroxene and orthopyroxene are identical indicating recent (within 40 Ma) re-equilibration of Nd.Sr and Nd concentrations in the whole rocks and clinopyroxenes show an excellent positive correlation, and have an average Sr/Nd ratio of 15. This ratio is similar to the primitive mantle value, as well as that found in primitive MORBs and OIBs, but is much lower than that measured in island arc basalts and what might be predicted for a subduction zone-derived fluid. This indicates that a significant proportion of the Sr and Nd in these peridotites is introduced as a basaltic melt with intraplate chemical characteristics.The isotopic compositions of the peridotites reflect long-term, small-scale heterogeneities in the continental lithospheric mantle, and are in marked contrast to the near uniform isotopic compositions of the host alkali basalts (⁸⁷Sr/⁸⁶Sr= 0.7038–0.7041andε_Nd = +3.6 to +2.9). A minimum of three evolutionary stages are identified in the growth of the continental lithospheric mantle: an early basalt depletion event, recording the initial development and stabilization of the lithospheric mantle, followed by at least two enrichment episodes. These observations are consistent with continental lithospheric mantle growth involving the underplating of refractory peridotite diapirs.     

Phase equilibrium modeling of zircon stability in mantle peridotite: Implication for crust-mantle interaction

Zircon is a common accessory mineral in various rocks,especially in the crustal ones.It is the best mineral for U-Pb dating.Meanwhile,trace elements and isotopes of the mineral can also provide much information concerning the formation and evolution of rocks.There are a growing number of reports of zircon existing in mantle peridotite.However,it is generally considered that zircon is unlikely crystallized in ultrabasic rocks due to SiO₂-unsaturation.In this paper,the SiO₂ activity and zircon/baddeleyite transition curve at different conditions were calculated through thermodynamic phase equilibrium modeling,to reveal the main factors affecting the SiO₂ activity and the stability of zircon/baddeleyite in ultrabasic and basic rocks,especially in mantle peridotite.These results provide a thermodynamic basis for interpreting the genesis and significance of zircon in mantle rocks.That is,the SiO₂ activity is mainly controlled by stable mineral assemblages and temperature-pressure conditions.The orthopyroxene+olivine assemblage in peridotite as an effective buffer restricts the SiO₂ activity in a relatively high range with a small variation.The upper temperature limit of zircon can reach more than 1500℃ with this mineral assemblage.During the low-temperature serpentinization of peridotite,the replacement of olivine and pyroxene by serpentine can result in a significant decrease of SiO₂ activity,and baddeleyite can be stabilized at<530℃ and<2.7 GPa.When peridotite is strongly metasomatized by the SiO₂-bearing fluid,the addition of SiO₂ can increase its activity and make zircon stable at low temperatures.The SiO₂ activity in ultrabasic-basic rocks is not only positively correlated with the SiO₂ content but also negatively correlated with the Ca and Na contents of rocks.This is because Ca and Na preferentially combine with Si and Al to form Si-rich minerals,such as clinopyroxene and feldspar.This process will consume excessive SiO₂,decreasing the SiO₂ activity.This may be the reason why zircon can be found in ultrabasic rocks,while baddeleyite can exist in some basic and alkaline rocks.The thermodynamic modeling can also reasonably explain the mutual transformation between zircon and baddeleyite in ultrabasic-basic rocks.Our results indicate that zircon can exist stably in mantle peridotite in a wide range of temperature-pressure conditions and its formation is related to melt/fluid metasomatism.That is,the presence of zircon in mantle peridotite is an important information carrier of crust-mantle interaction for deep material cycling.