Implications of melting for stabilisation of the lithosphere and heat loss in the Archaean

The increased depth and volume of melting induced in a higher temperature Archaean mantle controls the stability of the lithosphere, heat loss rates and the thickness of the oceanic crust. The relationship between density distributions in oceanic lithosphere and the depth of melting at spreading centres is investigated by calculating the mineral proportions and densities of residual mantle depleted by extraction of melt fractions. The density changes related to compositional gradients are comparable to those produced by thermal effects for lithosphere formed from a mantle which is 200°C or more hotter than modern upper mantle. If Archaean continental crust formed initially above oceanic lithosphere, the compositional density gradients may be sufficient to preserve a thick Archaean continental lithosphere within which the Archaean age diamonds are preserved. The amount of heat advected by melts at mid-ocean ridges today is small but heat advected by melting becomes proportionally more important as higher mantle temperatures lead to a greater volume of melt and as the rate of production of oceanic plates increases. Archaean tectonics could have been dominated by spreading rates 2–3 times greater than now and with mantle temperatures between ca. 1600°C and 1800°C at the depth of the solidus. Mid-ocean ridge melting would produce a relatively thick but light refractory lithosphere on which continents could form, protected from copious volcanism and high mantle temperatures.     

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

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

Thermal evolution of the Earth: Secular changes and fluctuations of plate characteristics

The average secular cooling rate of the Earth can be deduced from compositional variations of mantle melts through time and from rheological conditions at the onset of sub-solidus convection at the end of the initial magma ocean phase. The constraint that this places on the characteristics of mantle convection in the past are investigated using the global heat balance equation and a simple parameterization for the heat loss of the Earth. All heat loss parameterization schemes depend on a closure equation for the maximum age of oceanic plates. We use a scheme that accounts for the present-day distribution of heat flux at Earth's surface and that does not depend on any assumption about the dynamics of convection with rigid plates, which remain poorly understood. We show that heat supply to the base of continents and transient continental thermal regimes cannot be ignored. We find that the maximum sea floor age has not changed by large amounts over the last 3 Ga. Calculations lead to a maximum temperature at an age of about 3 Ga and cannot be extrapolated further back in time. By construction, these calculations are based on the present-day tectonic regime characterized by the subduction of large oceanic plates and hence indicate that this regime did not prevail until an age of about 3 Ga. According to this interpretation, the onset of rapid continental growth occurred when the current plate regime became stable.     

Core cooling by subsolidus mantle convection

Although vigorous mantle convection early in the thermal history of the Earth is shown to be capable of removing several times the latent heat content of the core, we are able to construct a thermal evolution model of the Earth in which the core does not solidify. The large amount of energy removed from the model Earth's core by mantle convection is supplied by the internal energy of the core which is assumed to cool from an initial high temperature given by the silicate melting temperature at the core-mantle boundary. For the smaller terrestrial planets, the iron and silicate melting temperatures at the core-mantle boundaries are more comparable than for the Earth, and the cores of these planets may not possess enough internal energy to prevent core solidification by mantle convection. Our models incorporate temperature-dependent mantle viscosity and radiogenic heat sources in the mantle. The Earth models are constrained by the present surface heat flux and mantle viscosity. Internal heat sources produce only about 55% of the Earth model's present surface heat flow.     

大陆岩石圈、地幔底部异常体与地幔对流相互作用的数值模拟

在地球表层存在着占地表面积约30%的具有低固有密度、高黏度的大陆岩石圈.由于其特殊的物理化学性质,大陆岩石圈通常不直接参与下方的地幔对流,但其与地幔对流格局有着重要的相互影响.大量研究显示,在中太平洋和非洲的下地幔底部,存在着两块占核幔边界（CMB）面积约20%的高密度热化学异常体（由于其剪切波速度较低,常称作低剪切波速度省（LSVPs））.LSVPs的演化既受地幔对流的影响,同时也影响地幔物质运动的格局和动力学过程.本文系统研究了存在大陆岩石圈,下地幔LSVPs的地幔对流模型.模拟结果显示：（1）当大陆体积较小时,其边缘常伴随着俯冲,大陆区域地幔常处于下涌状态,其上地幔温度较低,大陆岩石圈在水平方向处于压应力状态.随着大陆体积的增大,大陆边缘的俯冲逐渐减弱,大陆区域地幔由下涌转为上涌,其上地幔温度较高,大陆岩石圈水平方向处于拉应力状态.（2） 岩石圈与软流圈边界（LAB）在大陆下方较深,温度较低;在海洋区域较浅,温度较高.随着大陆体积的增大,陆洋之间LAB深度、温度的差异逐渐减小.（3）大陆区域地幔底部LSVPs物质的丰度与大陆的体积呈正相关.当大陆体积较小时,大陆下方的LSVPs丰度比海洋区域少.随着大陆体积的增大,大陆下方LSVPs的丰度逐渐增大.（4）海洋地区地表热流高,且随时间波动大,大陆地区地表热流低,随时间波动较小;LSVPs区域的核幔边界热流低.     

Continental growth and volatile exchange during Earth''s evolution

We investigate the thermal and degassing history of the Earth with the help of a parameterized mantle convection model including the volatile exchange between mantle and surface reservoirs. The weakening of mantle silicates by dissolved volatiles is described by a functional relationship between creep rate and water fugacity. We use flow law parameters of diffusion creep in olivine under dry and wet conditions. The mantle degassing rate is considered as directly proportional to the seafloor spreading rate, which is also dependent on the mantle heat flow and the continental area. To calculate the spreading rate, we assume three different continental growth models: constant growth, delayed growth, and the one proposed by Reymer and Schubert (1984, Tectonics, 3: 63–77). The rate of regassing also depends on the seafloor spreading rate, as well as on other factors. Both mechanisms (degassing and regassing) are coupled self-consistently with the help of a parameterized convection model under implementation of a temperature and volatile-content dependent mantle viscosity. We calculate time series for the Earth's evolution over 4.6 Gyr for the average mantle temperature, the mantle heat flow, the mantle viscosity, the Rayleigh number, the Urey ratio, the volatile loss, and the seafloor spreading rate. In those numerical simulations with continental growth from the beginning and a high initial average mantle temperature water is outgassed rapidly. In the delayed continental growth model there is a very early outgassing event and the delayed continental growth has no remarkable influence on the thermal and outgassing history. A similar situation is found for the linear continental growth model but not for the Reymer and Schubert (1984) model.     

The structure of convection in the spherical mantle with internal heating

Thermal convection in the mantle is caused by the heat transported upwards from the core and by the heat produced by the internal radioactive sources. According to the data on the heat transfer by the mantle plumes and geochemical evidence, only 20% of the total heat of the Earth is supplied to the mantle from the core, whereas most of the heat is generated by the internal sources. Along with the models that correctly allow for the internal heat sources, there are also many publications (including monographs) on the models of mantle convection that completely ignore the internal heating or the heat flux from below. In this study, we analyze to what extent these approximations could be correct. The analytical distributions of temperature and heat flux in the case of internal heating without convection and the results of the numerical modeling for convection with different intensity are presented. It is shown that the structure of thermal convection is governed by the distribution of the heat flux in the mantle but not by the heat balance, as it is typically implicitly assumed in most works. Heat production by the internal sources causes the growth of the heat flux as a function of radius. However, in the spherical mantle of the Earth, the heat flux decreases with radius due to the geometry. It turned out that with the parameters of the present Earth, both these effects compensate each other to a considerable extent, and the resulting heat flux turns out to be nearly constant as a function of radius. Since the structure of the convective flows in the mantle is determined by the distributions of heat flux and total heat flux, in the Cartesian models of the mantle convection the effective contribution of internal heating is small, and ignoring the heat flux from the core significantly distorts the structure of the convective currents and temperature distributions in the mantle.     

Temperature and distribution of radioactive matter in the upper mantle

Temperatures of the Earth’s upper mantle, derived from Gutenberg’s seismic velocities, have been revised, using some recent determinations of the clastic constants of dunites. For depths greater than about 50 km the temperature obtained is sufficient to melt the basaltic fraction of a periodotitic mantle (if this low-melting fraction still exists at those depths). The distribution of the radioactive heat sources appears to be fairly uniform down to about 110 km; below this level the radioactive matter is probably absent. This is what would be expected for the upper mantle below the oceans; therefore, Gutenberg’s velocities seem to correspond to oceanic areas, rather than to the mantle below a continental shield.     

Chemical differentiation of the Earth: the relationship between mantle, continental crust, and oceanic crust

The average chemical compositions of the continental crust and the oceanic crust (represented by MORB), normalized to primitive mantle values and plotted as functions of the apparent bulk partition coefficient of each element, form surprisingly simple, complementary concentration patterns. In the continental crust, the maximum concentrations are on the order of 50 to 100 times the primitive-mantle values, and these are attained by the most highly incompatible elements Cs, Rb, Ba, and Th. In the average oceanic crust, the maximum concentrations are only about 10 times the primitive mantle values, and they are attained by the moderately incompatible elements Na, Ti, Zr, Hf, Y and the intermediate to heavy REE.This relationship is explained by a simple, two-stage model of extracting first continental and then oceanic crust from the initially primitive mantle. This model reproduces the characteristic concentration maximum in MORB. It yields quantitative constraints about the effective aggregate melt fractions extracted during both stages. These amount to about 1.5% for the continental crust and about 8–10% for the oceanic crust.The comparatively low degrees of melting inferred for average MORB are consistent with the correlation of Na₂O concentration with depth of extrusion [1], and with the normalized concentrations of Ca, Sc, and Al ( 3) in MORB, which are much lower than those of Zr, Hf, and the HREE ( 10). Ca, Al and Sc are compatible with clinopyroxene and are preferentially retained in the residual mantle by this mineral. This is possible only if the aggregate melt fraction is low enough for the clinopyroxene not to be consumed.A sequence of increasing compatibility of lithophile elements may be defined in two independent ways: (1) the order of decreasing normalized concentrations in the continental crust; or (2) by concentration correlations in oceanic basalts. The results are surprisingly similar except for Nb, Ta, and Pb, which yield inconsistent bulk partition coefficients as well as anomalous concentrations and standard deviations.The anomalies can be explained if Nb and Ta have relatively large partition coefficients during continental crust production and smaller coefficients during oceanic crust production. In contrast, Pb has a very small coefficient during continental crust production and a larger coefficient during oceanic crust production. This is the reason why these elements are useful in geochemical discrimination diagrams for distinguishing MORB and OIB on the one hand from island arc and most intracontinental volcanics on the other.The results are consistent with the crust-mantle differentiation model proposed previously [2]. Nb and Ta are preferentially retained and enriched in the residual mantle during formation of continental crust. After separation of the bulk of the continental crust, the residual portion of the mantle was rehomogenized, and the present-day internal heterogeneities between MORB and OIB sources were generated subsequently by processes involving only oceanic crust and mantle. During this second stage, Nb and Ta are highly incompatible, and their abundances are anomalously high in both OIB and MORB.The anomalous behavior of Pb causes the so-called “lead paradox”, namely the elevated U/Pb and Th/Pb ratios (inferred from Pb isotopes) in the present-day, depleted mantle, even though U and Th are more incompatible than Pb in oceanic basalts. This is explained if Pb is in fact more incompatible than U and Th during formation of the continental crust, and less incompatible than U and Th during formation of oceanic crust.     

海洋板块俯冲作用下上覆大陆岩石层减薄机制的动力学模拟

几乎所有大陆岩石层的减薄现象,可能都与海洋板块的俯冲作用相关,但是两者之间的内在联系迄今仍不十分明确,为此,我们设计了一系列包含洋-陆俯冲系统的二维数值模型,来探讨海洋板块的俯冲作用对上覆大陆岩石层变形行为的影响,尤其对大陆岩石层减薄效应的制约.模型结果表明,海洋板块俯冲过程中的地幔楔熔体对大陆岩石层地幔的热侵蚀以及由熔体上升所诱发的地幔局部对流的强烈扰动会导致上覆大陆岩石层的减薄效应.这种效应不仅表现在横向上的向陆内蔓延,还表现在垂向上的向浅部发展.且多类动力学参数都能制约大陆岩石层的减薄效应.具体地,随着汇聚速率和洋壳厚度的增加,上覆大陆岩石层在横向上的减薄范围越大,在垂向上的减薄程度也越深;而随着俯冲海洋板块年龄的增加,上覆大陆岩石层在横向上的减薄范围增大,但在垂向上的减薄程度会减小;随着上覆大陆岩石层厚度的增加,其横向减薄范围会减小,但在垂向上的减薄程度会加深.本文研究成果能为揭示华北克拉通减薄/破坏的动力学过程提供一定的理论参考依据.     

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.     

Isotopic inferences on the chemical structure of the mantle

Variations in the isotopic composition of rocks derived from the upper mantle can be used to infer the chemical history and structure of the Earth's interior. The most prominent material in the upper mantle is the source of mid-ocean ridge basalts (MORB). The MORB source is characterized by a general depletion in incompatible elements caused by the extraction of the continental crust from the mantle. At least three other isotopically distinct components are recognized in the suboceanic mantle. All three could be generated by the recycling of near surface materials (oceanic crust, pelagic sediments, continental lithospheric mantle) into the mantle by subduction. Therefore, the isotope data do not require a compositionally layered mantle, but neither do they deny the existence of such layering. Correlations between the volumetric output of plume volcanism with the reversal frequency of the Earth's magnetic field, and between the geographic distribution of isotopic variability in oceanic volcanism with seismic tomography suggest input of deep mantle material to surface volcanism in the form of deep mantle plumes. Volcanism on the continents shows a much wider range in isotopic composition than does oceanic volcanism. The extreme isotopic compositions observed for some continental magmas and mantle xenoliths indicate long-term (up to 3.3 Gyr) preservation of compositionally distinct material in thick (>200 km) sections of continental lithospheric mantle.     

Heat loss from the earth: A constraint on Archaean tectonics from the relation between geothermal gradients and the rate of plate production

The models suggested for the oceanic lithosphere which best predict oceanic heat flow and depth profiles are the constant thickness model and a model in which the lithosphere thickens away from the ridge with a heat source at its base. The latter is considered to be more physically realistic. Such a model, constrained by the observed oceanic heat flow and depth profiles and a temperature at the ridge crest of between 1100°C and 1300°C, requires a heat source at the base of the lithosphere of between 0.5 and 0.9 h.f.u., thermal conductivities for the mantle between 0.005 and 0.0095 cal cm⁻¹ °C⁻¹ s⁻¹ and a coefficient of thermal expansion at 840°C between 4.1 × 10⁻⁵ and 5.1 × 10⁻⁵ °C⁻¹. Plate creation and subduction are calculated to dissipate about 45% of the total earth heat loss for this model. The efficiency of this mechanism of heat loss is shown to be strongly dependent on the magnitude of the basal heat source. A relation is derived for total earth heat loss as a function of the rate of plate creation and the amount of heat transported to the base of plates. The estimated heat transport to the base of the oceanic lithosphere is similar to estimates of mantle heat flow into the base of the continental lithosphere. If this relation existed in the past and if metamorphic conditions in late Archaean high-grade terrains can be used to provide a maximum constraint on equilibrium Archaean continental thermal gradients, heat flow into the base of the lithosphere in the late Archaean must have been less than about 1.2–1.5 h.f.u. The relation between earth heat loss, the rate of plate creation and the rate of heat transport to the base of the lithosphere suggests that a significant proportion of the heat loss in the Archaean must have taken place by the processes of plate creation and subduction. The Archaean plate processes may have involved much more rapid production of plates only slightly thinner than at present.     

Seismic tomography and continental drift

Based on data of seismic tomography, the structure of the mantle flows of the contemporary Earth and the continental drift are calculated. Results of calculation of the contemporary motion of continents and their future drift for 150 Myr are presented. The present-day positions of six continents and the nine largest islands are taken as an initial state. The contemporary temperature distribution in the mantle is calculated according to the data of seismic tomography. The 3-D distribution of seismic wave velocities is converted into the density distribution and then into the temperature distribution. The Stokes equation is numerically solved for flows in a viscous mantle with floating continents for the given initial temperature distribution. In this way, the velocities of convective flows are determined in the entire present-day mantle and the surface distribution for the Earth’s heat flux is obtained. The reliability of the calculated flows in the mantle is estimated by the comparison of the calculated velocities of the contemporary continents and oceanic lithosphere with data of satellite measurements. Further, evolutionary equations of convection with floating continents were numerically solved. The calculated structure of mantle flows, temperature distribution, and position of continents are presented for a time moment 150 Myr in the future. The resulting successive changes in the position of continents in time show how islands (in particular, Japan and Indonesia) will be attached to continents and how continents will converge, exhibiting a tendency toward the formation of a new supercontinent in the southern hemisphere of the Earth.     

The structure of thermal convection and geotherms beneath a continent and ocean

The formation of the thermal cross section of the lithosphere and mantle upon the interaction between the mantle convection and the immobile continent surrounded by the oceanic lithosphere is studied by numerical modeling. The convective temperature and velocity fields and then the averaged geotherms for subcontinental and suboceanic regions up to the boundary with the core are calculated from the solution of convection equations with a jump in viscosity in the continental zone. Using the experimental data on the solidus temperature in the rocks of the upper mantle, the average thickness of the continental and oceanic lithosphere is estimated at 190 and 30 km, respectively. The effect of a hot spot formed in the subcontinental upper mantle at a depth of 250–500 km, which has not been previously noted, is revealed. Although the temperature in this zone is typically assumed to be close to adiabatic, the calculations show that it is actually higher than adiabatic by up to 200°C. The physical mechanism responsible for this effect is associated with the accumulation of convective heat beneath the thermally insulating layer of the continental lithosphere. The revealed anomalies can be important in studying the phase and mineral transformations at the base of the lithosphere and in the regional geodynamical reconstructions.     

Petrogenetic model of the Permian Tarim Large Igneous Province

Over the last two decades great strides have been made in characterizing the spatial distribution, time sequence,geochemical characteristics, mantle sources, and magma evolution processes for various igneous rocks in the Early Permian Tarim Large Igneous Province(TLIP). This work has laid a solid foundation for revealing the evolutionary processes and genetic models of large igneous provinces(LIPs). This study systematically demonstrates the two-stage melting model for the TLIP based on our previous research work and predecessor achievements, and highlights the two types of magmatic rocks within the TLIP.The two-stage melting model suggests that the formation of the TLIP is mantle plume related. The early hot mantle plume caused the low-degree partial melting of the lithosphere mantle, while in the later stage, the plume partially melted due to adiabatic uplift and decompression. Therefore, this model carries signatures of both the "Parana" and "Deccan" models in terms of mantle plume activity. During the early stage, the mantle plume provided the heat required for partial melting of sub-continental lithosphere mantle(SCLM), similar to the "Parana Model", while later the plume acted as the main avenue for melting, as in the "Deccan Model". Basalts that erupted in the first stage have higher 87Sr/86 Sr, lower 143Nd/144 Nd ratios, and are enriched in large ion lithophile elements and high field strength elements, indicating a possible origin from the enriched continental lithosphere mantle,similar to the Parana type geochemical features. The basic-ultrabasic intrusive rocks in the second stage exhibit lower 87Sr/86 Sr,higher 143Nd/144 Nd ratios relative to the basalts, consistent with the involvement of a more depleted asthenospheric material,such as a mantle plume, similar to the Deccan type geochemical features. The first stage basalts can be further subdivided into two categories, i.e., Group 1 and Group 2 basalts. Group 2 basalts have lower 87Sr/86 Sr and higher 143Nd/144 Nd ratios than Group 1 basalts, and lie between compositions of the Group 1 basalts and second stage magmatism. Group 2 basalts may be the intermediate component of the TLIP, and the whole TLIP is the result of plume and lithosphere interaction. Developing this petrogenetic model for the TLIP aids in comprehensively understanding its magmatism and deep geological and geodynamic processes. Furthermore, this work enriches the theories describing the origin of large igneous province and mantle plume activity.     

The snowball Earth aftermath: Exploring the limits of continental weathering processes

Carbonates capping Neoproterozoic glacial deposits contain peculiar sedimentological features and geochemical anomalies ascribed to extraordinary environmental conditions in the snowball Earth aftermath. It is commonly assumed that post-snowball climate dominated by CO₂ partial pressures several hundred times greater than modern levels, would be characterized by extreme temperatures, a vigorous hydrological cycle, and associated high continental weathering rates. However, the climate in the aftermath of a global glaciation has never been rigorously modelled. Here, we use a hierarchy of numerical models, from an atmospheric general circulation model to a mechanistic model describing continental weathering processes, to explore characteristics of the Earth system during the supergreenhouse climate following a snowball glaciation. These models suggest that the hydrological cycle intensifies only moderately in response to the elevated greenhouse. Indeed, constraints imposed by the surface energy budget sharply limit global mean evaporation once the temperature has warmed sufficiently that the evaporation approaches the total absorbed solar radiation. Even at 400 times the present day pressure of atmospheric CO₂, continental runoff is only 1.2 times the modern runoff. Under these conditions and accounting for the grinding of the continental surface by the ice sheet during the snowball event, the simulated maximum discharge of dissolved elements from continental weathering into the ocean is approximately 10 times greater than the modern flux. Consequently, it takes millions of years for the silicate weathering cycle to reduce post-snowball CO₂ levels to background Neoproterozoic levels. Regarding the origin of the cap dolostones, we show that continental weathering alone does not supply enough cations during the snowball melting phase to account for their observed volume.     

Stored mafic/ultramafic crust and early Archean mantle depletion

Both early and late Archean rocks from greenstone belts and felsic gneiss complexes exhibit positive ε_Nd values of +1 to +5 by 3.5 Ga, demonstrating that a depleted mantle reservoir existed very early. The amount of preserved pre-3.0 Ga continental crust cannot explain such high ε values in the depleted residue unless the volume of residual mantle was very small: a layer less than 70 km thick by 3.0 Ga. Repeated and exclusive sampling of such a thin layer, especially in forming the felsic gneiss complexes, is implausible. Extraction of enough continental crust to deplete the early mantle and its destructive recycling before 3.0 Ga ago requires another implausibility, that the sites of crustal generation and of recycling were substantially distinct. In contrast, formation of mafic or ultramafic crust analogous to present-day oceanic crust was continuous from very early times. Recycled subducted oceanic lithosphere is a likely contributor to present-day hotspot magmas, and forms a reservoir at least comparable in volume to continental crust. Subduction of an early mafic/ultramafic “oceanic” crust and temporary storage rather than immediate mixing back into undifferentiated mantle may be responsible for the depletion and high ε_Nd values of the Archean upper mantle. Using oceanic crustal production proportional to heat productivity, we show that temporary storage in the mantle of that crust, whether basaltic as formed by 5–20% partial melting, or partly komatiitic and formed by higher extents of melting is sufficient to balance an early depleted mantle of significant volume with ε_Nd at least +3.0.     

Anisotropic models of the upper mantle

Long period Rayleigh wave and Love wave dispersion data, particularly for oceanic areas, have not been simultaneously satisfied by an isotropic structure. In this paper available phase and group velocity data are inverted by a procedure which includes the effects of transverse anisotropy, anelastic dispersion, sphericity, and gravity. We assume that the surface wave data represents an azimuthal average of actual velocities. Thus, we can treat the mantle as transversely isotropic. The resulting models for average Earth, average ocean, and oceanic regions divided according to the age of the ocean floor, are quite different from previous results which ignore the above effects. The models show a low-velocity zone with age dependent anisotropy and velocities higher than derived in previous surface wave studies. The correspondence between the anisotropy variation with age and a physical model based on flow aligned olivine is suggestive. For most of the Earth SH > SV in the vicinity of the low-velocity zone. Neat the East Pacific Rise, however, SV > SH at depth, consistent with ascending flow. Anisotropy is as important as temperature in causing radial and lateral variations in velocity. The models have a high velocity nearly isotropic layer at the top of the mantle that thickens with age. This layer defines the LID, or seismic lithosphere. In the Pacific, the LID thickens with age to a maximum thickness of ～50 km. This thickness is comparable to the thickness of the elastic lithosphere. The LID thickness is thinner than derived using isotropic or pseudo-isotropic procedures. A new model for average Earth is obtained which includes a thin LID. This model extends the fit of a PREM, type model to shorter period surface waves.     

Anomalously nucleogenic neon in North Chile Ridge basalt glasses suggesting a previously degassed mantle source

Fresh basalt glasses from the North Chile Ridge (NCR) in the southeastern Pacific have Ne isotopic compositions distinctly different from typical mid-ocean ridge basalts (MORB). In a three-isotope plot of ²⁰Ne/²²Ne vs. ²¹Ne/²²Ne, the NCR data define a correlation line with a slope smaller than that of the MORB correlation line, i.e. their Ne composition is more nucleogenic than that of MORB. ³He/⁴He ratios are slightly lower than the MORB average, whereas in a few stepwise heating fractions very high ⁴⁰Ar/³⁶Ar ratios up to 28,000 are found. One model to explain the data assumes contamination of the NCR mantle source by material from the continental or oceanic crust, but in addition to difficulties with quantitatively reconciling the noble gas patterns with such a model it seems unable to account for some geochemical characteristics of NCR basalts reported earlier [Bach et al., Earth Planet. Sci. Lett. 142 (1996) 223–240], such as depletions in highly incompatible elements and unradiogenic Sr isotope compositions. Therefore we favor the scenario of a mantle source which was depleted and degassed previously, possibly as a residue from mantle melting beneath the southern East Pacific Rise that was transported to the NCR and melted again. The time during which such a depleted reservoir would have to be separated from the MORB mantle is estimated at 10–100 Ma based on U/Th–Ne systematics, in reasonable agreement with the time scale deduced from the formation history of the NCR and the temporal evolution of the southeast Pacific.