Geology topics

Faults & earthquakes

Explore USGS earthquake and fault publications. Search fault names, locations and seismic hazards; follow the original report for maps, methods and publication context.

Follow this topic with RSS

At least 811 records · Page 45Linked to original sources

USGS2018

Aeromagnetic data reveal potentially seismogenic basement faults in the induced seismicity setting of Oklahoma

New aeromagnetic survey data collected over north central Oklahoma image possible seismogenic faults in the crystalline basement. Linear earthquake sequences associated with induced seismicity suggest the reactivation of ancient basement faults, but few of these sequences are aligned with mapped faults. The new data show many earthquake sequences aligned with linear magnetic gradients or offsets between anomalies, while mapped faults, which mainly describe sedimentary cover, show limited correspondence with either. This strongly suggests significant structural differences between the crystalline basement and sedimentary cover. Furthermore, while the earthquakes are occurring on reactivated ancient faults, most of these faults have likely been inactive for millions of years. The magnetic data exhibit many gradient lineaments that are optimally oriented for fault slip, and the earthquake data suggest additional optimally oriented faults. Together these data suggest the presence of potentially numerous seismogenic faults throughout the region, which may contribute to high levels of induced seismicity.

Oklahoma↗

USGS1996

Type of faulting and orientation of stress and strain as a function of space and time in Kilauea's south flank, Hawaii

Earthquake focal mechanisms of events occurring between 1972 and 1992 in the south flank of Kilauea volcano, Hawaii, are used to infer the state of stress and strain as a function of time and space. We have determined 870 fault plane solutions from P wave first motion polarities for events with magnitudes M L ≥ 2.5 and depth ranging between 6 and 12 km. Faulting is characterized by a mixture of decollement, reverse, and normal faults. Most large earthquakes with magnitude M < 7 slip on reverse faults striking NE at 40° and dipping SE between 60° and 70°. In Hawaii, the earthquakes with M > 7 rupture the decollement plane, since it is the only surface large enough to generate magnitude 7 or larger earthquakes. The percentage of reverse faulting events is high compared to the decollement and normal faulting mechanisms for the period 1972–1983. The percentage of decollement type focal mechanisms becomes dominant after 1983. This pattern of faulting activity suggests that pressure was building up within Kilauea's rift zone prior to the 1983 Puu'Oo eruption. Overall, a single stress orientation with the maximum compressive stress oriented SE perpendicular to the rift and dipping at 45° is compatible with the coeval existence of decollement, reverse, and normal faults. However, in a crustal volume east of longitude 155°10′W, we find a change of the orientation of σ 1 from nearly horizontal to plunging 45° SE occurring in 1979. This stress rotation suggests magma movements within the aseismic part of Kilauea's east rift zone. The strain and stress orientations are coaxial in the south flank except within the volume where the stress rotation is observed. We observe a change in the relationship between stress and strain directions caused either by the shifting of seismic activity from reverse faults to decollements, while stress stays constant, or by a rotation of stress, while strain remains constant. Assuming that the model of a noncohesive Coulomb wedge is appropriate for Kilauea's south flank, we find that high pore pressures are prevalent along the decollement and within the wedge for a coefficient of friction equal to 0.85.

Hawaii↗

USGS1996

Metastable mantle phase transformations and deep earthquakes in subducting oceanic lithosphere

Earth's deepest earthquakes occur as a population in subducting or previously subducted lithosphere at depths ranging from about 325 to 690 km. This depth interval closely brackets the mantle transition zone, characterized by rapid seismic velocity increases resulting from the transformation of upper mantle minerals to higher-pressure phases. Deep earthquakes thus provide the primary direct evidence for subduction of the lithosphere to these depths and allow us to investigate the deep thermal, thermodynamic, and mechanical ferment inside slabs. Numerical simulations of reaction rates show that the olivine → spinel transformation should be kinetically hindered in old, cold slabs descending into the transition zone. Thus wedge-shaped zones of metastable peridotite probably persist to depths of more than 600 km. Laboratory deformation experiments on some metastable minerals display a shear instability called transformational faulting. This instability involves sudden failure by localized superplasticity in thin shear zones where the metastable host mineral transforms to a denser, finer-grained phase. Hence in cold slabs, such faulting is expected for the polymorphic reactions in which olivine transforms to the spinel structure and clinoenstatite transforms to ilmenite. It is thus natural to hypothesize that deep earthquakes result from transformational faulting in metastable peridotite wedges within cold slabs. This consideration of the mineralogical states of slabs augments the traditional largely thermal view of slab processes and explains some previously enigmatic slab features. It explains why deep seismicity occurs only in the approximate depth range of the mantle transition zone, where minerals in downgoing slabs should transform to spinel and ilmenite structures. The onset of deep shocks at about 325 km is consistent with the onset of metastability near the equilibrium phase boundary in the slab. Even if a slab penetrates into the lower mantle, earthquakes should cease at depths near 700 km, because the seismogenic phase transformations in the slab are completed or can no longer occur. Substantial metastability is expected only in old, cold slabs, consistent with the observed restriction of deep earthquakes to those settings. Earthquakes should be restricted to the cold cores of slabs, as in any model in which the seismicity is temperature controlled, via the distribution of metastability. However, the geometries of recent large deep earthquakes pose a challenge for any such models. Transformational faulting may give insight into why deep shocks lack appreciable aftershocks and why their source characteristics, including focal mechanisms indicating localized shear failure rather than implosive deformation, are so similar to those of shallow earthquakes. Finally, metastable phase changes in slabs would produce an internal source of stress in addition to those due to the weight of the sinking slab. Such internal stresses may explain the occurrence of earthquakes in portions of lithosphere which have foundered to the bottom of the transition zone and/or are detached from subducting slabs. Metastability in downgoing slabs could have considerable geodynamic significance. Metastable wedges would reduce the negative buoyancy of slabs, decrease the driving force for subduction, and influence the state of stress in slabs. Heat released by metastable phase changes would raise temperatures within slabs and facilitate the transformation of spinel to the lower mantle mineral assemblage, causing slabs to equilibrate more rapidly with the ambient mantle and thus contribute to the cessation of deep seismicity. Because wedge formation should occur only for fast subducting slabs, it may act as a “parachute” and contribute to regulating plate speeds. Wedge formation would also have consequences for mantle evolution because the density of a slab stagnated near the bottom of the transition zone would increase as it heats up and the wedge transforms to denser spinel, favoring the subsequent sinking of the slab into the lower mantle.

Reviews of Geophysics↗

USGS2020

Relocated aftershocks and background seismicity in eastern Indonesia shed light on the 2018 Lombok and Palu earthquake sequences

High seismicity rates in eastern Indonesia occur due to the complex interaction of several tectonic plates which resulted in two deadly, destructive earthquake sequences that occurred in Lombok Island and the city of Palu, Sulawesi in 2018. The first sequence began in July with an M w 6.4 event near Lombok, culminating in an M w 7.0 event 8 d later. This was then followed by a nearby M w 6.9 event 12 d later. Approximately 1000 km to the northeast, a separate sequence began several weeks later near Palu where an M w 7.5 event occurred that triggered a tsunami. In this study, we present hypocentre relocations for both earthquake sequences as well as all other regional earthquakes in eastern Indonesia. The relocations were performed using a teleseismic double-difference relocation method and arrival times for P and S waves from stations at local, regional, and teleseismic distances. The catalogue and phase data were taken from the Agency for Meteorology, Climatology and Geophysics (BMKG) of Indonesia and the International Seismological Centre (ISC) for the period of April 2009 through November 2018. The relocated catalogue provides an improved view of seismicity in eastern Indonesia over the study period, sharpening locations and interpretations of seismogenic features throughout the region. In the Lombok area, the relocated earthquakes clearly show a backarc thrust to the north of the Sunda-Banda Arc transition zone. The relocated aftershocks show that the destructive M w 7.0 and M w 6.9 earthquakes of the Lombok sequence ruptured two different regions: The M w 7.0 earthquake propagated westward, whereas the M w 6.9 earthquake propagated eastward. The entire sequence of Lombok earthquakes was most likely started by the M w 6.4 event as the initial event or foreshock, which then triggered backarc thrusts on both sides. Several weeks later and far to the northeast, the M w 7.5 Palu earthquake occurred along the Palu-Koro Fault, filling a seismic gap that had not ruptured in an M w 6.0 event or larger since at least 1900. The distribution of aftershocks indicates that the northern part of the Palu-Koro Fault has lower relative seismicity rates than the southern part at shallow depths, and that off fault aftershocks are mostly located to the east of the Palu-Koro Fault.

Lombok, Palu↗

USGS1998

Near-surface structural model for deformation associated with the February 7, 1812, New Madrid, Missouri, earthquake

The February 7, 1812, New Madrid, Missouri, earthquake (M [moment magnitude] 8) was the third and final large-magnitude event to rock the northern Mississippi Embayment during the winter of 1811–1812. Although ground shaking was so strong that it rang church bells, stopped clocks, buckled pavement, and rocked buildings up and down the eastern seaboard, little coseismic surface deformation exists today in the New Madrid area. The fault(s) that ruptured during this event have remained enigmatic. We have integrated geomorphic data documenting differential surficial deformation (supplemented by historical accounts of surficial deformation and earthquake-induced Mississippi River waterfalls and rapids) with the interpretation of existing and recently acquired seismic reflection data, to develop a tectonic model of the near-surface structures in the New Madrid, Missouri, area. This model consists of two primary components: a north-northwest–trending thrust fault and a series of northeast-trending, strike-slip, tear faults. We conclude that the Reelfoot fault is a thrust fault that is at least 30 km long. We also infer that tear faults in the near surface partitioned the hanging wall into subparallel blocks that have undergone differential displacement during episodes of faulting. The northeast-trending tear faults bound an area documented to have been uplifted at least 0.5 m during the February 7, 1812, earthquake. These faults also appear to bound changes in the surface density of epicenters that are within the modern seismicity, which is occurring in the stepover zone of the left-stepping right-lateral strike-slip fault system of the modern New Madrid seismic zone.

Missouri↗

USGS2015

The 2011 Mineral, Virginia, earthquake and its significance for seismic hazards in eastern North America: overview and synthesis

The 23 August 2011 M w (moment magnitude) 5.7 ± 0.1, Mineral, Virginia, earthquake was the largest and most damaging in the central and eastern United States since the 1886 M w 6.8–7.0, Charleston, South Carolina, earthquake. Seismic data indicate that the earthquake rupture occurred on a southeast-dipping reverse fault and consisted of three subevents that progressed northeastward and updip. U.S. Geological Survey (USGS) "Did You Feel It?" intensity reports from across the eastern United States and southeastern Canada, rockfalls triggered at distances to 245 km, and regional groundwater-level changes are all consistent with efficient propagation of high-frequency seismic waves (&sim;1 Hz and higher) in eastern North America due to low attenuation. Reported damage included cracked or shifted foundations and broken walls or chimneys, notably in unreinforced masonry, and indicated intensities up to VIII in the epicentral area based on USGS "Did You Feel It?" reports. The earthquake triggered the first automatic shutdown of a U.S. nuclear power plant, located &sim;23 km northeast of the main shock epicenter. Although shaking exceeded the plant's design basis earthquake, the actual damage to safety-related structures, systems, and components was superficial. Damage to relatively tall masonry structures 130 km to the northeast in Washington, D.C., was consistent with source directivity, soft-soil ground-motion amplification, and anisotropic wave propagation with lower attenuation parallel to the northeast-trending Appalachian tectonic fabric. The earthquake and aftershocks occurred in crystalline rocks within Paleozoic thrust sheets of the Chopawamsic terrane. The main shock and majority of aftershocks delineated the newly named Quail fault zone in the subsurface, and shallow aftershocks defined outlying faults. The earthquake induced minor liquefaction sand boils, but notably there was no evidence of a surface fault rupture. Recurrence intervals, and evidence for larger earthquakes in the Quaternary in this area, remain important unknowns. This event, along with similar events during historical time, is a reminder that earthquakes of similar or larger magnitude pose a real hazard in eastern North America.

Virginia↗

USGS2022

The 8 April 1860 Jour de Pâques earthquake sequence in southern Haiti

The grave threat posed by the Enriquillo‐Plantain Garden fault zone (EPGFZ) and other fault systems on the Tiburon Peninsula in southern Haiti was highlighted by the catastrophic M 7.0 Léogâne earthquake on 12 January 2010 and again by the deadly M 7.2 Nippes earthquakes on 14 August 2021. Early Interferometric Synthetic Aperture Radar observations suggest the 2021 earthquake broke structures associated with this fault system farther west of the 2010 event, but the rupture zones of both events are separated by a ∼50 km gap. This sequence provided the impetus to reconsider a nineteenth century earthquake that may have occurred within this gap. Though previous studies identified a single moderately large event on 8 April 1860, original sources describe a complex and distributed seismic sequence to the west of Port‐au‐Prince. These provide evidence for an initial event to the west of Les Cayes, on the southern coast of the Tiburon Peninsula. This was followed on the morning of 8 April 1860 by a damaging earthquake near l’Anse‐à‐Veau along the northern coast of the peninsula, which was succeeded 14 hr later by a larger mainshock to the east. Although locations cannot be determined precisely from extant macroseismic data, our preferred scenario includes an intensity magnitude ( M I ) ∼ 6.2 "> ( M I ) ∼ 6.2 (MI)∼6.2 foreshock, followed by an M I ∼ 6.8 "> M I ∼ 6.8 MI∼6.8 mainshock wholly or partially on the EPGFZ to the east of the foreshock. Our results suggest that, although earthquakes tend to trigger subsequent earthquakes on immediately neighboring fault segments, sequences can include gaps due to strain release by poorly characterized earthquakes in the historic past. We also document 83 aftershocks to the 1860 earthquake sequence, including a previously unknown damaging shock, estimated M I ∼ 6.6 "> M I ∼ 6.6 MI∼6.6 ⁠ , on 10 April 1860 near Jérémie, farther west of l’Anse‐à‐Veau.

Tiburon Peninsula↗

USGS2016

Preliminary assessment of a previously unknown fault zone beneath the Daytona Beach sand blow cluster near Marianna, Arkansas

We collected new high‐resolution P ‐wave seismic‐reflection data to explore for possible faults beneath a roughly linear cluster of early to mid‐Holocene earthquake‐induced sand blows to the south of Marianna, Arkansas. The Daytona Beach sand blow deposits are located in east‐central Arkansas about 75 km southwest of Memphis, Tennessee, and about 80 km south of the southwestern end of the New Madrid seismic zone (NMSZ). Previous studies of these sand blows indicate that they were produced between 10,500 and 5350 yr B.P. (before A.D. 1950). The sand blows are large and similar in size to those in the heart of the NMSZ produced by the 1811–1812 earthquakes. The seismic‐reflection profiles reveal a previously unknown zone of near‐vertical faults imaged in the 100–1100‐m depth range that are approximately coincident with a cluster of earthquake‐induced sand blows and a near‐linear surface lineament composed of air photo tonal anomalies. These interpreted faults are expressed as vertical discontinuities with the largest displacement fault showing about 40 m of west‐side‐up displacement at the top of the Paleozoic section at about 1100 m depth. There are about 20 m of folding on reflections within the Eocene strata at 400 m depth. Increasing fault displacement with depth suggests long‐term recurrent faulting. The imaged faults within the vicinity of the numerous sand blow features could be a causative earthquake source, although it does not rule out the possibility of other seismic sources nearby. These newly located faults add to a growing list of potentially active Pleistocene–Holocene faults discovered over the last two decades that are within the Mississippi embayment region but outside of the historical NMSZ.

Arkansas↗

USGS1994

Magnetic field observations in the near-field the 28 June 1992 Mw 7.3 Landers, California, earthquake

Recent reports suggest that large magnetic field changes occur prior to, and during, large earthquakes. Two continuously operating proton magnetometers, LSBM and OCHM , at distances of 17.3 and 24.2 km, respectively, from the epicenter of the 28 June 1992 M w 7.3 Landers earthquake, recorded data through the earthquake and its aftershocks. These two stations are part of a differentially connected array of proton magnetometers that has been operated along the San Andreas fault since 1976. The instruments have a sensitivity of 0.25 nT or better and transmit data every 10 min through the GOES satellite to the USGS headquarters in Menlo Park, California. Seismomagnetic offsets of −1.2 ± 0.6 and −0.7 ± 0.7 nT were observed at these sites. In comparison, offsets of −0.3 ± 0.2 and −1.3 ± 0.2 nT were observed during the 8 July 1986 M L 5.9 North Palm Springs earthquake, which occurred directly beneath the OCHM magnetometer site. The observations are generally consistent with seismomagnetic models of the earthquake, in which fault geometry and slip have the same from as that determined by either inversion of the seismic data or inversion of geodetically determined ground displacements produced by the earthquake. In these models, right-lateral rupture occurs on connected fault segments in a homogeneous medium with average magnetization of 2 A/m. The fault-slip distribution has roughly the same form as the observed surface rupture, and the total moment release is 1.1 × 10 20 Nm. There is no indication of diffusion-like character to the magnetic field offsets that might indicate these effects result from fluid flow phenomena. It thus seems unlikely that these earthquake-generated offsets and those produced by the North Palm Springs earthquake were generated by electrokinetic effects. Also, there are no indications of enhanced low-frequency magnetic noise before the earthquake at frequencies below 0.001 Hz.

California↗

USGS2013

Fine-scale delineation of the location of and relative ground shaking within the San Andreas Fault zone at San Andreas Lake, San Mateo County, California

The San Francisco Public Utilities Commission is seismically retrofitting the water delivery system at San Andreas Lake, San Mateo County, California, where the reservoir intake system crosses the San Andreas Fault (SAF). The near-surface fault location and geometry are important considerations in the retrofit effort. Because the SAF trends through highly distorted Franciscan mélange and beneath much of the reservoir, the exact trace of the 1906 surface rupture is difficult to determine from surface mapping at San Andreas Lake. Based on surface mapping, it also is unclear if there are additional fault splays that extend northeast or southwest of the main surface rupture. To better understand the fault structure at San Andreas Lake, the U.S. Geological Survey acquired a series of seismic imaging profiles across the SAF at San Andreas Lake in 2008, 2009, and 2011, when the lake level was near historical lows and the surface traces of the SAF were exposed for the first time in decades. We used multiple seismic methods to locate the main 1906 rupture zone and fault splays within about 100 meters northeast of the main rupture zone. Our seismic observations are internally consistent, and our seismic indicators of faulting generally correlate with fault locations inferred from surface mapping. We also tested the accuracy of our seismic methods by comparing our seismically located faults with surface ruptures mapped by Schussler (1906) immediately after the April 18, 1906 San Francisco earthquake of approximate magnitude 7.9; our seismically determined fault locations were highly accurate. Near the reservoir intake facility at San Andreas Lake, our seismic data indicate the main 1906 surface rupture zone consists of at least three near-surface fault traces. Movement on multiple fault traces can have appreciable engineering significance because, unlike movement on a single strike-slip fault trace, differential movement on multiple fault traces may exert compressive and extensional stresses on built structures within the fault zone. Such differential movement and resulting distortion of built structures appear to have occurred between fault traces at the gatewell near the southern end of San Andreas Lake during the 1906 San Francisco earthquake (Schussler, 1906). In addition to the three fault traces within the main 1906 surface rupture zone, our data indicate at least one additional fault trace (or zone) about 80 meters northeast of the main 1906 surface rupture zone. Because ground shaking also can damage structures, we used fault-zone guided waves to investigate ground shaking within the fault zones relative to ground shaking outside the fault zones. Peak ground velocity (PGV) measurements from our guided-wave study indicate that ground shaking is greater at each of the surface fault traces, varying with the frequency of the seismic data and the wave type (P versus S). S-wave PGV increases by as much as 5–6 times at the fault traces relative to areas outside the fault zone, and P-wave PGV increases by as much as 3–10 times. Assuming shaking increases linearly with increasing earthquake magnitude, these data suggest strong shaking may pose a significant hazard to built structures that extend across the fault traces. Similarly complex fault structures likely underlie other strike-slip faults (such as the Hayward, Calaveras, and Silver Creek Faults) that intersect structures of the water delivery system, and these fault structures similarly should be investigated.

California↗

USGS1979

Trees as indicators of past movements on the San Andreas Fault

Trees are sources of information about fault movements that have occurred before the earliest historical reports. This kind of evidence can be used to improve estimates of when earthquakes will recur on faults known to be seismically active and to identify active faults that have no record of movement during recent history. The approach is not new. Robert Page of the U.S Geological Survey described the effects of the 1958 earthquake on trees along the Fairweather fault in Alaska. He showed that tree rings methods could have been used to identify and to closely date this event. We undertook a similar study of the northern part of the San Andreas fault, in part because there are no historic records prior to 1906 along this segment. Earthquakes and surface rupture along faults affect trees in several different ways. Direct effects include fracturing, twisting, and tilting of trees that grown on the surface of the break. In a much wider zone along the fault, trees may be felled or topped as a result of ground motion. Among the indirect effects are tilting, felling, or burial of trees in earthquake-triggered landslides. Long-term effects may include changes in growth rate due to local hydrologic and topographic changes as well as to biological effects such as the death of neighboring trees. Under favoralbe circumstances these can be dated by tree ring methods.

California↗

USGS1998

Multiscale seismic imaging of active fault zones for hazard assessment: A case study of the Santa Monica fault zone, Los Angeles, California

High‐resolution seismic reflection profiles at two different scales were acquired across the transpressional Santa Monica Fault of north Los Angeles as part of an integrated hazard assessment of the fault. The seismic data confirm the location of the fault and related shallow faulting seen in a trench to deeper structures known from regional studies. The trench shows a series of near‐vertical strike‐slip faults beneath a topographic scarp inferred to be caused by thrusting on the Santa Monica fault. Analysis of the disruption of soil horizons in the trench indicates multiple earthquakes have occurred on these strike‐slip faults within the past 50 000 years, with the latest being 1000 to 3000 years ago. A 3.8-km-long, high‐resolution seismic reflection profile shows reflector truncations that constrain the shallow portion of the Santa Monica Fault (upper 300 m) to dip northward between 30° and 55°, most likely 30° to 35°, in contrast to the 60° to 70° dip interpreted for the deeper portion of the fault. Prominent, nearly continuous reflectors on the profile are interpreted to be the erosional unconformity between the 1.2 Ma and older Pico Formation and the base of alluvial fan deposits. The unconformity lies at depths of 30–60 m north of the fault and 110–130 m south of the fault, with about 100 m of vertical displacement (180 m of dip‐slip motion on a 30°–35° dipping fault) across the fault since deposition of the upper Pico Formation. The continuity of the uncomformity on the seismic profile constrains the fault to lie in a relatively narrow (50 m) zone, and to project to the surface beneath Ohio Avenue immediately south of the trench. A very high‐resolution seismic profile adjacent to the trench images reflectors in the 15 to 60 m depth range that are arched slightly by folding just north of the fault. A disrupted zone on the profile beneath the south end of the trench is interpreted as being caused by the deeper portions of the trenched strike‐slip faults where they merge with the thrust fault.

Geophysics↗

USGS1991

Late Quaternary faulting along the Death Valley-Furnace Creek fault system, California and Nevada

The Death Valley-Furnace Creek fault system, in California and Nevada, has a variety of impressive late Quaternary neotectonic features that record a long history of recurrent earthquake-induced faulting. Although no neotectonic features of unequivocal historical age are known, paleoseismic features from multiple late Quaternary events of surface faulting are well developed throughout the length of the system. Comparison of scarp heights to amount of horizontal offset of stream channels and the relationships of both scarps and channels to the ages of different geomorphic surfaces demonstrate that Quaternary faulting along the northwest-trending Furnace Creek fault zone is predominantly right lateral, whereas that along the north-trending Death Valley fault zone is predominantly normal. These observations are compatible with tectonic models of Death Valley as a northwest-trending pull-apart basin. The largest late Quaternary scarps along the Furnace Creek fault zone, with vertical separation of late Pleistocene surfaces of as much as 64 m (meters), are in Fish Lake Valley. Despite the predominance of normal faulting along the Death Valley fault zone, vertical offset of late Pleistocene surfaces along the Death Valley fault zone apparently does not exceed about 15 m. Evidence for four to six separate late Holocene faulting events along the Furnace Creek fault zone and three or more late Holocene events along the Death Valley fault zone are indicated by rupturing of Q 1B (about 200-2,000 years old) geomorphic surfaces. Probably the youngest neotectonic feature observed along the Death Valley-Furnace Creek fault system, possibly historic in age, is vegetation lineaments in southernmost Fish Lake Valley. Near-historic faulting in Death Valley, within several kilometers south of Furnace Creek Ranch, is represented by (1) a 2,000-year-old lake shoreline that is cut by sinuous scarps, and (2) a system of young scarps with free-faceted faces (representing several faulting events) that cuts Q 1B surfaces.

California, Nevada↗

USGS2017

The finite, kinematic rupture properties of great-sized earthquakes since 1990

Here, I present a database of >160 finite fault models for all earthquakes of M 7.5 and above since 1990, created using a consistent modeling approach. The use of a common approach facilitates easier comparisons between models, and reduces uncertainties that arise when comparing models generated by different authors, data sets and modeling techniques. I use this database to verify published scaling relationships, and for the first time show a clear and intriguing relationship between maximum potency (the product of slip and area) and average potency for a given earthquake. This relationship implies that earthquakes do not reach the potential size given by the tectonic load of a fault (sometimes called “moment deficit,” calculated via a plate rate over time since the last earthquake, multiplied by geodetic fault coupling). Instead, average potency (or slip) scales with but is less than maximum potency (dictated by tectonic loading). Importantly, this relationship facilitates a more accurate assessment of maximum earthquake size for a given fault segment, and thus has implications for long-term hazard assessments. The relationship also suggests earthquake cycles may not completely reset after a large earthquake, and thus repeat rates of such events may appear shorter than is expected from tectonic loading. This in turn may help explain the phenomenon of “earthquake super-cycles” observed in some global subduction zones.

Earth and Planetary Science Letters↗

USGS2020

Plate boundary localization, slip-rates and rupture segmentation of the Queen Charlotte Fault based on submarine tectonic geomorphology

Linking fault behavior over many earthquake cycles to individual earthquake behavior is a primary goal in tectonic geomorphology, particularly across an entire plate boundary. Here, we examine the 1150-km-long, right-lateral Queen Charlotte-Fairweather fault system using comprehensive multibeam bathymetry data acquired along the Queen Charlotte Fault (QCF) offshore southeastern Alaska and western British Columbia. Fine-scale analysis of tectonic geomorphology allowed us to identify and reconstruct 184 strike-slip piercing points over a 630 km stretch of the QCF. Age constraints from glacial recession and offshore sedimentation patterns yield a consistent slip-rate of ∼50–57 mm/yr since ∼17–12 ka, the fastest rate for a continent-ocean strike-slip fault on Earth. These slip-rates equal or exceed estimates of Pacific-North America (PA-NA) relative motion from global plate reconstructions, indicating that PA-NA motion is highly localized. The QCF cuts the seafloor along a narrow and unusually straight trace for its entire length and multiple fault traces are observed only at local step-overs. The geometry and behavior of the QCF over many earthquake cycles is simple and typical of mature faults with relatively homogeneous stress fields. Since the QCF is the primary PA-NA plate boundary, we used the trace of the QCF to define the small circle path for relative plate motion and computed the associated Euler pole. Predicted along-strike obliquity variations based on the new pole agree with observed tectonic geomorphology and suggest that previous global plate reconstructions overestimated the degree of oblique convergence along the QCF. We also find that subtle, long-wavelength (75–150 km) bends and discrete step-overs appear to define the endpoints of M>7 earthquakes, suggesting that obliquity and resultant fault geometry may control rupture segmentation and asperity development. Lastly, the agreement between predicted obliquity and tectonic geomorphology along the entire length of QCF compelled a reevaluation of regional tectonic models. In the north, the eastern Yakatat Terrane appears to be translating northwest with the Pacific plate, and slip transferred from the QCF to the Fairweather Fault results in ∼20 mm/yr of convergence along the southern St. Elias mountains. In the south, we predict a reduced rate of convergence along the QCF west of Haida Gwaii (∼5–6 mm/yr of shortening, on average) relative to previous studies. Our results support a model for transpression and strike-slip partitioning along the edge of a hot and weak Pacific Plate, leading to crustal thickening and growth of the Queen Charlotte Terrace to the west of Haida Gwaii.

Alaska, British Columbia↗

USGS2001

Elastic rebound following the Kocaeli earthquake, Turkey, recorded using synthetic aperture radar interferometry

A basic model incorporating satellite synthetic aperture radar (SAR) interferometry of the fault rupture zone that formed during the Kocaeli earthquake of August 17, 1999, documents the elastic rebound that resulted from the concomitant elastic strain release along the North Anatolian fault. For pure strike-slip faults, the elastic rebound function derived from SAR interferometry is directly invertible from the distribution of elastic strain on the fault at criticality, just before the critical shear stress was exceeded and the fault ruptured. The Kocaeli earthquake, which was accompanied by as much as ∼5 m of surface displacement, distributed strain ∼110 km around the fault prior to faulting, although most of it was concentrated in a narrower and asymmetric 10-km-wide zone on either side of the fault. The use of SAR interferometry to document the distribution of elastic strain at the critical condition for faulting is clearly a valuable tool, both for scientific investigation and for the effective management of earthquake hazard.

Geology↗

USGS2007

The persistence of directivity in small earthquakes

We derive a simple inversion of peak ground acceleration (PGA) or peak ground velocity (PGV) for rupture direction and rupture velocity and then test this inversion on the peak motions obtained from seven 3.5≤ M ≤4.1 earthquakes that occurred in two clusters in November 2002 and February 2003 near San Ramon, California. These clusters were located on two orthogonal strike-slip faults so that the events share the same approximate focal mechanism but not the same fault plane. Three earthquakes exhibit strong directivity, but the other four earthquakes exhibit relatively weak directivity. We use the residual PGAs and PGVs from the other six events to determine station corrections for each earthquake. The inferred rupture directions unambiguously identify the fault plane for the three earthquakes with strong directivity and for three of the four earthquakes with weak directivity. The events with strong directivity have fast rupture velocities (0.63 β ≤ v ≤0.87 β ); the events with weak directivity either rupture more slowly (0.17 β ≤ v ≤0.35 β ) or bilaterally. The simple unilateral inversion cannot distinguish between slow and bilateral ruptures: adding a bilateral rupture component degrades the fit of the rupture directions to the fault planes. By comparing PGAs from the events with strong and weak directivity, we show how an up-dip rupture in small events can distort the attenuation of peak ground motion with distance. When we compare the rupture directions of the earthquakes to the location of aftershocks in the two clusters, we find than almost all the aftershocks of the three earthquakes with strong directivity occur within 70° of the direction of rupture.

Bulletin of the Seismological Society of America↗

USGS2020

The community code verification exercise for simulating sequences of earthquakes and aseismic slip (SEAS)

Numerical simulations of sequences of earthquakes and aseismic slip (SEAS) have made great progress over past decades to address important questions in earthquake physics. However, significant challenges in SEAS modeling remain in resolving multiscale interactions between earthquake nucleation, dynamic rupture, and aseismic slip, and understanding physical factors controlling observables such as seismicity and ground deformation. The increasing complexity of SEAS modeling calls for extensive efforts to verify codes and advance these simulations with rigor, reproducibility, and broadened impact. In 2018, we initiated a community code-verification exercise for SEAS simulations, supported by the Southern California Earthquake Center. Here, we report the findings from our first two benchmark problems (BP1 and BP2), designed to verify different computational methods in solving a mathematically well-defined, basic faulting problem. We consider a 2D antiplane problem, with a 1D planar vertical strike-slip fault obeying rate-and-state friction, embedded in a 2D homogeneous, linear elastic halfspace. Sequences of quasi-dynamic earthquakes with periodic occurrences (BP1) or bimodal sizes (BP2) and their interactions with aseismic slip are simulated. The comparison of results from 11 groups using different numerical methods show excellent agreements in long-term and coseismic fault behavior. In BP1, we found that truncated domain boundaries influence interseismic stressing, earthquake recurrence, and coseismic rupture, and that model agreement is only achieved with sufficiently large domain sizes. In BP2, we found that complexity of fault behavior depends on how well physical length scales related to spontaneous nucleation and rupture propagation are resolved. Poor numerical resolution can result in artificial complexity, impacting simulation results that are of potential interest for characterizing seismic hazard such as earthquake size distributions, moment release, and recurrence times. These results inform the development of more advanced SEAS models, contributing to our further understanding of earthquake system dynamics.

Seismological Research Letters↗