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

Late Pleistocene to Holocene slip rates for the Gurvan Bulag thrust fault (Gobi-Altay, Mongolia) estimated with 10Be dates

We surveyed morphotectonic markers along the central part of the Gurvan Bulag thrust, a fault that ruptured with the Bogd fault during the Gobi-Altay earthquake (1957, M 8.3), to document climatic and tectonic processes along the fault for the late Pleistocene- Holocene period. The markers were dated using 10Be produced in situ. Two major periods of alluviation ended at 131 ?? 20 and 16 ?? 4.8 ka. These appear to be contemporaneous with global climatic changes at the terminations of marine isotope stages (MIS) 6 and 2. The vertical slip rates, determined from offset measurements and surfaces ages, are 0.14 ?? 0.03 mm/yr over the late Pleistocene-Holocene and between 0.44 ?? 0.11 and 1.05 ?? 0.25 mm/yr since the end of the late Pleistocene. The higher of these slip rates for the last ???16 kyr is consistent with paleoseismic investigations along the fault [Prentice et al., 2002], and suggests that, at the end of late Pleistocene, the fault evolved from quiescence to having recurrence intervals of 4.0 ?? 1.2 kyr for surface ruptures with ???4 m vertical offset (similar to that of 1957). The inferred recurrence interval is comparable to that of the Bogd fault (3.7 ?? 1.3 kyr) suggesting that the two faults may have ruptured together also earlier during the last ???16 kyr.

Journal of Geophysical Research B: Solid Earth↗

High-resolution shallow reflection seismic image and surface evidence of the Upper Tiber Basin active faults (Northern Apennines, Italy)

Shallow seismic reflection prospecting has been carried out in order to investigate the faults that bound to the southwest and northeast the Quaternary Upper Tiber Basin (Northern Apennines, Italy). On the northeastern margin of the basin a ??? 1 km long reflection seismic profile images a fault segment and the associated up to 100 meters thick sediment wedge. Across the southwestern margin a 0.5 km-long seismic profile images a 50-55??-dipping extensional fault, that projects to the scarp at the base of the range-front, and against which a 100 m thick syn-tectonic sediment wedge has formed. The integration of surface and sub-surface data allows to estimate at least 190 meters of vertical displacement along the fault and a slip rate around 0.25 m/kyr. Southwestern fault might also be interpreted as the main splay structure of regional Alto Tiberina extensional fault. At last, the 1917 Monterchi earthquake (Imax=X, Boschi et alii, 2000) is correlable with an activation of the southwestern fault, and thus suggesting the seismogenic character of this latter.

Conference Paper↗

Rapid tsunami models and earthquake source parameters: Far-field and local applications

Rapid tsunami models have recently been developed to forecast far-field tsunami amplitudes from initial earthquake information (magnitude and hypocenter). Earthquake source parameters that directly affect tsunami generation as used in rapid tsunami models are examined, with particular attention to local versus far-field application of those models. First, validity of the assumption that the focal mechanism and type of faulting for tsunamigenic earthquakes is similar in a given region can be evaluated by measuring the seismic consistency of past events. Second, the assumption that slip occurs uniformly over an area of rupture will most often underestimate the amplitude and leading-wave steepness of the local tsunami. Third, sometimes large magnitude earthquakes will exhibit a high degree of spatial heterogeneity such that tsunami sources will be composed of distinct sub-events that can cause constructive and destructive interference in the wavefield away from the source. Using a stochastic source model, it is demonstrated that local tsunami amplitudes vary by as much as a factor of two or more, depending on the local bathymetry. If other earthquake source parameters such as focal depth or shear modulus are varied in addition to the slip distribution patterns, even greater uncertainty in local tsunami amplitude is expected for earthquakes of similar magnitude. Because of the short amount of time available to issue local warnings and because of the high degree of uncertainty associated with local, model-based forecasts as suggested by this study, direct wave height observations and a strong public education and preparedness program are critical for those regions near suspected tsunami sources.

ISET Journal of Earthquake Technology↗

Holocene geologic slip rate for the Banning strand of the southern San Andreas Fault, southern California

Northwest directed slip from the southern San Andreas Fault is transferred to the Mission Creek, Banning, and Garnet Hill fault strands in the northwestern Coachella Valley. How slip is partitioned between these three faults is critical to southern California seismic hazard estimates but is poorly understood. In this paper, we report the first slip rate measured for the Banning fault strand. We constrain the depositional age of an alluvial fan offset 25 ± 5 m from its source by the Banning strand to between 5.1 ± 0.4 ka (95% confidence interval (CI)) and 6.4 + 3.7/−2.1 ka (95% CI) using U-series dating of pedogenic carbonate clast coatings and 10 Be cosmogenic nuclide exposure dating of surface clasts. We calculate a Holocene geologic slip rate for the Banning strand of 3.9 + 2.3/−1.6 mm/yr (median, 95% CI) to 4.9 + 1.0/−0.9 mm/yr (median, 95% CI). This rate represents only 25–35% of the total slip accommodated by this section of the southern San Andreas Fault, suggesting a model in which slip is less concentrated on the Banning strand than previously thought. In rejecting the possibility that the Banning strand is the dominant structure, our results highlight an even greater need for slip rate and paleoseismic measurements along faults in the northwestern Coachella Valley in order to test the validity of current earthquake hazard models. In addition, our comparison of ages measured with U-series and 10 Be exposure dating demonstrates the importance of using multiple geochronometers when estimating the depositional age of alluvial landforms.

Journal of Geophysical Research B: Solid Earth↗

Lithospheric flexure under the Hawaiian volcanic load: Internal stresses and a broken plate revealed by earthquakes

Several lines of earthquake evidence indicate that the lithospheric plate is broken under the load of the island of Hawai`i, where the geometry of the lithosphere is circular with a central depression. The plate bends concave downward surrounding a stress-free hole, rather than bending concave upward as with past assumptions. Earthquake focal mechanisms show that the center of load stress and the weak hole is between the summits of Mauna Loa and Mauna Kea where the load is greatest. The earthquake gap at 21 km depth coincides with the predicted neutral plane of flexure where horizontal stress changes sign. Focal mechanism P axes below the neutral plane display a striking radial pattern pointing to the stress center. Earthquakes above the neutral plane in the north part of the island have opposite stress patterns; T axes tend to be radial. The M 6.2 Honomu and M 6.7 Kiholo main shocks (both at 39 km depth) are below the neutral plane and show radial compression, and the M 6.0 Kiholo aftershock above the neutral plane has tangential compression. Earthquakes deeper than 20 km define a donut of seismicity around the stress center where flexural bending is a maximum. The hole is interpreted as the soft center where the lithospheric plate is broken. Kilauea's deep conduit is seismically active because it is in the ring of maximum bending. A simplified two-dimensional stress model for a bending slab with a load at one end yields stress orientations that agree with earthquake stress axes and radial P axes below the neutral plane. A previous inversion of deep Hawaiian focal mechanisms found a circular solution around the stress center that agrees with the model. For horizontal faults, the shear stress within the bending slab matches the slip in the deep Kilauea seismic zone and enhances outward slip of active flanks.

Hawaii↗

Two models for earthquake forerunners

Similar precursory phenomena have been observed before earthquakes in the United States, the Soviet Union, Japan, and China. Two quite different physical models are used to explain these phenomena. According to a model developed by US seismologists, the so-called dilatancy diffusion model, the earthquake occurs near maximum stress, following a period of dilatant crack expansion. Diffusion of water in and out of the dilatant volume is required to explain the recovery of seismic velocity before the earthquake. According to a model developed by Soviet scientists growth of cracks is also involved but diffusion of water in and out of the focal region is not required. With this model, the earthquake is assumed to occur during a period of falling stress and recovery of velocity here is due to crack closure as stress relaxes. In general, the dilatancy diffusion model gives a peaked precursor form, whereas the dry model gives a bay form, in which recovery is well under way before the earthquake. A number of field observations should help to distinguish between the two models: study of post-earthquake recovery, time variation of stress and pore pressure in the focal region, the occurrence of pre-existing faults, and any changes in direction of precursory phenomena during the anomalous period.

Pure and Applied Geophysics PAGEOPH↗

High-frequency ground motion and source characteristics of the 2008 Wenchuan and 2013 Lushan, China, earthquakes

The 2008 M W 7.9 Wenchuan and the 2013 M W 6.6 Lushan earthquakes, which both occurred on the Longmen Shan thrust belt, show some interesting similarities and differences. Whereas the Wenchuan earthquake entailed a rupture zone that extended about 300 km northeastward, with fault slip extending to the surface, the Lushan earthquake was the result of a buried and much more compact zone of rupture. The high-frequency ground motions, however, for these two earthquakes, as measured by the peak ground acceleration, were evidently influenced by neither the extent of rupture nor the presence or absence of surface rupture. The source parameters for these two earthquakes tend to confirm the idea that high-frequency ground motion is controlled by stress changes in the rupture zone that give rise to the radiated ground acceleration. The apparent stresses for the Wenchuan and Lushan earthquakes are about 0.5 and 0.75 MPa, respectively, and the stress drops, in the same order, are about 2.5 and 3.5 MPa. The ratios of average stress drop to apparent stress are in the range 4.5–5 for both events, consistent with expectations based on the Brune (J Geophys Res 75(26):4997–5009, 1970) source model.

Lushan, Wenchuan↗

Seafloor gas seeps triggered by a northern California earthquake

A series of marine surveys using high-resolution seismic-reflection and side-scan sonar profiling was conducted on the Klamath River delta from 1977 to 1985. In November 1980, near the middle of this 8-year time period, a major earthquake (magnitude 7) occurred in the area. The earthquake caused liquefaction of shelf deposits off the Klamath River, and while mapping the sediment failures, we observed an inordinate number of gas vents and acoustic anomalies in the sediment (gas charging) and water column (gas clouds). Our analyses indicate that a change in gas seepage occurred after the earthquake. There was a marked increase in gas abundance immediately after the earthquake, a gradual decrease over the next year, and stabilization at pre-earthquake gas levels within five years or less. The acoustic anomalies were observed in pre-Quaternary rocks as well as in surficial Holocene deposits, and gas seeps overlay truncated anticlines, faults, and syncline exposures. Thus, the gas released may have been in part thermogenic in origin. A large area (20 km 2 ) on the delta failed during the earthquake, and the abundance of gas seeps and craters in the area suggests that sediment degassing may have been an important process contributing to failure.

Marine Geology↗

Spatial distribution of landslides triggered from the 2007 Niigata Chuetsu–Oki Japan Earthquake

Understanding the spatial distribution of earthquake-induced landslides from specific earthquakes provides an opportunity to recognize what to expect from future events. The July 16, 2007 M w 6.6 (M JMA 6.8) Niigata Chuetsu–Oki Japan earthquake triggered hundreds of landslides in the area surrounding the coastal city of Kashiwazaki and provides one such opportunity to evaluate the impacts of an offshore, magnitude 6 + earthquake on a steep coastal region. As part of a larger effort to document all forms of geotechnical damage from this earthquake, we performed landslide inventory mapping throughout the epicentral area and analyzed the resulting data for spatial, seismic-motion, and geologic correlations to describe the pattern of landsliding. Coupled with examination of a third-party, aerial-photo-based landslide inventory, our analyses reveal several areas of high landslide concentration that are not readily explained by either traditional epicentral and fault–plane-distance metrics or by recorded and inferred ground-motions. Whereas average landslide concentrations averaged less than 1 landslide per square kilometer (LS/km 2 ), some areas reached up to 2 LS/km 2 in the Nishiyama Hills to the northeast of Kashiwazaki and between 2 and 11 LS/km 2 in coastal areas to the north and south of the city. Correlation with seismometer-based and monument overturning back-calculated ground motions suggests that a minimum peak ground acceleration (PGA) of approximately 0.2 g was necessary for landsliding throughout the region, but does not explain the subregional areas of high landslide concentration. However, analysis of topographic slope and the distribution of generally weak, dip-slope, geologic units does sufficiently explain why, on a sub-regional scale, high landslide concentrations occurred where they did. These include: (1) an inland region of steep, dip-slope, anticlinal sedimentary strata with associated fold belt compression and uplift of the anticline and (2) coastal areas with generally weaker, weathered outcrop lithology and steeper slopes resulting from active marine and terrestrial cliff processes. The results offer lessons for understanding the effects of earthquakes on both regional and subregional scales with regard to the spatial distribution of landsliding.

Engineering Geology↗

Great-earthquake paleogeodesy and tsunamis of the past 2000 years at Alsea Bay, central Oregon coast, USA

The width of plate-boundary fault rupture at the Cascadia subduction zone, a dimension related to earthquake magnitude, remains uncertain because of the lack of quantitative information about land-level movements during past great-earthquake deformation cycles. Beneath a marsh at Alsea Bay, on the central Oregon coast, four sheets of tsunami-deposited sand blanket contacts between tidal mud and peat. Radiocarbon ages for the sheets match ages for similar evidence of regional coseismic subsidence and tsunamis during four of Cascadia's great earthquakes. Barring rapid, unrecorded postseismic uplift, reconstruction of changes in land level from core samples using diatom and foraminiferal transfer functions includes modest coseismic subsidence (0.4??0.2 m) during the four earthquakes. Interpretation is complicated, however, by the 30-38% of potentially unreliable transfer function values from samples with poor analogs in modern diatom and foraminiferal assemblages. Reconstructions of coseismic subsidence using good-analog samples range from 0.46??0.12 to 0.09??0.20 m showing greater variability than implied by sample-specific errors. From apparent high rates of land uplift following subsidence and tsunamis, we infer that postseismic rebound caused by slip on deep parts of the plate boundary and (or) viscoelastic stress relaxation in the upper plate may be almost as large as coseismic subsidence. Modest coseismic subsidence 100 km landward of the deformation front implies that plate-boundary ruptures in central Oregon were largely offshore. Ruptures may have been long and narrow during earthquakes near magnitude 9, as suggested for the AD 1700 earthquake, or of smaller and more variable dimensions and magnitudes. ?? 2008 Elsevier Ltd. All rights reserved.

Quaternary Science Reviews↗

Slip-parallel seismic lineations on the Northern Hayward Fault, California

A high-resolution relative earthquake location procedure is used to image the fine-scale seismicity structure of the northern Hayward fault, California. The seismicity defines a narrow, near-vertical fault zone containing horizontal alignments of hypocenters extending along the fault zone. The lineations persist over the 15-year observation interval, implying the localization of conditions on the fault where brittle failure conditions are met. The horizontal orientation of the lineations parallels the slip direction of the fault, suggesting that they are the result of the smearing of frictionally weak material along the fault plane over thousands of years.

Geophysical Research Letters↗

Rupture imaging of the M w 7.9 12 May 2008 Wenchuan earthquake from back projection of teleseismic P waves

[1] The M w 7.9 Wenchuan earthquake of 12 May 2008 was the most destructive Chinese earthquake since the 1976 Tangshan event. Tens of thousands of people were killed, hundreds of thousands were injured, and millions were left homeless. Here we infer the detailed rupture process of the Wenchuan earthquake by back-projecting teleseismic P energy from several arrays of seismometers. This technique has only recently become feasible and is potentially faster than traditional finite-fault inversion of teleseismic body waves; therefore, it may reduce the notification time to emergency response agencies. Using the IRIS DMC, we collected 255 vertical component broadband P waves at 30-95?? from the epicenter. We found that at periods of 5 s and greater, nearly all of these P waves were coherent enough to be used in a global array. We applied a simple down-sampling heuristic to define a global subarray of 70 stations that reduced the asymmetry and sidelobes of the array response function (ARF). We also considered three regional subarrays of seismometers in Alaska, Australia, and Europe that had apertures less than 30?? and P waves that were coherent to periods as short as 1 s. Individual ARFs for these subarrays were skewed toward the subarrays; however, the linear sum of the regional subarray beams at 1 s produced a symmetric ARF, similar to that of the groomed global subarray at 5 s. For both configurations we obtained the same rupture direction, rupture length, and rupture time. We found that the Wenchuan earthquake had three distinct pulses of high beam power at 0, 23, and 57 s after the origin time, with the pulse at 23 s being highest, and that it ruptured unilaterally to the northeast for about 300 km and 110 s, with an average speed of 2.8 km/s. It is possible that similar results can be determined for future large dip-slip earthquakes within 20-30 min of the origin time using relatively sparse global networks of seismometers such as those the USGS uses to locate earthquakes in near-real time. Copyright 2009 by the American Geophysical Union.

Geochemistry, Geophysics, Geosystems↗

Precise tremor source locations and amplitude variations along the lower-crustal central San Andreas Fault

We precisely locate 88 tremor families along the central San Andreas Fault using a 3D velocity model and numerous P and S wave arrival times estimated from seismogram stacks of up to 400 events per tremor family. Maximum tremor amplitudes vary along the fault by at least a factor of 7, with by far the strongest sources along a 25 km section of the fault southeast of Parkfield. We also identify many weaker tremor families, which have largely escaped prior detection. Together, these sources extend 150 km along the fault, beneath creeping, transitional, and locked sections of the upper crustal fault. Depths are mostly between 18 and 28 km, in the lower crust. Epicenters are concentrated within 3 km of the surface trace, implying a nearly vertical fault. A prominent gap in detectible activity is located directly beneath the region of maximum slip in the 2004 magnitude 6.0 Parkfield earthquake.

California↗

Plastic faulting in ice

Plastic faulting is a brittle‐like failure phenomenon exhibited by water ice and several other rock types under confinement. It is suspected to be the mechanism of deep earthquakes and extreme cases of shear localization in shallow rocks. Unlike ordinary Coulombic failure, plastic faulting is characterized by a pressure‐independent failure strength and fault plane oriented 45° to maximum principal stress. To research the question of how the instability initiates, we conducted over 50 constant‐displacement‐rate experiments on polycrystalline ice (phases Ih and II) near the brittle‐to‐ductile (B‐D) transition, at confining pressures P = 0–300 MPa, applied strain rates = 5 × 10 −5 – 7 × 10 −3 s −1 , temperatures T = 105–233 K, and mean grain sizes d = 0.25–1.18 mm. We find that (1) the width of the B‐D transition in variable space is vanishingly narrow, to the point of appearing as a crossover, (2) a plastic fault plane, once formed, is not a zone of subsequent weakness, (3) distributed ice I→II phase transformation in small amounts (<1 vol%) shows no causal relationship to subsequent failure, and (4) plastic faulting also occurs in ice II. We hypothesize that the elusive nucleating “trigger” parallels that of metals and ceramics undergoing severe plastic deformation, wherein transient local structural rearrangement occurs, in turn causing material strength to drop to a level sufficiently low, in a volume sufficiently large, that adiabatic instability is nucleated. Our results do not require and often are inconsistent with phase transformation. Plastic faulting may therefore be available to all solids undergoing severe deformation, and its appearance in so few is simply the result of insufficiently extreme conditions.

Journal of Geophysical Research- Solid Earth↗

Source parameters of the 1980 Mammoth Lakes, California, earthquake sequence

From the more than 1500 Mammoth Lakes earthquakes recorded on three-component digital seismographs (Spudich et al., 1981), 150 were used in an analysis of the locations, mechanism, and source parameters. A composite fault plane solution of nine earthquakes 3.9 ≤ M ≤ 5.1 defines a right-lateral strike slip mechanism on a steeply dipping nearly east-west plane striking S75°E or left-lateral strike slip on a nearly north-south plane striking N10°E. Vertical cross sections of well-located aftershocks indicate possibly three east-west planes that coincide with the locations of the four largest earthquakes with M L ≥ 6.0. Using the spectral analysis of S waves (Brune, 1970), source parameters for 67 earthquakes were determined. Forty-eight had magnitudes greater than or equal to 3.0. Seismic moments ranged from 9.20×10 18 dyn cm to 2.33×10 24 dyn cm. Earthquakes with seismic moment greater than about 1.0×10 21 dyn cm had nearly constant stress drops (≃ 50 bars); earthquakes with seismic moment less than about 1.0×10 21 dyn cm had stress drops that apparently decrease as seismic moment decreases.

California↗

Variations in fault slip and strain accumulation at Parkfield, California: Initial results using two-color geodimeter measurements 1984-1988

Repeated length measurements of several geodetic baselines near Parkfield, California, have revealed significant variations in the local rates of shallow fault slip and strain. This network of baselines, surveyed several times each week, straddles the San Andreas fault in the transition zone between the creeping section to the northwest and the locked section to the southeast. The length measurements, characterized by a precision approaching 0.1 ppm, reveal large fluctuations in the rates of baseline extension. Principal mode analysis of the length change data indicates that the two largest components of the signal are (1) secular extensions and contractions consistent with surface slip on the main strand of the San Andreas fault, and (2) a large seasonal oscillation with no obvious spatial coherence. On most of these baselines, the second component appears to be in phase with seasonal rainfall. When data from the baselines with the largest amplitude of the seasonal signal are excluded, the remaining data can be modeled in terms of both spatial and temporal variations in surface slip, variations in the components of the changes in uniform strain, and the possible displacement of the central monument in this radial network. In parameterizing this model, the spatial variation of slip beneath the near surface is reflected by changes in shear strain. Although the computed secular shear is highly dependent upon the specified parameterization of surface slip, the data are consistent with the hypothesis that slip at intermediate depths lags behind the surface slip rate. However, the range in models that fit the data does not necessarily imply that there is a deficit in slip at depth relative to the surface. Comparison of the inferred values of surface slip from the model with the observed fault slip measured by very short baseline creep meters indicates close agreement in secular rates, but the short-term variations observed with the creep meters are either highly attenuated or nonexistent in the modeled slip since the modeled slip is a spatial average which smooths out possible short-wavelength variations in the surface slip for which the creep instruments are most sensitive. An interesting conclusion from the two-color data is that surface slip on the San Andreas fault appears to be spread over a 2-km-wide zone on the south flank of Middle Mountain but is confined to a very narrow zone to the south as the fault passes through the center of the network. This conclusion is dependent upon the assumption that a few critical monuments are stable and track tectonic displacements in the long term. Finally, the largest observed strain change is an extensional strain coincident with the Kettleman Hills earthquake M 5.5 in August 1985.

Journal of Geophysical Research Solid Earth↗

Earthquake recurrence and risk assessment in circum-Pacific seismic gaps

THE development of the concept of seismic gaps, regions of low earthquake activity where large events are expected, has been one of the notable achievements of seismology and plate tectonics. Its application to long-term earthquake hazard assessment continues to be an active field of seismological research. Here I have surveyed well documented case histories of repeated rupture of the same segment of circum-Pacific plate boundary and characterized their general features. I find that variability in fault slip and spatial extent of great earthquakes rupturing the same plate boundary segment is typical rather than exceptional but sequences of major events fill identified seismic gaps with remarkable order. Earthquakes are concentrated late in the seismic cycle and occur with increasing size and magnitude. Furthermore, earthquake rupture starts near zones of concentrated moment release, suggesting that high-slip regions control the timing of recurrent events. The absence of major earthquakes early in the seismic cycle indicates a more complex behaviour for lower-slip regions, which may explain the observed cycle-to-cycle diversity of gap-filling sequences.

Nature↗

Book review: Earthquakes and water

It is really nice to see assembled in one place a discussion of the documented and hypothesized hydrologic effects of earthquakes. The book is divided into chapters focusing on particular hydrologic phenomena including liquefaction, mud volcanism, stream discharge increases, groundwater level, temperature and chemical changes, and geyser period changes. These hydrologic effects are inherently fascinating, and the large number of relevant publications in the past decade makes this summary a useful milepost. The book also covers hydrologic precursors and earthquake triggering by pore pressure. A natural need to limit the topics covered resulted in the omission of tsunamis and the vast literature on the role of fluids and pore pressure in frictional strength of faults. Regardless of whether research on earthquake-triggered hydrologic effects ultimately provides insight into the physics of earthquakes, the text provides welcome common ground for interdisciplinary collaborations between hydrologists and seismologists. Such collaborations continue to be crucial for investigating hypotheses about the role of fluids in earthquakes and slow slip. Review info: Earthquakes and Water. By Wang, C.-Y. and Manga, M., 2010. ISBN: 9783642008092, 218 pp.

Geofluids↗