Ultra-low velocity zones (ULVZs) have been identified as regions of extremely low velocity anomalies in the Earth's lowermost mantle using seismic observations from reflected, refracted and diffracted arrivals along the mantle side of the core–mantle boundary. Estimation of ULVZ geometrical (i.e. shape and size) and elastic (i.e. velocity and density) parameters with uncertainties is crucial in understanding the role of ULVZs in the ongoing dynamic processes within the Earth's mantle; however, these parameters are still poorly known due to uncertainties and tradeoffs of the seismically resolved ULVZ geometries and elastic parameters. Computation of synthetic waveforms for 2-D and 3-D ULVZs shapes is currently computationally feasible, but past studies utilize higher dimensional waveform modelling of mostly only low-frequency diffracted waves. Most studies focusing on high-frequency core-reflected waveforms (e.g. ScP) still use 1-D modelling approaches to determine ULVZ properties. This approach might lead to wrong results if the imaged structures have inherently 3-D geometries. This study investigates high-frequency synthetic ScP waveforms for various 2.5-D ULVZ geometries showing that additional seismic arrivals are generated even when the ScP geometrical ray path does not directly strike the location of the ULVZ. The largest amplitude additional phases in the 2.5-D models are post-cursor arrivals that are generated at the edges of the finite-length ULVZs. These newly identified ScP post-cursors can arrive within the ScsP post-cursor time window traditionally analysed in 1-D ULVZ studies. These post-cursors might then be misidentified or constructively/destructively interfere with the ScsP post-cursor, leading to incorrect estimation of ULVZ parameters. In this study we investigate the bias introduced by the 2.5-D morphologies on the 1-D estimated ULVZ elastic properties in a Bayesian waveform inversion scheme. We further expand the Bayesian method by including the data noise covariance matrix in the inversion and compare it to an autoregressive noise model that was utilized in previous studies. From the application to the observed ScP data, we find that the new approach converges faster, particularly for the inversion of data from multiple events, and the new algorithm retrieves ULVZ parameters with more realistic uncertainties. The inversion of 2.5-D synthetic ScP waveforms suggests that the retrieved ULVZ parameters can be misleading with unrealistically high confidence if we do not consider the data noise covariance matrix in the inversion. Our new approach can also retrieve the shape of a multidimensional Gaussian ULVZ if its length is 12° or longer in the great circle arc direction. However, 2.5-D synthetic waveforms show additional waveform complexity which can constructively interfere with the ScP wavefield. Hence, in many cases the estimation of ULVZ properties using 1-D forward modelling can provide incorrect ULVZ parameters. Hence previous ULVZ modelling efforts using 1-D parameter estimation methods may be incorrect.
Within the last decade, thin ultra-low velocity zone (ULVZ) layering, sitting directly on top of the core–mantle boundary (CMB), has begun to be investigated using the flip-reverse-stack (FRS) method. In this method, pre- and post-cursor arrivals that are symmetrical in time about the ScS arrival, but with opposite polarities, are stacked. This same methodology has also been applied to high velocity layering, with indications that ultra-high velocity zones (UHVZs) may also exist. Thus far, studies using the FRS technique have relied on 1-D synthetic predictions to infer material properties of ULVZs. 1-D ULVZ models predominantly show a SdS precursor that reflects off the top of the ULVZ and an ScscS reverberation within the ULVZ that arrives as a post-cursor. 1-D UHVZ models are more complex and have a different number of arrivals depending on a variety of factors including UHVZ thickness, velocity contrast, and lateral extent. 1-D modelling approaches assume that lower mantle heterogeneity is constant and continuous everywhere across the lower mantle. However, lower mantle features display lateral heterogeneity and are either finite in extent or display local thickness variations. We examine the interaction of the ScS wavefield with ULVZs and UHVZs in 2.5-D geometries of finite extent. We show that multiple additional arrivals exist that are not present in 1-D predictions. In particular, multipath ScS arrivals as well as additional post-cursor arrivals are generated. Subsequent processing by the FRS method generates complicated FRS traces with multiple peaks. Furthermore, post-cursor arrivals can be generated even when the ScS ray path does not directly strike the heterogeneity from above. Analysing these predictions for 2.5-D models using 1-D modelling techniques demonstrates that a cautious approach must be adopted in utilization and interpretation of FRS traces to determine if the ScS wavefield is interacting with a ULVZ or UHVZ through a direct strike on the top of the feature. In particular, traveltime delays or advances of the ScS arrival should be documented and symmetrical opposite polarity arrivals should be demonstrated to exist around ScS. The latter can be quantified by calculation of a time domain multiplication trace. Because multiple post-cursor arrivals are generated by finite length heterogeneities, interpretation should be confined to single layer models rather than to interpret the additional peaks as internal layering. Furthermore, strong trade-offs exist between S-wave velocity perturbation and thickness making estimations of ULVZ or UHVZ elastic parameters highly uncertain. We test our analysis methods using data from an event occurring in the Fiji-Tonga region recorded in North America. The ScS bounce points for this event sample the CMB region to the southeast of Hawaii, in a region where ULVZs have been identified in several recent studies. We see additional evidence for a ULVZ in this region centred at 14°N and 153°W with a lateral scale of at least 250 km × 360 km. Assuming a constant S-wave velocity decrease of −10 or −20 per cent with respect to the PREM model implies a ULVZ thickness of up to 16 or 9 km, respectively.
more » « less- PAR ID:
- 10542692
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 239
- Issue:
- 2
- ISSN:
- 0956-540X
- Format(s):
- Medium: X Size: p. 1038-1062
- Size(s):
- p. 1038-1062
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
SUMMARY -
Qualitative and quantitative analysis of seismic waveforms sensitive to the core–mantle boundary (CMB) region reveal the presence of ultralow-velocity zones (ULVZs) that have a strong decrease in compressional (P) and shear (S) wave velocity, and an increase in density within thin structures. However, understanding their physical origin and relation to the other large-scale structures in the lowermost mantle are limited due to an incomplete mapping of ULVZs at the CMB. The SKS and SPdKS seismic waveforms is routinely used to infer ULVZ presence, but has thus far only been used in a limited epicentral distance range. As the SKS/SPdKS wavefield interacts with a ULVZ it generates additional seismic arrivals, thus increasing the complexity of the recorded wavefield. Here, we explore utilization of the multi-scale sample entropy method to search for ULVZ structures. We investigate the feasibility of this approach through analysis of synthetic seismograms computed for PREM, 1-, 2.5-, and 3-D ULVZs as well as heterogeneous structures with a strong increase in velocity in the lowermost mantle in 1- and 2.5-D. We find that the sample entropy technique may be useful across a wide range of epicentral distances from 100° to 130°. Such an analysis, when applied to real waveforms, could provide coverage of roughly 85% by surface area of the CMB.more » « less
-
more » « less
Abstract The locations of ultralow‐velocity zones (ULVZs) at the core‐mantle boundary (CMB) have been linked to a variety of features including hot spot volcanoes and large low‐velocity province (LLVP) boundaries, yet only a small portion of the CMB region has been probed for ULVZ existence. Here we present a new map of lower mantle heterogeneity locations using a global collection of highly anomalous SPdKS recordings based on a dataset of more than 58,000 radial component seismograms, which sample 56.9% of the CMB by surface area. The inference of heterogeneity location using the SPdKS seismic phase is challenging due to source‐versus receiver‐side ambiguity. Due to this ambiguity, we conducted an inversion using the principle of parsimony. The inversion is conducted using a genetic algorithm which is repeated several thousand times in order to construct heterogeneity probability maps. This analysis reveals that at probabilities
0.5, 0.25, and 0.125 up to 1.3%, 8.2%, or 19.7% of the CMB may contain ULVZ‐like heterogeneities. These heterogeneities exist in all lower mantle settings, including both high‐ and low‐velocity regions. Additionally, we present evidence that the Samoan ULVZ may be twice as large as previously estimated, and also present evidence for the existence of additional mega‐sized ULVZs, such as a newly discovered ULVZ located to the east of the Philippines. We provide new evidence for the ULVZ east of the Philippines through an analysis of ScP records. -
The core mantle boundary (CMB), where the solid silicate mantle meets the liquid iron-nickel outer core, represents the largest density contrast on our planet, and it has long been recognized that the CMB is associated with significant structural heterogeneities. One CMB structure of particular interest are ultra low-velocity zones (ULVZs), laterally-varying, 5-50 km thick isolated patches seen in some locations just above the CMB that are associated with increased density and reduced seismic wave velocities. These variable characteristics have led to many questions regarding ULVZ origins, but less than 20% of the CMB has been surveyed for the presence of ULVZs given limited seismic coverage of the lowermost mantle. Therefore, investigations that sample the CMB with new geometries are critical to further our understanding of ULVZs and their potential connection to other deep Earth processes. The Transantarctic Mountains Northern Network (TAMNNET), a 15-station seismic array that was recently deployed in Antarctica, provides a unique dataset to further study ULVZ structure with new and unique path geometry. Core-reflected ScP phases from the TAMNNET dataset well sample the CMB in the vicinity of New Zealand in the southwestern Pacific, providing coverage between an area to the north where ULVZ structure has been previously identified and another region to the south, which shows no ULVZ evidence. This area is of particular interest because the data points sample across the boundary of the Pacific large low shear velocity province (LLSVP). The Weddell Sea region in Antarctica is also well sampled, providing new information on this area that has not been previously studied. By identifying and modeling pre- and post-cursor ScP energy, we have explored new portions of the CMB and found evidence for ULVZs in both regions. Given that ULVZs are detected within, along the edge of, and far from the Pacific LLSVP, our results may support that ULVZs are actually present everywhere along the CMB but that they are sometimes undetectable given associated methodology resolution.more » « less
-
Given limited seismic coverage of the lowermost mantle, less than one-fourth of the core-mantle boundary (CMB) has been surveyed for the presence of ultra-low velocity zones (ULVZs). Investigations that sample the CMB with new geometries are therefore important to further our understanding of ULVZ origins and their potential connection to other deep Earth processes. Using core-reflected ScP waves recorded by the recently deployed Transantarctic Mountains Northern Network in Antarctica, our study aims to expand ULVZ investigations in the southern hemisphere. Our dataset samples the CMB in the vicinity of New Zealand, providing coverage between an area to the northeast, where ULVZ structure has been previously identified, and another region to the south, where prior evidence for an ULVZ was inconclusive. This area is of particular interest because the data sample across the boundary of the Pacific Large Low Shear Velocity Province (LLSVP). The Weddell Sea region near Antarctica is also well sampled, providing new information on a region that has not been previously studied. A correlative scheme between 1-D synthetic seismograms and the observed ScP data demonstrates that ULVZs are required in both study regions. Modeling uncertainties limit our ability to definitively define ULVZ characteristics but also likely indicate more complex 3-D structure. Given that ULVZs are detected within, along the edge of, and far from the Pacific LLSVP, our results support the hypothesis that ULVZs are compositionally distinct from the surrounding mantle. ULVZs may be ubiquitous along the CMB; however, they may be thinner in many regions than can be resolved by current methods. Mantle convection currents may sweep the ULVZs into thicker piles in some areas, pushing these anomalies toward the boundaries of LLSVPs.more » « less
