Robust estimates of the ratio between S- and P-wave velocity anomalies in the Earth's mantle using normal modes

Conference item

Seismic tomography allows us to image the interior of the Earth. In general, to determine the nature of seismic anomalies, constraints on more than one seismic parameter are required. For example, the ratio R between perturbations in v_s and v_p (dlnv_s and dlnv_p, respectively) is studied extensively in the lowermost mantle and interpreted in terms of thermal and/or chemical anomalies. However, to jointly interpret tomographic models of variations in v_s and v_p or their ratio R, it is essential for them to share the same local resolution. Most existing models do not provide resolution information, and thus cannot guarantee to honour this condition. In addition, uncertainties are typically not provided, making it difficult to robustly interpret the ratio R = dlnv_s / dlnv_p. To overcome these issues, we utilise the recently developed SOLA tomographic method, a variant of the linear Backus–Gilbert inversion scheme. SOLA retrieves local-average model estimates, together with information on their uncertainties, whilst it also provides direct control on model resolution through target kernels. In this contribution, we apply SOLA to normal-mode data with sensitivity to both v_s and v_p, as well as density throughout the mantle. Specifically, we aim to develop models of both v_s and v_p with the same local resolution. We test our methodology and approach using synthetic tests for various noise cases (random noise, data noise or also additional 3D Earth noise due to variations in other physical parameters than the one of interest). We find that the addition of the 3D noise increases the uncertainties in our model estimates significantly, only allowing us to find model estimates in six or four layers for v_s and v_p, respectively. While the synthetic tests indicate that no satisfactory density models can be obtained, we easily manage to construct models of dlnv_s and dlnv_p with almost identical resolution, from which the ratio R can be robustly inferred. The obtained values of R in our synthetic experiments significantly depend on the noise case considered and the method used to calculate it, with the addition of 3D noise always leading to an overestimate of R. When applying our approach to real data, we obtain values of R in the range of 2.5–4.0 in the lowest 600 km of the mantle, which are consistent with previous studies. Our model estimates with related resolving kernels and uncertainties can be used to test geodynamic model predictions to provide further insights into the temperature and composition of the mantle.

Authors: Restelli, F; Zaroli, C; Koelemeijer, P

• Publication status: Published
• Peer review status: Peer reviewed
• Publisher: Elsevier
• Host title: Proceedings of the 17th Symposium of Study of the Earth’s Deep Interior (SEDI)
• Journal: Physics of the Earth and Planetary Interiors
• Volume: 347
• Article number: 107135
• Publication date: 2023-12-28
• Acceptance date: 2023-12-22
• Event title: 17th Study of the Earth's Deep Interior (SEDI 2022)
• Event location: Zurich, Switzerland
• Event website: https://sedi2022.earth.sinica.edu.tw/
• Event start date: 2022-07-11
• Event end date: 2024-07-15
• DOI: 10.1016/j.pepi.2023.107135
• EISSN: 1872-7395
• ISSN: 0031-9201
• Language: English
• Keywords: composition and structure of the mantle; seismic tomography; free oscillations; uncertainties