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 1,009 records · Page 56Linked to original sources

USGS2020

Coexisting seismic behavior of transform faults revealed by high-resolution bathymetry

Transform faults are known to have anomalously low rates of seismicity, but no direct observations reveal why this is the case. We use new, autonomous underwater vehicle high-resolution seafloor mapping to image the morphology of and offsets along transform fault segments in the Gulf of California. Fault splays display a varied history of activation and deactivation of individual fault strands over time, not unlike those mapped onshore or imaged within the bathymetry of the Queen Charlotte-Fairweather and the Palos Verdes faults of offshore western Canada and Southern California. A series of six identically offset depositional fans evidence 21–23 meters of slip along the main transform fault, which could not have been produced by a single earthquake. Rather, the lack of smaller-magnitude offsets indicates synchronous deposition and an absence of multiple slope failure-inducing earthquakes, thus providing the first direct evidence that creep and earthquakes occur at different times in the slip history of a given transform fault segment.

Geology↗

USGS2022

Climatic influence on the expression of strike-slip faulting

Earthquakes on strike-slip faults are preserved in the geomorphic record by offset landforms that span a range of displacements, from small offsets created in the most recent earthquake (MRE) to large offsets that record cumulative slip from multiple prior events. An exponential decay in the number of large cumulative offsets has been observed on many faults, and a leading hypothesis is that climate controls the rate of decay. We present offset measurements compiled from 31 studies of strike-slip faults with evidence of multiple paleoearthquakes and corresponding climatic and tectonic information to test this hypothesis. Both the global compilation and numerical landscape evolution modeling reveal that the decay rate in large offsets is negatively correlated with mean annual precipitation. Faults in dry regions with high drainage density more commonly preserve small MRE offsets, and faults in wet regions with lower drainage density more commonly preserve a mix of small MRE and large cumulative offsets. Geomorphology of faults in different climates supports this result and illustrates precipitation’s effect on the development and preservation of offset channels. Our findings imply that current and past climate affect how displacement on strike-slip faults is recorded and interpreted to inform earthquake history.

Geology↗

USGS2003

Gravity and magnetic expression of the San Leandro gabbro with implications for the geometry and evolution of the Hayward Fault zone, northern California

The Hayward Fault, one of the most hazardous faults in northern California, trends north-northwest and extends for about 90 km along the eastern San Francisco Bay region. At numerous locations along its length, distinct and elongate gravity and magnetic anomalies correlate with mapped mafic and ultramafic rocks. The most prominent of these anomalies reflects the 16-km-long San Leandro gabbroic block. Inversion of magnetic and gravity data constrained with physical property measurements is used to define the subsurface extent of the San Leandro gabbro body and to speculate on its origin and relationship to the Hayward Fault Zone. Modeling indicates that the San Leandro gabbro body is about 3 km wide, dips about 75??-80?? northeast, and extends to a depth of at least 6 km. One of the most striking results of the modeling, which was performed independently of seismicity data, is that accurately relocated seismicity is concentrated along the western edge or stratigraphically lower bounding surface of the San Leandro gabbro. The western boundary of the San Leandro gabbro block is the base of an incomplete ophiolite sequence and represented at one time, a low-angle roof thrust related to the tectonic wedging of the Franciscan Complex. After repeated episodes of extension and attenuation, the roof thrust of this tectonic wedge was rotated to near vertical, and in places, the strike-slip Hayward Fault probably reactivated or preferentially followed this pre-existing feature. Because earthquakes concentrate near the edge of the San Leandro gabbro but tend to avoid its interior, we qualitatively explore mechanical models to explain how this massive igneous block may influence the distribution of stress. The microseismicity cluster along the western flank of the San Leandro gabbro leads us to suggest that this stressed volume may be the site of future moderate to large earthquakes. Improved understanding of the three-dimensional geometry and physical properties along the Hayward Fault will provide additional constraints on seismic hazard probability, earthquake modeling, and fault interactions that are applicable to other major strike-slip faults around the world.

California↗

USGS2007

Late quaternary paleoseismology of the southern Steens fault zone, northern Nevada

The 192-km-long Steens fault zone is the most prominent normal fault system in the northern Basin and Range province of western North America. We use trench mapping and radiometric dating to estimate displacements and timing of the last three surface-rupturing earthquakes (E1-E3) on the southern part of the fault south of Denio, Nevada. Coseismic displacements range from 1.1 to 2.2 ?? 0.5 m, and radiometric ages indicate earthquake times of 11.5 ?? 2.0 ka (E3), 6.1 ?? 0.5 ka (E2), and 4.6 ?? 1.0 ka (E1). These data yield recurrence intervals of 5.4 ?? 2.1 k.y. between E3 and E2, 1.5 ?? 1.1 k.y. between E2 and E1, and an elapsed time of 4.6 ?? 1.0 k.y. since E1. The recurrence data yield variable interval slip rates (between 0.2 ?? 0.22 and 1.5 ?? 2.3 mm/yr), but slip rates averaged over the past ???18 k.y. (0.24 ?? 0.06 mm/year) are similar to long-term (8.5-12.5 Ma) slip rates (0.2 ?? 0.1 mm /yr) measured a few kilometers to the north. We infer from the lack of significant topographic relief across the fault in Bog Hot Valley that the fault zone is propagating southward and may now be connected with a fault at the northwestern end of the Pine Forest Range. Displacements documented in the trench and a rupture length of 37 km indicate a history of three latest Quaternary earthquakes with magnitudes of M 6.6-7.1 on the southern part of the Steens fault zone.

Bulletin of the Seismological Society of America↗

USGS2020

Estimating rupture dimensions of three major earthquakes in Sichuan, China, for early warning and rapid loss estimates

Large earthquakes like in Wenchuan in 2008, MW 7.9, Sichuan, China, provide opportunity for earthquake early warning (EEW) as many heavily shaken areas are far (~50 km) from the epicenter and warning time could be long enough (≥ 5 s) to take effective preventative action. On the other hand, earthquakes with magnitudes larger than ~M 6.5 are challenging for EEW since source dimensions need to be defined in order to adequately estimate shaking. The Finite-Fault Rupture Detector (FinDer) is an approach to identify fault rupture extents from real-time strong motion and/or broadband records. In this study, we playback local and regional on-scale strong motion waveforms recorded during the 2008 MW 7.9 Wenchuan, 2013 MW 6.6 Lushan, and 2017 MW 6.5 Jiuzhaigou earthquakes to study the performance of FinDer for the current layout of the China Strong Motion Network. Overall, the FinDer line-source models agree well with the observed spatial distribution of aftershocks and fault models determined from waveform inversion. However, since FinDer models are constructed to characterize seismic ground motions (as needed for EEW) instead of source parameters, the rupture length can be overestimated for events radiating high levels of high-frequency motions, as is the case in the Lushan earthquake. If the set of strong motion data used had been available in real-time, 50% to 80% of sites experiencing shaking of intensity MMI IV-VII (light to very strong) and 30% experiencing VIII-IX (severe to violent) could have been issued a warning with 10 s and 5 s, respectively, before the arrival of the destructive S-wave. We also show that loss estimates after devastating earthquakes based on the FinDer line-source are more accurate compared to a point-source model. For the Wenchuan earthquake, for example, they predict a four to six times larger number of fatalities and injured, which is consistent with official reports. At the same time, these losses could be provided 1/2~3 hours faster than if based on more complex inversion rupture models.

Sichuan↗

USGS2007

3D Visualization of Earthquake Focal Mechanisms Using ArcScene

We created a new tool, 3D Focal Mechanisms (3DFM), for viewing earthquake focal mechanism symbols three dimensionally. This tool operates within the Environmental Systems Research Institute (ESRI®) GIS software ArcScene® 9.x. The program requires as input a GIS point dataset of earthquake locations containing strike, dip, and rake values for a nodal plane of each earthquake. Other information, such as depth and magnitude of the earthquake, may also be included in the dataset. By default for each focal point, 3DFM will create a black and white sphere or “beach ball” that is oriented based on the strike, dip, and rake values. If depth values for each earthquake are included, the focal symbol will also be placed at its appropriate location beneath the Earth's surface. In addition to the default settings, there are several other options in 3DFM that can be adjusted. The appearance of the symbols can be changed by (1) creating rings around the fault planes that are colored based on magnitude, (2) showing only the fault planes instead of a sphere, (3) drawing a flat disc that identifies the primary nodal plane, (4) or by displaying the null, pressure, and tension axes. The size of the symbols can be changed by adjusting their diameter, scaling them based on the magnitude of the earthquake, or scaling them by the estimated size of the rupture patch based on earthquake magnitude. It is also possible to filter the data using any combination of the strike, dip, rake, magnitude, depth, null axis plunge, pressure axis plunge, tension axis plunge, or fault type values of the points. For a large dataset, these filters can be used to create different subsets of symbols. Symbols created by 3DFM are stored in graphics layers that appear in the ArcScene® table of contents. Multiple graphics layers can be created and saved to preserve the output from different symbol options.

Data Series↗

USGS2017

Shallow-depth location and geometry of the Piedmont Reverse splay of the Hayward Fault, Oakland, California

The Piedmont Thrust Fault, herein referred to as the Piedmont Reverse Fault (PRF), is a splay of the Hayward Fault that trends through a highly populated area of the City of Oakland, California (fig. 1A). Although the PRF is unlikely to generate a large-magnitude earthquake, slip on the PRF or high-amplitude seismic energy traveling along the PRF may cause considerable damage during a large earthquake on the Hayward Fault. Thus, it is important to determine the exact location, geometry (particularly dip), and lateral extent of the PRF within the densely populated Oakland area. In the near surface, the PRF juxtaposes Late Cretaceous sandstone (of the Franciscan Complex Novato Quarry terrane of Blake and others, 1984) and an older Pleistocene alluvial fan unit along much of its mapped length (fig. 1B; Graymer and others, 1995). The strata of the Novato Quarry unit vary greatly in strike (NW, NE, and E), dip direction (NE, SW, E, and NW), dip angle (15° to 85°), and lithology (shale and sandstone), and the unit has been intruded by quartz diorite in places. Thus, it is difficult to infer the structure of the fault, particularly at depth, with conventional seismic reflection imaging methods. To better determine the location and shallow-depth geometry of the PRF, we used high-resolution seismic imaging methods described by Catchings and others (2014). These methods involve the use of coincident P-wave (compressional wave) and S-wave (shear wave) refraction tomography and reflection data, from which tomographic models of P- and S-wave velocity and P-wave reflection images are developed. In addition, the coincident P-wave velocity (VP) and S-wave velocity (VS) data are used to develop tomographic models of VP/VS ratios and Poisson’s ratio, which are sensitive to shallow-depth faulting and groundwater. In this study, we also compare measurements of Swave velocities determined from surface waves with those determined from refraction tomography. We use the combination of seismic methods to infer the fault location, dip, and the National Earthquake Hazards Reduction Program (NEHRP) site classification along the seismic profile. Our seismic study is a smaller part of a larger study of the PRF by Trench and others (2016).

California↗

USGS1994

The Loma Prieta, California, earthquake of October 17, 1989: Tectonic processes and models

If there is a single theme that unifies the diverse papers in this chapter, it is the attempt to understand the role of the Loma Prieta earthquake in the context of the earthquake 'machine' in northern California: as the latest event in a long history of shocks in the San Francisco Bay region, as an incremental contributor to the regional deformation pattern, and as a possible harbinger of future large earthquakes. One of the surprises generated by the earthquake was the rather large amount of uplift that occurred as a result of the reverse component of slip on the southwest-dipping fault plane. Preearthquake conventional wisdom had been that large earthquakes in the region would probably be caused by horizontal, right-lateral, strike-slip motion on vertical fault planes. In retrospect, the high topography of the Santa Cruz Mountains and the elevated marine terraces along the coast should have provided some clues. With the observed ocean retreat and the obvious uplift of the coast near Santa Cruz that accompanied the earthquake, Mother Nature was finally caught in the act. Several investigators quickly saw the connection between the earthquake uplift and the long-term evolution of the Santa Cruz Mountains and realized that important insights were to be gained by attempting to quantify the process of crustal deformation in terms of Loma Prieta-type increments of northward transport and fault-normal shortening.

California↗

USGS1985

Earthquake hazards to domestic water distribution systems in Salt Lake County, Utah

A magnitude-7. 5 earthquake occurring along the central portion of the Wasatch Fault, Utah, may cause significant damage to Salt Lake County's domestic water system. This system is composed of water treatment plants, aqueducts, distribution mains, and other facilities that are vulnerable to ground shaking, liquefaction, fault movement, and slope failures. Recent investigations into surface faulting, landslide potential, and earthquake intensity provide basic data for evaluating the potential earthquake hazards to water-distribution systems in the event of a large earthquake. Water supply system components may be vulnerable to one or more earthquake-related effects, depending on site geology and topography. Case studies of water-system damage by recent large earthquakes in Utah and in other regions of the United States offer valuable insights in evaluating water system vulnerability to earthquakes.

Conference Paper↗

USGS1986

Earthquakes

An earthquake is the motion or trembling of the ground produced by sudden displacement of rock in the Earth's crust. Earthquakes result from crustal strain, volcanism, landslides, and collapse of caverns. Stress accumulates in response to tectonic forces until it exceeds the strength of the rock. The rock then breaks along a preexisting or new fracture called a fault. The rupture extends outward in all directions along the fault plane from its point of origin (focus). The rupture travels in an irregular manner until the stress is relatively equalized. If the rupture disturbs the surface, it produces a visible fault. Earthquakes can affect hundreds of thousands of square kilometers; cause damage to property measured in the tens of billions of dollars; result in loss of life and injury to hundreds of thousands of persons; and disrupt the social and economic functioning of the affected area. Although earthquakes in the United States occur most frequently in states west of the Rocky Mountains, devastating earthquakes have also occurred in the Midwest and East. All 50 states have some degree of risk from earthquakes. Many of these earthquake effects are depicted in the slides included in this set.

Report↗

USGS2014

Stress-based aftershock forecasts made within 24h post mainshock: Expected north San Francisco Bay area seismicity changes after the 2014 M=6.0 West Napa earthquake

We calculate stress changes resulting from the M = 6.0 West Napa earthquake on north San Francisco Bay area faults. The earthquake ruptured within a series of long faults that pose significant hazard to the Bay area, and we are thus concerned with potential increases in the probability of a large earthquake through stress transfer. We conduct this exercise as a prospective test because the skill of stress‐based aftershock forecasting methodology is inconclusive. We apply three methods: (1) generalized mapping of regional Coulomb stress change, (2) stress changes resolved on Uniform California Earthquake Rupture Forecast faults, and (3) a mapped rate/state aftershock forecast. All calculations were completed within 24 h after the main shock and were made without benefit of known aftershocks, which will be used to evaluative the prospective forecast. All methods suggest that we should expect heightened seismicity on parts of the southern Rodgers Creek, northern Hayward, and Green Valley faults.

California↗

USGS2004

Crustal structure of the northern margin of the eastern Tien Shan, China, and its tectonic implications for the 1906 M~7.7 Manas earthquake

The Tien Shan orogenic belt is the most active intracontinental mountain belt in the world. We describe an 86-km-long N–S-trending deep seismic reflection profile (which passes through the southern Junggar basin) located on the northeastern Tien Shan piedmont. Two distinct anticlines beneath the northern margin of the Tien Shan are clearly imaged in the seismic section. In addition, we have imaged two detachment surfaces at depths of ∼7 and ∼16 km. The detachment surface at 16-km depth corresponds to the main detachment that converges with the steep angle reverse fault (the Junggar Southern Marginal Fault) on which the 1906 M~7.7 Manas earthquake occurred. A 12–14-km-thick sedimentary basin is imaged beneath the southern Junggar basin near Shihezi. The crust beneath the northern margin of the Tien Shan is 50–55-km thick, and decreases beneath the Junggar basin to 40–45-km thick. The crustal image of the deep seismic reflection profile is consistent with models derived from nearby seismic refraction data and Bouguer gravity anomalies in the same region. The faulting associated with the 1906 Manas earthquake also fits within the structural framework imaged by the seismic reflection profile. Present-day micro-seismicity shows a hypocentral depth-distribution between 5 and 35 km, with a peak at 20 km. We hypothesize that the 1906 Manas earthquake initiated at a depth of ∼20 km and propagated upwards, causing northward slip on the sub-horizontal detachments beneath the southern Junggar basin. Thus, in accord with regional geological mapping, the current shortening within the eastern Tien Shan is accommodated both by high-angle reverse faulting and detachment faulting that can be clearly imaged at depth in seismic reflection data.

Earth and Planetary Science Letters↗

USGS2020

Mysterious tsunami in the Caribbean Sea following the 2010 Haiti earthquake possibly generated by dynamically triggered early aftershocks

Dynamically triggered offshore aftershocks, caused by passing seismic waves from main shocks located on land, are currently not considered in tsunami warnings. The M7.0 2010 Haiti earthquake epicenter was located on land 27 km north of the Caribbean Sea and its focal mechanism was oblique strike-slip. Nevertheless, a tsunami recorded on a Caribbean Deep-Ocean Assessment and Reporting of Tsunami (DART) buoy and a tide gauge produced runup heights of 1–3 m along Haiti southeast coast. Earthquake finite-fault model inversions of the DART waveform suggest that a reverse fault doublet with magnitudes of M6.8 and M6.5 located 85 km southwest of the epicenter may have excited the tsunami. This doublet collocates with dynamically triggered aftershocks, derived from back-projection analysis, that occurred 20-60 s after the main shock of the Haiti earthquake. The aftershocks are within a region of maximum dynamic strain predicted by the main shock, on a possibly tectonically active submarine ridge southwest of Haiti's Southern Peninsula. The agreement between the tsunami finite-fault source models and the seismic and tectonic evidence suggests that earthquakes on land, even strike-slip faults, can generate tsunamis by dynamically triggering offshore aftershocks.

Earth and Planetary Science Letters↗

USGS2007

Microearthquake streaks and seismicity triggered by slow earthquakes on the mobile south flank of Kilauea Volcano, Hawai'i

We perform waveform cross correlation and high precision relocation of both background seismicity and seismicity triggered by periodic slow earthquakes at Kilauea Volcano's mobile south flank. We demonstrate that the triggered seismicity dominantly occurs on several preexisting fault zones at the Hilina region. Regardless of the velocity model employed, the relocated earthquake epicenters and triggered seismicity localize onto distinct fault zones that form streaks aligned with the slow earthquake surface displacements determined from GPS. Due to the unknown effects of velocity heterogeneity and nonideal station coverage, our relocation analyses cannot distinguish whether some of these fault zones occur within the volcanic crust at shallow depths or whether all occur on the decollement between the volcano and preexisting oceanic crust at depths of ∼8 km. Nonetheless, these Hilina fault zones consistently respond to stress perturbations from nearby slow earthquakes.

Hawaii↗

USGS2009

Moderate rates of late Quaternary slip along the northwestern margin of the Basin and Range Province, Surprise Valley fault, northeastern California

The 86-km-long Surprise Valley normal fault forms part of the active northwestern margin of the Basin and Range province in northeastern California. We use trench mapping and radiocarbon, luminescence, and tephra dating to estimate displacements and timing of the past five surface-rupturing earthquakes on the central part of the fault near Cedarville. A Bayesian OxCal analysis of timing constraints indicates earthquake times of 18.2 ± 2.6, 10.9 ± 3.2, 8.5 ± 0.5, 5.8 ± 1.5, and 1.2 ± 0.1 ka. These data yield recurrence intervals of 7.3 ± 4.1, 2.5 ± 3.2, 2.7 ± 1.6, and 4.5 ± 1.5 ka and an elapsed time of 1.2 ± 0.1 ka since the latest surface-rupturing earthquake. Our best estimate of latest Quaternary vertical slip rate is 0.6 ?? 0.1 mm/a. This late Quaternary rate is remarkably similar to long-term (8-14 Ma) minimum vertical slip rates (>0.4-0.5 ± 0.3 mm/a) calculated from recently acquired seismic reflection and chronologic and structural data in Surprise Valley and the adjacent Warner Mountains. However, our slip rate yields estimates of extension that are lower than recent campaign GPS determinations by factors of 1.5-4 unless the fault has an unusually shallow (30°-35°) dip as suggested by recently acquired seismic reflection data. Coseismic displacements of 2-4.5 ± 1 m documented in the trench and probable rupture lengths of 53-65 km indicate a history of latest Quaternary earthquakes of M 6.8-7.3 on the central part of the. Surprise Valley fault.

California↗

USGS2019

Three-dimensional basin and fault structure from a detailed seismic velocity model of Coachella Valley, Southern California

The Coachella Valley in the northern Salton Trough is known to produce destructive earthquakes, making it a high seismic hazard area. Knowledge of the seismic velocity structure and geometry of the sedimentary basins and fault zones is required to improve earthquake hazard estimates in this region. We simultaneously inverted first P wave travel times from the Southern California Seismic Network (39,998 local earthquakes) and explosions (251 land/sea shots) from the 2011 Salton Seismic Imaging Project to obtain a 3-D seismic velocity model. Earthquakes with focal depths ≤10 km were selected to focus on the upper crustal structure. Strong lateral velocity contrasts in the top ~3 km correlate well with the surface geology, including the low-velocity (<5 km/s) sedimentary basin and the high-velocity crystalline basement rocks outside the valley. Sediment thickness is ~4 km in the southeastern valley near the Salton Sea and decreases to <2 km at the northwestern end of the valley. Eastward thickening of sediments toward the San Andreas fault within the valley defines Coachella Valley basin asymmetry. In the Peninsular Ranges, zones of relatively high seismic velocities (~6.4 km/s) between 2 to 4 km depth may be related to Late Cretaceous mylonite rocks or older inherited basement structures. Other high-velocity domains exist in the model down to 9 km depth and help define crustal heterogeneity. We identify a potential fault zone in Lost Horse Valley unassociated with mapped faults in Southern California from the combined interpretation of surface geology, seismicity, and lateral velocity changes in the model.

California↗

USGS2020

Near-fault velocity spectra from laboratory failures and their relation to natural ground motion

We compared near-fault velocity spectra recorded during laboratory experiments to that of natural earthquakes. We fractured crystalline rock samples at room temperature and intermediate confining pressure (50 MPa). Subsequent slip events were generated on the fracture surfaces under higher confinement (300 MPa). Velocity spectra from rock fracture resemble the inverse frequency (1/ f ) decay of natural earthquake velocity. This spectrum can be attributed to fault creation via seismic fracturing over a wide range of spatial scales. In contrast, subsequent slips on the rough fracture surfaces are depleted in high frequency energy and falloff approximately as 1/ f 2 . The 1/ f 2 spectrum is more consistent with a slider-block model obeying static-kinetic friction than a natural earthquake. The depleted high frequency content precludes the rough fault experiments from being directly analogous to natural sources. The suppression of high frequencies may have resulted from two possible factors: (1) the presence of a well-developed shear zone and coseismic damping of the fault motion by dissipation within it or, in our favored interpretation, (2) a smaller amount of energy dissipated by shearing relative to the total energy release at elevated confining pressure. In context of the latter explanation, a unifying concept that applies to these experiments, earthquakes, ground motion, and models of complex radiated motion is that high frequency radiated energy is relatively enhanced when total energy release is nearly balanced within the source region by dissipative processes. This near-critical energy release condition can be accessed at low normal stress in laboratory experiments.

Journal of Geophysical Research↗

USGS2021

Coseismic surface displacement in the 2019 ridgecrest earthquakes: Comparison of field measurements and optical image correlation results

A fundamental topic in earthquake studies is understanding the extent to which fault rupture at the surface is localized on primary fault strands as opposed to distributed tens to hundreds of meters away from primary ruptures through off‐fault deformation (OFD) via a combination of discrete secondary faulting and bulk deformation. The 2019 Ridgecrest, CA M w 6.4 and M w 7.1 earthquakes provide an opportunity to explore this problem via comparison of published field‐based and mostly on‐fault offset measurements, with new lateral displacement measurements made over length scales of hundreds of meters across the primary rupture using WorldView satellite images collected before and after the earthquakes. A mean fit to the field observations underestimates net slip (pixel‐correlation results) by an average of 0.4 m along the M w 6.4 rupture (∼41% of net left‐lateral displacement) and 0.7 m along the M w 7.1 rupture (∼65% of net right‐lateral displacement). We attribute these differences to substantial OFD along the lengths of the M w 6.4 (∼59%) and M w 7.1 (∼35%) ruptures. A maximum fit to the field observations provides an improved match to the pixel correlation results (difference of 0.1–0.2 m or ∼76%–98% of net displacement). OFD may in part depend on displacement magnitude, where locations along the rupture with smaller displacements (<2 m) are characterized by a greater percentage of OFD. Furthermore, OFD increases with rupture zone width, especially at the scale of individual rupture strands (hundreds of meters). These findings contribute to a growing understanding of the important role of OFD in accommodating deformation near the surface.

California↗