The joint work A Bayesian Formulation for Estimating the Composition of Earth’s Crust with the Dartmouth Earth science department (G Pease, B. Keller) along with A. Gelb (Dartmouth Math) has been accepted and published in JGR Solid Earth.
Abstract: Due to the inaccessibility of Earth’s deep interior, geologists have long attempted to estimate the composition of the continental crust from its seismic properties. Despite numerous sources of error including nonuniqueness in the mapping between composition and seismic properties, the corresponding uncertainties have typically been estimated qualitatively at best. We propose a Bayesian approach that uses mineralogical modeling to combine prior knowledge about the composition of the crust with seismic data to give a posterior distribution of the predicted composition at any location, combined with a Monte Carlo simulation to estimate the average composition of the Earth’s crust. Our approach yields an estimated composition of 59.5% silica in the upper crust (90% credible interval 58.9 %–60.1%), 57.9% in the middle crust (90% credible interval 57.2%–58.6%), and 53.6% in the lower crust (90% credible interval 53.0%–54.2%). Our estimate exhibits less compositional stratification over depth and a more intermediate composition in the upper and middle crust than previous estimates. Testing our approach on a simulated crust reveals the importance of prior assumptions in estimating the composition of the crust from its seismic properties, and suggests that future work should focus on quantifying those assumptions.