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 649 records · Page 36Linked to original sources

USGS2006

Quaternary tectonic faulting in the Eastern United States

Paleoseismological study of geologic features thought to result from Quaternary tectonic faulting can characterize the frequencies and sizes of large prehistoric and historical earthquakes, thereby improving the accuracy and precision of seismic-hazard assessments. Greater accuracy and precision can reduce the likelihood of both underprotection and unnecessary design and construction costs. Published studies proposed Quaternary tectonic faulting at 31 faults, folds, seismic zones, and fields of earthquake-induced liquefaction phenomena in the Appalachian Mountains and Coastal Plain. Of the 31 features, seven are of known origin. Four of the seven have nontectonic origins and the other three features are liquefaction fields caused by moderate to large historical and Holocene earthquakes in coastal South Carolina, including Charleston; the Central Virginia Seismic Zone; and the Newbury, Massachusetts, area. However, the causal faults of the three liquefaction fields remain unclear. Charleston has the highest hazard because of large Holocene earthquakes in that area, but the hazard is highly uncertain because the earthquakes are uncertainly located. Of the 31 features, the remaining 24 are of uncertain origin. They require additional work before they can be clearly attributed either to Quaternary tectonic faulting or to nontectonic causes. Of these 24, 14 features, most of them faults, have little or no published geologic evidence of Quaternary tectonic faulting that could indicate the likely occurrence of earthquakes larger than those observed historically. Three more features of the 24 were suggested to have had Quaternary tectonic faulting, but paleoseismological and other studies of them found no evidence of large prehistoric earthquakes. The final seven features of uncertain origin require further examination because all seven are in or near urban areas. They are the Moodus Seismic Zone (Hartford, Connecticut), Dobbs Ferry fault zone and Mosholu fault (New York City), Lancaster Seismic Zone and the epicenter of the shallow Cacoosing Valley earthquake (Lancaster and Reading, Pennsylvania), Kingston fault (central New Jersey between New York and Philadelphia), and Everona fault-Mountain Run fault zone (Washington, D.C., and Arlington and Alexandria, Virginia). ?? 2005 Elsevier B.V. All rights reserved.

Engineering Geology↗

USGS2008

Compilation of surface creep on California faults and comparison of WGCEP 2007 deformation model to Pacific-North American plate Mmtion

This Appendix contains 3 sections that 1) documents published observations of surface creep on California faults, 2) constructs line integrals across the WG-07 deformation model to compare to the Pacific - North America plate motion, and 3) constructs strain tensors of volumes across the WG-07 deformation model to compare to the Pacific - North America plate motion. Observation of creep on faults is a critical part of our earthquake rupture model because if a fault is observed to creep the moment released as earthquakes is reduced from what would be inferred directly from the fault's slip rate. There is considerable debate about how representative creep measured at the surface during a short time period is of the whole fault surface through the entire seismic cycle (e.g. Hudnut and Clark, 1989). Observationally, it is clear that the amount of creep varies spatially and temporally on a fault. However, from a practical point of view a single creep rate is associated with a fault section and the reduction in seismic moment generated by the fault is accommodated in seismic hazard models by reducing the surface area that generates earthquakes or by reducing the slip rate that is converted into seismic energy. WG-07 decided to follow the practice of past Working Groups and the National Seismic Hazard Map and used creep rate (where it was judged to be interseismic, see Table P1) to reduce the area of the fault surface that generates seismic events. In addition to following past practice, this decision allowed the Working Group to use a reduction of slip rate as a separate factor to accommodate aftershocks, post seismic slip, possible aseismic permanent deformation along fault zones and other processes that are inferred to affect the entire surface area of a fault, and thus are better modeled as a reduction in slip rate. C-zones are also handled by a reduction in slip rate, because they are inferred to include regions of widely distributed shear that is not completely expressed as earthquakes large enough to model. Because the ratio of the rate of creep relative to the total slip rate is often used to infer the average depth of creep, the depth of creep can be calculated and used to reduce the surface area of a fault that generates earthquakes in our model. This reduction of surface area of rupture is described by an aseismicity factor, assigned to each creeping fault in Appendix A. An aseismicity factor of less than 1 is only assigned to faults that are inferred to creep during the entire interseismic period. A single aseismicity factor was chosen for each section of the fault that creeps by expert opinion from the observations documented here. Uncertainties were not determined for the aseismicity factor, and thus it represents an unmodeled (and difficult to model) source of error. This Appendix simply provides the documentation of known creep, the type and precision of its measurement, and attempts to characterize the creep as interseismic, afterslip, transient or triggered. Parts 2 and 3 of this Appendix compare the WG-07 deformation model and the seismic source model it generates to the strain generated by the Pacific - North American plate motion. The concept is that plate motion generates essentially all of the elastic strain in the vicinity of the plate boundary that can be released as earthquakes. Adding up the slip rates on faults and all others sources of deformation (such as C-zones and distributed background seismicity) should approximately yield the plate motion. This addition is usually accomplished by one of four approaches: 1) line integrals that sum deformation along discrete paths through the deforming zone between the two plates, 2) seismic moment tensors that add up seismic moment of a representative set of earthquakes generated by a crustal volume spanning the plate boundary, 3) strain tensors generated by adding up the strain associated with all of the faults in a crustal volume spanning the plate

Open-File Report↗

USGS2002

Geologic and paleoseismic study of the Lavic Lake fault at Lavic Lake Playa, Mojave Desert, southern California

Paleoseismic investigations of the Lavic Lake fault at Lavic Lake playa place constraints on the timing of a possible earlier earthquake along the 1999 Hector Mine rupture trace and reveal evidence of the timing of the penultimate earthquake on a strand of the Lavic Lake fault that did not rupture in 1999. Three of our four trenches, trenches A, B, and C, were excavated across the 1999 Hector Mine rupture; a fourth trench, D, was excavated across a vegetation lineament that had only minor slip at its southern end in 1999. Trenches A-C exposed strata that are broken only by the 1999 rupture; trench D exposed horizontal bedding that is locally warped and offset by faults. Stratigraphic evidence for the timing of an earlier earthquake along the 1999 rupture across Lavic Lake playa was not exposed. Thus, an earlier event, if there was one along that rupture trace, predates the lowest stratigraphic level exposed in our trenches. Radiocarbon dating of strata near the bottom of trenches constrains a possible earlier event to some time earlier than about 4950 B.C. Buried faults revealed in trench D are below a vegetation lineament at the ground surface. A depositional contact about 80 cm below the ground surface acts as the upward termination of fault breaks in trench D. Thus, this contact may be the event horizon for a surface-rupturing earthquake prior to 1999-the penultimate earthquake on the Lavic Lake fault. Radiocarbon ages of detrital charcoal samples from immediately below the event horizon indicate that the earthquake associated with the faulting occurred later than A.D. 260. An approximately 1300-year age difference between two samples at about the same stratigraphic level below the event horizon suggests the potential for a long residence time of detrital charcoal in the area. Coupled with a lack of bioturbation that could introduce young organic material into the stratigraphic section, the charcoal ages provide only a maximum bounding age; thus, the recognized event may be younger. There is abundant, subtle evidence for pre-1999 activity of the Lavic Lake fault in the playa area, even though the fault was not mapped near the playa prior to the Hector Mine earthquake. The most notable indicators for long-term presence of the fault are pronounced, persistent vegetation lineaments and uplifted basalt exposures. Primary and secondary slip occurred in 1999 on two southern vegetation lineaments, and minor slip locally formed on a northern lineament; trench exposures across the northern vegetation lineament revealed the post-A.D. 260 earthquake, and a geomorphic trough extends northward into alluvial fan deposits in line with this lineament. The presence of two basalt exposures in Lavic lake playa indicates the presence of persistent compressional steps and uplift along the fault. Fault-line scarps are additional geomorphic markers of repeated slip events in basalt exposures.

California↗

USGS2004

The Loma Prieta, California, Earthquake of October 17, 1989 - Geologic setting and crustal structure

Although some scientists considered the Ms=7.1 Loma Prieta, Calif., earthquake of 1989 to be an anticipated event, some aspects of the earthquake were surprising. It occurred 17 km beneath the Santa Cruz Mountains along a left-stepping restraining bend in the San Andreas fault system. Rupture on the southwest-dipping fault plane consisted of subequal amounts of right-lateral and reverse motion but did not reach the surface. In the area of maximum uplift, severe shaking and numerous ground cracks occurred along Summit Road and Skyland Ridge, several kilometers south of the main trace of the San Andreas fault. The relatively deep focus of the earthquake, the distribution of ground failure, the absence of throughgoing surface rupture on the San Andreas fault, and the large component of uplift raised several questions about the relation of the 1989 Loma Prieta earthquake to the San Andreas fault: Did the earthquake actually occur on the San Andreas fault? Where exactly is the San Andreas fault in the heavily forested Santa Cruz Mountains, and how does the fault relate to ground ruptures that occurred there in 1989 and 1906? What is the geometry of the San Andreas fault system at depth, and how does it relate to the major crustal blocks identified by geologic mapping? Subsequent geophysical and geologic investigations of crustal structure in the Loma Prieta region have addressed these and other questions about the relation of the earthquake to geologic structures observed in the southern Santa Cruz Mountains. The diverse papers in this chapter cover several topics: geologic mapping of the region, potential- field and electromagnetic modeling of crustal structure, and the velocity structure of the crust and mantle in and below the source region for the earthquake. Although these papers were mostly completed between 1992 and 1997, they provide critical documentation of the crustal structure of the Loma Prieta region. Together, they present a remarkably coherent, three-dimensional picture of the earthquake source region--a geologically complex volume of crust with a long history of both right-lateral faulting and fault-normal compression, thrusting, and uplift.

California↗

USGS2012

Historic creep rate and potential for seismic slip along the Hayward Fault, California

The Hayward fault is considered the most likely source of one or more major earthquakes in the San Francisco Bay area in the next few decades. Historically, at least one, and probably two, major earthquakes (about M 6.8) occurred along the Hayward fault, one in 1836 and another in 1868. Little is known about the 1836 event, but the 1868 earthquake was accompanied by a surface rupture that extended as much as 41 km along the southern part of the fault. Although the amount of surface slip in 1868 is uncertain, right slip (including afterslip) reached at least several centimeters, and possibly several decimeters in places. This paper documents the spatial variation of creep rate along the Hayward fault since the 1868 earthquake. Creep (aseismic fault slip) occurs over at least 66 km and may extend over the fault's entire 82-km length, of which about 13 km lies underwater. Creep rate seems nearly constant over decades, but short-term variations occur. We derive creep rate mainly from our own systematic surveying of offset cultural features (curbs, fences, and buildings). On each feature we solve directly for accumulated creep by using multiple linear regression. Creep rate mostly falls in the range of 3.5–6.5 mm/yr; but systematic variation occurs along strike. Fault segments with distinctly higher and lower rates generally correspond to parts of the fault most salient from the overall average alinement of the fault. Most distinctive is a 4-km-long section near the south end of the fault that creeps at about 9 mm/yr. Such a high rate has occurred there at least since the 1920s and probably since the 1868 earthquake, as indicated by an offset railroad track built in 1869. We suggest that this 9 mm/yr slip rate may approach the long-term or deep slip rate that controls average recurrence interval between major earthquakes. If so, assuming an elastic rebound model, the potential for slip in large earthquakes below the surficial creeping zone is now &sim;1.1 m in the southern (1868) segment of the fault and ≥ 1.4 m in the northern (1836?) segment. Subtracting surface creep rates from a long-term slip rate of 9 mm/yr gives present potential for surface slip in large earthquakes of up to 0.8 m, with an average of 0.6 m in the northern segment and 0.4 m in the southern segment. We present a simple hypothesis for rupture potential that is compatible with historic creep rate, microseismicity distribution, and geodetic data. If seismic rupture occurs on segments 41 km long by 10 km deep (7 km fully locked, 3 km creeping), today's potential for seismic moment release is 1.4 × 10 19 and 1.1 × 10 19 N m for both 1836? and 1868 segments, respectively, and 2.5 × 10 19 N m for both segments jointly. Converting moment to magnitude gives M L 6.8 in the northern segment, M L 6.7 in the southern segment, and M L 7.0 for simultaneous rupture of both.

California↗

USGS2009

The SCEC/USGS dynamic earthquake rupture code verification exercise

Numerical simulations of earthquake rupture dynamics are now common, yet it has been difficult to test the validity of these simulations because there have been few field observations and no analytic solutions with which to compare the results. This paper describes the Southern California Earthquake Center/U.S. Geological Survey (SCEC/USGS) Dynamic Earthquake Rupture Code Verification Exercise, where codes that simulate spontaneous rupture dynamics in three dimensions are evaluated and the results produced by these codes are compared using Web-based tools. This is the first time that a broad and rigorous examination of numerous spontaneous rupture codes has been performed—a significant advance in this science. The automated process developed to attain this achievement provides for a future where testing of codes is easily accomplished. Scientists who use computer simulations to understand earthquakes utilize a range of techniques. Most of these assume that earthquakes are caused by slip at depth on faults in the Earth, but hereafter the strategies vary. Among the methods used in earthquake mechanics studies are kinematic approaches and dynamic approaches. The kinematic approach uses a computer code that prescribes the spatial and temporal evolution of slip on the causative fault (or faults). These types of simulations are very helpful, especially since they can be used in seismic data inversions to relate the ground motions recorded in the field to slip on the fault(s) at depth. However, these kinematic solutions generally provide no insight into the physics driving the fault slip or information about why the involved fault(s) slipped that much (or that little). In other words, these kinematic solutions may lack information about the physical dynamics of earthquake rupture that will be most helpful in forecasting future events. To help address this issue, some researchers use computer codes to numerically simulate earthquakes and construct dynamic, spontaneous rupture (hereafter called “spontaneous rupture”) solutions. For these types of numerical simulations, rather than prescribing the slip function at each location on the fault(s), just the friction constitutive properties and initial stress conditions are prescribed. The subsequent stresses and fault slip spontaneously evolve over time as part of the elasto-dynamic solution. Therefore, spontaneous rupture computer simulations of earthquakes allow us to include everything that we know, or think that we know, about earthquake dynamics and to test these ideas against earthquake observations.

Seismological Research Letters↗

USGS2007

A global search inversion for earthquake kinematic rupture history: Application to the 2000 western Tottori, Japan earthquake

[1] We present a two-stage nonlinear technique to invert strong motions records and geodetic data to retrieve the rupture history of an earthquake on a finite fault. To account for the actual rupture complexity, the fault parameters are spatially variable peak slip velocity, slip direction, rupture time and risetime. The unknown parameters are given at the nodes of the subfaults, whereas the parameters within a subfault are allowed to vary through a bilinear interpolation of the nodal values. The forward modeling is performed with a discrete wave number technique, whose Green's functions include the complete response of the vertically varying Earth structure. During the first stage, an algorithm based on the heat-bath simulated annealing generates an ensemble of models that efficiently sample the good data-fitting regions of parameter space. In the second stage (appraisal), the algorithm performs a statistical analysis of the model ensemble and computes a weighted mean model and its standard deviation. This technique, rather than simply looking at the best model, extracts the most stable features of the earthquake rupture that are consistent with the data and gives an estimate of the variability of each model parameter. We present some synthetic tests to show the effectiveness of the method and its robustness to uncertainty of the adopted crustal model. Finally, we apply this inverse technique to the well recorded 2000 western Tottori, Japan, earthquake ( Mw 6.6); we confirm that the rupture process is characterized by large slip (3-4 m) at very shallow depths but, differently from previous studies, we imaged a new slip patch (2-2.5 m) located deeper, between 14 and 18 km depth.

Tottori↗

USGS2009

Precise location of San Andreas Fault tremors near Cholame, California using seismometer clusters: Slip on the deep extension of the fault?

We examine a 24-hour period of active San Andreas Fault (SAF) tremor and show that this tremor is largely composed of repeated similar events. Utilizing this similarity, we locate the subset of the tremor with waveforms similar to an identified low frequency earthquake (LFE) "master template," located using P and S wave arrivals to be ???26 km deep. To compensate for low signal-to-noise, we estimate event-pair differential times at "clusters" of nearby stations rather than at single stations. We find that the locations form a near-linear structure in map view, striking parallel to the SAF and near the surface trace. Therefore, we suggest that at least a portion of the tremor occurs on the deep extension of the fault, likely reflecting shear slip, similar to subduction zone tremor. If so, the SAF may extend to the base of the crust, ???10 km below the deepest regular earthquakes on the fault. ?? 2009 by the American Geophysical Union.

Geophysical Research Letters↗

USGS2008

Potential effects of a scenario earthquake on the economy of southern California: Labor market exposure and sensitivity analysis to a magnitude 7.8 earthquake

The Multi-Hazards Demonstration Project (MHDP) is a collaboration between the U.S. Geological Survey (USGS) and various partners from the public and private sectors and academia, meant to improve Southern California's resiliency to natural hazards (Jones and others, 2007). In support of the MHDP objectives, the ShakeOut Scenario was developed. It describes a magnitude 7.8 (M7.8) earthquake along the southernmost 300 kilometers (200 miles) of the San Andreas Fault, identified by geoscientists as a plausible event that will cause moderate to strong shaking over much of the eight-county (Imperial, Kern, Los Angeles, Orange, Riverside, San Bernardino, San Diego, and Ventura) Southern California region. This report contains an exposure and sensitivity analysis of economic Super Sectors in terms of labor and employment statistics. Exposure is measured as the absolute counts of labor market variables anticipated to experience each level of Instrumental Intensity (a proxy measure of damage). Sensitivity is the percentage of the exposure of each Super Sector to each Instrumental Intensity level. The analysis concerns the direct effect of the scenario earthquake on economic sectors and provides a baseline for the indirect and interactive analysis of an input-output model of the regional economy. The analysis is inspired by the Bureau of Labor Statistics (BLS) report that analyzed the labor market losses (exposure) of a M6.9 earthquake on the Hayward fault by overlaying geocoded labor market data on Instrumental Intensity values. The method used here is influenced by the ZIP-code-level data provided by the California Employment Development Department (CA EDD), which requires the assignment of Instrumental Intensities to ZIP codes. The ZIP-code-level labor market data includes the number of business establishments, employees, and quarterly payroll categorized by the North American Industry Classification System. According to the analysis results, nearly 225,000 business establishments, or 44 percent of all establishments, would experience Instrumental Intensities between VII (7) and X (10). This represents more than 4 million employees earning over $45 billion in quarterly payroll. Over 57,000 of these establishments, employing over 1 million employees earning over $10 billion in quarterly payroll, would experience Instrumental Intensities of IX (9) or X (10). Based upon absolute counts and percentages, the Trade, Transportation, and Utilities Super Sector and the Manufacturing Super Sector are estimated to have the greatest exposure and sensitivity respectively. The Information and the Natural Resources and Mining Super Sectors are estimated to be the least impacted. Areas estimated to experience an Instrumental Intensity of X (10) account for approximately 3 percent of the region's labor market.

California↗

USGS2008

Potential effects of a scenario earthquake on the economy of southern California: Small business exposure and sensitivity analysis to a magnitude 7.8 earthquake

California↗

USGS2004

Earthquake scenario and probabilistic ground-shaking hazard maps for the Albuquerque-Belen-Santa Fe, New Mexico, corridor

New Mexico's population is concentrated along the corridor that extends from Belen in the south to Española in the north and includes Albuquerque and Santa Fe. The Rio Grande rift, which encompasses the corridor, is a major tectonically, volcanically, and seismically active continental rift in the western U.S. Although only one large earthquake (moment magnitude (M) ≥ 6) has possibly occurred in the New Mexico portion of the rift since 1849, paleoseismic data indicate that prehistoric surface-faulting earthquakes of M 6.5 and greater have occurred on aver- age every 400 yrs on many faults throughout the Rio Grande rift. We have developed a series of nine scenario and probabilistic hazard maps that portray the ground shaking that could occur in the Albuquerque-Belen-Santa Fe corridor from future earthquakes in New Mexico. These maps, at a scale of 1:500,000, display color-contoured ground-motion values in terms of the parameters of peak horizontal acceleration and horizontal spectral accelerations at 0.2 and 1.0 second (sec) periods. The maps depict surficial ground shaking and incorporate the site-response effects at locations underlain by unconsolidated sediments. The scenario maps are for a M 7.0 earthquake rupturing the Sandia-Rincon faults, which are adjacent to and dip west beneath Albuquerque. The probabilistic maps are for the two annual exceedance probabilities of building code relevance, 10% and 2% exceedance probabilities in 50 yrs (corresponding to return periods of 500 and 2,500 yrs, respectively). We included 57 Quaternary faults, all located within the Rio Grande rift, in the probabilistic seismic hazard analysis. These faults were characterized in terms of their geometry, rupture behavior (including possible segmentation), maximum expected earthquake magnitude, recurrence model, probability of activity, and slip rate. Preferred maximum magnitude values for these faults ranged from M 6.1 to 7.4 and preferred slip rates from 0.01 to approximately 0.12 mm/yr. Regional source zones and Gaussian smoothing of the historical seismicity were also included in the probabilistic hazard analysis to account for the hazard from background earthquakes (M ≤ 6.5). A numerical ground-motion modeling approach and empirical attenuation relation- ships appropriate for extensional tectonic regimes were used to compute the scenario earthquake and probabilistic ground motions on rock. Amplification factors were then used to modify the rock motions and hence to incorporate site response into the hazard maps. These factors were based on three generalized geologic site-response categories (hard rock, soft rock, and firm/stiff soil) and were adopted from similar California-based categories because insufficient subsurface geologic and geotechnical data are available for the map area. The resulting hazard maps indicate that from both scenario and probabilistic perspectives, the ground-shaking hazard in the Albuquerque–Belen–Santa Fe corridor from future earthquakes could be severe, damaging, and potentially disastrous. In the event of a M 7.0 earthquake occurring on the Sandia–Rincon faults, ground shaking as characterized by peak ground acceleration could reach 0.7 g in much of the eastern half of the Albuquerque metropolitan area. (1 g = 980 cm/sec, the rate of gravitational acceleration.) These high ground motions will be attributable to the city’s location directly over the Sandia–Rincon faults and the amplifying effect of the unconsolidated sediments within the Albuquerque Basin. These levels of ground shaking will probably result in severe damage to traditional adobe construction and even to modern buildings. Long- period ground motions (> 1.0 sec), which are significant to long and tall structures (e.g., tall buildings, long bridges, and highway overpasses), will also be high (> 1.0 g). Injuries and loss of life will be likely. For the 500- and 2,500-yr return period maps, the highest peak accelerations are predicted to be at the damaging levels of 0.3 g and 0.6 g, respectively. All maps show dramatically the frequency-dependent amplification of unconsolidated sediments in the basins along the Rio Grande valley (e.g., Albuquerque Basin). The pattern of amplification and deamplification is clearly a function of the distribution of unconsolidated sediments. These maps are not intended to be a substitute for site-specific studies for engineering design nor to replace standard maps commonly referenced in building codes. Rather, we hope that these maps will be used as a guide by government agencies; the engineering, urban planning, emergency preparedness, and response communities; and the general public as part of an overall program to reduce earthquake risk and losses in New Mexico.

New Mexico↗

USGS2024

Mesoproterozoic to Paleozoic tectonics, Pleistocene landforms, and Holocene seismicity in the Blue Ridge: Results from integrated studies of the 9 August 2020, Mw 5.1 earthquake area near Sparta, North Carolina, USA

This field trip examines the results of integrated geologic studies of the 9 August 2020, M w 5.1 earthquake near Sparta, North Carolina, USA. The earthquake generated ~4 km of coseismic surface rupture of the Little River fault and uplifted a surface area of ~11 km 2 . The Little River fault is a thrust fault oriented 110–130°/45–70°SW, and mapped fault segments are en echelon with scarp heights from <5–30 cm. The epicenter is in polydeformed rocks of the Ashe and Alligator Back Metamorphic Suites in the eastern Blue Ridge. Bedrock structure formed during multiple Paleozoic orogenies; the regional foliation strikes NE-SW and dips SE (mean orientation 063°/52°SE). Mapping identified late Paleozoic veins and shear zones, a regional joint set striking 330–340° and 250–240°, and brittle faults that cut the Paleozoic foliation. Brittle faults oriented similar to the Little River fault are mapped up to 4 km along strike from the coseismic rupture along Bledsoe Creek valley, and the combined length of the Little River fault system is ~8 km. Paleoseismic trenches across the Little River fault corroborate the reactivation of an older fault by the 2020 earthquake and reveal two events during late Pleistocene (<50 ka). Surficial mapping identified several terrace deposits, including a deposit along Bledsoe Creek that yielded a 26 Al/ 10 Be isochron burial age of 0.46 ± 0.13 Ma and overlies a brittle fault, thus constraining the timing of movement of the fault at that location. Paleoliquefaction studies document soft-sediment deformation features in alluvium that may represent paleoseismic events. Collectively, these results highlight long-lived paleoseismicity of the Blue Ridge and that the 9 August 2020 earthquake reactivated an older, suitably oriented brittle fault in the bedrock. The Little River fault is an example of a previously unknown but active fault lying outside of known seismic zones with demonstrated recurrence of paleo-ruptures, raising questions about the assumption that damaging earthquakes are limited to areas of ongoing background seismicity, which is counter to seismic hazard assessments in the eastern United States. Bedrock mapping separates eastern Blue Ridge lithostratigraphy of the Lynchburg Group and Ashe and Alligator Back Metamorphic Suites into separate fault-bound packages juxtaposed over various 1.3–1.0 Ga basement rocks of the northern French Broad massif by the Gossan Lead fault.

North Carolina↗

USGS2016

Statistical tests of simple earthquake cycle models

A central goal of observing and modeling the earthquake cycle is to forecast when a particular fault may generate an earthquake: a fault late in its earthquake cycle may be more likely to generate an earthquake than a fault early in its earthquake cycle. Models that can explain geodetic observations throughout the entire earthquake cycle may be required to gain a more complete understanding of relevant physics and phenomenology. Previous efforts to develop unified earthquake models for strike-slip faults have largely focused on explaining both preseismic and postseismic geodetic observations available across a few faults in California, Turkey, and Tibet. An alternative approach leverages the global distribution of geodetic and geologic slip rate estimates on strike-slip faults worldwide. Here we use the Kolmogorov-Smirnov test for similarity of distributions to infer, in a statistically rigorous manner, viscoelastic earthquake cycle models that are inconsistent with 15 sets of observations across major strike-slip faults. We reject a large subset of two-layer models incorporating Burgers rheologies at a significance level of α = 0.05 (those with long-term Maxwell viscosities η M <~ 4.0 × 10 19 Pa s and η M >~ 4.6 × 10 20 Pa s) but cannot reject models on the basis of transient Kelvin viscosity η K . Finally, we examine the implications of these results for the predicted earthquake cycle timing of the 15 faults considered and compare these predictions to the geologic and historical record.

Geophysical Research Letters↗

USGS2013

InSAR Evidence for an active shallow thrust fault beneath the city of Spokane Washington, USA

In 2001, a nearly five month long sequence of shallow, mostly small magnitude earthquakes occurred beneath the city of Spokane, a city with a population of about 200,000, in the state of Washington. During most of the sequence, the earthquakes were not well located because seismic instrumentation was sparse. Despite poor-quality locations, the earthquake hypocenters were likely very shallow, because residents near the city center both heard and felt many of the earthquakes. The combination of poor earthquake locations and a lack of known surface faults with recent movement make assessing the seismic hazards related to the earthquake swarm difficult. However, the potential for destruction from a shallow moderate-sized earthquake is high, for example Christchurch New Zealand in 2011, so assessing the hazard potential of a seismic structure involved in the Spokane earthquake sequence is important. Using interferometric synthetic aperture radar (InSAR) data from the European Space Agency ERS2 and ENVISAT satellites and the Canadian Space Agency RADARSAT-1, satellite we are able to show that slip on a shallow previously unknown thrust fault, which we name the Spokane Fault, is the source of the earthquake sequence. The part of the Spokane Fault that slipped during the 2001 earthquake sequence underlies the north part of the city, and slip on the fault was concentrated between ~0.3 and 2 km depth. Projecting the buried fault plane to the surface gives a possible surface trace for the Spokane Fault that strikes northeast from the city center into north Spokane.

Washington↗

USGS2009

The 2006-2007 Kuril Islands great earthquake sequence

[1] The southwestern half of a ∼500 km long seismic gap in the central Kuril Island arc subduction zone experienced two great earthquakes with extensive preshock and aftershock sequences in late 2006 to early 2007. The nature of seismic coupling in the gap had been uncertain due to the limited historical record of prior large events and the presence of distinctive upper plate, trench and outer rise structures relative to adjacent regions along the arc that have experienced repeated great interplate earthquakes in the last few centuries. The intraplate region seaward of the seismic gap had several shallow compressional events during the preceding decades (notably an M S 7.2 event on 16 March 1963), leading to speculation that the interplate fault was seismically coupled. This issue was partly resolved by failure of the shallow portion of the interplate megathrust in an M W = 8.3 thrust event on 15 November 2006. This event ruptured ∼250 km along the seismic gap, just northeast of the great 1963 Kuril Island (M w = 8.5) earthquake rupture zone. Within minutes of the thrust event, intense earthquake activity commenced beneath the outer wall of the trench seaward of the interplate rupture, with the larger events having normal-faulting mechanisms. An unusual double band of interplate and intraplate aftershocks developed. On 13 January 2007, an M W = 8.1 extensional earthquake ruptured within the Pacific plate beneath the seaward edge of the Kuril trench. This event is the third largest normal-faulting earthquake seaward of a subduction zone on record, and its rupture zone extended to at least 33 km depth and paralleled most of the length of the 2006 rupture. The 13 January 2007 event produced stronger shaking in Japan than the larger thrust event, as a consequence of higher short-period energy radiation from the source. The great event aftershock sequences were dominated by the expected faulting geometries; thrust faulting for the 2006 rupture zone, and normal faulting for the 2007 rupture zone. A large intraplate compressional event occurred on 15 January 2009 (M w = 7.4) near 45 km depth, below the rupture zone of the 2007 event and in the vicinity of the 16 March 1963 compressional event. The fault geometry, rupture process and slip distributions of the two great events are estimated using very broadband teleseismic body and surface wave observations. The occurrence of the thrust event in the shallowest portion of the interplate fault in a region with a paucity of large thrust events at greater depths suggests that the event removed most of the slip deficit on this portion of the interplate fault. This great earthquake doublet demonstrates the heightened seismic hazard posed by induced intraplate faulting following large interplate thrust events. Future seismic failure of the remainder of the seismic gap appears viable, with the northeastern region that has also experienced compressional activity seaward of the megathrust warranting particular attention.

Journal of Geophysical Research B: Solid Earth↗

USGS2002

The 1999 Izmit, Turkey, earthquake: A 3D dynamic stress transfer model of intraearthquake triggering

Before the August 1999 Izmit (Kocaeli), Turkey, earthquake, theoretical studies of earthquake ruptures and geological observations had provided estimates of how far an earthquake might jump to get to a neighboring fault. Both numerical simulations and geological observations suggested that 5 km might be the upper limit if there were no transfer faults. The Izmit earthquake appears to have followed these expectations. It did not jump across any step-over wider than 5 km and was instead stopped by a narrower step-over at its eastern end and possibly by a stress shadow caused by a historic large earthquake at its western end. Our 3D spontaneous rupture simulations of the 1999 Izmit earthquake provide two new insights: (1) the west- to east-striking fault segments of this part of the North Anatolian fault are oriented so as to be low-stress faults and (2) the easternmost segment involved in the August 1999 rupture may be dipping. An interesting feature of the Izmit earthquake is that a 5-km-long gap in surface rupture and an adjacent 25° restraining bend in the fault zone did not stop the earthquake. The latter observation is a warning that significant fault bends in strike-slip faults may not arrest future earthquakes.

Izmit↗

USGS2023

Combinatorial optimization of earthquake spatial distributions under minimum cumulative stress constraints

We determine optimal on‐fault earthquake spatial distributions using a combinatorial method that minimizes the long‐term cumulative stress resolved on the fault. An integer‐programming framework was previously developed to determine the optimal arrangement of a millennia‐scale earthquake sample that minimizes the misfit to a target slip rate determined from geodetic data. The resulting cumulative stress from just slip‐rate optimization, however, can greatly exceed fault strength estimates. Therefore, we add an objective function that minimizes cumulative stress and broad stress constraints to limit the solution space. We find that there is a trade‐off in the two objectives: minimizing the cumulative stress on a fault within fault strength limits concentrates earthquakes in specific areas of the fault and results in excursions from the target slip rate. Both slip‐rate and stress objectives can be combined in either a weighted or lexicographic (hierarchical) method. Using a combination of objectives, we demonstrate that a Gutenberg–Richter sample of earthquakes can be arranged on a constant slip‐rate finite fault with minimal stress and slip‐rate residuals. We apply this method to determine the optimal arrangement of earthquakes on the variable slip‐rate Nankai megathrust over 5000 yr. The sharp decrease in slip rate at the Tokai section of the fault results in surplus cumulative stress under all scenarios. Using stress optimization alone restricts this stress surplus to the northeast end of the fault at the expense of decreasing the slip rate away from the target slip rate at the southwest end of the fault. A combination of both slip‐rate and stress objectives provides an adequate fit to the data, although alternate model formulations for the fault are needed at the Tokai section to explain persistent excess cumulative stress. In general, incorporating stress objectives and constraints into the integer‐programming framework adds an important aspect of fault physics to the resulting earthquake rupture forecasts.

Bulletin of the Seismological Society of America↗

USGS1999

Abrupt along-strike change in tectonic style: San Andreas fault zone, San Francisco Peninsula

Seismicity and high-resolution aeromagnetic data are used to define an abrupt change from compressional to extensional tectonism within a 10- to 15-km-wide zone along the San Andreas fault on the San Francisco Peninsula and offshore from the Golden Gate. This 100-km-long section of the San Andreas fault includes the hypocenter of the M w = 7.8 1906 San Francisco earthquake as well as the highest level of persistent microseismicity along that ∼470-km-long rupture. We define two distinct zones of deformation along this stretch of the fault using well-constrained relocations of all post-1969 earthquakes based a joint one-dimensional velocity/hypocenter inversion and a redetermination of focal mechanisms. The southern zone is characterized by thrust- and reverse-faulting focal mechanisms with NE trending P axes that indicate “fault-normal” compression in 7- to 10-km-wide zones of deformation on both sides of the San Andreas fault. A 1- to 2-km-wide vertical zone beneath the surface trace of the San Andreas is characterized by its almost complete lack of seismicity. The compressional deformation is consistent with the young, high topography of the Santa Cruz Mountains/Coast Ranges as the San Andreas fault makes a broad restraining left bend (∼10°) through the southernmost peninsula. A zone of seismic quiescence ∼15 km long separates this compressional zone to the south from a zone of combined normal-faulting and strike-slip-faulting focal mechanisms (including a M L = 5.3 earthquake in 1957) on the northernmost peninsula and offshore on the Golden Gate platform. Both linear pseudogravity gradients, calculated from the aeromagnetic data, and seismic reflection data indicate that the San Andreas fault makes an abrupt ∼3-km right step less than 5 km offshore in this northern zone. A similar right-stepping (dilatational) geometry is also observed for the subparallel San Gregorio fault offshore. Persistent seismicity and extensional tectonism occur within the San Andreas right stepover region and at least 15 km along-strike both to the SE and NW. The 1906 San Francisco earthquake may have nucleated within the San Andreas right stepover, which may help explain the bilateral nature of rupture of this event. Our analysis suggests two seismic hazards for the San Francisco Peninsula in addition to the hazard associated with a M = 7 to 8 strike-slip earthquake along the San Andreas fault: the potential for a M ≈ 6 normal-faulting earthquake just 5–8 km west of San Francisco and a M = 6+ thrust faulting event in the southern peninsula.

Journal of Geophysical Research B: Solid Earth↗