next up previous
Next: The Core-Mantle Boundary Up: Mineral physics of iron Previous: Anisotropy of the

Impurities in the Outer and Inner Core

The experimental results shown in Figure 6 of Mao et al. [1990] demonstrate that at room temperature there is a 100 GPa difference between the measured of iron and the PREM at inner core densities. We can expect the thermal pressure, , to be about at inner core temperatures at those densities [Anderson, 1995, Page 269]. So there is about a 40 GPa discrepancy to account for by light impurities. This problem has been examined by Jephcoat and Olson [1987], who proposed that a substantial amount of FeO would be needed to satisfy PREM's density at the inner core. Their analysis was based on the phase as the stable phase at core pressures. As described above, in the section on phase diagrams, the phase may be the dominant phase in the inner core. In that case, the quantity of impurities in the inner core will almost certainly be less than calculated by Jephcoat and Olson, [ Isaak et al., 1994].

Poirier [1994a] has critically reviewed the literature regarding the types and quantities of light impurities in the outer core. He concludes that there is no reason to consider that only one light element is in the outer core. Several light elements (Si, S, O, H, C) seem certain to be in the outer core. Poirier concludes, ``it is not even obvious that one element should be particularly dominant.''

Although not a light element, nickel is one of the elements we should expect to find in the inner and outer core. Mao et al. [1990] showed experimentally that an iron-nickel alloy, FeNi, has virtually the same relationship as pure iron up to 300 GPa at 300 K. Therefore the presence of nickel will not disturb the calculation of density effects due to other impurities, and so the geophysical model may safely ignore, at least as a first approximation, the presence of nickel as other alloys are considered.

From the shock wave experiments of Jeanloz and Ahrens [1980] on Fe and FeO, we know that about 45% FeO by weight is required in combination with 55% pure iron by weight to satisfy the outer core density data. Similarly, using the shock wave data of Ahrens [1979] on FeS and FeS, an even larger amount of sulfur can be mixed with iron to satisfy the outer core density. Therefore oxygen and/or sulfur have been assumed by several [ Ringwood and Hibberson, 1991; Knittle and Jeanloz, 1991; Usselman, 1975] as likely candidates for impurities of the inner core in order to satisfy the density deficit. Poirier's [1994a] calculations show that if silicon, oxygen, and sulfur were simultaneously present, a number of combinations of mass fractions would give the required 10% density deficit.

It is of interest to consider how these light elements, when mixed with iron, might affect the melting temperature of the resulting alloy. Boehler [1992] measured melting of the Fe-FeS and Fe-FeO systems up to 50 GPa. His 's were substantially lower than those for these same systems as determined by Knittle and Jeanloz [1991]. For Fe0, Boehler's curve passed through the point at 17 GPa measured by Ringwood and Hibberson [1990]. Boehler found that for FeO, increases with faster than it does for pure iron itself, at least up to 80 GPa, the limit of the experiment. For FeS, the rate of increase of with is less than for pure iron. Therefore, it is possible that FeO in the inner core will not diminish the value of significantly, although it will diminish the density (see also a similar conclusion by Knittle and Jeanloz, 1991).

However, other abundant elements may cause the lowering of . Fukai [1992] reported that the presence of iron hydride substantially lowered the melting temperature at high pressure. Yagi et al. [1994] reported that when Fe, MgSiO, and water are mixed at high pressure, both FeH and FeO are formed, and the resulting FeH has a much lower and a lower density than pure iron. Similarly, Wood [1993], who presented a detailed analysis of Fe-C and the Fe-C-S systems at high pressure, showed that carbon is likely present in the inner core. Its effect would be to depress the melting point significantly.

In trying to estimate the melting temperature depression of iron by the concentration of minor elements, we have no compelling theory to help us. We cannot even usefully predict whether a given system will have a eutectic or act as a solid solution at very high pressures. There are a number of experiments that show a miscibility gap between the metallic F and the ionic oxide (such as FeO). In general, the experiments show that the miscibility gap decreases as the pressure increases, indicating that the eutectic behavior may disappear at sufficiently high P. For example, Wood [1993] found experimentally that the eutectic for the Fe-S-C system rises with pressure from to 15 GPa, and the miscibility gap decreases. He suggested that Fe-S-C liquids could form one stable liquid at core temperatures. The experimental data are incomplete, and such as exist are limited to pressures up to about 50 GPa.

If the eutectic type phase diagrams are valid at core pressure, then calculations and measurements of described in the previous section may not help very much in determining core temperatures. If, on the other hand, eutectic phase diagrams do not persist up to core pressures and impurities go into the iron lattice as solutes, making solid solutions, then the resulting can be calculated by finding the freezing point depression based on the melting curve of pure iron.

There is an experimental result suggesting that solid solutions are favored at very high pressures for iron systems. W.W. Anderson and Ahrens [1994] reported that the bulk modulus of pure liquid iron at outer core conditions is virtually identical to the seismic bulk modulus of the outer core. Thus of the outer core is independent of the concentration of the impurities in the iron. D.L. Anderson [1976] showed empirically that for silicate solid solutions (garnets, oxides, and pyroxenes), the bulk modulus is virtually independent of iron concentrations. This was verified in detail experimentally for olivine, pyroxenes, and garnets [ Isaak, 1992; Weidner et al., 1982]. The theory for the - relationship in solid solutions was established by Shankland [1972], who showed that the law holds, where is the bulk sound velocity for solid solutions. This law is equivalent to saying that is independent of for solid solutions at constant crystal structure. It is derived from Debye theory for a solid solution, where the crystal structure is held constant, and the concentration of impurities is varied. Iron fulfills the criteria of a Debye solid. Therefore, the W.W. Anderson and Ahrens [1994] report on bulk modulus supports the assumption that iron, under high compression and with a few percent impurities, is a solid solution, and eutectics need not be considered.

In addition, there are theoretical results indicating the onset of solubility conditions at high pressure. Boness and Brown [1990] and Sherman [1991] suggested that Fe-S exhibits continuous solution at core pressures. Using the van der Laar equation, Boness and Brown [1990] estimated that the melting depression drop from pure iron to an outer core composition is about 1000 K.

A solid solution is favored for an oxide in the core if the bonding is metallic. There has been a persistent effort to see whether FeO becomes more metallic at the pressure of its transition at about 80 GPa. The polymorphic wüstite transition found by Jeanloz and Ahrens [1980] was described as probably from to (NaCl to CsCl), implying ionic bonding in the high pressure phase. Jackson and Ringwood [1981] presented arguments that the transition should be to a nickel arsenide structure with a more covalent bond involving spin pairing.

Isaak et al. [1993], in a first principles calculation, (LSDA) investigated the phase stability of wüstite at high pressure. They predicted a transition to a metallic structure at 500 GPa. Sherman [1994], following the path of Isaak et al. [1993] but using the more general model called full potential linearized augmented wave (FLAPW) method of calculation, found the Ni-As structure at 125 GPa, which is metallic.

Fei et al. [1994] determined experimentally that wüstite transforms to the Ni-As-type structure at 93 GPa and 800 K, indicating a substantial increase in the metallicity of the structure.

Thus the evidence, both experimental and theoretical, has accumulated in favor of a solid solution of iron and other light impurities for the composition of the core.

Poirier and Shankland [1994] speculated that the temperature of the inner core at the ICB is the of pure iron less the freezing point depression due to about 10% light impurities, that is, a drop of 500--1000 degrees. Braginsky and Roberts [1995] adopted a melting point depression of 700 K. Stacey [1995], following Poirier's [1986] suggestion on melting depression, found the value of at the ICB to be about 5000 K. To quote Stacey, ``clearer evidence is needed.''

Boehler [1993] estimated a value close to that of Poirier, , using a different approach. Using his experimental results on the variation of with pressure for the Fe-(Fe-FeO) system, he noted that is linear in density (a variation of the Kraut-Kennedy law), so a simple extrapolation of his measured results to the ICB pressure gave him 4640 K for the Fe-FeO system. He assumed no triple points or eutectics above his measurements terminating at 60 GPa. It might be that triple points occurring in pure iron as seen in the phase diagram are absent in the Fe-FeO system. In the Poirier calculation, light impurities are not specified, except as a concentration value, whereas in the Boehler extrapolation, an assumed light impurity is used.

These estimates are based on incomplete data, however. To find the quantity of impurities, we must know the crystallographic structure of pure iron that exists at inner core conditions. Four years ago, it was universally assumed that the phase would be hcp, and the high pressure density was promptly determined [ Mao et al., 1990]. But now we see that the structure might be fcc. The arguments in favor of this are as follows: at ambient conditions, pure iron (bcc) tends to change to fcc when alloyed with small amounts of nickel, copper, palladium, or platinum; the fcc phase might be the phase found in pure iron at high conditions; or the high P, high phase could be an unforeseen phase, such as a simple cubic cell with a large number of atoms per cell, as found, for example, in manganese at high and high [ Young, 1991]. Also we must know more about the kinds of impurities and their relative concentrations in the inner core.

Current thinking of geophysicists ignores potassium as a significant impurity, although it was considered a viable alternative two decades ago. (A good discussion of the feasibility of potassium in the core is found in Jacobs [1987]).

To sum up, in spite of reservations, these various estimates of melting point depression lead to a value of near 5000 K. Braginsky and Roberts [1995] adopted .



next up previous
Next: The Core-Mantle Boundary Up: Mineral physics of iron Previous: Anisotropy of the



U.S. National Report to IUGG, 1991-1994
Rev. Geophys. Vol. 33 Suppl., © 1995 American Geophysical Union