An official website of the United States government
SUMMARY SS-precursor imaging is used to image sharp interfaces within Earth’s mantle. Current SS-precursor techniques require tightly bandpassed signals (e.g. 0.02–0.05 Hz), limiting both vertical and horizontal resolutions. Higher frequency content would allow for the detection of finer structure in and around the mantle transition zone (MTZ). Here, we present a new SS-precursor deconvolution technique based on multiple-taper correlation (MTC). We show that applying MTC to SS-precursor deconvolution can increase the frequency cut-off up to 0.5 Hz, which potentially sharpens vertical resolution to ∼10 km. Furthermore, the high-pass frequency can be lowered (≪ 0.01 Hz), allowing more long-period energy to be included in the calculation, to better constrain the signal and reduce side lobes. Our method is benchmarked on full-waveform synthetic seismograms computed via AxiSEM3D for the PREM 1-D Earth model. We apply our novel MTC-SS-precursor deconvolution to ∼7000 seismograms recorded at broad-band borehole sensors of the Global Seismographic Network with source–receiver bounce points in the North-Central Pacific Ocean. The MTZ in this region appears to be thin, which agrees with previous results. We do not observe the 520-km discontinuity in our SS-precursor estimates. Additionally, we detect a low-velocity zone above the MTZ to the north of the Hawaiian Islands that has previously been inferred from asymmetry in side lobe amplitudes. Our high-frequency analysis demonstrates this feature to be a sharp interface (≤ 10-km thickness), rather than a thick wave speed gradient.
more »
« less
Mantle Q structure from S , SS , SSS and SSSS amplitude measurements
SUMMARY The seismic quality factor (Q) of the Earth’s mantle is of great importance for the understanding of the physical and chemical properties that control mantle anelasticity. The radial structure of the Earth’s Q is less well resolved compared to its wave speed structure, and large discrepancies exist among global 1-D Q models. In this study, we build a global data set of amplitude measurements of S, SS, SSS and SSSS waves using earthquakes that occurred between 2009 and 2017 with moment magnitudes ranging from 6.5 to 8.0. Synthetic seismograms for those events are computed in a 1-D reference model PREM, and amplitude ratios between observed and synthetic seismograms are calculated in the frequency domain by spectra division, with measurement windows determined based on visual inspection of seismograms. We simulate wave propagation in a global velocity model S40RTS based on SPECFEM3D and show that the average amplitude ratio as a function of epicentral distance is not sensitive to 3-D focusing and defocusing for the source–receiver configuration of the data set. This data set includes about 5500 S and SS measurements that are not affected by mantle transition zone triplications (multiple ray paths), and those measurements are applied in linear inversions to obtain a preliminary 1-D Q model QMSI. This model reveals a high Q region in the uppermost lower mantle. While model QMSI improves the overall datafit of the entire data set, it does not fully explain SS amplitudes at short epicentral distances or the amplitudes of the SSS and SSSS waves. Using forward modelling, we modify the 1-D model QMSI iteratively to reduce the overall amplitude misfit of the entire data set. The final Q model QMSF requires a stronger and thicker high Q region at depths between 600 and 900 km. This anelastic structure indicates possible viscosity layering in the mid mantle.
more »
« less
- Award ID(s):
- 2017218
- PAR ID:
- 10368615
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 231
- Issue:
- 1
- ISSN:
- 0956-540X
- Page Range / eLocation ID:
- p. 703-716
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
SUMMARY Long-period underside SS wave reflections have been widely used to furnish global constraints on the presence and depth of mantle discontinuities and to document evidence for their origins, for example, mineral phase-transformations in the transition zone, compositional changes in the mid-mantle and dehydration-induced melting above and below the transition zone. For higher-resolution imaging, it is necessary to separate the signature of the source wavelet (SS arrival) from that of the distortion caused by the mantle reflectivity (SS precursors). Classical solutions to the general deconvolution problem include frequency-domain or time-domain deconvolution. However, these algorithms do not easily generalize when (1) the reflectivity series is of a much shorter period compared to the source wavelet, (2) the bounce point sampling is sparse or (3) the source wavelet is noisy or hard to estimate. To address these problems, we propose a new technique called SHARP-SS: Sparse High-Resolution Algorithm for Reflection Profiling with SS waves. SHARP-SS is a Bayesian deconvolution algorithm that makes minimal a-priori assumptions on the noise model, source signature and reflectivity structure. We test SHARP-SS using real data examples beneath the NoMelt Pacific Ocean region. We recover a low-velocity discontinuity at a depth of $$\sim 69 \pm 4$$ km which marks the base of the oceanic lithosphere, consistent with previous work derived from surface waves, body wave conversions, and ScS reverberations. We anticipate high-resolution fine mantle stratification imaging using SHARP-SS at locations where seismic stations are sparsely distributed.more » « less
-
SUMMARY Global variations in the propagation of fundamental-mode and overtone surface waves provide unique constraints on the low-frequency source properties and structure of the Earth’s upper mantle, transition zone and mid mantle. We construct a reference data set of multimode dispersion measurements by reconciling large and diverse catalogues of Love-wave (49.65 million) and Rayleigh-wave dispersion (177.66 million) from eight groups worldwide. The reference data set summarizes measurements of dispersion of fundamental-mode surface waves and up to six overtone branches from 44 871 earthquakes recorded on 12 222 globally distributed seismographic stations. Dispersion curves are specified at a set of reference periods between 25 and 250 s to determine propagation-phase anomalies with respect to a reference Earth model. Our procedures for reconciling data sets include: (1) controlling quality and salvaging missing metadata; (2) identifying discrepant measurements and reasons for discrepancies; (3) equalizing geographic coverage by constructing summary rays for travel-time observations and (4) constructing phase velocity maps at various wavelengths with combination of data types to evaluate inter-dataset consistency. We retrieved missing station and earthquake metadata in several legacy compilations and codified scalable formats to facilitate reproducibility, easy storage and fast input/output on high-performance-computing systems. Outliers can be attributed to cycle skipping, station polarity issues or overtone interference at specific epicentral distances. By assessing inter-dataset consistency across similar paths, we empirically quantified uncertainties in traveltime measurements. More than 95 per cent measurements of fundamental-mode dispersion are internally consistent, but agreement deteriorates for overtones especially branches 5 and 6. Systematic discrepancies between raw phase anomalies from various techniques can be attributed to discrepant theoretical approximations, reference Earth models and processing schemes. Phase-velocity variations yielded by the inversion of the summary data set are highly correlated (R ≥ 0.8) with those from the quality-controlled contributing data sets. Long-wavelength variations in fundamental-mode dispersion (50–100 s) are largely independent of the measurement technique with high correlations extending up to degree ∼25. Agreement degrades with increasing branch number and period; highly correlated structure is found only up to degree ∼10 at longer periods (T > 150 s) and up to degree ∼8 for overtones. Only 2ζ azimuthal variations in phase velocity of fundamental-mode Rayleigh waves were required by the reference data set; maps of 2ζ azimuthal variations are highly consistent between catalogues ( R = 0.6–0.8). Reference data with uncertainties are useful for improving existing measurement techniques, validating models of interior structure, calculating teleseismic data corrections in local or multiscale investigations and developing a 3-D reference Earth model.more » « less
-
SUMMARY Seismic anisotropy has been detected at many depths of the Earth, including its upper layers, the lowermost mantle and the inner core. While upper mantle seismic anisotropy is relatively straightforward to resolve, lowermost mantle anisotropy has proven to be more complicated to measure. Due to their long, horizontal ray paths along the core–mantle boundary (CMB), S waves diffracted along the CMB (Sdiff) are potentially strongly influenced by lowermost mantle anisotropy. Sdiff waves can be recorded over a large epicentral distance range and thus sample the lowermost mantle everywhere around the globe. Sdiff therefore represents a promising phase for studying lowermost mantle anisotropy; however, previous studies have pointed out some difficulties with the interpretation of differential SHdiff–SVdiff traveltimes in terms of seismic anisotropy. Here, we provide a new, comprehensive assessment of the usability of Sdiff waves to infer lowermost mantle anisotropy. Using both axisymmetric and fully 3-D global wavefield simulations, we show that there are cases in which Sdiff can reliably detect and characterize deep mantle anisotropy when measuring traditional splitting parameters (as opposed to differential traveltimes). First, we analyze isotropic effects on Sdiff polarizations, including the influence of realistic velocity structure (such as 3-D velocity heterogeneity and ultra-low velocity zones), the character of the lowermost mantle velocity gradient, mantle attenuation structure, and Earth’s Coriolis force. Secondly, we evaluate effects of seismic anisotropy in both the upper and the lowermost mantle on SHdiff waves. In particular, we investigate how SHdiff waves are split by seismic anisotropy in the upper mantle near the source and how this anisotropic signature propagates to the receiver for a variety of lowermost mantle models. We demonstrate that, in particular and predictable cases, anisotropy leads to Sdiff splitting that can be clearly distinguished from other waveform effects. These results enable us to lay out a strategy for the analysis of Sdiff splitting due to anisotropy at the base of the mantle, which includes steps to help avoid potential pitfalls, with attention paid to the initial polarization of Sdiff and the influence of source-side anisotropy. We demonstrate our Sdiff splitting method using three earthquakes that occurred beneath the Celebes Sea, measured at many transportable array stations at a suitable epicentral distance. We resolve consistent and well-constrained Sdiff splitting parameters due to lowermost mantle anisotropy beneath the northeastern Pacific Ocean.more » « less
-
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
