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

USGS2007

Investigation of the M6.6 Niigata-Chuetsu Oki, Japan, earthquake of July 16, 2007

The M6.6 mainshock of the Niigata Chuetsu Oki (offshore) earthquake occurred at 10:13 a.m. local time on July 16, 2007, and was followed by a sequence of aftershocks that were felt during the entire time of the reconnaissance effort. The mainshock had an estimated focal depth of 10 km and struck in the Japan Sea offshore Kariwa. Analysis of waveforms from source inversion studies indicates that the event occurred along a thrust fault with a NE trend. The fault plane is either a strike of 34 degrees with a dip of 51 degrees or a strike of 238 degrees with a dip of 41 degrees. Which of these two planes is associated with the mainshock rupture is unresolved, although attenuation relationship analysis indicates that the northwest-dipping fault is favored. The quake affected an approximately 100-km-wide area along the coastal areas of southwestern Niigata prefecture. The event triggered ground failures as far as the Unouma Hills, located in central Niigata approximately 50 km from the shore and the source area of the 2004 Niigata Chuetsu earthquake. The primary event produced tsunami run-ups that reached maximum runup heights of about 20 centimeters along the shoreline of southern Niigata Prrefecture.

Niigata↗

USGS2008

Potential effects of a scenario earthquake on the economy of southern California: Intraregional commuter, worker, and earnings flow analysis

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 (Jones and others, 2008). This report uses selected datasets from the U.S. Census Bureau and the State of California's Employment Development Department to develop preliminary estimates of the number and spatial distribution of commuters who cross the San Andreas Fault and to characterize these commuters by the industries in which they work and their total earnings. The analysis concerns the relative exposure of the region's economy to the effects of the earthquake as described by the location, volume, and earnings of those commuters who work in each of the region's economic sectors. It is anticipated that damage to transportation corridors traversing the fault would lead to at least short-term disruptions in the ability of commuters to travel between their places of residence and work.

California↗

USGS2014

Key recovery factors for the August 24, 2014, South Napa Earthquake

Through discussions between the Federal Emergency Management Agency (FEMA) and the U.S. Geological Survey (USGS) following the South Napa earthquake, it was determined that several key decision points would be faced by FEMA for which additional information should be sought and provided by USGS and its partners. This report addresses the four tasks that were agreed to. These tasks are (1) assessment of ongoing fault movement (called afterslip) especially in the Browns Valley residential neighborhood, (2) assessment of the shaking pattern in the downtown area of the City of Napa, (3) improvement of information on the fault hazards posed by the West Napa Fault System (record of past earthquakes and slip rate, for example), and (4) imagery acquisition and data processing to provide overall geospatial information support to FEMA.

California↗

USGS1974

Strain release mechanism of the 1906 San Francisco earthquake

Reexamination of geodetic data has shown that aseismic slip occurred on or near the San Andreas fault in the period of about 20 years after 1906. The inferred displacements are comparable to but at greater depths than the sudden slip that occurred at the time of the earthquake. The postseismic slip is constrained only between late 1906 and 1925, and data are insufficient to determine the movements, if any, below about 20 kilometers on the fault. Two independent observations also indicate substantial anomalous crustal deformation away from the fault at least 30 years before the earthquake.

Science↗

USGS1992

The San Andreas Fault System, California, USA

Geologists, seismologists, and geophysicists have intensively studied the San Andreas fault system for the past 20 to 30 years. Their goals were to learn more about damaging earthquakes, the behavior of major stirke-slip faults, and methods of reducing earthquake hazards in populated areas. Field geologic investigations, seismic networks, post-earthquake studies, precision geodetic surveys, and reflection and refraction seismic surveys are among the methods used to decipher the history, geometry, and mechanics of the system. -from Authors

Annales Tectonicae↗

USGS2004

Mitigation of earthquake damage

The article describes the use of a geologic map to help mitigate earthquake damage along the Denali Fault where the Trans-Alaska Pipeline crosses. Geologic mapping of bedrock and unconsolidated deposits reveals a history of horizontal right-lateral slip and local vertical separations at the fault. It was determined that the eastern 220 mile of the Denali and Totschunda fault system was the most likely segment to generate an 8+ magnitude earthquake.

Alaska↗

USGS2013

A fluid-driven earthquake swarm on the margin of the Yellowstone caldera

Over the past several decades, the Yellowstone caldera has experienced frequent earthquake swarms and repeated cycles of uplift and subsidence, reflecting dynamic volcanic and tectonic processes. Here, we examine the detailed spatial-temporal evolution of the 2010 Madison Plateau swarm, which occurred near the northwest boundary of the Yellowstone caldera. To fully explore the evolution of the swarm, we integrated procedures for seismic waveform-based earthquake detection with precise double-difference relative relocation. Using cross-correlation of continuous seismic data and waveform templates constructed from cataloged events, we detected and precisely located 8710 earthquakes during the three-week swarm, nearly four times the number of events included in the standard catalog. This high-resolution analysis reveals distinct migration of earthquake activity over the course of the swarm. The swarm initiated abruptly on January 17, 2010 at about 10 km depth and expanded dramatically outward (both shallower and deeper) over time, primarily along a NNW-striking, ~55º ENE-dipping structure. To explain these characteristics, we hypothesize that the swarm was triggered by the rupture of a zone of confined high-pressure aqueous fluids into a pre-existing crustal fault system, prompting release of accumulated stress. The high-pressure fluid injection may have been accommodated by hybrid shear and dilatational failure, as is commonly observed in exhumed hydrothermally affected fault zones. This process has likely occurred repeatedly in Yellowstone as aqueous fluids exsolved from magma migrate into the brittle crust, and it may be a key element in the observed cycles of caldera uplift and subsidence.

Wyoming↗

USGS1985

Prediction of central California earthquakes from soil-gas helium fluctuations

The observations of short-term decreases in helium soil-gas concentrations along the San Andreas Fault in central California have been correlated with subsequent earthquake activity. The area of study is elliptical in shape with radii approximately 160??80 km, centered near San Benito, and with the major axis parallel to the Fault. For 83 percent of the M>4 earthquakes in this area a helium decrease preceded seismic activity by 1.5 to 6.5 weeks. There were several earthquakes without a decrease and several decreases without a corresponding earthquake. Owing to complex and unresolved interaction of many geophysical and geochemical parameters, no suitable model is yet developed to explain the observations. ?? 1985 Birkha??user Verlag.

Pure and Applied Geophysics PAGEOPH↗

USGS2013

What do data used to develop ground-motion prediction equations tell us about motions near faults?

A large database of ground motions from shallow earthquakes occurring in active tectonic regions around the world, recently developed in the Pacific Earthquake Engineering Center’s NGA-West2 project, has been used to investigate what such a database can say about the properties and processes of crustal fault zones. There are a relatively small number of near-rupture records, implying that few recordings in the database are within crustal fault zones, but the records that do exist emphasize the complexity of ground-motion amplitudes and polarization close to individual faults. On average over the whole data set, however, the scaling of ground motions with magnitude at a fixed distance, and the distance dependence of the ground motions, seem to be largely consistent with simple seismological models of source scaling, path propagation effects, and local site amplification. The data show that ground motions close to large faults, as measured by elastic response spectra, tend to saturate and become essentially constant for short periods. This saturation seems to be primarily a geometrical effect, due to the increasing size of the rupture surface with magnitude, and not due to a breakdown in self similarity.

Pure and Applied Geophysics↗

USGS2017

Performance of Irikura recipe rupture model generator in earthquake ground motion simulations with Graves and Pitarka hybrid approach

We analyzed the performance of the Irikura and Miyake (Pure and Applied Geophysics 168(2011):85–104, 2011) (IM2011) asperity-based kinematic rupture model generator, as implemented in the hybrid broadband ground motion simulation methodology of Graves and Pitarka (Bulletin of the Seismological Society of America 100(5A):2095–2123, 2010), for simulating ground motion from crustal earthquakes of intermediate size. The primary objective of our study is to investigate the transportability of IM2011 into the framework used by the Southern California Earthquake Center broadband simulation platform. In our analysis, we performed broadband (0–20 Hz) ground motion simulations for a suite of M6.7 crustal scenario earthquakes in a hard rock seismic velocity structure using rupture models produced with both IM2011 and the rupture generation method of Graves and Pitarka (Bulletin of the Seismological Society of America, 2016) (GP2016). The level of simulated ground motions for the two approaches compare favorably with median estimates obtained from the 2014 Next Generation Attenuation-West2 Project (NGA-West2) ground motion prediction equations (GMPEs) over the frequency band 0.1–10 Hz and for distances out to 22 km from the fault. We also found that, compared to GP2016, IM2011 generates ground motion with larger variability, particularly at near-fault distances (<12 km) and at long periods (>1 s). For this specific scenario, the largest systematic difference in ground motion level for the two approaches occurs in the period band 1–3 s where the IM2011 motions are about 20–30% lower than those for GP2016. We found that increasing the rupture speed by 20% on the asperities in IM2011 produced ground motions in the 1–3 s bandwidth that are in much closer agreement with the GMPE medians and similar to those obtained with GP2016. The potential implications of this modification for other rupture mechanisms and magnitudes are not yet fully understood, and this topic is the subject of ongoing study. We concluded that IM2011 rupture generator performs well in ground motion simulations using Graves and Pitarka hybrid method. Therefore, we recommend it to be considered for inclusion into the framework used by the Southern California Earthquake Center broadband simulation platform.

Pure and Applied Geophysics↗

USGS1975

Vertical crustal movements preceding and accompanying the San Fernando earthquake of February 9, 1971: A summary

Comparative elevations referred chiefly to a tidal bench mark with a history of relatively positive movement show that much of the Transverse Ranges of southern California sustained major changes in elevation both before and in association with the M L 6.4 San Fernando earthquake of February 9, 1971. Preseismic changes in elevation in the southern Transverse Ranges were almost uniformly positive and generally episodic. Maximum uplift measured between 1960/61 and 1968/69 was 0.205 m (observed) or 0.200 m (adjusted); it occurred about 30 km northeast of the 1971 epicenter. As much as 0.129 m (observed) of preseismic uplift was also measured about 30 km northwest of the epicenter between 1964 and 1968. A comparison between 1968 and 1969 elevation data revealed 0.078 m (observed) of differential uplift centered about 10 km west of the epicenter; earlier leveling indicates that this uplift began between February 1967 and May 1968. Changes in elevation measured during the interval 1968/69-71 (postearthquake) are interpreted as chiefly coseismic and were concentrated along a 15-km ruptured segment of the east-trending, north-dipping San Fernando fault. These changes were characterized by subsidence of as much as 0.111 m (observed) or 0.092 m (adjusted) south of the fault and a ridge of uplift of up to at least 2.195 m (observed) or 2.196 m (adjusted) immediately north of the fault. The more episodic preseismic movements are interpretable as deep-seated creep events on the San Fernando fault. Alternatively, these episodic movements may be due partly to the operation of dilatancy; the onset of the 196849 epicentral uplift accords closely with the initiation of a V p / V s anomaly recognized in this same area. The 1968/69-71 elevation changes can be attributed almost entirely to slip on the San Fernando fault.

California↗

USGS2015

Simulation of acceleration field of the Lushan earthquake (Ms7.0, April 20, 2013, China)

The acceleration field of the Lushan earthquake (Ms7.0, April 20, 2013, China) is simulated using a new modified version of the stochastic finite-fault method (EXSIM) based on a dynamic corner frequency approach. To incorporate the effect of heterogeneous slip distribution on the variation of source spectrum, we adopt an empirical source spectral model and derive the corresponding dynamic parameters, which vary with the cumulative seismic moment of the ruptured area. The new modified method is validated by: 1) comparison of the simulation results with those obtained from the EXSIM method using near-fault ground motion data of the 1994 Northridge earthquake; 2) comparison of simulated PGA contour map inferred from synthetic time histories at 315 grid locations with the observed PGA shakemap for the 2013 Lushan earthquake; 3) comparison of simulated PGA with those predicted by ground-motion prediction equations (GMPEs); and 4) comparison of simulated time histories with observed acceleration records at six strong motion stations during the mainshock of the Lushan earthquake, in which local site response is considered in the simulation. These comparisons confirm the validity of the new simulation procedure for purposes of regional strong ground motion estimation. Limitations of the procedure in modeling the phasing of different arrivals in the seismic signal and near-surface response of geologic deposits are discussed.

Lushan↗

USGS2014

Ball-and-socket tectonic rotation during the 2013 M w 7.7 Balochistan earthquake

The September 2013 M w 7.7 Balochistan earthquake ruptured a ∼200-km-long segment of the curved Hoshab fault in southern Pakistan with 10±0.2 m of peak sinistral and ∼1.7±0.8 m of dip slip. This rupture is unusual because the fault dips 60±15° towards the focus of a small circle centered in northwest Pakistan, and, despite a 30° increase in obliquity along strike, the ratios of strike and dip slip remain relatively uniform. Surface displacements and geodetic and teleseismic source inversions quantify a bilateral rupture that propagated rapidly at shallow depths from a transtensional jog near the northern end of the rupture. Static friction prior to rupture was unusually weak (μ<0.05), and friction may have approached zero during dynamic rupture. Here we show that the inward-dipping Hoshab fault defines the northern rim of a structural unit in southeast Makran that rotates – akin to a 2-D ball-and-socket joint – counter-clockwise in response to India's penetration into the Eurasian plate. This rotation accounts for complexity in the Chaman fault system and, in principle, reduces seismic potential near Karachi; nonetheless, these findings highlight deficiencies in strong ground motion equations and tectonic models that invoke Anderson–Byerlee faulting predictions.

Balochistan↗

USGS2001

Source parameters for the 1952 Kern County earthquake, California: A joint inversion of leveling and triangulation observations

Coseismic leveling and triangulation observations are used to determine the faulting geometry and slip distribution of the July 21, 1952, M w 7.3 Kern County earthquake on the White Wolf fault. A singular value decomposition inversion is used to assess the ability of the geodetic network to resolve slip along a multisegment fault and shows that the network is sufficient to resolve slip along the surface rupture to a depth of 10 km. Below 10 km, the network can only resolve dip slip near the fault ends. The preferred source model is a two-segment right-stepping fault with a strike of 51° and a dip of 75° SW. The epicentral patch has deep (6–27 km) left-lateral oblique slip, while the northeastern patch has shallow (1–12.5 km) reverse slip. There is nearly uniform reverse slip (epicentral, 1.6 m; northeast, 1.9 m), with 3.6 m of left-lateral strike slip limited to the epicentral patch. The seismic moment is M 0 = 9.2 ± 0.5 × 10 19 N m ( M w = 7.2). The signal-to-noise ratio of the leveling and triangulation data is reduced by 96% and 49%, respectively. The slip distribution from the preferred model matches regional geomorphic features and may provide a driving mechanism for regional shortening across the Comanche thrust and structural continuity with the Scodie seismic lineament to the northeast.

California↗

USGS2004

Kinematic and dynamic rupture models of the November 3, 2002 Mw7.9 Denali, Alaska, earthquake

Regional seismic waveforms, continuous and campaign-mode GPS data, and surface slip measurements were used to obtain a kinematic model of the rupture process of the November 3, 2002 Mw 7.9 Denali, Alaska, earthquake. The event initiated as a Mw 7.0 reverse slip event on the north-dipping Susitna Glacier fault with subsequent right-lateral slip distributed over approximately 300 km of the Denali fault system. Near-shear rupture velocity is inferred from the kinematic modeling. The average and maximum slips were found to be 2.14 in and 10.3 m. Static stress drop varies from 1.3 to 5.0 MPa over the 5-segment fault model. Dynamic modeling shows the rupture propagated along the Susitna Glacier and Denali faults, then transferred to the Totschunda fault before stopping, largely due to the Totschunda's more favorable orientation with respect to the regional stress field.

Alaska↗

USGS2007

Stress orientations of Taiwan Chelungpu-Fault Drilling Project (TCDP) hole-A as observed from geophysical logs

[1] The Taiwan Chelungpu-fault Drilling Project (TCDP) drilled a 2-km-deep research borehole to investigate the structure and mechanics of the Chelungpu Fault that ruptured in the 1999 M w 7.6 Chi-Chi earthquake. Geophysical logs of the TCDP were carried out over depths of 500–1900 m, including Dipole Sonic Imager (DSI) logs and Formation Micro Imager (FMI) logs in order to identify bedding planes, fractures and shear zones. From the continuous core obtained from the borehole, a shear zone at a depth of 1110 meters is interpreted to be the Chelungpu fault, located within the Chinshui Shale, which extends from 1013 to 1300 meters depth. Stress-induced borehole breakouts were observed over nearly the entire length of the wellbore. These data show an overall stress direction (∼N115°E) that is essentially parallel to the regional stress field and parallel to the convergence direction of the Philippine Sea plate with respect to the Eurasian plate. Variability in the average stress direction is seen at various depths. In particular there is a major stress orientation anomaly in the vicinity of the Chelungpu fault. Abrupt stress rotations at depths of 1000 m and 1310 m are close to the Chinshui Shale's upper and lower boundaries, suggesting the possibility that bedding plane slip occurred during the Chi-Chi earthquake.

Geophysical Research Letters↗

USGS2009

Tsunamis and splay fault dynamics

The geometry of a fault system can have significant effects on tsunami generation, but most tsunami models to date have not investigated the dynamic processes that determine which path rupture will take in a complex fault system. To gain insight into this problem, we use the 3D finite element method to model the dynamics of a plate boundary/splay fault system. We use the resulting ground deformation as a time-dependent boundary condition for a 2D shallow-water hydrodynamic tsunami calculation. We find that if me stress distribution is homogeneous, rupture remains on the plate boundary thrust. When a barrier is introduced along the strike of the plate boundary thrust, rupture propagates to the splay faults, and produces a significantly larger tsunami man in the homogeneous case. The results have implications for the dynamics of megathrust earthquakes, and also suggest mat dynamic earthquake modeling may be a useful tool in tsunami researcn. Copyright 2009 by the American Geophysical Union.

Geophysical Research Letters↗

USGS2010

Postseismic relaxation following the 1994 M w 6.7 Northridge earthquake, southern California

We have reexamined the postearthquake deformation of a 65 km long linear array of 11 geodetic monuments extending north–south across the rupture (reverse slip on a blind thrust dipping 40°S–20°W) associated with the 1994 M w 6.7 Northridge earthquake. That array was surveyed frequently in the interval from 4 to 2650 days after the earthquake. The velocity of each of the monuments over the interval 100–2650 days postearthquake appears to be constant. Moreover, the profile of those velocities along the length of the array is very similar to a preearthquake velocity profile for a nearby, similarly oriented array. We take this to indicate that significant postseismic relaxation is evident only in the first 100 days postseismic and that the subsequent linear trend is typical of the interseismic interval. The postseismic relaxation (postseismic displacement less displacement that would have occurred at the preseismic velocity) is found to be almost wholly parallel (N70°W) to the nearby (40 km) San Andreas Fault with only negligible relaxation in the direction of coseismic slip (N20°E) on the Northridge rupture. We suggest that the N70°W relaxation is caused by aseismic, right-lateral slip at depth on the San Andreas Fault, excess slip presumably triggered by the Northridge rupture. Finally, using the Dieterich (1994) stress-seismicity relation, we show that return to the preseismic deformation rate within 100 days following the earthquake could be consistent with the cumulative number of M > 2.5 earthquakes observed following the main shock.

California↗