skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on September 10, 2026

Title: Paleoseismic evidence of directivity for the 1976 Mw 7.5 Motagua earthquake, Guatemala
Abstract On 4 February 1976, a Mw 7.5 earthquake along the Motagua fault, Guatemala, ruptured ~230 km of the North American and Caribbean plate boundary. Today, the plate boundary remains poorly monitored, and the 1976 earthquake is still not fully understood. Here, we present seismic reflection profiles and radiometrically dated sediment core data from six lakes around the Motagua fault, together with reports of destruction and a quasi-dynamic rupture model, which show that the 1976 earthquake experienced strong directivity that impacted the distribution of shaking. The earthquake left behind a detailed record of event deposits (EDs) in five of the six study lakes. Thicker EDs are present in Lake Atitlán, near the terminus of the earthquake rupture, whereas thinner EDs were found in lakes off-axis of the rupture direction. We argue that EDs can be utilized to constrain asymmetrical distribution of shaking during earthquakes and that paleoseismic studies should consider directivity as a factor controlling the thickness of EDs.  more » « less
Award ID(s):
2325891 2038179
PAR ID:
10650357
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ;
Publisher / Repository:
Geological Society of America
Date Published:
Journal Name:
Geology
Volume:
53
Issue:
11
ISSN:
0091-7613
Page Range / eLocation ID:
971 to 976
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract A devastating magnitude 7.2 earthquake struck Southern Haiti on 14 August 2021. The earthquake caused severe damage and over 2000 casualties. Resolving the earthquake rupture process can provide critical insights into hazard mitigation. Here we use integrated seismological analyses to obtain the rupture history of the 2021 earthquake. We find the earthquake first broke a blind thrust fault and then jumped to a disconnected strike‐slip fault. Neither of the fault configurations aligns with the left‐lateral tectonic boundary between the Caribbean and North American plates. The complex multi‐fault rupture may result from the oblique plate convergence in the region, so that the initial thrust rupture is due to the boundary‐normal compression and the following strike‐slip faulting originates from the Gonâve microplate block movement, orienting SW‐NE direction. The complex rupture development of the earthquake suggests that the regional deformation is accommodated by a network of segmented faults with diverse faulting conditions. 
    more » « less
  2. Abstract The 2023 Turkey earthquake sequence involved unexpected ruptures across numerous fault segments. We present 3D dynamic rupture simulations to illuminate the complex dynamics of the earthquake doublet. Our models are constrained by observations available within days of the sequence and deliver timely, mechanically consistent explanations of the unforeseen rupture paths, diverse rupture speeds, multiple slip episodes, heterogeneous fault offsets, locally strong shaking, and fault system interactions. Our simulations link both earthquakes, matching geodetic and seismic observations and reconciling regional seismotectonics, rupture dynamics, and ground motions of a fault system represented by 10 curved dipping segments and embedded in a heterogeneous stress field. The Mw 7.8 earthquake features delayed backward branching from a steeply branching splay fault, not requiring supershear speeds. The asymmetrical dynamics of the distinct, bilateral Mw 7.7 earthquake are explained by heterogeneous fault strength, prestress orientation, fracture energy, and static stress changes from the previous earthquake. Our models explain the northward deviation of its eastern rupture and the minimal slip observed on the Sürgü fault. 3D dynamic rupture scenarios can elucidate unexpected observations shortly after major earthquakes, providing timely insights for data-driven analysis and hazard assessment toward a comprehensive, physically consistent understanding of the mechanics of multifault systems. 
    more » « less
  3. Abstract The spectra of earthquake waveforms can provide important insight into rupture processes, but the analysis and interpretation of these spectra is rarely straightforward. Here we develop a Bayesian framework that embraces the inherent data and modeling uncertainties of spectral analysis to infer key source properties. The method uses a spectral ratio approach to correct the observedS‐wave spectra of nearby earthquakes for path and site attenuation. The objective then is to solve for a joint posterior probability distribution of three source parameters—seismic moment, corner frequency, and high‐frequency falloff rate—for each earthquake in the sequence, as well as a measure of rupture directivity for select target events with good azimuthal station coverage. While computationally intensive, this technique provides a quantitative understanding of parameter tradeoffs and uncertainties and allows one to impose physical constraints through prior distributions on all source parameters, which guide the inversion when data is limited. We demonstrate the method by analyzing in detail the source properties of 14 different target events of magnitude M5 in southern California that span a wide range of tectonic regimes and fault systems. These prominent earthquakes, while comparable in size, exhibit marked diversity in their source properties and directivity, with clear spatial patterns, depth‐dependent trends, and a preference for unilateral directivity. These coherent spatial variations source properties suggest that regional differences in tectonic setting, hypocentral depth or fault zone characteristics may drive variability in rupture processes, with important implications for our understanding of earthquake physics and its relation to hazard. 
    more » « less
  4. Abstract Earthquake early warning (EEW) systems aim to forecast the shaking intensity rapidly after an earthquake occurs and send warnings to affected areas before the onset of strong shaking. The system relies on rapid and accurate estimation of earthquake source parameters. However, it is known that source estimation for large ruptures in real‐time is challenging, and it often leads to magnitude underestimation. In a previous study, we showed that machine learning, HR‐GNSS, and realistic rupture synthetics can be used to reliably predict earthquake magnitude. This model, called Machine‐Learning Assessed Rapid Geodetic Earthquake model (M‐LARGE), can rapidly forecast large earthquake magnitudes with an accuracy of 99%. Here, we expand M‐LARGE to predict centroid location and fault size, enabling the construction of the fault rupture extent for forecasting shaking intensity using existing ground motion models. We test our model in the Chilean Subduction Zone with thousands of simulated and five real large earthquakes. The result achieves an average warning time of 40.5 s for shaking intensity MMI4+, surpassing the 34 s obtained by a similar GNSS EEW model. Our approach addresses a critical gap in existing EEW systems for large earthquakes by demonstrating real‐time fault tracking feasibility without saturation issues. This capability leads to timely and accurate ground motion forecasts and can support other methods, enhancing the overall effectiveness of EEW systems. Additionally, the ability to predict source parameters for real Chilean earthquakes implies that synthetic data, governed by our understanding of earthquake scaling, is consistent with the actual rupture processes. 
    more » « less
  5. ABSTRACT Directivity, or the focusing of energy along the direction of an earthquake rupture, is a common property of earthquakes of all sizes and can cause increased hazard due to azimuthally dependent ground-motion amplification. For small earthquakes, the effects of directivity are generally less pronounced due to reduced rupture size, yet the directivity in small events can bias source property estimates and provide important insights into general regional faulting patterns. However, due to observational limitations, directivity is usually only measured and modeled for large events. As such, many studies of small earthquakes either ignore directivity altogether or assume a constant rupture direction for all events in a cluster. In our study, we apply a refined directivity fitting method constrained with two separate methods of source deconvolution to the dataset of aftershocks of the 2019 Ridgecrest earthquakes, which contain a large number of well-recorded small-to-mid sized earthquakes occurring in close proximity to each other. The revealed directivity of 100+ small (M 2.4–5) earthquakes is highly heterogeneous and primarily oblique to and away from the main fault strike, suggesting a complex postseismic stress redistribution. In addition, the energy focusing effect of directivity appears to bias the selection of high-quality data from stations in the direction of rupture, leading to average stress-drop increases of 50% if directivity is not accounted for. 
    more » « less