Compositional constraints on the equation of state and thermal properties of the lower mantle

Citation
Fd. Stacey et Dg. Isaak, Compositional constraints on the equation of state and thermal properties of the lower mantle, GEOPHYS J I, 146(1), 2001, pp. 143-154
Citations number
42
Categorie Soggetti
Earth Sciences
Journal title
GEOPHYSICAL JOURNAL INTERNATIONAL
ISSN journal
0956540X → ACNP
Volume
146
Issue
1
Year of publication
2001
Pages
143 - 154
Database
ISI
SICI code
0956-540X(200107)146:1<143:CCOTEO>2.0.ZU;2-X
Abstract
By extrapolating the lower mantle equation of state (EoS) to P = 0, T = 290 K, we determine the EoS parameters that are compatible with a mixture of ( Mg, Fe)SiO3 perovskite (with a small admixture of Al2O3), (Mg, Fe)O magnesi owustite and CaSiO3 perovskite in arbitrary proportions and with arbitrary Fe/(Fe + Mg) ratio. The parameters fitted are density, rho, adiabatic incom pressibility, K-S, and its pressure derivative, K'(S) =(partial derivativeK (S)/partial derivativeP)(S). The first stage is adiabatic extrapolation to P = 0, T = T-0, that is, to the foot of the lower mantle adiabat, at which k'(0)(T-0) is allowed to have any value between 3.8 and 4.6, and 1500 K les s than or equal to T-0 less than or equal to 2000 K. It is important to use an equation for which the lower mantle fitting does not prescribe K'(0)(T- 0) and this rules out the third-order Birch theory, which gives a seriously wrong value. The further extrapolation to 290 K at P = 0 uses thermodynami c relationships with maximum generality, allowing all of the following ther moelastic parameters to be arbitrary functions of temperature: K; rho; Grun eisen parameter, gamma; q = (partial derivative 1n gamma/partial derivative 1n V)(T), where V is volume; volume coefficient of thermal expansion, alph a; adiabatic Anderson-Gruneisen parameter, delta (S) = (1/alpha) (partial d erivative 1n K-S/partial derivativeT)(P); and the mixed P, T derivative (pa rtial derivativeK'(S)/partial derivativeT)(P). The heat capacity at constan t volume, C-V, is assumed to follow the Debye function, so alpha is control led by that. The temperature dependences of the dimensionless parameters ga mma, q and delta (S) at P = 0 are slight. We find gamma to be precisely ind ependent of T at constant V. The parameter dK'(0)/dT increases strongly wit h T, as well as with the assumed value of K'(0)(T-0), where K'(0) is K'(S) at P = 0. The fitting disallows significant parameter ranges. In particular , we find solutions only if K'(0)(T0)greater than or equal to4.2 and the 29 0 K value of K'(0) for Mg perovskite is less than 3.8. Conclusions about co mposition are less secure, partly because of doubt about individual mineral properties. The volume of magnesiowustite is found to be between 10 and 25 per cent for respective T-0 values of 2000 and 1500 K, but the Ca-perovski te volume is no more than 6 per cent and has little influence on the other conclusions. The resulting overall Fe/(Fe + Mg) ratio is 0.12 to 0.15. Alth ough this ratio is higher than expected for a pyrolite composition, the rat io depends critically on the assumed mineral densities; some adjustment of the mineral mix may need to be considered.