FeO and H2O and the homogeneous accretion of the earth

We present new shock devolatilization recovery data for brucite (Mg(OH)₂) shocked to 13 and 23 GPa. These data combined with previous data for serpentine (Mg₃Si₂O₅(OH)₄) are used to constrain the minimum size terrestrial planet for which planetesimal infall will result in an impact-generated water atmosphere. Assuming a chondritic abundance of minerals including 3–6%, by mass water, in hydrous phyllosilicates, we carried out model calculations simulating the interaction of metallic iron with impact-released free water on the surface of the accreting Earth. We assume that the reaction of water with iron in the presence of enstatite is the prime source of the terrestrial FeO component of silicates and oxides. Lower and upper bounds on the terrestrial FeO budget are based on mantle FeO content and possible incorporation of FeO in the outer core. We demonstrate that the iron-water reaction would result in the absence of atmospheric/hydrospheric water, if homogeneous accretion is assumed. In order to obtain10²⁵g of atmospheric water by the end of accretion, slightly heterogeneous accretion with initially 36% by mass iron planetesimals, as compared to a homogeneous value of 34% is required. Such models yield final FeO budgets, which either require a higher FeO content of the mantle (17 wt.%) or oxygen as a light element in the outer core of the Earth.     

Composition of the core,I. Solubility of oxygen in molten iron at high temperatures

Experiments on the solubility of FeO in molten iron have been carried out at temperatures between 2100 and 2550°C. The results show that liquid FeO is extensively soluble in molten iron at 2500°C and indicate that they probably become completely miscible above 2800°C. Liquid iron in equilibrium with crystalline magnesiowüstite (Mg_0.8Fe_0.2)O which is believed to be an important mineral in the lower mantle, would dissolve about 14 mol.% of FeO at 2500°C and 40 mol.% of FeO at 2800°C. The geochemical implications of these results are discussed. It is concluded that the outer core probably contains a large amount of dissolved FeO and that oxygen is probably the principal light element in the outer core.     

Heterogeneous accretion,composition and core–mantle differentiation of the Earth

A model of core formation is presented that involves the Earth accreting heterogeneously through a series of impacts with smaller differentiated bodies. Each collision results in the impactor's metallic core reacting with a magma ocean before merging with the Earth's proto-core. The bulk compositions of accreting planetesimals are represented by average solar system abundances of non-volatile elements (i.e. CI-chondritic), with 22% enhancement of refractory elements and oxygen contents that are defined mainly by the Fe metal/FeO silicate ratio. Based on an anhydrous bulk chemistry, the compositions of coexisting core-forming metallic liquid and peridotitic silicate liquid are calculated by mass balance using experimentally-determined metal/silicate partition coefficients for the elements Fe, Si, O, Ni, Co, W, Nb, V, Ta and Cr. Oxygen fugacity is fixed by the partitioning of Fe between metal and silicate and depends on temperature, pressure and the oxygen content of the starting composition. Model parameters are determined by fitting the calculated mantle composition to the primitive mantle composition using least squares minimization. Models that involve homogeneous accretion or single-stage core formation do not provide acceptable fits. In the most successful models, involving 24 impacting bodies, the initial 60–70% (by mass) of the Earth accretes from highly-reduced material with the final 30–40% of accreted mass being more oxidised, which is consistent with results of dynamical accretion simulations. In order to obtain satisfactory fits for Ni, Co and W, it is required that the larger (and later) impactor cores fail to equilibrate completely before merging with the Earth's proto-core, as proposed previously on the basis of Hf-W isotopic studies. Estimated equilibration conditions may be consistent with magma oceans extending to the core–mantle boundary, thus making core formation extremely efficient. The model enables the compositional evolution of the Earth's mantle and core to be predicted throughout the course of accretion. The results are consistent with the late accretion of the Earth's water inventory, possibly with a late veneer after core formation was complete. Finally, the core is predicted to contain ~ 5 wt.% Ni, ~ 8 wt.% Si, ~ 2 wt.% S and ~ 0.5 wt.% O.     

Volume-temperature relationship for iron at 330 GPa and the Earth's core density deficit

Estimates of core density deficit (cdd) of the Earth's outer core recently reported by Anderson and Isaak [Another look at the core density deficit of Earth's outer core, Phys. Earth Planet Int. 131 (2002) 19-27] are questionable in view of the serious errors in the pressure-volume and bulk modulus data due to an inadequacy in the calibration process used by Mao et al. [Static compression of iron to 300 GPa and Fe_0.8Ni_0.2 alloy to 200 GPa: implications for the core, J. Geophys. Res. 94 (1990) 21737-21742]. The data used by Anderson and Isaak deviate significantly from the corresponding values derived from seismology. In the present study we have used the input data on density, isothermal bulk modulus and its pressure derivative from Stacey and Davis [High pressure equations of state with application to lower mantle and core, Phys. Earth Planet Int. 142 (2004) 137-184] which are consistent with the seismological data. Volumes of hexagonal close-packed iron have been calculated at different temperatures under isobaric conditions at P = 330 GPa, the inner core boundary (ICB) pressure using the relationship between thermal pressure and volume expansion based on the lattice potential theory originally due to Born and Huang [Dynamical Theory of Crystal Lattices, Oxford University Press, Oxford, 1954, p. 50]. The formulation for thermal pressure used by Anderson and Isaak has been modified by taking into account the variations of thermal expansivity α and isothermal bulk modulus K_T with temperature. Values of cdd are then estimated corresponding to different temperatures ranging from 4000 to 8000 K. The results for cdd at different temperatures obtained in the present study are significantly higher than those estimated by Anderson and Isaak suggesting that the cdd for the Earth's outer core is nearly 10%. The effects of nickel when an Fe-Ni alloy replaces Fe are estimated and found to be insignificant.     

Temporal variation in primary magma compositions in the northeast Japan Arc

Abstract A remarkable temporal variation in primary magma compositions has been found in the Northeast Japan arc. The trench-side magmas have become more enriched in FeO* and the backarc-side magmas have become more depleted in FeO* while retaining almost constant SiO, levels for the last ∼20 million years. In order to understand the origin of the temporal variation, FeO* and SiO, contents in partial melts are modeled for an adiabatically-rising mantle as a function of potential temperature and original composition of the mantle material. The result demonstrates that the primary magmas that are more depleted in FeO* were derived from the mantle materials either at lower potential temperatures or with compositions more depleted in basaltic components. A possible mechanism for the inferred primary magma variability is the change in depth intervals with time of magma production in a compositionally-layered mantle wedge; greater degrees of depletion at a greater depth is reconciled with a probable thermal regime in the mantle wedge.     

The system iron-enstatite-water at high pressures and temperatures—formation of iron hydride and some geophysical implications

The system iron-enstatite-water was investigated at pressures around 5 GPa and at temperatures ranging from 1000 to 1200°C, using several different kinds of starting materials. Quenched samples showed the coexistence of iron, olivine and pyroxene. Synthesis of the Fe-containing olivine in the run products proves that a series of reactions, Fe + H₂O → FeH_x + FeO and FeO + MgSiO₃ → (Mg, Fe)₂SiO₄, have taken place. Spherical “balls of iron” were observed in the 1200°C run. This strongly indicates that the melting temperature of iron decreased by ～ 500 K by the possible dissolution of hydrogen. Following geophysical implications are derived from these experimental results. If water was retained in the hydrous minerals in the primordial material, the iron-water reaction is expected to occur throughout the core-formation process. The reaction product FeH_x will melt and then sink to form a proto-core and iron oxide will be dissolved in the Earth's mantle. The dissolution of hydrogen in the Earth's core is a natural consequence of the core-formation process.     

Coupled Os–Os enrichments in the Earth’s mantle – core–mantle interaction or recycling of ferromanganese crusts and nodules?

¹⁸⁶Os enrichments in volcanic rocks and peridotite-derived iridosmine grains have been attributed to contributions from Earth’s outer core to the mantle, and apparently constrain the scale of mantle convection and an early timing for inner–outer core segregation more than 3.5 Gyr ago. Here, we highlight that marine ferromanganese crusts and nodules are characterised by high Pt/Os ratios and Pt–Os contents that develop much larger ¹⁸⁶Os excesses over geological time (≥0.2%/Gyr) than those hypothesised for Earth’s outer core (<0.005–0.01%/Gyr). ¹⁸⁷Os/¹⁸⁸Os ratios in ferromanganese crusts are radiogenic due to sequestering of continental Os from seawater. Similarly, ancient ferromanganese materials may have had ¹⁸⁶Os excesses (>0.1%) as a result of high Pt/Os ratios in continental crust, even prior to in-growth of ¹⁸⁶Os after formation due to their high Pt/Os ratios. Past recycling of small amounts of these materials into the Earth’s mantle will produce coupled ¹⁸⁷Os–¹⁸⁶Os excesses and little change in Re and platinum-group-element concentrations, as observed in Hawaiian picrites, and in contrast to the predicted result of outer core addition to the mantle. ¹⁸⁷Os and ¹⁸⁶Os enrichments in the Hawaiian mantle source are potentially consistent with it comprising recycled oceanic lithosphere, pelagic sediments and ferromanganese materials, and questions the notion that Os isotopes can be used to uniquely identify core–mantle interactions and the depth at which mantle sources for volcanism originate.     

甜水海钻孔TS95若干地球化学指标的非线性分析与意义

应用关联维分析和R／S分析对西昆伦甜水海孔TS95的4项指标（FeO,Fe2O3,Fe2O3/FeO和有机碳）进行分析,发现它们存在明显的混沌特征和Hurst现象,4项指标的混沌吸引子分别为：FeO－2．8;Fe2O3－3．2;Fe2O／FeO－2．9;有机碳2．5,饱和嵌入维数为5到6,表明控制该地区气候环境演化的动力系统是由有限维数确定的混沌系统,构筑该系统所需变量至少3到4个,最多5到6个,这种体上与深海氧同位素比较接近,与黄土剖面一些指标差异较大,4项指标折Hurst指数（FeO－0．85;Fe2O3－0．76;Fe2O3／FeO－0．65;有机碳－0．74）显示该地区气候环境演化存在明显的持续性成分,这与早期的研究一致,很可能与青藏高原的构造抬升有关,表明除了全球性气候环境波动背景外,局域性因素,如肝藏高原构造抬升和未海地区水系变迁,都会在各种气候环境演化的替代性指标上有所反映。4项指标Hurst指数的差异可能反映FeO容易迁移,受到了流域水系变化的影响,因而有机碳和Fe2O3可能是反映这种环境演化持续性分更好的指标。     

Metal/silicate fractionation in the solar system

Fractionation between the metal and silicate components of objects in the inner solar system has long been recognized as a necessity in order to explain the observed density variations of the terrestrial planets and the H-group, L-group dichotomy of the ordinary chondrites. This paper discusses the densities of the terrestrial planets in light of current physical and chemical models of processes in the solar nebula. It is shown that the observed density trends in the inner solar system need not be the result of special fractionation processes, and that the densities of the planets may be direct results of simultaneous application of both physical and chemical restraints on the structure of the nebula, most notably the variation of temperature with heliocentric distance. The density of Mercury is easily attributed to accretion at temperatures so high that MgSiO₃ is only partially retained but Fe metal is condensed. The densities of the other terrestrial planets are shown to be due to different degrees of retention of S, O and H as FeS, FeO and hydrous silicates produced in chemical equilibrium between condensates and solar-composition gases. It is proposed that Mercury and Venus Have cores of Fe⁰, Earth has a core of Fe⁰ containing substantial amounts of FeS, and Mars has a quite small core of FeS with more FeO in its mantle than in Earth's. Geophysical and geochemical consequences of these conclusions are discussed.     

Gravitational energy of core evolution: implications for thermal history and geodynamo power

Using density–pressure relationships for mantle silicate and core alloy closely matching PREM we have constructed six models of the Earth in different evolutionary states. Gravitational energies and elastic strain energies are calculated for models with homogeneous composition, separated mantle and liquid core, separated inner and outer cores with the inner core either liquid or solid and models with increased densities, representing cooling of either the mantle or core. In this way we have isolated the gravitational energy released by each of several evolutionary processes and subtracted the consequent increase in strain energy to obtain the net energy released as heat or geodynamo power. Radiogenic heat (∼7.8×10³⁰ J) is found to contribute only about 25% of the total heat budget, the balance originating as residual gravitational energy from the original accretion and from core separation (14×10³⁰ J). The total energy of compositional convection, driven by inner core formation, is 3.68×10²⁸ J and this is the most important (or even the only) energy source for the dynamo for the most recent 2 billion years. It appears unlikely that the inner core existed much before that time. The total net (gravitational minus strain) energy released in the core by the process of inner core formation, 11.92×10²⁸ J, is not much less than the thermal energy released in this process, 15.1×10²⁸ J. In the mantle the net (gravitational minus strain) energy released by thermal contraction is about 20% of the heat release. All of the numerical results are presented in a manner that allows simple rescaling to any revised density estimates.     

Chemical record of the projectile in the graded fall-back sedimentary unit from the Ries Crater, Germany

The so-called “Graded Unit” comprises 17.2 m of the core from the research borehole Nördlingen 1973 and is thought to represent fall-back material, which has been sorted similar to tephra. Samples from six levels of these air-borne impact debris sediments in the Ries Crater have been analyzed for Ir, Ni, Co, Cr and Fe. The concentrations of these elements are negatively correlated with average grain size and accordingly the depth of the samples indicating the presence of a surface correlated component and are hence consistent with vapour deposition of meteoritic material on small dust grains. The largest Ir concentration found is 230 pg/g, a factor of 28 higher than the indigenous concentration. Enrichment factors in the uppermost layer of the “Graded Unit” are about 28 for Ir, approximately 2 for Fe, Co, Ni and 1.5 for Cr. If the enrichment of all elements determined were of meteoritic origin, then the “net meteoritic composition” is most closely matched by aubrites as suggested previously. However, it seems more likely that the observed concentration trends of Cr, Fe, Co and Ni are due to mineral sorting. Therefore the only as yet certain indication of meteoritic material is the presence of Ir, which is not sufficient to further constrain the nature of the projectile.     

Geochemistry of the oldest MORB and OIB in the Isua Supracrustal Belt, southern West Greenland: Implications for the composition and temperature of early Archean upper mantle

Abstract Recent geological investigations of the Isua Supracrustal Belt (3.8 Ga), southern West Greenland, have suggested that it is the oldest accretionary complex on earth, defined by an oceanic plate‐type stratigraphy and a duplex structure. Plate history from mid‐oceanic ridge through plume magmatism to subduction zone has been postulated from analysis of the reconstructed oceanic plate stratigraphy in the accretionary complex. Comparison between field occurrence of greenstones in modern and ancient accretionary complexes reveals that two types of tholeiitic basalt from different tectonic settings, mid‐oceanic ridge basalt (MORB) and oceanic island basalt (OIB), occur. This work presents major, trace and rare earth element (REE) compositions of greenstones derived from Isua MORB and OIB, and of extremely rare relict igneous clinopyroxene in Isua MORB. The Isua clinopyroxenes (Cpx) have compositional variations equivalent to those of Cpx in modern MORB; in particular, low TiO₂ and Na₂O contents. The Isua Cpx show slightly light (L)REE‐depleted REE patterns, and the calculated REE pattern of the host magma is in agreement with that of Isua MORB. Analyses of 49 least‐altered greenstones carefully selected from approximately 1200 samples indicate that Isua MORB are enriched in Al₂O₃, and depleted in TiO₂, FeO*, Y and Zr at the given MgO content, compared with Isua OIB. In addition, Isua MORB show an LREE‐depleted pattern, whereas Isua OIB forms a flat REE pattern. Such differences suggest that the Early Archean mantle had already become heterogeneous, depending on the tectonic environment. Isua MORB are enriched in FeO compared with modern MORB. Comparison of Isua MORB with recent melting experiments shows that the source mantle had 85–87 in Mg? and was enriched in FeO. Potential mantle temperature is estimated to be approximately 1480°C, indicating that the Early Archean mantle was hotter by at most approximately 150°C than the modern mantle.     

Inversion of seismic and gravity data for the composition and core sizes of the Moon

We model the internal structure of the Moon, initially homogeneous and later differentiated due to partial melting. The chemical composition and the internal structure of the Moon are retrieved by the Monte-Carlo inversion of the gravity (the mass and the moment of inertia), seismic (compressional and shear velocities), and petrological (balance equations) data. For the computation of phase equilibrium relations and physical properties, we have used a method of minimization of the Gibbs free energy combined with a Mie-Gr@uneisen equation of state within the CaO-FeO-MgO-Al₂O₃-SiO₂ system. The lunar models with a different degree of constraints on the solution are considered. For all models, the geophysically and geochemically permissible ranges of seismic velocities and concentrations in three mantle zones and the sizes of Fe-10%S core are estimated. The lunar mantle is chemically stratified; different mantle zones, where orthopyroxene is the dominant phase, have different concentrations of FeO, Al₂O₃, and CaO. The silicate portion of the Moon (crust + mantle) may contain 3.5–5.5% Al₂O₃ and 10.5–12.5% FeO. The chemical boundary between the middle and the lower mantle lies at a depth of 620–750 km. The lunar models with and without a chemical boundary at a depth of 250–300 km are both possible. The main parameters of the crust, the mantle, and the core of the Moon are estimated. At the depths of the lower mantle, the P and S velocities range from 7.88 to 8.10 km/s and from 4.40 to 4.55 km/s, respectively. The radius of a Fe-10%S core is 340 ± 30 km.     

A new high-pressure phase of Ca2Al2SiO7 and implications for the earth''s interior

Gehlenite (Ca₂Al₂SiO₇) has been found to transform to a new phase at pressures greater than 100 kbar and at about 1000°C, using a diamond-anvil pressure cell coupled with laser heating. The atoms of the new phase appear to be arranged in a perovskite-related structure similar to that described for Na₂Ti₃O₇. The structure probably consists of layers of (Al₂SiO₇)⁴⁻, which are built up from blocks of edge-sharing (Al, Si)O₆ octahedra and these blocks are joined by common octahedra corners. A small cubic unit cell with a = 3.719 ± 0.004 Å indexes completely the strong lines of the powder diffraction pattern, and a superlattice with a = 14.88 ± 0.02 Å satisfies all the observed weak lines in addition to the strong ones. However, the cell may be pseudocubic. The small cell contains a half of the gehlenite formula while the large cell contains 32 gehlenite formulae. Hence the molar volume for the new phase of Ca₂Al₂SiO₇ is calculated to be 61.96 ± 0.20 cm³ at atmospheric pressure and room temperature. The new sodium titanate-type structure is probably more closely packed than an ordinary perovskite-type structure in which all octahedral corners are shared. This view is strongly supported by the very great density of this new phase, which is about 8% denser than the equivalent mixture of CaAl₂O₄ (calcium ferrite type) plus CaSiO₃ (cubic perovskite type). The new phase is probably the most closely packed silicate known. Mg₂SiO₄ (spinel) was found to transform to an assemblage containing MgSiO₃ (perovskite) plus MgO (periclase) at P-T conditions equivalent to the upper part of the lower mantle. By reacting with MgO, the perovskite modification of both MgSiO₃ and MgSiO₃ · xAl₂O₃ may adopt the sodium titanate structure at the still greater depths of the lower mantle. If the sodium titanate structures of Mg₂(Al₂Si)O₇ and Mg₂(MgSi₂)O₇ are present in the deep part of the lower mantle, MgO does not exist as a separate phase at the mantle-core boundary. This might be an obstacle to the possibility of dissolving these oxides (specifically the FeO component) in the molten Fe in the outer core as suggested by geophysical and geochemical studies of the earth's interior. The mechanism for developing the chemical plumes in the deep mantle proposed by Anderson does not appear to be consistent with studies of phase transformations in Ca-Al-rich compounds as outlined in this paper.     

Thermal histories of the core and mantle

Recognition that the cooling of the core is accomplished by conduction of heat into a thermal boundary layer (D″) at the base of the mantle, partly decouples calculations of the thermal histories of the core and mantle. Both are controlled by the temperature-dependent rheology of the mantle, but in different ways. Thermal parameters of the Earth are more tightly constrained than hitherto by demanding that they satisfy both core and mantle histories. We require evolution from an early state, in which the temperatures of the top of the core and the base of the mantle were both very close to the mantle solidus, to the present state in which a temperature increment, estimated to be ～ 800 K, has developed across D″. The thermal history is not very dependent upon the assumption of Newtonian or non-Newtonian mantle rheology. The thermal boundary layer at the base of the mantle (i.e., D″) developed within the first few hundred million years and the temperature increment across it is still increasing slowly. In our preferred model the present temperature at the top of the core is 3800 K and the mantle temperature, extrapolated to the core boundary without the thermal boundary layer, is 3000 K. The mantle solidus is 3860 K. These temperatures could be varied within quite wide limits without seriously affecting our conclusions. Core gravitational energy release is found to have been remarkably constant at ～ 3 × 10¹¹ W. nearly 20% of the core heat flux, for the past 3 × 10⁹ y, although the total terrestrial heat flux has decreased by a factor of 2 or 3 in that time. This gravitational energy can power the “chemical” dynamo in spite of a core heat flux that is less than that required by conduction down an adiabatic gradient in the outer core; part of the gravitational energy is used to redistribute the excess heat back into the core, leaving 1.8 × 10¹¹ W to drive the dynamo. At no time was the dynamo thermally driven and the present radioactive heating in the core is negligibly small. The dynamo can persist indefinitely into the future; available power 10¹⁰ y from now is estimated to be 0.3 × 10¹¹ W if linear mantle rheology is assumed or more if mantle rheology is non-linear. The assumption that the gravitational constant decreases with time imposes an implausible rate of decrease in dynamo energy. With conventional thermodynamics it also requires radiogenic heating of the mantle considerably in excess of the likely content of radioactive elements.     

Equation of state fits to the lower mantle and outer core

The lower mantle and outer core are subjected to tests for homogeneity and adiabaticity. An earth model is used which is based on the inversion of body waves and Q-corrected normal-mode data. Homogeneous regions are found at radii between 5125 and 4825 km, 4600 and 3850 km, and 3200 and 2200 km. The lower mantle and outer core are inhomogeneous on the whole and are only homogeneous in the above local regions.Finite-strain and atomistic equations of state are fit to the homogeneous regions. The apparent convergence of the finite-strain relations is examined to judge their applicability to a given region. In some cases the observed pressure derivatives of the elastic moduli are used as additional constraints. The effect of minor deviations from adiabaticity on the extrapolations is also considered. An ensemble of zero-pressure values of the density and seismic velocities are found for these regions. The range of extrapolated values from these several approaches provides a measure of uncertainties involved.     

The amplitude ratioPcP/P and the core-mantle boundary

Summary Using the Haskell matrix formulation, theoretical reflection coefficient curves have been calculated for a multi-layered core-mantle boundary for comparison with observational data. Two cases are considered, first when the shear velocity in the core is equal to zero and second when the core has a finite rigidity. If the velocity contrast is large between the imbedded layer and the mantle, the reflection coefficient curves for the multi-layered medium are irregular in shape as compared to those for two half-spaces, representing the core and the mantle, respectively. The reflection coefficient curves show an oscillatory character if the imbedded layer is thick and has a high velocity contrast.The observational data consist of short-period vertical-component seismograph records ofP andPcP from nuclear explosions in the Aleutian chain, Nevada, Novaya Zemlya, Kazakh and Sahara. Attenuation and geometrical spreading are taken into consideration. Four different models for the quality factorQ are applied to the observational data. The data are found to be much affected by theQ-model used for the corrections.Based on proposedQ-values, a model for the core-mantle boundary is found, characterized by two low-velocity layers at the bottom of the mantle. The thicknesses are 16.10 km (outer layer) and 19.96 km (inner layer), the compressional wave velocities 12.17 km/sec and 10.94 km/sec and the shear wave velocities are 6.29 km/sec and 5.33 km/sec, respectively. A better fit to this model is found when in addition the shear velocity in the outer core is 2.20 km/sec and the density ratio at the core-mantle boundary is 1.07. In other words, the observations favour a layer of finite rigidity in the outer core rather than a fluid one.     

Chemical composition of Archaean komatiites: Implications for early history of the earth and mantle evolution

Estimates of the chemical composition of the Archaean mantle, derived from elemental abundance ratios in komatiites combined with ultramafic xenolith data, support a model of a multistage heterogeneous accretion history of the Earth and synchronous core formation, 4.6 Ga ago.Most refractory lithophile element abundance ratios in komatiites and xenoliths are close to chondritic except for V/Ti and Ca/Al. Depletion of vanadium is likely due to its partial incorporation into the iron core; whereas fractionation of Ca/Al observed in Archaean Al-undepleted komatiites (1.20 times chondrites) and in some modern fertile spinel lherzolite xenoliths (1.15 times chondrites) could be due to small amounts of garnet (rich in Al but poor in Ca) segregation into the lower mantle during partial or complete melting of the upper mantle in the very early history of the Earth. However, this process may have had only a small effect on the overall chemical composition of the upper mantle.Simultaneous occurrence of early Archaean Al-undepleted (Al/Ti chondrites) and Al-depleted (Al/Ti 0.5 chondrites, and depletion of Sc and heavy REE) peridotitic komatiites in the Barberton area, S. Africa, and late Archaean Newton Township, Canada, argue against derivation of peridotitic komatiites from a circum-global magma ocean. Garnet separation from a mantle diapir which intersects the solidus at great depth ( 200 km) in a hotter early Archaean mantle could explain the chemical characteristics of Al-depleted komatiites. Alternatively, these two types of komatiites could have been derived from different layers in a fractionated mantle. A limited amount of Hf isotope data for Archaean komatiites seems to suggest that both mechanisms are important. This chemically and minerallogically layered mantle, if it existed, was homogenized by mantle convection after early Archaean times.Constant P₂O₅/TiO₂, Ni/Mg, Co/Mg, Fe/Mg ratios (siderophile/lithophile) and PGE abundances, estimated for the mantle sources of komatiites from Archaean to modern times, strongly argue against continuous growth of the Earth's core since the early Archaean.Extensive crustal contamination might have been involved in the generation of Archaean-early Proterozoic siliceous high magnesian basalts with “boninite affinity”. However, involvement of chemically modified ancient continental lithosphere may also be important in the generation of these basalts.     

A multi-layered kinematic dynamo model: implications of a stratified upper layer in the Earth's core

It has been suggested that there exists a stably stratified electrically conducting layer at the top of the Earth's outer fluid core and that lateral temperature gradients in the lower mantle is capable of a driving thermal-wind-type flow near the core–mantle boundary. We investigate how such a flow in a stable layer could influence the geomagnetic field and the geodynamo using a very simple two-dimensional kinematic dynamo model in Cartesian geometry. The dynamo has four layers representing the inner core, convecting lower outer core, stable upper core, and insulating mantle. An α² dynamo operates in the convecting outer core and a horizontal shear flow is imposed in the stable layer. Exact dynamo solutions are obtained for a range of parameters, including different conductivities for the stable layer and inner core. This allows us to connect our solutions with known, simpler solutions of a single-layer α² dynamo, and thereby assess the effects of the extra layers. We confirm earlier results that a stable, static layer can enhance dynamo action. We find that shear flows produce dynamo wave solutions with a different spatial structure from the steady α² dynamos solutions. The stable layer controls the behavior of the dynamo system through the interface conditions, providing a new means whereby lateral variations on the boundary can influence the geomagnetic field.     

Earth tides and polar motion

We establish a general theory that describes the rotational motion of a layered, oblate, elastic Earth under the influence of tidal forces when account is taken of the liquid outer core. We obtain a linearized version of the Navier-Stokes equation; within it not only have we retained the Coriolis and centrifugal acceleration terms, but also have included the nutational terms. We also make use of the Euler equation for angular momentum to analytically relate the nutational motion of the rotational axis with the oscillations of the liquid core and obtain a constraint for the nutational amplitude. Consideration of the Poisson equation for density variation completes our analytical model.We primarily discuss the equations of motion for the liquid core and present the solution as the sum of two terms: one being a component of the spheroidal displacement field, the other of the toroidal field. We also formulate the equations valid for the solid mantle when rotational effects are included, and establish the boundary conditions that must hold at the various interfaces in order that a complete integration of the differential system of equations be accomplished.We assume that the outer core consists of an inviscid fluid and ignore the existence of any boundary layer. We do not impose, however, any restriction on the stratification of the fluid. The dynamical coupling between liquid core and solid mantle is represented by a torque which is generated by the forced oscillations within the liquid core; these oscillations are in turn triggered by the diurnal tides.The expected influence of the liquid core/solid mantle boundary on the nutational motion is discussed in view of Poincare's results concerning a liquid core surrounded by a rigid shell. Comparison is finally made of our model with Molodenskii's 1961 theory for a neutral core and the 1976 Shen-Mansinha nutational theory for an unrestricted core.