Title: New constraints on the 1922 Atacama, Chile, earthquake from Historical seismograms
SUMMARY We recently found the original Omori seismograms recorded at Hongo, Tokyo, of the 1922 Atacama, Chile, earthquake (MS = 8.3) in the historical seismogram archive of the Earthquake Research Institute (ERI) of the University of Tokyo. These recordings enable a quantitative investigation of long-period seismic radiation from the 1922 earthquake. We document and provide interpretation of these seismograms together with a few other seismograms from Mizusawa, Japan, Uppsala, Sweden, Strasbourg, France, Zi-ka-wei, China and De Bilt, Netherlands. The 1922 event is of significant historical interest concerning the cause of tsunami, discovery of G wave, and study of various seismic phase and first-motion data. Also, because of its spatial proximity to the 1943, 1995 and 2015 great earthquakes in Chile, the 1922 event provides useful information on similarity and variability of great earthquakes on a subduction-zone boundary. The 1922 source region, having previously ruptured in 1796 and 1819, is considered to have significant seismic hazard. The focus of this paper is to document the 1922 seismograms so that they can be used for further seismological studies on global subduction zones. Since the instrument constants of the Omori seismographs were only incompletely documented, we estimate them using the waveforms of the observed records, a calibration pulse recorded on the seismogram and the waveforms of better calibrated Uppsala Wiechert seismograms. Comparison of the Hongo Omori seismograms with those of the 1995 Antofagasta, Chile, earthquake (Mw = 8.0) and the 2015 Illapel, Chile, earthquake (Mw = 8.3) suggests that the 1922 event is similar to the 1995 and 2015 events in mechanism (i.e. on the plate boundary megathrust) and rupture characteristics (i.e. not a tsunami earthquake) with Mw = 8.6 ± 0.25. However, the initial fine scale rupture process varies significantly from event to event. The G1 and G2, and R1 and R2 of the 1922 event are comparable in amplitude, suggesting a bilateral rupture, which is uncommon for large megathrust earthquakes.  more » « less
Award ID(s):
1802364
PAR ID:
10145603
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Geophysical Journal International
Volume:
219
Issue:
1
ISSN:
0956-540X
Page Range / eLocation ID:
645 to 661
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. On 29 July 2025, an Mw 8.8 earthquake struck off Kamchatka, Russia, generating a Pacific-wide tsunami and marking the largest earthquake since the launch of the surface water and ocean topography (SWOT) satellite in 2022. We analyze tsunami observations from SWOT together with three nearby deep-ocean assessment and reporting of tsunamis (DART) buoys to resolve the source of the event. SWOT provided the first high-resolution spaceborne track of a great subduction-zone tsunami, capturing waveforms that reveal complex propagation, dispersion, and scattering. Inversion of the DART time series using Gaussian unit sources shows that the rupture extended ∼400 km along strike, with peak uplift of ∼4 m, significantly different from the published finite-fault model. A blended source that combines the DART-inverted uplift with subsidence from the seismic–geodetic model best matches both datasets and reproduces the SWOT observations. Comparison with reconstructions of the 1952 Mw 9.0 Kamchatka earthquake indicates that the 2025 rupture likely reactivated significant portions of the megathrust that broke in 1952 but occurred farther down-dip and with little to no near-trench slip, consistent with its smaller tsunami impact. These findings highlight the hazard implications of short recurrence intervals of great earthquakes and show how rupture style governs tsunami severity. They also demonstrate the value of satellite altimetry for improving tsunami source characterization, post-event forecasting, and understanding of hydrodynamic processes. 
    more » « less
  2. SUMMARY Observations of historical tsunami earthquakes reveal that ruptures of these earthquakes propagate slowly at shallow depth with longer duration, depletion in high-frequency radiation and larger discrepancy of Mw–Ms than ordinary megathrust earthquakes. They can effectively generate tsunami and lead to huge damage to regional populated areas near the coast. In this study, we use a recently developed dynamic earthquake simulator to explore tsunami earthquake generation from a physics-based modelling point of view. We build a shallow-dipping subduction zone model in which locally locked, unstable patches (asperities) are distributed on a conditionally stable subduction interface at shallow depth. The dynamic earthquake simulator captures both quasi-static and dynamic processes of earthquake cycles. We find that earthquakes can nucleate on these asperities and propagate into the surrounding conditionally stable zone at slow speeds, generating tsunami earthquakes. A high normal stress asperity, representing a subducted seamount, can act as an asperity in some events but as a barrier in other events over multiple earthquake cycles. Low normal stress asperities typically act as asperities in tsunami earthquakes. The degree of velocity-weakening in the conditionally stable zone, which may sustain rupture at different speeds or stop rupture, is critical for tsunami earthquake generation and affects its recurrence interval. Distributed asperities may rupture in isolated events separated by tens of years, or in a sequence of events separated by hours to days, or in one large event in a cascade fashion, demonstrating complex interactions among them. The recurrence interval on a high normal stress asperity is much larger than that on low normal stress asperities. These modelling results shed lights on the observations from historical tsunami earthquakes, including the 1994 and 2006 Java tsunami earthquakes and 2010 Mentawai tsunami earthquake. 
    more » « less
  3. So far in this century, six very large–magnitude earthquakes ( M W ≥ 7.8) have ruptured separate portions of the subduction zone plate boundary of western South America along Ecuador, Peru, and Chile. Each source region had last experienced a very large earthquake from 74 to 261 y earlier. This history led to their designation in advance as seismic gaps with potential to host future large earthquakes. Deployments of geodetic and seismic monitoring instruments in several of the seismic gaps enhanced resolution of the subsequent faulting processes, revealing preevent patterns of geodetic slip deficit accumulation and heterogeneous coseismic slip on the megathrust fault. Localized regions of large slip, or asperities, appear to have influenced variability in how each source region ruptured relative to prior events, as repeated ruptures have had similar, but not identical slip distributions. We consider updated perspectives of seismic gaps, asperities, and geodetic locking to assess current very large earthquake hazard along the South American subduction zone, noting regions of particular concern in northern Ecuador and Colombia (1958/1906 rupture zone), southeastern Peru (southeasternmost 1868 rupture zone), north Chile (1877 rupture zone), and north-central Chile (1922 rupture zone) that have large geodetic slip deficit measurements and long intervals (from 64 to 154 y) since prior large events have struck those regions. Expanded geophysical measurements onshore and offshore in these seismic gaps may provide critical information about the strain cycle and fault stress buildup late in the seismic cycle in advance of the future great earthquakes that will eventually strike each region. 
    more » « less
  4. Long‐period ground motions from large (Mw≥7.0) subduction‐zone earthquakes are a real threat for large‐scale human‐made structures. The Nankai subduction zone, Japan, is expected to host a major megathrust earthquake in the near future and has therefore been instrumented with offshore and onshore permanent seismic networks. We use the ambient seismic field continuously recorded at these stations to simulate the long‐period (4–10 s) ground motions from past and future potential offshore earthquakes. First, we compute impulse response functions (IRFs) between an ocean‐bottom seismometer of the Dense Oceanfloor Network System for Earthquakes and Tsunamis (DONET) network, which is located offshore on the accretionary wedge, and 60 onshore Hi‐net stations using seismic interferometry by deconvolution. As this technique only preserves the relative amplitude information of the IRFs, we use a moderate Mw 5.5 event to calibrate the amplitudes to absolute levels. After calibration, the IRFs are used together with a uniform stress‐drop source model to simulate the long‐period ground motions of the 2004 Mw 7.2 intraplate earthquake. For both events, the residuals of the 5\\% damped spectral acceleration (SA) computed from the horizontal and vertical components of the observed and simulated waveforms exhibit almost no bias and acceptable uncertainties. We also compare the observed SA values of the Mw 7.2 event to those from the subduction‐zone BC Hydro ground‐motion model (GMM) and find that our simulations perform better than the model. Finally, we simulate the long‐period ground motions of a hypothetical Mw 8.0 subduction earthquake that could occur along the Nankai trough. For this event, our simulations generally exhibit stronger long‐period ground motions than those predicted by the BC Hydro GMM. This study suggests that the ambient seismic field recorded by the ever‐increasing number of ocean‐bottom seismometers can be used to simulate the long‐period ground motions from large megathrust earthquakes. 
    more » « less
  5. Abstract On 19 September 2022, a major earthquake struck the northwestern Michoacán segment along the Mexican subduction zone. A slip model is obtained that satisfactorily explains geodetic, teleseismic, and tsunami observations of the 2022 event. The preferred model has a compact large-slip patch that extends up-dip and northwestward from the hypocenter and directly overlaps a 1973 Mw 7.6 rupture. Slip is concentrated offshore and below the coast at depths from 10 to 30 km with a peak value of ∼2.9 m, and there is no detected coseismic slip near the trench. The total seismic moment is 3.1×1020  N·m (Mw 7.6), 72% of which is concentrated in the first 30 s. Most aftershocks are distributed in an up-dip area of the mainshock that has small coseismic slip, suggesting near-complete strain release in the large-slip patch. Teleseismic P waveforms of the 2022 and 1973 earthquakes are similar in duration and complexity with high cross-correlation coefficients of 0.68–0.98 for long P to PP signal time windows, indicating that the 2022 earthquake is a quasi-repeat of the 1973 earthquake, possibly indicating persistent frictional properties. Both the events produced more complex P waveforms than comparable size events along Guerrero and Oaxaca, reflecting differences in patchy locking of the Mexican megathrust. 
    more » « less