Applications of liquid state physics to the Earth's core

By use of the modern theory of liquids and some guidance from the hard-sphere model of liquid structure, the following new results have been derived for application to the Earth's outer core. (1) dK/$\text{d}\text{P ? 5 ? 5. 6P}$/K, where K is the incompressibility and P the pressure. This is valid for a high-pressure liquid near its melting point, provided that the pressure is derived primarily from a strongly repulsive pair potential φ. This result is consistent with seismic data, except possibly in the lowermost region of the outer core, and demonstrates the approximate universality of dK/dP proposed by Birch (1939) and Bullen (1949). (2) dlnT_M/dlnρ = (γC_V ? 1)/($\text{C}{}_{\mathrm{<mtext>V</mtext>}}\text{?}\text{3}\text{2}\text{),}\text{where}\text{T}{}_{\mathrm{<mtext>M</mtext>}}$ is the melting point, ρ the density, γ the atomic thermodynamic Grüneisen parameter and C_V the atomic contribution to the specific heat in units of Boltzmann's constant per atom. This reduces to Lindemann's law for C_V = 3 and provides further support for the approximate validity of this law. (3) It follows that the “core paradox” of Higgins and Kennedy can only occur if $\text{γ <}\text{2}\text{3}$. However, it is shown that $\text{γ <}\text{2}\text{3}\text{? ∫}{}^{\mathrm{\infty }}{}_0\text{(?g/?T)}{}_{\mathrm{\rho }}\text{r(}\text{d}\text{/}\text{d}\text{r)(r}{}^2\text{φ)}\text{d}\text{r > 0}$, which cannot be achieved for any strongly repulsive pair potential φ and the corresponding pair distribution function g. It is concluded that $\text{γ >}\text{2}\text{3}$ and that the core paradox is almost certainly impossible for any conceivable core composition. Approximate calculations suggest that γ ～ 1.3–1.5 in the core. Further work on the thermodynamics of the liquid core must await development of a physically realistic pair potential, since existing pair potentials may be unsatisfactory.     

Thermal equation of state of magnesiowüstite (Mg0.6Fe0.4)O

Volume measurements for magnesiowüstite (Mg_0.6Fe_0.4)O, were carried out up to pressures of 10.1 GPa in the temperature range 300–1273 K, using energy-dispersive synchrotron X-ray diffraction. These data allow reliable determination of the temperature dependence of the bulk modulus and good constraint on the thermal expansitivity at ambient pressure which was previously not known for magnesiowüstite. From these data, thermal and elastic parameters were derived from various approaches based on the Birch–Murnaghan equation of state (EOS) and on the relevant thermodynamic relations. The results from three different equations of state are remarkably consistent. With (∂K_T/∂P)_T fixed at 4, we obtained K₀=158(2) GPa, (∂K_T/∂T)_P=−0.029(3) GPa K⁻¹, (∂K_T/∂T)_V=−3.9(±2.3)×10⁻³ GPa K⁻¹, and α_T=3.45(18)×10⁻⁵+1.14(28)×10⁻⁸T. The K₀, (∂K_T/∂T)_P, and (∂K_T/∂T)_V values are in agreement with those of Fei et al. (1992) and are similar to previously determined values for MgO. The zero pressure thermal expansitivity of (Mg_0.6Fe_0.4)O is found to be similar to that for MgO (Suzuki, 1975). These results indicate that, for the compositional range x=0–0.4 in (Mg_1−xFe_x)O, the thermal and elastic properties of magnesiowüstite exhibit a dependence on the iron content that is negligibly small, within uncertainties of the experiments. They are consequently insensitive to the Fe–Mg partitioning between (Mg, Fe)SiO₃ perovskite and magnesiowüstite when applied to compositional models of the lower mantle. With the assumption that (Mg_0.6Fe_0.4)O is a Debye-like solid, a modified equation of heat capacity at constant pressure is proposed and thermodynamic properties of geophysically importance are calculated and tabulated at high temperatures.     

A universal thermal equation-of-state

Evidence has steadily accumulated to show that at high temperatures (above the Debye temperature, θ) the thermal pressure, P_TH, of solids, is linear with T to a close approximation. This empirical finding yields a simple relationship between P, V, and T quite useful for the computation of the equation-of-state (EOS). For geophysical applications, the empirical data is, so far, limited to a few minerals, all of which are important to our geophysical models of the Earth. The same results have been found for a variety of types of solids, including alkali metals, noble gas solids, alkali halides and metals in addition to minerals. It is argued that the linearity between P_TH and T is a general high-temperature property of solids. This includes minerals. Thus it is proposed that there exists a common thermal EOS which transcends the chemical bonding type and crystallographic class.     

Are anharmonicity corrections needed for temperature-profile calculations of interiors of terrestrial planets?

It is shown that there is linearity between the thermal pressure P_TH and T between the Debye temperature θ and some high temperature $\text{T}{}^1$. $\text{T}{}^1$ has been measured at 1 atm and is reported for several minerals including, for example, MgO (1300 K) and forsterite (1200 K). The change in thermal pressure from room temperature for five solids, so far measured, indicate striking linearity with T at high temperatures.It is further shown that the value of $\text{T}{}^1$ increases greatly as the pressure increases. It is therefore concluded that P_TH is probably linear with T for mantle minerals under mantle conditions. The proportionality constant is derived from the measurements of thermal expansivity and bulk modulus at high temperature and zero pressure.The argument is then reversed. Assuming that the thermal pressure is in fact linear with T for the various shells in a planet, the resulting density and temperature profile of the planet is derived. The resulting density profile of the Earth compares favorably with corresponding values of recent seismic profiles.     

The high-temperature acoustic Grüneisen parameter in the earth's interior

The value of the acoustic Grüneisen parameter, γ_a, in the earth's interior has been computed using data from recent models obtained by inversion of normal data. In this paper we emphasize the data from the PEM model of the earth because there has been sufficient smoothing of the seismic data so that the derivatives d ln ν_s/d ? and d ln ν_p/d ? can be well defined at all depths.The results for the lower mantle show that γ_a decreases exponentially from 1.3 to 1.0, and there are several consistent cross-checks of the limiting values. We find γ_a is about 1.5 for the inner core and outer core. These results confirm, in broad outline, the results of others who computed γ for the core by entirely different methods. They also confirm a higher value of γ in the inner core. The value of γ_a in the lower mantle follows a ρ^?1.35 law, which is reminiscent of the expirical law γρ = constant, commonly used in shock-wave analyses.     

Effects of chemical,mineralogical, and temperature variations on the density and bulk sound velocity of the bottom lower mantle

In this study, we have modeled the density(ρ) and bulk sound velocity(VΦ) profiles of the bottom lower mantle using the experimental thermal equation of state(EoS) parameters of lower-mantle minerals, including bridgmanite, ferropericlase,CaSiO3-perovskite, and post-perovskite. We re-evaluated the literature pressure-volume-temperature relationships of these minerals using a self-consistent pressure scale in order to avoid the long-standing pressure scale problem and to provide more reliable constraints on the thermal EoS parameters. With the obtained thermal EoS parameters, we have constructed the ρ and VΦ profiles of the bottom lower mantle in different composition, mineralogy, and temperature models. Our modelling results show that the variations of chemistry, mineralogy, and temperature have different seismic signatures from each other. The Fe and Al enrichment at the bottom lower mantle can cause an increase in ρ but greatly lower VΦ. A change in mineralogy needs to be considered with the lateral variation in temperature. The cold slabs will be shown as denser regions compared to the normal mantle because of the combined effect of a lower temperature and the presence of a denser post-perovskite at a shallower depth,whereas the hot regions will have a 1–2% lower ρ than the normal mantle. VΦ of both cold slabs and hot regions will be lower than the normal mantle when bridgmanite is the dominant phase in the normal mantle, yet they will be greater once bridgmanite transforms into post-perovskite in the normal mantle. Our modeling also shows that the presence of a(Fe, Al)-enriched bridgmanite thermal pile above the core-mantle boundary will exhibit a seismic signature of enhanced ρ and VΦ, but a reduced VS,which is consistent with the observed seismic anomalies in the large-low-shear-velocity-provinces(LLSVPs). The existence of such a(Fe, Al)-enriched bridgmanite thermal pile thus can help to understand the origin of the LLSVPs. These results provide new insights for the chemical and structure of the deepest lower mantle.     

Thermal equation of state of majorite with MORB composition

In situ synchrotron X-ray diffraction experiments were conducted using the SPEED-1500 multi-anvil press at SPring-8 on majoritic garnet synthesized from natural mid-ocean ridge basalt (MORB), whose chemical composition is close to the average of oceanic crust, at 19 GPa and 2200 K. Pressure-volume-temperature data were collected using a newly developed high-pressure cell assembly to 21 GPa and 1273 K. Data were fit to the high-temperature Birch-Murnaghan equation of state, with fixed values for the ambient cell volume (V₀ = 1574.14(4) Å³) and the pressure derivative of the isothermal bulk modulus (K^′T = 4). This yielded an isothermal bulk modulus of K_T0 = 173(1) GPa, a temperature derivative of the bulk modulus (∂K_T/∂T)_P = −0.022(5) GPa K⁻¹, and a volumetric coefficient of thermal expansivity α = a + bT with values of a = 2.0(3) × 10⁻⁵ K⁻¹ and b = 1.0(5) × 10⁻⁸ K⁻². The derived thermoelastic parameters are very similar to those of pyrope. The density of subducted oceanic crust compared to pyrolitic mantle at the conditions in Earth's transition zone (410-660 km depth) was calculated using these results and previously reported thermoelastic parameters for MORB and pyrolite mineral assembledges. These calculations show that oceanic crust is denser than pyrolitic mantle throughout the mantle transition zone along a normal geotherm, and the density difference is insensitive to temperature at the pressures in lower part of the transition zone.     

The thermal boundary-layer interpretation of D″ and its role as a plume source

The “anomalous” layer in the lowermost mantle, identified as D″ in the notation of K.E. Bullen, appears in the PREM Earth model as a 150 km-thick zone in which the gradient of incompressibility with pressure, $\text{d}\text{K}\text{d}\text{P}$, is almost 1.6, instead of 3.2 as in the overlying mantle. Since PREM shows no accompanying change in density or density gradient, we identify D″ as a thermal boundary layer and not as a chemically distinct zone. The anomaly in $\text{d}\text{K}\text{d}\text{P}$ is related to the temperature gradient by the temperature dependence of K_s, for which we present a thermodynamic identity in terms of accessible quantities. This gives the numerical result (?K_s/?T)_P=?1.6×10⁷ Pa K^?1 for D″ material. The corresponding temperature increment over the D″ range is 840 K. Such a layer cannot be a static feature, but must be maintained by a downward motion of the lower mantle toward the core-mantle boundary with a strong horizontal flow near the base of D″. Assuming a core heat flux of 1.6 × 10¹² W, the downward speed is 0.07 mm y^?1 and the temperature profile in D″, scaled to match PREM data, is approximately exponential with a scale height of 73 km. The inferred thermal conductivity is 1.2 W m^?1 K^?1. Using these values we develop a new analytical model of D″ which is dynamically and thermally consistent. In this model, the lower-mantle material is heated and softened as it moves down into D″ where the strong temperature dependence of viscosity concentrates the horizontal flow in a layer ～ 12 km thick and similarly ensures its removal via narrow plumes.     

Thermal equation of state of magnesite to 32 GPa and 2073 K

Pressure–volume–temperature relations have been measured to 32 GPa and 2073 K for natural magnesite (Mg_0.975Fe_0.015Mn_0.006Ca_0.004CO₃) using synchrotron X-ray diffraction with a multianvil apparatus at the SPring-8 facility. A least-squares fit of the room-temperature compression data to a third-order Birch–Murnaghan equation of state (EOS) yielded K₀ = 97.1 ± 0.5 GPa and K′ = 5.44 ± 0.07, with fixed V₀ = 279.55 ± 0.02 Å³. Further analysis of the high-temperature compression data yielded the temperature derivative of the bulk modulus (∂K_T/∂T)_P = −0.013 ± 0.001 GPa/K and zero-pressure thermal expansion α = a₀ + a₁T with a₀ = 4.03 (7) × 10⁻⁵ K⁻¹ and a₁ = 0.49 (10) × 10⁻⁸ K⁻². The Anderson–Grüneisen parameter is estimated to be δ_T = 3.3. The analysis of axial compressibility and thermal expansivity indicates that the c-axis is over three times more compressible (K_Tc = 47 ± 1 GPa) than the a-axis (K_Tc = 157 ± 1 GPa), whereas the thermal expansion of the c-axis (a₀ = 6.8 (2) × 10⁻⁵ K⁻¹ and a₁ = 2.2 (4) × 10⁻⁸ K⁻²) is greater than that of the a-axis (a₀ = 2.7 (4) × 10⁻⁵ K⁻¹ and a₁ = −0.2 (2) × 10⁻⁸ K⁻²). The present thermal EOS enables us to accurately calculate the density of magnesite to the deep mantle conditions. Decarbonation of a subducting oceanic crust containing 2 wt.% magnesite would result in a 0.6% density reduction at 30 GPa and 1273 K. Using the new EOS parameters we performed thermodynamic calculations for magnesite decarbonation reactions at pressures to 20 GPa. We also estimated stability of magnesite-bearing assemblages in the lower mantle.     

Compositional effect on the pressure derivatives of bulk modulus of silicate melts

Although the bulk moduli (K_T0) of silicate melts have a relatively narrow range of values, the pressure derivatives of the isothermal bulk modulus (K_T0^′) can assume a broad range of values and have an important influence on the compositional dependence of the melt compressibility at high pressure. Based on the melt density data from sink/float experiments at high pressures in the literature, we calculate K_T0^′ using an isothermal equation of state (EOS) (e.g., Birch–Murnaghan EOS and Vinet EOS) with the previously determined values of room-pressure density (ρ₀) and room-pressure bulk modulus (K_T0). The results show that best estimates of K_T0^′ vary considerably from ~ 3 to ~ 7 for different compositions. K_T0^′ is nearly independent of Mg # (molar Mg/(Mg + Fe)), but decreases with SiO₂ content. Hydrous melts have anomalously small K_T0^′ leading to a high degree of compression at high pressures. For anhydrous melts, K_T0^′ is ~ 7 for peridotitic melts, ~ 6 for picritic melts, ~ 5 for komatiitic melts, and ~ 4 for basaltic melts.     

Evaporation from three water bodies of different sizes and climates: Measurements and scaling analysis

Evaporation from small reservoirs, wetlands, and lakes continues to be a theoretical and practical problem in surface hydrology and micrometeorology because atmospheric flows above such systems can rarely be approximated as stationary and planar-homogeneous with no mean subsidence (hereafter referred to as idealized flow state). Here, the turbulence statistics of temperature (T) and water vapor (q) most pertinent to lake evaporation measurements over three water bodies differing in climate, thermal inertia and degree of advective conditions are explored. The three systems included Lac Léman in Switzerland (high thermal inertia, near homogeneous conditions with no appreciable advection due to long upwind fetch), Eshkol reservoir in Israel (intermediate thermal inertia, frequent strong advective conditions) and Tilopozo wetland in Chile (low thermal inertia, frequent but moderate advection). The data analysis focused on how similarity constants for the flux-variance approach, C_T/C_q, and relative transport efficiencies R_wT/R_wq, are perturbed from unity with increased advection or the active role of temperature. When advection is small and thermal inertia is large, C_T/C_q < 1 (or R_wT/R_wq > 1) primarily due to the active role of temperature, which is consistent with a large number of studies conducted over bare soil and vegetated surfaces. However, when advection is significantly large, then C_T/C_q > 1 (orR_wT/R_wq < 1). When advection is moderate and thermal inertia is low, then C_T/C_q ∼ 1. This latter equality, while consistent with Monin–Obukhov similarity theory (MOST), is due to the fact that advection tends to increase C_T/C_q above unity while the active role of temperature tends to decrease C_T/C_q below unity. A simplified scaling analysis derived from the scalar variance budget equation, explained qualitatively how advection could perturb MOST scaling (assumed to represent the idealized flow state).     

The transport of water in subduction zones

The transport of water from subducting crust into the mantle is mainly dictated by the stability of hydrous minerals in subduction zones. The thermal structure of subduction zones is a key to dehydration of the subducting crust at different depths. Oceanic subduction zones show a large variation in the geotherm, but seismicity and arc volcanism are only prominent in cold subduction zones where geothermal gradients are low. In contrast, continental subduction zones have low geothermal gradients, resulting in metamorphism in cold subduction zones and the absence of arc volcanism during subduction. In very cold subduction zone where the geothermal gradient is very low(?5?C/km), lawsonite may carry water into great depths of ?300 km. In the hot subduction zone where the geothermal gradient is high(25?C/km), the subducting crust dehydrates significantly at shallow depths and may partially melt at depths of 80 km to form felsic melts, into which water is highly dissolved. In this case, only a minor amount of water can be transported into great depths. A number of intermediate modes are present between these two end-member dehydration modes, making subduction-zone dehydration various. Low-T/low-P hydrous minerals are not stable in warm subduction zones with increasing subduction depths and thus break down at forearc depths of ?60–80 km to release large amounts of water. In contrast, the low-T/low-P hydrous minerals are replaced by low-T/high-P hydrous minerals in cold subduction zones with increasing subduction depths, allowing the water to be transported to subarc depths of 80–160 km. In either case, dehydration reactions not only trigger seismicity in the subducting crust but also cause hydration of the mantle wedge. Nevertheless, there are still minor amounts of water to be transported by ultrahigh-pressure hydrous minerals and nominally anhydrous minerals into the deeper mantle. The mantle wedge overlying the subducting slab does not partially melt upon water influx for volcanic arc magmatism, but it is hydrated at first with the lowest temperature at the slab-mantle interface, several hundreds of degree lower than the wet solidus of hydrated peridotites. The hydrated peridotites may undergo partial melting upon heating at a later time. Therefore, the water flux from the subducting crust into the overlying mantle wedge does not trigger the volcanic arc magmatism immediately.     

Equation of state of gold and its application to the phase boundaries near 660 km depth in Earth’s mantle

The pressure-volume-temperature equation of state (EOS) of gold is fundamental to high-pressure science because of its widespread use as an internal pressure standard. In particular, the EOS of gold has been used in recent in situ multi-anvil press studies for determination of phase boundaries related to the 660-km seismic discontinuity. These studies show that the boundaries are lower by 2 GPa than expected from the depth of the 660-km discontinuity. Here we report a new P-V-T EOS of gold based on the inversion of quasi-hydrostatic compression and shock wave data using the Mie-Grüneisen relation and the Birch-Murnaghan-Debye equation. The previously poorly constrained pressure derivative of isothermal bulk modulus and the volume dependence of Grüneisen parameter (q=d lnγ/d ln V) are determined by including both phonon and electron effects implicitly: K′_0T=5.0±0.2 and q=1.0±0.1. This combined with other accurately measured parameters enables us to calculate pressure at a given volume and temperature. At 660-km depth conditions, this new EOS yields 1.0±0.2 GPa higher pressure than Anderson et al.’s EOS which has been used in the multi-anvil experiments. However, after the correction, there still exists a 1.5-GPa discrepancy between the post-spinel boundary measured by multi-anvil studies and the 660-km discontinuity. Other potential error sources, such as thermocouple emf dependence on pressure or systematic errors in spectroradiometry, should be investigated. Theoretical and experimental studies to better understand electronic and anharmonic effects in gold at high P-T are also needed.     

Melting temperature distribution and fractionation in the lower mantle

The melting curve of perovskite MgSiO₃ and the liquidus and solidus curves of the lower mantle were estimated from thermodynamic data and the results of experiments on phase changes and melting in silicates.The initial slope of the melting curve of perovskite MgSiO₃ was obtained as $\text{d}\text{T}{}_{\mathrm{<mtext>m</mtext>}}\text{/}\text{d}\text{P?77}\text{K}\text{GPa}{}^{\mathrm{?1}}$ at 23 GPa. The melting curve of perovskite was expressed by the Kraut-Kennedy equation as $\text{T}{}_{\mathrm{<mtext>m</mtext>}}\text{(}\text{K}\text{)=917(}\text{1+29.6ΔV}\text{V}{}_0\text{)}$, where T_m?2900 K and P?23 GPa; and by the Simon equation, $\text{P(}\text{GPa}\text{)?23=21.2[}\text{(T}{}_{\mathrm{<mtext>m</mtext>}}\text{(}\text{K}\text{)}\text{2900)}{}^{1.75}\text{?1}\text{]}$.The liquidus curve of the lower mantle was estimated as $\text{T}{}_{\mathrm{<mtext>liq</mtext>}}\text{? 0.9 T}{}_{\mathrm{<mtext>m</mtext>}}$ (perovskite) and this gives the liquidus temperature T_liq=7000 ±500 K at the mantle-core boundary. The solidus curve of the lower mantle was also estimated by extrapolating the solidus curve of dry peridotite using the slope of the solidus curve of magnesiowüstite at high pressures. The solidus temperature is ～ 5000 K at the base of the lower mantle. If the temperature distribution of the mantle was 1.5 times higher than that given by the present geotherm in the early stage of the Earth's history, partial melting would have proceeded into the deep interior of the lower mantle.Estimation of the density of melts in the MgOFeOSiO₂ system for lower mantle conditions indicates that the initial melt formed by partial fusion of the lower mantle would be denser than the residual solid because of high concentration of iron into the melt. Thus, the melt generated in the lower mantle would tend to move downward toward the mantle-core boundary. This downward transportation of the melt in the lower mantle might have affected the chemistry of the lower mantle, such as in the D″ layer, and the distribution of the radioactive elements between mantle and core.     

Pressure dependence of the thermal Grüneisen parameter,with application to the earth's lower mantle and outer core

By considering high-temperature (classical) thermal oscillations of atoms in certain simple crystal structures with purely central interatomic forces, the treatment of anharmonic oscillations is generalised to random three-dimensional motion, yielding the Vashchenko and Zubarev relationship for the Grüneisen ratio γ at any pressure. If one-dimensional atomic oscillations only are considered the equation reduces to the Dugdale-MacDonald expression. To account for non-central forces additional terms must be introduced, giving:

$\text{γ=}\text{1}\text{2}\text{d}\text{K}\text{d}\text{P}\text{?}\text{5}\text{6}\text{+}\text{2}\text{9}\text{P}\text{K}\text{?}\text{f}\text{18K}\text{+}\text{1}\text{6}\text{d}\text{f}\text{d}\text{P}\text{1?}\text{4}\text{3}\text{P}\text{K}\text{+}\text{f}\text{3K}$

where f = 0 for purely central forces. Calculations of f in terms of the Poisson ratio for different crystal structures have not been made, but for many materials the central-force approximation suffices. This is believed to be true both for the outer core $\text{(}\text{γ}\text{≈1.4)}$ and for the close-packed structures of the lower mantle $\text{(}\text{γ}\text{≈1.0)}$. For the upper mantle non-central atomic forces are important and we have no estimate of ($\text{γ}$ independently of laboratory values for plausible minerals which suggest γ ≈ 0.8.     

Approximate solutions for the formation of the lithosphere

The lithosphere is interpreted as a thermal boundary layer. Approximate solutions of the boundary layer cooling problem are developed which include mantle radioactivity, partial melt in the asthenosphere, a temperature gradient in the asthenosphere, and a non-zero lithospheric thickness at the ridge crests. The cooling history of oceanic lithosphere is found to be remarkably insensitive to assumptions about the amount of radioactivity in the upper mantle and the extent of melting in the asthenosphere. Determinations of the thickness of oceanic lithosphere and the depths of oceans as a function of age are in excellent agreement with boundary layer predictions which include a heat flux from the asthenosphere. However, the determinations do not resolve how much of the total asthenospheric heat flux might be caused by a temperature gradient in the asthenosphere. Simple thermal arguments indicate that the initial lithospheric thickness, L₀, at ridge crests should depend on the local half-spreading rate, V, as L₀ = 3 km/V(cm/year).     

The effect of ferromagnetism on the equation of state of Fe3C studied by first-principles calculations

First-principles calculations have been used to determine the equation of state of Fe₃C in both its low-pressure magnetically ordered and high-pressure non-magnetically ordered states; at 0 K the ferromagnetic transition was found to occur at about 60 GPa. In the high pressure, non-magnetically ordered regime at 0 K the material may be described by a Birch-Murnaghan third-order equation of state with V₀=8.968(7) ų per atom, K₀=316.62(2) GPa and K′=4.30(2). At atmospheric pressure the ferromagnetic phase transition in Fe₃C occurs at ∼483 K; preliminary measurements of the thermal expansion by powder neutron diffraction show that this transition produces a large effect on thermoelastic properties. The volumetric thermal expansion coefficient in the paramagnetic phase was found to be 4.34×10⁻⁵ K⁻¹ at T∼550 K. By applying a thermal expansion correction to the calculated equation of state at 0 K, predicted values for the density and adiabatic incompressibility of this material at core pressures and temperatures were obtained. These results appear to be sufficiently different from seismological data so as to preclude Fe₃C as the major inner core-forming phase.     

On the origin of El Chichón volcano and subduction of Tehuantepec Ridge: A geodynamical perspective

The origin of El Chichón volcano is poorly understood, and we attempt in this study to demonstrate that the Tehuantepec Ridge (TR), a major tectonic discontinuity on the Cocos plate, plays a key role in determining the location of the volcano by enhancing the slab dehydration budget beneath it. Using marine magnetic anomalies we show that the upper mantle beneath TR undergoes strong serpentinization, carrying significant amounts of water into subduction. Another key aspect of the magnetic anomaly over southern Mexico is a long-wavelength (∼ 150 km) high amplitude (∼ 500 nT) magnetic anomaly located between the trench and the coast. Using a 2D joint magnetic-gravity forward model, constrained by the subduction P–T structure, slab geometry and seismicity, we find a highly magnetic and low-density source located at 40–80 km depth that we interpret as a partially serpentinized mantle wedge formed by fluids expelled from the subducting Cocos plate. Using phase diagrams for sediments, basalt and peridotite, and the thermal structure of the subduction zone beneath El Chichón we find that ∼ 40% of sediments and basalt dehydrate at depths corresponding with the location of the serpentinized mantle wedge, whereas the serpentinized root beneath TR strongly dehydrates (∼90%) at depths of 180-200 km comparable with the slab depths beneath El Chichón (200-220 km). We conclude that this strong deserpentinization pulse of mantle lithosphere beneath TR at great depths is responsible for the unusual location, singularity and, probably, the geochemically distinct signature (adakitic-like) of El Chichón volcano.     

Thermodynamic properties of calcium ferrite-type MgAl2O4: A first principles study

The diamond anvil cell experiments have revealed that the calcium ferrite(CF)-type aluminous phase is probably an important component of subducted mid-oceanic ridge basalt(MORB) in the lower mantle. In this study, we have performed first principles lattice dynamics calculations for the Mg Al_2O_4 end-member of the aluminous phase based on density functional perturbation theory using two functionals within the local density approximation(LDA) and generalized gradient approximation(GGA) for bracketing the calculated properties at their lower and upper limits, respectively. A simple empirical pressure correction at zero temperature has been applied to both LDA and GGA. The results of room-temperature equation of state(EOS) and zero-pressure thermal expansion calculated by GGA with pressure correction have shown the best agreement with available experimental data. The high-pressure and temperature thermodynamic properties have been obtained using the GGA with correction method. The pressure-volume relations are fitted with a third-order high-temperature Birch-Murnaghan EOS. The isobaric heat capacity, the coefficient of thermal expansion and isothermal bulk modulus are fitted with polynomials and their coefficients are reported in the range of 0–40 GPa and 300–2000 K. The density profile of MORB estimated using the computational thermo-elastic constants supports the hypothesis that the subducted oceanic slabs could gain enough downwelling forces into the lower mantle.