The rupture process of the Peru intermediate and deep earthquakes

Author(s):  
Carolina López Sánchez ◽  
Elisa Buforn ◽  
Maurizio Mattesini ◽  
Hernando Tavera

<p>The seismicity of Peru is associated with the subduction process of the Nazca plate under South America and characterized by the occurrence of shallow, intermediate and deep earthquakes. In this study, we focus our attention on the rupture process of earthquakes (Mw>6.0 ) that occurred during the period 2018-2019 at intermediate depth (50<h<200 km) and deep depth (500<h<700 km). Focal mechanisms have been estimated from slip inversion of body waves at teleseismic distances (Kikuchi and Kanamori, 1991). We investigate possible differences in the moment rate functions at different focal depths using our results and those provided by SCARDEC database. Furthermore, an estimation of the radiated seismic energy (E<sub>R</sub>) was provided from the direct integration of the velocity P wave recorded at teleseismic and regional distances, getting values between 10<sup>15 </sup>to 10<sup>16</sup> J. The data were corrected by geometrical spreading, anelastic attenuation, and free surface effect. These results are interpreted in terms of the seismotectonics of the region</p>

2021 ◽  
Author(s):  
Elisa Buforn ◽  
Carmen Pro ◽  
Hernando Tavera ◽  
Agustin Udias ◽  
Maurizio Mattesini

<p>We analyze the differences in the rupture process for twelve very deep earthquakes (h>500 km) at the Peruvian-Brazilian subduction zone. These earthquakes are produced by the contact between the Nazca and the South America Plates. We have estimated the focal mechanism from teleseismic waveforms, using the slip inversion over the rupture plane, testing rupture velocities ranging from 2.5 km/s to 4.5 km/s, and analyzing the slip distribution for each  rupture velocity. The selected 12 earthquakes have occurred in the period 1994- 2016, with magnitudes between 5.9 and 8.2 and focal depth 500- 700 km. They can be separated in two groups attending to their epicentral location. The first group is formed by 9 events located, in the Peruvian-Brazil border, with epicenters following a NNW-SSE alignment, parallel to the trench. Their focal mechanisms present solutions of normal faulting with planes oriented in NS direction, dipping about 45 degrees and with vertical pressure axis. The second group is formed by three earthquakes located to the south of the first group in northern Bolivia. Their mechanisms show dip-slip motion with a near vertical plane, oriented in NW-SE direction and the pressure axis dipping 45º to the NE. The moment rate functions correspond to single ruptures with time durations from 6s to 12s, with the exception of the large 1994 Bolivian earthquake (Mw = 8.2) which presents a complex and longer STF. The different mechanisms for the two groups of earthquakes confirm the different dip of the subducting Nazca plate at the two areas, with the steeper part at the southern one.  </p>


2019 ◽  
Vol 5 (5) ◽  
pp. eaav2032 ◽  
Author(s):  
Diego Melgar ◽  
Gavin P. Hayes

Whether earthquakes of different sizes are distinguishable early in their rupture process is a subject of debate. Studies have shown that the frequency content of radiated seismic energy in the first seconds of earthquakes scales with magnitude, implying determinism. Other studies have shown that recordings of ground displacement from small to moderate-sized earthquakes are indistinguishable, implying a universal early rupture process. Regardless of how earthquakes start, events of different sizes must be distinguishable at some point. If that difference occurs before the rupture duration of the smaller event, this implies some level of determinism. We show through analysis of a database of source time functions and near-source displacement records that, after an initiation phase, ruptures of M7 to M9 earthquakes organize into a slip pulse, the kinematic properties of which scale with magnitude. Hence, early in the rupture process—after about 10 s—large and very large earthquakes can be distinguished.


1995 ◽  
Vol 85 (3) ◽  
pp. 716-735 ◽  
Author(s):  
John F. Cassidy ◽  
Garry C. Rogers

Abstract On 6 April 1992, a magnitude 6.8 (MS) earthquake occurred in the triple-junction region at the northern end of the Cascadia subduction zone. This was the largest earthquake in at least 75 yr to occur along the 110-km-long Revere-Dellwood-Wilson (RDW) transform fault and the first large earthquake in this region recorded by modern broadband digital seismic networks. It thus provides an opportunity to examine the rupture process along a young (<2 Ma) oceanic transform fault and to gain better insight into the tectonics of this triple-junction region. We have investigated the source parameters and the rupture process of this earthquake by modeling broadband body waves and long-period surface waves and by accurately locating the mainshock and the first 10 days of aftershocks using a well-located “calibration” event recorded during an ocean-bottom seismometer survey. Analysis of P and SH waveforms reveals that this was a complex rupture sequence consisting of three strike-slip subevents in 12 sec. The initial rupture occurred 5 to 6 km to the SW of the seafloor trace of the RDW fault at 50.55° N, 130.46° W. The dominant subevent occurred 2 to 3 sec later and 4.3 km beneath the seafloor trace of the RDW fault, and a third subevent occurred 5 sec later, 18 km to the NNW, suggesting a northwestward propagating rupture. The aftershock sequence extended along a 60- to 70-km-long segment of the RDW fault, with the bulk of the activity concentrated ∼30 to 40 km to the NNW of the epicenter, consistent with this interpretation. The well-constrained mechanism of the initial rupture (strike/dip/slip 339°/90°/−168°) and of the largest aftershock (165°/80°/170°) are rotated 15° to 20° clockwise relative to the seafloor trace of the RDW fault but are parallel to the Pacific/North America relative plate motion vector. In contrast, the mechanisms of the dominant subevent (326°/87°/−172°), and the long-period solution derived from surface waves aligns with the RDW fault. This suggests that small earthquakes (M < 6) in this area occur along faults that are optimally aligned with respect to the regional stress field, whereas large earthquakes, involving tens of kilometers of rupture, activate the RDW fault. For the mainshock, we estimate a seismic moment (from surface waves) of 1.0 × 1026 dyne-cm, a stress drop of 60 bars, and an average slip of 1.2 m. This represents only 21 yr of strain accumulation, implying that there is either a significant amount of aseismic slip along the RDW fault or that much of the strain accumulation manifests itself as deformation within the Dellwood and Winona blocks or along the continental margin.


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.


1983 ◽  
Vol 73 (6A) ◽  
pp. 1835-1850
Author(s):  
Robert B. Herrmann ◽  
Andrzej Kijko

Abstract The applicaton of the Nutli (1973) definition of the mbLg magnitude to instruments and wave periods other than the short-period WWSSN seismograph is examined. The basic conclusion is that the Nuttli (1973) definition is applicable to a wider range of seismic instruments if the log10(A/T) term is replaced by log10A. For consistency and precision, the notation mbLg should be applied only to magnitudes based upon 1.0 Hz observations. The mbLg magnitude definition was constrained to be consistent with teleseismic P-wave mb estimates from four Central United States earthquakes. In general, for measurements made at a frequency f, the notation mLg(f) should be used, where m L g ( f ) = 2.94 + 0.833 log ⁡ 10 ( r / 10 ) + 0.4342 γ r + log ⁡ 10 A , and r is the epicentral distance in kilometers, γ is the coefficient of anelastic attenuation, and A is the reduced ground amplitude in microns. Given its stability when estimated from different instruments, the mLg(f) magnitude is an optimum choice for an easily applied, standard magnitude scale for use in regional seismic studies.


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.


1997 ◽  
Vol 40 (1) ◽  
Author(s):  
G. Drakatos ◽  
G. Karantonis ◽  
G. N. Stavrakakis

The three-dimensional velocity structure of the crust in the Aegean sea and the surrounding regions (34.0º-42.OºN, 19.0ºE-29.0ºE) is investigated by inversion of about 10000 residuals of arrival times of P-wave from local events. The resulting velocity structure shows strong horizontal variations due to the complicated crustal structure and the variations of crustal thickness. The northern part of the region generally shows high velocities. In the inner part of the volcanic arc (Southern Aegean area), relatively low velocities are observed, suggesting a large-scale absorption of seismic energy as confirmed by the low seismicity of the region. A low velocity zone was observed along the subduction zone of the region, up to a depth of 4 km. The existence of such a zone could be due to granitic or other intrusions in the crust during the uplift of the region during Alpidic orogenesis.


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.


1973 ◽  
Vol 63 (2) ◽  
pp. 587-597
Author(s):  
Ta-Liang Teng ◽  
James P. Tung

abstract Recent observations of P′P′ and its precursors, identified as reflections from within the Earth's upper mantle, are used to examine the structure of the uppermantle discontinuities with specific reference to the density, the S velocity, and the Q variations. The Haskell-Thomson matrix method is used to generate the complex reflection spectrum, which is then Fourier synthesized for a variety of upper-mantle velocity-density and Q models. Surface displacements are obtained for the appropriate recording instrument, permitting a direct comparison with the actual seismograms. If the identifications of the P′P′ precursors are correct, our proposed method yields the following: (1) a structure of Gutenberg-Bullen A type is not likely to produce observable P′P′ upper-mantle reflections, (2) in order that a P′P′ upper-mantle reflection is strong enough to be observed, first-order density and S-velocity discontinuities together with a P-wave discontinuity are needed at a depth of about 650 km, and (3) corresponding to a given uppermantle velocity-density model, an estimate can be made of the Q in the upper mantle for short-period seismic body waves.


2021 ◽  
Author(s):  
Thomas Möller ◽  
Wolfgang Friederich

<p>Modeling waveforms of teleseismic body waves requires the solution of the seismic wave equation in the entire Earth. Since fully-numerical 3D simulations on a global scale with periods of a few seconds are far too computationally expensive, we resort to a hybrid approach in which fully-numerical 3D simulations are performed only within the target region and wave propagation through the rest of the Earth is modeled using methods that are much faster but apply only to spherically symmetric Earth models.</p><p>We present a hybrid method that uses GEMINI to compute wave fields for a spherically symmetric Earth model up to the boundaries of a regional box. The wavefield is injected at the boundaries, where wave propagation is continued using SPECFEM-Cartesian. Inside the box, local heterogeneities in the velocity distribution are allowed, which can cause scattered and reflected waves. To prevent these waves from reflecting off the edges of the box absorbing boundary conditions are specifically applied to these parts of the wavefields. They are identified as the difference between the wavefield calculated with SPECFEM at the edges and the incident wavefield.</p><p>The hybrid method is applied to a target region in and around the Alps as a test case. The region covers an area of 1800 by 1350 km centered at 46.2°N and 10.87°E and includes crust and mantle to a depth of 600 km. We compare seismograms with a period of up to ten seconds calculated with the hybrid method to those calculated using GEMINI only for identical 1D earth models. The comparison of the seismograms shows only very small differences and thus validates the hybrid method. In addition, we demonstrate the potential of the method by calculating seismograms where the 1D velocity model inside the box is replaced by a velocity model generated using P-wave traveltime tomography.</p>


Sign in / Sign up

Export Citation Format

Share Document