Welcome to the new version of CaltechAUTHORS. Login is currently restricted to library staff. If you notice any issues, please email coda@library.caltech.edu
Published November 2015 | Published
Journal Article Open

A method of spherical harmonic analysis in the geosciences via hierarchical Bayesian inference

Abstract

The problem of decomposing irregular data on the sphere into a set of spherical harmonics is common in many fields of geosciences where it is necessary to build a quantitative understanding of a globally varying field. For example, in global seismology, a compressional or shear wave speed that emerges from tomographic images is used to interpret current state and composition of the mantle, and in geomagnetism, secular variation of magnetic field intensity measured at the surface is studied to better understand the changes in the Earth's core. Optimization methods are widely used for spherical harmonic analysis of irregular data, but they typically do not treat the dependence of the uncertainty estimates on the imposed regularization. This can cause significant difficulties in interpretation, especially when the best-fit model requires more variables as a result of underestimating data noise. Here, with the above limitations in mind, the problem of spherical harmonic expansion of irregular data is treated within the hierarchical Bayesian framework. The hierarchical approach significantly simplifies the problem by removing the need for regularization terms and user-supplied noise estimates. The use of the corrected Akaike Information Criterion for picking the optimal maximum degree of spherical harmonic expansion and the resulting spherical harmonic analyses are first illustrated on a noisy synthetic data set. Subsequently, the method is applied to two global data sets sensitive to the Earth's inner core and lowermost mantle, consisting of PKPab-df and PcP-P differential traveltime residuals relative to a spherically symmetric Earth model. The posterior probability distributions for each spherical harmonic coefficient are calculated via Markov Chain Monte Carlo sampling; the uncertainty obtained for the coefficients thus reflects the noise present in the real data and the imperfections in the spherical harmonic expansion.

Additional Information

© 2015 The Authors. Published by Oxford University Press on behalf of The Royal Astronomical Society. Accepted 2015 August 27. Received 2015 August 26; in original form 2015 August 25. Jack Muir acknowledges the financial support of the Research School of Physics at the Australian National University, and the General Sir John Monash Foundation during this project.

Attached Files

Published - Geophys._J._Int.-2015-Muir-1164-71.pdf

Files

Geophys._J._Int.-2015-Muir-1164-71.pdf
Files (3.1 MB)
Name Size Download all
md5:1ebff1fd6406a114bcdb292821ffaa86
3.1 MB Preview Download

Additional details

Created:
August 22, 2023
Modified:
October 17, 2023