skip to main content

Title: FastTrip: A Fast MPI‐Accelerated 1D Triplication Waveform Inversion Package for Constraining Mantle Transition Zone Discontinuities
The 410‐ and 660‐km discontinuities define the top and bottom of the mantle transition zone (MTZ). The properties of these mineralogical phase transformation interfaces provide critical constraints on the dynamics, temperature, and composition of the MTZ. Triplicated body waves that bottom near these discontinuities carry rich information about them. To streamline the modeling of upper‐mantle triplications recorded at regional distances (13°–30°), we have developed a (Fast) Message Passing Interface (MPI)‐accelerated 1D (Tr)iplication Waveform (I)nversion (P)ackage (FastTrip). With triplication waveform data as input, FastTrip uses a global search method to output a set of acceptable 1D velocity models. Quantitative estimation of the model uncertainties can be further derived based on the range of acceptable models. FastTrip supports central processing unit (CPU) parallel acceleration (15,000 models within 2 hr with 100 CPUs) and is portable to other inversion problems that can be described by a relatively small number of model parameters.
; ; ; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Seismological research letters
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  1. SUMMARY The method of ScS reverberation migration is based on a ‘common reflection point’ analysis of multiple ScS reflections in the mantle transition zone (MTZ). We examine whether ray-theoretical traveltimes, slownesses and reflection points are sufficiently accurate for estimating the thickness H of the MTZ, defined by the distance between the 410- and 660-km phase transitions. First, we analyse ScS reverberations generated by 35 earthquakes and recorded at hundreds of seismic stations from the combined Arrays in China, Hi-NET in Japan and the Global Seismic Network. This analysis suggests that H varies by about 30 km and therefore that dynamic processes have modified the large-scale structure of the MTZ in eastern Asia and the western Pacific region. Second, we apply the same procedure to spectral-element synthetics for PREM and two 3-D models. One 3-D model incorporates degree-20 topography on the 410 and 660 discontinuities, otherwise preserving the PREM velocity model. The other model incorporates the degree-20 velocity heterogeneity of S20RTS and leaves the 410 and 660 flat. To optimize reflection point coverage, our synthetics were computed assuming a homogeneous grid of stations using 16 events, four of which are fictional. The resolved image using PREM synthetics resembles the PREM structure andmore »indicates that the migration approach is correct. However, ScS reverberations are not as strongly sensitive to H as predicted ray-theoretically because the migration of synthetics for a model with degree-20 topography on the 410 and 660: H varies by less than 5 km in the resolved image but 10 km in the original model. In addition, the relatively strong influence of whole-mantle shear-velocity heterogeneity is evident from the migration of synthetics for the S20RTS velocity model and the broad sensitivity kernels of ScS reverberations at a period of 15 s. A ray-theoretical approach to modelling long-period ScS traveltimes appears inaccurate, at least for continental-scale regions with relatively sparse earthquake coverage. Additional modelling and comparisons with SS precursor and receiver function results should rely on 3-D waveform simulations for a variety of structures and ultimately the implementation of full wave theory.« less

    Long-period (T > 10 s) shear wave reflections between the surface and reflecting boundaries below seismic stations are useful for studying phase transitions in the mantle transition zone (MTZ) but shear-velocity heterogeneity and finite-frequency effects complicate the interpretation of waveform stacks. We follow up on a recent study by Shearer & Buehler (hereafter SB19) of the top-side shear wave reflection Ssds as a probe for mapping the depths of the 410-km and 660-km discontinuities beneath the USArray. Like SB19, we observe that the recorded Ss410s-S and Ss660s-S traveltime differences are longer at stations in the western United States than in the central-eastern United States. The 410-km and 660-km discontinuities are about 40–50 km deeper beneath the western United States than the central-eastern United States if Ss410s-S and Ss660s-S traveltime differences are transformed to depth using a common-reflection point (CRP) mapping approach based on a 1-D seismic model (PREM in our case). However, the east-to-west deepening of the MTZ disappears in the CRP image if we account for 3-D shear wave velocity variations in the mantle according to global tomography. In addition, from spectral-element method synthetics, we find that ray theory overpredicts the traveltime delays of the reverberations. Undulations of the 410-kmmore »and 660-km discontinuities are underestimated when their wavelengths are smaller than the Fresnel zones of the wave reverberations in the MTZ. Therefore, modelling of layering in the upper mantle must be based on 3-D reference structures and accurate calculations of reverberation traveltimes.

    « less
  3. SUMMARY The detailed structure near the 410-km discontinuity provides key constraints of the dynamic interactions between the upper mantle and the lower mantle through the mantle transition zone (MTZ) via mass and heat exchange. Meanwhile, the temperature of the subducting slab, which can be derived from its fast wave speed perturbation, is critical for understanding the mantle dynamics in subduction zones where the slab enters the MTZ. Multipathing, i.e. triplicated, body waves that bottom near the MTZ carry rich information of the 410-km discontinuity structure and can be used to constrain the discontinuity depth and radial variations of wave speeds across it. In this study, we systematically analysed the trade-off between model parameters in triplication studies using synthetic examples. Specifically, we illustrated the necessity of using array-normalized amplitude. Two 1-D depth profiles of the wave speed below the Tatar Strait of Russia in the Kuril subduction zone are obtained. We have observed triplications due to both the 410-km discontinuity and the slab upper surface. And, seismic structures for these two interfaces are simultaneously inverted. Our derived 410-km discontinuity depths for the northern and southern regions are at 420$\pm $15 and 425$\pm $15 km, respectively, with no observable uplift. The slab upper surface is inverted tomore »be located about 50–70 km below the 410-km discontinuity. This location is between the depths of the 1 and 2 per cent P-wave speed perturbation contours of a regional 3-D full-waveform inversion (FWI) model, but we found twice the wave speed perturbation amplitude. A wave speed increase of 3.9–4.6 per cent within the slab, compared to 2.0–2.4 per cent from the 3-D FWI model, is necessary to fit the waveforms with the shortest period of 2 s, indicating that high-frequency waves are required to accurately resolve the detailed structures near the MTZ.« less
  4. Summary To provide constraints on a number of significant controversial issues related to the structure and dynamics of the Australian continent, we utilize P-to-S receiver functions (RFs) recorded by 182 stations to map the 410 and 660 km discontinuities (d410 and d660, respectively) bordering the mantle transition zone (MTZ). The RFs are stacked in successive circular bins with a radius of 1o under a non-plane wavefront assumption. The d410 and d660 depths obtained using the 1-D IASP91 Earth model show a systematic apparent uplifting of about 15 km for both discontinuities in central and western Australia relative to eastern Australia, as the result of higher seismic wavespeeds in the upper mantle beneath the former area. After correcting the apparent depths using the Australian Seismological Reference Model, the d410 depths beneath the West Australia Craton are depressed by ∼10 km on average relative to the normal depth of 410 km, indicating a positive thermal anomaly of 100 K at the top of the MTZ which could represent a transition from a thinner than normal MTZ beneath the Indian ocean and the normal MTZ beneath central Australia. The abnormally thick MTZ beneath eastern Australia can be adequately explained by subducted cold slabs in the MTZ. A localizedmore »normal thickness of the MTZ beneath the Newer Volcanics Province provides supporting evidence of non-mantle-plume mechanism for intraplate volcanic activities in the Australian continent.« less
  5. With the ongoing discussion of Earth structure under West Antarctica and how it relates to the extension and volcanism of the area, we explore the possibility of a hydrated or thermally perturbed mantle underneath the region. Using P-wave receiver functions, we focus on the Mantle Transition Zone (MTZ) and how its thickness fluctuates from the global average (240-260 km). Prior studies have explored the West Antarctic regions of Marie Byrd Land and the West Antarctic Rift, but we expand this to include ~3-5 years of recent, additional seismic data from the Amundsen Sea and Pine Island Bay regions. Several years of additional data from the Ronne-Fichtner Ice Shelf, Ellsworth Land, and Marie Byrd Land regions will help provide a more complete picture of the mantle transition zone. Data for this study was obtained from IRIS for earthquakes of a 5.5 magnitude or greater. We use an iterative, time domain deconvolution method, filtered with Gaussian widths of 0.5, 0.75, and 1.0. All events within their respective Gaussian filter have undergone quality check by removing waveforms that have lower than 85% fit and visually checking for clear outliers. We migrate the receiver functions to depth and stack, using both single station stackingmore »and common conversion point (CCP) stacking. We migrate the CCP stacks assuming both 1D (AK-135) and 3D velocity models throughout the region. Preliminary results from single-station stacks beneath the Thurston Island and Amundsen Sea regions suggest that the MTZ thickness is similar to the global average and the depth to the transition zone appears to be depressed, with average transition zone boundaries appearing around 430 and 680 km. If the MTZ is thinner than the global average, it may be an indication for high temperature thermal anomalies or a plume under West Antarctica that may help explain the history of extension and uplift there. These results could be useful for glacial isostatic adjustment and/or geothermal heat flux models that attempt to understand ice sheet history and stability.« less