Array-based iterative measurements of SmKS travel times and their constraints on outermost core structure

Vigorous convection in Earth's outer core led to the suggestion that it is chemically homogeneous. However, there is increasing seismic evidence for structural complexities close to the outer core's upper and lower boundaries. Both body waves and normal mode data have been used to estimate a P wave velocity, V_p, at the top of the outer core (the E’ layer), which is lower than that in the Preliminary Reference Earth Model. However, these low V_p models do not agree on the form of this velocity anomaly. One reason for this is the difficulty in retrieving and measuring SmKS arrival times. To address this issue, we propose a novel approach using data from seismic arrays to iteratively measure SmKS-SKKS-differential travel times. This approach extracts individual SmKS signal from mixed waveforms of the SmKS series, allowing us to reliably measure differential travel times. We successfully use this method to measure SmKS time delays from earthquakes in the Fiji‐Tonga and Vanuatu subduction zones. SmKS time delays are measured by waveform cross correlation between SmKS and SKKS, and the cross‐correlation coefficient allows us to access measurement quality. We also apply this iterative scheme to synthetic SmKS seismograms to investigate the 3‐D mantle structure's effects. The mantle structure corrections are not negligible for our data, and neglecting them could bias the V_p estimation of uppermost outer core. After mantle structure corrections, we can still see substantial time delays of S3KS, S4KS, and S5KS, supporting a low V_p at the top of Earth's outer core

Mantle Transition Zone Receiver Functions for Bermuda:Automation, Quality Control, and Interpretation

The origin of the Bermuda rise remains ambiguous, despite, or perhaps because of, the existence of sometimes incongruous seismic wave‐speed and discontinuity models in the sub‐Bermudian mantle. Hence, whether Bermuda is the surface manifestation of a mantle plume remains in question. Using the largest data set of seismic records from Bermuda to date, we estimate radial receiver functions at the Global Seismographic Network station BBSR in multiple frequency bands, using iterative time‐domain deconvolution. Motivated by synthetic experiments using axisymmetric spectral‐element forward waveform modeling, we devise a quality metric for our receiver functions to aid in the automation and reproduction of mantle transition zone discontinuity studies. We interpret the complex signals we observe by considering the mineralogical controls on mantle transition zone discontinuity structure, and conclude that our results are likely to be indicative of a thicker than average mantle transition zone. Our result is incompatible with the canonical model of a whole mantle plume in an olivine dominated mantle; however, considerations of phase transitions in the garnet system would allow us to reconcile our observations with the possible presence of a through‐going hot thermal anomaly beneath Bermuda

Evidence for a kilometre-scale seismically slow layer atop the core-mantle boundary from normal modes

Geodynamic modelling and seismic studies have highlighted the possibility that a thin layer of low seismic velocities, potentially molten, may sit atop the core-mantle boundary but has thus far eluded detection. In this study we employ normal modes, an independent data type to body waves, to assess the visibility of a seismically slow layer atop the core-mantle boundary to normal mode centre frequencies. Using forward modelling and a dataset of 353 normal mode observations we find that some centre frequencies are sensitive to one-dimensional kilometre-scale structure at the core-mantle boundary. Furthermore, a global slow and dense layer 1 - 3 km thick is better-fitting than no layer. The well-fitting parameter space is broad with a wide range of possible seismic parameters, which precludes inferring a possible composition or phase. Our methodology cannot uniquely detect a layer in the Earth but one should be considered possible and accounted for in future studies

Seismically determined elastic parameters for Earth’s outer core

Turbulent convection of the liquid iron alloy outer core generates Earth’s magnetic field and supplies heat to the mantle. The exact composition of the iron alloy is fundamentally linked to the processes powering the convection and can be constrained by its seismic properties. Discrepancies between seismic models determined using body waves and normal modes show that these properties are not yet fully agreed upon. In addition, technical challenges in experimentally measuring the equation-of-state (EoS) parameters of liquid iron alloys at high pressures and temperatures further complicate compositional inferences. We directly infer EoS parameters describing Earth’s outer core from normal mode center frequency observations and present the resulting Elastic Parameters of the Outer Core (EPOC) seismic model. Unlike alternative seismic models, ours requires only three parameters and guarantees physically realistic behavior with increasing pressure for a well-mixed homogeneous material along an isentrope, consistent with the outer core’s condition. We show that EPOC predicts available normal mode frequencies better than the Preliminary Reference Earth Model (PREM) while also being more consistent with body wave–derived models, eliminating a long-standing discrepancy. The velocity at the top of the outer core is lower, and increases with depth more steeply, in EPOC than in PREM, while the density in EPOC is higher than that in PREM across the outer core. The steeper profiles and higher density imply that the outer core comprises a lighter but more compressible alloy than that inferred for PREM. Furthermore, EPOC’s steeper velocity gradient explains differential SmKS body wave travel times better than previous one-dimensional global models, without requiring an anomalously slow ~90- to 450-km-thick layer at the top of the outer core