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.
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.
more » « less- PAR ID:
- 10537647
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 239
- Issue:
- 1
- ISSN:
- 0956-540X
- Format(s):
- Medium: X Size: p. 591-620
- Size(s):
- p. 591-620
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
SUMMARY -
more » « less
Abstract Ultralow‐velocity zones (ULVZs) have been studied using a variety of seismic phases; however, their physical origin is still poorly understood. Short period ScP waveforms are extensively used to infer ULVZ properties because they may be sensitive to all ULVZ elastic moduli and thickness. However, ScP waveforms are additionally complicated by the effects of path attenuation, coherent noise, and source complexity. To address these complications, we developed a hierarchical Bayesian inversion method that allows us to invert ScP waveforms from multiple events simultaneously and accounts for path attenuation and correlated noise. The inversion method is tested with synthetic predictions which show that the inclusion of attenuation is imperative to recover ULVZ parameters accurately and that the ULVZ thickness and S‐wave velocity decrease are most reliably recovered. Utilizing multiple events simultaneously reduces the effects of coherent noise and source time function complexity, which in turn allows for the inclusion of more data to be used in the analyses. We next applied the method to ScP data recorded in Australia for 291 events that sample the core‐mantle boundary beneath the Coral Sea. Our results indicate, on average, ∼12‐km thick ULVZ with ∼14% reduction in S‐wave velocity across the region, but there is a greater variability in ULVZ properties in the south than that in the north of the sampled region. P‐wave velocity reductions and density perturbations are mostly below 10%. These ScP data show more than one ScP post‐cursor in some areas which may indicate complex 3‐D ULVZ structures.
-
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
-
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
-
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
