Abstract

Seismic free oscillations, or normal modes, provide a convenient tool to calculate low-frequency seismograms in heterogeneous Earth models. A procedure called ‘full mode coupling’ allows the seismic response of the Earth to be computed. However, in order to be theoretically exact, such calculations must involve an infinite set of modes. In practice, only a finite subset of modes can be used, introducing an error into the seismograms. By systematically increasing the number of modes beyond the highest frequency of interest in the seismograms, we investigate the convergence of full-coupling calculations. As a rule-of-thumb, it is necessary to couple modes 1–2 mHz above the highest frequency of interest, although results depend upon the details of the Earth model. This is significantly higher than has previously been assumed. Observations of free oscillations also provide important constraints on the heterogeneous structure of the Earth. Historically, this inference problem has been addressed by the measurement and interpretation of splitting functions. These can be seen as secondary data extracted from low frequency seismograms. The measurement step necessitates the calculation of synthetic seismograms, but current implementations rely on approximations referred to as self- or group-coupling and do not use fully accurate seismograms. We therefore also investigate whether a systematic error might be present in currently published splitting functions. We find no evidence for any systematic bias, but published uncertainties must be doubled to properly account for the errors due to theoretical omissions and regularization in the measurement process. Correspondingly, uncertainties in results derived from splitting functions must also be increased. As is well known, density has only a weak signal in low-frequency seismograms. Our results suggest this signal is of similar scale to the true uncertainties associated with currently published splitting functions. Thus, it seems that great care must be taken in any attempt to robustly infer details of Earth's density structure using current splitting functions.

Highlights

  • Our understanding of the Earth’s large-scale interior structure and dynamics draws heavily on observations of seismic free oscillations (‘normal modes’)

  • It is not clear how to define a meaningful prior for the splitting function inversion, and in any case such analysis would neglect the inherent non-linearity of the measurement process

  • By systematically increasing the number of modes in full-coupling calculations, we have demonstrated that accurate full-coupling requires inclusion of modes significantly beyond the highest frequency range of interest in the seismograms

Read more

Summary

INTRODUCTION

Our understanding of the Earth’s large-scale interior structure and dynamics draws heavily on observations of seismic free oscillations (‘normal modes’). Al-Attar et al (2012) highlighted that use of either the self- or group-coupling approximations when computing synthetic spectra leads to errors of similar magnitude to the signal likely to be attributable to lateral variations in density within the Earth. Building on the existing body of work, the first aim of the present paper is to identify more precisely the conditions under which full-coupling can be regarded as ‘sufficiently accurate’ This information permits us to perform high-quality synthetic experiments to address the second question: can robust information about earth structure be inferred from measurements of ‘splitting functions’ (Woodhouse & Giardini 1985; Giardini et al 1987, 1988)? Repeating the experiment for models where density perturbations are set to zero, we conclude that splitting functions are contaminated by a significant approximation error, a similar regularization error and contain a density signal of the level of the total measurement uncertainty

THEORETICALB AC KG RO U N D
Reference eigenfunctions
Computing synthetic seismograms
Full-coupling
The self- and group-coupling approximations
Normal modes of an aspherical earth model
From spectra to structure: the splitting functions
Splitting functions for isolated groups of multiplets
Splitting functions in the real Earth
T RU N C AT IONER RO R S I N FULL-COUPLING SEISMOGRAMS
The scale of acceptable error
Observations using S20RTS
Model-dependence of the truncation error
Summary
THE INTERPRETABILITY OF SPLITTING FUNCTION MEASUREMENTS
Measurement of synthetic splitting functions
Observations from synthetic splitting functions
CONCLUSIONS
The seismic wave equation
Spherical-earth eigenfunctions
Completeness and orthogonality
Seismograms in a spherically symmetric model
Seismograms in an aspherical earth model
Self- and group-coupling
Findings
Normal modes in aspherical models
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call