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,783 records · Page 99Linked to original sources

USGS2025

Reply to, “Comment on ‘The 1886 Charleston, South Carolina, earthquake: Relic railroad offset reveals rupture,’ by Roger Bilham and Susan E. Hough”

We welcome this opportunity to respond to Pratt et al. (2024) (hereinafter P24). Bilham and Hough (2023) proposed a “first-cut” elastic deformation model for the 1886 earthquake, a quantitative source model constrained by identified coseismic constraints. A key observation was the measurement of a lateral offset of a railroad line south of Summerville, leading to a model with predominately dextral slip and minor convergence, from which we concluded that active faulting had raised the Penholoway Marine Terrace >6 m since ∼770 ka. P24 questioned these constraints and proposed an alternative rupture model with predominantly reverse slip. This alternative model is neither consistent with coseismic constraints nor with other geophysical data. In a revised model presented here, we recognize that uplift of the Penholoway Terrace is confined to the eastern edge of the terrace, which we conclude results from active folding and tectonic transpression centered on the dextral fault that offset the railroad in 1886.

South Carolina↗

USGS1968

Earthquake of August 19, 1966, Varto Area, eastern Turkey

The earthquake of August 19, 1966, centered near Varto in eastern Turkey. It occurred at 12:22:09.6 G. M. T. and had a Richter magnitude of about 7. Intensities reached IX. Reportedly, 2,529 people were killed, 1,500 injured, and 19,013 buildings were demolished or heavily damaged. Surface fractures in three parallel fracture zones trending N65-70°W display both right-lateral strike slip and extension normal to the fractures. The trend and sense of displacement relate the fracturing to the North Anatolia fault system. Examples of geologic factors influencing damage include: (1) degree of saturation of the alluvium, (2) proximity to landslides, and (3) relation to surface faulting.

Near Varto in eastern Turkey↗

USGS1996

Directional topographic site response at Tarzana observed in aftershocks of the 1994 Northridge, California, earthquake: Implications for mainshock motions

The Northridge earthquake caused 1.78 g acceleration in the east-west direction at a site in Tarzana, California, located about 6 km south of the mainshock epicenter. The accelerograph was located atop a hill about 15-m high, 500-m long, and 130-m wide, striking about N78°E. During the aftershock sequence, a temporary array of 21 three-component geophones was deployed in six radial lines centered on the accelerograph, with an average sensor spacing of 35 m. Station C00 was located about 2 m from the accelerograph. We inverted aftershock spectra to obtain average relative site response at each station as a function of direction of ground motion. We identified a 3.2-Hz resonance that is a transverse oscillation of the hill (a directional topographic effect). The top/base amplification ratio at 3.2 Hz is about 4.5 for horizontal ground motions oriented approximately perpendicular to the long axis of the hill and about 2 for motions parallel to the hill. This resonance is seen most strongly within 50 m of C00. Other resonant frequencies were also observed. A strong lateral variation in attenuation, probably associated with a fault, caused substantially lower motion at frequencies above 6 Hz at the east end of the hill. There may be some additional scattered waves associated with the fault zone and seen at both the base and top of the hill, causing particle motions (not spectral ratios) at the top of the hill to be rotated about 20° away from the direction transverse to the hill. The resonant frequency, but not the amplitude, of our observed topographic resonance agrees well with theory, even for such a low hill. Comparisons of our observations with theoretical results indicate that the 3D shape of the hill and its internal structure are important factors affecting its response. The strong transverse resonance of the hill does not account for the large east-west mainshock motions. Assuming linear soil response, mainshock east-west motions at the Tarzana accelerograph were amplified by a factor of about 2 or less compared with sites at the base of the hill. Probable variations in surficial shear-wave velocity do not account for the observed differences among mainshock acceleration observed at Tarzana and at two different sites within 2 km of Tarzana.

California↗

USGS1986

Effects of the 1906 Earthquake on the Bald Hill Outlet System, San Mateo County, California

Following the earthquake of April 18, 1906, it was discovered that a brick forebay and other parts of the reservoir outlet system were in the slip zone of the San Andreas fault. The original outlet through which water was directed to San Francisco consisted of two tunnels joined at the brick forebay; one tunnel extends 2,820 ft to the east under Bald Hill on Buri Buri Ridge, and the other tunnel intersects the lake bottom about 250 ft west of the forebay. In 1897 a second intake was added to the system, also joining the original forebay. During the present study the accessible parts of this original outlet system were examined with the hope of learning how the system had been affected by fault slip in 1906.

California↗

USGS1974

Faults and their potential hazards in Santa Cruz County, California

This report and map area designed to provide Santa Cruz County with basic data on the location, pattern, recency of movement, potential for future surface rupture, and anticipated earthquake magnitudes and recurrence intervals for several of the faults located there for consideration in the county's land-use planning program and the preparation of its seismic safety element. This report does not consider other seismic hazards such as ground shaking, seismically induced ground failure, or tsunamis (seismically generated sea waves) that can accompany offshore earthquakes of large magnitude.

California↗

USGS1981

Digital seismograms of aftershocks of the Imperial Valley, California, earthquake of October 15, 1979

This report presents the digital seismograms collected for the aftershocksof the Imperial Valley, California earthquake (M L = 6.6) that occurred onL October 15, 1979. A short summary follows that describes the field procedures, instrumentation, and routine computer analysis that was used to produce the ground motion time histories. Boore and Fletcher (1981) may be referenced for selected fault plane solutions, spectra and moments. Requests for digital time histories on magnetic tape may be made to the Branch of Ground Motion and Faulting, Menlo Park, California.

California↗

USGS2014

Surface deformation associated with the March 1996 earthquake swarm at Akutan Island, Alaska, revealed by C-band ERS and L-band JERS radar interferometry

In March 1996, an intense earthquake swarm beneath Akutan Island, Alaska, was accompanied by extensive ground cracking but no eruption of Akutan volcano. Radar interferograms produced from L-band JERS-1 and C-band ERS-1/2 images show uplift associated with the swarm by as much as 60 cm on the western part of the island. The JERS-1 interferogram has greater coherence, especially in areas with loose surface material or thick vegetation. It also shows subsidence of similar magnitude on the eastern part of the island and displacements along faults reactivated during the swarm. The axis of uplift and subsidence strikes about N70°W, which is roughly parallel to a zone of fresh cracks on the northwest flank of the volcano, to normal faults that cut the island and to the inferred maximum compressive stress direction. A common feature of models that fit the deformation is the emplacement of a shallow dike along this trend beneath the northwest flank of the volcano. Both before and after the swarm, the northwest flank was uplifted 5–20 mm/year relative to the southwest flank, probably by magma intrusion. The zone of fresh cracks subsided about 20 mm during 1996–1997 and at lesser rates thereafter, possibly because of cooling and degassing of the intrusion.

Alaska↗

USGS1996

Rupture directivity and slip distribution of the M 4.3 foreshock to the 1992 Joshua Tree earthquake, Southern California

Details of the M 4.3 foreshock to the Joshua Tree earthquake were studied using P waves recorded on the Southern California Seismic Network and the Anza network. Deconvolution, using an M 2.4 event as an empirical Green's function, corrected for complicated path and site effects in the seismograms and produced simple far-field displacement pulses that were inverted for a slip distribution. Both possible fault planes, north-south and east-west, for the focal mechanism were tested by a least-squares inversion procedure with a range of rupture velocities. The results showed that the foreshock ruptured the north-south plane, similar to the mainshock. The foreshock initiated a few hundred meters south of the mainshock and ruptured to the north, toward the mainshock hypocenter. The mainshock (M 6.1) initiated near the northern edge of the foreshock rupture 2 hr later. The foreshock had a high stress drop (320 to 800 bars) and broke a small portion of the fault adjacent to the mainshock but was not able to immediately initiate the mainshock rupture.

Bulletin of the Seismological Society of America↗

USGS2016

Earthquake probabilities for the Wasatch front region in Utah, Idaho, and Wyoming

In a letter to The Salt Lake Daily Tribune in September 1883, U.S. Geological Survey (USGS) geologist G.K. Gilbert warned local residents about the implications of observable fault scarps along the western base of the Wasatch Range. The scarps were evidence that large surface-rupturing earthquakes had occurred in the past and more would likely occur in the future. The main actor in this drama is the 350-km-long Wasatch fault zone (WFZ), which extends from central Utah to southernmost Idaho. The modern Wasatch Front urban corridor, which follows the valleys on the WFZ’s hanging wall between Brigham City and Nephi, is home to nearly 80% of Utah’s population of 3 million. Adding to this circumstance of “lots of eggs in one basket,” more than 75% of Utah’s economy is concentrated along the Wasatch Front in Utah’s four largest counties, literally astride the five central and most active segments of the WFZ.

Utah, Idaho, Wyoming↗

USGS2016

Effective stress, friction and deep crustal faulting

Studies of crustal faulting and rock friction invariably assume the effective normal stress that determines fault shear resistance during frictional sliding is the applied normal stress minus the pore pressure. Here we propose an expression for the effective stress coefficient α f at temperatures and stresses near the brittle-ductile transition (BDT) that depends on the percentage of solid-solid contact area across the fault. α f varies with depth and is only near 1 when the yield strength of asperity contacts greatly exceeds the applied normal stress. For a vertical strike-slip quartz fault zone at hydrostatic pore pressure and assuming 1 mm and 1 km shear zone widths for friction and ductile shear, respectively, the BDT is at ~13 km. α f near 1 is restricted to depths where the shear zone is narrow. Below the BDT α f = 0 is due to a dramatically decreased strain rate. Under these circumstances friction cannot be reactivated below the BDT by increasing the pore pressure alone and requires localization. If pore pressure increases and the fault localizes back to 1 mm, then brittle behavior can occur to a depth of around 35 km. The interdependencies among effective stress, contact-scale strain rate, and pore pressure allow estimates of the conditions necessary for deep low-frequency seismicity seen on the San Andreas near Parkfield and in some subduction zones. Among the implications are that shear in the region separating shallow earthquakes and deep low-frequency seismicity is distributed and that the deeper zone involves both elevated pore fluid pressure and localization.

Journal of Geophysical Research B: Solid Earth↗

USGS2016

Coseismic slip and early afterslip of the 2015 Illapel, Chile, earthquake: Implications for frictional heterogeneity and coastal uplift

Great subduction earthquakes are thought to rupture portions of the megathrust, where interseismic coupling is high and velocity-weakening frictional behavior is dominant, releasing elastic deformation accrued over a seismic cycle. Conversely, postseismic afterslip is assumed to occur primarily in regions of velocity-strengthening frictional characteristics that may correlate with lower interseismic coupling. However, it remains unclear if fixed frictional properties of the subduction interface, coseismic or aftershock-induced stress redistribution, or other factors control the spatial distribution of afterslip. Here we use interferometric synthetic aperture radar and Global Position System observations to map the distribution of coseismic slip of the 2015 M w 8.3 Illapel, Chile, earthquake and afterslip within the first 38 days following the earthquake. We find that afterslip overlaps the coseismic slip area and propagates along-strike into regions of both high and moderate interseismic coupling. The significance of these observations, however, is tempered by the limited resolution of geodetic inversions for both slip and coupling. Additional afterslip imaged deeper on the fault surface bounds a discrete region of deep coseismic slip, and both contribute to net uplift of the Chilean Coastal Cordillera. A simple partitioning of the subduction interface into regions of fixed frictional properties cannot reconcile our geodetic observations. Instead, stress heterogeneities, either preexisting or induced by the earthquake, likely provide the primary control on the afterslip distribution for this subduction zone earthquake. We also explore the occurrence of coseismic and postseismic coastal uplift in this sequence and its implications for recent hypotheses concerning the source of permanent coastal uplift along subduction zones.

Journal of Geophysical Research B: Solid Earth↗

USGS2016

A study of the 2015 Mw 8.3 Illapel earthquake and tsunami: Numerical and analytical approaches

The September 16, 2015 Illapel, Chile earthquake triggered a large tsunami, causing both economic losses and fatalities. To study the coastal effects of this earthquake, and to understand how such hazards might be accurately modeled in the future, different finite fault models of the Illapel rupture are used to define the initial condition for tsunami simulation. The numerical code Non-hydrostatic Evolution of Ocean WAVEs (NEOWAVE) is employed to model the tsunami evolution through the Pacific Ocean. Because only a short time is available for emergency response, and since the earthquake and tsunami sources are close to the coast, gaining a rapid understanding of the near-field run-up behavior is highly relevant to Chile. Therefore, an analytical solution of the 2 ? 1 D shallow water wave equations is considered. With this solution, we show that we can quickly estimate the run-up distribution along the coastline, to first order. After the earthquake and tsunami, field observations were measured in the surrounded coastal region, where the tsunami resulted in significant run-up. First, we compare the analytical and numerical solutions to test the accuracy of the analytical approach and the field observations, implying the analytic approach can accurately model tsunami run-up after an earthquake, without sacrificing the time necessary for a full numerical inversion. Then, we compare both with field run-up measurements. We observe the consistency between the two approaches. To complete the analysis, a tsunami source inversion is performed using run-up field measurements only. These inversion results are compared with seismic models, and are shown to capture the broad-scale details of those models, without the necessity of the detailed data sets they invert.

Illapel↗

USGS2001

Variations in creep rate along the Hayward Fault, California, interpreted as changes in depth of creep

Variations ill surface creep rate along the Hayward fault are modeled as changes in locking depth using 3D boundary elements. Model creep is driven by screw dislocations at 12 km depth under the Hayward and other regional faults. Inferred depth to locking varies along strike from 4-12 km. (12 km implies no locking.) Our models require locked patches under the central Hayward fault, consistent with a M6.8 earthquake in 1868, but the geometry and extent of locking under the north and south ends depend critically on assumptions regarding continuity and creep behavior of the fault at its ends. For the northern onshore part of the fault, our models contain 1.4-1.7 times more stored moment than the model of Bu??rgmann et al. [2000]; 45-57% of this stored moment resides in creeping areas. It is important for seismic hazard estimation to know how much of this moment is released coseismically or as aseismic afterslip.

Geophysical Research Letters↗

USGS2008

Constraining fault constitutive behavior with slip and stress heterogeneity

We study how enforcing self-consistency in the statistical properties of the preshear and postshear stress on a fault can be used to constrain fault constitutive behavior beyond that required to produce a desired spatial and temporal evolution of slip in a single event. We explore features of rupture dynamics that (1) lead to slip heterogeneity in earthquake ruptures and (2) maintain these conditions following rupture, so that the stress field is compatible with the generation of aftershocks and facilitates heterogeneous slip in subsequent events. Our three-dimensional fmite element simulations of magnitude 7 events on a vertical, planar strike-slip fault show that the conditions that lead to slip heterogeneity remain in place after large events when the dynamic stress drop (initial shear stress) and breakdown work (fracture energy) are spatially heterogeneous. In these models the breakdown work is on the order of MJ/m2, which is comparable to the radiated energy. These conditions producing slip heterogeneity also tend to produce narrower slip pulses independent of a slip rate dependence in the fault constitutive model. An alternative mechanism for generating these confined slip pulses appears to be fault constitutive models that have a stronger rate dependence, which also makes them difficult to implement in numerical models. We hypothesize that self-consistent ruptures could also be produced by very narrow slip pulses propagating in a self-sustaining heterogeneous stress field with breakdown work comparable to fracture energy estimates of kJ/M2. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2011

Testing earthquake source inversion methodologies

Source Inversion Validation Workshop; Palm Springs, California, 11-12 September 2010; Nowadays earthquake source inversions are routinely performed after large earthquakes and represent a key connection between recorded seismic and geodetic data and the complex rupture process at depth. The resulting earthquake source models quantify the spatiotemporal evolution of ruptures. They are also used to provide a rapid assessment of the severity of an earthquake and to estimate losses. However, because of uncertainties in the data, assumed fault geometry and velocity structure, and chosen rupture parameterization, it is not clear which features of these source models are robust. Improved understanding of the uncertainty and reliability of earthquake source inversions will allow the scientific community to use the robust features of kinematic inversions to more thoroughly investigate the complexity of the rupture process and to better constrain other earthquakerelated computations, such as ground motion simulations and static stress change calculations.

Eos, Transactions, American Geophysical Union↗

USGS2020

Semiautomated estimates of directivity and related source properties of small to moderate southern California earthquakes using second seismic moments

We develop a semiautomated method for estimating with second seismic moments the directivity, rupture area, duration, and centroid velocity of earthquakes. The method is applied to 41 southern California earthquakes with magnitude in the range 3.5–5.2 and provides stable results for 28 events. Apparent source time functions (ASTFs) of P and S phases are derived using deconvolution with three stacked empirical Green's functions (seGf). The use of seGf suppresses nongeneric source effects, improves the focal mechanism correspondence to the analyzed earthquakes, and typically allows inclusion of 5 to 15 more ASTFs compared with analysis using a single eGf. Most analyzed earthquakes in the Trifurcation area of the San Jacinto Fault have directivities toward the northwest, while events around Cajon Pass and San Gabriel Mountain tend to propagate toward the southeast. These results are generally consistent with predictions for dynamic rupture on bimaterial interfaces associated with the imaged velocity contrasts in the area. The second moment inversions also provide constraints on the upper and lower bounds of rupture areas in our data set. Stress drops and uncertainties are estimated for elliptical ruptures using the derived characteristic rupture length and width. The semiautomated second moment method with seGfs can be used for routine application to moderate earthquakes in locations with good station coverage.

California↗

USGS2021

Can modeling the geologic record contribute to constraining the tectonic source of the AD 1755 Great Lisbon earthquake?

The precise location of the seismic source of 1755 CE Great Lisbon earthquake is still uncertain. The aim of this work is to use an onland sedimentary record in southern Portugal to test and validate seismic sources for the earthquake. To achieve this, tsunami deposit thicknesses from over 150 cores collected at Salgados in southern Portugal were compared to the results of a tsunami sediment transport model (Delft3D-FLOW) that simulates tsunami propagation, inundation, erosion, and deposition. Five different hypothetical seismic sources were modeled with varying bed roughness coefficients to assess how well they reproduced observed patterns of tsunami deposit thicknesses and dune. Modeled and observed historical tsunami arrival times were also used to test different earthquake sources. Based on these comparisons, three modeled earthquake sources were able to reproduce the observed data, suggesting they should be regarded as somewhat more likely sources for the 1755 earthquake in contrast to four other modeled sources. The fault closest to shore (Marquês de Pombal) yielded the best correlations between model and observations.

Lisbon↗

USGS2025

Unveiling coseismic deformation from differenced legacy aerial photography and modern lidar topography: The 1983 M6.9 Borah Peak earthquake, Idaho, USA

The 1983 M6.9 Borah Peak, Idaho, earthquake is one of the largest historical normal fault earthquakes in the western United States. We quantified meter-scale vertical change along the 35 km-long rupture using topographic differencing of 1966 aerial imagery and 2019 lidar-derived data. The initial differencing results are largely obscured by horizontal and vertical georeferencing errors and flight-line stripes. Our error corrections are designed to be insensitive to the coseismic deformation and reduced error by 50%. We calculated vertical separation and resolved a maximum of 2.02 ± 0.46 m at Doublespring Pass. Our vertical separation measurements are generally consistent with those from prior studies using field data and post-earthquake topographic data. However, the differencing measurements are a few decimeters lower than these prior measurements, indicating that differencing can isolate historical from prehistoric earthquake deformation. Our study demonstrates that revisiting historical earthquakes can provide new insights into the magnitude and patterns of coseismic deformation.

Idaho↗