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 379 records · Page 21Linked to original sources

Slip triggered on southern California faults by the 1992 Joshua Tree, Landers, and big bear earthquakes

Five out of six functioning creepmeters on southern California faults recorded slip triggered at the time of some or all of the three largest events of the 1992 Landers earthquake sequence. Digital creep data indicate that dextral slip was triggered within 1 min of each mainshock and that maximum slip velocities occurred 2 to 3 min later. The duration of triggered slip events ranged from a few hours to several weeks. We note that triggered slip occurs commonly on faults that exhibit fault creep. To account for the observation that slip can be triggered repeatedly on a fault, we propose that the amplitude of triggered slip may be proportional to the depth of slip in the creep event and to the available near-surface tectonic strain that would otherwise eventually be released as fault creep. We advance the notion that seismic surface waves, perhaps amplified by sediments, generate transient local conditions that favor the release of tectonic strain to varying depths. Synthetic strain seismograms are presented that suggest increased pore pressure during periods of fault-normal contraction may be responsible for triggered slip, since maximum dextral shear strain transients correspond to times of maximum fault-normal contraction.

California↗

Characterizing the Los Angeles Aqueduct crossing of the San Andreas Fault for improved earthquake resilience

The five-mile-long Elizabeth Tunnel, which crosses the San Andreas fault (SAF) zone near Lake Hughes, California, is part of the Los Angeles Aqueduct (LAA) that delivers water from Owens Valley to the City of Los Angeles. Geologic characterization of the Elizabeth Tunnel alignment is focused on developing a better understanding of fault displacement hazards at the SAF crossing to support design of both short- and long-term strategies to increase the earthquake resilience of the City’s water supply system. The results of this study define a fault zone that is wider and more complex at the surface than at tunnel depth. A 750-ft transect of deep, angled core borings, located ~ 100 ft west of, and parallel to the tunnel alignment, penetrated a wide zone (>550 ft) of fault-damaged granitic and gneissic rocks containing a few thick fault zones that can be correlated confidently between borings as well as a multitude of very thin faults and shears throughout the rock mass. Fault rocks are predominantly defined as coherent cataclasites and ultracataclasites, as opposed to incoherent clay gouge, even at this relatively shallow sampling depth. Integrating geomorphic, geologic, petrographic, and geophysical data defines a steeply south-dipping SAF that splays upward and increases in complexity toward the surface, and provides the basis for engineering design considerations at this critical SAF crossing.

California↗

In search of earthquake-related hydrologic and chemical changes along Hayward Fault

Flow and chemical measurements have been made about once a month, and more frequently when required, since 1976 at two springs in Alum Rock Park in eastern San Jose, California, and since 1980 at two shallow wells in eastern Oakland in search of earthquake-related changes. All sites are on or near the Hayward Fault and are about 55 km apart. Temperature, electric conductivity, and water level or flow rate were measured in situ with portable instruments. Water samples were collected for later chemical and isotopic analyses in the laboratory. The measured flow rate at one of the springs showed a long-term decrease of about 40% since 1987, when a multi-year drought began in California. It also showed several increases that lasted a few days to a few months with amplitudes of 2.4 to 8.6 times the standard deviations above the background rate. Five of these increases were recorded shortly after nearby earthquakes of magnitude 5.0 or larger, and may have resulted from unclogging of the flow path and increase of permeability caused by strong seismic shaking. Two other flow increases were possibly induced by exceptionally heavy rainfalls. The water in both wells showed seasonal temperature and chemical variations, largely in response to rainfall. In 1980 the water also showed some clear chemical changes unrelated to rainfall that lasted a few months; these changes were followed by a magnitude 4 earthquake 37 km away. The chemical composition at one of the wells and at the springs also showed some longer-term variations that were not correlated with rainfall but possibly correlated with the five earthquakes mentioned above. These correlations suggest a common tectonic origin for the earthquakes and the anomalies. The last variation at the affected well occurred abruptly in 1989, shortly before a magnitude 5.0 earthquake 54 km away. ?? 1993.

California↗

A slow earthquake sequence on the San Andreas fault

EARTHQUAKES typically release stored strain energy on timescales of the order of seconds, limited by the velocity of sound in rock. Over the past 20 years, observations and laboratory experiments have indicated that capture can also occur more slowly, with durations up to hours. Such events may be important in earthquake nucleation and in accounting for the excess of plate convergence over seismic slip in subduction zones. The detection of events with larger timescales requires near-field deformation measurements. In December 1992, two borehole strainmeters close to the San Andreas fault in California recorded a slow strain event of about a week in duration, and we show here that the strain changes were produced by a slow earthquake sequence (equivalent magnitude 4.8) with complexity similar to that of regular earthquakes. The largest earthquakes associated with these slow events were small (local magnitude 3.7) and contributed negligible strain release. The importance of slow earthquakes in the seismogenic process remains an open question, but these observations extend the observed timescale for slow events by two orders of magnitude.

Nature↗

Implications of the earthquake cycle for inferring fault locking on the Cascadia megathrust

GPS velocity fields in the Western US have been interpreted with various physical models of the lithosphere-asthenosphere system: (1) time-independent block models; (2) time-dependent viscoelastic-cycle models, where deformation is driven by viscoelastic relaxation of the lower crust and upper mantle from past faulting events; (3) viscoelastic block models, a time-dependent variation of the block model. All three models are generally driven by a combination of loading on locked faults and (aseismic) fault creep. Here we construct viscoelastic block models and viscoelastic-cycle models for the Western US, focusing on the Pacific Northwest and the earthquake cycle on the Cascadia megathrust. In the viscoelastic block model, the western US is divided into blocks selected from an initial set of 137 microplates using the method of Total Variation Regularization, allowing potential trade-offs between faulting and megathrust coupling to be determined algorithmically from GPS observations. Fault geometry, slip rate, and locking rates (i.e. the locking fraction times the long term slip rate) are estimated simultaneously within the TVR block model. For a range of mantle asthenosphere viscosity (4.4 × 10 18 to 3.6 × 10 20 Pa s) we find that fault locking on the megathrust is concentrated in the uppermost 20 km in depth, and a locking rate contour line of 30 mm yr −1 extends deepest beneath the Olympic Peninsula, characteristics similar to previous time-independent block model results. These results are corroborated by viscoelastic-cycle modelling. The average locking rate required to fit the GPS velocity field depends on mantle viscosity, being higher the lower the viscosity. Moreover, for viscosity ≲ 10 20 Pa s, the amount of inferred locking is higher than that obtained using a time-independent block model. This suggests that time-dependent models for a range of admissible viscosity structures could refine our knowledge of the locking distribution and its epistemic uncertainty.

Geophysical Journal International↗

Late Holocene slip rate of the San Andreas fault and its accommodation by creep and moderate-magnitude earthquakes at Parkfield, California

Investigation of a right-laterally offset channel at the Miller's Field paleoseismic site yields a late Holocene slip rate of 26.2 +6.4/−4.3 mm/yr (1σ) for the main trace of the San Andreas fault at Parkfield, California. This is the first well-documented geologic slip rate between the Carrizo and creeping sections of the San Andreas fault. This rate is lower than Holocene measurements along the Carrizo Plain and rates implied by far-field geodetic measurements (∼35 mm/yr). However, the rate is consistent with historical slip rates, measured to the northwest, along the creeping section of the San Andreas fault (<30 mm/yr). The paleoseismic exposures at the Miller's Field site reveal a pervasive fabric of clay shear bands, oriented clockwise oblique to the San Andreas fault strike and extending into the uppermost stratigraphy. This fabric is consistent with dextral aseismic creep and observations of surface slip from the 28 September 2004 M6 Parkfield earthquake. Together, this slip rate and deformation fabric suggest that the historically observed San Andreas fault slip behavior along the Parkfield section has persisted for at least a millennium, and that significant slip is accommodated by structures in a zone beyond the main San Andreas fault trace.

California↗

Slip rate on the San Diego trough fault zone, inner California Borderland, and the 1986 Oceanside earthquake swarm revisited

The San Diego trough fault zone (SDTFZ) is part of a 90-km-wide zone of faults within the inner California Borderland that accommodates motion between the Pacific and North American plates. Along with most faults offshore southern California, the slip rate and paleoseismic history of the SDTFZ are unknown. We present new seismic reflection data that show that the fault zone steps across a 5-km-wide stepover to continue for an additional 60 km north of its previously mapped extent. The 1986 Oceanside earthquake swarm is located within the 20-km-long restraining stepover. Farther north, at the latitude of Santa Catalina Island, the SDTFZ bends 20&deg; to the west and may be linked via a complex zone of folds with the San Pedro basin fault zone (SPBFZ). In a cooperative program between the U.S. Geological Survey (USGS) and the Monterey Bay Aquarium Research Institute (MBARI), we measure and date the coseismic offset of a submarine channel that intersects the fault zone near the SDTFZ&ndash;SPBFZ junction. We estimate a horizontal slip rate of about 1:5  0:3 mm=yr over the past 12,270 yr.

California↗

S -wave triggering of tremor beneath the Parkfield, California, section of the San Andreas fault by the 2011 Tohoku, Japan earthquake: observations and theory

The dynamic stresses that are associated with the energetic seismic waves generated by the M w 9.0 Tohoku earthquake off the northeast coast of Japan triggered bursts of tectonic tremor beneath the Parkfield section of the San Andreas fault (SAF) at an epicentral distance of ∼8200 km. The onset of tremor begins midway through the ∼100‐s‐period S‐wave arrival, with a minor burst coinciding with the SHSH arrival, as recorded on the nearby broadband seismic station PKD. A more pronounced burst coincides with the Love arrival, followed by a series of impulsive tremor bursts apparently modulated by the 20‐ to 30‐s‐period Rayleigh wave. The triggered tremor was located at depths between 20 and 30 km beneath the surface trace of the fault, with the burst coincident with the S wave centered beneath the fault 30 km northwest of Parkfield. Most of the subsequent activity, including the tremor coincident with the SHSH arrival, was concentrated beneath a stretch of the fault extending from 10 to 40 km southeast of Parkfield. The seismic waves from the Tohoku epicenter form a horizontal incidence angle of ∼14°, with respect to the local strike of the SAF. Computed peak dynamic Coulomb stresses on the fault at tremor depths are in the 0.7–10 kPa range. The apparent modulation of tremor bursts by the small, strike‐parallel Rayleigh‐wave stresses (∼0.7 kPa) is likely enabled by pore pressure variations driven by the Rayleigh‐wave dilatational stress. These results are consistent with the strike‐parallel dynamic stresses (δτ s ) associated with the S, SHSH, and surface‐wave phases triggering small increments of dextral slip on the fault with a low friction (μ∼0.2). The vertical dynamic stresses δτ d do not trigger tremor with vertical or oblique slip under this simple Coulomb failure model.

California↗

Basement and regional structure along strike of the Queen Charlotte Fault in the context of modern and historical earthquake ruptures

The Queen Charlotte fault (QCF) is a dextral transform system located offshore of southeastern Alaska and western Canada, accommodating &sim;4.4&thinsp;&thinsp;cm/yr of relative motion between the Pacific and North American plates. Oblique convergence along the fault increases southward, and how this convergence is accommodated is still debated. Using seismic reflection data, we interpret offshore basement structure, faulting, and stratigraphy to provide a geological context for two recent earthquakes, an M w 7.5 strike‐slip event near Craig, Alaska, and an M w 7.8 thrust event near Haida Gwaii, Canada. We map downwarped Pacific oceanic crust near 54&deg; N, between the two rupture zones. Observed downwarping decreases north and south of 54&deg; N, parallel to the strike of the QCF. Bending of the Pacific plate here may have initiated with increased convergence rates due to a plate motion change at &sim;6&thinsp;&thinsp;Ma. Tectonic reconstruction implies convergence‐driven Pacific plate flexure, beginning at 6 Ma south of a 10&deg; bend the QCF (which is currently at 53.2&deg; N) and lasting until the plate translated past the bend by &sim;2&thinsp;&thinsp;Ma. Normal‐faulted approximately late Miocene sediment above the deep flexural depression at 54&deg; N, topped by relatively undeformed Pleistocene and younger sediment, supports this model. Aftershocks of the Haida Gwaii event indicate a normal‐faulting stress regime, suggesting present‐day plate flexure and underthrusting, which is also consistent with reconstruction of past conditions. We thus favor a Pacific plate underthrusting model to initiate flexure and accommodation space for sediment loading. In addition, mapped structures indicate two possible fault segment boundaries along the QCF at 53.2&deg; N and at 56&deg; N.

Alaska↗

The Catfish Lake Scarp, Allyn, Washington: Preliminary field data and implications for earthquake hazards posed by the Tacoma fault

The Tacoma fault bounds gravity and aeromagnetic anomalies for 50 km across central Puget lowland from Tacoma to western Kitsap County. Tomography implies at least 6 km of post-Eocene uplift to the north of the fault relative to basinal sedimentary rocks to the south. Coastlines north of the Tacoma fault rose about 1100 years ago during a large earthquake. Abrupt uplift up to several meters caused tidal flats at Lynch Cove, North Bay, and Burley Lagoon to turn into forested wetlands and freshwater marshes. South of the fault at Wollochet Bay, Douglas-fir forests sank into the intertidal zone and changed into saltmarsh. Liquefaction features found beneath the marsh at Burley Lagoon point to strong ground shaking at the time of uplift. Recent lidar maps of the area southwest of Allyn, Washington revealed a 4 km long scarp, or two closely spaced en-echelon scarps, which correspond closely to the Tacoma fault gravity and aeromagnetic anomalies. The scarp, named the Catfish Lake scarp, is north-side-up, trends east-west, and clearly displace striae left by a Vashon-age glacier. A trench across the scarp exposed evidence for postglacial folding and reverse slip. No organic material for radiocarbon dating was recovered from the trench. However, relationships in the trench suggest that the folding and faulting is postglacial in age.

Washington↗

New fault picture points toward San Francisco Bay area earthquakes

Recent earthquakes and a new way of looking at faults suggest that damaging earthquakes are closing in on the San Francisco area. Earthquakes Awareness Week 1989 in northern California started off with a bang on Monday, 3 April, when a magnitude 4.8 earthquake struck 15 kilometers northeast of San Jose. The relatively small shock-its primary damage was the shattering of an air-control tower window-got the immediate attention of three U.S Geological Survey seismologists in Menlo Park near San Francisco. David Oppenheimer, William Bakun, and Allan Lindh had forecast a nearby earthquake in a just completed report, and this, they thought, might be it.

California↗

A 15-year catalog of more than 1 million low-frequency earthquakes: Tracking tremor and slip along the deep San Andreas Fault

Low-frequency earthquakes (LFEs) are small, rapidly recurring slip events that occur on the deep extensions of some major faults. Their collective activation is often observed as a semi-continuous signal known as tectonic (or non-volcanic) tremor. This manuscript presents a catalog of more than 1 million LFEs detected along the central San Andreas Fault from 2001-2016. These events have been detected via a multi-channel matched-filter search, cross-correlating waveform templates representing 88 different LFE families with continuous seismic data. Together, these source locations span nearly 150 km along the central San Andreas Fault, ranging in depth from ~16-30 km. This accumulating catalog has been the source of numerous studies examining the behavior of these LFE sources and the inferred slip behavior of the deep fault. The relatively high temporal and spatial resolution of the catalog has provided new insights into properties such as tremor migration, recurrence, and triggering by static and dynamic stress perturbations. Collectively, these characteristics are inferred to reflect a very weak fault likely under near-lithostatic fluid pressure, yet the physical processes controlling the stuttering rupture observed as tremor and LFE signals remain poorly understood. This paper aims to document the LFE catalog assembly process and associated caveats, while also updating earlier observations and inferred physical constraints. The catalog itself accompanies this manuscript as part of the electronic supplement, with the goal of providing a useful resource for continued future investigations.

California↗

Constraints on friction, dilatancy, diffusivity, and effective stress from low-frequency earthquake rates on the deep San Andreas Fault

Families of recurring low-frequency earthquakes (LFEs) within nonvolcanic tremor on the San Andreas Fault in central California are sensitive to tidal stresses. LFEs occur at all levels of the tides, are strongly correlated and in phase with the ~200 Pa shear stresses, and weakly and not systematically correlated with the ~2 kPa tidal normal stresses. We assume that LFEs are small sources that repeatedly fail during shear within a much larger scale, aseismically slipping fault zone and consider two different models of the fault slip: (1) modulation of the fault slip rate by the tidal stresses or (2) episodic slip, triggered by the tides. LFEs are strongly clustered with duration much shorter than the semidiurnal tide; they cannot be significantly modulated on that time scale. The recurrence times of clusters, however, are many times longer than the semidiurnal, leading to an appearance of tidal triggering. In this context we examine the predictions of laboratory-observed triggered frictional (dilatant) fault slip. The undrained end-member model produces no sensitivity to the tidal normal stress, and slip onsets are in phase with the tidal shear stress. The tidal correlation constrains the diffusivity to be less than ~1 × 10 −6 /s and the product of the friction and dilatancy coefficients to be at most 5 × 10 −7 , orders of magnitude smaller than observed at room temperature. In the absence of dilatancy the effective normal stress at failure would be about ~55 kPa. For this model the observations require intrinsic weakness, low dilatancy, and lithostatic pore fluid.

Journal of Geophysical Research↗

Aftershocks of the June 20, 1978, Greece earthquake: A multimode faulting sequence

A 10-station portable seismograph network was deployed in northern Greece to study aftershocks of the magnitude (m b ) 6.4 earthquake of June 20, 1978. The main shock occurred (in a graben) about 25 km northeast of the city of Thessaloniki and caused an east-west zone of surface rupturing 14 km long that splayed to 7 km wide at the west end. The hypocenters for 116 aftershocks in the magnitude range from 2.5 to 4.5 were determined. The epicenters for these events cover an area 30 km (east-west) by 18 km (north-south), and focal depths ranges from 4 to 12 km. Most of the aftershocks in the east half of the aftershock zone are north of the surface rupture and north of the graben. Those in the west half are located within the boundaries of the graben. Composite focalmechanism solutions for selected aftershocks indicate reactivation of geologically mapped normal faults in the area. Also, strike-slip and dip-slip faults that splay off the western end of the zone of surface ruptures may have been activated. The epicenters for four large (M ≯ 4.8) foreshocks and the main shock were relocated using the method of joint epicenter determination. Collectively, those five epicenters form an arcuate pattern convex southward, that is north of and 5 km distant from the surface rupturing. The 5-km separation, along with a focal depth of 8 km (average aftershock depth) or 16 km (NEIS main-shock depth), implies that the fault plane dips northward 58° or 73°, respectively. A preferred nodal-plane dip of 36° was determined by B.C. Papazachos and his colleagues in 1979 from a focal-mechanism solution for the main shock. If this dip is valid for the causal fault and that fault projects to the zone of surface rupturing, a decrease of dip with depth is required.

Tectonophysics↗

A reevaluation of the Pallett Creek earthquake chronology based on new AMS radiocarbon dates, San Andreas fault, California

The Pallett Creek paleoseismic record occupies a keystone position in most attempts to develop rupture histories for the southern San Andreas fault. Previous estimates of earthquake ages at Pallett Creek were determined by decay counting radiocarbon methods. That method requires large samples which can lead to unaccounted sources of uncertainty in radiocarbon ages because of the heterogeneous composition of organic layers. In contrast, accelerator mass spectrometry (AMS) radiocarbon dates may be obtained from small samples that have known carbon sources and also allow for a more complete sampling of the section. We present 65 new AMS radiocarbon dates that span nine ground‐rupturing earthquakes at Pallett Creek. Overall, the AMS dates are similar to and reveal no dramatic bias in the conventional dates. For many layers, however, individual charcoal samples were younger than the conventional dates, leading to earthquake ages that are overall slightly younger than previously reported. New earthquake ages are determined by Bayesian refinement of the layer ages based on stratigraphic ordering and sedimentological constraints. The new chronology is more regular than previously published records in large part due to new samples constraining the age of event R. The closed interval from event C to 1857 has a mean recurrence of 135 years ( σ = 83.2 years) and a quasiperiodic coefficient of variation (COV) of 0.61. We show that the new dates and resultant earthquake chronology have a stronger effect on COV than the specific membership of this long series and dating precision improvements from sedimentation rates.

California↗

Periodic, chaotic, and doubled earthquake recurrence intervals on the deep San Andreas Fault

Earthquake recurrence histories may provide clues to the timing of future events, but long intervals between large events obscure full recurrence variability. In contrast, small earthquakes occur frequently, and recurrence intervals are quantifiable on a much shorter time scale. In this work, I examine an 8.5-year sequence of more than 900 recurring low-frequency earthquake bursts composing tremor beneath the San Andreas fault near Parkfield, California. These events exhibit tightly clustered recurrence intervals that, at times, oscillate between ~3 and ~6 days, but the patterns sometimes change abruptly. Although the environments of large and low-frequency earthquakes are different, these observations suggest that similar complexity might underlie sequences of large earthquakes.

California↗

Multiple seismogenic processes for high-frequency earthquakes at Katmai National Park, Alaska: Evidence from stress tensor inversions of fault-plane solutions

The volcanological significance of seismicity within Katmai National Park has been debated since the first seismograph was installed in 1963, in part because Katmai seismicity consists almost entirely of high-frequency earthquakes that can be caused by a wide range of processes. I investigate this issue by determining 140 well-constrained first-motion fault-plane solutions for shallow (depth < 9 km) earthquakes occurring between 1995 and 2001 and inverting these solutions for the stress tensor in different regions within the park. Earthquakes removed by several kilometers from the volcanic axis occur in a stress field characterized by horizontally oriented σ 1 and σ 3 axes, with σ 1 rotated slightly (12°) relative to the NUVEL1A subduction vector, indicating that these earthquakes are occurring in response to regional tectonic forces. On the other hand, stress tensors for earthquake clusters beneath several Katmai cluster volcanoes have vertically oriented σ 1 axes, indicating that these events are occurring in response to local, not regional, processes. At Martin-Mageik, vertically oriented σ 1 is most consistent with failure under edifice loading conditions in conjunction with localized pore pressure increases associated with hydrothermal circulation cells. At Trident-Novarupta, it is consistent with a number of possible models, including occurrence along fractures formed during the 1912 eruption that now serve as horizontal conduits for migrating fluids and/or volatiles from nearby degassing and cooling magma bodies. At Mount Katmai, it is most consistent with continued seismicity along ring-fracture systems created in the 1912 eruption, perhaps enhanced by circulating hydrothermal fluids and/or seepage from the caldera-filling lake.

Alaska↗

Kinematics of fault slip associated with the July 4-6 2019 Ridgecrest, Californai earthquakes sequence

The 2019 Ridgecrest, California, earthquake sequence produced observable crustal deformation over much of central and southern California, as well as surface rupture over several tens of kilometers. To obtain a detailed picture of the fault slip involved in the 4 July M 6.4 foreshock and 6 July M 7.1 mainshock, we combine strong‐motion seismic waveforms with crustal deformation observations to obtain kinematic and static slip models of both events. We sample the regional seismic wavefield for both the foreshock and mainshock with three‐component records from 31 stations of the California Integrated Seismic Network. The deformation observations include Global Positioning System (GPS), Interferometric Synthetic Aperture Radar (InSAR), and borehole strainmeter recordings of the dynamic strain field. These data collectively constrain the kinematic coseismic slip distributions of the events, with measurements variously observing coseismic slip from one event (e.g., seismic waveforms, kinematic solutions from continuous GPS, and strainmeter time series) or coseismic slip from both events combined (InSAR). We find that the foreshock ruptured two separate faults, one with left‐lateral strike slip on a northeast–southwest‐trending fault and the other with right‐lateral strike slip on an orthogonal fault, with unilateral rupture propagation along both. The mainshock ruptured a series of northwest–southeast‐trending faults with right‐lateral strike slip concentrated in the uppermost 6 km with exceptionally low‐rupture velocity averaging 1.0 &#x2013; 1.5 &#x2009;&#x2009; km / s ">

California↗