USGS Science⌕ Search

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.

At least 1,747 records · Page 97Linked to original sources

Depositional history and fault-related studies, Bolinas Lagoon, California

Studies of core sediments and seismic reflection profiles elucidate the structure and depositional history of Bolinas Lagoon, Calif., which covers 4.4 km 2 and lies in the San Andreas fault zone at the southeast corner of the Point Reyes Peninsula 20 km northwest of San Francisco. The 1906 trace of the San Andreas fault crosses the west side of the lagoon and was determined from (1) tectonically caused salt-marsh destruction indicated by comparison of 1854 and 1929 U.S. Coast and Geodetic Survey (U.S.C. & G.S.) topographic surveys, (2) formation of a tidal channel along the border of destroyed salt marshes, and (3) azimuths of the trend of the fault measured in 1907. Subsidence in the lagoon of 30 cm occurred east of the San Andreas fault in 1906. Near the east shore, seismic-reflection profiling indicates the existence of a graben fault that may connect to a graben fault on the Golden Gate Platform. Comparison of radiocarbon dates on shells and plant debris from boreholes drilled on Stinson Beach spit with a relative sea-level curve constructed for southern San Francisco Bay indicates 5.8 to more than 17.9 m of tectonic subsidence of sediments now located 33 m below mean sea level. Cored sediments indicate a marine transgression dated at 7770?65 yrs B.P. overlying freshwater organic-rich lake deposits. Fossil pollen including 2 to 8 percent Picea (spruce) indicate a late Pleistocene (?)-Early Holocene climate, cooler, wetter, and foggier than at present. Above the transgression are discontinuous and interfingering sequences of transgressive-regressive marine, estuarine, and barrier sediments that reflect rapid lateral and vertical shifts of successive depositional environments. Fossil megafauna indicate (1) accumulation in a protected, shallow-water estuary or bay, and (2) that the lagoon was probably continuously shallow and never a deep-water embayment. Analysis of grain-size parameters, pollen frequencies, and organic remains from a core near the north end of the lagoon indicates (1) that mid-nineteenth-century redwood logging correlates with rates of sediment accumulation of l.3 to 1.9 cm/yr that are three to 6 times higher than post-1906 rates of 0.3 to 0.4 cm/yr, (2) accumulation of up to 115 cm of sediment since 1849, and (3) an anomalously coarse-grained sediment that may correlate with the 1906 earthquake.

Open-File Report↗

Chapter B. The Loma Prieta, California, Earthquake of October 17, 1989 - Highway Systems

This paper summarizes the impact of the Loma Prieta earthquake on highway systems. City streets, urban freeways, county roads, state routes, and the national highway system were all affected. There was damage to bridges, roads, tunnels, and other highway structures. The most serious damage occurred in the cities of San Francisco and Oakland, 60 miles from the fault rupture. The cost to repair and replace highways damaged by this earthquake was $2 billion. About half of this cost was to replace the Cypress Viaduct, a long, elevated double-deck expressway that had a devastating collapse which resulted in 42 deaths and 108 injuries. The earthquake also resulted in some positive changes for highway systems. Research on bridges and earthquakes began to be funded at a much higher level. Retrofit programs were started to upgrade the seismic performance of the nation's highways. The Loma Prieta earthquake changed earthquake policy and engineering practice for highway departments not only in California, but all over the world.

Professional Paper↗

Using seismic b-values to interpret seismicity rates and physical processes during the preeruptive earthquake swarm at Augustine Volcano 2005-2006

We use seismic b -values to explore physical processes during the Augustine Volcano 2005–6 preeruptive earthquake swarm. The preeruptive earthquake swarm was divided into two parts: the “long swarm,” which extended from April 30, 2005, to January 10, 2006; and the "short swarm," which started 13 hours before the onset of explosive activity on January 11, 2006. Calculations of b -value for each of these swarms and for a background period were performed. The short swarm, directly preceding the eruption, had the lowest calculated b -value. In addition to the low value, the shape of the b -value plot for the short swarm appears to have two separate slopes, a shallower slope for magnitudes as great as 1.2 and a steeper slope for magnitudes greater than 1.2. Calculations of b were also run for three precursory deformation stages suggested by a separate investigation of deformation at Augustine Volcano. The highest b -value, found in stage 2, may indicate an increase in pore pressure and in thermal gradient, which matches the geodetic interpretation of a proposed dike intrusion. Finer resolution changes of b are explored through calculations of b -value versus time. An initial drop in b -value in late 2004 preceded the onset of increased seismicity. The temporal nature of this change and its timing are corroborated by atmospheric temperature data recorded on the summit of the volcano, which increased at approximately the same time. Stress at Augustine Volcano was also studied using 79 earthquakes that returned acceptable focal mechanisms between January 1, 2002, and January 10, 2006. These mechanisms and an attempted stress-tensor inversion imply that stresses within the Augustine edifice are highly variable and do not display a dominant faulting style. A population of high-frequency volcano-tectonic earthquakes during the short swarm is found to have accompanying very-long-period (20 seconds and greater) energy. Statistical analysis indicates that these earthquakes are a separate population of events. We interpret this population of earthquakes to represent a separate and distinct physical process that was not seen before the 13 hours preceding the eruption. The b -value time series also indicates that when changes in stress, pore pressure, and thermal gradient occur simultaneously, that stress effects dominate the observed b -value.

Alaska↗

Earthquakes and tsunami

Earthquakes occur as a burst of sudden ground shaking created by the release of accumulated stress along a fault, often influenced by movement of the world’s tectonic plates. Ground shaking from an earthquake can generate additional hazards, including landslides, liquefaction, and tsunami. According to the 2019 “Global Assessment Report on Disaster Risk Reduction”, earthquakes combined with tsunami are the most damaging environmental hazards globally. Impacts of earthquakes and tsunami on people have increased around the world as human development of built infrastructure continues to expand. Adverse earthquake and tsunami impacts can be reduced through strategies including land-use planning, engineering, mitigation and preparedness, emergency planning, warnings, and exercises. The specific disaster risk reduction approaches taken will depend on the country, considering the geography, built environment, and social and cultural contexts. Wherever the location, it is important that such measures are considered and implemented holistically, as singular approaches may not be effective in addressing earthquake and tsunami challenges.

Book chapter↗

Preliminary analysis of the earthquake (MW 8.1) and tsunami of April 1, 2007, in the Solomon Islands, southwestern Pacific Ocean

On April 1, 2007, a destructive earthquake (Mw 8.1) and tsunami struck the central Solomon Islands arc in the southwestern Pacific Ocean. The earthquake had a thrust-fault focal mechanism and occurred at shallow depth (between 15 km and 25 km) beneath the island arc. The combined effects of the earthquake and tsunami caused dozens of fatalities and thousands remain without shelter. We present a preliminary analysis of the Mw-8.1 earthquake and resulting tsunami. Multichannel seismic-reflection data collected during 1984 show the geologic structure of the arc's frontal prism within the earthquake's rupture zone. Modeling tsunami-wave propagation indicates that some of the islands are so close to the earthquake epicenter that they were hard hit by tsunami waves as soon as 5 min. after shaking began, allowing people scant time to react.

Science of Tsunami Hazards↗

Photointerpretation of Alaskan post-earthquake photography

Aerial photographs taken after the March 27, 1964, Good Friday, Alaskan earthquake were examined stereoscopically to determine effects of the earthquake in areas remote from the towns, highways, and the railroad. The two thousand black and white photographs used in this study were taking in April, after the earthquake, by the U. S. Coast & Geodetic Survey and were generously supplied to the U. S. Geological Survey. Part of the photographs, at a scale of 1/24,000, provide blanket coverage of approximately 2,000 square miles of land area north and west of Prince William Sound, including parts of the mainland and some of the adjacent islands. The epicenter of the earthquake, near the head of Unakwik Inlet, is located in this area. The rest of the photographs, at scales ranging from 1/17,000 to 1/40,000, cover isolated strips of the coastline of the mainland and nearby islands in the general area of Prince William Sound. Figure 1 shows the area of new photo coverage used in this study. The objective of the study was to determine quickly whether geological features resulting from the earthquake, such as faults, changes in shoreline, cracks in surficial material, pressure ridges in lake ice, fractures in glaciers and lake ice, and rock slides and avalanches, might be identifiable by photointerpretation. The study was made without benefit of comparisons with older, or pre-earthquake photography, which was not readily available for immediate use.

Alaska↗

Morphology and age of fault scarps in the Rio Grande Rift, south-central Colorado

Fault scarps in the Rio Grande rift of Colorado provide most of the evidence of paleoseismicity in the state, and are thus a major focus of assessments of earthquake hazards. Critical components of such assessments are the ages of past faulting events; age control is scarce and generally coarse for the Rio Grande rift in Colorado. Ages estimated from fault-scarp morphology are thus valuable, partly because they are directly related to the fault eVents, and partly for want of alternative methods. The processes responsible for scarp degradation are poorly understood and the many variables that affect the rates and relative importance of those processes limits the precision and accuracy of age estimates derived from scarp morphology. Both analytical and empirical methods are available for relating scarp degradation processes to scarp morphology and age. Diffusion-equation models that have been derived for scarp degradation are useful for their predictive value and for the insight they provide about the degradation processes. However, because of the limitations mentioned above and because of the multiple-event nature of the scarps in the Rio Grande rift, the simpler empirical approach involving the relation between scarp angle and the logarithm of scarp height is used here. The scarp-morphology data from the Rio Grande rift suggest that the youngest movements on all of the scarps examined occurred less than 15,000 years ago, with the exception of two scarps that are formed on lithologically dissimilar deposits. This age estimate is consistent with the observation that all of the scarps occur at least in part on deposits of Pinedale age. The scarp-morphology data suggest that some of the scarps are closer to 5,000 years than to 15,000 years in age; several of these scarps occur on deposits of early Holocene age and are associated with early Holocene radiocarbon dates. One scarp occurs on deposits with radiocarbon dates of only a few thousand years, although the morphology of this scarp is unusual for its age. Three scarps have morphologies that suggest that they may represent fault events younger than 5,000 years old. Although no corroboration is available, these scarps may represent some of the youngest fault movements in Colorado.

Colorado↗

Stochastic simulation of earthquake ground motions for the seismic assessment of monumental masonry structures: Source-based vs site-based approaches

Earthquakes are among the most destructive natural disasters and have resulted in a massive number of fatalities and economic losses all over the world. Simulated ground motion records are valuable, particularly for regions lacking seismic stations or with a limited history of large-magnitude earthquakes. Notably, a significant percentage of monumental masonry buildings are located in regions with limited access to real records; hence, simulated records play a paramount role in their seismic protection. However, few studies have investigated the structural response of heritage buildings via response history analyses to assess the performance of simulated earthquakes against real ones. To accomplish this, this study simulates the recorded time-series of the 9th of July 1998 Faial earthquake in the Azores (M w = 6.2) at four available stations, using two different simulation approaches, that is, a source-based stochastic finite-fault method and a site-based broadband stochastic method. First, two masonry facades with sidewalls characterized by different slenderness levels are adopted to conduct this research. Moreover, the proposed approach is also applied to an existing monumental structure, that is, São Francisco Church, located at Horta, which was affected by damage during the Faial earthquake. Results demonstrate that both simulation approaches provide similar results in terms of structural response prediction. The proposed framework also demonstrates that a small mismatch in terms of predicted damage patterns can result in a significant relative error in terms of displacement predictions.

Earthquake Engineering and Structural Dynamics↗

The October 1992 Parkfield, California, earthquake prediction

A magnitude 4.7 earthquake occurred near Parkfield, California, on October 20, 992, at 05:28 UTC (October 19 at 10:28 p.m. local or Pacific Daylight Time).This moderate shock, interpreted as the potential foreshock of a damaging earthquake on the San Andreas fault, triggered long-standing federal, state and local government plans to issue a public warning of an imminent magnitude 6 earthquake near Parkfield. Although the predicted earthquake did not take place, sophisticated suites of instruments deployed as part of the Parkfield Earthquake Prediction Experiment recorded valuable data associated with an unusual series of events. this article describes the geological aspects of these events, which occurred near Parkfield in October 1992. The accompnaying article, an edited version of a press conference b Richard Andrews, the Director of the California Office of Emergency Service (OES), describes governmental response to the prediction.

California↗

Integrated geophysical characteristics of the 2015 Illapel, Chile, earthquake

On 16 September 2015, a Mw 8.3 earthquake ruptured the subduction zone offshore of Illapel, Chile, generating an aftershock sequence with 14 Mw 6.0–7.0 events. A double source W phase moment tensor inversion consists of a Mw 7.2 subevent and the main Mw 8.2 phase. We determine two slip models for the mainshock, one using teleseismic broadband waveforms and the other using static GPS and InSAR surface displacements, which indicate high slip north of the epicenter and west-northwest of the epicenter near the oceanic trench. These models and slip distributions published in other studies suggest spatial slip uncertainties of ~25 km and have peak slip values that vary by a factor of 2. We relocate aftershock hypocenters using a Bayesian multiple-event relocation algorithm, revealing a cluster of aftershocks under the Chilean coast associated with deep (20–45 km depth) mainshock slip. Less vigorous aftershock activity also occurred near the trench and along strike of the main aftershock region. Most aftershocks are thrust-faulting events, except for normal-faulting events near the trench. Coulomb failure stress change amplitudes and signs are uncertain for aftershocks collocated with deeper mainshock slip; other aftershocks are more clearly associated with loading from the mainshock. These observations reveal a frictionally heterogeneous interface that ruptured in patches at seismogenic depths (associated with many aftershocks) and with homogeneous slip (and few aftershocks) up to the trench. This event likely triggered seismicity separate from the main slip region, including along-strike events on the megathrust and intraplate extensional events.

Illapel↗

Strain pattern represented by scarps formed during the earthquakes of October 2, 1915, Pleasant Valley, Nevada

The pattern of scarps developed during the earthquakes of October 2, 1915, in Pleasant Valley, Nevada, may have formed as a result of a modern stress system acting on a set of fractures produced by an earlier stress system which was oriented differently. Four major scarps developed in a right-stepping, en-echelon pattern suggestive of left-lateral slip across the zone and an extension axis oriented approximately S85°W. The trend of the zone is N25°E. However, the orientation of simple dip-slip on most segments trending approximately N20—40° E and a right-lateral component of displacement on several N- and NW-trending segments of the scarps indicate that the axis of regional extension was oriented between N50° and 70° W, normal to the zone. The cumulative length of the scarps is 60 km, average vertical displacement 2 m, and the maximum vertical displacement near the Pearce School site 5.8 m. Almost everywhere the 1915 scarps formed along an older scarp line, and in some places older scarps represent multiple previous events. The most recent displacement event prior to 1915 is interpreted to have occurred more than 6600 years ago, but possibly less than 20,000 years ago. Some faults expressed by older scarps that trend northwest were not reactivated in 1915, possibly because they are oriented at a low angle with respect to the axis of modern regional extension. The 1915 event occurred in an area of overlap of three regional fault trends oriented northwest, north, and northeast and referred to, respectively, as the Oregon—Nevada, Northwest Nevada, and Midas—Battle Moutain trends. Each of these trends may have developed at a different time; the Oregon—Nevada trend was possibly the earliest and developed in Late Miocene time (Stewart et al. 1975). Segments of the 1915 scarps ar

Nevada↗

36 - Implications of crustal strain during conventional, slow, and silent earthquakes

Uniform block-slip motion consistent with simple shear on locked fault segments is the primary feature apparent in geodetic measurements of strain accumulation along plate boundaries (Savage, 1983). However, almost every aspect of fault failure is nonlinear in character. This premise derives from theoretical models (Kostrov, 1966; Richards, 1976; Andrews, 1976; Freund, 1979; Rice and Rudnicki, 1979; Rice, 1983, 1992; Stuart, 1979, Stuart and Mavko, 1979; Das and Scholz, 1981; Rundle et al., 1984; Tse and Rice, 1986; Lorenzetti and Tullis, 1989; Segall and Rice, 1995; Shaw, 1997; Miller et al., 1999) and laboratory-generated frictional failure of crustal materials (Dieterich, 1979, 1981; Mogi, 1981; Mogi et al., 1982) which predict accelerating deformation will occur before dynamic slip instabilities, better known as earthquakes.

International Geophysics↗

On the modified Mercalli intensities and magnitudes of the 1811-1812 New Madrid earthquakes

We reexamine original felt reports from the 1811–1812 New Madrid earthquakes and determine revised isoseismal maps for the three principal mainshocks. In many cases we interpret lower values than those assigned by earlier studies. In some cases the revisions result from an interpretation of original felt reports with an appreciation for site response issues. Additionally, earlier studies had assigned modified Mercalli intensity (MMI) values of V-VII to a substantial number of reports that we conclude do not describe damage commensurate with intensities this high. We investigate several approaches to contouring the MMI values using both analytical and subjective methods. For the first mainshock on 02∶15 LT December 16, 1811, our preferred contouring yields M ω7.2–7.3 using the area-moment regressions of Johnston [1996]. For the 08∶00 LT on January 23, 1812, and 03∶45 LT on February 7, 1812, mainshocks, we obtain M ω7.0 and M ω7.4–7.5, respectively. Our magnitude for the February mainshock is consistent with the established geometry of the Reelfoot fault, which all evidence suggests to have been the causative structure for this event. We note that the inference of lower magnitudes for the New Madrid events implies that site response plays a significant role in controlling seismic hazard at alluvial sites in the central and eastern United States. We also note that our results suggest that thrusting may have been the dominant mechanism of faulting associated with the 1811–1812 sequence.

Arkansas, Illinois, Kentucky, Missouri, Tennessee↗

Ground deformation associated with the March 1996 earthquake swarm at Akutan volcano, Alaska, revealed by satellite radar interferometry

In March 1996 an intense swarm of volcano-tectonic earthquakes (???3000 felt by local residents, Mmax = 5.1, cumulative moment of 2.7 ??1018 N m) beneath Akutan Island in the Aleutian volcanic arc, Alaska, produced extensive ground cracks but no eruption of Akutan volcano. Synthetic aperture radar interferograms that span the time of the swarm reveal complex island-wide deformation: the western part of the island including Akutan volcano moved upward, while the eastern part moved downward. The axis of the deformation approximately aligns with new ground cracks on the western part of the island and with Holocene normal faults that were reactivated during the swarm on the eastern part of the island. The axis is also roughly parallel to the direction of greatest compressional stress in the region. No ground movements greater than 2.83 cm were observed outside the volcano's summit caldera for periods of 4 years before or 2 years after the swarm. We modeled the deformation primarily as the emplacement of a shallow, east-west trending, north dipping dike plus inflation of a deep, Mogi-type magma body beneath the volcano. The pattern of subsidence on the eastern part of the island is poorly constrained. It might have been produced by extensional tectonic strain that both reactivated preexisting faults on the eastern part of the island and facilitated magma movement beneath the western part. Alternatively, magma intrusion beneath the volcano might have been the cause of extension and subsidence in the eastern part of the island. We attribute localized subsidence in an area of active fumaroles within the Akutan caldera, by as much as 10 cm during 1992-1993 and 1996-1998, to fluid withdrawal or depressurization of the shallow hydrothermal system. Copyright 2000 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Response of seismicity to Coulomb stress triggers and shadows of the 1999 Mw=7.6 Chi-Chi, Taiwan, earthquake

The correlation between static Coulomb stress increases and aftershocks has thus far provided the strongest evidence that stress changes promote seismicity, a correlation that the Chi-Chi earthquake well exhibits. Several studies have deepened the argument by resolving stress changes on aftershock focal mechanisms, which removes the assumption that the aftershocks are optimally oriented for failure. Here one compares the percentage of planes on which failure is promoted after the main shock relative to the percentage beforehand. For Chi-Chi we find a 28% increase for thrust and an 18% increase for strike-slip mechanisms, commensurate with increases reported for other large main shocks. However, perhaps the chief criticism of static stress triggering is the difficulty in observing predicted seismicity rate decreases in the stress shadows, or sites of Coulomb stress decrease. Detection of sustained drops in seismicity rate demands a long catalog with a low magnitude of completeness and a high seismicity rate, conditions that are met at Chi-Chi. We find four lobes with statistically significant seismicity rate declines of 40-90% for 50 months, and they coincide with the stress shadows calculated for strike-slip faults, the dominant faulting mechanism. The rate drops are evident in uniform cell calculations, 100-month time series, and by visual inspection of the M ??? 3 seismicity. An additional reason why detection of such declines has proven so rare emerges from this study: there is a widespread increase in seismicity rate during the first 3 months after Chi-Chi, and perhaps many other main shocks, that might be associated with a different mechanism. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Postearthquake relaxation after the 2004 M6 Parkfield, California, earthquake and rate-and-state friction

An unusually complete set of measurements (including rapid rate GPS over the first 10 days) of postseismic deformation is available at 12 continuous GPS stations located close to the epicenter of the 2004 M6.0 Parkfield earthquake. The principal component modes for the relaxation of the ensemble of those 12 GPS stations were determined. The first mode alone furnishes an adequate approximation to the data. Thus, the relaxation at all stations can be represented by the product of a common temporal function and distinct amplitudes for each component (north or east) of relaxation at each station. The distribution in space of the amplitudes indicates that the relaxation is dominantly strike slip. The temporal function, which spans times from about 5 min to 900 days postearthquake, can be fit by a superposition of three creep terms, each of the form ??l loge(1 + t/??l), with characteristic times ??, = 4.06, 0.11, and 0.0001 days. It seems likely that what is actually involved is a broad spectrum of characteristic times, the individual components of which arise from afterslip on different fault patches. Perfettini and Avouac (2004) have shown that an individual creep term can be explained by the spring-slider model with rate-dependent (no state variable) friction. The observed temporal function can also be explained using a single spring-slider model (i.e., single fault patch) that includes rate-and-state-dependent friction, a single-state variable, and either of the two commonly used (aging and slip) state evolution laws. In the latter fits, the rate-and-state friction parameter b is negative.

Journal of Geophysical Research B: Solid Earth↗

Aftershock seismicity of the 2010 Maule Mw=8.8 Chile, earthquake: Correlation between co-seismic slip models and aftershock distribution?

The 27 February 2010 Maule, Chile (Mw=8.8) earthquake is one of the best instrumentally observed subduction zone megathrust events. Here we present locations, magnitudes and cumulative equivalent moment of the first -2 months of aftershocks, recorded on a temporary network deployed within 2 weeks of the occurrence of the mainshock. Using automatically-determined onset times and a back projection approach for event association, we are able to detect over 30,000 events in the time period analyzed. To further increase the location accuracy, we systematically searched for potential S-wave arrivals and events were located in a regional 2D velocity model. Additionally, we calculated regional moment tensors to gain insight into the deformation history of the aftershock sequence. We find that the aftershock seismicity is concentrated between 40 and 140 km distance from the trench over a depth range of 10 to 35 km. Focal mechanisms indicate a predominance of thrust faulting, with occasional normal faulting events. Increased activity is seen in the outer-rise region of the Nazca plate, predominantly in the northern part of the rupture area. Further down-dip, a second band of clustered seismicity, showing mainly thrust motion, is located at depths of 40–45 km. By comparing recent published mainshock source inversions with our aftershock distribution, we discriminate slip models based on the assumption that aftershocks occur in areas of rapid transition between high and low slip, surrounding high-slip regions of the mainshock.

Geophysical Research Letters↗

Rupture process of the Mw7.0 December 5, 2024 Offshore Cape Mendocino earthquake

The Mw7.0 December 5, 2024 Offshore Cape Mendocino earthquake ruptured a km long portion of the east-west trending Mendocino fault zone (MFZ). In order to clarify the rupture process, we assemble three-component seismograms from regional seismic stations, horizontal coseismic displacement vectors derived from Global Navigation Satellite System (GNSS) time series, and a Sentinel-1 ascending interferogram. These data are interpreted with a model of slip distributed on two vertical fault planes representative of the eastern MFZ and spanning the ~70 km length of the aftershock zone. Assuming right-lateral strike slip, we find that the rupture initiates in the oceanic mantle at 20-30 km depth and proceeds unilaterally updip and toward the east. Early aftershocks locate adjacent to the peak slip areas, tracking the coseismic rupture propagation from oceanic mantle to shallower depth and implying a significant role of static stress transfer in driving aftershocks in an ocean plate environment.

California↗