On the bulk composition of the lower mantle: Predictions and limitations from generalized inversion of radial seismic profiles

J. Matas, J. Bass, Y. Ricard, E. Mattern, M. S.T. Bukowinski

Research output: Contribution to journalArticle

Abstract

We examine the problem of obtaining the thermal structure and bulk chemical composition of the lower mantle from its seismologically determined velocity and density profiles, and the most recent results on the elastic properties of the relevant phases (including, of particular importance, shear moduli). A novel aspect of this paper is the application of an iterative technique solving generalized non-linear inverse problem, which allows us to simultaneously consider a complex chemical system (the MgO-FeO-SiO2-Al2O3 -CaO system, which includes all major components in the lower mantle), and to rigorously evaluate the full covariance and resolution matrices. The effects of experimental uncertainties in the shear moduli are carefully accounted for. We show that although the a posteriori uncertainties in the results for lower-mantle compositions are relatively large, the averaged lower-mantle Mg/Si ratio should be lower than 1.3 in order to satisfactorily fit the 1-D seismic profiles. Two distinct families of best-fitting models are determined. The first is based upon a value for the pressure derivative of the perovskite shear modulus that is representative of various existing experimental measurements (μ′0 = 1.8). Under this assumption, it is not possible to match the lower mantle seismic properties with an adiabatic geotherm and uniform chemical composition. Instead, this family of solutions is characterized by a geotherm with large temperature gradients (dT/dz increases from 0.5 to 0.9 K km-1 between 800 and 2700 km and the temperature reaches 3400 K at the depth of 2700 km), and a depth dependent bulk composition with an Mg/Si ratio decreasing from 1.18 ± 0.14 to 1.03 ± 0.16 between 800 and 2700 km. The second family of solutions is obtained when we attempt to fit the lower mantle with a simpler compositional and thermal structure. This can only be done when the pressure derivative of the shear modulus for perovskite is close to the most recent values obtained by Brillouin spectroscopy, that is, with a μ′0 close to 1.6 instead of 1.8. The resulting temperature gradient is 0.25 K km-1 in the upper part of the lower mantle and 0.5 K km-1 below 1700 km depth; the geotherm reaches 2800 K at a depth of 2700 km. Corresponding Mg/Si ratio remains rather constant and close to 1.16 throughout the lower mantle. We show that the temperature gradient is strongly correlated with the pressure derivative μ′0 of the shear modulus of perovskite: lower values of μ′0 imply lower thermal gradients. We also discuss the importance of the Bullen parameter as an additional constraint. In order to refine conclusions on the lower-mantle structure, additional independent observables, such as accurate observations on electrical conductivity and 1-D Q profiles, are necessary.

Original languageEnglish (US)
Pages (from-to)764-780
Number of pages17
JournalGeophysical Journal International
Volume170
Issue number2
DOIs
StatePublished - Aug 1 2007

Fingerprint

lower mantle
Earth mantle
Elastic moduli
Thermal gradients
inversions
shear modulus
profiles
prediction
predictions
Chemical analysis
Derivatives
shear
perovskite
temperature gradient
temperature gradients
thermal structure
Inverse problems
chemical composition
inversion
Spectroscopy

Keywords

  • Density
  • Inverse problem
  • Mantle
  • Mineralogy
  • Seismic velocities

ASJC Scopus subject areas

  • Geophysics
  • Geochemistry and Petrology

Cite this

On the bulk composition of the lower mantle : Predictions and limitations from generalized inversion of radial seismic profiles. / Matas, J.; Bass, J.; Ricard, Y.; Mattern, E.; Bukowinski, M. S.T.

In: Geophysical Journal International, Vol. 170, No. 2, 01.08.2007, p. 764-780.

Research output: Contribution to journalArticle

@article{12e57328eb374382a51a559fc352bb2a,
title = "On the bulk composition of the lower mantle: Predictions and limitations from generalized inversion of radial seismic profiles",
abstract = "We examine the problem of obtaining the thermal structure and bulk chemical composition of the lower mantle from its seismologically determined velocity and density profiles, and the most recent results on the elastic properties of the relevant phases (including, of particular importance, shear moduli). A novel aspect of this paper is the application of an iterative technique solving generalized non-linear inverse problem, which allows us to simultaneously consider a complex chemical system (the MgO-FeO-SiO2-Al2O3 -CaO system, which includes all major components in the lower mantle), and to rigorously evaluate the full covariance and resolution matrices. The effects of experimental uncertainties in the shear moduli are carefully accounted for. We show that although the a posteriori uncertainties in the results for lower-mantle compositions are relatively large, the averaged lower-mantle Mg/Si ratio should be lower than 1.3 in order to satisfactorily fit the 1-D seismic profiles. Two distinct families of best-fitting models are determined. The first is based upon a value for the pressure derivative of the perovskite shear modulus that is representative of various existing experimental measurements (μ′0 = 1.8). Under this assumption, it is not possible to match the lower mantle seismic properties with an adiabatic geotherm and uniform chemical composition. Instead, this family of solutions is characterized by a geotherm with large temperature gradients (dT/dz increases from 0.5 to 0.9 K km-1 between 800 and 2700 km and the temperature reaches 3400 K at the depth of 2700 km), and a depth dependent bulk composition with an Mg/Si ratio decreasing from 1.18 ± 0.14 to 1.03 ± 0.16 between 800 and 2700 km. The second family of solutions is obtained when we attempt to fit the lower mantle with a simpler compositional and thermal structure. This can only be done when the pressure derivative of the shear modulus for perovskite is close to the most recent values obtained by Brillouin spectroscopy, that is, with a μ′0 close to 1.6 instead of 1.8. The resulting temperature gradient is 0.25 K km-1 in the upper part of the lower mantle and 0.5 K km-1 below 1700 km depth; the geotherm reaches 2800 K at a depth of 2700 km. Corresponding Mg/Si ratio remains rather constant and close to 1.16 throughout the lower mantle. We show that the temperature gradient is strongly correlated with the pressure derivative μ′0 of the shear modulus of perovskite: lower values of μ′0 imply lower thermal gradients. We also discuss the importance of the Bullen parameter as an additional constraint. In order to refine conclusions on the lower-mantle structure, additional independent observables, such as accurate observations on electrical conductivity and 1-D Q profiles, are necessary.",
keywords = "Density, Inverse problem, Mantle, Mineralogy, Seismic velocities",
author = "J. Matas and J. Bass and Y. Ricard and E. Mattern and Bukowinski, {M. S.T.}",
year = "2007",
month = "8",
day = "1",
doi = "10.1111/j.1365-246X.2007.03454.x",
language = "English (US)",
volume = "170",
pages = "764--780",
journal = "Geophysical Journal International",
issn = "0956-540X",
publisher = "Wiley-Blackwell",
number = "2",

}

TY - JOUR

T1 - On the bulk composition of the lower mantle

T2 - Predictions and limitations from generalized inversion of radial seismic profiles

AU - Matas, J.

AU - Bass, J.

AU - Ricard, Y.

AU - Mattern, E.

AU - Bukowinski, M. S.T.

PY - 2007/8/1

Y1 - 2007/8/1

N2 - We examine the problem of obtaining the thermal structure and bulk chemical composition of the lower mantle from its seismologically determined velocity and density profiles, and the most recent results on the elastic properties of the relevant phases (including, of particular importance, shear moduli). A novel aspect of this paper is the application of an iterative technique solving generalized non-linear inverse problem, which allows us to simultaneously consider a complex chemical system (the MgO-FeO-SiO2-Al2O3 -CaO system, which includes all major components in the lower mantle), and to rigorously evaluate the full covariance and resolution matrices. The effects of experimental uncertainties in the shear moduli are carefully accounted for. We show that although the a posteriori uncertainties in the results for lower-mantle compositions are relatively large, the averaged lower-mantle Mg/Si ratio should be lower than 1.3 in order to satisfactorily fit the 1-D seismic profiles. Two distinct families of best-fitting models are determined. The first is based upon a value for the pressure derivative of the perovskite shear modulus that is representative of various existing experimental measurements (μ′0 = 1.8). Under this assumption, it is not possible to match the lower mantle seismic properties with an adiabatic geotherm and uniform chemical composition. Instead, this family of solutions is characterized by a geotherm with large temperature gradients (dT/dz increases from 0.5 to 0.9 K km-1 between 800 and 2700 km and the temperature reaches 3400 K at the depth of 2700 km), and a depth dependent bulk composition with an Mg/Si ratio decreasing from 1.18 ± 0.14 to 1.03 ± 0.16 between 800 and 2700 km. The second family of solutions is obtained when we attempt to fit the lower mantle with a simpler compositional and thermal structure. This can only be done when the pressure derivative of the shear modulus for perovskite is close to the most recent values obtained by Brillouin spectroscopy, that is, with a μ′0 close to 1.6 instead of 1.8. The resulting temperature gradient is 0.25 K km-1 in the upper part of the lower mantle and 0.5 K km-1 below 1700 km depth; the geotherm reaches 2800 K at a depth of 2700 km. Corresponding Mg/Si ratio remains rather constant and close to 1.16 throughout the lower mantle. We show that the temperature gradient is strongly correlated with the pressure derivative μ′0 of the shear modulus of perovskite: lower values of μ′0 imply lower thermal gradients. We also discuss the importance of the Bullen parameter as an additional constraint. In order to refine conclusions on the lower-mantle structure, additional independent observables, such as accurate observations on electrical conductivity and 1-D Q profiles, are necessary.

AB - We examine the problem of obtaining the thermal structure and bulk chemical composition of the lower mantle from its seismologically determined velocity and density profiles, and the most recent results on the elastic properties of the relevant phases (including, of particular importance, shear moduli). A novel aspect of this paper is the application of an iterative technique solving generalized non-linear inverse problem, which allows us to simultaneously consider a complex chemical system (the MgO-FeO-SiO2-Al2O3 -CaO system, which includes all major components in the lower mantle), and to rigorously evaluate the full covariance and resolution matrices. The effects of experimental uncertainties in the shear moduli are carefully accounted for. We show that although the a posteriori uncertainties in the results for lower-mantle compositions are relatively large, the averaged lower-mantle Mg/Si ratio should be lower than 1.3 in order to satisfactorily fit the 1-D seismic profiles. Two distinct families of best-fitting models are determined. The first is based upon a value for the pressure derivative of the perovskite shear modulus that is representative of various existing experimental measurements (μ′0 = 1.8). Under this assumption, it is not possible to match the lower mantle seismic properties with an adiabatic geotherm and uniform chemical composition. Instead, this family of solutions is characterized by a geotherm with large temperature gradients (dT/dz increases from 0.5 to 0.9 K km-1 between 800 and 2700 km and the temperature reaches 3400 K at the depth of 2700 km), and a depth dependent bulk composition with an Mg/Si ratio decreasing from 1.18 ± 0.14 to 1.03 ± 0.16 between 800 and 2700 km. The second family of solutions is obtained when we attempt to fit the lower mantle with a simpler compositional and thermal structure. This can only be done when the pressure derivative of the shear modulus for perovskite is close to the most recent values obtained by Brillouin spectroscopy, that is, with a μ′0 close to 1.6 instead of 1.8. The resulting temperature gradient is 0.25 K km-1 in the upper part of the lower mantle and 0.5 K km-1 below 1700 km depth; the geotherm reaches 2800 K at a depth of 2700 km. Corresponding Mg/Si ratio remains rather constant and close to 1.16 throughout the lower mantle. We show that the temperature gradient is strongly correlated with the pressure derivative μ′0 of the shear modulus of perovskite: lower values of μ′0 imply lower thermal gradients. We also discuss the importance of the Bullen parameter as an additional constraint. In order to refine conclusions on the lower-mantle structure, additional independent observables, such as accurate observations on electrical conductivity and 1-D Q profiles, are necessary.

KW - Density

KW - Inverse problem

KW - Mantle

KW - Mineralogy

KW - Seismic velocities

UR - http://www.scopus.com/inward/record.url?scp=34547171158&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=34547171158&partnerID=8YFLogxK

U2 - 10.1111/j.1365-246X.2007.03454.x

DO - 10.1111/j.1365-246X.2007.03454.x

M3 - Article

AN - SCOPUS:34547171158

VL - 170

SP - 764

EP - 780

JO - Geophysical Journal International

JF - Geophysical Journal International

SN - 0956-540X

IS - 2

ER -