skip to main content

Title: Mantle Q structure from SSSSSS and SSSS amplitude measurements

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 more » 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.

« less
; ; ;
Publication Date:
Journal Name:
Geophysical Journal International
Page Range or eLocation-ID:
p. 703-716
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this

    Precisely constraining the source parameters of large earthquakes is one of the primary objectives of seismology. However, the quality of the results relies on the quality of synthetic earth response. Although earth structure is laterally heterogeneous, particularly at shallow depth, most earthquake source studies at the global scale rely on the Green's functions calculated with radially symmetric (1-D) earth structure. To avoid the impact of inaccurate Green's functions, these conventional source studies use a limited set of seismic phases, such as long-period seismic waves, broad-band P and S waves in teleseismic distances (30° < ∆ < 90°), and strong ground motion records at close-fault stations. The enriched information embedded in the broad-band seismograms recorded by global and regional networks is largely ignored, limiting the spatiotemporal resolution. Here we calculate 3-D strain Green's functions at 30 GSN stations for source regions of 9 selected global earthquakes and one earthquake-prone area (California), with frequency up to 67 mHz (15 s), using SPECFEM3D_GLOBE and the reciprocity theorem. The 3-D SEM mesh model is composed of mantle model S40RTS, crustal model CRUST2.0 and surface topography ETOPO2. We surround each target event with grids in horizontal spacing of 5 km and vertical spacing of 2.0–3.0 km, allowing usmore »to investigate not only the main shock but also the background seismicity. In total, the response at over 210 000 source points is calculated in simulation. The number of earthquakes, including different focal mechanisms, centroid depth range and tectonic background, could further increase without additional computational cost if they were properly selected to avoid overloading individual CPUs. The storage requirement can be reduced by two orders of magnitude if the output strain Green's functions are stored for periods over 15 s. We quantitatively evaluate the quality of these 3-D synthetic seismograms, which are frequency and phase dependent, for each source region using nearby aftershocks, before using them to constrain the focal mechanisms and slip distribution. Case studies show that using a 3-D earth model significantly improves the waveform similarity, agreement in amplitude and arrival time of seismic phases with the observations. The limitations of current 3-D models are still notable, dependent on seismic phases and frequency range. The 3-D synthetic seismograms cannot well match the high frequency (>40 mHz) S wave and (>20 mHz) Rayleigh wave yet. Though the mean time-shifts are close to zero, the standard deviations are notable. Careful calibration using the records of nearby better located earthquakes is still recommended to take full advantage of better waveform similarity due to the use of 3-D models. Our results indicate that it is now feasible to systematically study global large earthquakes using full 3-D earth response in a global scale.

    « less
  2. SUMMARY A fleet of autonomously drifting profiling floats equipped with hydrophones, known by their acronym mermaid, monitors worldwide seismic activity from inside the oceans. The instruments are programmed to detect and transmit acoustic pressure conversions from teleseismic P wave arrivals for use in mantle tomography. Reporting seismograms in near-real time, within hours or days after they were recorded, the instruments are not usually recovered, but if and when they are, their memory buffers can be read out. We present a unique 1-yr-long data set of sound recorded at frequencies between 0.1 and 20 Hz in the South Pacific around French Polynesia by a mermaid float that was, in fact, recovered. Using time-domain, frequency-domain and time-frequency-domain techniques to comb through the time-series, we identified signals from 213 global earthquakes known to published catalogues, with magnitudes 4.6–8.0, and at epicentral distances between 24° and 168°. The observed signals contain seismoacoustic conversions of compressional and shear waves travelling through crust, mantle and core, including P, S, Pdif, Sdif, PKIKP, SKIKS, surface waves and hydroacoustic T phases. Only 10 earthquake records had been automatically reported by the instrument—the others were deemed low-priority by the onboard processing algorithm. After removing all seismic signals from the record,more »and also those from other transient, dominantly non-seismic, sources, we are left with the infrasonic ambient noise field recorded at 1500 m depth. We relate the temporally varying noise spectral density to a time-resolved ocean-wave model, WAVEWATCH III. The noise record is extremely well explained, both in spectral shape and in temporal variability, by the interaction of oceanic surface gravity waves. These produce secondary microseisms at acoustic frequencies between 0.1 and 1 Hz according to the well-known frequency-doubling mechanism.« less

    Accurate synthetic seismic wavefields can now be computed in 3-D earth models using the spectral element method (SEM), which helps improve resolution in full waveform global tomography. However, computational costs are still a challenge. These costs can be reduced by implementing a source stacking method, in which multiple earthquake sources are simultaneously triggered in only one teleseismic SEM simulation. One drawback of this approach is the perceived loss of resolution at depth, in particular because high-amplitude fundamental mode surface waves dominate the summed waveforms, without the possibility of windowing and weighting as in conventional waveform tomography.

    This can be addressed by redefining the cost-function and computing the cross-correlation wavefield between pairs of stations before each inversion iteration. While the Green’s function between the two stations is not reconstructed as well as in the case of ambient noise tomography, where sources are distributed more uniformly around the globe, this is not a drawback, since the same processing is applied to the 3-D synthetics and to the data, and the source parameters are known to a good approximation. By doing so, we can separate time windows with large energy arrivals corresponding to fundamental mode surface waves. This opens the possibility ofmore »designing a weighting scheme to bring out the contribution of overtones and body waves. It also makes it possible to balance the contributions of frequently sampled paths versus rarely sampled ones, as in more conventional tomography.

    Here we present the results of proof of concept testing of such an approach for a synthetic 3-component long period waveform data set (periods longer than 60 s), computed for 273 globally distributed events in a simple toy 3-D radially anisotropic upper mantle model which contains shear wave anomalies at different scales. We compare the results of inversion of 10 000 s long stacked time-series, starting from a 1-D model, using source stacked waveforms and station-pair cross-correlations of these stacked waveforms in the definition of the cost function. We compute the gradient and the Hessian using normal mode perturbation theory, which avoids the problem of cross-talk encountered when forming the gradient using an adjoint approach. We perform inversions with and without realistic noise added and show that the model can be recovered equally well using one or the other cost function.

    The proposed approach is computationally very efficient. While application to more realistic synthetic data sets is beyond the scope of this paper, as well as to real data, since that requires additional steps to account for such issues as missing data, we illustrate how this methodology can help inform first order questions such as model resolution in the presence of noise, and trade-offs between different physical parameters (anisotropy, attenuation, crustal structure, etc.) that would be computationally very costly to address adequately, when using conventional full waveform tomography based on single-event wavefield computations.

    « less
  4. 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,more »could provide coverage of roughly 85% by surface area of the CMB.« less

    Determinations of seismic anisotropy, or the dependence of seismic wave velocities on the polarization or propagation direction of the wave, can allow for inferences on the style of deformation and the patterns of flow in the Earth’s interior. While it is relatively straightforward to resolve seismic anisotropy in the uppermost mantle directly beneath a seismic station, measurements of deep mantle anisotropy are more challenging. This is due in large part to the fact that measurements of anisotropy in the deep mantle are typically blurred by the potential influence of upper mantle and/or crustal anisotropy beneath a seismic station. Several shear wave splitting techniques are commonly used that attempt resolve seismic anisotropy in deep mantle by considering the presence of multiple anisotropic layers along a raypath. Examples include source-side S-wave splitting, which is used to characterize anisotropy in the deep upper mantle and mantle transition zone beneath subduction zones, and differential S-ScS and differential SKS-SKKS splitting, which are used to study anisotropy in the D″ layer at the base of the mantle. Each of these methods has a series of assumptions built into them that allow for the consideration of multiple regions of anisotropy. In this work, we systematically assess themore »accuracy of these assumptions. To do this, we conduct global wavefield modelling using the spectral element solver AxiSEM3D. We compute synthetic seismograms for earth models that include seismic anisotropy at the periods relevant for shear wave splitting measurements (down to 5 s). We apply shear wave splitting algorithms to our synthetic seismograms and analyse whether the assumptions that underpin common measurement techniques are adequate, and whether these techniques can correctly resolve the anisotropy incorporated in our models. Our simulations reveal some inaccuracies and limitations of reliability in various methods. Specifically, explicit corrections for upper mantle anisotropy, which are often used in source-side direct S splitting and S-ScS differential splitting, are typically reliable for the fast polarization direction ϕ but not always for the time lag δt, and their accuracy depends on the details of the upper mantle elastic tensor. We find that several of the assumptions that underpin the S-ScS differential splitting technique are inaccurate under certain conditions, and we suggest modifications to traditional S-ScS differential splitting approaches that lead to improved reliability. We investigate the reliability of differential SKS-SKKS splitting intensity measurements as an indicator for lowermost mantle anisotropy and find that the assumptions built into the splitting intensity formula can break down for strong splitting cases. We suggest some guidelines to ensure the accuracy of SKS-SKKS splitting intensity comparisons that are often used to infer lowermost mantle anisotropy. Finally, we suggest a new strategy to detect lowermost mantle anisotropy which does not rely on explicit upper mantle corrections and use this method to analyse the lowermost mantle beneath east Asia.

    « less