One of most universal statistical properties of earthquakes is the tendency to cluster in space and time. Yet while clustering is pervasive, individual earthquake sequences can vary markedly in duration, spatial extent, and time evolution. In July 2014, a prolific earthquake sequence initiated within the Sheldon Wildlife Refuge in northwest Nevada, USA. The sequence produced 26 M4 earthquakes and several hundred M3s, with no clear mainshock or obvious driving force. Here we combine a suite of seismological analysis techniques to better characterize this unusual earthquake sequence. High-precision relocations reveal a clear, east dipping normal fault as the dominant structure that intersects with a secondary, subvertical cross fault. Seismicity occurs in burst of activity along these two structures before eventually transitioning to shallower structures to the east. Inversion of hundreds of moment tensors constrain the overall normal faulting stress regime. Source spectral analysis suggests that the stress drops and rupture properties of these events are typical for tectonic earthquakes in the western US. While regional station coverage is sparse in this remote study region, the timely installation of a temporary seismometer allows us to detect nearly 70,000 earthquakes over a 40-month time period when the seismic activity is highest. Such immense productivity is difficult to reconcile with current understanding of crustal deformation in the region and may be facilitated by local hydrothermal processes and earthquake triggering at the transitional intersection of subparallel fault systems.
more »
« less
Coherent Spatial Variations in the Productivity of Earthquake Sequences in California and Nevada
Abstract Earthquakes are clustered in space and time, with individual sequences composed of events linked by stress transfer and triggering mechanisms. On a global scale, variations in the productivity of earthquake sequences—a normalized measure of the number of triggered events—have been observed and associated with regional variations in tectonic setting. Here, we focus on resolving systematic variations in the productivity of crustal earthquake sequences in California and Nevada—the two most seismically active states in the western United States. We apply a well-tested nearest-neighbor algorithm to automatically extract earthquake sequence statistics from a unified 40 yr compilation of regional earthquake catalogs that is complete to M ∼ 2.5. We then compare earthquake sequence productivity to geophysical parameters that may influence earthquake processes, including heat flow, temperature at seismogenic depth, complexity of quaternary faulting, geodetic strain rates, depth to crystalline basement, and faulting style. We observe coherent spatial variations in sequence productivity, with higher values in the Walker Lane of eastern California and Nevada than along the San Andreas fault system in western California. The results illuminate significant correlations between productivity and heat flow, temperature, and faulting that contribute to the understanding and ability to forecast crustal earthquake sequences in the area.
more »
« less
- PAR ID:
- 10505735
- Publisher / Repository:
- Seismological Society of America (Open Access)
- Date Published:
- Journal Name:
- The Seismic Record
- Volume:
- 3
- Issue:
- 4
- ISSN:
- 2694-4006
- Page Range / eLocation ID:
- 322 to 331
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The 2010MW7.2 El Mayor‐Cucapah earthquake ruptured a zone of ~120 km in length in northern Baja California. The geographic distribution of this earthquake sequence was well constrained by waveform relocation. The depth distribution, however, was poorly determined as it is near the edge of, or outside, the Southern California Seismic Network. Here we use two complementary methods to constrain the focal depths of moderate‐sized events (M≥ 4.0) in this sequence. We first determine the absolute earthquake depth by modeling the regional depth phases at high frequencies (~1 Hz). We mainly focus onPnand its depth phasespPnandsPn, which arrive early at regional distance and are less contaminated by crustal multiples. To facilitate depth phase identification and to improve signal‐to‐noise ratio, we take advantage of the dense Southern California Seismic Network and use array analysis to align and stackPnwaveforms. For events without clear depth phases, we further determine their relative depths with respect to those with known depths using differential travel times of thePn, directP, and directSphases recorded for event pairs. Focal depths of 93 out of 122M≥ 4.0 events are tightly constrained with absolute uncertainty of about 1 km. Aftershocks are clustered in the depth range of 3–10 km, suggesting a relatively shallow seismogenic zone, consistent with high surface heat flow in this region. Most aftershocks are located outside or near the lower terminus of coseismic high‐slip patches of the main shock, which may be governed by residual strains, local stress concentration, or postseismic slip.more » « less
-
The state of Nevada is home to one of the most seismically active regions in the world, with crustal deformation associated with the Walker Lane transitioning into Basin and Range tectonics as one traverses from west to east across the state. Despite hosting numerous prominent earthquake sequences over the past century and beyond, at present, there exists no unified research-quality earthquake catalog for the state and its surrounding region. Here, we present a newly compiled, high-precision catalog of more than 180,000 earthquakes occurring around Nevada from 2008 to 2023. The data processing workflow to create this catalog includes an absolute location step that accounts for topography and 3D variations in subsurface wavespeed, and a relative relocation step that refines event positions using differential times measured from waveform cross-correlation. We also provide an update to the local magnitude scale that better accounts for the observed distance attenuation of waveform amplitudes as well as local site effects. We describe some fundamental insights that can be derived from the new catalog, including regional variations in event depth distributions and sequence clustering statistics, and publish the catalog to the wider community to facilitate future research efforts.more » « less
-
Abstract The role of upper‐plate faulting in the seismic cycle of large megathrust earthquakes remains poorly understood. We use quasi‐dynamic numerical simulations of seismic cycles to analyze the interaction between crustal faulting and the foreshock sequence of the 2014 Iquique (Mw 8.2) earthquake in Northern Chile. Multi‐cycle models incorporating upper‐plate faulting align better with coseismic displacements, replicating events akin to the Iquique earthquake. Upper‐plate faulting significantly influences foreshock seismicity and deformation patterns. By calibrating the average hydraulic state—varying the effective normal stress—along the megathrust with pre‐earthquake seismicity, we find that lower pore pressure ratios result in more seismicity before the mainshock. This implies that the hydraulic state of the megathrust is critical for foreshock activity. This comprehensive modeling approach underscores the importance of the mechanical interplay between the megathrust and upper‐plate faults in precursory sequences of large subduction zone earthquakes.more » « less
-
Abstract Foreshocks are the only currently widely identified precursory seismic behavior, yet their utility and even identifiability are problematic, in part because of extreme variation in behavior. Here, we establish some global trends that help identify the expected frequency of foreshocks as well the type of earthquake most prone to foreshocks. We establish these tendencies using the global earthquake catalog of the U.S. Geological Survey National Earthquake Information Center with a completeness level of magnitude 5 and mainshocks with Mw≥7.0. Foreshocks are identified using three clustering algorithms to address the challenge of distinguishing foreshocks from background activity. The methods give a range of 15%–43% of large mainshocks having at least one foreshock but a narrower range of 13%–26% having at least one foreshock with magnitude within two units of the mainshock magnitude. These observed global foreshock rates are similar to regional values for a completeness level of magnitude 3 using the same detection conditions. The foreshock sequences have distinctive characteristics with the global composite population b-values being lower for foreshocks than for aftershocks, an attribute that is also manifested in synthetic catalogs computed by epidemic-type aftershock sequences, which intrinsically involves only cascading processes. Focal mechanism similarity of foreshocks relative to mainshocks is more pronounced than for aftershocks. Despite these distinguishing characteristics of foreshock sequences, the conditions that promote high foreshock productivity are similar to those that promote high aftershock productivity. For instance, a modestly higher percentage of interplate mainshocks have foreshocks than intraplate mainshocks, and reverse faulting events slightly more commonly have foreshocks than normal or strike-slip-faulting mainshocks. The western circum-Pacific is prone to having slightly more foreshock activity than the eastern circum-Pacific.more » « less
An official website of the United States government

