Abstract
Metallic melt containing iron (Fe) and carbon (C) may be present at depths greater than 250 km inside the Earth. Depending on its wetting behavior, such dense melt may be trapped locally or drain into deep mantle and core. Here, we report experimental data on the wetting behavior of Fe–C melt in silicates at the conditions of Earth’s mid-mantle between 10 and 23 GPa and 1600 and 1800°C. The measured dihedral angles of Fe–C melt in olivine, ringwoodite or bridgmanite and ferropericlase matrixes are 117 ± 14°, 120 ± 14°, and 107 ± 16°, respectively, well above the critical value of 60° for complete wetting. The estimated percolation thresholds are at least 7% in volume, far exceeding the amount of metal in the mantle. Consequently, slab-derived Fe–C melt in the mid-mantle is expected to occur as isolated pockets and would not percolate through its silicate matrix.
Introduction
Iron–carbon alloys represent reduced forms of carbon and are important components of Earth’s long-term carbon cycle. Metallic iron may have been actively produced in the lower mantle (e.g., ), the mantle transition zone, and deep part of the upper mantle (e.g., ), due to disproportionation of ferrous iron (Fe2+) in silicates into iron metal (Fe0) and ferric iron (Fe3+). At depths greater than ∼250 km, metallic iron may reduce subducted carbonates to produce elemental carbon or carbide (MgCO3 + 2 Fe0 = 3 (Fe, Mg)O + C, e.g., ). Studies of diamond inclusions show that metallic melt containing Fe and C is present in the deep mantle (e.g., ; ). Depending on its wetting behavior, such melt may be trapped locally or drain into greater depths. Furthermore, carbon is a candidate for the light alloying element for the Earth’s iron-rich core, which likely constitutes the predominant carbon reservoir of the planet (e.g., Wood et al., 2013; ). If the Fe–C melt could form an interconnected network, it may eventually segregate to the core-mantle boundary (CMB), possibly merging into the liquid outer core. It has been proposed that slab-derived Fe–C melt in isolated patches can match seismically observed density and velocity features of some ultralow velocity zones (ULVZs) (). For these reasons, knowledge of the wetting behavior of Fe–C is crucial for assessing the fate of slab-derived metallic liquid in deep Earth.
The wetting behavior of iron-rich melt containing sulfur, oxygen, and silicon at high pressures has been studied previously in order to understand the mechanisms of core formation (e.g., ; ; ). However, that of Fe–C melt in silicate matrix remains poorly understood. Experimental data on the wetting behavior of Fe–C melt at mid-mantle conditions is scarce and limited to <6 GPa (; Zhimulev et al., 2018). In this study, we conduct experiments to investigate the wetting behavior of Fe–C melt in a matrix of polycrystalline olivine polymorphs and bridgmanite/ferropericlase at the conditions of Earth’s mid-mantle between 10 and 23 GPa and 1600 and 1800°C. The results are applied to assess the feasibility of percolative drainage of iron–carbon melt through the mantle.
Materials and Methods
Experimental
High pressure experiments were carried out using the 1000-ton Walker-type multi-anvil press at the University of Michigan. Experiments at pressures up to 17 GPa used Toshiba-Tungaloy F-grade tungsten carbide (WC) cubes with 5 mm truncation edge length (TEL) and those at pressures up to 23 GPa used similar WC cubes with 3 mm TEL (Figure 1). Temperature was monitored using the power curves determined by previous multi-anvil experiments with type-C (W95Re5-W74Re26) thermocouples. The pressure uncertainty is estimated to be ±2 GPa and the temperature uncertainty is estimated to be ±100 K, on the basis of duplication experiments ().
FIGURE 1
The starting material of the experiments consists of a powder mixture of San Carlos olivine, iron, and graphite. The olivine powder was prepared by grinding natural grains under ethanol in an agate mortar for 1 h and then dried at 110°C. Its composition is approximately (Mg0.9Fe0.1)2SiO4 with 2900 ppm Ni (). An Fe–C mixture with 96 wt% Fe and 4 wt% graphite, corresponding to the eutectic composition of the Fe–C binary system at 10 GPa () was homogenized and then mixed with the olivine power with a ratio of approximately 97:3 in volume. The volume fraction of Fe–C in the starting material was increased above 7% in one run for comparison. The sample was packed into a MgO capsule, which was dried at 950°C for approximately 2 h before loading to minimize trapped moisture. COMPRES 10/5 and 8/3 assemblies were used for experiments up to 17 GPa and those above 20 GPa, respectively ().
In each experiment, the sample was first pressurized to 5 GPa and then sintered at 800–1000°C for 8–12 h during further compression to the target pressure. Target temperature was set below the melting temperature of the silicate phase and above the melting temperature of the iron-rich mixture. The sample was held at target pressure and temperature for approximately 30 min and then quenched by cutting power to the heater and slowly decompressed to ambient pressure.
Analytical
Recovered experimental products were examined for texture, structure, and composition using optical microscopes, micro-Raman spectroscope, and scanning electron microscope (SEM) with energy dispersive spectrometer (EDS). The sample was first mounted in epoxy, sectioned and polished for analyses using optical microscopes. Raman analysis was then performed using a Renishaw inVia confocal Raman Microscope with a 532 nm continuous wave laser. The input laser power varied between 4 and 40 mW and the laser beam was focused through a 50× Mitutoyo objective lens to a diameter of 3–5 microns or through a 20× Olympus objective lens to a beam diameter of 7–12 microns. Raman signal was dispersed with an 1800 grooves/cm grating to achieve spectral resolution of 1–2 wave number (wn). The spectral window was centered at 520 wn. The signal was recorded with a Peltier cooled CCD. Acquisition time is typically 1 s per CCD window. Replicate measurements were conducted to evaluate heterogeneity. For composition and further texture analyses, the polished sample was aluminum-coated and examined using a JOEL 700FLV field emission SEM in the Central Campus Electron Microbeam Analysis Laboratory (EMAL) at the University of Michigan.
Apparent dihedral angles were measured from back scattered electron (BSE) images of the recovered samples by using the thresholding and protractor functions of the image processing software Fiji (). To establish a consistent criterion and hence to increase reproducibility of dihedral angle measurement, the BSE images of the measured areas were converted into binary format. Each individual dihedral angle was measured by the “protractor” function in Fiji. The actual fraction of Fe–C in each experiment was estimated from the polished cross-sections of recovered samples using the “Analyze Particles…” function in Fiji.
Results and Discussion
Four experiments at pressures between 10 and 23 GPa and temperatures at 1600°C or 1800°C produced domains of quenched Fe–C melt dispersed in solid matrix (Table 1 and Figure 2). The molten state of the Fe–C phase is inferred from the spherical shape and/or the dendritic texture of interspersed carbon-rich and carbon-poor domains. The solid state of the silicate or oxide is inferred from the euhedral or sub-hedral shapes. The coexistence of Fe–C melt with solid silicate and oxide at the experimental pressures and temperatures are consistent with known melting behavior of silicates and iron–carbon binary system (Figure 3).
TABLE 1
| Exp. # | Pressure (GPa) | Temperature (°C) | Melt fraction (%) | Matrix | # of angles measured | Dihedral angles (°) |
| M053016 | 10 | 1600 | 2.6 | Olivine | 208 | 117 ± 14 |
| M052416 | 17 | 1600 | 0.5 | Ringwoodite | 91 | 120 ± 14 |
| M013017 | 23 | 1800 | 1.2 | Bridgmanite/Ferropericlase | 257 | 107 ± 16 |
| M060316 | 23 | 1800 | 7.3 | Bridgmanite/Ferropericlase | – | – |
Summary of experimental conditions and results.
Pressure uncertainty is ±2 GPa and temperature uncertainty is ±100 K. In M060316, no clear triple junction is available for dihedral measurement. Melt fraction is estimated from a cross-section of the recovered sample for each run.
FIGURE 2
FIGURE 3
In all experiments, the Fe–C melt after quench formed isolated pockets at the grain boundaries of silicate or silicate and oxide matrix. The volume fraction of the Fe–C melt ranged from 0.5 to 7.3%. At small melt fractions between 0.5 and 2.6%, the Fe–C melt assumed nearly spherical shape with a radius of R = 0.2–1 micron (R, radius of an equal area circle) and showed no connectivity at all (Figure 2, upper). The Fe–C sphere appear somewhat distorted at 17 and 23 GPa. Although some large melt domains (R ≥ 5 micron) were formed at a melt fraction of 7.3%, the dihedral angles of the small isolated melt pockets (R ≤ 1 micron) remained large and full connectivity was not developed (Figure 2, lower). All the textures suggest that Fe–C melt does not wet silicate or silicate plus oxide at the conditions of the mid-mantle.
Raman analyses showed that the olivine converted into ringwoodite at 17 GPa, and dissociated into bridgmanite and ferropericlase at 23 GPa (Figure 4). The carbon contents of the Fe–C melt phases were only estimated semi-quantitatively from EDS data because the domains are typically sub-micron in size. They are consistent with the expected 4 ± 2 wt%. The Fe–C melt phases sometimes contained a small amount of nickel, likely originating from San Carlos olivine in the starting materials.
FIGURE 4
Dihedral Angle of Fe–C Melt in Mantle Silicate
The apparent dihedral angles of Fe–C melt in silicate matrix at mid-mantle conditions are measured in cross-sections of the experimental products. Dihedral angle is defined as the angle between the junction of two solid grains and melt pocket in the plane normal to the junction of the two solid grains (e.g., ). Previous studies suggested that a large number of measurements (50–100) are adequate to determine the effective dihedral angles reliably (). Here, between 91 and 257 measurements yielded distributions close to Gaussian (Table 1 and Figure 5). The mean value of the apparent dihedral angle distribution is taken as the effective dihedral angle, with one standard deviation as the associated uncertainty. The dihedral angles of Fe–C melt in olivine, ringwoodite, and bridgmenite-ferropericlase are 117 ± 14° at 10 GPa, 120 ± 14° at 17 GPa, and 107 ± 16° at 23 GPa, respectively. In previous studies, the median is sometimes taken as the effective dihedral angle of the system when the distribution is not near Gaussian. For comparison, we also report the medians of the measured dihedral angles in each experiment with a 95% confidence interval around the medians (), which are 116+4–2° at 10 GPa, 119+4–4° at 17 GPa, and 108+2–2° at 23 GPa, respectively.
FIGURE 5
Textural Equilibrium
The measured dihedral angle is accurate only when the system has reached textural equilibrium so that with time the microstructures will only change in size but not in type or shape. At textural equilibrium, the distribution of the apparent dihedral angels formed by the liquid phase assumes a Gaussian pattern, with a single peak that drops off rapidly on both sides (). Moreover, the distribution of grain sizes approaches the theoretically expected pattern, and the apparent angles formed by the solid grains at triple junction center around 120°, also dropping off rapidly above or below the mean value. Assuming that the effect of solid silicate matrixes is negligible, the dihedral angle measured in an equilibrated phase assembly is characteristic of the composition at specific pressure and temperature conditions.
Time series studies showed that iron-rich melt in solid silicate matrix reached steady-state texture within 12 min at 25 GPa and 2500 K () in laser-heated diamond anvil cells, and within 50 min at 14 GPa and 1700°C in multi-anvil apparatus (). In this study, the experiments were sintered at 5 GPa and 800–1000°C for 8–12 h and then held at 10–23 GPa and 1600–1800°C for about 30 min. The pressure and temperature conditions and experimental durations are comparable to the studies with demonstrated textural equilibrium. Moreover, the distributions of the apparent angles can be fitted by Gaussian curves (Figure 5), indicating sufficient approach to texture equilibrium. Equilibrium is further supported by the nearly 120° of dihedral angles at triple junctions of solid grains (Figure 2). Although the equilibrated triple junctions of solid grains may have formed at the pre-annealing or sintering stage, textural re-equilibration between solid grains at a slightly higher temperature would take shorter time than starting from ambient temperature, and the Fe–C melt is expected to establish textural equilibrium with the solid phases much faster than within the solid grains.
Interfacial Energy and Surface Tension Between Fe–C Melt and Mantle Silicate
The measured dihedral angles allow us to estimate the interfacial energy and surface tension between Fe–C melt and various mantle silicate phases at high pressure and high temperature conditions. The dihedral angle of a liquid in a solid matrix is governed by the interfacial energy per unit area between the liquid and the solid, γsl, and that between the solid phases, γss. In an isotropic system, force balance at the triple junction can be described by the Cahn’s equation (e.g., ):
The interfacial energy between solid silicate grains (olivine) is estimated at about 0.9 J/m2 with a relatively large uncertainty of ±30% (). If this value is applicable to silicates at high pressures and high temperatures, then the measured dihedral angles of 107–120° yield the interfacial energy between Fe–C melt and major mantle silicates at approximately 0.76–0.9 J/m2 at 10–23 GPa. The surface tension of Fe–C melt on mantle silicate, which has the identical numerical value to the interfacial energy between liquid and solid (), is 0.76–0.9 N/m.
The relation 2⋅cos θ/2 = γss/γsl is derived for simple systems containing a single isotropic solid phase, and therefore it does not hold accurately for anisotropic mantle silicates such as olivine and bridgmanite. showed that with sufficiently high anisotropy the melt may form isolated plane-faced pockets at grain corners. A few experimental studies examined the anisotropy of γsl at olivine-basalt interface at low pressures (e.g., ), but the influence of anisotropy on the wetting behavior of iron-rich melt system remains to be explored.
According to the relationship between dihedral angle and interfacial energies, θ = 120° when γss = γsl. Below 120°, θ decreases rapidly with the ratio and reaches the critical value of 60° when or γsl∼0.6γss. This relationship can be applied to understand the effects of composition, pressure, temperature on the wetting properties of Fe–C melt, as will be discussed below.
Composition Effect on the Dihedral Angles of Iron-Rich Melt
Our results show that the measured dihedral angles of Fe–C melt in silicates remained near 120°, indicating that the Fe–C melt does not wet major silicate phases in the mantle. The surface-avoiding or surface-phobic behavior of carbon is similar to that observed in metallurgical studies at ambient pressure (; ). It has been shown that at ambient pressure adding carbon to an iron liquid has little effect on the surface tension of liquid iron, which remains high even with more than 10 wt% of carbon (Figure 6A). In contrast, metallurgical studies showed that the surface tension of iron-rich liquid on aluminum oxide is sensitive to the presence of alloying elements such as O and S. At 1 bar and 1843 K, adding 0.1 wt% oxygen rapidly reduces the surface tension of liquid iron on alumina to by about 35%, whereas adding the same amount of sulfur to liquid iron reduces the surface tension by 30% (Figure 6A). These elements are therefore known as surface-loving elements or surfactants.
FIGURE 6
At high pressures, the dihedral angle of iron-rich melt in silicate is also influenced by the presence of alloying elements such as O and S (Figure 6B). For example, adding oxygen to an Fe–S melt lowers its dihedral angle in silicate matrix, which could fall below the critical value of 60° at sufficiently high oxygen concentrations (, , ). Other non-metal elements, such as S (e.g., ), Si (e.g., ), and P (e.g., ) can also reduce the dihedral angle of iron-rich melt in silicate, although the effect is weaker than oxygen.
Previous studies found that oxidized starting material or moistures trapped in sample capsule may introduce oxygen to the melt, as revealed by the presence of tiny FeO blobs exsolved from iron-rich melt upon quench (e.g., ). Carbon and oxygen are likely incompatible in iron-rich melt (). In this study, no tiny blob was observed in the quenched Fe–C melt, suggesting little influence of oxygen on the measured dihedral angles.
Effect of Pressure on the Dihedral Angle of Fe–C Melt in Silicate
The results from this study show little pressure dependence of the dihedral angle of Fe–C melt in mantle silicate between 10 and 23 GPa, with the variations from 107 to 120° falling within the measurement uncertainty of ±14–16° (Figure 7). The measured angles are also comparable to the dihedral angle of 115° between Fe–C melt with 4.3 wt% carbon and olivine at 3 GPa and 1400°C ().
FIGURE 7
Pressure can influence the dihedral angle of Fe-rich melt in silicate in a number of ways (). Misalignment of crystal lattices produces excess volume or gap, which contributes to interfacial energy γss. The excess volume between melt and crystal is likely smaller because atoms in the melt are not constrained by a lattice and hence moves more freely to fill space. Upon compression, the amount of lattice misalignment may change as a result of pressure-induced preferred orientation, polymorphic structural transformation, or even chemical reaction such as the breakdown of ringwoodite into bridgmanite and ferropericlase. Compression likely reduces the difference between the two excess volumes, leading to a smaller γss/γsl, hence larger dihedral angle.
Compression may also affect the chemical properties of the interface such as the partitioning of surface-active elements between the coexisting melt and solid. For instance, significant chemical change such as enhanced solubility of oxygen in iron-rich alloy may have reduced γsl, thus lowering the dihedral angle (). Even if the solubilities of surface-active elements remain constant, their surface activities may change with pressure (). The same amount of sulfur affects dihedral angle less significantly at high pressure () than at 1 bar (). Furthermore, segregation of certain species into the grain boundary may occur as a result of compression (). These changes could affect γss and γsl and thus the dihedral angle.
The absence of pressure effect on the dihedral angle of Fe–C melt in mantle silicates suggests that ratio of interfacial energies remain nearly unchanged between 3 and 23 GPa, unless various pressure-induced effects on the ratio fortuitously offset each other. Within the pressure range of the upper mantle and transition zone, the dihedral angle of Fe–S melt with 12–14 wt% S in silicate also falls between 100 and 120°, showing little pressure dependence (Figure 7; ). For Fe–Si melt in silicate, available data at pressure below 5 GPa do not show any pressure dependence either (Figure 7; e.g., ).
At higher pressures, diamond-anvil cell experiments on Fe–Si melt up to 47 GPa and 3000 K () and on Fe–S melt up to 64 GPa and 3300 K () showed that the dihedral angles of molten iron-alloys in bridgmanite matrix decreased rapidly with increasing pressure and temperature to below the critical value of 60° (Figure 7). The large negative pressure effect may be related to large deviatoric stress in solid phases, which would increase γss significantly and therefore reducing θ (Figure 5; ). Alternatively, significant chemical change such as enhanced solubility of oxygen in the iron-rich alloy may have reduced γsl, thus reducing the dihedral angle (). It is also possible that temperature effect plays an important in the observed reduction of dihedral angle, as discussed below.
Effect of Temperature on the Wetting Behavior of Fe–C Melt in Mantle Silicate
No temperature effect on the dihedral angle of Fe–C melt in mantle silicate was observed between 1600 and 1800°C, at 10–23 GPa. Strong temperature dependence of surface tension was reported in metallurgical studies on iron alloys at 1 bar (e.g., ). At 25 GPa, the dihedral angle of Fe–Si melt in bridgmanite decreases from 94° at 2400 K to 70° at 2800 K (), also suggesting potentially significant temperature effect on the surface tension of compressed iron-rich melt, especially at temperatures approaching the melting point of the solid phase. In this study, the experimental temperatures were near the eutectic temperatures of Fe–C binary systems and far below the melting point of the silicate matrixes. The relatively low temperatures and small temperature range may explain the absence of temperature effect in our data.
At deep Earth conditions, it is often difficult to separate the effects of pressure, temperature, and composition on the dihedral angle. The temperature of a dihedral angle experiment must fall above the liquidus of iron-rich alloy and below the solidus of silicates, leaving a relatively narrow temperature range at which iron-rich melt coexists with solid silicate. This range usually shifts upward as pressure increases (Figure 3), and therefore it may be impossible to isolate the pressure and temperature effects. For Fe–S, the pressure rise from 25 to 64 GPa is accompanied by an increase of temperature from 2300 to 3300 K (); thus it remains unclear if temperature, pressure, or both are responsible for the decrease in the dihedral angle of Fe–S in silicates.
Effect of Crystalline Structure on the Dihedral Angle of Fe–C Melt in Silicate
Our results show that the dihedral angle of Fe–C melt is not affected by pressure-induced polymorphic transition or decomposition of silicate at mid-mantle conditions. In the experiments, olivine powder grew in crystal size at 10 GPa, and transformed into ringwoodite at 17 GPa. At 23 GPa, ringwoodite broke down to form bridgmanite and ferropericlase. The measured dihedral angles of Fe–C melt remained near 115°. In contrast, the wetting behavior of Fe–S changed abruptly upon the break-down of ringwoodite, with the dihedral angle of Fe–S melt with 12–14 wt% S dropping from approximately 108° at upper mantle conditions to 71° at 25 GPa (). It is possible that the presence of ferropericlase increased the solubility of oxygen in the Fe–S melt, thus lowering its surface tension, whereas oxygen remains highly insoluble in the Fe–C melt, which is more reducing than the Fe–S melt.
Melt Connectivity and Percolation Threshold
The dihedral angles of Fe–C melt in silicate matrix at 3–23 GPa are well above the critical value of 60°, implying that Fe–C melt can only flow through mantle silicate when its volume fraction exceeds ∼7%. Dihedral angle measures the ability of liquid to wet a solid matrix. A dihedral angle of 60° or smaller implies complete wetting, where an interconnected network of liquid forms regardless of the melt fraction. A dihedral angle that is greater than 60° implies that a fraction of melt may be stranded in the solid matrix as isolated pockets. The minimum volume fraction of liquid to form a network that enables passage is known as the percolation threshold, fp. For homogeneously distributed melt in solid matrix, fp can be estimated from θ through the following relationship:
With θ∼ 120°, the percolation threshold is estimated at ∼7 vol% (e.g., ; ). It should be noted that dihedral angle is a function of interfacial energy ratio (Eq. 1) while connectivity can be affected by both dihedral angle and melt fraction. Full connectivity of Fe-rich melt in silicate will develop at a wetting angle of >60° if the melt fraction is sufficiently large. However, at very low melt fractions, as relevant to the amount of metallic melt in the Earth’s mantle, percolative drainage (forming an interconnected network) can only occur when the wetting angle of Fe-rich melt is below 60°.
Below the threshold, percolation of iron-rich melt through silicate may still occur if the silicate matrix becomes partially molten (e.g., ) or deforms under shear stress (e.g., ). Although the presence of partial melt of silicate does not change the dihedral angles of iron-rich melt (e.g., ), droplets of dense iron-rich melt can sink through molten silicate, especially when the amount of silicate melt is large enough to create channels (). The presence of non-wetting silicate melt may also raise the percolation threshold of iron-rich melt (Yoshino et al., 2004).
A number of studies have investigated if shear deformation can enhance connectivity of iron-rich melt and facilitate its percolation through silicate, but the effect remains poorly constrained. Walte et al. (2011) suggested that deformation only promote Fe–S melt segregation at high strain rates, with no influence at geologically relevant strain rates. also showed that the ability of Fe–S melt to flow through olivine matrix did not change with the amount of shear deformation. further suggested that the onset of silicate melting nullifies any effect of deformation on melt segregation and may hinter segregation of Fe-rich melt by reducing its connectivity. In contrast, suggested that deformation aided percolation at high pressure over a very large range of strain rates.
Implications for Earth’s Deep Carbon Cycle
The measured large dihedral angles and estimated high percolation thresholds imply that percolation may not be a viable mechanism to segregate Fe–C melt from the mantle into the core. The amount of Fe–C melt in the mid-mantle was estimated at no more than 1 wt% (). This is well below the percolation threshold of ∼7 vol%, and therefore the Fe–C melt is expected to be stranded as isolated pockets at the mid-mantle pressures.
The presence of other alloying elements such as sulfur, hydrogen, and/or oxygen may allow Fe–C melt to percolate through surrounding silicates toward greater depths. Metallic phases in diamond inclusions consist not only Fe and C but also S and H (). Adding sulfur may reduce the interfacial energy of Fe–C melt and lower the percolation threshold. The effect of hydrogen is not known. While oxygen is effective in lowering the surface tension of iron-rich melt, its solubility in Fe–C melt is likely negligible because of the reducing condition.
Partial melt of silicates may also facilitate percolation of iron-rich melt (). However, the present-day mantle temperature is too low to cause extensive melting of the silicates, thus segregation of Fe–C melt through channels of molten silicate is unlikely.
Given that the percolation threshold of Fe–C binary melt is much higher than the amount of Fe–C melt in the mid-mantle, and that partial melting of silicate matrix is unlikely with the present-day geotherm, we conclude that the slab-derived Fe–C melt is most likely to be retained within the subducting slab at mid-mantle depths during diamond growth, unless its wetting behavior is sufficiently modified by the presence of other alloying elements to allow percolation.
Although slab-derived Fe–C melt may not be able to percolate through the mantle, it may still reach the CMB if the sinking slab materials can survive viscous stirring associated with mantle plumes and remain largely unmixed with the ambient mantle. proposed that near the base of the mantle, the presence of slab-derived iron–carbon melt in isolated patches may explain seismically observed density excess and velocity features of some ULVZs. Here, we found that the dihedral angle of Fe–C in mid-mantle is larger than 100°. If applicable at the conditions near the CMB, then more than 10 vol% Fe–C is needed to match the ULVZ observations (). This amount is, however, larger than the percolation threshold, implying that not enough Fe–C can be retained in silicate matrix to produce long-lived ULVZs, unless its drainage is hindered by the presence of silicate partial melt (). Alternatively, metallic melt may be replenished through on-going subduction or infiltration of core melt back into ferropericlase matrix due to morphological instability ().
Once at the CMB, Fe–C melt may segregate into the core if the high temperature at the base of the lower mantle lowers its dihedral angle and percolation threshold sufficiently. With an estimated 1 wt% Fe–C in the mantle (e.g., Wood et al., 2013), loss of Fe–C from 40% of the mantle corresponds to 0.8% growth in core mass at the expense of 0.4% reduction in mantle mass, which is equivalent to about 10-km expansion in core radius over Earth’s history. This is insignificant for the iron-rich core but would be a major leak for carbon in the silicate Earth.
Conclusion
The wetting behavior of Fe–C melt in silicates at mid-mantle conditions have been studied using multi-anvil experiments. The dihedral angles of Fe–C melt in olivine, ringwoodite, and bridgmanite/ferropericlase matrixes are measured at 117° at 10 GPa, 120 ± 14° at 17 GPa, and 107 ± 16° at 23 GPa, respectively. Carbon is found to behave as a surface-avoiding element in silicates at pressures up to 23 GPa, consistent with the high surface tension of Fe–C melt on alumina at 1 bar. Within the experimental ranges, no effect of pressure or temperature was observed.
The measure dihedral angles and estimated percolation thresholds are applied to assess the behavior of Fe–C alloys at deep Earth conditions. The percolation threshold of Fe–C melt in silicates is estimated at 7 vol%, much larger than the amount produced in the mantle. Consequently, Fe–C melt is expected to be retained within subducting slabs at mid-mantle depths (Figure 8), unless its surface tension is sufficiently lowered by the presence of other alloying elements such as sulfur. The amount of Fe–C melt needed to match the low shear wave velocity of the ULVZs, however, exceeds the percolation threshold. As a result, continuous replenishment, dynamic stirring, or the presence of non-wetting silicate melt is needed to generate ULVZ through slab-derived metallic melt. The non-wetting behavior of Fe–C implies that active core growth through percolation is unlikely, even though limited draining of Fe–C to the core remains feasible and may cause significant carbon loss from the silicate Earth.
FIGURE 8
Further investigations are necessary to fully assess the wetting behavior of Fe–C melt during core formation and deep carbon cycle. This study only reached the pressure at the uppermost lower mantle. The wetting behavior of Fe–C melt in main stretch of the lower mantle remains to be studied. Additional work is also needed to demonstrate equilibrium texture through time-series experiments, quantify the composition of the Fe–C melt, evaluate the effect of temperature, and to explore the influence of anisotropy on the wetting behavior of iron-rich melts.
Statements
Data availability statement
All datasets generated for this study are included in the manuscript/supplementary files.
Author contributions
JL and JD conceived the idea and designed the project. JD performed the experiments and analyzed the data, as part of his undergraduate honors thesis under the supervision of JL. JD and JL wrote the manuscript. FZ helped with the experiments.
Funding
This work was partially supported by the National Science Foundation Grants AST 1344133 and EAR 1763189, and the Alfred P. Sloan Foundation Deep Carbon Observatory Grant G-2017-9954 to JL; and a Turner Undergraduate Award from the Department of Earth and Environmental Sciences at University of Michigan, a James Mills Peirce Fellowship from the Graduate School of Arts and Sciences at Harvard University, open access publication funding from the Harvard Open-Access Publishing Equity (HOPE) office, and the Department of Earth and Planetary Sciences (EPS) at Harvard University to JD.
Acknowledgments
We thank the reviewers, GB and VC, for their reviews that helped to improve and clarify the manuscript; Hidenori Terasaki, Shun-ichiro Karato, Rebecca A. Fischer, Wendy L. Mao, Carl B. Agee, and David (Dave) Walker for discussions and comments; and Paul Kelley, the Harvard Open-Access Publishing Equity (HOPE) office, and the Frontiers Wavier Office for providing open-access funding opportunities.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
BergM. T.BromileyG. D.ButlerI. B.FrostM.BradleyR.CarrJ.et al (2017). Deformation-aided segregation of Fe-S liquid from olivine under deep earth conditions: implications for core formation in the early solar system.Phys. Earth Planet. Inter.26338–54. 10.1016/j.pepi.2017.01.004
2
CerantolaV.WalteN. P.RubieD. C. (2015). Deformation of a crystalline olivine aggregate containing two immiscible liquids: implications for early core–mantle differentiation.Earth Planet. Sci. Lett.41767–77. 10.1016/j.epsl.2015.02.014
3
DuncanM. S.FeiY. (2017). “Experimental constraints on metal percolation through silicate: implications for core formation on asteroids and planetesimals,” in Lunar and Planetary Science Conference, Vol. 48 (Houston, TX).
4
FeiY.BroshE. (2014). Experimental study and thermodynamic calculations of phase relations in the Fe–C system at high pressure.Earth Planet. Sci. Lett.408155–162. 10.1016/j.epsl.2014.09.044
5
FischerR. A.NakajimaY.CampbellA. J.FrostD. J.HarriesD.LangenhorstF.et al (2015). High pressure metal–silicate partitioning of Ni, Co, V, Cr, Si, and O.Geochim. Cosmochim. Acta167177–194. 10.1073/pnas.1108544108
6
FrostD. J.LiebskeC.LangenhorstF.McCammonC. A.TrønnesR. G.RubieD. C. (2004). Experimental evidence for the existence of iron-rich metal in the Earth’s lower mantle.Nature428409–412. 10.1038/nature02413
7
HaldenF. A.KingeryW. D. (1955). Surface tension at elevated temperatures. II. Effect of C, N, O and S on liquid iron surface tension and interfacial energy with Al2O3.J. Phys. Chem.59557–559. 10.1021/j150528a018
8
HerzbergC.ViditoC.StarkeyN. A. (2016). Nickel–cobalt contents of olivine record origins of mantle peridotite and related rocks.Am. Mineral.1011952–1966. 10.2138/am-2016-5538
9
HerzbergC.RaterronP.ZhangJ. (2000). New experimental observations on the anhydrous solidus for peridotite KLB−1.Geochem., Geophys., Geosyst.1:1051. 10.1029/2000GC000089
10
HirayamaY.FujiiT.KuritaK. (1993). The melting relation of the system, iron and carbon at high pressure and its bearing on the early stage of the Earth.Geophys. Res. Let.202095–2098. 10.1029/93gl02131
11
HolnessM. B.LewisS. (1997). The structure of the halite-brine interface inferred from pressure and temperature variations of equilibrium dihedral angles in the halite-H2O-CO2 system.Geochim. Cosmochim. Acta61795–804. 10.1016/s0016-7037(96)00370-5
12
HolzheidA. (2013). Sulphide melt distribution in partially molten silicate aggregates: implications to core formation scenarios in terrestrial planets.Eur. J. Mineral.25267–277. 10.1127/0935-1221/2013/0025-2264
13
HolzheidA.SchmitzM. D.GroveT. L. (2000). Textural equilibria of iron sulfide liquids in partly molten silicate aggregates and their relevance to core formation scenarios.J. Geophys. Res.: Solid Earth10513555–13567. 10.1029/2000jb900046
14
HultgrenR.OrrR. L.AndersonP. D.KelleyK. K. (1963). Selected Values of Thermodynamic Properties of Metals and Alloys.New York, NY: John Wiley and Sons, Inc, 107.
15
IidaT.GuthrieR. I. (2015). The Thermophysical Properties of Metallic Liquids: Fundamentals, Vol. 1. Oxford: Oxford University Press.
16
KaminskyF. V.WirthR. (2011). Iron carbide inclusions in lower-mantle diamond from Juina. Brazil.Can. Mineral.49555–572. 10.3749/canmin.49.2.555
17
LaporteD.WatsonE. B. (1995). Experimental and theoretical constraints on melt distribution in crustal sources: the effect of crystalline anisotropy on melt interconnectivity.Chem. Geol.124161–184. 10.1016/0009-2541(95)00052-n
18
LeinenweberK. D.TyburczyJ. A.SharpT. G.SoignardE.DiedrichT.PetuskeyW. B.et al (2012). Cell assemblies for reproducible multi-anvil experiments (the COMPRES assemblies).Am. Mineral.97353–368. 10.2138/am.2012.3844
19
LiJ.FeiY. (2014). Experimental constraints on core composition.Treatise Geochem.2568–546.
20
LiZ.LiJ. (2015). Melting curve of NaCl to 20 GPa from electrical measurements of capacitive current.Am. Mineral.1001892–1898. 10.2138/am-2015-5248
21
LiuJ.LiJ.HrubiakR.SmithJ. S. (2016). Origins of ultralow velocity zones through slab-derived metallic melt.Proc. Natl. Acad. Sci. U.S.A.1135547–5551. 10.1073/pnas.1519540113
22
LordO. T.WalterM. J.DasguptaR.WalkerD.ClarkS. M. (2009). Melting in the Fe–C system to 70 GPa.Earth Planet. Sci. Lett.284157–167. 10.1016/j.epsl.2009.04.017
23
MannU.FrostD. J.RubieD. C. (2008). The wetting ability of Si-bearing liquid Fe-alloys in a solid silicate matrix-percolation during core formation under reducing conditions?Phys. Earth Planet. Int.1671–7. 10.1016/j.pepi.2007.12.002
24
OtsukaK.KaratoS.-I. (2012). Deep penetration of molten iron into the mantle caused by a morphological instability.Nature492243–246. 10.1038/nature11663
25
PasseroneA.SangiorgiR. (1985). Solid-liquid interfacial tensions by the dihedral angle method. A mathematical approach.Acta Metallurgica33771–776. 10.1016/0001-6160(85)90100-2
26
PorterD. A.EasterlingK. E.SherifM. (2009). Phase Transformations in Metals and Alloys, (Revised Reprint).CRC press.
27
RohrbachA.BallhausC.Golla–SchindlerU.UlmerP.KamenetskyV. S.KuzminD. V. (2007). Metal saturation in the upper mantle.Nature449456–458. 10.1038/nature06183
28
RohrbachA.SchmidtM. W. (2011). Redox freezing and melting in the Earth’s deep mantle resulting from carbon–iron redox coupling.Nature472209–212. 10.1038/nature09899
29
RubieD. C.JacobsonS. A. (2016). Mechanisms and geochemical models of core formation.Deep Earth: Physics and Chemistry of the Lower Mantle and Core217181–190. 10.1002/9781118992487
30
SchäferF. N.FoleyS. F. (2002). The effect of crystal orientation on the wetting behaviour of silicate melts on the surfaces of spinel peridotite minerals.Contrib. Mineral. Petrol.143254–262. 10.1007/s00410-001-0339-0
31
SchindelinJ.Arganda-CarrerasI.FriseE.KaynigV.LongairM.PietzschT.et al (2012). Fiji: an open-source platform for biological-image analysis.Nat. Methods9676–682. 10.1038/nmeth.2019
32
ShannonM. C. (1998). Percolation of Iron Melts Through Mantle Minerals at High Pressure: Implications for Core Formation.Ph.D. Dissertation, Harvard University, Cambridge, MA.
33
ShannonM. C.AgeeC. B. (1996). High pressure constraints on percolative core formation.Geophys. Res. Lett.232717–2720. 10.1029/96gl02817
34
ShannonM. C.AgeeC. B. (1998). Percolation of core melts at lower mantle conditions.Science2801059–1061. 10.1126/science.280.5366.1059
35
ShiC. Y.ZhangL.YangW.LiuY.WangJ.MengY.et al (2013). Formation of an interconnected network of iron melt at Earth’s lower mantle conditions.Nat. Geosci.6971. 10.1038/ngeo1956
36
SmithE. M.ShireyS. B.NestolaF.BullockE. S.WangJ.RichardsonS. H.et al (2016). Large gem diamonds from metallic liquid in Earth’s deep mantle.Science3541403–1405. 10.1126/science.aal1303
37
StevensonD. J. (1990). “Fluid Dynamics of Core Formation,” in Origin of the Earth.Oxford: Oxford University Press, 231–249.
38
StickelsC. A.HuckeE. E. (1964). Measurement of dihedral angles.Trans. Metallurgical Soc. AIME230:795.
39
StixrudeL.Lithgow-BertelloniC. (2011). Thermodynamics of mantle minerals-II. Phase equilibria.Geophys. J. Int.1841180–1213. 10.1111/j.1365-246x.2010.04890.x
40
TakafujiN.HiroseK.OnoS.XuF.MitomeM.BandoY. (2004). Segregation of core melts by permeable flow in the lower mantle.Earth Planet. Sci. Lett.224249–257. 10.1016/j.epsl.2004.05.016
41
TerasakiH.FrostD. J.RubieD. C.LangenhorstF. (2005). The effect of oxygen and sulphur on the dihedral angle between Fe–O–S melt and silicate minerals at high pressure: implications for Martian core formation.Earth Planet. Sci. Lett.232379–392. 10.1016/j.epsl.2005.01.030
42
TerasakiH.FrostD. J.RubieD. C.LangenhorstF. (2007). Interconnectivity of Fe–O–S liquid in polycrystalline silicate perovskite at lower mantle conditions.Phys. Earth Planet. Int.161170–176. 10.1016/j.pepi.2007.01.011
43
TerasakiH.FrostD. J.RubieD. C.LangenhorstF. (2008). Percolative core formation in planetesimals.Earth Planet. Sci. Lett.273132–137. 10.1073/pnas.1707580114
44
ToddK. A.WatsonH. C.YuT.WangY. (2016). The effects of shear deformation on planetesimal core segregation: results from in-situ X-ray micro-tomography.Am. Mineral.1011996–2004. 10.2138/am-2016-5474
45
von BargenN.WaffH. S. (1986). Permeabilities, interfacial areas and curvatures of partially molten systems: results of numerical computations of equilibrium microstructures.Journal of Geophysical Research: Solid Earth, 919261–9276. 10.1029/JB091iB09p09261
46
WalteN. P.RubieD. C.BonsP. D.FrostD. J. (2011). Deformation of a crystalline aggregate with a small percentage of high-dihedral-angle liquid: implications for core–mantle differentiation during planetary formation.Earth Planet. Sci. Lett.305124–134. 10.1016/j.epsl.2011.02.049
47
WilliamsQ.JeanlozR.BassJ.SvendsenB.AhrensT. J. (1987). The melting curve of iron to 250 gigapascals: a constraint on the temperature at Earth’s center.Science236181–182. 10.1126/science.236.4798.181
48
WoodB. J.LiJ.ShaharA. (2013). Carbon in the core: its influence on the properties of core and mantle.Rev. Mineral. Geochem.75231–250. 10.1126/sciadv.1701840
49
YoshinoT.WalterM. J.KatsuraT. (2004). Connectivity of molten Fe alloy in peridotite based on in situ electrical conductivity measurements: implications for core formation in terrestrial planets.Earth Planet. Sci. Lett.222625–643. 10.1016/s0012-821x(04)00182-7
50
ZhimulevE. I.ChepurovA. I.SoninV. M.LitasovK. D.ChepurovA. A. (2018). Experimental modeling of percolation of molten iron through polycrystalline olivine matrix at 2.0–5.5 GPa and 1600°C.High Press. Res.38153–164. 10.1080/08957959.2018.1458847
Summary
Keywords
dihedral angle, metallic melt, percolation, ultra-low velocity zone, core growth
Citation
Dong J, Li J and Zhu F (2019) Wetting Behavior of Iron–Carbon Melt in Silicates at Mid-Mantle Pressures With Implications for Earth’s Deep Carbon Cycle. Front. Earth Sci. 7:268. doi: 10.3389/feart.2019.00268
Received
31 May 2019
Accepted
30 September 2019
Published
22 October 2019
Volume
7 - 2019
Edited by
Sami Mikhail, University of St Andrews, United Kingdom
Reviewed by
Geoffrey David Bromiley, University of Edinburgh, United Kingdom; Valerio Cerantola, European X-Ray Free Electron Laser, Germany
Updates
Copyright
© 2019 Dong, Li and Zhu.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jie Li, jackieli@umich.edu
†Present address: Feng Zhu, Hawai‘i Institute of Geophysics and Planetology, University of Hawai‘i at Mānoa, Honolulu, HI, United States
This article was submitted to Earth and Planetary Materials, a section of the journal Frontiers in Earth Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.