The January 2019 (Mw 6.7) Coquimbo Earthquake: Insights from a Seismic Sequence within the Nazca Plate

2019 ◽  
Author(s):  
Sergio Ruiz ◽  
Jean‐Baptiste Ammirati ◽  
Felipe Leyton ◽  
Leoncio Cabrera ◽  
Bertrand Potin ◽  
...  
Keyword(s):  
Moment Tensor ◽  
High Stress ◽  
Normal Fault ◽  
Seismic Source ◽  
Velocity Model ◽  
Seismic Sequence ◽  
Ground Acceleration ◽  
Nazca Plate ◽  
3D Velocity Model ◽  
Macroseismic Intensities

ABSTRACT On 20 January 2019, the Chilean cities of Coquimbo and La Serena were shaken by an intraplate earthquake of Mw 6.7 located at 70 km depth. High peak ground acceleration values and macroseismic intensities were reported. The mainshock was followed by more than 150 aftershocks higher than ML 2.5, a seismic sequence completely recorded by local stations. Using a 3D velocity model, we precisely located the seismicity. The aftershocks were located some 20 km above and shifted from the mainshock but still inside the Nazca plate. We also performed moment tensor inversion of nine events obtaining mostly normal‐fault focal mechanisms and kinematic inversions using the elliptical‐patch approach. We found that the mainshock broke an approximated zone of 6 km by 8 km, propagated upward in the northwest direction and away from the aftershock area. The rupture inverted from accelerograms containing up to 1 Hz was characterized with a high stress drop of 7.51 MPa and a short seismic source time function of only 3 s duration.

Aftershock Analysis of the 2018 Mw 7.1 Anchorage, Alaska, Earthquake: Relocations and Regional Moment Tensors

2019 ◽  
Vol 91 (1) ◽  
pp. 114-125 ◽  
Author(s):  
Natalia A. Ruppert ◽  
Avinash Nayak ◽  
Clifford Thurber ◽  
Cole Richards
Keyword(s):  
Moment Tensor ◽  
Hanging Wall ◽  
Velocity Model ◽  
Aftershock Sequence ◽  
Fault Rupture ◽  
Moment Tensors ◽  
3D Velocity Model ◽  
The Moment ◽  
The Relationship ◽  
Alaska Earthquake

Abstract The 30 November 2018 magnitude 7.1 Anchorage earthquake occurred as a result of normal faulting within the lithosphere of subducted Yakutat slab. It was followed by a vigorous aftershock sequence with over 10,000 aftershocks reported through the end of July 2019. The Alaska Earthquake Center produced a reviewed aftershock catalog with a magnitude of completeness of 1.3. This well‐recorded dataset provides a rare opportunity to study the relationship between the aftershocks and fault rupture of a major intraslab event. We use tomoDD algorithm to relocate 2038 M≥2 aftershocks with a regional 3D velocity model. The relocated aftershocks extend over a 20 km long zone between 47 and 57 km depth and are primarily confined to a high VP/VS region. Aftershocks form two clusters, a diffuse southern cluster and a steeply west‐dipping northern cluster with a gap in between where maximum slip has been inferred. We compute moment tensors for the Mw>4 aftershocks using a cut‐and‐paste method and careful selection of regional broadband stations. The moment tensor solutions do not exhibit significant variability or systematic differences between the northern and southern clusters and, on average, agree well with the mainshock fault‐plane parameters. We propose that the mainshock rupture initiated in the Yakutat lower crust or uppermost mantle and propagated both upward into the crust to near its top and downward into the mantle. The majority of the aftershocks are confined to the seismically active Yakutat crust and located both on and in the hanging wall of the mainshock fault rupture.

scholarly journals A repeatable seismic source for tomography at volcanoes

1999 ◽  
Vol 42 (3) ◽  
Author(s):  
U. Wegler ◽  
B. G. Lühr ◽  
A. Ratdomopurbo
Keyword(s):  
Seismic Tomography ◽  
Velocity Structure ◽  
Seismic Source ◽  
Velocity Model ◽  
Seismic Events ◽  
Seismic Signals ◽  
Improvement Project ◽  
Site Survey ◽  
3D Velocity Model ◽  
Source Signal

One major problem associated with the interpretation of seismic signals on active volcanoes is the lack of knowledge about the internal structure of the volcano. Assuming a 1D or a homogeneous instead of a 3D velocity structure leads to an erroneous localization of seismic events. In order to derive a high resolution 3D velocity model of<r>Mt. Merapi (Java) a seismic tomography experiment using active sources is planned as a part of the MERAPI (Mechanism Evaluation, Risk Assessment and Prediction Improvement) project. During a pre-site survey in August 1996 we tested a seismic source consisting of a 2.5 l airgun shot in water basins that were constructed in different flanks of the volcano. This special source, which in our case can be fired every two minutes, produces a repeatable, identical source signal. Using this source the number of receiver locations is not limited by the number of seismometers. The seismometers can be moved to various receiver locations while the source reproduces the same source signal. Additionally, at each receiver location we are able to record the identical source signal several times so that the disadvantage of the lower energy compared to an explosion source can be reduced by skipping disturbed signals and stacking several recordings.

Probabilistic Moment Tensor Inversion for Hydrocarbon-Induced Seismicity in the Groningen Gas Field, The Netherlands, Part 1: Testing

2020 ◽  
Vol 110 (5) ◽  
pp. 2095-2111 ◽  
Author(s):  
Daniela Kühn ◽  
Sebastian Heimann ◽  
Marius P. Isken ◽  
Elmer Ruigrok ◽  
Bernard Dost
Keyword(s):  
Input Data ◽  
Moment Tensor ◽  
Good Practice ◽  
Velocity Model ◽  
Moment Tensor Inversion ◽  
Phase Identification ◽  
Moment Tensors ◽  
Velocity Models ◽  
3D Velocity Model ◽  
1D Velocity Model

ABSTRACT Since 1991, induced earthquakes have been observed and linked to gas production in the Groningen field. Recorded waveforms are complex, resulting partly from a Zechstein salt layer overlying the reservoir and partly from free-surface reverberations, internal multiples, interface conversions, guided waves, and waves diving below the reservoir. Therefore, picking of polarities or amplitudes for use in moment tensor inversion is problematic, whereas phase identification may be circumvented employing full waveform techniques. Although moment tensors have become a basic tool to analyze earthquake sources, their uncertainties are rarely reported. We introduce a method for probabilistic moment tensor estimation and demonstrate its use on the basis of a single event within the Groningen field, concentrating on detailed tests of input data and inversion parameters to derive rules of good practice for moment tensor estimation of events recorded in the Groningen field. In addition to the moment tensor, event locations are provided. Hypocenters estimated simultaneously with moment tensors are often less sensitive to uncertainties in crustal structure, which is pertinent for the application to the Groningen field, because the task of relating earthquakes to specific faults hitherto suffers from a limited resolution of earthquake locations. Because of the probabilistic approach, parameter trade-offs, uncertainties, and ambiguities are mapped. In addition, the implemented bootstrap method implicitly accounts for modeling errors affecting every station and phase differently. A local 1D velocity model extracted from a full 3D velocity model yields more consistent results than other models applied previously. For all velocity models and combinations of input data tested, a shift in location of 1 km to the south is observed for the test event compared to the public catalog. A full moment tensor computed employing the local 1D velocity model features negative isotropic components and may be interpreted as normal fault and collapse at reservoir level.

The Crustal Seismicity of the Western Andean Thrust (Central Chile, 33°–34° S): Implications for Regional Tectonics and Seismic Hazard in the Santiago Area

2019 ◽  
Vol 109 (5) ◽  
pp. 1985-1999
Author(s):  
Jean-Baptiste Ammirati ◽  
Gabriel Vargas ◽  
Sofía Rebolledo ◽  
Rachel Abrahami ◽  
Bertrand Potin ◽  
...  
Keyword(s):  
Velocity Model ◽  
Fault Scarp ◽  
S Wave ◽  
The West ◽  
Nazca Plate ◽  
Study Region ◽  
Central Chile ◽  
3D Velocity Model ◽  
Wave Travel ◽  
Regional Tectonics

Abstract Most of the recorded seismicity in central Chile can be linked to the subduction of the Nazca plate. To the east, a much smaller fraction is observed at 0–30 km depths beneath the western Andean thrust. Paleoseismic studies evidenced the occurrence of at least two major earthquakes (M&gt;7) over the past 17 ka, associated with the San Ramón fault (SRF): an important tectonic feature characterizing the west Andean thrust, close the Santiago metropolitan area. To better constrain the crustal seismicity in this area, the Chilean Seismological Center (CSN) extended its permanent seismic network with seven new broadband seismometers deployed around the scarp of the SRF and farther east. The improved azimuthal distribution and reduced station spacing allowed to complete the CSN catalog with more than 900 smaller magnitude earthquakes (ML&lt;2.5) detected and located within the study region. The use of a 3D velocity model derived from P- and S-wave travel-time tomography considerably lowered the uncertainties associated with hypocentral locations. Our results show an important seismicity beneath the Principal Cordillera located at a depth of ∼10  km, and a deeper seismicity (~15 km) aligned with the main Andean thrust more to the west, parallel to the scarp of the SRF. Regional stress inversion results suggest that the seismicity of the west Andean thrust is accommodating northeast–southwest compressional stress, consistent with the convergence of the Nazca plate. Based on our improved crustal seismicity, combined with observations from previous studies, we have been able to refine the scenario of an Mw 7.5 earthquake rupturing the SRF. Ground-motion prediction results show peak ground accelerations of ∼0.8g close to the fault scarp.

scholarly journals Earthquake location in tectonic structures of the Alpine Chain: the case of the Constance Lake (Central Europe) seismic sequence

2021 ◽  
Author(s):  
Francesca D’Ajello Caracciolo ◽  
Rodolfo Console
Keyword(s):  
Central Europe ◽  
Velocity Model ◽  
Moho Depth ◽  
Seismic Sequence ◽  
Tectonic Structures ◽  
Study Region ◽  
Waveform Data ◽  
Seismic Stations ◽  
Velocity Models ◽  
Master Event

AbstractA set of four magnitude Ml ≥ 3.0 earthquakes including the magnitude Ml = 3.7 mainshock of the seismic sequence hitting the Lake Constance, Southern Germany, area in July–August 2019 was studied by means of bulletin and waveform data collected from 86 seismic stations of the Central Europe-Alpine region. The first single-event locations obtained using a uniform 1-D velocity model, and both fixed and free depths, showed residuals of the order of up ± 2.0 s, systematically affecting stations located in different areas of the study region. Namely, German stations to the northeast of the epicenters and French stations to the west exhibit negative residuals, while Italian stations located to the southeast are characterized by similarly large positive residuals. As a consequence, the epicentral coordinates were affected by a significant bias of the order of 4–5 km to the NNE. The locations were repeated applying a method that uses different velocity models for three groups of stations situated in different geological environments, obtaining more accurate locations. Moreover, the application of two methods of relative locations and joint hypocentral determination, without improving the absolute location of the master event, has shown that the sources of the four considered events are separated by distances of the order of one km both in horizontal coordinates and in depths. A particular attention has been paid to the geographical positions of the seismic stations used in the locations and their relationship with the known crustal features, such as the Moho depth and velocity anomalies in the studied region. Significant correlations between the observed travel time residuals and the crustal structure were obtained.

scholarly journals Fully probabilistic seismic source inversion – Part 1: Efficient parameterisation

2013 ◽  
Vol 5 (2) ◽  
pp. 1125-1162 ◽  
Author(s):  
S. C. Stähler ◽  
K. Sigloch
Keyword(s):  
Moment Tensor ◽  
Source Parameters ◽  
Practical Importance ◽  
Principal Component ◽  
Seismic Source ◽  
Empirical Orthogonal Functions ◽  
Time Function ◽  
Source Inversion ◽  
Source Time Function ◽  
Earthquake Parameters

Abstract. Seismic source inversion is a non-linear problem in seismology where not just the earthquake parameters themselves, but also estimates of their uncertainties are of great practical importance. Probabilistic source inversion (Bayesian inference) is very adapted to this challenge, provided that the parameter space can be chosen small enough to make Bayesian sampling computationally feasible. We propose a framework for PRobabilistic Inference of Source Mechanisms (PRISM) that parameterises and samples earthquake depth, moment tensor, and source time function efficiently by using information from previous non-Bayesian inversions. The source time function is expressed as a weighted sum of a small number of empirical orthogonal functions, which were derived from a catalogue of >1000 STFs by a principal component analysis. We use a likelihood model based on the cross-correlation misfit between observed and predicted waveforms. The resulting ensemble of solutions provides full uncertainty and covariance information for the source parameters, and permits to propagate these source uncertainties into travel time estimates used for seismic tomography. The computational effort is such that routine, global estimation of earthquake mechanisms and source time functions from teleseismic broadband waveforms is feasible.

Gaussian beam migration

Geophysics ◽  
1990 ◽  
Vol 55 (11) ◽  
pp. 1416-1428 ◽  
Author(s):  
N. Ross Hill
Keyword(s):  
Gaussian Beam ◽  
Wave Field ◽  
Seismic Source ◽  
Velocity Model ◽  
Downward Continuation ◽  
Gaussian Beams ◽  
Seismic Velocities ◽  
Migration Method ◽  
Gaussian Beam Migration ◽  
Lateral Variations

Just as synthetic seismic data can be created by expressing the wave field radiating from a seismic source as a set of Gaussian beams, recorded data can be downward continued by expressing the recorded wave field as a set of Gaussian beams emerging at the earth’s surface. In both cases, the Gaussian beam description of the seismic‐wave propagation can be advantageous when there are lateral variations in the seismic velocities. Gaussian‐beam downward continuation enables wave‐equation calculation of seismic propagation, while it retains the interpretive raypath description of this propagation. This paper describes a zero‐offset depth migration method that employs Gaussian beam downward continuation of the recorded wave field. The Gaussian‐beam migration method has advantages for imaging complex structures. Like finite‐difference migration, it is especially compatible with lateral variations in velocity, but Gaussian beam migration can image steeply dipping reflectors and will not produce unwanted reflections from structure in the velocity model. Unlike other raypath methods, Gaussian beam migration has guaranteed regular behavior at caustics and shadows. In addition, the method determines the beam spacing that ensures efficient, accurate calculations. The images produced by Gaussian beam migration are usually stable with respect to changes in beam parameters.

Strong ground motion of mine tremors: Some implications for near-source ground motion parameters

1981 ◽  
Vol 71 (1) ◽  
pp. 295-319
Author(s):  
A. McGarr ◽  
R. W. E. Green ◽  
S. M. Spottiswoode
Keyword(s):  
Ground Motion ◽  
Source Parameters ◽  
Seismic Source ◽  
Shear Velocity ◽  
Recording System ◽  
Ground Acceleration ◽  
Hypocentral Distance ◽  
Earthquake Size ◽  
Mine Tremors ◽  
Stress Drops

abstract Ground acceleration was recorded at a depth of about 3 km in the East Rand Proprietary Mines, South Africa, for tremors with −1 ≦ ML ≦ 2.6 in the hypocentral distance range 50 m &lt; R ≦ 1.6 km. The accelerograms typically had predominant frequencies of several hundred Hertz and peak accelerations, a, as high as 12 g. The peak accelerations show a dependence on magnitude, especially when expressed as dynamic shear-stress differences, defined as σ˜ = ρRa, where ρ is density. For the mine tremors, σ˜ varies from 2 to 500 bars and depends on magnitude according to log σ˜ = 1.40 + 0.38 · ML. Accelerograms for 12 events were digitized and then processed to determine velocity and, for seven events with especially good S/N, displacement and seismic source parameters. Peak ground velocities v ranged up to 6 cm/sec and show a well-defined dependence one earthquake size as measured by ML or by seismic moment, Mo. On the basis of regression fits to the mine data, with −0.76 ≦ ML ≦ 1.45, log Rv = 3.95 + 0.57 ML, where Rv is in cm2/sec, and log Rv = −4.68 + 0.49 log Mo. These regression lines agree excellently with the corresponding data for earthquakes of ML up to 6.4 or Mo to 1.4 × 1026 dyne-cm. At a given value of ML or Mo, a, at fixed R, shows considerably greater variation than v and appears to depend on the bandwidth of the recording system. The peak acceleration at small hypocentral distances is broadly consistent with ρRa = 1.14 Δτrofs/β, where Δτ is stress drop, ro is the source radius, β is shear velocity, and fs is the bandwidth of the recording system. The peak velocity data agree well with Rv = 0.57 βΔτro/μ, where μ is the modulus of rigidity; both expressions follow from Brune's model of the seismic source and were compared with data for events in the size range 5 × 1016 ≦ Mo ≦ 1.4 × 1026 dyne-cm. Measurements of the source parameters indicated that, as for earthquakes, the stress drops for the tremors range from 1 to 100 bars and show no consistent dependence on Mo down to Mo = 5 × 1016 dyne-cm.

Sign in / Sign up

Export Citation Format

Share Document