scholarly journals Triplicated P-wave measurements for waveform tomography of the mantle transition zone

Solid Earth ◽  
2012 ◽  
Vol 3 (2) ◽  
pp. 339-354 ◽  
Author(s):  
S. C. Stähler ◽  
K. Sigloch ◽  
T. Nissen-Meyer

Abstract. Triplicated body waves sample the mantle transition zone more extensively than any other wave type, and interact strongly with the discontinuities at 410 km and 660 km. Since the seismograms bear a strong imprint of these geodynamically interesting features, it is highly desirable to invert them for structure of the transition zone. This has rarely been attempted, due to a mismatch between the complex and band-limited data and the (ray-theoretical) modelling methods. Here we present a data processing and modelling strategy to harness such broadband seismograms for finite-frequency tomography. We include triplicated P-waves (epicentral distance range between 14 and 30°) across their entire broadband frequency range, for both deep and shallow sources. We show that is it possible to predict the complex sequence of arrivals in these seismograms, but only after a careful effort to estimate source time functions and other source parameters from data, variables that strongly influence the waveforms. Modelled and observed waveforms then yield decent cross-correlation fits, from which we measure finite-frequency traveltime anomalies. We discuss two such data sets, for North America and Europe, and conclude that their signal quality and azimuthal coverage should be adequate for tomographic inversion. In order to compute sensitivity kernels at the pertinent high body wave frequencies, we use fully numerical forward modelling of the seismic wavefield through a spherically symmetric Earth.

2012 ◽  
Vol 4 (2) ◽  
pp. 783-821 ◽  
Author(s):  
S. C. Stähler ◽  
K. Sigloch ◽  
T. Nissen-Meyer

Abstract. Triplicated body waves sample the mantle transition zone more extensively than any other wave type, and interact strongly with the discontinuities at 410 km and 660 km. Since the seismograms bear a strong imprint of these geodynamically interesting features, it is highly desirable to invert them for structure of the transition zone. This has rarely been attemped, due to the mismatch between the complex and bandlimited data and the (ray-theoretical) modeling methods. Here we present a data processing and modeling strategy to harness such broadband seismograms for finite-frequency tomography. We include triplicated P-waves (epicentral distance range between 14 and 30°) across their entire broadband frequency range, for both deep and shallow sources. We show that it is possible to predict the complex sequence of arrivals in these seismograms, but only after a careful effort to estimate source time functions and other source parameters from data, variables that strongly influence the waveforms. Modeled and observed waveforms then yield decent cross-correlation fits, from which we measure finite-frequency traveltime anomalies. We discuss two such data sets, for North America and Europe, and conclude that their signal quality and azimuthal coverage should be adequate for tomographic inversion. In order to compute sensitivity kernels at the pertinent high body-wave frequencies, we use fully numerical forward modelling of the seismic wavefield through a spherically symmetric earth.


2021 ◽  
Author(s):  
Felix Bissig ◽  
Amir Khan ◽  
Domenico Giardini

<p>The mantle transition zone (MTZ) is bounded by seismic discontinuities at average depths of 410 km and 660 km, which are generally associated with major mantle mineral transformations. A body wave impinging from above on these discontinuities develops a refracted and reflected branch, leading to multiple arrivals of the same wavetype within a short time window. These so-called triplicated body waves are observed at regional epicentral distances (15-30°) and carry information on MTZ structure due to their strong interaction with the 410 km and 660 km discontinuities. Careful data selection and processing as well as the assessment of source parameters are necessary steps in obtaining a high quality triplication data set. In this study, we consider recordings of events in Central America at permanent and transportable USArray stations, which are inverted for mantle structure. Our methodology is based on a joint consideration of mineral physics and seismic data in a probabilistic inversion framework and allows for determination of mantle thermo-chemical and seismic velocity structure. We present constraints on the mantle structure underneath the Gulf of Mexico.</p>


1983 ◽  
Vol 73 (2) ◽  
pp. 419-434
Author(s):  
Jeffery S. Barker ◽  
Charles A. Langston

abstract Teleseismic P-wave first motions for the M ≧ 6 earthquakes near Mammoth Lakes, California, are inconsistent with the vertical strike-slip mechanisms determined from local and regional P-wave first motions. Combining these data sets allows three possible mechanisms: a north-striking, east-dipping strike-slip fault; a NE-striking oblique fault; and a NNW-striking normal fault. Inversion of long-period teleseismic P and SH waves for the events of 25 May 1980 (1633 UTC) and 27 May 1980 (1450 UTC) yields moment tensors with large non-double-couple components. The moment tensor for the first event may be decomposed into a major double couple with strike = 18°, dip = 61°, and rake = −15°, and a minor double couple with strike = 303°, dip = 43°, and rake = 224°. A similar decomposition for the last event yields strike = 25°, dip = 65°, rake = −6°, and strike = 312°, dip = 37°, and rake = 232°. Although the inversions were performed on only a few teleseismic body waves, the radiation patterns of the moment tensors are consistent with most of the P-wave first motion polarities at local, regional, and teleseismic distances. The stress axes inferred from the moment tensors are consistent with N65°E extension determined by geodetic measurements by Savage et al. (1981). Seismic moments computed from the moment tensors are 1.87 × 1025 dyne-cm for the 25 May 1980 (1633 UTC) event and 1.03 × 1025 dyne-cm for the 27 May 1980 (1450 UTC) event. The non-double-couple aspect of the moment tensors and the inability to obtain a convergent solution for the 25 May 1980 (1944 UTC) event may indicate that the assumptions of a point source and plane-layered structure implicit in the moment tensor inversion are not entirely valid for the Mammoth Lakes earthquakes.


Geophysics ◽  
2014 ◽  
Vol 79 (1) ◽  
pp. R1-R12 ◽  
Author(s):  
E. Diego Mercerat ◽  
Guust Nolet ◽  
Christophe Zaroli

We evaluated a comprehensive numerical experiment of finite-frequency tomography with ray-based (“banana-doughnut”) kernels that tested all aspects of this method, starting from the generation of seismograms in a 3D model, the window selection, and the crosscorrelation with seismograms predicted for a background model, to the final regularized inversion. In particular, we tested if the quasilinearity of crosscorrelation delays allowed us to forego multiple (linearized) iterations in the case of strong reverberations characterizing multiple scattering and the gain in resolution that can be obtained by observing body-wave dispersion. Contrary to onset times, traveltimes observed by crosscorrelation allowed us to exploit energy arriving later in the time window centered in the P-wave or any other indentifiable ray arrival, either scattered from, or diffracted around, lateral heterogeneities. We tested using seismograms calculated by the spectral element method in a cross-borehole experiment conducted in a 3D checkerboard cube. The use of multiple frequency bands allowed us to estimate body-wave dispersion caused by diffraction effects. The large velocity contrast (10%) and the regularity of the checkerboard pattern caused severe reverberations that arrived late in the crosscorrelation windows. Nevertheless, the model resulting from the inversion with a data fit with reduced [Formula: see text] resulted in an excellent correspondence with the input model and allowed for a complete validation of the linearizations that lay at the basis of the theory. The use of multiple frequencies led to a significant increase in resolution. Moreover, we evaluated a case in which the sign of the anomalies in the checkerboard was systematically reversed in the ray-theoretical solution, a clear demonstration of the reality of the “doughnut-hole” effect. The experiment validated finite-frequency theory and disqualified ray-theoretical inversions of crosscorrelation delay times.


1974 ◽  
Vol 64 (6) ◽  
pp. 1887-1899
Author(s):  
George A. McMechan ◽  
Warren G. Workman

abstract The observed behavior of P-wave relative amplitudes, as a function of epicentral distance, between 10° and 35°, is controlled primarily by the velocity-depth structure of the upper mantle. P-wave synthetic seismograms calculated by the new quantized ray theory technique are used to determine theoretical log (A/T) versus log Δ curves from a number of upper mantle models. Maximum amplitude arrivals show less model dependence than the first arrivals in the same wave trains, and hence are more consistent magnitude indicators for regions where the upper mantle structure is poorly known. Log (A/T) versus log Δ curves vary considerably, but predictably, from model to model. This model-dependent variation can account for a major part of the large standard deviations usually associated with the calculation of magnitudes from body waves.


1965 ◽  
Vol 55 (2) ◽  
pp. 203-235 ◽  
Author(s):  
Ari Ben-Menahem ◽  
Stewart W. Smith ◽  
Ta-Liang Teng

Abstract The well-known first motion method of Nakano and Byerly is extended, generalized and combined with recent new ideas in body wave theory in order to set up a routine procedure for extracting source parameters from spectral analysis of isolated P and S pulses recorded at a net of standardized stations around a non-shallow source. The method consists of compensating the observed spectrums for instrumental and propagational effects. A combined study of the resulting radiation patterns, initial phases, and the initial amplitudes will render information regarding the spatial and temporal nature of deep and intermediate earthquake sources as seen through the spectral window of 10-100 seconds. The shorter periods can be used for source studies only if an accurate station correction is available.


1974 ◽  
Vol 64 (3-1) ◽  
pp. 571-579
Author(s):  
R. K. Dube ◽  
J. C. Bhayana

abstract Crustal structure in the Gangetic Plains of India has been investigated using body-wave data of earthquakes. A three-layered crust has been interpreted, consisting of a top layer of 2.7-km/sec P-wave velocity and 3.7-km thickness, an intermediate layer of 5.64-km/sec velocity and 15.2-km thickness and a bottom layer of 6.49-km/sec velocity and 21.4-km thickness. The average depth of the M discontinuity obtained is 40.3 km. The shear-wave velocities for the Sg, S*, Sn phases are 3.45, 3.85 and 4.61 km/sec, respectively. The velocities of both P and S waves are lower than those obtained for the Peninsular Shield of India.


1984 ◽  
Vol 74 (2) ◽  
pp. 361-376
Author(s):  
John Boatwright ◽  
Jon B. Fletcher

Abstract Seventy-three digitally recorded body waves from nine multiply recorded small earthquakes in Monticello, South Carolina, are analyzed to estimate the energy radiated in P and S waves. Assuming Qα = Qβ = 300, the body-wave spectra are corrected for attenuation in the frequency domain, and the velocity power spectra are integrated over frequency to estimate the radiated energy flux. Focal mechanisms determined for the events by fitting the observed displacement pulse areas are used to correct for the radiation patterns. Averaging the results from the nine events gives 27.3 ± 3.3 for the ratio of the S-wave energy to the P-wave energy using 0.5 〈Fi〉 as a lower bound for the radiation pattern corrections, and 23.7 ± 3.0 using no correction for the focal mechanisms. The average shift between the P-wave corner frequency and the S-wave corner frequency, 1.24 ± 0.22, gives the ratio 13.7 ± 7.3. The substantially higher values obtained from the integral technique implies that the P waves in this data set are depleted in energy relative to the S waves. Cursory inspection of the body-wave arrivals suggests that this enervation results from an anomalous site response at two of the stations. Using the ratio of the P-wave moments to the S-wave moments to correct the two integral estimates gives 16.7 and 14.4 for the ratio of the S-wave energy to the P-wave energy.


Sign in / Sign up

Export Citation Format

Share Document